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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Testing coupled dark energy with large scale structure observation

The coupling between the dark components provides a new approach to mitigate the coincidence problem of cosmological standard model. In this paper, dark energy is treated as a fluid with a constant equation of state, whose coupling with dark matter is Q-bar =3Hξ{sub x}ρ-bar {sub x}. In the frame of dark energy, we derive the evolution equations for the density and velocity perturbations. According to the Markov Chain Monte Carlo method, we constrain the model by currently available cosmic observations which include cosmic microwave background radiation, baryon acoustic oscillation, type Ia supernovae, and fσ{sub 8}(z) data points from redshift-space distortion. The results show the interaction rate in σ regions: ξ{sub x} = 0.00328{sub -0.00328-0.00328-0.00328}{sup +0.000736+0.00549+0.00816}, which means that the recently cosmic observations favor a small interaction rate which is up to the order of 10{sup -2}, meanwhile, the measurement of redshift-space distortion could rule out the large interaction rate in the σ region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Adaptive hybrid simulations for multiscale stochastic reaction networks

The probability distribution describing the state of a Stochastic Reaction Network (SRN) evolves according to the Chemical Master Equation (CME). It is common to estimate its solution using Monte Carlo methods such as the Stochastic Simulation Algorithm (SSA). In many cases, these simulations can take an impractical amount of computational time. Therefore, many methods have been developed that approximate sample paths of the underlying stochastic process and estimate the solution of the CME. A prominent class of these methods include hybrid methods that partition the set of species and the set of reactions into discrete and continuous subsets. Such a partition separates the dynamics into a discrete and a continuous part. Simulating such a stochastic process can be computationally much easier than simulating the exact discrete stochastic process with SSA. Moreover, the quasi-stationary assumption to approximate the dynamics of fast subnetworks can be applied for certain classes of networks. However, as the dynamics of a SRN evolves, these partitions may have to be adapted during the simulation. We develop a hybrid method that approximates the solution of a CME by automatically partitioning the reactions and species sets into discrete and continuous components and applying the quasi-stationary assumption on identifiable fast subnetworks. Our method does not require any user intervention and it adapts to exploit the changing timescale separation between reactions and/or changing magnitudes of copy-numbers of constituent species. We demonstrate the efficiency of the proposed method by considering examples from systems biology and showing that very good approximations to the exact probability distributions can be achieved in significantly less computational time. This is especially the case for systems with oscillatory dynamics, where the system dynamics change considerably throughout the time-period of interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Modeling the spectral modification of lower hybrid wave in the presence of drift-wave type density fluctuation in the scrape-off-layer of the EAST tokamak

In this work, the spectrum change of the lower hybrid (LH) waves caused by a low-frequency density fluctuation in the scrape-off-layer is studied by applying the wave scattering model developed by Bonoli and Ott [Phys. Fluids 25, 359 (1982)] via a Monte Carlo method. A ray-tracing model specific to this purpose is developed to evaluate the probability distribution of both the poloidal refractive index (N θ ) and the parallel refractive index (N ∥ ) of the LH wave at the last closed flux surface (LCFS).The probability distributions of N ∥ and N θ at the LCFS are studied using the EAST parameters as a function of wave frequency, the initial N ∥ , and the polar injection position, which may influence the lower hybrid current drive efficiency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Effects of magnetic drift tangential to magnetic surfaces on neoclassical transport in non-axisymmetric plasmas

In evaluating neoclassical transport by radially local simulations, the magnetic drift tangential to a flux surface is usually ignored in order to keep the phase-space volume conservation. In this paper, effect of the tangential magnetic drift on the local neoclassical transport is investigated. To retain the effect of the tangential magnetic drift in the local treatment of neoclassical transport, a new local formulation for the drift kinetic simulation is developed. The compressibility of the phase-space volume caused by the tangential magnetic drift is regarded as a source term for the drift kinetic equation, which is solved by using a two-weight δf Monte Carlo method for non-Hamiltonian system [G. Hu and J. A. Krommes, Phys. Plasmas 1, 863 (1994)]. It is demonstrated that the effect of the drift is negligible for the neoclassical transport in tokamaks. In non-axisymmetric systems, however, the tangential magnetic drift substantially changes the dependence of the neoclassical transport on the radial electric field E{sub r}. The peaked behavior of the neoclassical radial fluxes around E{sub r }={sub }0 observed in conventional local neoclassical transport simulations is removed by taking the tangential magnetic drift into account.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Particle Acceleration in Two Converging Shocks

Observations by spacecraft such as ACE , STEREO , and others show that there are proton spectral “breaks” with energy E {sub br} at 1–10 MeV in some large CME-driven shocks. Generally, a single shock with the diffusive acceleration mechanism would not predict the “broken” energy spectrum. The present paper focuses on two converging shocks to identify this energy spectral feature. In this case, the converging shocks comprise one forward CME-driven shock on 2006 December 13 and another backward Earth bow shock. We simulate the detailed particle acceleration processes in the region of the converging shocks using the Monte Carlo method. As a result, we not only obtain an extended energy spectrum with an energy “tail” up to a few 10 MeV higher than that in previous single shock model, but also we find an energy spectral “break” occurring on ∼5.5 MeV. The predicted energy spectral shape is consistent with observations from multiple spacecraft. The spectral “break,” then, in this case is caused by the interaction between the CME shock and Earth’s bow shock, and otherwise would not be present if Earth were not in the path of the CME.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Effect of size on electronic states in a strained pyramidal InAs-GaAs quantum dot system

The effect of size on electronic states in a strained pyramidal InAs-GaAs quantum dot system was studied. A comparison was made between two InAs quantum pyramids of different sizes embedded inside a cubic GaAs susbtrate material. Strain relaxation was carried out via the Metropolis Monte Carlo method and the calculated local strain tensors were then included to solve the energy values and the wave functions of the electronic states inside the two simulation cube. The 3D finite difference scheme was employed to solve the time independent Schrödinger equation based on the decoupled electron-hole model. Calculated energy values of the four lowest electronic states showed that the transitions between the electron and hole states widen as the size of the dot becomes smaller especially between the ground states. The confinement of electrons and holes become weaker as the size of the dot reduces.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

New capabilities of the MORET 6 Monte Carlo neutron transport code

The MORET code is a simulation tool that solves the transport equation for neutrons using the Monte Carlo method. It allows users to model complex three-dimensional geometrical configurations in a user-friendly way. New features have been introduced to extend the application field of MORET beyond the usual criticality calculations for which it has been initially designed. The most important change is the addition of an analog fixed source mode which allows studies of systems of any reactivity combined with very flexible outputs. Other useful improvements have been added concerning the geometric part, the fission matrix, the multigroup sensitivity coefficients and the outputs. This paper presents an overview of these new features. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

A partially reflecting random walk on spheres algorithm for electrical impedance tomography

In this work, we develop a probabilistic estimator for the voltage-to-current map arising in electrical impedance tomography. This novel so-called partially reflecting random walk on spheres estimator enables Monte Carlo methods to compute the voltage-to-current map in an embarrassingly parallel manner, which is an important issue with regard to the corresponding inverse problem. Our method uses the well-known random walk on spheres algorithm inside subdomains where the diffusion coefficient is constant and employs replacement techniques motivated by finite difference discretization to deal with both mixed boundary conditions and interface transmission conditions. We analyze the global bias and the variance of the new estimator both theoretically and experimentally. Subsequently, the variance of the new estimator is considerably reduced via a novel control variate conditional sampling technique which yields a highly efficient hybrid forward solver coupling probabilistic and deterministic algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Optical fields in porous polylactide matrices

Optical fields induced by laser light inside porous polymer matrices (scaffolds) under stimulating effect on cells are investigated. The spectra of light scattering coefficients and anisotropy factors of porous polylactide scaffolds in the visible and near-IR ranges are determined by numerical simulation based on experimental data obtained using a double integrating sphere. Studies are conducted for dry and wet scaffolds fabricated by supercritical fluid foaming and surface-selective laser sintering. The Monte Carlo method is used to calculate the fluence rate distribution under the assumptions of the Henyey – Greenstein phase scattering function and the model phase function of scattering on spherical cavities. It is shown that th e fluence rate distribution mainly depends on the scattering coefficient and the mean cosine of the deflection angle (g-factor). The results obtained are of great practical importance for tissue engineering based on the use of polymeric porous scaffolds with cells to replace tissue defects; they can be used to determine the dose range for laser stimulation of cells. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Steady state in a gas of inelastic rough spheres heated by a uniform stochastic force

We study here the steady state attained in a granular gas of inelastic rough spheres that is subject to a spatially uniform random volume force. The stochastic force has the form of the so-called white noise and acts by adding impulse to the particle translational velocities. We work out an analytical solution of the corresponding velocity distribution function from a Sonine polynomial expansion that displays energy non-equipartition between the translational and rotational modes, translational and rotational kurtoses, and translational-rotational velocity correlations. By comparison with a numerical solution of the Boltzmann kinetic equation (by means of the direct simulation Monte Carlo method), we show that our analytical solution provides a good description that is quantitatively very accurate in certain ranges of inelasticity and roughness. We also find three important features that make the forced granular gas steady state very different from the homogeneous cooling state (attained by an unforced granular gas). First, the marginal velocity distributions are always close to a Maxwellian. Second, there is a continuous transition to the purely smooth limit (where the effects of particle rotations are ignored). And third, the angular translational-rotational velocity correlations show a preference for a quasiperpendicular mutual orientation (which is called “lifted-tennis-ball” behavior)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

The use of Compton scattering in detecting anomaly in soil-possible use in pyromaterial detection

The Compton scattering is able to determine the signature of land mine detection based on dependency of density anomaly and energy change of scattered photons. In this study, 4.43 MeV gamma of the Am-Be source was used to perform Compton scattering. Two detectors were placed between source with distance of 8 cm and radius of 1.9 cm. Detectors of thallium-doped sodium iodide NaI(TI) was used for detecting gamma ray. There are 9 anomalies used in this simulation. The physical of anomaly is in cylinder form with radius of 10 cm and 8.9 cm height. The anomaly is buried 5 cm deep in the bed soil measured 80 cm radius and 53.5 cm height. Monte Carlo methods indicated the scattering of photons is directly proportional to density of anomalies. The difference between detector response with anomaly and without anomaly namely contrast ratio values are in a linear relationship with density of anomalies. Anomalies of air, wood and water give positive contrast ratio values whereas explosive, sand, concrete, graphite, limestone and polyethylene give negative contrast ratio values. Overall, the contrast ratio values are greater than 2 % for all anomalies. The strong contrast ratios result a good detection capability and distinction between anomalies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

THE EFFECT OF MAGNETIC FIELDS AND AMBIPOLAR DIFFUSION ON CORE MASS FUNCTIONS

Linear analysis of the formation of protostellar cores in planar magnetic interstellar clouds yields information about length scales involved in star formation. Combining these length scales with various distributions of other environmental variables (i.e., column density and mass-to-flux ratio) and applying Monte Carlo methods allow us to produce synthetic core mass functions (CMFs) for different environmental conditions. Our analysis shows that the shape of the CMF is directly dependent on the physical conditions of the cloud. Specifically, magnetic fields act to broaden the mass function and develop a high-mass tail while ambipolar diffusion will truncate this high-mass tail. In addition, we analyze the effect of small number statistics on the shape and high-mass slope of the synthetic CMFs. We find that observed CMFs are severely statistically limited, which has a profound effect on the derived slope for the high-mass tail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

On the evaporation of solar dark matter: spin-independent effective operators

As a part of the effort to investigate the implications of dark matter (DM)-nucleon effective interactions on the solar DM detection, in this paper we focus on the evaporation of the solar DM for a set of the DM-nucleon spin-independent (SI) effective operators. In order to put the evaluation of the evaporation rate on a more reliable ground, we calculate the non-thermal distribution of the solar DM using the Monte Carlo methods, rather than adopting the Maxwellian approximation. We then specify relevant signal parameter spaces for the solar DM detection for various SI effective operators. Based on the analysis, we determine the minimum DM masses for which the DM-nucleon coupling strengths can be probed from the solar neutrino observations. As an interesting application, our investigation also shows that evaporation effect can not be neglectd in a recent proposal aiming to solve the solar abundance problem by invoking the momentum-dependent asymmetric DM in the Sun.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

COMPARISON OF DIVERSITY OF TYPE IIB SUPERNOVAE WITH ASYMMETRY IN CASSIOPEIA A USING LIGHT ECHOES

We compare the diversity of spectral line velocities in a large sample of type IIb supernovae (SNe IIb) with the expected asphericity in the explosion, as measured from the light echoes (LEs) of Cassiopeia A (Cas A), which was a historical galactic SN IIb. We revisit the results of Rest et al., who used LEs to observe Cas A from multiple lines of sight and hence determine its asphericity, as seen in the velocity of three spectral lines (He i λ 5876, H α , and the Ca ii near-infrared (NIR) triplet). We confirm and improve on this measurement by reproducing the effect of the LEs in the spectra of several extragalactic SNe IIb found in the literature as well as mean SN IIb spectra recently created by Liu et al. and comparing these to the observed light echo spectra of Cas A, including their associated uncertainties. In order to quantify the accuracy of this comparison, we smooth the light echo spectra of Cas A using Gaussian processes and use a Monte Carlo method to measure the absorption velocities of these three features in the spectra. We then test the hypothesis that the diversity of ejecta velocities seen in SNe IIb can be explained by asphericity. We do this by comparing the range of velocities seen in the different LEs, and hence different lines of sight, of Cas A to that seen in the population of SNe IIb. We conclude that these two ranges are of the same order and thus asphericity could be enough to explain the diversity in the expansion velocity alone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Numerical simulation of ion charge breeding in electron beam ion source

The Electron Beam Ion Source particle-in-cell code (EBIS-PIC) tracks ions in an EBIS electron beam while updating electric potential self-consistently and atomic processes by the Monte Carlo method. Recent improvements to the code are reported in this paper. The ionization module has been improved by using experimental ionization energies and shell effects. The acceptance of injected ions and the emittance of extracted ion beam are calculated by extending EBIS-PIC to the beam line transport region. An EBIS-PIC simulation is performed for a Cs charge-breeding experiment at BNL. The charge state distribution agrees well with experiments, and additional simulation results of radial profiles and velocity space distributions of the trapped ions are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Computational study of trimer self-assembly and fluid phase behavior

The fluid phase diagram of trimer particles composed of one central attractive bead and two repulsive beads was determined as a function of simple geometric parameters using flat-histogram Monte Carlo methods. A variety of self-assembled structures were obtained including spherical micelle-like clusters, elongated clusters, and densely packed cylinders, depending on both the state conditions and shape of the trimer. Advanced simulation techniques were employed to determine transitions between self-assembled structures and macroscopic phases using thermodynamic and structural definitions. Simple changes in particle geometry yield dramatic changes in phase behavior, ranging from macroscopic fluid phase separation to molecular-scale self-assembly. In special cases, both self-assembled, elongated clusters and bulk fluid phase separation occur simultaneously. Our work suggests that tuning particle shape and interactions can yield superstructures with controlled architecture.

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

SU-E-T-598: Parametric Equation for Quick and Reliable Estimate of Stray Neutron Doses in Proton Therapy and Application for Intracranial Tumor Treatments

Purpose: Development of a parametric equation suitable for a daily use in routine clinic to provide estimates of stray neutron doses in proton therapy. Methods: Monte Carlo (MC) calculations using the UF-NCI 1-year-old phantom were exercised to determine the variation of stray neutron doses as a function of irradiation parameters while performing intracranial treatments. This was done by individually changing the proton beam energy, modulation width, collimator aperture and thickness, compensator thickness and the air gap size while their impact on neutron doses were put into a single equation. The variation of neutron doses with distance from the target volume was also included in it. Then, a first step consisted in establishing the fitting coefficients by using 221 learning data which were neutron absorbed doses obtained with MC simulations while a second step consisted in validating the final equation. Results: The variation of stray neutron doses with irradiation parameters were fitted with linear, polynomial, etc. model while a power-law model was used to fit the variation of stray neutron doses with the distance from the target volume. The parametric equation fitted well MC simulations while establishing fitting coefficients as the discrepancies on the estimate of neutron absorbed doses were within 10%. The discrepancy can reach ∼25% for the bladder, the farthest organ from the target volume. Finally, the validation showed results in compliance with MC calculations since the discrepancies were also within 10% for head-and-neck and thoracic organs while they can reach ∼25%, again for pelvic organs. Conclusion: The parametric equation presents promising results and will be validated for other target sites as well as other facilities to go towards a universal method.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Conformation and elasticity of a charged polymer chain bridging two nanoparticles

A complex composed of a charged flexible polymer chain irreversibly attached with its ends to surfaces of two nanoparticles was investigated using the Metropolis Monte Carlo method on a simple cubic lattice. The simulations were performed in the presence of explicit ions. The bridging chain and the nanoparticles bearing the same and the opposite sign charges were considered. Changes in the free energy of the complex upon its stretching or compression, together with the magnitude of the elastic force, were examined. The relative roles of energetic and entropic effects in determining the properties of the complex were identified. Also, the adsorption of charged monomers on the opposite-sign charged nanoparticles and its influence on the examined quantities was studied. Moreover, a simple semi-analytical approach to the thermodynamics of the polymer bridge was derived.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗