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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Lunar Relay Onboard Navigation Performance and Effects on Lander Descent to Surface

A system of Lunar relay satellites has been proposed to address communication and navigation needs and ensure robustness for the variety of upcoming robotic and human exploration missions to the Moon. The relays are envisioned to estimate self position and time knowledge onboard, allowing the system to provide in-situ navigation services to missions in the Lunar and cis-Lunar environment. The quality and accuracy of those services are highly dependent on the navigation performance of the relay itself. To assess the Lunar relay navigation performance, a series of orbit determination (OD) Monte Carlo (MC) simulations are run using Lunar gravity modeling up to a degree and order of 250 and a variety of onboard clocks and measurement types including weak-signal GNSS, Ground Network (GN) pseudorange (PR) and Doppler, and optical navigation (OpNav) center-finding (CF). The estimated trajectories produced by these Lunar relay MC simulations, along with the associated errors, and transmitted navigation reference signal parameters, are used to evaluate the expected navigation performance of a user on a descent trajectory to the Lunar surface. The scenario features a lander system performing onboard navigation relying on one-way range and Doppler measurements from reference signals emitted by the Lunar relay. This paper can be used as a reference in determining the onboard clock and measurement types necessary to obtain acceptable navigation performance for the Lunar relay, and as a baseline for Lunar lander navigation performance using accurate measurements from a relay reference signal.

Jeffrey L. Small↗
NASA NTRSDate not supplied

Lunar Relay Onboard Navigation Performance and Effects on Lander Descent to Surface

A system of Lunar relay satellites has been proposed to address communication and navigation needs and ensure robustness for the variety of upcoming robotic and human exploration missions to the Moon. The relays are envisioned to estimate self-position and time knowledge onboard, allowing the system to provide in-situ navigation services to missions in the Lunar and cis-Lunar environment. The quality and accuracy of those services are highly dependent on the navigation performance of the relay itself. To assess the Lunar relay navigation performance, a series of orbit determination (OD) Monte Carlo (MC) simulations are run using Lunar gravity modeling up to a degree and order of 250 and a variety of onboard clocks and measurement types including weak-signal GNSS, Ground Network (GN) pseudorange (PR) and Doppler, and optical navigation (OpNav) center-finding (CF). The estimated trajectories produced by these Lunar relay MC simulations, along with the associated errors, and transmitted navigation reference signal parameters, are used to evaluate the expected navigation performance of a user on a descent trajectory to the Lunar surface. The scenario features a lander system performing onboard navigation relying on one-way range and Doppler measurements from reference signals emitted by the Lunar relay. This paper can be used as a reference in determining the onboard clock and measurement types necessary to obtain acceptable navigation performance for the Lunar relay, and as a baseline for Lunar lander navigation performance using accurate measurements from a relay reference signal.

Jeffrey L. Small↗
NASA NTRS2021

Variability of Disk Emission in Pre-main Sequence and Related Stars. V. OccultationEvents from the Innermost Disk Region of the Herbig Ae Star HD 163296

HD 163296 is a Herbig Ae star that underwent a dramatic∼0.8 magnitude drop in brightness in the V photometric band in 2001 and a brightening in the near-IR in 2002. Because the star possesses Herbig–Haro objects traveling in out flowing bipolar jets, it was suggested that the drop in brightness was due to a clump of dust entrained in a disk wind, blocking the line of sight toward the star. In order to quantify this hypothesis, we investigated the brightness drop at visible wavelengths and the brightening at near-IR wavelengths of HD 163296 using the Monte Carlo Radiative Transfer Code,HOCHUNK3D. We created three models to understand the events. Model 1 describes the quiescent state of the system. Model 2 describes the change in structure that led to the drop in brightness in 2001.Model 3 describes the structure needed to produce the observed 2002 brightening of the near-IR wavelengths. Models 2 and 3 utilize a combination of a disk wind and central bipolar flow. By introducing a filled bipolar cavity in Models 2 and 3, we were able to successfully simulate a jet-like structure for the star with a disk wind and created the drop and subsequent increase in brightness of the system. On the other hand, when the bipolar cavity is not filled, Model 1 replicates the quiescent state of the system

Monika Pikhartova↗
NASA NTRSDate not supplied

Lunar Lander Landing Modeling and Simulation in Simscape

The landing gear is a critical subsystem of a lander for a successful landing on the moon or a planet. During a landing event, structural loading, energy absorption and toppling stability are important factors for the landing gear development. These aspects can be studied using simulations. The landing gear comprising the primary struct, secondary structs, footpad, and joints, is easily modeled using Simscape Multibody toolbox. The energy absorption mechanisms in the structs is readily modeled using Simulink and Multibody blocks based on the relative motions of the inner and outer cylinder of the structures. During touchdown, the footpads strike the lunar surface, and the contact forces can be modeled including the soil mechanical properties. Slosh dynamics are included using a mechanical pendulum model. Tip-over controllers are developed using Simulink blocks and toolboxes. This model was intensively tested manually then additionally with Monte Carlo simulation. It offers the opportunity for landing stability study, i.e. the max slope a lander can land without tipping over in 3D realistic landing situations.

Wei Wu↗
DOE OSTI2024

Data-Driven Computation of Probabilistic Marching Cubes for Efficient Visualization of Level-Set Uncertainty

Uncertainty visualization is an important emerging research area. Being able to visualize data uncertainty can help scientists improve trust in analysis and decision-making. However, visualizing uncertainty can add computational overhead, which can hinder the efficiency of analysis. In this paper, we propose novel data-driven techniques to reduce the computational requirements of the probabilistic marching cubes (PMC) algorithm. PMC is an uncertainty visualization technique that studies how uncertainty in data affects level-set positions. However, the algorithm relies on expensive Monte Carlo (MC) sampling for the multivariate Gaussian uncertainty model because no closed-form solution exists for the integration of multivariate Gaussian. In this work, we propose the eigenvalue decomposition and adaptive probability model techniques that reduce the amount of MC sampling in the original PMC algorithm and hence speed up the computations. Our proposed methods produce results that show negligible differences compared with the original PMC algorithm demonstrated through metrics, including root mean squared error, maximum error, and difference images. We demonstrate the performance and accuracy evaluations of our data-driven methods through experiments on synthetic and real datasets.

Athawale, Tushar↗
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 OSTI2024

Posterior Covariance Matrix Approximations

Here, the Davis equation of state (EOS) is commonly used to model thermodynamic relationships for high explosive (HE) reactants. Typically, the parameters in the EOS are calibrated, with uncertainty, using a Bayesian framework and Markov Chain Monte Carlo (MCMC) methods. However, MCMC methods are computationally expensive, especially for complex models with many parameters. This paper provides a comparison between MCMC and less computationally expensive Variational methods (Variational Bayesian and Hessian Variational Bayesian) for computing the posterior distribution and approximating the posterior covariance matrix based on heterogeneous experimental data. All three methods recover similar posterior distributions and posterior covariance matrices. This study demonstrates that for this EOS parameter calibration application, the assumptions made in the two Variational methods significantly reduce the computational cost but do not substantially change the results compared to MCMC.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Comparisons of the Maxwell and CLL Gas/Surface Interaction Models Using DSMC

Two contrasting models of gas-surface interactions are studied using the Direct Simulation Monte Carlo (DSMC) method. The DSMC calculations examine differences in predictions of aerodynamic forces and heat transfer between the Maxwell and Cercignani-Lampis-Lord (CLL) models for flat plate configurations at freestream conditions corresponding to a 140 km orbit around Venus. The size of the flat plate is that of one of the solar panels on the Magellan spacecraft, and the freestream conditions are one of those experienced during aerobraking maneuvers. Results are presented for both a single flat plate and a two-plate configuration as a function of angle of attack and gas-surface accommodation coefficients. The two plate system is not representative of the Magellan geometry, but is studied to explore possible experiments that might be used to differentiate between the two gas surface interaction models.

Hedahl, Marc O.↗
NASA NTRS1993

Investigation of radiative interaction in laminar flows using Monte Carlo simulation

The Monte Carlo method (MCM) is employed to study the radiative interactions in fully developed laminar flow between two parallel plates. Taking advantage of the characteristics of easy mathematical treatment of the MCM, a general numerical procedure is developed for nongray radiative interaction. The nongray model is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. To validate the Monte Carlo simulation for nongray radiation problems, the results of radiative dissipation from the MCM are compared with two available solutions for a given temperature profile between two plates. After this validation, the MCM is employed to solve the present physical problem and results for the bulk temperature are compared with available solutions. In general, good agreement is noted and reasons for some discrepancies in certain ranges of parameters are explained.

Liu, Jiwen↗
NASA NTRSDate not supplied

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗
NASA NTRS1995

Monte Carlo Simulation of the Rapid Crystallization of Bismuth-Doped Silicon

In this Letter we report Ising model simulations of the growth of alloys which predict quite different behavior near and far from equilibrium. Our simulations reproduce the phenomenon which has been termed 'solute trapping,' where concentrations of solute, which are far in excess of the equilibrium concentrations, are observed in the crystal after rapid crystallization. This phenomenon plays an important role in many processes which involve first order phase changes which take place under conditions far from equilibrium. The underlying physical basis for it has not been understood, but these Monte Carlo simulations provide a powerful means for investigating it.

Jackson, Kenneth A.↗
DOE OSTI2020

Learning Hidden Structure in Multi-Fidelity Information Sources for Efficient Uncertainty Quantification (LDRD 218317)

This report summarizes the work done under the Laboratory Directed Research and Development (LDRD) project entitled "Learning Hidden Structure in Multi-Fidelity Information Sources for Efficient Uncertainty Quantification". In this project we investigated multi-fidelity strategies for fusing data from information sources of varying cost and accuracy. Most existing strategies exploit hierarchical relationships between models, for example that occur when different models are generated by refining a numerical discretization parameter. In this work we focused on encoding the relationships between information sources using directed acyclic graphs. The multi-fidelity networks can have general structure and represent a significantly greater variety of modeling relationships than recursive networks used in the current state literature. Numerical results show that a non-hierarchical multi-fidelity Monte Carlo strategy can reduce the cost of estimating uncertainty in predictions of a model of plasma expanding in a vacuum by almost two orders of magnitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1997

CGRO Guest Investigator Program

The following are highlights from the research supported by this grant: (1) Theory of gamma-ray blazars: We studied the theory of gamma-ray blazars, being among the first investigators to propose that the GeV emission arises from Comptonization of diffuse radiation surrounding the jet, rather than from the synchrotron-self-Compton mechanism. In related work, we uncovered possible connections between the mechanisms of gamma-ray blazars and those of intraday radio variability, and have conducted a general study of the role of Compton radiation drag on the dynamics of relativistic jets. (2) A Nonlinear Monte Carlo code for gamma-ray spectrum formation: We developed, tested, and applied the first Nonlinear Monte Carlo (NLMC) code for simulating gamma-ray production and transfer under much more general (and realistic) conditions than are accessible with other techniques. The present version of the code is designed to simulate conditions thought to be present in active galactic nuclei and certain types of X-ray binaries, and includes the physics needed to model thermal and nonthermal electron-positron pair cascades. Unlike traditional Monte-Carlo techniques, our method can accurately handle highly non-linear systems in which the radiation and particle backgrounds must be determined self-consistently and in which the particle energies span many orders of magnitude. Unlike models based on kinetic equations, our code can handle arbitrary source geometries and relativistic kinematic effects In its first important application following testing, we showed that popular semi-analytic accretion disk corona models for Seyfert spectra are seriously in error, and demonstrated how the spectra can be simulated if the disk is sparsely covered by localized 'flares'.

Begelman, Mitchell C.↗
NASA NTRS1994

Probability density function approach for compressible turbulent reacting flows

The objective of the present work is to extend the probability density function (PDF) tubulence model to compressible reacting flows. The proability density function of the species mass fractions and enthalpy are obtained by solving a PDF evolution equation using a Monte Carlo scheme. The PDF solution procedure is coupled with a compression finite-volume flow solver which provides the velocity and pressure fields. A modeled PDF equation for compressible flows, capable of treating flows with shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed. Two super sonic diffusion flames are studied using the proposed PDF model and the results are compared with experimental data; marked improvements over solutions without PDF are observed.

Hsu, A. T.↗
DOE OSTI2020

Multilevel and multifidelity uncertainty quantification for cardiovascular hemodynamics

Standard approaches for uncertainty quantification in cardiovascular modeling pose challenges due to the large number of uncertain inputs and the significant computational cost of realistic three-dimensional simulations. Here, we propose an efficient uncertainty quantification framework utilizing a multilevel multifidelity Monte Carlo (MLMF) estimator to improve the accuracy of hemodynamic quantities of interest while maintaining reasonable computational cost. This is achieved by leveraging three cardiovascular model fidelities, each with varying spatial resolution to rigorously quantify the variability in hemodynamic outputs. We employ two low-fidelity models (zero- and one-dimensional) to construct several different estimators. Our goal is to investigate and compare the efficiency of estimators built from combinations of these two low-fidelity model alternatives and our high-fidelity three-dimensional models. We demonstrate this framework on healthy and diseased models of aortic and coronary anatomy, including uncertainties in material property and boundary condition parameters. Our goal is to demonstrate that for this application it is possible to accelerate the convergence of the estimators by utilizing a MLMF paradigm. Therefore, we compare our approach to single fidelity Monte Carlo estimators and to a multilevel Monte Carlo approach based only on three-dimensional simulations, but leveraging multiple spatial resolutions. We demonstrate significant, on the order of 10 to 100 times, reduction in total computational cost with the MLMF estimators. We also examine the differing properties of the MLMF estimators in healthy versus diseased models, as well as global versus local quantities of interest. As expected, global quantities such as outlet pressure and flow show larger reductions than local quantities, such as those relating to wall shear stress, as the latter rely more heavily on the highest fidelity model evaluations. Similarly, healthy models show larger reductions than diseased models. In all cases, our workflow coupling Dakota’s MLMF estimators with the SimVascular cardiovascular modeling framework makes uncertainty quantification feasible for constrained computational budgets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Enforcing detailed balance in the Borgnakke–Larsen redistribution method with temperature dependent relaxation models

For decades, it has been observed that the commonly used Borgnakke–Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke–Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke–Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke–Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). Here, the consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Calculation of Groundwater Pathway Radiological Dose for the Hanford Site Composite Analysis Null-Space Model Carlo Flow Model Set

The purpose of this environmental calculation file (ECF) is to present the results of the exposure route-specific and total radiological dose assessments for the groundwater pathway based on the null space Monte Carlo (NSMC) groundwater concentrations as a part an uncertainty analysis for the updated Hanford Site Composite Analysis (CA). The Plateau-to River (P2R) Groundwater Model (CP-57037, Model Package Report: Plateau to River Groundwater Model Version 8.3) is the CA base case that simulates the fate and transport of radiological contaminants within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. An NSMC analysis was performed to identify and quantify the potential uncertainties associated with the P2R Model. The result of the NSMC analysis is a set of flow and transport simulations that provide an estimated range of possible outcomes that are used to quantify the uncertainty associated with the simulated base case concentrations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

Predicting the performance of solar illuminated thermal louvers for the GOES I, J, K instruments

The new series of Geostationary Operational Environmental Satellites (GOES I,J,K,L/M) will have two earth viewing optical instruments which use passive, bi-metallic actuated, blade-type louvers for their thermal control. GOES is a three-axis stabilized spacecraft in geostationary orbit, which means that sun can enter the instrument apertures for several hours before and after local midnight. The solar heat absorbed in the instruments causes the louvers to open thus allowing the radiators to reject the heat. These louvered radiators are located on the north face of the instruments. To perform thermal analyses of the instruments, a modeling method was needed for predicting the thermal characteristics of the louvers. This task was complicated by the fact that sun shines on the louvers during the summer season. A Monte Carlo analytical technique was used to develop absorbed power and effective emissivity tables for the louvers. These tables could be used with the instrument math models to predict the thermal behavior of the instruments. Data from the Monte Carlo analysis showed that solar entrapment could have a significant effect on the heat rejection capability of the louvered radiators.

Mengers, David↗