Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Modelling”

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 127 records · Page 7

Radiation modeling with direct simulation Monte Carlo

Improvements in the modeling of radiation in low density shock waves with direct simulation Monte Carlo (DSMC) are the subject of this study. A new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. This scheme attempts to move the DSMC programs toward a more detailed modeling of the physics and more reliance on available experimental data. The new method is compared with the current modeling technique and both techniques are compared with available data. The differences in the results are evaluated. The test case is based on an AVCO-Everett shock tube experiment, a 10-km/s standing shock wave in air at 0.1 Torr. The new method agrees with the available data as well as the results from the earlier scheme and is more easily extrapolated to different flow conditions.

Carlson, Ann B.↗

Radiation Modeling with Direct Simulation Monte Carlo

Improvements in the modeling of radiation in low density shock waves with direct simulation Monte Carlo (DSMC) are the subject of this study. A new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. This scheme attempts to move the DSMC programs toward a more detailed modeling of the physics and more reliance on available rate data. The new method is compared with the current modeling technique and both techniques are compared with available experimental data. The differences in the results are evaluated. The test case is based on experimental measurements from the AVCO-Everett Research Laboratory electric arc-driven shock tube of a normal shock wave in air at 10 km/s and .1 Torr. The new method agrees with the available data as well as the results from the earlier scheme and is more easily extrapolated to di erent ow conditions.

Carlson, Ann B.↗

Strategies for Automation of Model Tuning in Multifidelity Trajectory Uncertainty Propagation

Multi-model Monte Carlo methods are efficient strategies to perform forward uncertainty quantification studies in entry, descent, and landing (EDL) applications. These multi-model methods are based on the classical Monte Carlo estimator, but fuse predictions from several low-fidelity models to obtain estimators with greater precision given a prescribed computational budget. The effectiveness of these approaches relies on the magnitudes of correlations between the low-fidelity models and the high-fidelity model, as well as the relative computational costs of all models. Identifying and exploiting the best trade-off between correlation and cost, which ultimately depends on the selection of hyperparameters in the low-fidelity models, is a task often performed by hand or simply inspired by the deterministic understanding available for a specific application. This work extends a preliminary effort,

Marten Thompson↗

Computation of Domain-Averaged Shortwave Irradiance by a One-Dimensional Algorithm Incorporating Correlations between Optical Thickness and Direct Incident Radiation

A one-dimensional radiative transfer algorithm that accounts for correlations between the optical thickness and the incident direct solar radiation is developed to compute the domain-averaged shortwave irradiance profile. It divides the direct irradiance into four components and treats the direct irradiance in two separate, clear and cloudy columns to account for the fact that clouds attenuate the direct irradiance more than clear-sky. The horizontal inhomogeneity of clouds in the cloudy column is treated by the gamma weighted two-stream approximation, which assumes that the optical thickness of clouds follows a gamma distribution. The algorithm inputs the cloud fraction, cumulative cloud fraction as a function of height, and a parameter expressing the shape of the probability density function of the cloud optical thickness distribution in addition to inputs required for a two-stream radiative transfer model. These cloud property inputs can be obtained using ground- and satellite-based instruments. Therefore, the algorithm can treat realistic cloud overlap features and horizontal inhomogeneity of clouds in a framework of one- dimensional radiative transfer. Heating rates computed by the algorithm using cloud fields generated by cloud resolving models agree with those computed with a Monte Carlo model. If optical properties in computational layers that divide a vertically extensive cloud are correlated, the irradiance profile computed by the algorithm further improves.

Kato, S.↗

Modeling of thin-film GaAs growth

A solid Monte Carlo model is constructed for the simulation of crystal growth. The model assumes thermally accommodated adatoms impinge upon the surface during a delta time interval. The surface adatoms are assigned a random energy from a Boltzmann distribution, and this energy determines whether the adatoms evaporate, migrate, or remain stationary during the delta time interval. For each addition or migration of an adatom, potential wells are adjusted to reflect the absorption, migration, or desorption potential changes.

Heinbockel, J. H.↗

Modeling of thin film GaAs growth

A potential scaling Monte Carlo model of crystal growth is developed. The model is a modification of the solid-on-solid method for studying crystal growth in that potentials at surface sites are continuously updated on a time scale reflecting the surface events of migration, incorporation and evaporation. The model allows for B on A type of crystal growth and lattice disregistry by the assignment of potential values at various surface sites. The surface adatoms are periodically assigned a random energy from a Boltzmann distribution and this energy determines whether the adatoms evaporate, migrate or remain stationary during the sampling interval. For each addition or migration of an adatom, the surface potentials are adjusted to reflect the adsorption, migration or desorption potential changes.

Heinbockel, J. H.↗

Space-based solar power conversion and delivery systems study. Volume 3: Economic analysis of space-based solar power systems

A variety of economic and programmatic issues are discussed concerning the development and deployment of a fleet of space-based solar power satellites (SSPS). The costs, uncertainties and risks associated with the current photovoltaic SSPS configuration, and with issues affecting the development of an economically viable SSPS development program are analyzed. The desirability of a low earth orbit (LEO) demonstration satellite and a geosynchronous (GEO) pilot satellite is examined and critical technology areas are identified. In addition, a preliminary examination of utility interface issues is reported. The main focus of the effort reported is the development of SSPS unit production, and operation and maintenance cost models suitable for incorporation into a risk assessment (Monte Carlo) model (RAM). It is shown that the key technology area deals with the productivity of man in space, not, as might be expected, with some hardware component technology.

Hazelrigg, G. A., Jr.↗

Space-based solar power conversion and delivery systems study. Volume 5: Economic analysis

Space-based solar power conversion and delivery systems are studied along with a variety of economic and programmatic issues relevant to their development and deployment. The costs, uncertainties and risks associated with the current photovoltaic Satellite Solar Power System (SSPS) configuration, and issues affecting the development of an economically viable SSPS development program are addressed. In particular, the desirability of low earth orbit (LEO) and geosynchronous (GEO) test satellites is examined and critical technology areas are identified. The development of SSPS unit production (nth item), and operation and maintenance cost models suitable for incorporation into a risk assessment (Monte Carlo) model (RAM) are reported. The RAM was then used to evaluate the current SSPS configuration expected costs and cost-risk associated with this configuration. By examining differential costs and cost-risk as a function of postulated technology developments, the critical technologies, that is, those which drive costs and/or cost-risk, are identified. It is shown that the key technology area deals with productivity in space, that is, the ability to fabricate and assemble large structures in space, not, as might be expected, with some hardware component technology.

Source record↗

Hydrogen Balmer alpha intensity distributions and line profiles from multiple scattering theory using realistic geocoronal models

The H Balmer alpha nightglow is investigated by using Monte Carlo models of asymmetric geocoronal atomic hydrogen distributions as input to a radiative transfer model of solar Lyman-beta radiation in the thermosphere and atmosphere. It is shown that it is essential to include multiple scattering of Lyman-beta radiation in the interpretation of Balmer alpha airglow data. Observations of diurnal variation in the Balmer alpha airglow showing slightly greater intensities in the morning relative to evening are consistent with theory. No evidence is found for anything other than a single sinusoidal diurnal variation of exobase density. Dramatic changes in effective temperature derived from the observed Balmer alpha line profiles are expected on the basis of changing illumination conditions in the thermosphere and exosphere as different regions of the sky are scanned.

Anderson, D. E., Jr.↗

Kinetic Monte Carlo Tool for Kinetic Modeling of Linear Step-Growth Polymerization: Insight into Recycling of Polyurethanes

A kinetic Monte Carlo model of polyurethane polymerization which explicitly tracks the polymer sequences is developed and shared. Here this model is benchmarked against theoretical and experimental polyurethane data and used to investigate the effect on oligomer distributions of unequal reactivity of the first and second isocyanate to react. The reverse reactions using thermodynamic consistency are then added to the framework, and analogous to the addition polymerization concept of ceiling temperature, equilibrium chain length distributions at various temperatures are calculated. For a mixture of three monomers AA, BB, and CC, where BB and CC do not react with one another, are present in stoichiometric proportions, and have different enthalpies of reaction with AA, an odd-even effect emerges. Odd length chains are more likely than even length chains for temperatures at which BB and CC have significantly different equilibrium conversions. The concept of ceiling temperature that is typically cited for addition polymers is extended here to provide a measure of conditions under which depolymerization for recycling is favored.

condensation polymerization↗

The effect of finite geometry on the three-dimensional transfer of solar irradiance in clouds

Results are presented for a Monte Carlo model applied to a wide range of cloud widths and heights, and for an analytical model restricted in its application to cuboidally shaped clouds whose length, breadth, and depth may be varied independently; the clouds must be internally homogeneous with respect to their intrinsic radiative properties. Comparative results from the Monte Carlo method and the derived analytical model are presented for a wide range of cloud sizes, with special emphasis on the effects of varying the single scatter albedo, the solar zenith angle, and the scattering phase angle.

Davies, R.↗

Investigation of Magmatic Activities on Early Mars Using Igneous Mineral Chemistry in Gale Crater, Mars

One objective of rover missions is exploring the geological context of the surroundings. Over the years, igneous petrology and sedimentology have been disconnected, the first investigating magmatic processes and volcanic activities, and the second seeking environmental conditions in the past and assessing the habitability of the planet. Although different, one is related to the other: igneous rocks are altered and broken down, leading to the formation of sedimentary rocks, which can in turn be used to back out the nature of their magmatic source. The Curiosity rover that landed in the 3.7 Gyr old impact crater Gale is traveling through sedimentary rocks. About fifty float rocks have been observed, and several of them with ambiguous texture and composition have been classified as igneous or sedimentary depending on studies such as Jake_M. The composition of several unambiguous igneous rocks has been analyzed [4- 6] but their heterogeneity at a larger instrumental (measurement size < 2 cm) scale prevents the measurement of a bulk composition as performed on Earth. An original approach avoiding these two last issues is to consider igneous mineral chemistry analyzed within igneous and sedimentary rocks to assess magmatic processes that could have formed them. Most Curiosity data are used to explore ancient environmental conditions, and a significant number of compositional analyses are under-explored for constraining magmatic activities. We will present how we can make use of sedimentary data for investigating igneous processes in the vicinity of Gale crater. Geological Context: We focus on the first 750 martian days, corresponding to measurements in a coherent lacustrine sedimentary unit called Bradbury, because all sedimentary rocks were sourced from the same watershed and appear to have a consistent source with minimal alteration [2-3]. Igneous detrital minerals including feldspar and pyroxene, are observed in sedimentary rocks. Monte Carlo models showed that minimal cation loss is observed based on the composition of all Bradbury rocks, implying negligible weathering [3]. Although clay minerals are detected in few rocks [7], chemical compositions of rocks can be explained by a mixture of primary igneous minerals [3]. Variation of composition within Bradbury rocks can be explained by mineral sorting and one distinct source component. While a common magmatic source is suggested, Bradbury sediments likely come from several volcanic eruptions from a single magmatic chamber [9- 10]. The occurrence of alkali minerals like sanidine and K-rich rocks throughout Bradbury supports the presence of a potassic component, likely trachytic, while plagioclase and a mafic composition suggest a basaltic component [8-9]. Instruments: Mineral chemistry can be estimated by three instruments onboard Curiosity. The CheMin instrument enables detection of mineral assemblages using X-ray diffraction (XRD). Using Rietveld refinement, each mineral is identified according to their 1D XRD pattern [11]. Note that distinction between pyroxene minerals is challenging with the CheMin instrument due to overlapping peaks on XRD patterns and low angular resolution of the instrument [12]. Then, using least square regression and optimization algorithms based on unit-cell parameters, mineral chemistry has been estimated by [11]. Plagioclase composition has been estimated using the NaAlSi3O8- CaAl2Si2O8 system and alkali feldspar is based on the NaAlSi3O8-KAlSi3O8 system (stars in Fig.1). Two mudstone samples (John Klein and Cumberland) and one sandstone sample (Windjana) were analyzed by CheMin at Yellowknife Bay and Kimberley, respectively. Figure 1. Ternary diagrams of feldspars (top) and pyroxene (bottom) quadrilateral. Stars correspond to CheMin composition and the gray patches to ChemCam composition. The colored dots are the composition of feldspar and pyroxene that crystallized during fractional crystallization at FMQ+1 of a melt extracted at distinct melting degree during the adiabatic ascent of a primitive mantle composition, without any water (left panels) and with 0.5 wt.% of water (right panels) at distinct pressure. The ChemCam instrument enables the analysis of the chemical compositions of rocks at hundreds of micrometer scale (350-550 μm) using laser induced breakdown spectroscopy (LIBS), which may provide the composition of minerals when they are larger than the beam spot (>550 μm) [13]. Within >5000 LIBS points, we performed a typical stoichiometric filtering allowing us to distinguish 56 feldspar and 10 pyroxene mineral compositions (grey patches in Fig. 1). Finally, the Alpha Particle X-ray Spectrometer (APXS) analyzes the composition of rocks with a 1.6 cm diameter spot size. Monte Carlo mass balance modeling allowed [3] to decipher a feldspar range varying between An30 and An40 (Fig. 1). Discussion: Although there could be a more complex history and other ways to form the whole compositional range of igneous minerals analyzed within the Bradbury formation, we are presenting here simple magmatic pathways commonly occurring on Earth using the thermodynamical softwares pMELTS and rhyoliteMELTS [14]. The objective is to find reasonable igneous processes that produce minerals that parallel the compositions of feldspar and pyroxene analyzed by the Curiosity rover. As commonly observed for mid-ocean ridge basalts, the adiabatic ascent of a primitive mantle composition [15] partially melting at 2 GPa has been modeled, followed by the extraction of a liquid at distinct degrees of partial melting, which undergoes fractional crystallization at an oxygen fugacity +1 log unit above the fayalite-magnetite-quartz (FMQ) buffer within the crust (0.02-0.4 GPa) with H2O = 0- 0.5 wt. %. These latter conditions correspond to those recorded within igneous clasts from the Noachian martian breccia NWA 7034 and paired and within Gale igneous rocks (colored dots in Fig. 1) [16-17]. To check the reliability of these 2-step models, we also tested fractional crystallization at similar conditions (FMQ+1; P=0.02-0.4 GPa; H2O = 0-0.5 wt. %) of starting compositions corresponding to that of magmas with distinct melting degrees obtained from isobaric experiments at 2 GPa [18]. Mineral compositions obtained from both models are similar. As shown on Fig. 1, the whole range of observed feldspar compositions cannot be reproduced by fractionation of one magma only. Indeed, while alkali feldspar and Na-plagioclase likely crystallized from fractional crystallization of a low-degree melt (here <15%), plagioclase and pyroxene can only be formed by fractional crystallization of a higher degree melt (here >19%). The corresponding liquid descent lines are broadly in agreement with compositions estimated by ChemCam corresponding to float igneous rocks (Fig. 2) [4-6]. Figure 2. Silica versus alkali content. Lines show the liquid lines of descent from magmas with distinct degrees of melting. Gray patches represents the composition of Gale igneous rocks [4-6]. Trachytic to rhyolitic magmas crystallize alkali feldspar, and andesite to dacite magmas likely form plagioclase. Therefore, at least two starting magmas at distinct melting degrees, which could easily come from a single mantle source, are necessary to explain the whole compositional range of feldspar and pyroxene analyzed within Bradbury rocks. Conclusion: Because rocks from the Bradbury formation are likely originating from the same magmatic source with minimal weathering as supported by several studies using different approaches, igneous mineral chemistry analyzed by CheMin and ChemCam allows us to back out reasonable magmatic pathways that could have crystallized them. Fractional crystallization of at least two starting magmas originating from distinct melting degrees of a single mantle source can explain the whole range of feldspar and pyroxene composition. Both alkaline and sub-alkaline liquids can be produced, with compositions corresponding to those of the igneous rocks analyzed by ChemCam within the Bradbury formation, highlighting the complexity of Mars magmatism.

Payre, V.↗

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ozone measurement systems improvements studies

Results are summarized of an initial study of techniques for measuring atmospheric ozone, carried out as the first phase of a program to improve ozone measurement techniques. The study concentrated on two measurement systems, the electro chemical cell (ECC) ozonesonde and the Dobson ozone spectrophotometer, and consisted of two tasks. The first task consisted of error modeling and system error analysis of the two measurement systems. Under the second task a Monte-Carlo model of the Dobson ozone measurement technique was developed and programmed for computer operation.

Thomas, R. W.↗

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made substantial progress toward the goal of producing a unified model of the basic physical processes responsible for solar wind acceleration. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a 1D model of plasma heating and acceleration. The accomplishments in Year 2 are divided into these two categories: 1a. Focused Study of Kinetic Magnetohydrodynamic (MHD) Turbulence. lb. Focused Study of Non - WKB Alfven Wave Rejection. and 2. The Unified Model Code. We have continued the development of the computational model of a time-study open flux tube in the extended corona. The proton-electron Monte Carlo model is being tested, and collisionless wave-particle interactions are being included. In order to better understand how to easily incorporate various kinds of wave-particle processes into the code, the PI performed a detailed study of the so-called "Ito Calculus", i.e., the mathematical theory of how to update the positions of particles in a probabilistic manner when their motions are governed by diffusion in velocity space.

Cranmer, Steven R.↗

Investigation of the significance of excited species in wave-heated dense non-ideal microplasmas using particle-in-cell Monte Carlo-collision modeling

A computational model for the interaction of a high-pressure microplasma with an electromagnetic wave is presented. A one-dimensional particle-in-cell Monte Carlo collision model is used to investigate the plasma non-ideality effects in a second-stage laser-heated xenon plasma with a comprehensive chemistry mechanism, including excited species. A xenon microplasma at a temperature of 300 K and pressure of 10 bar becomes non-ideal after the laser heating stage with lower ionization and electronic excitation barriers, which makes the presence of excited species important to the plasma generation process. For these conditions, two-step ionization via excited species becomes a major contributor to plasma ionization. We find that full ionization of the plasma is obtained in ~2 ps when excited species are included in the chemistry model, as opposed to ~2.5 ps for a chemistry mechanism that does not consider their effect. With excited species, the ionization mode transitions from direct ionization via electron-neutral collisions to direct and stepwise ionization as the plasma generation progresses.

36 MATERIALS SCIENCE↗

Studying strangeness and baryon production mechanisms through angular correlations between charged Ξ baryons and identified hadrons in pp collisions at $\sqrt{s}$ = 13 TeV

The angular correlations between charged Ξ baryons and associated identified hadrons (pions, kaons, protons, Λ baryons, and Ξ baryons) are measured in pp collisions at $\sqrt{s}$ = 13 TeV with the ALICE detector to give insight into the particle production mechanisms and balancing of quantum numbers on the microscopic level. In particular, the distribution of strangeness is investigated in the correlations between the doubly-strange Ξ baryon and mesons and baryons that contain a single strange quark, K and Λ. As a reference, the results are compared to Ξπ and Ξp correlations, where the associated mesons and baryons do not contain a strange valence quark. These measurements are expected to be sensitive to whether strangeness is produced through string breaking or in a thermal production scenario. Furthermore, the multiplicity dependence of the correlation functions is measured to look for the turn-on of additional particle production mechanisms with event activity. The results are compared to predictions from the string-breaking model P YTHIA 8, including tunes with baryon junctions and rope hadronisation enabled, the cluster hadronisation model H ERWIG 7, and the core-corona model E POS-LHC . While some aspects of the experimental data are described quantitatively or qualitatively by the Monte Carlo models, no model can match all features of the data. These results provide stringent constraints on the strangeness and baryon number production mechanisms in pp collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Forward and inverse predictive transient models of TREAT using surrogate reactivity models

In this work, we present two novel approaches for predicting the power evolution and control rod height of the Transient Reactor Test Facility (TREAT) to support experiment modeling; specifically, transient analysis of the NASA-sponsored Sirius series of experiments. These approaches utilize steady-state Monte Carlo model, point kinetics model, and surrogate models to predict power evolution and control rod axial position during transient experiments. Both approaches were tested and validated against several Sirius experiments that were performed in the TREAT facility at different power levels. The validation test results show very good agreement with the experimental data, and the models were able to accurately predict the power evolution and the axial control rod position with an average error within 3.0%. Finally, this indicates that these approaches will help the reactor engineering team of the TREAT facility in preparing and predicting the power and temperature of the experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗