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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A general CFD framework for fault-resilient simulations based on multi-resolution information fusion

We develop a general CFD framework for multi-resolution simulations to target multiscale problems but also resilience in exascale simulations, where faulty processors may lead to gappy, in space-time, simulated fields. We combine approximation theory and domain decomposition together with statistical learning techniques, e.g. coKriging, to estimate boundary conditions and minimize communications by performing independent parallel runs. To demonstrate this new simulation approach, we consider two benchmark problems. First, we solve the heat equation (a) on a small number of spatial “patches” distributed across the domain, simulated by finite differences at fine resolution and (b) on the entire domain simulated at very low resolution, thus fusing multi-resolution models to obtain the final answer. Second, we simulate the flow in a lid-driven cavity in an analogous fashion, by fusing finite difference solutions obtained with fine and low resolution assuming gappy data sets. We investigate the influence of various parameters for this framework, including the correlation kernel, the size of a buffer employed in estimating boundary conditions, the coarseness of the resolution of auxiliary data, and the communication frequency across different patches in fusing the information at different resolution levels. In addition to its robustness and resilience, the new framework can be employed to generalize previous multiscale approaches involving heterogeneous discretizations or even fundamentally different flow descriptions, e.g. in continuum-atomistic simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Flow-induced vibration analysis of a helical coil steam generator experiment using large eddy simulation

This paper describes a numerical study of flow-induced vibration in a helical coil steam generator experiment conducted at Argonne National Laboratory in the 1980s. In the experiment, a half-scale sector model of a steam generator helical coil tube bank was subjected to still and flowing air and water, and the vibrational characteristics were recorded. The research detailed in this document utilizes the multi-physics simulation toolkit SHARP developed at Argonne National Laboratory, in cooperation with Lawrence Livermore National Laboratory, to simulate the experiment. SHARP uses the spectral element code Nek5000 for fluid dynamics analysis and the finite element code DIABLO for structural analysis. The flow around the coil tubes is modeled in Nek5000 by using a large eddy simulation turbulence model. Transient pressure data on the tube surfaces is sampled and transferred to DIABLO for the structural simulation. The structural response is simulated in DIABLO via an implicit time-marching algorithm and a combination of continuum elements and structural shells. Tube vibration data (acceleration and frequency) are sampled and compared with the experimental data. Currently, only one-way coupling is used, which means that pressure loads from the fluid simulation are transferred to the structural simulation but the resulting structural displacements are not fed back to the fluid simulation

CFD↗
DOE OSTI2025

Lowering entry barriers to developing custom simulators of distributed applications and platforms with SimGrid

Researchers in parallel and distributed computing (PDC) often resort to simulation because experiments conducted using a simulator can be for arbitrary experimental scenarios, are less resource-, labor-, and time-consuming than their real-world counterparts, and are perfectly repeatable and observable. Many frameworks have been developed to ease the development of PDC simulators, and these frameworks provide different levels of accuracy, scalability, versatility, extensibility, and usability. Further, the SimGrid framework has been used by many PDC researchers to produce a wide range of simulators for over two decades. Its popularity is due to a large emphasis placed on accuracy, scalability, and versatility, and is in spite of shortcomings in terms of extensibility and usability. Although SimGrid provides sensible simulation models for the common case, it was difficult for users to extend these models to meet domain-specific needs. Furthermore, SimGrid only provided relatively low-level simulation abstractions, making the implementation of a simulator of a complex system a labor-intensive undertaking. In this work we describe developments in the last decade that have contributed to vastly improving extensibility and usability, thus lowering or removing entry barriers for users to develop custom SimGrid simulators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

MS–STEM–FEM: A parallelized multi-slice fluctuation TEM simulation tool

Atomic configurations of glassy or amorphous materials containing medium-range order (MRO) may be identified by comparing fluctuation transmission electron microscopy (FTEM) measurements to FTEM simulations obtained using model configurations. Candidate model sizes have traditionally been much thinner than the samples measured experimentally, and publicly available FTEM simulation software has until now omitted microscope parameters, dynamical scattering, and the phase of the diffracted electron wave. We introduce MS–STEM–FEM, an open-source software package for simulating FTEM experiments using established multi-slice TEM simulation techniques to emulate experiment more closely by incorporating microscope parameters and simulating electron scattering and propagation as a complex valued wave. Simulations using established models are compared with results of experimental STEM–FEM to validate the software. Several statistical measures of diffraction are implemented and their responses to model features are compared. Dynamical scattering is found to be less influential than the variety of crystallite orientations which occur in thicker models. Simulations of variable resolution microscopy confirm that cumulative intensity of the FTEM signal decreases with reduced model MRO and increased coherence volume. In conclusion, advantageous model scaling characteristics and efficient processor performance scaling are demonstrated, along with a study of convergence with respect to pertinent simulation parameters to identify accuracy requirements.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Structure and Evolution of Mesoscale Convective Systems: Sensitivity to Cloud Microphysics in Convection‐Permitting Simulations Over the United States

Abstract Regional climate simulations over the continental United States were conducted for the 2011 warm season using the Weather Research and Forecasting model at convection‐permitting resolution (4 km) with two commonly used microphysics parameterizations (Thompson and Morrison). Sensitivities of the simulated mesoscale convective system (MCS) properties and feedbacks to large‐scale environments are systematically examined against high‐resolution geostationary satellite and 3‐D mosaic radar observations. MCS precipitation including precipitation amount, diurnal cycle, and distribution of hourly precipitation intensity are reasonably captured by the two simulations despite significant differences in their simulated MCS properties. In general, the Thompson simulation produces better agreement with observations for MCS upper level cloud shield and precipitation area, convective feature horizontal and vertical extents, and partitioning between convective and stratiform precipitation. More importantly, Thompson simulates more stratiform rainfall, which agrees better with observations and results in top‐heavier heating profiles from robust MCSs compared to Morrison. A stronger dynamical feedback to the large‐scale environment is therefore seen in Thompson, wherein an enhanced mesoscale vortex behind the MCS strengthens the synoptic‐scale trough and promotes advection of cool and dry air into the rear of the MCS region. The latter prolongs the MCS lifetimes in the Thompson relative to the Morrison simulations. Hence, different treatment of cloud microphysics not only alters MCS convective‐scale dynamics but also has significant impacts on their macrophysical properties such as lifetime and precipitation. As long‐lived MCSs produced 2–3 times the amount of rainfall compared to short‐lived ones, cloud microphysics parameterizations have profound impact in simulating extreme precipitation and the hydrologic cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Refining fast calorimeter simulations with a Schrödinger Bridge

Machine learning-based simulations, especially calorimeter simulations, are promising tools for approximating the precision of classical high energy physics simulations with a fraction of the generation time. Nearly all methods proposed so far learn neural networks that map a random variable with a known probability density, like a Gaussian, to realistic-looking events. In many cases, physics events are not close to Gaussian and so these neural networks have to learn a highly complex function. We study an alternative approach: Schrödinger bridge Quality Improvement via Refinement of Existing Lightweight Simulations (SQuIRELS). SQuIRELS leverages the power of diffusion-based neural networks and Schrödinger bridges to map between samples where the probability density is not known explicitly. We apply SQuIRELS to the task of refining a classical fast simulation to approximate a full classical simulation. On simulated calorimeter events, we find that SQuIRELS is able to reproduce highly non-trivial features of the full simulation with a fraction of the generation time.

Calorimeter methods↗
DOE OSTI2025

Lie-algebraic classical simulations for quantum computing

The classical simulation of quantum dynamics plays an important role in our understanding of quantum complexity and in the development of quantum technologies. Efficient techniques such as those based on the Gottesman-Knill theorem for Clifford circuits, tensor networks for low entanglement-generating circuits, or Wick's theorem for fermionic Gaussian states have become central tools in quantum computing. In this work, we contribute to this body of knowledge by presenting a framework for classical simulations, dubbed “𝔤-sim”, which is based on the underlying Lie algebraic structure of the dynamical process. When the dimension of the algebra grows at most polynomially in the system size, there exist observables for which the simulation is efficient. Indeed, we show that 𝔤-sim enables new regimes for classical simulations, is able to deal with certain forms of noise in the evolution, as well as can be used to tackle several paradigmatic variational and nonvariational quantum computing tasks. For the former, we perform Lie-algebraic simulations to train and optimize parametrized quantum circuits (thus effectively showing that some variational models can be dequantized), design enhanced parameter initialization strategies, solve tasks of quantum circuit synthesis, and train a quantum-phase classifier. For the latter, we report large-scale noiseless and noisy simulations on benchmark problems. By comparing the limitations of 𝔤-sim and certain Wick's theorem-based simulations, we find that the two methods become inefficient for different types of states or observables, hinting at the existence of distinct, nonequivalent resources for classical simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Evaluation of the Inertial Response of Variable-Speed Wind Turbines Using Advanced Simulation

In this paper, we focus on the temporary frequency support effect provided by wind turbine generators (WTGs) through the inertial response. With the implemented inertial control methods, the WTG is capable of increasing its active power output by releasing parts of the stored kinetic energy when the frequency excursion occurs. The active power can be boosted temporarily above the maximum power points, but the rotor speed deceleration follows and an active power output deficiency occurs during the restoration of rotor kinetic energy. We evaluate and compare the inertial response induced by two distinct inertial control methods using advanced simulation. In the first stage, the proposed inertial control methods are analyzed in offline simulation. Using an advanced wind turbine simulation program, FAST with TurbSim, the response of the researched wind turbine is comprehensively evaluated under turbulent wind conditions, and the impact on the turbine mechanical components are assessed. In the second stage, the inertial control is deployed on a real 600kW wind turbine - Controls Advanced Research Turbine, 3-bladed (CART3), which further verifies the inertial control through a hardware-in-the-loop (HIL) simulation. Various inertial control methods can be effectively evaluated based on the proposed two-stage simulation platform, which combines the offline simulation and real-time HIL simulation. The simulation results also provide insights in designing inertial control for WTGs.

17 WIND ENERGY↗
DOE OSTI2021

Quasi-Static Time Series Fatigue Simulation for PV Inverter Semiconductors with Long-Term Solar Profile

Power system simulations with long-term data tend to have large time steps varying from one second to a few minutes. However, for PV inverter semiconductors, the minimum thermal stresses cycle is with line frequency. This requires the time step of the fatigue simulation to be much smaller than the line period. This small time step results in poor simulation speed, especially for long-term simulations. This paper proposes a fast fatigue simulation for inverter semiconductors using the quasi-static time series (QSTS) simulation concept. The fatigue analysis typically focuses on the peak and valley values of a strain and neglect the transients from peaks to valleys. The proposed simulation utilizes this property of fatigue analysis and calculates the steady state of the semiconductor junction temperature only. The resulting time step of the fatigue simulation is 15 minutes, which is consistent with the solar dataset without losing accuracy.

fast simulation↗
DOE OSTI2022

Electromagnetic Transient (EMT) Simulation Algorithms for Evaluation of Large-Scale Extreme Fast Charging Systems (T&D Models)

Simulation of high-fidelity models of extreme fast charging (XFC) systems and large-area power grids with many XFCs can be time consuming in traditional simulators. Traditional simulators use a single method of discretization for all the components that results in imposing a large computational burden of inverting a large matrix as well as increased computations related to single method of discretization (that is typically a trapezoidal method). To overcome the problem of simulating large-area power grids with many XFCs, in this paper, advanced numerical simulation algorithms are applied for the first time together to reduce the dimension of matrix inversion. Here, the algorithms include numerical stiffness-based segregation, time constant-based segregation, clustering and aggregation on differential algebraic equations (DAEs), and multi-order integration approaches. These algorithms apply multiple discretization algorithms rather than a single discretization algorithm that further reduces the computational burden. The approaches mentioned here have resulted in speed-up of up to 18x in the simulation of a single distribution system with 15 XFCs and of up to 271x in the simulation of a transmission-distribution system with 300 XFCs in multiple distribution feeders with respect to conventional simulators (like power systems computer aided design [PSCAD]).

42 ENGINEERING↗
DOE OSTI2022

Evaluation of the Fitch Wind-Farm Wake Parameterization with Large-Eddy Simulations of Wakes Using the Weather Research and Forecasting Model

Abstract Wind-farm parameterizations in weather models can be used to predict both the power output and farm effects on the flow; however, their correctness has not been thoroughly assessed. We evaluate the wind-farm parameterization of the Weather Research and Forecasting Model with large-eddy simulations (LES) of the wake performed with the same model. We study the impact on the velocity and turbulence kinetic energy (TKE) of inflow velocity, roughness, resolution, number of turbines (one or two), and inversion height and strength. We compare the mesoscale with the LES by spatially averaging the LES within areas correspondent to the mesoscale horizontal spacing: one covering the turbine area and two downwind. We find an excellent agreement of the velocity within the turbine area between the two types of simulations. However, within the same area, we find the largest TKE discrepancies because in mesoscale simulations, the turbine-added TKE has to be highest at the turbine position to be advected downwind. Within the downwind areas, differences between velocities increase as the wake recovers faster in the LES, whereas for the TKE both types of simulations show similar levels. From the various configurations, the impact of inversion height and strength is small for these heights and inversion levels. The highest impact for the one-turbine simulations appears under the low-speed case due to the higher thrust, whereas the impact of resolution is low for the large-eddy simulations but high for the mesoscale simulations. Our findings demonstrate that higher-fidelity simulations are needed to validate wind-farm parameterizations.

17 WIND ENERGY↗
DOE OSTI2020

Simulation as a sustainable trading zone: Aiming at intergenerational justice

The paper, drawing on the example of simulation codes used in nuclear physics and high-energy physics, seeks to call attention to the ethical implications of discontinuing support for simulation codes and the loss of knowledge embodied in them. Predicated on the concept of trading zones and actor network models, the paper addresses the problem of extinction of simulation codes and attempts to understand their evolution and development within those frameworks. We show that simulation codes of closed type develop to the level of creoles, becoming local languages and standards of scientific centers and disappearing as their few main developers leave, whereas codes of open types become universal languages, imposing problemsolving patterns on the entire community and crowding out other codes. The paper suggests that because of simulations’ reliance on tacit knowledge, practices entrenched in codes cannot be exhaustively explicated or transmitted through writing alone; on the contrary, the life cycle of a simulation code is determined by the life cycle of its trading zone. We examine the extent to which both of these phenomena pose a risk to the preservation of knowledge. Bearing upon intergenerational ethics, we draw analogies between the pure intergenerational problem (PIP) and the problem of preserving the knowledge implemented in simulation codes and transmitting it to future generations. We argue that for the complete transfer of knowledge, it is necessary to develop and maintain inhabitability and sustainability of simulation trading zones in a controllable way, at least until the demand for these codes is warranted to cease in the future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Ground-motions site and event specificity: Insights from assessing a suite of simulated ground motions in the San Francisco Bay Area

This article presents the results of a research that is part of a larger collaborative effort between the Lawrence Berkeley National Laboratory and the Pacific Earthquake Engineering Research Center, funded by the US Department of Energy Office of Cybersecurity, Energy Security and Emergency Response. The main objective of this study is to assess a suite of near and far-field simulated ground motions obtained from 20 realizations of an M7 Hayward Fault earthquake in the San Francisco Bay Area, California USA, and inform the selection of rupture simulation parameters leading to strong motions. To this aim, comparisons are conducted with NGA-W2 and directivity ground-motion models and a selected population of records. An archetypal steel moment-resisting frame is utilized to assess infrastructure response distributions. The analyses carried out for each simulated event and subdomain with consistent properties in terms of shallow shear-wave velocity proved to be instrumental for better interpreting the differences between simulated motions and empirical models. The main reasons identified for variances between simulations and empirical relationships included (1) directivity effects fully captured by the simulations across the full breadth of rupture models; (2) site vicinity to ruptures that incorporate large-slip patches, particularly if these are in the forward-directivity direction; and (3) presence of geologic structures that can “trap” seismic waves and produce ground motions with large amplitude and long signal duration. The analyses carried out in this work provide a path for interpreting ground-motion site and event specificity obtained from a suite of physics-based simulations, differing only in the rupture model characterization, to inform the selection of simulation scenarios for site-specific engineering analyses under strong excitations. Evidence from this work points to the possibility that current hazard models may underestimate ground-motion intensities in areas where the combined effect of directivity and site conditions results in large ground-motion amplitudes.

58 GEOSCIENCES↗
DOE OSTI2022

Simulation of SPND Beta-Current in AP1000 core

The Self-Powered Neutron Detectors (SPNDs) are widely used in Generation-III commercial PWR due to its small size and external-power free character. It can measure the in-core neutron flux density immediately to obtain the power profile in the core. To better understand the actual state of the simulated core, a SPND beta-current simulation module has been established in the NECP-Bamboo code. Firstly, neutron simulation is carried out in the lattice code Bamboo-Lattice together with nuclide depletion simulation to obtain the depleted composition of the entire lattice including the SPND, and the few-group parameters including density and cross section of emitter material. Secondly, the depleted lattice can then be piped into a Monte-Carlo code for neutron-photon-electron simulation to obtain the electron escape efficiency from the emitter of the SPND to the corresponding collector. Thirdly, the electron escape efficiency can be parameterized into the function of the lattice state parameters similarly to the common few-group constants. Finally, the beta-current was calculated by the core code Bamboo-Core simultaneously with the simulation of the core operation process. The beta-current of the Vanadium SPND in the AP1000 core was simulated and analyzed in this paper. Comparison between simulation and measurement of the beta-current distribution in the specific core state is shown in this paper to verify the effectiveness of the method and the code. The radial normalized current distribution is in good agreement with the measurement. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations

We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Latent Space Simulation for Carbon Capture Design Optimization

The CO2 capture efficiency in solvent-based carbon capture systems (CCSs) critically depends on the gas-solvent interfacial area (IA), making maximization of IA a foundational challenge in CCS design. While the IA associated with a particular CCS design can be estimated via a computational fluid dynamics (CFD) simulation, using CFD to derive the IAs associated with numerous CCS designs is prohibitively costly. However, previous works such as Deep Fluids (Kim et al., 2019) show that large simulation speed-ups and low error are achievable by replacing CFD simulators with neural-network (NN) surrogates that mimic the CFD simulation process. This raises the possibility for a fast, accurate replacement for a CFD simulator, and thus computationally feasible IA-based CCS-design optimization. As such, here, we explore whether an existing NN-surrogate approach (and variants we develop) can successfully be applied to our complex carbon-capture CFD simulations, with the ultimate goal of obtaining a fast and accurate simulator for our CCS-design application. Our experiments build on the Deep Fluids approach and find that resulting surrogates can produce large speed ups (4000x) while maintaining IA relative errors as low as 4% on unseen CCS configurations (interpolating between configurations seen during training). Thus, despite less faithfulness to the underlying physics of the problem, NN surrogates may be a promising tool for our CCS design optimization problem. Notably, though, the Deep Fluids approach has limitations for our application (such as model non-transferability to CCS-packing changes), which we discuss. We conclude with potential directions for future work that may, like innovations we introduced here (e.g., transformer-based dynamics prediction), improve performance on our complicated dataset of CCS CFD simulations.

AI Surrogate Simulation, Carbon Capture, design op↗
DOE OSTI2019

Blade-Resolved, Single-Turbine Simulations Under Atmospheric Flow

The ultimate goal of the ExaWind project is to enable scientific discovery through predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The objective of this milestone is to demonstrate a 'production run' of a megawatt-scale turbine to exercise and demonstrate the turbine simulation capabilities that have been developed during the previous two years of the ExaWind project. Demonstrating this simulation capability is the critical first step in enabling Nalu-Wind to be used as a predictive modeling tool that can advance the scientific understanding of wind turbine performance. In Section 2, the milestone description is provided. Next, in Section 3, we present the NREL 5-megawatt (MW) reference turbine that was used to demonstrate Nalu-Wind's capabilities and describe the four computational meshes that were generated. In Section 4, we present results from simulations of the turbine operating under uniform inflow and homogeneous isotropic turbulence (HIT) conditions. Finally, in Section 5, notable achievements, conclusions, and a path forward are described.

17 WIND ENERGY↗
DOE OSTI2024

Integrated Transmission and Distribution Co-Simulation Platform for Demonstration of Bulk Grid Services Using Distributed Energy Resources

In September 2020, the Federal Energy Regulatory Commission (FERC) released Order 2222, which opens wholesale markets to small-capacity distributed energy resources (DERs), recognizing their potential in improving operational efficiency by providing bulk grid services. Therefore, a co-simulation capability that can connect transmission and distribution (T&D) simulations and evaluate the impacts of DER provision of bulk grid services is needed. In this paper, we present a new integrated T&D co-simulation platform that incorporates T&D system simulators, DER aggregator/group strategies, and a co-simulation coordinator. Industry-standard communication protocols are employed to mimic real-world conditions. Secondary frequency regulation is selected as the representative bulk grid service, and we simulate the responses of DERs to the frequency regulation signals. The simulation results for a solar-rich distribution feeder in Colorado, USA, demonstrate how the T&D co-simulation setup is used to evaluate the contributions of DERs to minimize the bulk grid frequency deviation.

24 POWER TRANSMISSION AND DISTRIBUTION↗