SEARCH · Search NASA
Results for “Constant Number Monte Carlo”
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.
Nanoparticle dynamics in the spatial afterglows of nonthermal plasma synthesis reactors
Nonthermal plasma flow tube reactors are industrially scalable systems for the production of nanocrystal (NC) based materials and coatings. One key advantage of nonthermal plasma synthesis is the ability to both synthesize NCs and deposit films in a single reactor, as at the reactor outlet, NCs can be inertially deposited onto a target substrate. The size and morphology of deposited particles can substantially influence the film structure and function. Though NCs are typically near-spherical and monodispersed as-produced in plasma synthesis reactors, NC charge and growth dynamics can be altered substantially when NCs are sampled out of the plasma and through the spatial afterglow region, affecting deposition. Experiments have demonstrated changes of NC size and charge in the spatial afterglow; however, these dynamics remain unexplored and unexplained via theory and simulation. To address this, we developed a constant number Monte Carlo (CNMC) simulation model to examine the mechanisms of NC decharging and growth in the spatial afterglow of plasma flow tube reactors. Collisions between NC and plasma species, diffusive deposition, and electron desorption from NCs are incorporated in the CNMC simulation. The simulation results are specifically compared with previous experiments on Si NCs synthesized from a low pressure Ar-SiH 4 nonthermal plasma reactor. Furthermore, the experiment-model comparison shows that CNMC models can be implemented which accurately model NC size distribution evolution in a spatial afterglow. Simultaneously, results show that improved collision models, energetic species diffusion models, and electron desorption models will be necessary to accurately depict NC dynamics in spatial afterglows.
Particle charge distributions in the effluent of a flow-through atmospheric pressure low temperature plasma
Atmospheric pressure low-temperature plasmas are often utilized to perform particle synthesis, treatment, and removal. It is well-known that dust particles are highly negatively charged in these plasmas; however, little is known about dust particle charging behavior as particles leave the plasma volume and pass through the spatial afterglow region. In this work, monodisperse particles of various sizes and work functions were introduced into an atmospheric pressure radiofrequency capacitively coupled flow-through plasma. Dust particle electrical mobility distributions downstream of the flow-through plasma were measured utilizing a differential mobility analyzer in conjunction with a condensation particle counter at various gas flow velocities. Charge distributions were determined from the measured electrical mobility distributions. Experiments confirm that particles become less negatively charged, and even net-positively charged after leaving the plasma volume, with a distribution that follows a shifted Boltzmann charge distribution. Additionally, particle charge in the effluent of the flow-through plasma is negligibly dependent on work function but highly size and flow velocity dependent. Larger particles were shown to have a higher magnitude of charge under all studied conditions; however, particle polarity was switchable by varying gas flow velocity. The charging dynamics were simulated utilizing a constant number Monte Carlo model that accounts for electron temperature decay and the transition from ambipolar to free diffusion of electrons and ions in the spatial afterglow. Simulation results also suggest that, at the same flow velocity, larger particles obtain a greater magnitude of charge, negative or positive. Furthermore, the decrease in electron mobility and the difference between ion and electron convective loss rates create an ion-rich region in the plasma effluent that promotes ion–particle collisions and drives particle charge removal and even reversal of polarity. Larger particles more favorably collide with energetic species in these environments, which results in higher charge states.
Understanding the depletion of electrons in dusty plasmas at atmospheric pressure
The nucleation and growth of NPs in the gas phase using atmospheric-pressure plasma systems is an important approach to synthesizing novel dimensionally-controlled materials. In this paper, we investigated the effect of the NPs on a typical type of continuous-flow, substrate-free plasma at atmospheric pressure to understand their potential contribution to electron density changes. A tandem plasma system was set up consisting of one plasma reactor that produced carbonaceous NPs from mixtures of argon and hexane, and another identical plasma reactor where the as-grown particles were injected and non-intrusive electrical and optical measurements were performed. The electron densities obtained from conductivity measurements and a plasma fluid model were found to decrease in the presence of NPs. However, control experiments revealed that the main source of the electron depletion was residual vapor or reaction byproducts in the form of molecular species or nanoclusters and not the particles themselves. These results were validated by constant number Monte Carlo simulations which showed that at the experimentally-measured conditions, the NPs were not of sufficiently high enough concentration to reduce the electron density; however, if molecules or clusters are ionizable, they remain in sufficient concentration to deplete electron densities. Our study shows that at atmospheric pressure, because of their typically larger electron density values, particle-producing plasmas are distinct from those at low pressure, and nanoparticle formation does not have the same impact while molecular-scale species may be a more important consideration.
Collaborative Research: Understanding Nanoparticle-Plasma Interactions in Dusty Non-Thermal Plasmas by Nanoparticle Probes and Online Aerosol Characterization (Final Technical Report)
This report describes results from the University of Minnesota’s contribution to studies performed under the award DE-SC0018202 (principal Investigator: Christopher J. Hogan Jr.. This award is a collaborative research grant with support equally provided to Case Western Reserve University and to the University of Illinois (with principal investigator R. Mohan Sankaran). To date ten peer-reviewed published papers have resulted from the studies performed with University of Minnesota personnel contributing; this report describes the results of these studies as well as several unpublished studies, grouped into seven separate study types. The overarching goal of the carried out studies was to understand nanoparticle behavior and dynamics in non-thermal plasma synthesis reactors, specifically adapting aerosol measurement technology to provide online measurements of particle size, charge, and structure in plasma reactor systems.
Direct Monte Carlo Simulations of Chemical Equilibrium Composition of Molecular Fluid Mixtures Under Shock Conditions
A new Monte Carlo simulation method is implemented which holds the number and kind of atoms constant rather than the number and kind of molecules for a molecular fluid mixture. Chemical reactions are allowed as a correlated move of atoms such as an interchange of atoms between molecules. Equilibrium chemical composition is determined as an average over the simulation along with equation of state quantities. Results are presented for N 2 + O 2 → 2NO at 30GPa and 3000K. The equilibrium composition is shown to be sensitive to the potentials between unlike species.
Defect energetics of cubic hafnia from quantum Monte Carlo simulations
Cubic hafnia (HfO 2 ) is of great interest for a number of applications in electronics because of its high dielectric constant. However, common defects in such applications degrade the properties of hafina. We have investigated the electronic properties of oxygen vacancies and nitrogen substitution in cubic HfO 2 using first-principles calculations based on density functional theory (DFT) and many-body diffusion Monte Carlo (DMC) methods. We investigate five different charge defect states of oxygen vacancies, as well as substitutional N defects which can lead to local magnetic moments. Both DMC and DFT calculations shows that an oxygen vacancy induces strong lattice relaxations around the defect. Finally, we compare defect formation energies, charge and spin densities obtained from DMC with results obtained using DFT. While the obtained formation energies from DMC are 0.6–1.5 eV larger than those from GGA+ U , the agreement for the most important defects, neutral and positively charged oxygen vacancies, and nitrogen substitutional defect, under oxygen-poor conditions are in reasonably good agreement. Lastly, our work confirms that nitrogen can act to passivate cubic hafnia for applications in electronics.
Probability techniques for reliability analysis of composite materials
Traditional design approaches for composite materials have employed deterministic criteria for failure analysis. New approaches are required to predict the reliability of composite structures since strengths and stresses may be random variables. This report will examine and compare methods used to evaluate the reliability of composite laminae. The two types of methods that will be evaluated are fast probability integration (FPI) methods and Monte Carlo methods. In these methods, reliability is formulated as the probability that an explicit function of random variables is less than a given constant. Using failure criteria developed for composite materials, a function of design variables can be generated which defines a 'failure surface' in probability space. A number of methods are available to evaluate the integration over the probability space bounded by this surface; this integration delivers the required reliability. The methods which will be evaluated are: the first order, second moment FPI methods; second order, second moment FPI methods; the simple Monte Carlo; and an advanced Monte Carlo technique which utilizes importance sampling. The methods are compared for accuracy, efficiency, and for the conservativism of the reliability estimation. The methodology involved in determining the sensitivity of the reliability estimate to the design variables (strength distributions) and importance factors is also presented.
A quasi-static Monte Carlo algorithm for the simulation of sub-prompt critical transients
A pseudo-time-dependent Monte Carlo algorithm suitable for simulating thermal reactor transients with timescales on the order of seconds is presented. An initial experimental implementation has been written into the Monte Carlo Application ToolKit (MCATK) developed at Los Alamos National Laboratory. The algorithm is essentially a combination of the fully time-dependent Monte Carlo algorithm, including the explicit simulation of delayed neutron precursors, and a quasi-static algorithm that does not use the point-kinetics equations. The fully time-dependent calculation is performed for some number of inactive and active micro time steps. The dynamic reactor time constant and delayed neutron precursor production rates are calculated directly from the active cycles, rather than from point-kinetics parameters. These quantities are then used to advance the neutron population and delayed neutron precursors forward in time across a macro time step. As a result, the method is applied to the C5G7-TD benchmark operating in steady state and in transient mode.
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.
Combining translational and rotational seismic motions to invert local-scale seismic data for time-variable moment tensors: do rotational motions help for high-frequency seismic data produced by underground explosions?
SUMMARY We present an analysis of combining translational and rotational seismic data in an inversion for the time-variable source time functions corresponding to the components of the seismic moment tensor. We conduct a series of numerical experiments where the data are simulated by a combination of an underground explosion and a co-located double couple shear source and recorded on surface-mounted seismometers within 1–2 km of the source. The experiments are designed to mimic explosion seismology experiments, and thus the data are in the 1–10 Hz frequency range and contain very few surface waves. We use a Monte Carlo method to propagate Earth model uncertainty into the estimates of seismic source parameters. In our experiments, we find that the uncertainty of the estimated seismic source parameters increases when we add rotational seismic motions to the inversion when using a constant number of data channels. In this case, the increased degree of uncertainty in the final results is most likely due to the near-surface Earth model uncertainty that we introduce in our simulations. However, for a fixed number of seismic stations, adding rotational seismic motions to the inversion acts to decrease the uncertainty of the estimated seismic source parameters, most likely due to the increase in the number of data channels used in the inversion.
Computational Materials Research
Computational Materials aims to model and predict thermodynamic, mechanical, and transport properties of polymer matrix composites. This workshop, the second coordinated by NASA Langley, reports progress in measurements and modeling at a number of length scales: atomic, molecular, nano, and continuum. Assembled here are presentations on quantum calculations for force field development, molecular mechanics of interfaces, molecular weight effects on mechanical properties, molecular dynamics applied to poling of polymers for electrets, Monte Carlo simulation of aromatic thermoplastics, thermal pressure coefficients of liquids, ultrasonic elastic constants, group additivity predictions, bulk constitutive models, and viscoplasticity characterization.
Charge conservation and higher moments of charge fluctuations
Higher moments of distributions of net charge and baryon number in heavy-ion collisions have been proposed as signals of fundamental QCD phase transitions. In order to better understand background processes for these observables, models are presented which enable one to gauge the effects of local charge conservation, decays of resonances and clusters, Bose symmetrization, and volume fluctuations. Monte Carlo methods for generating samplings of particles consistent with local charge conservation are presented and are followed by a review of simple analytic models involving a single type of charge with a constant experimental efficiency. The main model consists of thermal emission superimposed onto a simple parametrization of collective flow, known as a blast wave, with emission being consistent with individual canonical ensembles. The spatial extent of local charge conservation is parameterized by the size and extent over which charge is conserved. Here, the sensitivity of third- and fourth-order moments, skewness and kurtosis, to these parameters, and to beam energy and baryon density is explored. Comparisons with STAR data show that a significant part of the observed non-Poissonian fluctuations in net-proton fluctuations are explained by charge and baryon-number conservation, but that measurements of the STAR collaboration for fluctuations of net electric charge significantly differ from expectations of the models presented here.
Measuring signatures in photon angular spectra to distinguish nonlinear Compton scattering models
The collision of a high-energy electron beam with a laser pulse may be used to study radiation reaction and nonlinear Compton scattering among many other processes in strong-field quantum electrodynamics. Predictions from simulation and theory for these interactions rely on a number of approximations and assumptions that have not been experimentally tested. Here, experimentally measurable signatures are identified that might be able to distinguish between radiation reaction models, i.e., classical or quantum, or between the local constant field and local monochromatic approximations used to calculate the properties of the nonlinear Compton process. These signatures are considered through Monte Carlo simulations of various experimental conditions that are relevant to today's laser facilities. Potential detection schemes for measuring the signatures are proposed. We find that single-photon counting of keV photons to resolve harmonics and scintillator-based detection of MeV photons may allow us to validate nonlinear Compton scattering models and radiation reaction models respectively. This will require electron beams with divergence angles less than 2 mrad and less than 20% energy spread.
Geant4 Monte-Carlo (GEMC) A database-driven simulation program
GEMC[1] is an application that harnesses the power of databases to execute Geant4 Monte-Carlo simulations. The databases (MYSQL, CSQL, TEXT) define the geometry, materials, digitization algorithms, readout electronics and output formats. Implemented in C++, GEMC also boasts a user-friendly Python API that facilitates detector construction and database population. GEMC can handle real-life scenarios such as geometry variations and the run number-dependent calibration constants and digitization parameters. This abstract provides an overview of GEMC, accompanied by examples that showcase its versatility. We delve into the practical application of GEMC within the the CLAS12 experimental program at Jefferson Lab.
Evaluation of the spatial self-shielding impact for TRISO-based nuclear fuel depletion
Reactor physics analyses of nuclear cores with nuclear fuel concepts containing tristructural isotropic (TRISO) particles, such as pebbles or compact fuel elements, rely on various degrees of simplification to keep these highly heterogeneous problems computationally tractable. One such limitation regards the level of spatial discretization employed during burnup calculations, where traditionally only a limited number of spatial zones are modeled at the full core level and assume that the spectrum is constant within the fuel elements and TRISO particles in this depletion zone. This type of assumption neglects the impact of spatial self-shielding effect within the kernels (microscale level) as well as within the compact or pebbles (mesoscale level). Furthermore, the Monte Carlo code Serpent 2 contains many relevant features for efficiently modeling this type of geometry, including a collision-based domain decomposition intended for very large burnup calculations, which we leveraged for this work to quantify the impact of capturing neutron flux variations occurring at the micro- and mesoscale level on a series of high-temperature gas-cooled reactor fuel element depletion problems. While spatial self-shielding is observed at both scales, with differences from a volume-averaged burnup of ±7% within the kernels and ±2% between TRISO particles within the fuel element, the conjugated effect on nuclide inventories and multiplication factor are negligible, hence confirming that assuming a single average spectrum value may be sufficient for most applications.
Dissipative ground state preparation in ab initio electronic structure theory
Dissipative engineering is a powerful tool for quantum state preparation, and has drawn significant attention in quantum algorithms and quantum many-body physics in recent years. In this work, we introduce a novel approach using the Lindblad dynamics to efficiently prepare the ground state for general ab initio electronic structure problems on quantum computers, without variational parameters. These problems often involve Hamiltonians that lack geometric locality or sparsity structures, which we address by proposing two generic types of jump operators for the Lindblad dynamics. Type-I jump operators break the particle number symmetry and should be simulated in the Fock space. Type-II jump operators preserves the particle number symmetry and can be simulated more efficiently in the full configuration interaction space. For both types of jump operators, we prove that in a simplified Hartree-Fock framework, the spectral gap of our Lindbladian is lower bounded by a universal constant. For physical observables such as energy and reduced density matrices, the convergence rate of our Lindblad dynamics with Type-I jump operators remains universal, while the convergence rate with Type-II jump operators only depends on coarse grained information such as the number of orbitals and the number of electrons. To validate our approach, we employ a Monte Carlo trajectory-based algorithm for simulating the Lindblad dynamics for full ab initio Hamiltonians, demonstrating its effectiveness on molecular systems amenable to exact wavefunction treatment.
(U) Correlated Sampling Using Batch Statistics to Reduce the Uncertainty of Combinations of KSEN Outputs with MCNP6
The relative sensitivity of k eff to the densities of nuclides in a material are combined to compute relative sensitivities to other inputs. When computed in a single Monte Carlo run, the nuclide density sensitivities are correlated, and the statistical uncertainties propagated to other inputs will be incorrect unless those correlations are accounted for. Equations are presented to apply correlated sampling using batch statistics for the sum of an arbitrary number of random tallies and the difference of two random tallies when each is multiplied by a different constant. When correlated sampling is used on a recent benchmark evaluation, the correct statistical uncertainties for certain combinations of sensitivities are dramatically smaller than the incorrect uncertainties. A user-controlled, regular output of MCNP6’s KSEN sensitivities that allows batch statistics to be applied to combinations would be extremely valuable.