Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation techniques”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

2D-to-3D image translation of complex nanoporous volumes using generative networks

Image-based characterization offers a powerful approach to studying geological porous media at the nanoscale and images are critical to understanding reactive transport mechanisms in reservoirs relevant to energy and sustainability technologies such as carbon sequestration, subsurface hydrogen storage, and natural gas recovery. Nanoimaging presents a trade off, however, between higher-contrast sample-destructive and lower-contrast sample-preserving imaging modalities. Furthermore, high-contrast imaging modalities often acquire only 2D images, while 3D volumes are needed to characterize fully a source rock sample. In this work, we present deep learning image translation models to predict high-contrast focused ion beam-scanning electron microscopy (FIB-SEM) image volumes from transmission X-ray microscopy (TXM) images when only 2D paired training data is available. We introduce a regularization method for improving 3D volume generation from 2D-to-2D deep learning image models and apply this approach to translate 3D TXM volumes to FIB-SEM fidelity. We then segment a predicted FIB-SEM volume into a flow simulation domain and calculate the sample apparent permeability using a lattice Boltzmann method (LBM) technique. Results show that our image translation approach produces simulation domains suitable for flow visualization and allows for accurate characterization of petrophysical properties from non-destructive imaging data.

58 GEOSCIENCES↗
DOE OSTI2023

Statistical Generation of Ocean Forcing With Spatiotemporal Variability for Ice Sheet Models

Melting of ice at the base of floating ice shelves that fringe the Antarctic ice sheet has been identified as a significant source of uncertainty in sea level rise projections. Part of this uncertainty derives from chaotic internal variability of the coupled ocean-atmosphere system. For numerical ice sheet model projections, this uncertainty has not previously been quantified because of the prohibitive computational expense of running large climate model ensembles. Here, we develop and demonstrate a technique that generates independent realizations of internal climate variability from a single climate model simulation. Building on previous developments in model emulation, this technique uses empirical orthogonal function decomposition and Fourier-phase randomization to generate statistically consistent realizations of spatiotemporal variability fields for the target climate variable. The method facilitates efficient sampling of a wide range of climate trajectories, which can also be incorporated within ice sheet or other physical models to represent feedback processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Validation of the stochastic inversion algorithm for acoustic travel-time tomography: a large eddy simulation study

Acoustic tomography (AT) is explored as a remote sensing technique to obtain instantaneous snapshots of temperature and velocity fluctuations for wind energy applications. This study integrates Large Eddy Simulation (LES) with the Stochastic Inversion (SI) method to validate the algorithm’s capacity for accurate reconstruction of atmospheric fluctuations. The initial findings demonstrate the efficacy of the method in accurately capturing the predominant flow structures. Normalized L2 error evaluations further inform the algorithm’s precision, with errors accentuated in less sampled peripheral regions. The results underscore the method’s promise as a non-intrusive observational tool, with ongoing development poised to improve its precision and reliability.

17 WIND ENERGY↗
DOE OSTI2024

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial-flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that are neglected for wind turbines. These include buoyancy, added mass and inertial loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbines at a speed comparable to real time. While efficient for large sets of load cases and design studies, mid-fidelity codes make simplifying assumptions that may impact their accuracy. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system was modeled both with OpenFAST and with the commercial CFD code STAR-CCM+. The CFD model used a three-dimensional unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modeled with a custom coupling of the CFD solver with the open-source code MoorDyn. This improves upon the existing quasi-static catenary solver in STAR-CCM+, which lacks seabed contact or line-to-line connections. Simulation results for a floating marine turbine are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗
DOE OSTI2018

In Situ Analysis and Visualization of Fusion Simulations: Lessons Learned

The trends in high performance computing, where far more data can be computed that can ever be stored, have made in situ techniques an important area of research and development. Simulation campaigns, where domain scientists work with computer scientists to run a simulation and perform in situ analysis and visualization are important, and complex undertakings. In this paper we report our experiences performing in situ analysis and visualization on two campaigns. The two campaigns were related, but had important differences in terms of the codes that were used, the types of analysis and visualization required, and the visualization tools used. Further, we report the lessons learned from each campaign.

Kim, Mark B.↗
DOE OSTI2024

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System: Preprint

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory (NREL) has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that can be neglected for wind turbines. These include buoyant loads, added mass and inertia loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbine systems at a speed comparable to real time. While efficient for long simulations, large sets of load cases, and design studies, mid-fidelity codes cannot capture all of the potentially important physical phenomenon impacting marine turbine systems. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system designed by NREL was modeled both with OpenFAST and with the commercial CFD code, STARCCM+. The CFD model used a 3-D unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modelled with a custom coupling of the CFD solver with the open-source code, MoorDyn. This improves upon the existing quasi-static catenary solver in STARCCM+, which lacks seabed contact or line-to-line connections. Spatial and temporal convergence studies were conducted. The simulation results for a combined current and wave condition are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗
DOE OSTI2023

ReaxFF reactive force field model enables accurate prediction of physiochemical and mechanical properties of crystalline and amorphous shape‐memory polyurethane

Abstract We investigate reactive force field (ReaxFF) prediction of physiochemical and mechanical properties, and stimuli response behavior of a shape‐memory polyurethane (SMPU) rigid segment. We used SMPUs as a platform to link the rigid/hard segment domain crystallinity with properties that can also be experimentally measured. Specifically, we developed models for crystalline and amorphous forms of a common rigid segment of SMPUs, for example, 4,4′‐diphenylmethane diisocyanate (MDI) with n‐butanediol (BDO), denoted as MDI‐BDO. ReaxFF molecular dynamics simulations with equilibrium dynamics are performed by controlling mass, temperature, and pressure or volume to predict the structural and physiochemical properties of crystalline and amorphous rigid segments of SMPU systems, followed by simulating their XRD and FTIR patterns, respectively. Non‐equilibrium ReaxFF simulations are implemented by using uniaxial box deformation technique to study response‐behavior of crystalline SMPU under tensile loading and during stress‐relaxation. The model and simulation scheme presented in this study demonstrates a robust theoretical pathway toward application‐dictated, on‐demand design of high‐performance SMPUs in applications such as soft robotics.

Afroz, Mohammad M.↗
DOE OSTI2020

Simulation of the Charge Measurements for the SuperCDMS Soudan Experiment

Astronomical observations lead us to believe that dark matter constitutes the majority of the mass of the Universe, yet it has never been directly observed experimentally. The Super Cryogenic Dark Matter Search Experiment (SuperCDMS) aims to discover the particle nature of dark matter, by developing and using some of the most sensitive detectors ever built. The SuperCDMS detectors are capable of rejecting known particles to search for dark matter interactions, but the discrimination is a complicated analysis task. Previous analyses with SuperCDMS set world leading sensitivity limits to dark matter interactions, but the analysis techniques have not been based on simulations. Now we have developed an efficient simulations infrastructure that will help better understand the detectors, develop a better background model, and optimize the next generation dark matter searches. In this thesis we present the fully-functional simulations infrastructure to produce high quality samples a s well as results for the charge system readout simulations, and comparisons to expectations from data. We show that our simulations well-reproduce dominant features of data, bringing SuperCDMS closer to a robust simulation-based analysis method, much needed for the upcoming SuperCDMS dark matter searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Scientific Data Compression for Large Scale Computational Fluid Dynamics (CFD) Simulations

This Cooperative Research and Development Agreement (CRADA) between Oak Ridge National Laboratory (ORNL) and General Electric (GE) investigated methods for reducing the size of large computational fluid dynamics (CFD) simulation datasets using scientific data compression techniques. The work focused on adapting the MultiGrid Adaptive Reduction of Data (MGARD) compression framework and integrating it with high-performance I/O and visualization tools used in CFD workflows. MGARD uses hierarchical multilevel decomposition to enable error-controlled compression of floating-point scientific data while preserving quantities of interest. During the project, MGARD compression was integrated with the ADIOS I/O framework and visualization tools such as ParaView to enable efficient storage, transfer, and analysis of simulation data. The collaboration also explored approaches for improving compression performance for CFD data defined on unstructured meshes. Results demonstrate that scientific data compression can significantly reduce storage requirements and improve data management for large-scale CFD simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Review of Lattice-Boltzmann Models Coupled with Geochemical Modeling Applied for Simulation of Advanced Waterflooding and Enhanced Oil Recovery Processes

To maintain economic profit and improve the oil production efficiency after the primary and secondary production phase, advanced waterflooding techniques such as low salinity waterflooding in carbonate reservoirs have been investigated in numerical simulations, laboratory experiments, and field pilot tests. Multiple underlying mechanisms have been proposed based on these studies, and they are still under debate. Various numerical modeling approaches are introduced, but there exists a lack of a pore-scale comprehensive modeling scheme to fully understand the processes. Lattice-Boltzmann method (LBM) is a type of numerical fluid flow modeling technique that shows capabilities and flexibilities in modeling pore-scale fluid flow to integrate physical–chemical processes within complex structures. The intrinsic feature of LBM makes it a promising framework for simulating advanced waterflooding due to its flexibility, accuracy, and parallel efficiency. LBM works either by itself for solving reactive transport problems or by coupling with a third-party reaction solver. This review mainly introduces the LBM fluid flow and reactive transport capabilities and the concept and modeling approaches to simulate advanced waterflooding techniques. Meanwhile, an evaluation of the coupled LBM models for enhanced oil recovery (EOR) simulations is discussed with future research challenges and directions concluded.

02 PETROLEUM↗
DOE OSTI2025

Partnership Center for High-Fidelity Boundary Plasma Simulation (Final Report)

Within the Partnership Center for High-Fidelity Boundary Plasma Simulation (HBPS), work at UT-Austin was aimed at improved verification, validation, and uncertainty quantification (VVUQ) for edge plasma simulations and on performing gyrokinetics simulations of pedestal instabilities and turbulence in order to expand foundational understanding of pedestal transport. Regarding VVUQ, the accomplishments can be summarized as follows. First, it was shown that the Moment Preserving Constrained Resampling technique, when applied periodically in particle-in-cell simulations in the XGC code, can dramatically improve the accuracy of the simulation at essentially equivalent computational cost. Second, a technique for estimating model correlations, which are required to solve the model selection and sample allocation problem in multifidelity UQ techniques, without sampling the highest fidelity, most computationally expensive model, was developed and demonstrated. Third, previously developed methods for estimating statistical and discretization errors were applied to numerical methods relevant to edge plasma simulations, namely in particle-in-cell-based approaches, and shown to work. Finally, benchmark studies for comparing gyrokinetic codes were developed and performed, leading to reasonable agreement between four commonly used codes. Regarding physics studies, gyrokinetic simulations to investigate microtearing modes in the DIII-D pedestal were performed using the GENE code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling Dust Production, Growth, and Destruction in Reionization-era Galaxies with the CROC Simulations: Methods and Parameter Exploration

We introduce a model for the explicit evolution of interstellar dust in a cosmological galaxy formation simulation. We post-process a simulation from the Cosmic Reionization on Computers project, integrating an ordinary differential equation for the evolution of the dust-to-gas ratio along pathlines in the simulation sampled with a tracer particle technique. This model incorporates the effects of dust grain production in asymptotic giant branch star winds and supernovae, grain growth due to the accretion of heavy elements from the gas phase of the interstellar medium, and grain destruction due to thermal sputtering in the high-temperature gas of supernova remnants due to thermal sputtering. A main conclusion of our analysis is the importance of a carefully chosen dust destruction model, for which different reasonable parameterizations can predict very different values at the ∼100 pc resolution of the interstellar medium (ISM) in our simulations. We run this dust model on the single most massive galaxy in a 10 h$^{−1}$ comoving megaparsec box, which attains a stellar mass of ∼2 × 10$^{9}$ M $_{⊙}$ by z = 5. We find that the model is capable of reproducing dust masses and dust-sensitive observable quantities broadly consistent with existing data from high-redshift galaxies. The total dust mass in the simulated galaxy is somewhat sensitive to parameter choices for the dust model, especially the timescale for grain growth due to accretion in the ISM. Consequently, observable quantity observations that can constrain galaxy dust masses at these epochs are potentially useful for placing constraints on dust physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Ising meson spectroscopy on a noisy digital quantum simulator

Abstract Quantum simulation has the potential to be an indispensable technique for the investigation of non-perturbative phenomena in strongly-interacting quantum field theories (QFTs). In the modern quantum era, with Noisy Intermediate Scale Quantum (NISQ) simulators widely available and larger-scale quantum machines on the horizon, it is natural to ask: what non-perturbative QFT problems can be solved with the existing quantum hardware? We show that existing noisy quantum machines can be used to analyze the energy spectrum of several strongly-interacting 1+1D QFTs, which exhibit non-perturbative effects like ‘quark confinement’ and ‘false vacuum decay’. We perform quench experiments on IBM’s quantum simulators to compute the energy spectrum of 1+1D quantum Ising model with a longitudinal field. Our results demonstrate that digital quantum simulation in the NISQ era has the potential to be a viable alternative to numerical techniques such as density matrix renormalization group or the truncated conformal space methods for analyzing QFTs.

Science & Technology - Other Topics↗
DOE OSTI2018

A processing diagram for high-density Ti-6Al-4V by selective laser melting

Density optimization is the first critical step in building additively manufactured parts with highquality and good mechanical properties. In earlier work, we developed an approach that combines simulations and experiments to identify processing parameters for high-density parts. Here we present the results for Ti-6Al-4V using the laser powder-bed-fusion technique. We show how simple, but approximate, simulations are used to design parameters for a limited set of single track experiments. The resulting melt-pool characteristics are then used to identify processing parameters for high-density pillars. From these experimental data, a processing diagram based on the normalized energy density concept is constructed, illustrating an optimized processing window for high- or low-density samples. Excellent mechanical properties are obtained for Ti-6Al-4V samples built from the optimized window. We discuss the limitations and implications of the processing diagram on the microstructure and mechanical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Supporting Co-Design of Extreme-Scale Systems with In Situ Visual Analysis of Event-Driven Simulations (Final Report)

Extreme-scale systems for compute- and data-centric science are pushing the boundaries of system architectures in order to achieve ambitious goals in terms of computing capability, resilience, and power efficiency. Future systems necessitate emerging designs that incorporate new technologies, system organizations, and computer science techniques, and these changes are not yet well understood. Discrete-event-driven simulation of computing system architectures and subsystems has emerged as a productive and cost-effective means to evaluating potential designs, along with capabilities for executing these simulations of extreme-scale systems. However, techniques for analyzing the behavior of these simulations have not kept pace with our ability to perform the simulations, preventing us from extracting the most value from this promising approach. The objective of this project is to support the codesign of extreme-scale system architectures for compute and data-centric science through research and development of novel methods for analysis and visualization of large-scale event-driven simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Convective and turbulent motions in non-precipitating Cu. Part 1: Method of separation of convective and turbulent motions

Atmospheric motions in clouds and cloud surroundings have a wide range of scales, from several kilometers to centimeters. These motions have different impacts on cloud dynamics and microphysics. Larger-scale motions (hereafter referred to as convective motions) are responsible for mass transport over distances comparable with cloud scale, while motions of smaller scales (hereafter referred to as turbulent motions) are stochastic and responsible for mixing and cloud dilution. This distinction substantially simplifies the analysis of dynamic and microphysical processes in clouds. The present research is Part I of the study aimed at describing the method for separating the motion scale into a convective component and a turbulent component. An idealized flow is constructed, which is a sum of an initially prescribed field of the convective velocity with updrafts in the cloud core and downdrafts outside the core, and a stochastic turbulent velocity field obeying the turbulent properties, including the -5/3 law and the 2/3 structure function law. A wavelet method is developed allowing separation of the velocity field into the convective and turbulent components, with parameter values being in a good agreement with those prescribed initially. The efficiency of the method is demonstrated by an example of a vertical velocity field of a cumulus cloud simulated using the System for Atmospheric Modeling (SAM) with bin microphysics and resolution of 10 m. It is shown that vertical velocity in clouds indeed can be represented as a sum of convective velocity (forming zone of cloud updrafts and subsiding shell) and a stochastic velocity obeying laws of homogeneous and isotropic turbulence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An Evaluation of Dynamical Downscaling Methods Used to Project Regional Climate Change

In the past decade, dynamical downscaling using “pseudo-global-warming” (PGW) techniques has been applied frequently to project regional climate change. Such techniques generate signals by adding mean global climate model (GCM)-simulated climate change signals in temperature, moisture, and circulation to lateral and surface boundary conditions derived from reanalysis. An alternative to PGW is to downscale GCM data directly. This technique should be advantageous, especially for simulation of extremes, since it incorporates the GCM's full spectrum of changing synoptic-scale dynamics in the regional solution. Here, we test this assumption, by comparing simulations in Europe and Western North America. We find that for warming and changes in temperature extremes, PGW often produces similar results to direct downscaling in both regions. For mean and extreme precipitation changes, PGW generally also performs surprisingly well in many cases. Moisture budget analysis in the Western North America domain reveals why. Large fractions of the downscaled hydroclimate changes arise from mean changes in large-scale thermodynamics and circulation, that is, increases in temperature, moisture, and winds, included in PGW by design. The one component PGW may have difficulty with is the contribution from changes in synoptic-scale variability. When this component is large, PGW performance could be degraded. Global analysis of GCM data shows there are regions where it is large or dominant. Hence, our results provide a road map to identify, through GCM analyses, the circumstances when PGW would not be expected to accurately regionalize GCM climate signals.

54 ENVIRONMENTAL SCIENCES↗