SEARCH · Search NASA

Results for “quick 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Storage Sizing and Placement Simulation: Quick-Start Case Study User’s Guide

The Storage Sizing and Placement Simulation (SSIM) application allows a user to define the possible sizes and locations of energy storage elements on an existing grid model defined in OpenDSS. Given these possibilities, the software will automatically search through them and attempt to determine which configurations result in the best overall grid performance. This quick-start guide will go through, in detail, the creation of an SSIM model based on a modified version of the IEEE 34 bus test feeder system. There are two primary parts of this document. The first is a complete list of instructions with little-to-no explanation of the meanings of the actions requested. The second is a detailed description of each input and action stating the intent and effect of each. There are links between the two sections.

25 ENERGY STORAGE↗

DOE OSTI2021

Deeper Fluids

A framework for building and training neural network surrogates for CFD simulation, with a focus on quickly simulating volume fraction dynamics for carbon capture efficiency prediction. Augments the Deep Fluids (Kim et al., 2019) approach as described in our paper "Latent Space Simulation for Carbon Capture Design Optimization" (Bartoldson et al., 2021).

Nguyen, PhuocChanh N↗

DOE OSTI2024

Deep Generative Models for Fast Photon Shower Simulation in ATLAS

The need for large-scale production of highly accurate simulated event samples for the extensive physics programme of the ATLAS experiment at the Large Hadron Collider motivates the development of new simulation techniques. Building on the recent success of deep learning algorithms, variational autoencoders and generative adversarial networks are investigated for modelling the response of the central region of the ATLAS electromagnetic calorimeter to photons of various energies. The properties of synthesised showers are compared with showers from a full detector simulation using GEANT4 . Both variational autoencoders and generative adversarial networks are capable of quickly simulating electromagnetic showers with correct total energies and stochasticity, though the modelling of some shower shape distributions requires more refinement. This feasibility study demonstrates the potential of using such algorithms for ATLAS fast calorimeter simulation in the future and shows a possible way to complement current simulation techniques.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

The MFiX Particle-in-Cell Method (MFiX-PIC) Theory Guide

MFiX (Multiphase Flow with Interphase eXchanges) is an open-source multiphase flow solver developed at the National Energy Technology Laboratory. Within the code, users have access to a single phase or interpenetrating continua-based multiphase two-fluid model (TFM), a discrete element model (DEM), and a particle-in-cell model (PIC). TFM, DEM, and PIC can all be used to create multiphase simulations that include hydrodynamics, chemical reactions, and heat transfer.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2021

Solidification dynamics in metal additive manufacturing: analysis of model assumptions

Solidification characteristics in the meltpool drive the process-microstructure relationship which helps dictate the material properties of as-built parts in additive manufacturing; therefore, being able to accurately and quickly simulate the size, shape, and solidification characteristics in the melt pool is of great interest to the field. In this study, we investigate various important physical phenomena (dynamic material properties, fluid-flow, radiation and vaporization) which can either be included or neglected in a continuum finite volume model (FVM) and their effect on the solidification conditions. Additionally, since the simplest form of such a model (conduction only) has an analytic solution which is much faster, its viability is also considered. Since the inclusion of some of these physical phenomena will inherently change the net energy input as well as the amount of energy needed to achieve melting of a control volume, each set of included phenomena had an effective absorption efficiency which was calibrated to closely match the dimensions of the melt pool to that of the ground truth data. The ground truth data for this study was defined to be the output of the FVM which included all the physical phenomena (OF). This study then goes on to compare the effects on solidification conditions each of these calibrated models has. It was found that most of the change in solidification conditions comes from the inclusion of latent heat. A posterior correlation factor (PCF) is then introduced to enable an analytic model to predict similar solidification conditions to OF model.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Development of Building Design Optimization Methodology: Residential Building Applications

Building design optimization is a highly complex problem, requiring long computational running processes because of the many options that exist when a building is being designed. This paper introduces an integrated approach through which to perform this optimization within an acceptable time frame. The approach includes the methods of variable selection, model simplification, and a sequential optimization process. Using singular value decomposition, a large number of design variables is reduced to a smaller subset that can be solved more quickly through the optimization algorithm. To expedite the variable selection process, a modeling approach that quickly simulates annual energy consumption was developed to replace full annual energy simulations. The developed methodology was applied to two residential buildings in the US, and the results are discussed herein. To assess the accuracy of the integrated optimization methodology, the optimized life cycle costs are compaa variables demonstrating the strongest contributions in the optimization study were identified. The proposed methodology significantly shortened the time requirements for the optimization processes of the two case studies by 74% and 84%; the optimized life cycle costs were within 0.05% and 0.06%, respectively, of the optimum point.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Mahakala: A Python-based Modular Ray-tracing and Radiative Transfer Algorithm for Curved Spacetimes

We introduce Mahakala, a Python-based, modular, radiative ray-tracing code for curved spacetimes. We employ Google's JAX framework for accelerated automatic differentiation, which can efficiently compute Christoffel symbols directly from the metric, allowing the user to easily and quickly simulate photon trajectories through non-Kerr spacetimes. JAX also enables Mahakala to run in parallel on both CPUs and GPUs. Mahakala natively uses the Cartesian Kerr–Schild coordinate system, which avoids numerical issues caused by the pole in spherical coordinate systems. We demonstrate Mahakala's capabilities by simulating 1.3 mm wavelength images (the wavelength of Event Horizon Telescope observations) of general relativistic magnetohydrodynamic simulations of low-accretion rate supermassive black holes. The modular nature of Mahakala allows us to quantitatively explore how different regions of the flow influence different image features. We show that most of the emission seen in 1.3 mm images originates close to the black hole and peaks near the photon orbit. We also quantify the relative contribution of the disk, forward jet, and counterjet to 1.3 mm images.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Field Sensitivity Analysis of Turbulence Model Parameters for Flow Over a Wing

Reynolds-averaged-Navier-Stokes (RANS) turbulence models are a critical tool in computational-fluid-dynamics simulations of aerodynamic systems, but simulation results can be highly sensitive to RANS-model parameter choices. Sensitivity analysis can be used to quantify these impacts, and the objective of this study is to demonstrate field sensitivity analysis with respect to ten parameters in the 2003 Menter shear-stress-transport (SST) turbulence model. Here, the analysis is demonstrated for an application relevant to wind energy, namely, flow over a NACA 0015 wing at 12 degree angle of attack and a Reynolds number of 1.5 x 10 6 . We quantify sensitivity using Sobol indices and the mean-squared gradient, which are estimated using polynomial chaos and active subspace models, respectively. Our results indicate that there are substantial spatial variations in parameter sensitivities, with different sets of most-sensitive parameters near the wing, as well as in the downstream wake, consistent with the physical interpretations of the turbulence model inputs. We show that, for this particular turbulence model and flow, simultaneous dimension reduction is possible across all quantities of interest, enabling efficient exploration of model outcomes. Ultimately, this analysis provides new insights into turbulence model parameter sensitivities in incompressible flows, and also demonstrates the implementation of field sensitivity analysis for applications relevant to aerodynamics simulations.

17 WIND ENERGY↗

DOE OSTI2020

Dynamic Phasor Modeling of Various Multipulse Rectifiers and a VSI Fed by 18-Pulse Asymmetrical Autotransformer Rectifier Unit for Fast Transient Analysis

To study the fast modelling of dynamic characteristics caused by power electronic switching, the dynamic phasor (DP) theory based on the time-varying Fourier decomposition and frequency shifting is firstly applied to the modelling of various multipulse rectifiers. The DP model for symmetrical 12-pulse phase-shifting reactor rectifier unit (PSR-RU) is derived with model order reduction, relations between ac and dc terminals, and Taylor series expansion. The DP model of asymmetrical 18-pulse autotransformer rectifier unit (AT-RU) is proposed based on the switching functions expressed in the DP domain. Meanwhile, the DP model of voltage source inverter (VSI) fed by asymmetrical 18-pulse AT-RU is built with the harmonic state-space (HSS) equations. Under both balanced and unbalanced conditions, the good calculation accuracy and rapid simulation speed of the developed DP models are validated by the detailed time-domain (TD) simulation.

42 ENGINEERING↗

DOE OSTI2023

QES-Plume v1.0: a Lagrangian dispersion model

Low-cost simulations providing accurate predictions of transport of airborne material in urban areas, vegetative canopies, and complex terrain are demanding because of the small-scale heterogeneity of the features influencing the mean flow and turbulence fields. Common models used to predict turbulent transport of passive scalars are based on the Lagrangian stochastic dispersion model. The Quick Environmental Simulation (QES) tool is a low-computational-cost framework developed to provide high-resolution wind and concentration fields in a variety of complex atmospheric-boundary-layer environments. Part of the framework, QES-Plume, is a Lagrangian dispersion code that uses a time-implicit integration scheme to solve the generalized Langevin equations which require mean flow and turbulence fields. Here, QES-Plume is driven by QES-Winds, a 3D fast-response model that computes mass-consistent wind fields around buildings, vegetation, and hills using empirical parameterizations, and QES-Turb, a local-mixing-length turbulence model. In this paper, the particle dispersion model is presented and validated against analytical solutions to examine QES-Plume’s performance under idealized conditions. In particular, QES-Plume is evaluated against a classical Gaussian plume model for an elevated continuous point-source release in uniform flow, the Lagrangian scaling of dispersion in isotropic turbulence, and a non-Gaussian plume model for an elevated continuous point-source release in a power-law boundary-layer flow. In these cases, QES-Plume yields a maximum relative error below 6 % when compared with analytical solutions. In addition, the model is tested against wind-tunnel data for a uniform array of cubical buildings. QES-Plume exhibits good agreement with the experiment with 99 % of matched zeros and 59 % of the predicted concentrations falling within a factor of 2 of the experimental concentrations. Furthermore, results also emphasize the importance of using high-quality turbulence models for particle dispersion in complex environments. Finally, QES-Plume demonstrates excellent computational performance.

58 GEOSCIENCES↗

DOE OSTI2020

A Holistic Algorithmic Approach to Improving Accuracy, Robustness, and Computational Efficiency for Atmospheric Dynamics

Atmospheric weather and climate models must perform simulations very quickly to be useful. Therefore, modelers have traditionally focused on reducing computations as much as possible. However, in our new era of increasingly compute-capable hardware, data movement is now the prohibiting expense. This study examines the computational benefits of a new algorithmic approach to modeling atmospheric dynamics on scales relevant to weather and climate simulation. Rather than minimizing computations, this new approach considers the larger problem more holistically, including spatial accuracy, temporal accuracy, robustness (i.e., oscillations), on-node efficiency, and internode data transfers together at once. Numerical experiments demonstrate how computations can be strategically increased to simultaneously address each of these constraints while reducing data movement to adapt to modern accelerated hardware. The new algorithm can achieve at times up to 80% peak floating point throughput in single precision on the Nvidia Tesla V100 GPU, where the traditional approach is shown to only achieve single-digit floating point efficiency. Further, the new algorithm is twice as fast as a standard Runge--Kutta time integrator, and high-order accuracy with Weighted Essentially Non-Oscillatory (WENO) limiting came at less than 30% additional runtime cost on a GPU, thus increasing the accuracy per degree of freedom.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Progress on Optimizing Wind Farms and Rotor Designs Using Adjoints

Modern wind plants are increasingly tasked with multiple performance objectives. In addition to designing plants that maximize power output and minimize the levelized cost of energy (LCOE), the design and operation of wind plants is increasingly influenced by challenges regarding grid integration of variable generation renewables. This places a growing emphasis on making wind plants more controllable and predictable. WindSE is a Reynolds-averaged Navier-Stokes (RANS) model designed around analytical gradient and adjoint methods, with the ability to capture terrain-induced effects, as shown in Figure 1. The recent addition of an unsteady solver with an actuator line method (ALM) and ongoing work to enable massively parallel optimizations gives it a unique niche to explore coupled plant-level controls and design problems. This code is an open source python package built on the FEniCS framework that utilizes fast, parallel PETSc solvers to model fluid flow throughout wind-farm scale domains. Two recent studies performed using WindSE demonstrate the capability to optimize under a wide variety of flow conditions and objective functions. In the first, we present an optimization focused on modifying the layout of a wind farm with a fixed number of turbines for maximum total power output [1]. This study highlights the ability to quickly perform simulations using the steady Navier-Stokes solver combined with rotors represented as actuator disks while also stressing the importance of capturing terrain-induced effects. Gradient-based optimization using the RANS equations is viable due to the inclusion of efficiently computed adjoint derivatives. We interpret the physical results of the optimal layout and also discuss the computational cost of scaling to larger problems. In the second study, we present the capabilities of the unsteady Navier-Stokes solver, where rotor-blade profiles represented by actuator lines are optimized to enhance wake steering effects and overall power production [2]. We quantify the wind plant performance gains obtained from this type of simultaneous control co-design optimization as compared to optimizing the blade design and yaw independently. Figure 2 shows the differences between a baseline two-turbine system and an optimized system where we fine-tune the blade chord profile. Results and challenges from each study are quickly summarized and used to motivate the current development efforts within WindSE. Current and future work is focused on enabling higher-resolution studies with more degrees of freedom through parallelization of both the simulation and optimization algorithms. We present benchmarking results to show that WindSE performs well in both weak- and strong-scaling tests and further demonstrate that the optimizer obtains the same convergence rates in both shared- and distributed-memory environments. Using larger wind farms, we can study deep-array effects within an optimization context, allowing the use of objective functions that have been previously unstudied. As an example, we present ongoing work on a blockage metric which characterizes the loss of available kinetic energy due to wake effects from multiple upstream turbines.

adjoint optimization↗

DOE OSTI2023

Impacts of Permeability Heterogeneity and Background Flow on Supercritical CO 2 Dissolution in the Deep Subsurface

Motivated by CO 2 capture and sequestration (CCS) design considerations, we consider the coupled effects of permeability heterogeneity and background flow on the dissolution of a supercritical CO 2 lens into an underlying deep, confined aquifer. We present the results of a large-scale Monte Carlo simulation study examining the interaction of background flow rate and three parameters describing multi-Gaussian log-permeability fields: mean, variance, and correlation length. Hundreds of high-resolution simulations were performed using the PFLOTRAN finite volume software to model CO 2 dissolution in a kilometer-scale aquifer over 1,000 years. Predictive dimensionless scaling relationships relating CO 2 dissolution rate to heterogeneity statistics, Rayleigh (Ra) and Péclet (Pe) numbers were developed for both gravitationally dominated free convection to background flow-dominated forced convection regimes. An empirical criterion, Pe = Ra 3/4 , was discovered for regime transition. All simulations converged quickly to a quasi-steady, approximately linear dissolution rate. However, this rate displayed profound variability between permeability field realizations sharing the same heterogeneity statistics, even under mild permeability heterogeneity. In general, increased heterogeneity was associated with a lower mean and higher variance of dissolution rate, undesirable from a CCS design perspective. The relationship between dissolution rate and background flow was found to be complex and nonlinear. Dimensionless scaling relationships were uncovered for a number of special cases. Results call into question the validity of the Boussinesq approximation in the context of modest-to-high background flow rates and the general applicability of numerical simulations without background flow.

42 ENGINEERING↗

DOE OSTI2023

A riboswitch separated from its ribosome-binding site still regulates translation

Abstract Riboswitches regulate downstream gene expression by binding cellular metabolites. Regulation of translation initiation by riboswitches is posited to occur by metabolite-mediated sequestration of the Shine-Dalgarno sequence (SDS), causing bypass by the ribosome. Recently, we solved a co-crystal structure of a prequeuosine1-sensing riboswitch from Carnobacterium antarcticum that binds two metabolites in a single pocket. The structure revealed that the second nucleotide within the gene-regulatory SDS, G34, engages in a crystal contact, obscuring the molecular basis of gene regulation. Here, we report a co-crystal structure wherein C10 pairs with G34. However, molecular dynamics simulations reveal quick dissolution of the pair, which fails to reform. Functional and chemical probing assays inside live bacterial cells corroborate the dispensability of the C10–G34 pair in gene regulation, leading to the hypothesis that the compact pseudoknot fold is sufficient for translation attenuation. Remarkably, the C. antarcticum aptamer retained significant gene-regulatory activity when uncoupled from the SDS using unstructured spacers up to 10 nucleotides away from the riboswitch—akin to steric-blocking employed by sRNAs. Accordingly, our work reveals that the RNA fold regulates translation without SDS sequestration, expanding known riboswitch-mediated gene-regulatory mechanisms. The results infer that riboswitches exist wherein the SDS is not embedded inside a stable fold.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Transfer Operator Framework for Earth System Predictability and Water Cycle Extremes

For chaotic dynamical systems, nonlinear instabilities lead to exponentially divergent trajectories in the evolution of system states. Unless a simulation is initialized with an infinite-precision snapshot of the state of the true system and all known physical effects that go into its evolution are directly computed, the future state predicted by the simulation will quickly diverge from that of the true system. Moreover, the Earth system is highly structured and contains localized coherent structures that are particularly important to predict. Predicting extreme events associated with coherent structures, like hurricanes and blocking events, is crucial for understanding the effects of global warming on the water cycle.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Transfer Operator Framework for Earth System Predictability and Water Cycle Extremes

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Adaptive Machine Learning for Robust Diagnostics and Control of Time-Varying Particle Accelerator Components and Beams

Machine learning (ML) is growing in popularity for various particle accelerator applications including anomaly detection such as faulty beam position monitor or RF fault identification, for non-invasive diagnostics, and for creating surrogate models. ML methods such as neural networks (NN) are useful because they can learn input-output relationships in large complex systems based on large data sets. Once they are trained, methods such as NNs give instant predictions of complex phenomenon, which makes their use as surrogate models especially appealing for speeding up large parameter space searches which otherwise require computationally expensive simulations. However, quickly time varying systems are challenging for ML-based approaches because the actual system dynamics quickly drifts away from the description provided by any fixed data set, degrading the predictive power of any ML method, and limits their applicability for real time feedback control of quickly time-varying accelerator components and beams. In contrast to ML methods, adaptive model-independent feedback algorithms are by design robust to un-modeled changes and disturbances in dynamic systems, but are usually local in nature and susceptible to local extrema. In this work, we propose that the combination of adaptive feedback and machine learning, adaptive machine learning (AML), is a way to combine the global feature learning power of ML methods such as deep neural networks with the robustness of model-independent control. We present an overview of several ML and adaptive control methods, their strengths and limitations, and an overview of AML approaches.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

GADRAS Batch Inject Tool User Guide

Gamma Detector Response and Analysis Software (GADRAS) is used by the radiation detection and emergency response community to perform modeling and spectral analysis for gamma detector systems. Built into GADRAS is the ability to define a detector, geometry, background characteristics and source composition to generate synthetic spectra for drills and exercises (injects). Consequence Management is currently in development of a sample result data simulator tool in which a deposition model is probed for source conditions at moments in time and locations in space. These values are used to generate realistic sample results for use in drills and exercises. In addition to sample results, there is a need to simulate the actual spectra that would be observed in the field by downlooking HPGe instruments given a deposition activity. This way, the FRMAC Gamma Spectroscopist can practice their process of generating quantified results from spectra on realistic data as well. Recognizing the decades of work done in GADRAS to accurately generate synthetic spectra, this team decided to build a link between the new simulator and GADRAS to generate these spectra quickly and easily. The simulator tool will generate a file that specifies the name of the spectra, its location, date/time of measurement, duration of measurement, height off the ground, and the deposition activity and age for every radionuclide in the simulation. Then, a new tool within the Inject Tab of GADRAS was developed to read in this file given a detector selection and generate In-Situ spectra for each row in the file in any file format the user chooses. This way, simulation cell staff can take these files and then upload them to the appropriate data system (RAMS or RadResponder) for use during drills and exercises. An advanced feature of this tool allows for generating any spectra given an appropriate model and mapping of source to model layer in the batch inject tool. This way, spectra from field sample counts, mobile laboratories, or even fixed laboratories can be generated in bulk given an estimate of the radioactivity concentration or total radioactivity in an import file. This expands the capabilities of this tool a great deal and will make it a more useful tool for CM and others to help estimate detector response for nearly any situation. This user guide will explain the steps needed to perform a batch inject file generation.

61 RADIATION PROTECTION AND DOSIMETRY↗