Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation”

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,495 records · Page 83

Coupling reduced‐order modeling and coarse‐grained CFD‐DEM to accelerate coal gasifier simulation and optimization

Abstract High fidelity three‐dimensional (3‐D) numerical simulations of commercial‐scale coal gasifiers are very time consuming and expensive. This article proposes a reduced‐order modeling approach that uses quasi one‐dimensional (1‐D) CFD–DEM simulation results to serve as an accurate initial condition for the 3‐D simulation, which significantly shortens the physical time from several hours to minutes to achieve steady states. The 3‐D simulation employs a previously validated coarse‐grained CFD‐DEM approach, which further accelerates the simulation. Both 1‐D and 3‐D simulations adopt a particle shrinkage scheme to correctly account for the particle volume change due to chemical reaction and thus provides an accurate bed height. Comparison with available experimental measurements of the bed temperature and synthesis gas composition is used to validate the simulation. The final syngas composition and flowrates are strongly affected by the fuel or gasification agent rate and process operating conditions as expected. This approach is then used to optimize two important operating conditions (air and steam flows) to achieve specific gasifier performance goals for the production of Fischer–Tropsch synthesis gas and production of fuel gas for gas turbines.

Yu, Jia↗

The impact of coupled air–sea interaction on extreme East Asian summer monsoon simulation in CMIP5 models

Abstract In this study, the relationship between the ability to simulate air–sea interactions over the western North Pacific (WNP), and to reproduce the extreme East Asian summer monsoon (EASM), were investigated by comparing the performances of several global climate models (GCMs). High ranked in air–sea interaction simulation (HRA) and low ranked in air–sea interaction simulation (LRA) models were selected, according to their performance in simulating relations between sea surface temperature (SST) and precipitation over the WNP, from the ensemble of models that participated in the third and fifth phases of the Coupled Model Intercomparison Project (CMIP3, CMIP5). Compared with CMIP3 models, CMIP5 models exhibited improved simulations of the distinctive air–sea interaction over the WNP, namely, the strong atmospheric forcing on the ocean. Among CMIP5 models, HRA models, which reproduced intrinsic negative correlations between precipitation and SST over the WNP, could simulate the extreme EASM better than LRA models. In particular, HRA models generated a more realistic spatial distribution of the extreme EASM compared with LRA models. The defects of the LRA models resulted from distorted synoptic fields, including underestimated geopotential height and overestimated low‐level wind over the WNP, inducing unrealistic moisture supply and convection due to the exaggerated SST forcing. In contrast, reasonable air–sea interactions represented in HRA models lead to realistic synoptic fields over the WNP, and proper simulation of the extreme EASM.

Kim, Taehyung↗

Efficient Simulation of Open Quantum Systems on NISQ Trapped‐Ion Hardware

Abstract Simulating open quantum systems, which interact with external environments, presents significant challenges on noisy intermediate‐scale quantum (NISQ) devices due to limited qubit resources and noise. In this study, an efficient framework is proposed for simulating open quantum systems on NISQ hardware by leveraging a time‐perturbative Kraus operator representation of the system's dynamics. This approach avoids the computationally expensive Trotterization method and exploits the Lindblad master equation to represent time evolution in a compact form, particularly for systems satisfying specific commutation relations. The efficiency of this method is demonstrated by simulating quantum channels, such as the continuous‐time Pauli channel and damped harmonic oscillators, on NISQ trapped‐ion hardware, including IonQ Harmony and Quantinuum H1‐1. Additionally, hardware‐agnostic error mitigation techniques are introduced, including Pauli channel fitting and quantum depolarizing channel inversion, to enhance the fidelity of quantum simulations. These results show strong agreement between the simulations on real quantum hardware and exact solutions, highlighting the potential of Kraus‐based methods for scalable and accurate simulation of open quantum systems on NISQ devices. This framework opens pathways for simulating more complex systems under realistic conditions in the near term.

Burdine, Colin [Department of Electrical and Compu↗

Using AI Simulations to Dynamically Model Multi-agent Multi-team Energy Systems

The complexity of energy systems is well known as they are complex and intricate systems. As a result, many extant studies have used many simplifications or generalizations that do not accurately reflect the nature of this complex system. In particular, most HVAC systems are modeled as a single unit, or several large units, rather than as a hierarchical composite (e.g., as a floor rather than as a collection of disparate rooms). The net result of this is that the simulations are too generic to perform meaningful analysis, machine learning, or integrated simulation. We propose using a multi-agent multi-team strategic simulations framework called SiMAMT to better define, model, simulate, and learn the HVAC environment. SiMAMT allows us to create distinct models for each type of room, hierarchically aggregate them into units (like floors, or sections), and then into larger sets (like buildings or a campus), and then perform a simulation that interacts with each sub-element individually, the teams of sub-elements collectively, and the entire set in aggregation. Further, and most importantly, we additionally model another ‘team’ within the simulation framework - the users of the systems. Again, each individual is modeled distinctly, aggregated into sub-sets, then collected into large sets. Each user, or agent, is performing on their own but with respect to the larger team goals. This provides a simulation that has a much higher model fidelity and more applicable results that match the real-world.

Franklin, D. Michael↗

The Current State of High-Fidelity Simulations for Main Gas Path Turbomachinery Components and Their Industrial Impact

High-fidelity simulations have become feasible for most main gas path turbomachinery components, over the past two decades. This paper introduces the key challenges of simulating and modelling turbomachinery flows and presents an overview of possible simulation strategies. A comprehensive overview is given of which particular design challenges have to date been addressed by high-fidelity simulations, with a particular focus on high-pressure and centrifugal compressors, and high-pressure and low-pressure turbines. Recommendations are provided for how quality and accuracy can be ensured for high-fidelity simulations, using direct numerical simulations of a low-pressure turbine as a case study. It is argued that industrial impact from high-fidelity simulations can be achieved in two ways, either by conducting design-of-experiment-like studies that can provide designers with insight into certain physical mechanisms and phenomena, or by helping mature and improve lower-order models. The latter is discussed with particular emphasis on recent advances in machine learning for model assessment and improvement, and the potential of one selected data-driven approach is demonstrated on predictions of wake mixing for low-pressure and high-pressure turbines.

Mechanics↗

Neutrino transport in general relativistic neutron star merger simulations

Abstract Numerical simulations of neutron star–neutron star and neutron star–black hole binaries play an important role in our ability to model gravitational-wave and electromagnetic signals powered by these systems. These simulations have to take into account a wide range of physical processes including general relativity, magnetohydrodynamics, and neutrino radiation transport. The latter is particularly important in order to understand the properties of the matter ejected by many mergers, the optical/infrared signals powered by nuclear reactions in the ejecta, and the contribution of that ejecta to astrophysical nucleosynthesis. However, accurate evolutions of the neutrino transport equations that include all relevant physical processes remain beyond our current reach. In this review, I will discuss the current state of neutrino modeling in general relativistic simulations of neutron star mergers and of their post-merger remnants. I will focus on the three main types of algorithms used in simulations so far: leakage, moments, and Monte-Carlo scheme. I will review the advantages and limitations of each scheme, as well as the various neutrino–matter interactions that should be included in simulations. We will see that the quality of the treatment of neutrinos in merger simulations has greatly increased over the last decade, but also that many potentially important interactions remain difficult to take into account in simulations (pair annihilation, oscillations, inelastic scattering).

Foucart, Francois (ORCID:0000000346174738)↗

Numerical simulation of involute-plate research reactor flow behavior using RANS, LES and DNS

This paper investigates the flow behavior of involute-plate research reactors by performing Reynolds-Averaged Navier Stokes simulation (RANS), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) of the channel flow between fuel plates. By modeling turbulence with different numerical approaches, this study provides data with three levels of fidelity. For the RANS simulation, three widely used turbulence models, i.e., k-ε, k-ω, Reynolds Stress Turbulence model (RST) are applied by using the commercial CFD code STAR-CCM +. For LES and DNS, the open-source CFD code, Nek5000, is used given its outstanding scalability on High Performance Computer (HPC) and high-order technique. The results from RANS simulations are compared with that from LES and DNS for benchmarking. Both macroscale parameters and turbulence statistics, such as velocity magnitude, lateral velocity and turbulence kinetic energy, are presented and analyzed. The results from RANS simulation achieve good agreement with LES and DNS on velocity and turbulence kinetic energy prediction. The RST turbulence model predicts the most similar flow pattern of lateral velocity as compared to LES and DNS. The Lambda-2 (λ2) criterion with a reasonable threshold is used to demonstrate the instantaneous vortices distribution in the involute channel from both LES and DNS calculation. The DNS simulation captures more detailed turbulence especially near the corner, which explains the discrepancy between LES and DNS results near the corner. The normalized RMS error are defined and calculated to assess the performance of those turbulence models. The RST model captures the anisotropic feature of turbulence, which enable it to outperform other turbulence models for predicting the flow behavior in an involute channel. Although some discrepancies are found between LES and DNS results in the corner, the overall deviations between LES and DNS are found to be small. In conclusion, given that the computational cost of DNS calculation is an order of magnitude higher, using LES data for benchmarking RANS model is a cost-effective approach.

DNS↗

Stochastic simulation of occupant-driven energy use in a bottom-up residential building stock model

The residential buildings sector is one of the largest electricity consumers worldwide and contributes disproportionally to peak electricity demand in many regions. Strongly driven by occupant activities, household energy consumption is stochastic and heterogeneous in nature. However, most residential energy models applied by industry use homogeneous, deterministic activity schedules, which work well for predictions of annual energy consumption, but can result in unrealistic hourly or sub-hourly electric load profiles, with exaggerated or muted peaks. The increasing proportion of variable renewable energy generators means that representing the heterogeneity and stochasticity of occupant behavior is now crucial for reliable planning at both bulk-power and distribution-system scales. This work presents a novel and open-source occupancy simulation approach that can simulate a diverse set of individual occupant and household event schedules for all major electricity, fuel, and hot water end uses. To accomplish this, we evaluated three alternative occupant activity simulation approaches before selecting a hybrid combining time-inhomogeneous Markov chains and probability-sampling of event durations and magnitudes. Further, we integrated the stochastic occupancy simulation with an open-source bottom-up physics-simulation building stock model and published a set of 550,000 diverse household end-use activity schedules representing a national housing stock. The simulator was verified against time-use survey data, and simulation results were validated against measured end-use electricity data for accuracy and reliability. While we use data for the United States, our application demonstrates how similar approaches could be applied using the time-use survey data collected in many countries around the world.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Dynamics and lipid membrane coupling of the RAS-RAF complex revealed via multiscale simulations

To gain molecular and mechanistic insights into initiation of the RAS-RAF signaling cascade, we developed and used a combination of multiscale simulation and experimental approaches. The influence and impact of the membrane on RAS and RAF proteins is a factor we are just beginning to understand and appreciate in more detail. Molecular simulation is an ideal methodology to further study this complicated relationship between the membrane and associated proteins. Our previous work using Multiscale Machine-learned Modeling Infrastructure investigated different lipid compositions solely around the KRAS4b protein and the interplay between protein behavior and these membrane environments. Multiscale Machine-learned Modeling Infrastructure uses machine learning to couple adjacent simulation scales and has been efficiently scaled across some of the world’s largest high-performance computers. Recently, we have expanded this multiresolution framework to include the all-atom simulation scale and to incorporate the RAF RBDCRD domains. Here, we present the overall analysis results from this new simulation campaign comprising a mixture of RAS and RAF RBDCRD proteins. Approximately 35,000 coarse-grained and 10,000 all-atom molecular dynamics simulations were completed, sampled from a variety of protein/lipid composition configurations that were generated from a micron-scale continuum simulation containing hundreds of copies of the proteins. Our studies suggest that orientations of the RAS-RBDCRD complex on the membrane occupy distinct configurational states, and the spatial patterns of lipid arrangements around these different protein states are unique to each state. The extent and size of lipid “fingerprints” imposed on the membrane by the RAS-RBDCRD protein complex are significantly larger than observed for just the RAS protein on its own. These protein complexes strongly associate, but we do not observe statistically significant preferred protein-protein orientations. These observations indicate that spatial colocalization of RAS-RBDCRD proteins in the same vicinity may be assisted by specific membrane environments, acting to increase the probability of signaling complex formation.

Carpenter, Timothy S. [Lawrence Livermore National↗

A simulation approach to estimate energy savings potential of occupant behavior measures

Occupant behavior in buildings is a leading factor influencing energy use in buildings. Low-cost behavioral solutions have demonstrated significant potential energy savings. Estimating the behavioral savings potential is important for a more effective design of behavior change interventions, which in turn will support more effective energy-efficiency policies. This study introduces a simulation approach to estimate the energy savings potential of occupant behavior measures. First it defines five typical occupant behavior measures in office buildings, then simulates and analyzes their individual and integrated impact on energy use in buildings. The energy performance of the five behavior measures was evaluated using EnergyPlus simulation for a real office building across four typical U.S. climates and two vintages. The Occupancy Simulator was used to simulate the occupant movement in each zone with inputs from the site survey of the case building. Based on the simulation results, the occupant behavior measures can achieve overall site energy savings as high as 22.9% for individual measures and up to 41.0% for integrated measures. Although energy savings of behavior measures would vary depending upon many factors, the presented simulation approach is robust and can be adopted for other studies aiming to quantify occupant behavior impact on building performance.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

From models to reality: a systematic review on simulated and measured residential heat pump energy savings

High-performance HVAC solutions are central to residential energy management. A substantial share of these are electric, reversible-cycle systems, with heat pumps representing the largest portion of current and near-term adoption. This review synthesizes peer-reviewed and grey literature on residential space heating and cooling heat pumps. The academic literature is dominated by modeling (73.8%), with limited field measurement (13.1%). Grey literature from United States serve as a supplemental resource providing measured savings. Conversions from electric-resistance heating consistently show the largest site energy reductions, while oil/propane baselines yield moderate savings, and gas baseline scenario often deliver small and region-dependent savings. This study cross-checks the grey literature measured data with simulation data filtered from the ResStock dataset. The comparison indicates a discrepancy between simulations and measured data: simulated site EUIs are typically lower than measured EUIs, but percentage energy savings fall in similar ranges, implying simulations capture directional effects while underestimating energy use. Factors associated with variability and model–measurement differences include system characterization and control representation (e.g., backup heat engagement, thermostat/setpoint strategies, commissioning/installation quality), occupant behavior, weather normalization, metering scope, and envelope characterization. This paper also outlines the proposed methodology for comparing simulation and measured data for heat pumps. It emphasizes the metrics used for comparison and units harmonization, building characteristics matching, and compact metadata are needed for simulations to match measured data. The proposed methodology is expected to improve the credibility of simulated savings as measured evidence grows.

Yu, Lili↗

Enabling co-simulation of tokamak plant models and plasma control systems

A system for connecting the Plasma Control System and a model of the tokamak Plant in closed loop co-simulation for plasma control development has been in routine use at DIII-D for more than 20 years and at other fusion labs that use variants of the DIII-D PCS for approximately the last decade. Here, co-simulation refers to the simultaneous execution of two independent codes with the exchange of data - Plant actuator commands and tokamak diagnostic data - between them during execution. Interest in this type of PCS-Plant simulation technology has also been growing recently at other fusion facilities. In fact, use of such closed loop control simulations is assumed to play an even larger role in the development of both the ITER Plasma Control System (PCS) and the experimental operation of the ITER device, where they will be used to support verification/validation of the PCS and also for ITER pulse schedule development and validation. We describe the key use cases that motivate the co-simulation capability and the features that must be provided by the Plasma Control System to support it. These features could be provided by the PCS itself or by a model of the PCS. If the PCS itself is chosen to provide them, there are requirements imposed on its architecture. If a PCS model is chosen, there are requirements imposed on the initial implementation of this simulation as well as long-term consequences for its continued development and maintenance. We describe these issues for each use case and discuss the relative merits of the two choices. Several examples are given illustrating uses of the co-simulation method to address problems of plasma control during the operation of DIII-D and of other devices that use the DIII-D PCS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sensitivity of the simulation of passive neutron emission from UF 6 cylinders to the uncertainties in both 19 F(α,n) energy spectrum and thick target yield of 234 U in UF 6

Interest in safeguards verification measurements using passive thermal neutron counting to assay 235 U content in large 30B UF 6 canisters has grown in recent years. Here, the prohibitively high cost and impracticality of using reference 30B calibration cylinders extensively will likely make accurate simulations of increasing interest. Accuracy of the simulated response will define the confidence in the predicted response and the extent to which simulations can reasonably be relied upon. With 234 U driven 19 F(α, n) reactions being the main neutron source in low enriched UF 6 the uncertainties of the 19 F(α, n) energy spectrum and the thick target yield of 234 U in UF 6 propagate into the uncertainty in the predicted response and represent a major influence of basic nuclear data. Here sensitivity of the simulated total (Singles) and coincidence (Doubles) count rates are assessed for the Passive Neutron Enrichment Meter using six potential 19 F(α, n) neutron energy spectra over a range of enrichments and material distributions. The results indicate that variations in the Singles and Doubles due to simulated (α, n) neutron spectrum are less than 1.5% for this set of simulated neutron spectra, with dependence varying inversely with enrichment. Singles uncertainty is only slightly less than that of the thick target 19 F(α, n) yield corresponding to the primary neutron source, whereas the 19 F(α, n) yield dependence of the Doubles is reduced by the non-negligible 238 U(SF) coincident neutron emissions. Based on available thick target 19 F(α, n) yield estimates the uncertainty is on the order of 5%, establishing this as the main nuclear data limitation when simulating thermal neutron detectors response for 30B UF 6 storage cylinders. Based on these findings, it appears that the measurement and evaluation of the thick target 19 F(α, n) yield for uranium hexafluoride is due.

19F(α,n) neutron spectrum↗

Automation and optimization of stopping and range of ions in matter simulation runtime

Prior to every ion implantation experiment a simulation of the ion range and other relevant parameters is performed using Monte-Carlo based codes. Although increasing computing power has improved the speed of these calculations, the demands on Monte-Carlo codes are also increasing, requiring evaluation of the optimal number of simulations while ensuring accuracy within threshold bounds. We evaluate the “Stopping and Range of Ions in Matter” (SRIM) code due to its widespread usage. We show how dividing simulations into multiple parallel simulations with different random seeds can lead to calculation speedup and find lower bounds for the required number of ion traces simulated based on an exemplar system of a Ga focused ion beam and a high energy C beam as used in high linear energy transfer testing. Here our results indicate simulations can yield results within the underlying data accuracy of SRIM at 10X and 100X shorter simulation time than the SRIM default values.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

One-way coupled simulation of FIV in a 7-pin wire-wrapped fuel pin bundle

A normal fuel geometry found in sodium cooled reactor designs is the hexagonal lattice of wire-wrapped fuel pins. The wire-wrappers, which replace the more common spacer grids, act as a spacer between the fuel pins, while simultaneously increasing the fluid mixing between subchannels. The induced transverse flow, while good for mixing, increases the force load experienced by the fuel pins, increasing the possibility of problematic flow induced vibrations. The primary source of vibrations is turbulent buffeting, so very small in amplitude, but capable of rupturing the fuel cladding by fretting. The vast majority of research on fretting has been for bare rods with spacer grids. There are several significant ways that wire-wrapped fuel pins differ from bare rods that may have significant impact on the vibrations present. The increased transverse flow and asymmetrical flow pattern may result in larger and imbalanced force loads on the pins. The helical bracing of each pin may result in unique vibration patterns, and gaps may form between pins affecting how they are braced. Simulation allows for the observation of the forces present, how the flow field affects these forces, and how the structure responds and can be repeated for multiple cases for relatively low cost. Here, simulation of the flow induced vibrations for a 7-pin wire-wrapped fuel pin bundle is carried out by coupling Nek5000, a spectral element LES/DNS solver with DIABLO, a finite element CSM solver. The small amplitude vibrations allow for the assumption that the structure motion has negligible effect on the fluid flow, permitting the much less computationally expensive one-direction coupling. In order to gather significant force history data for the 2.63 helical pitch structure, a scheme was developed to artificially create force history from the power spectrum density of the one-pitch CFD simulation forces. The vibrations could then be simulated in the structure multiple times varying the contact and gaps between various pins. This provides for the search of worst-case scenarios and informs about the types of fretting that might be present. The simulation geometry is based on an experiment for comparison. However, the viscosity of the simulation was increased by a factor of four to make the simulation computationally feasible. While validation is not possible, comparisons are still valuable.

42 ENGINEERING↗

Comparison of volatiles evolving from selected highland and mare lunar regolith simulants during vacuum sintering

Volatiles evolving from JSC-1A, NU-LHT-4M and CSM-LHT-1G lunar regolith simulants during in vacuo thermal processing were analyzed using mass spectrometry as a function of temperature. Two high-fidelity simulants, JSC-1A (mare) and NU-LHT-4M (highland), were compared to a newly developed CSM-LHT-1G highland simulant, modified to closely match lunar geochemistry. Large autogenous gas loads were observed for all investigated materials. Mineralogical knowledge was used to identify and attribute individual volatile species to reacting, transforming, or decomposing constituents (hydrates, carbonates, sulfates, sulfides, clays, etc.) of the respective regolith simulant in the self-generated gas environment. Cumulative mass losses for individual simulant components as a function of temperature were quantified using mass spectrometry in conjunction with thermogravimetric analysis. Investigation of the four components of CSM-LHT-1G – anorthosite, basalt, augite, and glass – aided the attribution of volatile species to specific compounds and their respective sources. The results showed significant decomposition of non-lunar phases present in the man-made regolith simulants below the typical glass crystallization temperatures, which paves the way to devising methods for enhancing the fidelity of the simulants. Finally, high gas loads and corrosive gases (HF and HCl) were recognized as potential hazards, pertaining to the development of large testbed facilities.

42 ENGINEERING↗

Discrepancy quantification between experimental and simulated data of CO 2 adsorption isotherm using hierarchical Bayesian estimation

Here, to quantitatively analyze the inconsistencies commonly observed between experimental and simulated adsorption isotherms, parameter estimation of adsorption isotherm models was conducted by hierarchical Bayesian estimation with parameter uncertainties being quantified as probability distributions. The estimation method was implemented using Markov Chain Monte Carlo (MCMC) to analyze multiple data sets obtained from different sources, including a publicly available database. To describe the discrepancies of experimental and simulated adsorption data, the simulation data was set as the reference to which experimental measurements were compared. We applied the proposed approach to analyze CO 2 adsorption isotherms that are measured and simulated on zeolite 13X and MIL-101(Cr). In these case studies, the discrepancy of CO 2 adsorption isotherm was successfully quantified between experimental measurements and predictions given by molecular simulations using Grand Canonical Monte Carlo (GCMC), where uncertainties were quantified as probability distributions. Furthermore, experimental data sets that agree well with the GCMC simulation have been identified, providing insights into experimental and measurement methods as well as choosing the right assumptions in the molecular simulation.

42 ENGINEERING↗

$\mathrm{SIM}$-$\mathrm{STEM}$ $\mathrm{Lab}$: Incorporating Compressed Sensing Theory for Fast STEM Simulation

We report recently it has been shown that precise dose control and an increase in the overall acquisition speed of atomic resolution scanning transmission electron microscope (STEM) images can be achieved by acquiring only a small fraction of the pixels in the image experimentally and then reconstructing the full image using an inpainting algorithm. In this paper, we apply the same inpainting approach (a form of compressed sensing) to simulated, sub-sampled atomic resolution STEM images. We find that it is possible to significantly sub-sample the area that is simulated, the number of g-vectors contributing the image, and the number of frozen phonon configurations contributing to the final image while still producing an acceptable fit to a fully sampled simulation. Here we discuss the parameters that we use and how the resulting simulations can be quantifiably compared to the full simulations. As with any Compressed Sensing methodology, care must be taken to ensure that isolated events are not excluded from the process, but the observed increase in simulation speed provides significant opportunities for real time simulations, image classification and analytics to be performed as a supplement to experiments on a microscope to be developed in the future

47 OTHER INSTRUMENTATION↗