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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Development of the X-33 Aerodynamic Uncertainty Model

An aerodynamic uncertainty model for the X-33 single-stage-to-orbit demonstrator aircraft has been developed at NASA Dryden Flight Research Center. The model is based on comparisons of historical flight test estimates to preflight wind-tunnel and analysis code predictions of vehicle aerodynamics documented during six lifting-body aircraft and the Space Shuttle Orbiter flight programs. The lifting-body and Orbiter data were used to define an appropriate uncertainty magnitude in the subsonic and supersonic flight regions, and the Orbiter data were used to extend the database to hypersonic Mach numbers. The uncertainty data consist of increments or percentage variations in the important aerodynamic coefficients and derivatives as a function of Mach number along a nominal trajectory. The uncertainty models will be used to perform linear analysis of the X-33 flight control system and Monte Carlo mission simulation studies. Because the X-33 aerodynamic uncertainty model was developed exclusively using historical data rather than X-33 specific characteristics, the model may be useful for other lifting-body studies.

Cobleigh, Brent R.↗
NASA NTRS2011

The Distribution of Chromosomal Aberrations in Human Cells Predicted by a Generalized Time-Dependent Model of Radiation-Induced Formation of Aberrations

New experimental data show how chromosomal aberrations for low- and high-LET radiation are dependent on DSB repair deficiencies in wild-type, AT and NBS cells. We simulated the development of chromosomal aberrations in these cells lines in a stochastic track-structure-dependent model, in which different cells have different kinetics of DSB repair. We updated a previously formulated model of chromosomal aberrations, which was based on a stochastic Monte Carlo approach, to consider the time-dependence of DSB rejoining. The previous version of the model had an assumption that all DSBs would rejoin, and therefore we called it a time-independent model. The chromosomal-aberrations model takes into account the DNA and track structure for low- and high-LET radiations, and provides an explanation and prediction of the statistics of rare and more complex aberrations. We compared the program-simulated kinetics of DSB rejoining to the experimentally-derived bimodal exponential curves of the DSB kinetics. We scored the formation of translocations, dicentrics, acentric and centric rings, deletions, and inversions. The fraction of DSBs participating in aberrations was studied in relation to the rejoining time. Comparisons of simulated dose dependence for simple aberrations to the experimental dose-dependence for HF19, AT and NBS cells will be made.

Ponomarev, Artem L.↗
NASA NTRS2004

Mars Exploration Rover Terminal Descent Mission Modeling and Simulation

Because of NASA's added reliance on simulation for successful interplanetary missions, the MER mission has developed a detailed EDL trajectory modeling and simulation. This paper summarizes how the MER EDL sequence of events are modeled, verification of the methods used, and the inputs. This simulation is built upon a multibody parachute trajectory simulation tool that has been developed in POST I1 that accurately simulates the trajectory of multiple vehicles in flight with interacting forces. In this model the parachute and the suspended bodies are treated as 6 Degree-of-Freedom (6 DOF) bodies. The terminal descent phase of the mission consists of several Entry, Descent, Landing (EDL) events, such as parachute deployment, heatshield separation, deployment of the lander from the backshell, deployment of the airbags, RAD firings, TIRS firings, etc. For an accurate, reliable simulation these events need to be modeled seamlessly and robustly so that the simulations will remain numerically stable during Monte-Carlo simulations. This paper also summarizes how the events have been modeled, the numerical issues, and modeling challenges.

Raiszadeh, Behzad↗
DOE OSTI2024

Measurement of the Pion Exclusive Electro-Production Cross-Section in the E12-19-006 Experiment in Hall-C at Jefferson Lab

One of the most effective methods for exploring the transition from hadronic degrees of freedom to quark-gluon degrees of freedom in Quantum Chromodynamics (QCD) is through the investigation of \exclusive" pion and kaon electro-production reactions at various Q2 and ?t values. The E12-19-006 experiment is conducted within the confi?nes of experimental Hall C at the Thomas Jefferson National Accelerator Facility, USA, for such studies. The primary aim of the experiment is to ?first enhance our comprehension of the pion electro-production cross-section and its form factor at Q2 = 0.38 and 0.42 GeV2. This is the fi?rst run period of the E12-19-006 experiment which ran in summer 2019. A more profound understanding of the pion electro-production reaction, 1H(e,e'?+)n, at low Q2 is deemed essential to employ this electro-production reaction (an indirect technique) for the high Q2 studies, thereby delving deeper into the realm of QCD. Consequently, this dissertation presents a thorough analysis of the experimental data acquired in the ?first run period of the E12-19-006 experiment. In pursuit of precision, a series of systematic studies (target boiling correction study, the elastic reaction cross-section measurements, study for determining vari- ous kinematics o?sets, etc.) are conducted to discern the accuracy of the analyzed data, a prerequisite for the use of Rosenbluth separation technique to separate the pion electro-production cross-section terms in t bins. The separated pion electro-production cross-section through the Rosenbluth separation technique is then used to extract the pion electromagnetic form factor. In this dissertation, the pion electro-production cross-section is carefully dissected into its four constituent components: longitudinal (?L), transverse (?T ), longitudinal-transverse (?LT ), and transverse-transverse (?TT ), using the full version of Rosenbluth separation technique for the Q2 = 0.38 GeV2. The technique is simultaneously fi?tted to the unseparated pion electro-production cross-sections at the three values of polarization of the virtual photon (?), i.e., ? = 0.286, 0.629 and 0.781. An iterative process is applied to re?ne the parameters of the model cross-sections until the yield ratio of experimental and Monte Carlo simulation converges. In this study, 21 iterations are conducted to re?ne the model cross-section parameters. The fi?nal pion electro-production cross-section terms are then determined for 7 t bins using the optimized parameters of the model cross-sections.

Kumar, Vijay↗
NASA NTRS2012

Mapping the Milky Way Galaxy with LISA

Gravitational wave detectors in the mHz band (such as the Laser Interferometer Space Antenna, or LISA) will observe thousands of compact binaries in the galaxy which can be used to better understand the structure of the Milky Way. To test the effectiveness of LISA to measure the distribution of the galaxy, we simulated the Close White Dwarf Binary (CWDB) gravitational wave sky using different models for the Milky Way. To do so, we have developed a galaxy density distribution modeling code based on the Markov Chain Monte Carlo method. The code uses different distributions to construct realizations of the galaxy. We then use the Fisher Information Matrix to estimate the variance and covariance of the recovered parameters for each detected CWDB. This is the first step toward characterizing the capabilities of space-based gravitational wave detectors to constrain models for galactic structure, such as the size and orientation of the bar in the center of the Milky Way

McKinnon, Jose A.↗
NASA NTRS1983

A potential-energy scaling model to simulate the initial stages of thin-film growth

A solid on solid (SOS) Monte Carlo computer simulation employing a potential energy scaling technique was used to model the initial stages of thin film growth. The model monitors variations in the vertical interaction potential that occur due to the arrival or departure of selected adatoms or impurities at all sites in the 400 sq. ft. array. Boltzmann ordered statistics are used to simulate fluctuations in vibrational energy at each site in the array, and the resulting site energy is compared with threshold levels of possible atomic events. In addition to adsorption, desorption, and surface migration, adatom incorporation and diffusion of a substrate atom to the surface are also included. The lateral interaction of nearest, second nearest, and third nearest neighbors is also considered. A series of computer experiments are conducted to illustrate the behavior of the model.

Heinbockel, J. H.↗
DOE OSTI2025

Multiclass Classification Using Bayesian Multivariate Adaptive Regression Splines

We present a new Bayesian model for the problem of multiclass classification. In this model, the probabilities of class membership of a given observation are determined by the mean of a latent Gaussian distribution. The mean functions of this latent distribution consist of combinations of highly flexible basis functions of the inputs: multivariate adaptive regression splines (MARS), first developed for multiple regression. We use reversible jump Markov chain Monte Carlo to make inference on the classification model, including the number of basis functions. We compare the probabilistic classification performance of our proposed approach to existing methods on simulated and benchmark data, and compare uncertainty estimates on simulated data. Our proposed method compares favorably with existing Bayesian and frequentist multiclass classification methods in out-of-sample probabilistic classification, and uncertainty estimation of these probabilistic classifications. We examine the fit of the proposed method to a data set of hurricane storm surge levels near Delaware Bay, US, and conclude that sea level rise is a key contributor to damage delivered by storm surge.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Modeling Enhancements, Cross-Section Generation Updates, and Benchmarking with Shift

This technical report documents the modeling enhancements, cross-section generation updates, and bench marking with the Shift Monte Carlo code performed under the US Department of Energy Nuclear Energy Advanced Modeling and Simulation Program in FY 2024. The work performed included several modeling enhancements, such as integration of cross-section generation in Titan and the ability to produce microscopic multigroup cross sections with Shift. Benchmarking of the cross sections produced by Shift and the two-step workflow with Griffin was performed for three problems: the Advanced Breeder Test Reactor, a generic pebble bed reactor, and a TRISO heat pipe microreactor. Comparisons of results from these benchmark problems were done with Serpent, OpenMC, and Griffin. These enhancements provide a robust foundation for applying Shift for both reference and two-step neutronics analysis for advanced reactor simulation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2008

SEE Rate Estimation: Model Complexity and Data Requirements

Statistical Methods outlined in [Ladbury, TNS20071 can be generalized for Monte Carlo Rate Calculation Methods Two Monte Carlo Approaches: a) Rate based on vendor-supplied (or reverse-engineered) model SEE testing and statistical analysis performed to validate model; b) Rate calculated based on model fit to SEE data Statistical analysis very similar to case for CREME96. Information Theory allows simultaneous consideration of multiple models with different complexities: a) Model with lowest AIC usually has greatest predictive power; b) Model averaging using AIC weights may give better performance if several models have similar good performance; and c) Rates can be bounded for a given confidence level over multiple models, as well as over the parameter space of a model.

Ladbury, Ray↗
DOE OSTI2022

SERPENT 2 and TRIPOLI-4 neutronics benchmarking of the Jules Horowitz reactor

Jules Horowitz Reactor (JHR) is a future center for nuclear material testing in Europe, currently under construction at CEA Cadarache in southern France. It has been modelled extensively with CEA's HORUS3D/N neutron calculation package, that includes a Monte-Carlo neutron transport code TRIPOLI-4. In order to support the future computational analyses on the reactor core physics, VTT's Monte Carlo neutronics code SERPENT 2 has been used to create an alternative model of the JHR, taking advantage of the pre-existing Constructive Solid Geometry (CSG) structures on which the TRIPOLI-4 model was built. This work describes a benchmark, in which two JHR core configurations were analyzed using both of these models. The first configuration is a fresh fuel core at the beginning of life, in which aluminum based dummy-devices are installed in the reactor test positions. The second configuration represents the end-of-cycle equilibrium core with a heterogeneous burnup distribution, in which CALIPSO devices containing cladding material samples are installed at the in-core experimental positions, while the reflector positions are occupied by DEN-REP devices containing UO{sub 2} fuel. Neutron flux spectra and heat deposition are the parameters of interest in this study, due to their high importance in creating the desired experimental irradiation environment. The impact of delayed gamma energy deposition will be studied in detail, since the two models account for it with different methods. This work aims to establish a SERPENT 2 model of JHR, understanding its differences and agreements with TRIPOLI-4, and discussing its potential future applications. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Learning stochastic dynamics and predicting emergent behavior using transformers

We show that a neural network originally designed for language processing can learn the dynamical rules of a stochastic system by observation of a single dynamical trajectory of the system, and can accurately predict its emergent behavior under conditions not observed during training. We consider a lattice model of active matter undergoing continuous-time Monte Carlo dynamics, simulated at a density at which its steady state comprises small, dispersed clusters. We train a neural network called a transformer on a single trajectory of the model. The transformer, which we show has the capacity to represent dynamical rules that are numerous and nonlocal, learns that the dynamics of this model consists of a small number of processes. Forward-propagated trajectories of the trained transformer, at densities not encountered during training, exhibit motility-induced phase separation and so predict the existence of a nonequilibrium phase transition. Transformers have the flexibility to learn dynamical rules from observation without explicit enumeration of rates or coarse-graining of configuration space, and so the procedure used here can be applied to a wide range of physical systems, including those with large and complex dynamical generators.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2012

Modeling Phase-Aligned Gamma-Ray and Radio Millisecond Pulsar Light Curves

Since the discovery of the first eight gamma-ray millisecond pulsars (MSPs) by the Fermi Large Area Telescope, this population has been steadily expanding. Four of the more recent detections, PSR J00340534, PSR J1939+2134 (B1937+21; the first MSP ever discovered), PSR J1959+2048 (B1957+20; the first discovery of a black widow system), and PSR J2214+3000, exhibit a phenomenon not present in the original discoveries: nearly phase-aligned radio and gamma-ray light curves (LCs). To account for the phase alignment, we explore models where both the radio and gamma-ray emission originate either in the outer magnetosphere near the light cylinder or near the polar caps. Using a Markov Chain Monte Carlo technique to search for best-fit model parameters, we obtain reasonable LC fits for the first three of these MSPs in the context of altitude-limited outer gap (alOG) and two-pole caustic (alTPC) geometries (for both gamma-ray and radio emission). These models differ from the standard outer gap (OG)/two-pole caustic (TPC) models in two respects: the radio emission originates in caustics at relatively high altitudes compared to the usual conal radio beams, and we allow both the minimum and maximum altitudes of the gamma-ray and radio emission regions to vary within a limited range (excluding the minimum gamma-ray altitude of the alTPC model, which is kept constant at the stellar radius, and that of the alOG model, which is set to the position-dependent null charge surface altitude). Alternatively, phase-aligned solutions also exist for emission originating near the stellar surface in a slot gap scenario (low-altitude slot gap (laSG) models). We find that the alTPC models provide slightly better LC fits than the alOG models, and both of these give better fits than the laSG models (for the limited range of parameters considered in the case of the laSG models). Thus, our fits imply that the phase-aligned LCs are likely of caustic origin, produced in the outer magnetosphere, and that the radio emission for these pulsars may come from close to the light cylinder. In addition, we were able to constrain the minimum and maximum emission altitudes with typical uncertainties of 30% of the light cylinder radius. Our results therefore describe a third gamma-ray MSP subclass, in addition to the two previously found by Venter et al.: those with LCs fit by standard OG/TPC models and those with LCs fit by pair-starved polar cap models.

Venter, C.↗
NASA NTRS2017

A Quantitative Risk-Benefit Analysis of Prophylactic Surgery Prior to Extended-Duration Spaceflight

INTRODUCTION: Among otherwise healthy astronauts undertaking deep space missions, the risks for acute appendicitis (AA) and cholecystitis (AC) are not zero. If these conditions were to occur during spaceflight they may require surgery for definitive care. The proposed study quantifies and compares the risks of developing de novo AA and AC in-flight to the surgical risks of prophylactic laparoscopic appendectomy (LA) and cholecystectomy (LC) using NASA's Integrated Medical Model (IMM). METHODS: The IMM is a Monte Carlo simulation that forecasts medical events during spaceflight missions and estimates the impact of these medical events on crew health. In this study, four Design Reference Missions (DRMs) were created to assess the probability of an astronaut developing in-flight small-bowel obstruction (SBO) following prophylactic 1) LA, 2) LC, 3) LA and LC, or 4) neither surgery (SR# S-20160407-351). Model inputs were drawn from a large, population-based 2011 Swedish study that examined the incidence and risks of post-operative SBO over a 5-year follow-up period. The study group included 1,152 patients who underwent LA, and 16,371 who underwent LC. RESULTS: Preliminary results indicate that prophylactic LA may yield higher mission risks than the control DRM. Complete analyses are pending and will be subsequently available. DISCUSSION: The risk versus benefits of prophylactic surgery in astronauts to decrease the probability of acute surgical events during spaceflight has only been qualitatively examined in prior studies. Within the assumptions and limitations of the IMM, this work provides the first quantitative guidance that has previously been lacking to this important question for future deep space exploration missions.

Carroll, Danielle↗
DOE OSTI2025

Combined Influence of Rotation and Scrape-Off Layer Drifts on Recycling Asymmetries in Tokamak Plasmas

Coupled 2D fluid-kinetic simulations of a DIII-D high confinement tokamak plasma show that plasma rotation coupled with drift effects near the plasma edge play a significant role in the creation of the observed poloidal distribution of neutrals. It is observed that including either drift or rotation effects enhances particle flux at the inner target in the case of ion 𝐵×∇𝐵 drift toward the 𝑋-point. However, the particle flux asymmetry is significantly higher with the combination of drifts and rotation than either effect alone. The heightened particle flux asymmetry allows for improved simulation of the strong in-out asymmetry of the Lyman-𝛼 brightness profiles measured in the experiment. Enhancement of radial transport of parallel momentum changes the upstream scrape-off layer flow pattern, increasing the fraction of deuterium flux that reaches the inboard divertor entrance while lowering that which arrives at the outboard. In conclusion, this Letter indicates that by combining drifts, rotation, and viscous coupling, existing boundary plasma models can achieve a satisfactory agreement with experimentally measured neutral asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2006

Determination of the Cosmic Distance Scale from Sunyaev-Zel'dovich Effect and Chandra X-ray Measurements of High Redshift Galaxy Clusters

We determine the distance to 38 clusters of galaxies in the redshift range 0.14 less than or equal to z less than or equal to 0.89 using X-ray data from Chandra and Sunyaev-Zeldovich Effect data from the Owens Valley Radio Observatory and the Berkeley-Illinois-Maryland Association interferometric arrays. The cluster plasma and dark matter distributions are analyzed using a hydrostatic equilibrium model that accounts for radial variations in density, temperature and abundance, and, the statistical and systematic errors of this method are quantified. The analysis is performed via a Markov chain Monte Carlo technique that provides simultaneous estimation of all model parameters. W

Bonamente, Massimiliano↗
NASA NTRS2000

A Validation Summary of the NCC Turbulent Reacting/non-reacting Spray Computations

This pper provides a validation summary of the spray computations performed as a part of the NCC (National Combustion Code) development activity. NCC is being developed with the aim of advancing the current prediction tools used in the design of advanced technology combustors based on the multidimensional computational methods. The solution procedure combines the novelty of the application of the scalar Monte Carlo PDF (Probability Density Function) method to the modeling of turbulent spray flames with the ability to perform the computations on unstructured grids with parallel computing. The calculation procedure was applied to predict the flow properties of three different spray cases. One is a nonswirling unconfined reacting spray, the second is a nonswirling unconfined nonreacting spray, and the third is a confined swirl-stabilized spray flame. The comparisons involving both gas-phase and droplet velocities, droplet size distributions, and gas-phase temperatures show reasonable agreement with the available experimental data. The comparisons involve both the results obtained from the use of the Monte Carlo PDF method as well as those obtained from the conventional computational fluid dynamics (CFD) solution. Detailed comparisons in the case of a reacting nonswirling spray clearly highlight the importance of chemistry/turbulence interactions in the modeling of reacting sprays. The results from the PDF and non-PDF methods were found to be markedly different and the PDF solution is closer to the reported experimental data. The PDF computations predict that most of the combustion occurs in a predominantly diffusion-flame environment. However, the non-PDF solution predicts incorrectly that the combustion occurs in a predominantly vaporization-controlled regime. The Monte Carlo temperature distribution shows that the functional form of the PDF for the temperature fluctuations varies substantially from point to point. The results also bring to the fore some of the deficiencies associated with the use of assumed-shape PDF methods in spray computations.

Raju, M. S.↗
DOE OSTI2022

Chemical order transitions within extended interfacial segregation zones in NbMoTaW

Interfacial segregation and chemical short-range ordering influence the behavior of grain boundaries in complex concentrated alloys. In this study, we use atomistic modeling of a NbMoTaW refractory complex concentrated alloy to provide insight into the interplay between these two phenomena. Hybrid Monte Carlo and molecular dynamics simulations are performed on columnar grain models to identify equilibrium grain boundary structures. Here our results reveal extended near-boundary segregation zones that are much larger than traditional segregation regions, which also exhibit chemical patterning that bridges the interfacial and grain interior regions. Furthermore, structural transitions pertaining to an A2-to-B2 transformation are observed within these extended segregation zones. Both grain size and temperature are found to significantly alter the widths of these regions. An analysis of chemical short-range order indicates that not all pairwise elemental interactions are affected by the presence of a grain boundary equally, as only a subset of elemental clustering types are more likely to reside near certain boundaries. The results emphasize the increased chemical complexity that is associated with near-boundary segregation zones and demonstrate the unique nature of interfacial segregation in complex concentrated alloys.

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

Uncertainty quantification of bank vegetation impacts on the flood flow field in the American River, California, using large‐eddy simulations

Bank vegetation plays a key role in both hydrodynamics and morphodynamics of natural rivers; however, these effects are often unaccounted for in the computational flow dynamics of natural waterways. Recent studies using the large‐eddy simulation (LES), however, have attempted to gain insights into the impacts of bank vegetation on the mean flow field of the natural rivers using a vegetation model, which applies a sink term to the momentum equations of motion. This approach accounts for the effects of the vegetation and provides a practical approach to account for the complex patches of bank vegetation in large‐scale rivers. To implement the vegetation model, a drag coefficient reflecting the overall resistance of vegetal structures to the flow is needed, but due to the lack of calibrated data and range of size, density and type of vegetation, this parameter can be a significant source of uncertainty in the model results. Here, in this study, we use uncertainty quantification (UQ) to investigate the hydrodynamics and bed shear results when a bank vegetation is incorporated in an LES model. To this end, we used the polynomial chaos expansion and Monte Carlo sampling techniques to determine the uncertainties associated with the drag coefficient in the vegetation model and from uncertainties in the bed roughness and inflow discharge. The UQ analysis provided spatially varying confidence levels for the spanwise and vertical distribution of velocity magnitude and for the bed shear stress distributions. In addition, Sobol indices were computed to indicate the relative influence that each parameter had on the overall uncertainty. In general, it was found that uncertainty in flow discharge was the dominant source of uncertainty; however, the drag coefficient in the vegetation model and the bed roughness parameter also made significant contribution to the uncertainty near the banks and bed, respectively.

54 ENVIRONMENTAL SCIENCES↗