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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Studies of concentration and temperature dependences of precipitation kinetics in iron-copper alloys using kinetic Monte Carlo and stochastic statistical simulations

The previously developed ab initio model and the kinetic Monte Carlo method (KMCM) are used to simulate precipitation in a number of iron-copper alloys with different copper concentrations x and temperatures T. The same simulations are also made using an improved version of the previously suggested stochastic statistical method (SSM). The results obtained enable us to make a number of general conclusions about the dependences of the decomposition kinetics in Fe-Cu alloys on x and T. We also show that the SSM usually describes the precipitation kinetics in good agreement with the KMCM, and using the SSM in conjunction with the KMCM allows extending the KMC simulations to the longer evolution times. The results of simulations seem to agree with available experimental data for Fe-Cu alloys within statistical errors of simulations and the scatter of experimental results. Comparison of simulation results with experiments for some multicomponent Fe-Cu-based alloys allows making certain conclusions about the influence of alloying elements in these alloys on the precipitation kinetics at different stages of evolution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Dynamic heterogeneity in crossover spin facilitated model of supercooled liquid and fractional Stokes-Einstein relation

Kinetically constrained models have gained much interest as models that assign the origins of interesting dynamic properties of supercooled liquids to dynamical facilitation mechanisms that have been revealed in many experiments and numerical simulations. In this work, we investigate the dynamic heterogeneity in the fragile-to-strong liquid via Monte Carlo method using the model that linearly interpolates between the strong liquid-like behavior and the fragile liquid-like behavior by an asymmetry parameter b. When the asymmetry parameter is sufficiently small, smooth fragile-to-strong transition is observed both in the relaxation time and the diffusion constant. Using these physical quantities, we investigate fractional Stokes-Einstein relations observed in this model. When b is fixed, the system shows constant power law exponent under the temperature change, and the exponent has the value between that of the Frederickson-Andersen model and the East model. Furthermore, we investigate the dynamic length scale of our systems and also find the crossover relation between the relaxation time. We ascribe the competition between energetically favored symmetric relaxation mechanism and entropically favored asymmetric relaxation mechanism to the fragile-to-strong crossover behavior.

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

HIGH-TEMPERATURE PROCESSING OF SOLIDS THROUGH SOLAR NEBULAR BOW SHOCKS: 3D RADIATION HYDRODYNAMICS SIMULATIONS WITH PARTICLES

A fundamental, unsolved problem in solar system formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks have been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals and planetary embryos on eccentric orbits can produce bow shocks as they move supersonically through the disk gas, and are one possible source of chondrule-melting shocks. We investigate chondrule formation in bow shocks around planetoids through three-dimensional radiation hydrodynamics simulations. A new radiation transport algorithm that combines elements of flux-limited diffusion and Monte Carlo methods is used to capture the complexity of radiative transport around bow shocks. An equation of state that includes the rotational, vibrational, and dissociation modes of H{sub 2} is also used. Solids are followed directly in the simulations and their thermal histories are recorded. Adiabatic expansion creates rapid cooling of the gas, and tail shocks behind the embryo can cause secondary heating events. Radiative transport is efficient, and bow shocks around planetoids can have luminosities ∼few× 10{sup –8} L{sub ☉}. While barred and radial chondrule textures could be produced in the radiative shocks explored here, porphyritic chondrules may only be possible in the adiabatic limit. We present a series of predicted cooling curves that merit investigation in laboratory experiments to determine whether the solids produced by bow shocks are represented in the meteoritic record by chondrules or other solids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE INITIAL MASS FUNCTION AND THE SURFACE DENSITY PROFILE OF NGC 6231

We have performed new wide-field photometry of the young open cluster NGC 6231 to study the shape of the initial mass function (IMF) and mass segregation. We also investigated the reddening law toward NGC 6231 from optical to mid-infrared color excess ratios, and found that the total-to-selective extinction ratio is R{sub V} = 3.2, which is very close to the normal value. But many early-type stars in the cluster center show large color excess ratios. We derived the surface density profiles of four member groups, and found that they reach the surface density of field stars at about 10', regardless of stellar mass. The IMF of NGC 6231 is derived for the mass range 0.8-45 M{sub Sun }. The slope of the IMF of NGC 6231 ({Gamma} = -1.1 {+-} 0.1) is slightly shallower than the canonical value, but the difference is marginal. In addition, the mass function varies systematically, and is a strong function of radius-it is very shallow at the center, and very steep at the outer ring suggesting the cluster is mass segregated. We confirm the mass segregation for the massive stars (m {approx}> 8 M{sub Sun }) by a minimum spanning tree analysis. Using a Monte Carlo method, we estimate the total mass of NGC 6231 to be about 2.6 ({+-} 0.6) Multiplication-Sign 10{sup 3} M{sub Sun }. We constrain the age of NGC 6231 by comparison with evolutionary isochrones. The age of the low-mass stars ranges from 1 to 7 Myr with a slight peak at 3 Myr. However, the age of the high-mass stars depends on the adopted models and is 3.5 {+-} 0.5 Myr from the non-rotating or moderately rotating models of Brott et al. as well as the non-rotating models of Ekstroem et al. But the age is 4.0-7.0 Myr if the rotating models of Ekstroem et al. are adopted. This latter age is in excellent agreement with the timescale of ejection of the high-mass runaway star HD 153919 from NGC 6231, albeit the younger age cannot be entirely excluded.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The rigorous stochastic matrix multiplication scheme for the calculations of reduced equilibrium density matrices of open multilevel quantum systems

Understanding the roles of the temporary and spatial structures of quantum functional noise in open multilevel quantum molecular systems attracts a lot of theoretical interests. I want to establish a rigorous and general framework for functional quantum noises from the constructive and computational perspectives, i.e., how to generate the random trajectories to reproduce the kernel and path ordering of the influence functional with effective Monte Carlo methods for arbitrary spectral densities. This construction approach aims to unify the existing stochastic models to rigorously describe the temporary and spatial structure of Gaussian quantum noises. In this paper, I review the Euclidean imaginary time influence functional and propose the stochastic matrix multiplication scheme to calculate reduced equilibrium density matrices (REDM). In addition, I review and discuss the Feynman-Vernon influence functional according to the Gaussian quadratic integral, particularly its imaginary part which is critical to the rigorous description of the quantum detailed balance. As a result, I establish the conditions under which the influence functional can be interpreted as the average of exponential functional operator over real-valued Gaussian processes for open multilevel quantum systems. I also show the difference between the local and nonlocal phonons within this framework. With the stochastic matrix multiplication scheme, I compare the normalized REDM with the Boltzmann equilibrium distribution for open multilevel quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Susceptibility Critical Exponent of a 1D Ising Ring-Type Ferromagnetic

Using the Monte Carlo method, critical behavior of the one-dimensional ferromagnetic Ising model has been investigated with allowance for the interaction of the second and third neighbors and four-particle interaction. The obtained results on the critical temperature were compared with the critical temperature of the quasi-one-dimensional Ising magnetic [(CH{sub 3}){sub 3}NH] · FeCl{sub 3} · 2H{sub 2}O and with the magnitude of the exchange interaction J/k{sub B} = 17.4 K. Within the scope of the finite-dimensional scaling theory, the critical susceptibility exponent has been calculated. It has been shown that values of the susceptibility exponent for the one-dimensional Ising model with periodic boundary conditions are considerably less than the known values of the exponents for three-dimensional systems. The critical susceptibility exponent strongly depends on energy parameters; namely, it decreases with an increase in them.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Electron mobility calculation for graphene on substrates

By a semiclassical Monte Carlo method, the electron mobility in graphene is calculated for three different substrates: SiO{sub 2}, HfO{sub 2}, and hexagonal boron nitride (h-BN). The calculations account for polar and non-polar surface optical phonon (OP) scatterings induced by the substrates and charged impurity (CI) scattering, in addition to intrinsic phonon scattering in pristine graphene. It is found that HfO{sub 2} is unsuitable as a substrate, because the surface OP scattering of the substrate significantly degrades the electron mobility. The mobility on the SiO{sub 2} and h-BN substrates decreases due to CI scattering. However, the mobility on the h-BN substrate exhibits a high electron mobility of 170 000 cm{sup 2}/(V·s) for electron densities less than 10{sup 12 }cm{sup −2}. Therefore, h-BN should be an appealing substrate for graphene devices, as confirmed experimentally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

AUTONOMOUS GAUSSIAN DECOMPOSITION

We present a new algorithm, named Autonomous Gaussian Decomposition (AGD), for automatically decomposing spectra into Gaussian components. AGD uses derivative spectroscopy and machine learning to provide optimized guesses for the number of Gaussian components in the data, and also their locations, widths, and amplitudes. We test AGD and find that it produces results comparable to human-derived solutions on 21 cm absorption spectra from the 21 cm SPectral line Observations of Neutral Gas with the EVLA (21-SPONGE) survey. We use AGD with Monte Carlo methods to derive the H i line completeness as a function of peak optical depth and velocity width for the 21-SPONGE data, and also show that the results of AGD are stable against varying observational noise intensity. The autonomy and computational efficiency of the method over traditional manual Gaussian fits allow for truly unbiased comparisons between observations and simulations, and for the ability to scale up and interpret the very large data volumes from the upcoming Square Kilometer Array and pathfinder telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Charge exchange and ionization cross sections of H{sup +}+H collision in dense quantum plasmas

The plasma screening effects of dense quantum plasmas on H{sup +}+H charge exchange and ionization cross sections are calculated by the classical trajectory Monte Carlo method. For charge exchange cross sections, it is found that the screening effects reduce cross sections slightly in weak screening conditions. However, cross sections are reduced substantially in strong screening conditions. For ionization cross sections, with the increase of screening effects, cross sections for low energies increase more rapidly than those for high energies. When the screening effects are strong enough, it is found that ionization cross sections decrease with the increase of incident H{sup +} energy. In addition, the cross sections have been compared with those in weakly coupled plasmas. It is found that in weak screening conditions, plasma screening effects in the two plasmas are approximately the same, while in strong screening conditions, screening effects of dense quantum plasmas are stronger than those of weakly coupled plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Extragalactic synchrotron transients in the era of wide-field radio surveys. I. Detection rates and light curve characteristics

The impending era of wide-field radio surveys has the potential to revolutionize our understanding of astrophysical transients. Here we evaluate the prospects of a wide range of planned and hypothetical radio surveys using the properties and volumetric rates of known and hypothetical classes of extragalactic synchrotron radio transients (e.g., on-axis and off-axis gamma-ray bursts (GRBs), supernovae, tidal disruption events, compact object mergers). Utilizing these sources and physically motivated considerations we assess the allowed phase space of radio luminosity and peak timescale for extragalactic transients. We also include for the first time effects such as redshift evolution of the rates, K-corrections, and non-Euclidean luminosity distance, which affect the detection rates of the most sensitive surveys. The number of detected events is calculated by means of a Monte Carlo method, using the various survey properties (depth, cadence, area) and realistic detection criteria that include a cut on the minimum variability of the transients during the survey and an assessment of host galaxy contamination. We find that near-term GHz frequency surveys (ASKAP/VAST, Very Large Array Sky Survey) will detect few events: ≲30−50 on- and off-axis long GRBs (LGRBs) and off-axis tidal disruption events, and ∼50−100 neutron star binary mergers if ∼0.5% of the mergers result in a stable millisecond magnetar. Low-frequency surveys (e.g., LOFAR) are unlikely to detect any transients, while a hypothetical large-scale mm survey may detect ∼40 on-axis LGRBs. On the other hand, we find that SKA1 surveys at ∼0.1−1 GHz have the potential to uncover thousands of transients, mainly on-axis and off-axis LGRBs, on-axis short GRBs, off-axis TDEs, and neutron star binary mergers with magnetar remnants.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Monte Carlo study of vibrational relaxation processes

A new model is proposed for the computation of vibrational nonequilibrium in the direct simulation Monte Carlo method (DSMC). This model permits level to level vibrational transitions for the first time in a Monte Carlo flowfield simulation. The model follows the Landau-Teller theory for a harmonic oscillator in which the rates of transition are related to an experimental correlation for the vibrational relaxation time. The usual method for simulating such processes in the DSMC technique applies a constant exchange probability to each collision and the vibrational energy is treated as a continuum. A comparison of these two methods is made for the flow of nitrogen over a wedge. Significant differences exist for the vibrational temperatures computed. These arise as a consequence of the incorrect application of a constant exchange probability in the old method. It is found that the numerical performances of the two vibrational relaxation models are equal.

Boyd, Iain D.↗
DOE OSTI2018

Aging Effects in the Nonequilibrium Behavior of Magnetic Superstructures and Their Manifestation in Magnetoresistance

The nonequilibrium behavior of the magnetic superstructures consisting of alternating magnetic and nonmagnetic nanolayers is numerically simulated by Monte Carlo methods. An analysis of the calculated two-time dependence of an autocorrelation function during the evolution from a high-temperature initial state has revealed aging effects, which are characterized by slowing down of the correlation effects in the system when the waiting time increases. In contrast to bulk magnetic systems, the aging effects are shown to appear in magnetic superstructures both near the critical ferromagnetic ordering temperature T{sub c} in films and in a low-temperature phase at T ≤ T{sub c}. The aging effects in the correlation processes in a magnetic multilayer structure weaken when ferromagnetic film thickness N increases at T = T{sub c}(N), and these effects increase with film thickness N at temperatures T = T{sub c}(N)/2. When simulating the transport properties of the Co/Cu(001)/Co structure, we calculated the temperature dependence of equilibrium magnetoresistance and were the first to reveal the influence of nonequilibrium behavior of the structure on the magnetoresistance and the manifestation of the aging effects in it.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2016

Phase transitions in a two-dimensional antiferromagnetic Potts model on a triangular lattice with next-nearest neighbor interactions

Phase transitions (PTs) and frustrations in two-dimensional structures described by a three-vertex antiferromagnetic Potts model on a triangular lattice are investigated by the Monte Carlo method with regard to nearest and next-nearest neighbors with interaction constants J{sub 1} and J{sub 2}, respectively. PTs in these models are analyzed for the ratio r = J{sub 2}/J{sub 1} of next-nearest to nearest exchange interaction constants in the interval |r| = 0–1.0. On the basis of the analysis of the low-temperature entropy, the density of states function of the system, and the fourth-order Binder cumulants, it is shown that a Potts model with interaction constants J{sub 1} < 0 and J{sub 2} < 0 exhibits a first-order PT in the range of 0 ⩽ r < 0.2, whereas, in the interval 0.2 ⩽ r ⩽ 1.0, frustrations arise in the system. At the same time, for J{sub 1} > 0 and J{sub 2} < 0, frustrations arise in the range 0.5 < |r| < 1.0, while, in the interval 0 ⩽ |r| ⩽ 1/3, the model exhibits a second-order PT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

A comparative study of transport properties of monolayer graphene and AlGaN-GaN heterostructure

The electronic transport properties of monolayer graphene are presented with an Ensemble Monte Carlo method where a rejection technique is used to account for the occupancy of the final states after scattering. Acoustic and optic phonon scatterings are considered for intrinsic graphene and in addition, ionized impurity and surface roughness scatterings are considered for the case of dirty graphene. The effect of screening is considered in the ionized impurity scattering of electrons. The time dependence of drift velocity of carriers is obtained where overshoot and undershoot effects are observed for certain values of applied field and material parameters for intrinsic graphene. The field dependence of drift velocity of carriers showed negative differential resistance and disappeared as acoustic scattering becomes dominant for intrinsic graphene. The variation of electron mobility with temperature is calculated for intrinsic (suspended) and dirty monolayer graphene sheets separately and they are compared. These are also compared with the mobility of two dimensional electrons at an AlGaN/GaN heterostructure. It is observed that interface roughness may become very effective in limiting the mobility of electrons in graphene.

36 MATERIALS SCIENCE↗
NASA NTRS2016

A DSMC Surface Chemistry Model for Carbon-Based Ablators

A detailed molecular surface chemistry model for the DSMC (Direct Simulation Monte Carlo) method is proposed and implemented into the SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) DSMC solver. Molchanova et al. constructed a molecular model for surface recombination in DSMC that includes different surface processes (adsorption, desoprtion, Eley-Rideal and Langmuir-Hinshelwood). All surface processes can be divided into two groups: surface mechanisms, which involve only the particle adsorbed by the surface (desorption and Langmuir-Hinshelwood), and impact mechanisms, which also involve gas-phase particles (adsorption, Eley-Rideal). Using a similar approach, the 14-reaction kinetic model of oxygen-carbon interaction suggested by Zhlukhtov and Abe, as well as more recent models by Alba et al., Poovathinghal et al., and a new model developed in the scope of this work, are implemented in SPARTA. The computational results for the different oxidation models are compared with experimental results from Murray et al. (oxidation of a vitreous carbon surface due to a hyperthermal beam of O and O2), with a particular focus on fluxes, angular and Time-Of-Flight distributions of scattered particles.

oxidation↗
DOE OSTI2014

Dosimetric feasibility of real-time MRI-guided proton therapy

Purpose: Magnetic resonance imaging (MRI) is a prime candidate for image-guided radiotherapy. This study was designed to assess the feasibility of real-time MRI-guided proton therapy by quantifying the dosimetric effects induced by the magnetic field in patients’ plans and identifying the associated clinical consequences. Methods: Monte Carlo dose calculation was performed for nine patients of various treatment sites (lung, liver, prostate, brain, skull-base, and spine) and tissue homogeneities, in the presence of 0.5 and 1.5 T magnetic fields. Dose volume histogram (DVH) parameters such as D{sub 95}, D{sub 5}, and V{sub 20} as well as equivalent uniform dose were compared for the target and organs at risk, before and after applying the magnetic field. The authors further assessed whether the plans affected by clinically relevant dose distortions could be corrected independent of the planning system. Results: By comparing the resulting dose distributions and analyzing the respective DVHs, it was determined that despite the observed lateral beam deflection, for magnetic fields of up to 0.5 T, neither was the target coverage jeopardized nor was the dose to the nearby organs increased in all cases except for prostate. However, for a 1.5 T magnetic field, the dose distortions were more pronounced and of clinical concern in all cases except for spine. In such circumstances, the target was severely underdosed, as indicated by a decrease in D{sub 95} of up to 41% of the prescribed dose compared to the nominal situation (no magnetic field). Sites such as liver and spine were less affected due to higher tissue homogeneity, typically smaller beam range, and the choice of beam directions. Simulations revealed that small modifications to certain plan parameters such as beam isocenter (up to 19 mm) and gantry angle (up to 10°) are sufficient to compensate for the magnetic field-induced dose disturbances. The authors’ observations indicate that the degree of required corrections strongly depends on the beam range and direction relative to the magnetic field. This method was also applicable to more heterogeneous scenarios such as skull-base tumors. Conclusions: This study confirmed the dosimetric feasibility of real-time MRI-guided proton therapy and delivering a clinically acceptable dose to patients with various tumor locations within magnetic fields of up to 1.5 T. This work could serve as a guide and encouragement for further efforts toward clinical implementation of hybrid MRI–proton gantry systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Investigation of charge transfer and ionization in He-like systems (Li{sup +}, Be{sup 2+}, B{sup 3+}, C{sup 4+}, N{sup 5+}, O{sup 6+})-hydrogen atom collisions in Debye plasmas

The charge transfer and ionization cross sections have been calculated for He-like system (Li{sup +}, Be{sup 2+}, B{sup 3+}, C{sup 4+}, N{sup 5+}, O{sup 6+}) and hydrogen atom collisions in Debye plasmas for energies ranging from 1 to 500 keV using the classical trajectory Monte Carlo method. Interactions of the active electrons with the projectile ions have been described by model potentials. Cross sections are determined in both screening and unscreening environments, and a comparative study between both environments has been carried out. In particular, an interesting feature of sudden increase in the ionization cross sections at lower velocities is also observed in all the systems like pandey et al.[M. K. Pandey et al., Phys. Plasmas 19, 062104 (2012)] calculations for O{sup 8+}+ H collision. The feature of sudden increase in ionization cross sections at lower velocities and the effect of plasmas condition on it are explained in terms of the classical trajectory framework. We have found the cross sections for both capture and ionization are dependence on Debye screening lengths throughout the collision energies range, but are particularly pronounced at low projectile collisions energies. The calculated cross sections for the unscreened case are found in reasonable agreement with available experimental and theoretical results.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Composite Boson Description of a Low-Density Gas of Excitons

Ground-state properties of a fermionic Coulomb gas are calculated using the fixed-node diffusion Monte Carlo method. The validity of the composite boson description is tested for different densities. We extract the exciton–exciton s-wave scattering length by solving the four-body problem in a harmonic trap and mapping the energy to that of two trapped bosons. The equation of state is consistent with the Bogoliubov theory for composite bosons interacting with the obtained s-wave scattering length. The perturbative expansion at low density has contributions physically coming from (a) exciton binding energy, (b) mean-field Gross–Pitaevskii interaction between excitons, and (c) quantum depletion of the excitonic condensate (Lee–Huang–Yang terms for composite bosons). In addition, for low densities we find a good agreement with the Bogoliubov bosonic theory for the condensate fraction of excitons. The equation of state in the opposite limit of large density is found to be well described by the perturbative theory including (a) mixture of two ideal Fermi gases and (b) exchange energy. We find that for low densities both energetic and coherent properties are correctly described by the picture of composite bosons (excitons).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗