Search NASA⌕ Search

SEARCH · Search NASA

Results for “differentiable simulations”

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 361 records · Page 20

GrainGNN: A dynamic graph neural network for predicting 3D grain microstructure

We propose GrainGNN, a surrogate model for the evolution of polycrystalline grain structure under rapid solidification conditions in metal additive manufacturing. High fidelity simulations of solidification microstructures are typically performed using multicomponent partial differential equations (PDEs) with moving interfaces. The inherent randomness of the PDE initial conditions (grain seeds) necessitates ensemble simulations to predict microstructure statistics, e.g., grain size, aspect ratio, and crystallographic orientation. Here, currently such ensemble simulations are prohibitively expensive and surrogates are necessary.In GrainGNN, we use a dynamic graph to represent interface motion and topological changes due to grain coarsening. We use a reduced representation of the microstructure using hand-crafted features; we combine pattern finding and altering graph algorithms with two neural networks, a classifier (for topological changes) and a regressor (for interface motion). Both networks have an encoder-decoder architecture; the encoder has a multi-layer transformer long-short-term-memory architecture; the decoder is a single layer perceptron.We evaluate GrainGNN by comparing it to high-fidelity phase field simulations for in-distribution and out-of-distribution grain configurations for solidification under laser power bed fusion conditions. GrainGNN results in 80%–90% pointwise accuracy; and nearly identical distributions of scalar quantities of interest (QoI) between phase field and GrainGNN simulations compared using Kolmogorov-Smirnov test. GrainGNN's inference speedup (PyTorch on single x86 CPU) over a high-fidelity phase field simulation (CUDA on a single NVIDIA A100 GPU) is 150×–2000× for 100-initial grain problem. Further, using GrainGNN, we model the formation of 11,600 grains in 220 seconds on a single CPU core.

36 MATERIALS SCIENCE↗

Computational Methods to Characterize Panel Loading Conditions for Accelerated Testing

Panel cracking and degradation due to wind loading are known to have a detrimental effect on power output. Recreating these damaging conditions in a controlled experimental setting requires an understanding of the loads being generated at the panel surface as a function of both wind speed, wind direction, and panel orientation. To better understand these relationships, a computational fluid dynamics (CFD) simulation package was constructed using an open-source Python library for solving partial differential equations with the finite element method. This CFD package allows the simulated wind speed and panel orientation of a multi-panel array to be easily changed and provides the capability to measure the traction forces at discrete points along the sun-facing and ground-facing surface of an interior panel residing in the wake of one or more upstream panels. A parameter exploration was performed in which the panel angle was varied in even increments from -70 deg to 70 deg and the wind speed was varied from 2 to 30 m/s. During post processing, the measured traction along the panel surface was averaged spatially and interpreted as a time-varying signal, where further processing of this signal yielded the root mean square amplitude and a characteristic frequency associated with the loading. This study found that higher wind speeds are generally associated with increased amplitude of loading and that the panel orientation angle can significantly exacerbate or mitigate this loading. These outcomes are presented along with current work on higher-fidelity verification simulations and recommendations for performing accelerated experimental testing.

loading↗

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↗

Testing the Preliminary X-33 Navigation System

The X-33 Reusable Launch Vehicle (RLV) must meet the demanding requirements of landing autonomously on a narrow landing strip following a flight that reaches an altitude of up to 200,000 feet and a speed in excess of Mach 9 with significant in-flight energy bleed-off maneuvers. To execute this flight regimen a highly reliable avionics system has been designed that includes three LN-100G Inertial Navigation System/Global Positioning System (INS/GPS) units as the primary navigation system for the X-33. NASA's Marshall Space Flight Center (MSFC) tested an INS/GPS system in real-time simulations to determine the ability of this navigation suite to meet the in flight and autonomous landing requirements of the X-33 RLV. A total of sixty-one open loop tests were performed to characterize the navigation accuracy of the LN-100G. Twenty-seven closed-loop tests were also performed to evaluate the performance of the X-33 Guidance, Navigation and Control (GN&C) algorithms with the real navigation hardware. These closed-loop tests were also designed to expose any integration or operational issues with the real-time X-33 vehicle simulation. Dynamic road tests of the INS/GPS were conducted by Litton to assess the performance of differential and nondifferential INS/GPS hybrid navigation solutions. The results of the simulations and road testing demonstrate that this novel solution is capable of meeting the demanding requirements of take-off, in-flight navigation, and autonomous landing of the X-33 RLV. This paper describes the test environment developed to stimulate the LN-100G and discusses the results of this test effort. This paper also presents recommendations for a navigation system suitable to an operational RLV system.

Lomas, James J.↗

Simulation study of the lift roll coupling problem for hovering VTOL aircraft

The effects of lift-roll coupling on the handling qualities of hovering VTOL aircraft using differential thrust for roll control were assessed in a piloted simulation study in the Ames six-degrees-of-freedom motion simulator. The configuration tested has three vertical thrust sources, one on the roll axis and two laterally displaced from the roll axis, with a thrust distribution of 25%/50%/25%. The outboard thrust sources were modulated to provide roll control whereas all three provided height control. Maximum thrust-to-weight ratio was varied together with a coupling parameter that combined roll-inertia, weight, and engine location. Results showed that handling qualities are affected not only by the occurrence of lift-roll coupling (dependent on both variables) but also by the severity of the coupling (dependent on the coupling parameter alone). However, the advantages of differential thrust for control can be retained with careful design.

Jacobsen, R. A.↗

Evolution of planetesimals. I - Dynamics: Relaxation in a thin disk. II - Numerical simulations

The study examines the effects of density inhomogeneity and differential rotation as well as inelastic collisions on the dynamical evolution of planetesimals. Consideration is given to a three-step analysis: the dynamical evolution of the planetesimals, collisions and mass accumulation, and interaction with gas. It is shown that the velocity dispersion of a cold system of planetesimals increases rapidly due to elastic gravitational scattering. When the dispersion in the epicycle amplitude becomes comparable to the planetesimals' Roche radius, energy is transferred from the systematic Keplerian shear to the dispersive motion. With a numerical N-body scheme, gravitational scattering and physical collisions among a system of planetesimals is simulated. It is shown that dynamical equilibrium is attained with a velocity dispersion comparable to the surface escape velocity of those planetesimals which contribute most of the system mass.

Palmer, P. L.↗

Scientific machine learning for closure models in multiscale problems: A review

Here, closure problems are omnipresent when simulating multiscale systems, where some quantities and processes cannot be fully prescribed despite their effects on the simulation's accuracy. Recently, scientific machine learning approaches have been proposed as a way to tackle the closure problem, combining traditional (physics-based) modeling with data-driven (machine-learned) techniques, typically through enriching differential equations with neural networks. This paper reviews the different reduced model forms, distinguished by the degree to which they include known physics, and the different objectives of a priori and a posteriori learning. The importance of adhering to physical laws (such as symmetries and conservation laws) in choosing the reduced model form and choosing the learning method is discussed. The effect of spatial and temporal discretization and recent trends toward discretization-invariant models are reviewed. In addition, we make the connections between closure problems and several other research disciplines: inverse problems, Mori-Zwanzig theory, and multi-fidelity methods. In conclusion, much progress has been made with scientific machine learning approaches for solving closure problems, but many challenges remain. In particular, the generalizability and interpretability of learned models is a major issue that needs to be addressed further.

97 MATHEMATICS AND COMPUTING↗

A computer simulation of the midlatitude plasmasphere and ionosphere

A computer model has been developed to simulate species density, temperature, and plasma flow in the ionosphere and plasmasphere. A new approach, the flux preserving method involving a Newton iteration, is used to solve a system of governing equations comprising four second-order partial differentials. Tests are performed to demonstrate that the simulation converges to a stable steady state solution, then steady state ion fluxes are analyzed. Finally, simulations of the collapse of the sunset ionosphere are presented. Comparisons with satellite and radar data show good agreement in a number of cases.

Young, E. R.↗

Analyzing the Effects of Atmospheric Turbulent Fluctuations on the Wake Structure of Wind Turbines and Their Blade Vibrational Dynamics

In recent trends, a rising demand for renewable energy has driven wind turbines to larger proportions, where lighter blade designs are often adopted to reduce the costs associated with logistics and production. This causes modern utility-scale wind turbine blades to be inherently more flexible, and their amplified aeroelastic sensitivity results in complex multi-physics reactions to variant atmospheric conditions, including dynamic patterns of aerodynamic loading at the rotor and vortex structure evolutions within the wake. In this paper, we analyze the influence of inflow variance for wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex wake interactions for a range of flexible NRT blade variations, operating in differing conditions of variant atmospheric flow. CODEF solutions of turbine operation in Steady-In-The-Average (SITA) wind conditions are compared to SITA wind conditions featuring a controlled gust-like pulse overimposed, to isolate the effects of typical wind fluctuations. Finally, simulations of realistic time-varying wind conditions from SWiFT meteorological tower measurements are compared to the solutions of SITA wind conditions. These increasingly complex atmospheric inflow variations are tested to show the differing effects evoked by various patterns of spatiotemporal atmospheric flow fluctuations. An analysis is presented for solutions of wind turbine aeroelastic response and vortex wake evolution, to elucidate the consequences of variant inflow, which pertain to wind turbine dynamics at an individual and farm-collective scale. The comparisons of simulated farm flow for SITA and measured fluctuating wind conditions show that certain regions of the wake contain up to a 12% difference in normalized axial velocity, due to the introduction of wind fluctuations. The findings of this study prove valuable for practical applications in wind farm control and optimization strategies, with particular significance for modern utility-scale wind power plants operating in variant atmospheric conditions.

Farrell, Alayna (ORCID:000000023555720X)↗

Modeling and Analysis of Wind Turbine Wake Vortex Evolution Due to Time-Constant Spatial Variations in Atmospheric Flow

Modern utility-scale wind turbines are evolving toward larger, lighter, and more flexible designs to meet the growing demand for renewable energy while minimizing logistical costs. However, these advancements in lightweight design result in heightened aeroelastic sensitivity, leading to complex interactions which affect the rotor’s capacity to withstand aerodynamic loading and the cascading effects that manifest in the wake’s vortex-structure evolution under variable atmospheric conditions. In this paper, we analyze the influence of stream-wise fluctuating atmospheric flow conditions on wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex–wake interactions for a range of conditions with spatially variant atmospheric flow. CODEF solutions for turbine operation in wind conditions featuring only one parameter fluctuation are compared to wind conditions with several wind parameter variations in combination. By isolating individual inflow variations and comparing them to multi-parameter scenarios, we determine the contributions of each atmospheric factor to rotor dynamics, wake evolution, and downstream wind farm interactions. The purpose of this paper is to analyze the effects of spatial variations in atmospheric flow on the topological evolution of wind turbine vortex wakes, which constitutes a gap in the current understanding of wind turbine wake dynamics. The insights gained from this study are particularly valuable for the development of wind farm control strategies aimed at mitigating the adverse effects of wake interactions, enhancing energy capture, and improving the overall stability of wind farm operations. With these insights, we aim to contribute to the development of modeling and simulation tools to optimize utility-scale wind power plants operating in diverse atmospheric environments.

Farrell, Alayna (ORCID:000000023555720X)↗

Coronal mass ejections and the evolution of the large-scale corona

Coronal mass ejections (CMEs) are dynamic, large-scale events in the solar corona that expel plasma and magnetic fields into the solar wind. The structure of the large-scale corona influences CME occurrence; for example, CMEs are most frequently observed to be disruptions of coronal or helmet streamer configurations on the Sun. We have investigated the evolution of the large-scale corona in the presence of differential rotation, using time-dependent magnetohydrodynamic (MHD) simulations in two and three dimensions. An equilibrium coronal configuration is first computed by superimposing a Parker solar wind flow on an initial potential field, the equations are integrated in time until a steady-state is reached. When differential rotation is introduced, we find that the configuration does not reach a steady state; the closed field regions (helmet streamers) disrupt recurrently, ejecting plasmoids into the solar wind. Our results suggest that differential rotation may be one of the mechanisms by which mass ejections are initiated. We will describe the evolution of the polarization brightness (pB) for the resulting configurations, as well as the propagation of the ejected plasmoids in the solar wind, and we will compare the effects of differential rotation with other possible initiation mechanisms, such as the emergence of new magnetic flux.

Linker, J. A.↗

Organic Evaporation and Oxidation Testing in Support of Hanford Sample-and-Send

The Hanford site has approximately 56 million gallons of radioactive mixed waste stored in 177 underground storage tanks. The Hanford Waste Treatment and Immobilization Plant (WTP) is being built to treat and immobilize the tank waste. The baseline method for immobilization of Low Activity Waste (LAW) through the WTP is vitrification, but additional immobilization capacity is needed to supplement the initial LAW melters. An alternative cementitious waste form is being investigated for that future immobilization method to supplement vitrification. However, one impediment to a cementitious waste form is the presence of Land Disposal Restricted (LDR) organic chemicals in tank waste. This work evaluates potential avenues to eliminate that impediment to permit possible use of a cementitious waste form and work towards a decision whether additional LDR organic pretreatment would be required. Savannah River National Laboratory (SRNL) performed testing using simulants to examine evaporation as a method to remove some prevalent organics from LAW. Spiking the caustic LAW simulant with selected regulated organic chemicals found one that clearly decomposes because of caustic instability. Oxidation testing of other organic chemicals found some LDR organics degrade as desired and others are stable in the presence of peroxide and permanganate. In addition to studies with simulants, a literature review was performed to evaluate radiological stability of LDR organics. Descriptions of the experimental details, equipment, and results are included in this report. Evaporation testing consisted of preparing the LAW simulant, spiking that simulant with organic chemicals, and evaporating the mixture via differential distillation. The apparatus was a laboratory-scale vacuum evaporator operated at 60 ±5 torr absolute (vacuum evaporation) and also at atmospheric pressure. The LAW simulant represented the liquid expected to be retrieved from the Hanford tank farms at approximately 4.0 M [Na + ] total sodium ion concentration. The concentration of the organic chemicals added was significantly higher than typically found in the tank waste samples since the higher levels were necessary to assist in analytical measurement and tracking of the spiked species. Organic chemicals were chosen for the work with a consideration of how their volatility compares with that of methanol. This was done by comparing the ratio of the pure water Henry’s law coefficient (K h ) of methanol to that of the compound in question (hereafter termed the K h ratio), where ratios above unity indicated less volatility than methanol. Methanol was chosen because it is a common regulated chemical with relatively low volatility but which has been removed by evaporation in previous laboratory work. While organic separation results depend on evaporator design, laboratory experiments verified that organic partitioning to the overhead condensate stream by evaporation is a practical process. The work reported here found difficulties in quantitative analysis of the organic chemicals in aqueous samples. Most of the time there was insufficient analysis to close a mass balance for evaporator runs, but qualitative evidence of carryover was obtained. The methods were also able to show whether organic chemicals were susceptible or resistant to solution oxidation in permanganate or hydrogen peroxide tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Field Work Proposal ERKJ358: Black-box training for scientific machine learning models (Final Report)

The overarching goal of this project is to develop a scalable black-box training capability for scientific machine learning (SciML) problems that are non-trainable with existing automatic differentiation (AD)-based algorithms. AD assumes that a loss function can be decomposed into a sequence of elementary operations whose derivatives are known. This assumption is violated when the loss function includes a black-box physical model (e.g., a legacy simulator). The current strategy, converting a black-box simulator to an AD-enabled code via differential programming, is inflexible and time-, labor-consuming. Thus, black-box optimization is a main workhorse for training SciML models, e.g., in scientific reinforcement learning, hyper-parameter fine tuning, designing SciML models with adversarial robustness, etc.

97 MATHEMATICS AND COMPUTING↗

The condensational instability in the solar transition region and corona

The stability of plasmas at temperatures and densities typical of the solar transition region and corona is investigated using both a linear analysis and nonlinear time-dependent numerical simulations. Growth rates, decay rates, and oscillation frequencies of the perturbations determined from the linear analysis are in excellent agreement with the simulations. The nonlinear regime is characterized by a bifurcation of the plasma into a cool dense condensation surrounded by a hot tenuous corona. The condensation may then be accelerated by forces in the plasma such as those arising from gravity or differential heating. Finally, the results of the detailed simulation show that the transition region is a dynamically stable structure which is the result of the nonlinear evolution of the condensational instability.

Oran, E. S.↗

Phase-field model of precipitation processes with coherency loss

A phase-field model is proposed to simulate coherency loss coupled with microstructure evolution. A special field variable is employed to describe the degree of coherency loss of each particle and its evolution is governed by a Ginzburg-Landau type kinetic equation. For the sake of computational efficiency, a flood-fill algorithm is introduced that can drastically reduce the required number of field variables, which allows the model to efficiently simulate a large number of particles sufficient for characterizing their statistical features during Ostwald ripening. The model can incorporate size dependence of coherency loss, metastability of coherent particles, and effectively incorporate the underlying mechanisms of coherency loss by introducing a so-called differential energy criterion. The model is applied to simulate coarsening of Al 3 Sc precipitates in aluminum alloy and comprehensively compared with experiments. Our results clearly show how the particle size distribution is changed during coherency loss and affects the coarsening rate.

36 MATERIALS SCIENCE↗

Marmot V2

MARMOT is a robust numerical tool for mesoscale modeling of fuel performance developed under the NEAMS Fuels technical area to predict the coevolution of microstructure and properties in fuel and cladding materials. MARMOT accomplishes this using the phase field method coupled with finite strain mechanics and heat conduction. MARMOT is based on the open source Multiphysics Object-Oriented Simulation Environment (MOOSE) and solves the coupled partial differential equations defining the physics using the finite element method. MARMOT is being developed in order to facilitate the development of improved materials models for fuel performance, but it is also being developed as a powerful tool in and of itself for the simulation of mesoscale fuel performance.

Aagesen, LarryK.↗

A multi-dimensional parametric study of variability in multi-phase flow dynamics during geologic CO 2 sequestration accelerated with machine learning

Successful geologic CO 2 storage projects depend on numerical simulations to predict reservoir performance during site selection, injection verification, and post-injection monitoring phases of the project. These numerical simulations solve non-linear sets of coupled partial differential equations, while accounting for multi-phase fluid dynamics on the basis of constitutive equations that are embedded into the solution scheme. As a consequence, individual simulations often require tens to hundreds of hours to complete on high-performance computing clusters. Moreover, laboratory experiments reveal that parametric functions for capillary pressure and relative permeability exhibit substantial variability, even within the same rock type. This combination of computational expense and wide-ranging parametric variability means that there remains substantial uncertainty in the behavior of multi-phase CO 2 -water systems, particularly in the context of feedbacks between relative permeability and capillary pressure. To bridge this knowledge gap, here we develop a novel workflow that utilizes physics-based numerical simulation to train an artificial neural network (ANN) emulator for interrogating the multivariate parameter space that governs both capillary pressure and relative permeability. With this approach, the ANN is trained to emulate both fluid pressure distribution and CO 2 saturation, which are then interrogated quantitatively to generate parametric response surface mappings with high-fidelity resolution. Results from this study initially show that capillary entry pressure is the dominant control on both CO 2 plume geometry and fluid pressure propagation when considering the combined effects of capillary pressure and relative permeability, particularly when phase interference is low and residual CO 2 saturation is high. Moreover, the ANN emulator provides tremendous computational speed-up by computing 2691 individual simulations in several minutes; whereas, the same simulation ensemble would have required ~3 years of simulation time using only physics-based simulation methods (25,000 times speed up).

58 GEOSCIENCES↗

Study of fully coupled three-dimensional envelope instability using automatic differentiation

Automatic differentiation is a powerful tool for computing derivatives of simulation results with respect to given parameters. In this Letter, we have applied this tool to investigate the instability of a dynamical system governed by 21 ordinary differential equations. This second-order instability (named envelope instability) is driven by space-charge effects and has a significant impact on the operational regimes of particle accelerators. Our study delves into the three-dimensional envelope instability, incorporating both transverse and longitudinal coupling. Conventionally, analyzing this complex system would necessitate solving 441 ordinary differential equations, which is computationally intractable. However, by employing automatic differentiation, we were able to track only 21 equations. This approach allowed us to uncover an additional instability stopband, which arises from space-charge-induced coupling and has not been reported in previous studies. This research highlights the significant advantages of automatic differentiation in analyzing complicated dynamical systems involving a large number of ordinary differential equations.

Qiang, Ji [Lawrence Berkeley National Laboratory (↗