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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

On the uncertainty of temperature estimation in a rapid compression machine

Rapid compression machines (RCMs) have been widely used in the combustion literature to study the low-to-intermediate temperature ignition of many fuels. In a typical RCM, the pressure during and after the compression stroke is measured. Yet, measurement of the temperature history in the RCM reaction chamber is challenging. Thus, the temperature is generally calculated by the isentropic relations between pressure and temperature, assuming that the adiabatic core hypothesis holds. To estimate the uncertainty in the calculated temperature, an uncertainty propagation analysis must be carried out. Our previous analyses assumed that the uncertainties of the parameters in the equation to calculate the temperature were normally distributed and independent, but these assumptions do not hold for typical RCM operating procedures. Here, a Monte Carlo method is developed to estimate the uncertainty in the calculated temperature, while taking into account the correlation between parameters and the possibility of non-normal probability distributions. In addition, the Monte Carlo method is compared to an analysis that assumes normally distributed, independent parameters. Both analysis methods show that the magnitude of the initial pressure and the uncertainty of the initial temperature have strong influences on the magnitude of the uncertainty. Finally, the uncertainty estimation methods studied here provide a reference value for the uncertainty of the reference temperature in an RCM and can be generalized to other similar facilities.

42 ENGINEERING↗
DOE OSTI2021

Origin of the $ν$ = 1/2 fractional quantum Hall effect in wide quantum wells

The nature of the fractional quantum Hall effect at ν = 1/2, observed in wide quantum wells al- most three decades ago, is still under debate. Previous studies have investigated it by the variational Monte Carlo method, which assumes that the transverse wave function and the gap between the symmetric and antisymmetric subbands obtained in a local density approximation at zero magnetic field remain valid even at high perpendicular magnetic fields; this method also ignores the effect of Landau level mixing. We develop in this work a three-dimensional fixed phase diffusion Monte Carlo method, which gives, in a single framework, the total energies of various candidate states in a finite width quantum well, including Landau level mixing, directly in a large magnetic field. This method can be applied to one-component states, and also to two-component states in the limit where the symmetric and antisymmetric bands are nearly degenerate. Our three-dimensional fixed-phase diffusion Monte Carlo calculations find that the one-component composite-fermion Fermi sea and the one-component Pfaffian states are very close in energy for a range of quantum well widths and densities, suggesting that the observed 1/2 fractional quantum Hall state in wide quantum wells is likely to be the one-component Pfaffian state. Furthermore, we hope that this will motivate further experimental studies of this state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

SU-F-T-620: Development of a Convolution/Superposition Dose Engine for CyberKnife System

Purpose: Current CyberKnife treatment planning system (TPS) provided two dose calculation algorithms: Ray-tracing and Monte Carlo. Ray-tracing algorithm is fast, but less accurate, and also can’t handle irregular fields since a multi-leaf collimator system was recently introduced to CyberKnife M6 system. Monte Carlo method has well-known accuracy, but the current version still takes a long time to finish dose calculations. The purpose of this paper is to develop a GPU-based fast C/S dose engine for CyberKnife system to achieve both accuracy and efficiency. Methods: The TERMA distribution from a poly-energetic source was calculated based on beam’s eye view coordinate system, which is GPU friendly and has linear complexity. The dose distribution was then computed by inversely collecting the energy depositions from all TERMA points along 192 collapsed-cone directions. EGSnrc user code was used to pre-calculate energy deposition kernels (EDKs) for a series of mono-energy photons The energy spectrum was reconstructed based on measured tissue maximum ratio (TMR) curve, the TERMA averaged cumulative kernels was then calculated. Beam hardening parameters and intensity profiles were optimized based on measurement data from CyberKnife system. Results: The difference between measured and calculated TMR are less than 1% for all collimators except in the build-up regions. The calculated profiles also showed good agreements with the measured doses within 1% except in the penumbra regions. The developed C/S dose engine was also used to evaluate four clinical CyberKnife treatment plans, the results showed a better dose calculation accuracy than Ray-tracing algorithm compared with Monte Carlo method for heterogeneous cases. For the dose calculation time, it takes about several seconds for one beam depends on collimator size and dose calculation grids. Conclusion: A GPU-based C/S dose engine has been developed for CyberKnife system, which was proven to be efficient and accurate for clinical purpose, and can be easily implemented in TPS.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Estimating Large Global Significances with a New Monte Carlo Extrapolation Method

In particle physics, it is needed to evaluate the possibility that excesses of events in mass spectra are due to statistical fluctuations as quantified by the standards of local and global significances. Without prior knowledge of a particle’s mass, it is especially critical to estimate its global significance. The usual approach is to count the number of times a significance limit is exceeded in a collection of simulated Monte Carlo (MC) “toy experiments”. To demonstrate this conventional method for global significance, we performed simulation studies according to a recent Compact Muon Solenoid (CMS) result to show its effectiveness. However, this counting method is not practical for computing large global significances. To address this problem, we developed a new “extrapolation” method to evaluate the global significance. We compared the global significance estimated by our new method with that of the conventional approach, and verified its feasibility and effectiveness. This method is also applicable for cases where only small toy MC samples are available. In this approach, the significance is calculated based on p-values, assuming symmetrical Gaussian distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1994

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗
DOE OSTI2015

How large are nonadiabatic effects in atomic and diatomic systems?

With recent developments in simulating nonadiabatic systems to high accuracy, it has become possible to determine how much energy is attributed to nuclear quantum effects beyond zero-point energy. Here, we calculate the non-relativistic ground-state energies of atomic and molecular systems without the Born-Oppenheimer approximation. For this purpose, we utilize the fixed-node diffusion Monte Carlo method, in which the nodes depend on both the electronic and ionic positions. Our report shows the ground-state energies for all systems studied, ionization energies for the first-row atoms and atomization energies for the first-row hydrides. We find the ionization energies of the atoms to be nearly independent of the Born-Oppenheimer approximation, within the accuracy of our results. The atomization energies of molecular systems, however, show small effects of the nonadiabatic coupling between electrons and nuclei.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Correlated rattling-ion origins of dielectric properties in reentrant dipole glasses BaTiO{sub 3}-BiScO{sub 3}

The local structure of the pseudo-cubic solid solution 0.6BaTiO{sub 3}-0.4BiScO{sub 3}, which exhibits reentrant dipole-glass behavior, has been determined using the Reverse Monte Carlo method to simultaneously fit (1) neutron and X-ray total scattering data (including the corresponding real-space pair-distribution functions), (2) Bi and Sc extended X-ray absorption fine structure, and (3) patterns of diffuse scattering in electron diffraction. These structural refinements revealed the multi-site probability density distributions for both Bi (14-sites) and Ti (8 sites), whereas Ba and Sc featured normal unimodal distributions. Bi atoms are displaced along both the 〈111〉 and 〈100〉 directions, while Ti atoms are shifted along 〈111〉. Correlated dynamic hopping of Bi and Ti over their corresponding split sites combined with chemical disorder is proposed as the origin of the strong frequency dispersion observed in dielectric measurements. The existence of split sites also explains the reentrant dipole-glass behavior reported for this system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Impact of system parameters and geospatial variables on the reliability of residential systems with PV and energy storage

A reliable power supply is the foundation of modern society, enabling technologies used to function within a society. Residential systems are places where the end users directly consume power, enabling technologies to sustain life. With the emergence of behind-the-meter resources, the end-users have some control over power supply reliability. The intermittency and variability of these resources impact residential system reliability. In this work, we study the reliability performance of a grid-supplemented residential system with behind-the-meter Distributed Energy Resources(DER) subject to various system parameters and geospatial variables. We propose a multistate reliability model for the behind-the-meter microinverter-based Photo Voltaic (PV) system and integrated inverter-based energy storage (ES) system. A sequential Monte Carlo method is then presented to evaluate the reliability indices for the residential system with behind-the-meter DERs as the main supply and the grid as the backup. The sequential Monte Carlo method is used to analyze the reliability performance of the residential model at the top 100 populous counties of the United States, where actual load and solar irradiance data at the counties is used. In the analysis, the sensitivity of indices to system parameters such as DER size and the sensitivity of the indices to climate zone and gross horizontal irradiation (GHI), which affects the load conditions and the PV output of the residential system, is performed. The analysis results show that the system parameters and the geospatial variables significantly impact the residential system’s reliability. The insights from this analysis will be of immense value to the distribution system planners to provide zone-specific guidelines for DER system sizing and toward the evolution of utility business models. Finally, the methodology developed can be used to extend the analysis to other locations.

14 SOLAR ENERGY↗
DOE OSTI2015

Development and application of a standardized flow measurement uncertainty analysis framework to various low-head short-converging intake types across the United States federal hydropower fleet

Turbine discharges at low-head short converging intakes are difficult to measure accurately. The proximity of the measurement section to the intake entrance admits large uncertainties related to asymmetry of the velocity profile, swirl, and turbulence. Existing turbine performance codes [10, 24] do not address this special case and published literature is largely silent on rigorous evaluation of uncertainties associated with this measurement context. The American Society of Mechanical Engineers (ASME) Committee investigated the use of Acoustic transit time (ATT), Acoustic scintillation (AS), and Current meter (CM) in a short converging intake at the Kootenay Canal Generating Station in 2009. Based on their findings, a standardized uncertainty analysis (UA) framework for velocity-area method (specifically for CM measurements) is presented in this paper given the fact that CM is still the most fundamental and common type of measurement system. Typical sources of systematic and random errors associated with CM measurements are investigated, and the major sources of uncertainties associated with turbulence and velocity fluctuations, numerical velocity integration technique (bi-cubic spline), and the number and placement of current meters are being considered for an evaluation. Since the velocity measurements in a short converging intake are associated with complex nonlinear and time varying uncertainties (e.g., Reynolds stress in fluid dynamics), simply applying the law of propagation of uncertainty is known to overestimate the measurement variance while the Monte Carlo method does not. Therefore, a pseudo-Monte Carlo simulation method (random flow generation technique [8]) which was initially developed for the purpose of establishing upstream or initial conditions in the Large-Eddy Simulation (LES) and the Direct Numerical Simulation (DNS) is used to statistically determine uncertainties associated with turbulence and velocity fluctuations. This technique is then combined with a bi-cubic spline interpolation method which converts point velocities into a continuous velocity distribution over the measurement domain. Subsequently the number and placement of current meters are simulated to investigate the accuracy of the estimated flow rates using the numerical velocity-area integration method outlined in ISO 3354 [12]. The authors herein consider that statistics on generated flow rates processed with bi-cubic interpolation and sensor simulations are the combined uncertainties which already accounted for the effects of all those three uncertainty sources. A preliminary analysis based on the current meter data obtained through an upgrade acceptance test of a single unit located in a mainstem plant has been presented.

Smith, Brennan T↗
DOE OSTI2021

Stochastically estimated covariance matrices for independent and cumulative fission yields in the ENDF/B-VIII.0 and JEFF-3.3 evaluations

A Monte-Carlo method for the generation of correlation and covariance matrices for independent and cumulative fission yields has been developed. The method uses a constrained Monte-Carlo resampling structure in order to vary evaluated fission yield libraries in a way that meets basic conservation principles. This results in the generation of correlation/covariance matrices with limited model bias and uncertainty; the matrices are primarily reflective of the evaluated fission yield uncertainties and correlations that arise from the evaluation process. This method has been applied to generate correlation and covariance matrices for all of the fissioning systems of the ENDF/B-VIII.0 and JEFF-3.3 evaluations, marking the first time such matrices have been generated for all of these systems. These covariance matrices have been published online for immediate public use. These correlation and covariance matrices can be used to improve uncertainty estimation in calculations of reactor antineutrino emission rates, decay heat problems, and nuclear forensics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Generalizing the self-healing diffusion Monte Carlo approach to finite temperature: A path for the optimization of low-energy many-body bases

A statistical method is derived for the calculation of thermodynamic properties of many-body systems at low temperatures. This method is based on the self-healing diffusion Monte Carlo method for complex functions [F. A. Reboredo, J. Chem. Phys. 136, 204101 (2012)] and some ideas of the correlation function Monte Carlo approach [D. M. Ceperley and B. Bernu, J. Chem. Phys. 89, 6316 (1988)]. In order to allow the evolution in imaginary time to describe the density matrix, we remove the fixed-node restriction using complex antisymmetric guiding wave functions. In the process we obtain a parallel algorithm that optimizes a small subspace of the many-body Hilbert space to provide maximum overlap with the subspace spanned by the lowest-energy eigenstates of a many-body Hamiltonian. We show in a model system that the partition function is progressively maximized within this subspace. We show that the subspace spanned by the small basis systematically converges towards the subspace spanned by the lowest energy eigenstates. Possible applications of this method for calculating the thermodynamic properties of many-body systems near the ground state are discussed. The resulting basis can also be used to accelerate the calculation of the ground or excited states with quantum Monte Carlo.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

New method for level-lifetime measurements with thick scintillators

Level lifetimes provide key insight into the structure of atomic nuclei and serve as stringent tests of theoretical descriptions. Though several methods for determining level lifetimes exist for both reaction measurements and decay studies, here the focus is on techniques involving the direct measurement of time differences between population and subsequent depopulation of excited states. The techniques presented herein are broadly applicable across multiple timing ranges, but the approach is specifically described for the βγ timing method. A multi-step, amplitude-dependent time walk correction was employed to address the data analysis complications that arise from using thick scintillators for electron detection. Further, additional corrections for the depth of the interaction must also be performed when the parent isotope has been implanted into the detector. Techniques are presented for performing these time walk and depth of interaction corrections. Subsequently, a new Monte Carlo method utilizing measured detector responses obtained from the data, coupled with chi-square minimization, is presented for extracting excited state lifetimes ≳ 100 picoseconds. The framework of this Monte Carlo method is developed for the decay of a state in 68 Zn with a known 1.6 ps half life, which is considered prompt given the detection sensitivity of the technique, and then benchmarked using two other excited states in neutron-rich Ni isotopes with 120(34) ps and 1.05(3) ns half lives. Using this new method which takes into account the thick scintillator used, these same half lives were measured to be 135(10) ps and 1.04(24) ns, respectively. The overall good agreement demonstrates the validity of the technique.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SCoPE: an efficient method of Cosmological Parameter Estimation

Markov Chain Monte Carlo (MCMC) sampler is widely used for cosmological parameter estimation from CMB and other data. However, due to the intrinsic serial nature of the MCMC sampler, convergence is often very slow. Here we present a fast and independently written Monte Carlo method for cosmological parameter estimation named as Slick Cosmological Parameter Estimator (SCoPE), that employs delayed rejection to increase the acceptance rate of a chain, and pre-fetching that helps an individual chain to run on parallel CPUs. An inter-chain covariance update is also incorporated to prevent clustering of the chains allowing faster and better mixing of the chains. We use an adaptive method for covariance calculation to calculate and update the covariance automatically as the chains progress. Our analysis shows that the acceptance probability of each step in SCoPE is more than 95% and the convergence of the chains are faster. Using SCoPE, we carry out some cosmological parameter estimations with different cosmological models using WMAP-9 and Planck results. One of the current research interests in cosmology is quantifying the nature of dark energy. We analyze the cosmological parameters from two illustrative commonly used parameterisations of dark energy models. We also asses primordial helium fraction in the universe can be constrained by the present CMB data from WMAP-9 and Planck. The results from our MCMC analysis on the one hand helps us to understand the workability of the SCoPE better, on the other hand it provides a completely independent estimation of cosmological parameters from WMAP-9 and Planck data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Charge order–superfluidity transition in a two-dimensional system of hard-core bosons and emerging domain structures

We have used high-performance parallel computations by NVIDIA graphics cards applying the method of nonlinear conjugate gradients and Monte Carlo method to observe directly the developing ground state configuration of a two-dimensional hard-core boson system with decrease in temperature, and its evolution with deviation from a half-filling. This has allowed us to explore unconventional features of a charge order—superfluidity phase transition, specifically, formation of an irregular domain structure, emergence of a filamentary superfluid structure that condenses within of the charge-ordered phase domain antiphase boundaries, and formation and evolution of various topological structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Integral Criticality Estimators in MCATK

The Monte Carlo Application ToolKit (MCATK) is a component-based software toolset for delivering customized particle transport solutions using the Monte Carlo method. Currently under development in the XCP Monte Carlo group at Los Alamos National Laboratory, the toolkit has the ability to estimate the ke f f and a eigenvalues for static geometries. This paper presents a description of the estimators and variance reduction techniques available in the toolkit and includes a preview of those slated for future releases. Along with the description of the underlying algorithms is a description of the available user inputs for controlling the iterations. The paper concludes with a comparison of the MCATK results with those provided by analytic solutions. The results match within expected statistical uncertainties and demonstrate MCATK’s usefulness in estimating these important quantities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

On the Sensitivity of α/β Prediction to Dose Calculation Methodology in Prostate Brachytherapy

Purpose: To study the relationship between the accuracy of the dose calculation in brachytherapy and the estimations of the radiosensitivity parameter, α/β, for prostate cancer. Methods and Materials: In this study, Monte Carlo methods and more specifically the code ALGEBRA was used to produce accurate dose calculations in the case of prostate brachytherapy. Equivalent uniform biologically effective dose was calculated for these dose distributions and was used in an iso-effectiveness relationship with external beam radiation therapy. Results: By considering different levels of detail in the calculations, the estimation for the α/β parameter varied from 1.9 to 6.3 Gy, compared with a value of 3.0 Gy suggested by the American Association of Physicists in Medicine Task Group 137. Conclusions: Large variations of the α/β show the sensitivity of this parameter to dose calculation modality. The use of accurate dose calculation engines is critical for better evaluating the biological outcomes of treatments.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Temperature-extrapolation method for Implicit Monte Carlo - Radiation hydrodynamics calculations

We present a method for implementing temperature extrapolation in Implicit Monte Carlo solutions to radiation hydrodynamics problems. The method is based on a BDF-2 type integration to estimate a change in material temperature over a time step. We present results for radiation only problems in an infinite medium and for a 2-D Cartesian hohlraum problem. Additionally, radiation hydrodynamics simulations are presented for an RZ hohlraum problem and a related 3D problem. Our results indicate that improvements in noise and general behavior are possible. We present considerations for future investigations and implementations. (authors)

97 MATHEMATICS AND COMPUTING↗