Monte Carlo method with complex weights for neutron leakage-corrected calculations and anisotropic diffusion coefficient generations

Purpose This paper aims to present development and application of the Bayesian inverse approach for retrieving parameters of non-linear diffusion coefficient based on the integral information. Design/methodology/approach The Bayes formula was used to construct posterior distribution of the unknown parameters of non-linear diffusion coefficient. The resulting aposteriori distribution of sought parameters was integrated using Markov Chain Monte Carlo method to obtain expected values of estimated diffusivity parameters as well as their confidence intervals. Unsteady non-linear diffusion equation was discretised with the Global Radial Basis Function Collocation method and solved in time using Crank–Nicholson technique. Findings A number of manufactured analytical solutions of the non-linear diffusion problem was used to verify accuracy of the developed inverse approach. Reasonably good agreement, even for highly correlated parameters, was obtained. Therefore, the technique was used to compute concentration dependent diffusion coefficient of water in paper. Originality/value An original inverse technique, which couples efficiently meshless solution of the diffusion problem with the Bayesian inverse methodology, is presented in the paper. This methodology was extensively verified and applied to the real-life problem.

Download Full-text

Outbursts of comet Schwassmann-Wachmann 1, and the cloud of interplanetary boulders

International Astronomical Union Colloquium ◽

10.1017/s0252921100034035 ◽

1974 ◽

Vol 22 ◽

pp. 307 ◽

Cited By ~ 3

Author(s):

Zdenek Sekanina

Keyword(s):

Monte Carlo ◽

Monte Carlo Method ◽

Interplanetary Space ◽

Porous Matrix ◽

Maximum Diameter ◽

Expansion Velocity ◽

Solid Constituent ◽

Meteoric Material ◽

Periodic Comet

AbstractIt is suggested that the outbursts of Periodic Comet Schwassmann-Wachmann 1 are triggered by impacts of interplanetary boulders on the surface of the comet’s nucleus. The existence of a cloud of such boulders in interplanetary space was predicted by Harwit (1967). We have used the hypothesis to calculate the characteristics of the outbursts – such as their mean rate, optically important dimensions of ejected debris, expansion velocity of the ejecta, maximum diameter of the expanding cloud before it fades out, and the magnitude of the accompanying orbital impulse – and found them reasonably consistent with observations, if the solid constituent of the comet is assumed in the form of a porous matrix of lowstrength meteoric material. A Monte Carlo method was applied to simulate the distributions of impacts, their directions and impact velocities.

Download Full-text

Structures and crystallization of vacuum-deposited amorphous Se and Sb films

Proceedings, annual meeting, Electron Microscopy Society of America ◽

10.1017/s0424820100173832 ◽

1990 ◽

Vol 48 (4) ◽

pp. 140-141

Author(s):

Makoto Shiojiri ◽

Toshiyuki Isshiki ◽

Tetsuya Fudaba ◽

Yoshihiro Hirota

Keyword(s):

Monte Carlo ◽

Electron Microscope ◽

Monte Carlo Method ◽

Computer Program ◽

Radial Distribution ◽

Film Formation ◽

Amorphous State ◽

Two Dimensional ◽

Amorphous Films ◽

Binding Condition

In hexagonal Se crystal each atom is covalently bound to two others to form an endless spiral chain, and in Sb crystal each atom to three others to form an extended puckered sheet. Such chains and sheets may be regarded as one- and two- dimensional molecules, respectively. In this paper we investigate the structures in amorphous state of these elements and the crystallization.HRTEM and ED images of vacuum-deposited amorphous Se and Sb films were taken with a JEM-200CX electron microscope (Cs=1.2 mm). The structure models of amorphous films were constructed on a computer by Monte Carlo method. Generated atoms were subsequently deposited on a space of 2 nm×2 nm as they fulfiled the binding condition, to form a film 5 nm thick (Fig. 1a-1c). An improvement on a previous computer program has been made as to realize the actual film formation. Radial distribution fuction (RDF) curves, ED intensities and HRTEM images for the constructed structure models were calculated, and compared with the observed ones.

Download Full-text