Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

On-the-Fly Treatment of Temperature Dependent Cross Sections in the Unresolved Resonance Regions in RMC Code

Monte Carlo method can provide high fidelity neutronics analysis of different nuclear reactors, owing to its advantages of the flexible geometry modeling and the use of continuous-energy nuclear cross sections. It can also be coupled with thermal-hydraulics codes to consider the feedback. For thermal reactors such as PWR and HTGR, two main temperature effects should be considered. The first is the effect of thermal motion of target nuclei in resolved resonance energy regions, which is known as Doppler effect. The second is the thermal scattering and bound effect in thermal energy regions. The traditional approach of pre-generated cross sections has difficulty in memory footprint for detailed temperature modelling in multi-physics calculations. Recently, on-the-fly (OTF) technique has been proposed in order to reduce the memory consumption for both resolved resonance energy and thermal energy. For fast reactors calculations, the temperature dependence of cross sections in the unresolved resonance energy regions is also important. However, less attention was paid to the unresolved resonance regions compared with resolved resonance energy and thermal energy. In this paper, the on-the-fly temperature treatment of cross sections in unresolved resonance regions was proposed and developed in RMC codes. The results show that the on-the-fly treatment has high efficiency and satisfactory fidelity. The on-the-fly treatment of unresolved resonance region was also combined with the target motion sampling (TMS) method for resolved resonance energy. The formats of libraries were changed, in which the probability tables were separated from the cross sections. Then the interpolation of probability table was performed on-the-fly during the transport routine. The results show that the on-the-fly treatment has high efficiency and satisfactory fidelity. The on-the-fly treatment of unresolved resonance region was also combined with the target motion sampling (TMS) method for resolved resonance energy, showing that TMS method is effective and accurate for cross sections treatment in resolved resonance energy regions. The memory saving of the on-the-fly treatment is also considerable. For the high-fidelity 'transport-burnup-thermal-hydraulics' coupling calculation, the memory requirement is 60 Mb. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

THE FINAL-PARSEC PROBLEM IN THE COLLISIONLESS LIMIT

A binary supermassive black hole loses energy via ejection of stars in a galactic nucleus, until emission of gravitational waves becomes strong enough to induce rapid coalescence. Evolution via the gravitational slingshot requires that stars be continuously supplied to the binary, and it is known that in spherical galaxies the reservoir of such stars is quickly depleted, leading to stalling of the binary at parsec-scale separations. Recent N-body simulations of galaxy mergers and isolated nonspherical galaxies suggest that this stalling may not occur in less idealized systems. However, it remains unclear to what degree these conclusions are affected by collisional relaxation, which is much stronger in the numerical simulations than in real galaxies. In this study, we present a novel Monte Carlo method that can efficiently deal with both collisional and collisionless dynamics, and with galaxy models having arbitrary shapes. We show that without relaxation, the final-parsec problem may be overcome only in triaxial galaxies. Axisymmetry is not enough, but even a moderate departure from axisymmetry is sufficient to keep the binary shrinking. We find that the binary hardening rate is always substantially lower than the maximum possible, “full-loss-cone” rate, and that it decreases with time, but that stellar-dynamical interactions are nevertheless able to drive the binary to coalescence on a timescale ≲1 Gyr in any triaxial galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Particle-image Velocimetry Sensitivity through Automatic Differentiation

A commonly encountered difficulty in modern experimental fluid dynamics experiments is the need for full-flow visualization and measurement; traditional techniques such as hot-wire probes and laser Doppler velocimetry cannot provide such information without difficult scanning procedures. Particle-image velocimetry (PIV) is a measurement technique which can overcome this spatial limitation to provide real-time three-dimensional fluid velocity vectors for an entire plane within the flow, as opposed to a single point that is provided by other measurement techniques. To accomplish this, a PIV system seeds the fluid flow with tracer particles; these must have nearly the same density as that of the fluid and a high drag coefficient to ensure that they follow the fluid flow without deviation. A laser is sent through optics to expand its beam into a sheet which illuminates the plane of interest in the fluid. The particles reflect more light than the fluid, and their image is recorded using a high speed camera. Due to limitations in bandwidth, video cannot be taken fast enough to be useful, and instead pairs of images are captured with a known separation in time. The illuminated particles will move between the two images taken by the camera, and their relative motion can be used to compute the fluid velocity at that point. When a single camera is used, only two-dimensional velocity vectors can be computed, but with the addition of another camera this can be improved to full three-dimensional vectors. The flow's time evolution can be measured by taking series of image pairs and processing each pair to get the velocity field at the time between each image pair. Only the two-dimensional version will be considered in this work. To compute the velocity vectors of the flow, the images are broken up into interrogation regions, usually in a uniform grid. In the most basic PIV algorithm, corresponding regions in the image pairs are cross-correlated to compute the most likely displacement from one image to the next. This is used since the cross correlation will achieve a maximum value when the two interrogation regions line up the best. The location of the peak in the correlation map and scaling information for the optical setup is used to compute the particle displacement from one image to the next. This process is similar to holding two semi-transparent images and moving one until they match as well as possible. Sub-pixel displacements between images can be obtained by fitting a Gaussian function to the correlation peak. Knowing the displacement vector for each set of interrogation regions in the image pair, it is a simple matter to divide by the time gap between the two images to compute the fluid velocity at the point in time and space. Since there are many interrogation regions in each image pair, the velocity field for a sheet is sampled according to the arrangement of the interrogation regions. The sensitivity of the final computed vector to the input images is difficult to compute, and has led most considerations of uncertainty in the PIV field to apply Monte Carlo methods and correlations to quantify the systematic and random errors which result from application of various PIV algorithms. In this work, automatic differentiation (AD) has been used to trace the dependence of velocity vectors computed using the correlation-peak PIV algorithm back to the input values given to an artificial image generator. These initial results form the foundation for continued work in the application of AD to algorithms in PIV. They demonstrate that AD methods can be successfully applied to PIV algorithms without undue complexity. The derivatives of each individual vector with regard to the coefficients of the velocity field used to create the artificial images can be computed and the statistics of them considered collectively, or correlated between each other to quantify the effectiveness of a PIV algorithm. With the initial implementation complete, adding extra algorithms is relatively simple and considerable data can be generated for analysis. It is hoped that application of these methods can lead to improved PIV algorithms which demonstrate better sensitivity to small variations in particle displacement. While these methods can be used as an analysis tool for refinement of new algorithms, computation of actual resultant vector sensitivities to input images' individual pixels is the next logical step of this work. Such information will be analyzed with experimental data to produce a better understanding of how output vectors' uncertainty relates to the actual image pairs, and lead to automatic bad vector removal. (authors)

42 ENGINEERING↗
DOE OSTI2014

On geometric factors for neutral particle analyzers

Neutral particle analyzers (NPA) detect neutralized energetic particles that escape from plasmas. Geometric factors relate the counting rate of the detectors to the intensity of the particle source. Accurate geometric factors enable quick simulation of geometric effects without the need to resort to slower Monte Carlo methods. Previously derived expressions [G. R. Thomas and D. M. Willis, “Analytical derivation of the geometric factor of a particle detector having circular or rectangular geometry,” J. Phys. E: Sci. Instrum. 5(3), 260 (1972); J. D. Sullivan, “Geometric factor and directional response of single and multi-element particle telescopes,” Nucl. Instrum. Methods 95(1), 5–11 (1971)] for the geometric factor implicitly assume that the particle source is very far away from the detector (far-field); this excludes applications close to the detector (near-field). The far-field assumption does not hold in most fusion applications of NPA detectors. We derive, from probability theory, a generalized framework for deriving geometric factors that are valid for both near and far-field applications as well as for non-isotropic sources and nonlinear particle trajectories.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Improved modeling of two-dimensional transitions in dense phases on crystalline surfaces. Krypton–graphite system

This paper presents a refined technique to describe two-dimensional phase transitions in dense fluids adsorbed on a crystalline surface. Prediction of parameters of 2D liquid–solid equilibrium is known to be an extremely challenging problem, which is mainly due to a small difference in thermodynamic functions of coexisting phases and lack of accuracy of numerical experiments in case of their high density. This is a serious limitation of various attempts to circumvent this problem. To improve this situation, a new methodology based on the kinetic Monte Carlo method was applied. The methodology involves analysis of equilibrium gas–liquid and gas–solid systems undergoing an external potential, which allows gradual shifting parameters of the phase coexistence. The interrelation of the chemical potential and tangential pressure for each system is then treated with the Gibbs–Duhem equation to obtain the point of intersection corresponding to the liquid/solid–solid equilibrium coexistence. The methodology is demonstrated on the krypton–graphite system below and above the 2D critical temperature. Using experimental data on the liquid–solid and the commensurate–incommensurate transitions in the krypton monolayer derived from adsorption isotherms, the Kr–graphite Lennard–Jones parameters have been corrected resulting in a higher periodic potential modulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

First-principles equation of state and electronic properties of warm dense oxygen

We perform all-electron path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of oxygen. Our simulations cover a wide density-temperature range of 1–100 g cm{sup −3} and 10{sup 4}–10{sup 9} K. By combining results from PIMC and DFT-MD, we are able to compute pressures and internal energies from first-principles at all temperatures and provide a coherent equation of state. We compare our first-principles calculations with analytic equations of state, which tend to agree for temperatures above 8 × 10{sup 6} K. Pair-correlation functions and the electronic density of states reveal an evolving plasma structure and ionization process that is driven by temperature and density. As we increase the density at constant temperature, we find that the ionization fraction of the 1s state decreases while the other electronic states move towards the continuum. Finally, the computed shock Hugoniot curves show an increase in compression as the first and second shells are ionized.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Matrix-assisted relaxation in Fe(phen){sub 2}(NCS){sub 2} spin-crossover microparticles, experimental and theoretical investigations

In this study, we present the influence of the embedding matrix on the relaxation of Fe(phen){sub 2}(NCS){sub 2} (phen = 1,10-phenanthroline) spin-transition microparticles as revealed by experiments and provide an explanation within the framework of an elastic model based on a Monte-Carlo method. Experiments show that the shape of the high-spin → low-spin relaxation curves is drastically changed when the particles are dispersed in glycerol. This effect was considered in the model by means of interactions between the microparticles and the matrix. A faster start of the relaxation for microparticles embedded in glycerol is due to an initial positive local pressure acting on the edge spin-crossover molecules from the matrix side. This local pressure diminishes and eventually becomes negative during relaxation, as an effect of the decrease of the volume of spin-crossover microparticles from high-spin to low-spin.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗
DOE OSTI2018

Convergence and Efficiency of Adaptive Importance Sampling Techniques with Partial Biasing

We propose a new Monte Carlo method to efficiently sample a multimodal distribution (known up to a normalization constant). We consider a generalization of the discrete-time Self Healing Umbrella Sampling method, which can also be seen as a generalization of well-tempered metadynamics. The dynamics is based on an adaptive importance technique. The importance function relies on the weights (namely the relative probabilities) of disjoint sets which form a partition of the space. These weights are unknown but are learnt on the fly yielding an adaptive algorithm. In the context of computational statistical physics, the logarithm of these weights is, up to an additive constant, the free-energy, and the discrete valued function defining the partition is called the collective variable. The algorithm falls into the general class of Wang–Landau type methods, and is a generalization of the original Self Healing Umbrella Sampling method in two ways: (i) the updating strategy leads to a larger penalization strength of already visited sets in order to escape more quickly from metastable states, and (ii) the target distribution is biased using only a fraction of the free-energy, in order to increase the effective sample size and reduce the variance of importance sampling estimators. We prove the convergence of the algorithm and analyze numerically its efficiency on a toy example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Ion energy-angle distribution functions at the plasma-material interface in oblique magnetic fields

The ion energy-angle distribution (IEAD) at the wall of a magnetized plasma is of fundamental importance for the determination of the material processes occurring at the plasma-material interface, comprising secondary emissions and material sputtering. Here, we present a numerical characterization of the IEAD at the wall of a weakly collisional magnetized plasma with the magnetic field inclined at an arbitrary angle with respect to the wall. The analysis has been done using two different techniques: (1) a fluid-Monte Carlo method, and (2) particle-in-cell simulations, the former offering a fast but approximate method for the determination of the IEADs, the latter giving a computationally intensive but self-consistent treatment of the plasma behavior from the quasi-neutral region to the material boundary. The two models predict similar IEADs, whose similarities and differences are discussed. Data are presented for magnetic fields inclined at angles from normal to grazing incidence (0°–85°). We show the scaling factors of the average and peak ion energy and trends of the pitch angle at the wall as a function of the magnetic angle, for use in the correlation of fluid plasma models to material models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Phase transitions and critical properties in the antiferromagnetic Ising model on a layered triangular lattice with allowance for intralayer next-nearest-neighbor interactions

The phase transitions (PTs) and critical properties of the antiferromagnetic Ising model on a layered (stacked) triangular lattice have been studied by the Monte Carlo method using a replica algorithm with allowance for the next-nearest-neighbor interactions. The character of PTs is analyzed using the histogram technique and the method of Binder cumulants. It is established that the transition from the disordered to paramagnetic phase in the adopted model is a second-order PT. Static critical exponents of the heat capacity (α), susceptibility (γ), order parameter (β), and correlation radius (ν) and the Fischer exponent η are calculated using the finite-size scaling theory. It is shown that (i) the antiferromagnetic Ising model on a layered triangular lattice belongs to the XY universality class of critical behavior and (ii) allowance for the intralayer interactions of next-nearest neighbors in the adopted model leads to a change in the universality class of critical behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Historical remarks on exponential product and quantum analysis

The exponential product formula [1, 2] was substantially introduced in physics by the present author [2]. Its systematic applications to quantum Monte Carlo Methods [3] were preformed [4, 5] first in 1977. Many interesting applications [6] of the quantum-classical correspondence (namely S-T transformation) have been reported. Systematic higher-order decomposition formulae were also discovered by the present author [7-11], using the recursion scheme [7, 9]. Physically speaking, these exponential product formulae play a conceptual role of separation of procedures [3,14]. Mathematical aspects of these formulae have been integrated in quantum analysis [15], in which non-commutative differential calculus is formulated and a general quantum Taylor expansion formula is given. This yields many useful operator expansion formulae such as the Feynman expansion formula and the resolvent expansion. Irreversibility and entropy production are also studied using quantum analysis [15].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Metis aboard the Solar Orbiter space mission: Doses from galactic cosmic rays and solar energetic particles

The aim of this work is to calculate the dose released by galactic cosmic rays (GCRs) and solar energetic particles (SEPs) in the polarimeter of the Multi Element Telescope for Imaging and Spectroscopy (METIS) coronagraph [1] aboard the Solar Orbiter. This investigation is performed with a Monte Carlo method by considering the role of SEP events of proper intensity at a heliocentric distance from the Sun averaged along the spacecraft orbit. Our approach can be extended to other space missions reaching short distances from the Sun, such as Solar Probe Plus. This study indicates that the deposited dose on the whole set of polarimeter lenses and filters during ten years of the Solar Orbiter mission is of about 2000 Gy. For cerium treated lenses, a dose of 10{sup 6} Gy of gamma radiation from a {sup 60}Co source causes a few percent transmittance loss.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Rotational Brownian Dynamics simulations of clathrin cage formation

The self-assembly of nearly rigid proteins into ordered aggregates is well suited for modeling by the patchy particle approach. Patchy particles are traditionally simulated using Monte Carlo methods, to study the phase diagram, while Brownian Dynamics simulations would reveal insights into the assembly dynamics. However, Brownian Dynamics of rotating anisotropic particles gives rise to a number of complications not encountered in translational Brownian Dynamics. We thoroughly test the Rotational Brownian Dynamics scheme proposed by Naess and Elsgaeter [Macromol. Theory Simul. 13, 419 (2004); Naess and Elsgaeter Macromol. Theory Simul. 14, 300 (2005)], confirming its validity. We then apply the algorithm to simulate a patchy particle model of clathrin, a three-legged protein involved in vesicle production from lipid membranes during endocytosis. Using this algorithm we recover time scales for cage assembly comparable to those from experiments. We also briefly discuss the undulatory dynamics of the polyhedral cage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Crossover-Induced Spin Fluctuation and Electron Pairing in Strongly Correlated Electrons

We investigate the two-dimensional (2D) Hubbard model and d-p model as a model for cuprate high-temperature superconductors, by employing an optimized variational Monte Carlo method. We improved the many-body wave function starting from the simple Gutzwiller function. On the basis of optimized variational wave functions, we show that there is a superconducting phase in the ground state of the 2D Hubbard model. The on-site repulsive Coulomb interaction U can bring about high-temperature superconductivity because its characteristic energy is of the order of eV. The attractive interaction between electrons is induced by the crossover from weakly to strongly correlated regions. High-temperature superconductivity is possible in the strongly correlated region, where U is greater than the bandwidth.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Analysis of the Effect of Different Initial States and Structural Defects on the Characteristics of the Nonequilibrium Critical Behavior of the 3D Ising Model

The effect of different initial values m{sub 0} of magnetization and structural defects on the nonequilibrium critical behavior of the 3D Ising model have been analyzed numerically using the Monte Carlo method. Analysis of the two-time dependences of the autocorrelation function and dynamic susceptibility has revealed a substantial influence of the initial states on the aging effects that are characterized by anomalous retardation of relaxation and correlation in the system upon an increase in the waiting time. We have studied the violations of the fluctuation–dissipation theorem and calculated the limiting fluctuation–dissipation ratio. It is shown that in the nonequilibrium critical behavior of the 3D Ising model, two universality subclasses corresponding to the evolution of the system from the high-temperature (with m{sub 0} = 0) and low-temperature (with m{sub 0} = 1) initial states with the values of the limiting fluctuation–dissipation ratio typical of these states can be singled out.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Carbon Abundances in Starburst Galaxies of the Local Universe

The cosmological origin of carbon, the fourth most abundant element in the universe, is not well known and a matter of heavy debate. We investigate the behavior of C/O to O/H in order to constrain the production mechanism of carbon. We measured emission-line intensities in the spectral range from 1600 to 10000 Å on Space Telescope Imaging Spectrograph (STIS) long-slit spectra of 18 starburst galaxies in the local universe. We determined chemical abundances through traditional nebular analysis, and we used a Markov Chain Monte Carlo method to determine where our carbon and oxygen abundances lie in the parameter space. We conclude that our C and O abundance measurements are sensible. We analyzed the behavior of our sample in the [C/O] versus [O/H] diagram with respect to other objects such as DLAs, neutral ISM measurements, and disk and halo stars, finding that each type of object seems to be located in a specific region of the diagram. Our sample shows a steeper C/O versus O/H slope with respect to other samples, suggesting that massive stars contribute more to the production of C than N at higher metallicities, only for objects where massive stars are numerous; otherwise, intermediate-mass stars dominate the C and N production.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Variational Monte Carlo calculations of lithium atom in strong magnetic field

The variational Monte Carlo method is applied to investigate the ground state and some excited states of the lithium atom and its ions up to Z = 10 in the presence of an external magnetic field regime with γ = 0–100 arb. units. The effect of increasing field strength on the ground state energy is studied and precise values for the crossover field strengths were obtained. Our calculations are based on using accurate forms of trial wave functions, which were put forward in calculating energies in the absence of magnetic field. Furthermore, the value of Y at which ground-state energy of the lithium atom approaches to zero was calculated. The obtained results are in good agreement with the most recent values and also with the exact values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗