SEARCH · Engineering Papers

Results for “Simplified Simulator”

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 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Validation of wind resource and energy production simulations for small wind turbines in the United States

Abstract. Due to financial and temporal limitations, the small wind community relies upon simplified wind speed models and energy production simulation tools to assess site suitability and produce energy generation expectations. While efficient and user-friendly, these models and tools are subject to errors that have been insufficiently quantified at small wind turbine heights. This study leverages observations from meteorological towers and sodars across the United States to validate wind speed estimates from the Wind Integration National Dataset (WIND) Toolkit, the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis v5 (ERA5), and the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2), revealing average biases within ±0.5 m s−1 at small wind hub heights. Observations from small wind turbines across the United States provide references for validating energy production estimates from the System Advisor Model (SAM), Wind Report, MyWindTurbine.com, and Global Wind Atlas 3 (GWA3), which are seen to overestimate actual annual capacity factors by 2.5, 4.2, 11.5, and 7.3 percentage points, respectively. In addition to quantifying the error metrics, this paper identifies sources of model and tool discrepancies, noting that interannual fluctuation in the wind resource, wind speed class, and loss assumptions produces more variability in estimates than different horizontal and vertical interpolation techniques. The results of this study provide small wind installers and owners with information about these challenges to consider when making performance estimates and thus possible adjustments accordingly. Looking to the future, recognizing these error metrics and sources of discrepancies provides model and tool researchers and developers with opportunities for product improvement that could positively impact small wind customer confidence and the ability to finance small wind projects.

17 WIND ENERGY↗

DOE OSTI2026

Benchmarking of massively parallel phase-field codes for directional solidification

We present a detailed benchmark comparing two state-of-the-art phase-field implementations for simulating alloy solidification under experimentally relevant conditions. The study investigates the directional solidification of Al-3wt%Cu under high-velocity solidification conditions and SCN-0.46wt% camphor under microgravity conditions from National Aeronautics and Space Administration (NASA) DECLIC-DSI-R experiments. Both codes, one employing finite-difference discretization with uniform mesh and GPU-acceleration (GPU-PF) and the other one employing finite-element discretization with adaptive-mesh and CPU-parallelization (PRISMS-PF), solve the same quantitative phase-field formulation that incorporates an anti-trapping current for the solidification of dilute alloys. We evaluate the predictions of each code for dendritic morphology, primary spacing, and tip dynamics in both 2D and 3D, as well as their numerical convergence and computational performance. While existing benchmark problems have primarily focused on simplified or small-scale simulations, they do not reflect the computational and modeling challenges posed by employing experimentally relevant time and length scales. Our results provide a practical framework for assessing phase-field code performance as well as validating and facilitating their application in integrated computational materials engineering (ICME) workflows that require integration with realistic experimental data.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗

NASA NTRS2019

Hybrid Simulation of Solar-Wind-Like Turbulence

We present 2.5D hybrid simulations of the spectral and thermodynamic evolution of an initial state of magnetic field and plasma variables that in many ways represents solar wind fluctuations. In accordance with Helios near-Sun high-speed stream observations, we start with Alfvénic fluctuations along a mean magnetic field in which the fluctuations in the magnitude of the magnetic field are minimized. Since fluctuations in the radial flow speed are the dominant free energy in the observed fluctuations, we include a field-aligned v∥(k⊥) with an k(exp −1) spectrum of velocity fluctuations to drive the turbulent evolution. The flow rapidly distorts the Alfvénic fluctuations, yielding spectra (determined by spacecraft-like cuts) transverse to the field that become comparable to the k∥ fluctuations, as in spacecraft observations. The initial near constancy of the magnetic field is lost during the evolution; we show this also takes place observationally. We find some evolution in the anisotropy of the thermal fluctuations, consistent with expectations based on Helios data. We present 2D spectra of the fluctuations, showing the evolution of the power spectrum and cross-helicity. Despite simplifying assumptions, many aspects of simulations and observations agree. The greatly faster evolution in the simulations is at least in part due to the small scales being simulated, but also to the non-equilibrium initial conditions and the relatively low overall Alfvénicity of the initial fluctuations.

Roberts, D. Aaron↗

NASA NTRS2012

Soil Moisture Active/Passive (SMAP) Forward Brightness Temperature Simulator

The SMAP is one of four first-tier missions recommended by the US National Research Council's Committee on Earth Science and Applications from Space (Earth Science and Applications from Space: National Imperatives for the Next Decade and Beyond, Space Studies Board, National Academies Press, 2007) [1]. It is to measure the global soil moisture and freeze/thaw from space. One of the spaceborne instruments is an L-band radiometer with a shared single feedhorn and parabolic mesh reflector. While the radiometer measures the emission over a footprint of interest, unwanted emissions are also received by the antenna through the antenna sidelobes from the cosmic background and other error sources such as the Sun, the Moon and the galaxy. Their effects need to be considered accurately, and the analysis of the overall performance of the radiometer requires end-to-end performance simulation from Earth emission to antenna brightness temperature, such as the global simulation of L-band brightness temperature simulation over land and sea [2]. To assist with the SMAP radiometer level 1B algorithm development, the SMAP forward brightness temperature simulator is developed by adapting the Aquarius simulator [2] with necessary modifications. This poster presents the current status of the SMAP forward brightness simulator s development including incorporating the land microwave emission model and its input datasets, and a simplified atmospheric radiative transfer model. The latest simulation results are also presented to demonstrate the ability of supporting the SMAP L1B algorithm development.

Peng, Jinzheng↗

NASA NTRS2016

Simplified ISCCP Cloud Regimes for Evaluating Cloudiness in CMIP5 Models

We take advantage of the ISCCP simulator implemented in many CMIP5 models, in order to introduce a credible method to compare model cloud output with corresponding ISCCP observations based on the cloud regime (CR) concept. A novel simplified CR derivation method is introduced that exploits the co-variations of three variables, cloud optical thickness, cloud top pressure and cloud fraction (, pc, CF). Following evaluation criteria established in a previous paper of ours, we assess model cloud simulation performance based on how well the simplified CRs are simulated in terms of similarity of centroids, global values and map correlations of relative frequency-of-occurrence (RFO), and long term total clouds amounts (as opposed to average instantaneous cloud fractions). Mirroring prior results, modeled clouds tend to be too thick optically and not as extensive as in observations. Averaging across models does not help in this or other evaluation aspects, except for cloud geographical locations. A good model for one criterion is usually found to be also a good performer for other criteria despite plenty of exceptions for the criterion of centroid similarity. On the other hand, CRs that include convectively generated high-altitude clouds are not as well simulated here compared to the previous study, but other regimes containing low clouds with high CF improve. Models that have previously performed well distinguish themselves again here, but improvements from other models are also found. The simplified cloud evaluation method thus proves to be more forgiving while still good enough to reveal model weaknesses.

global climate models↗

DOE OSTI2022

LANL Contributions to pyiron: An integrated development environment for computational materials science

To support and accelerate the development of simulation protocols in atomistic modelling, we introduce an integrated development environment (IDE) for computational materials science called pyiron (http://pyiron.org). The pyiron IDE combines a web based source code editor, a job management system for build automation, and a hierarchical data management solution. The core components of the pyiron IDE are pyiron objects based on an abstract class, which links application structures such as atomistic structures, projects, jobs, simulation protocols and computing resources with persistent storage and an interactive user environment. The simulation protocols within the pyiron IDE are constructed using the Python programming language. To highlight key concepts of this tool as well as to demonstrate its ability to simplify the implementation and testing of simulation protocols we discuss two applications. In these examples we show how pyiron supports the whole life cycle of a typical simulation, seamlessly combines ab initio with empirical potential calculations, and how complex feedback loops can be implemented. While originally developed with focus on ab initio thermodynamics simulations, the concepts and implementation of pyiron are general thus allowing to employ it for a wide range of simulation topics.

Janssen, Jan↗

DOE OSTI2022

Role of electronic energy loss on defect production and interface stability: Comparison between ceramic materials and high-entropy alloys

High-entropy alloys (HEAs) and some complex alloys exhibit desirable properties and significant structural stability in harsh environments, including possible applications in advanced reactors. Energetic ion irradiation is often used as a surrogate for neutron irradiation; however, the impact of ion electronic energy deposition and dissipation is often neglected. Moreover, differences in recoil energy spectrum and density of cascade events on damage evolution must also be considered. In many chemically complex alloys, the mean free path of electrons is reduced significantly, thus their decreased thermal conductivity and slow dissipation of localized radiation energy can have noticeable effects on displacement cascade evolution that is greatly different from metals with high thermal conductivity. In this work, nanocrystalline HEAs of Ni 20 Fe 20 Co 20 Cr 20 Cu 20 and nonequiatomic (NiFeCoCr) 97 Cu 3 , both having much lower room-temperature thermal conductivity than pure Ni or Fe, are chosen as model HEAs to reveal the role that electronic energy loss during ion irradiation has in complex alloys. The response of nanocrystalline HEAs is investigated under irradiation at room temperature using MeV Ni and Au ions that have different ratios of electronic energy to damage energy, which is the energy dissipated in displacing atoms. Different from previously reported amorphization of nanocrystalline SiC, experimental results on these HEAs show that, similar to the process in nanocrystalline oxide materials, both inelastic thermal spikes via electron–phonon coupling and elastic thermal spikes via collisions among atomic nuclei contribute to the overall grain growth. The growth follows a power law dependence with the total deposited ion energy, and the derived value of the power-exponent suggests that the irradiation-induced instability at and near grain boundaries leads to local rapid atomic rearrangements and consequently grain growth. The high power-exponent value can be attributed to the sluggish diffusion and delayed defect evolution arising from the chemical complexity intrinsic to HEAs. Here, this work calls attention to quantified fundamental understanding of radiation damage processes beyond that of simplified displacement events, especially in simulating neutron environments.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Structural and Dynamical Roles of Bound Polymer Chains in Rubber Reinforcement

The addition of nanofillers to rubber matrices is a powerful route to improve the mechanical properties. Here, we focus on a molecular understanding of basic mechanisms that are important for the reinforcement in rubbers. The key role in this process is ascribed to bound rubber (BR) that engages with the matrix as well as with adjacent nanofillers. To date, this understanding has been impeded by the lack of experimental tools to directly probe the BR chains buried in a polymer matrix composed of the same polymer. To tackle this challenge, we combine neutron scattering/spectroscopy techniques with isotope-labeling and molecular dynamics simulations. The system is a simplified carbon-black-filled polybutadiene. The combined experimental and computational results provide new insights into the local structural and dynamical heterogeneities of BR chains and their interactions with the matrix polymer, highlighting (i) the structural partition of the bound chains into three components (i.e., trains, loops, and tails) and their fractions; (ii) their dynamical hierarchies, i.e., the trains that remain immobile on the filler surface, the loops that are fairly large and hence allow the interdigitation of matrix chains, and the tails with their unique characteristics to reach far out into the matrix and entangle with matrix chains. These multiple roles of the constituent components of the BR chains promote the formation of a well-developed adhesive polymer–filler interface, enhancing the elastic property of a filled rubber. Finally, the comprehensive understanding derived and validated by the model rubber will be translatable to many other polymer nanocomposites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Risk-Constrained Reinforcement Learning for Inverter-Dominated Power System Controls

Here, this paper develops a risk-aware controller for grid-forming inverters (GFMs) to minimize large frequency oscillations in GFM inverter-dominated power systems. To tackle the high variability from loads/renewables, we incorporate a mean-variance risk constraint into the classical linear quadratic regulator (LQR) formulation for this problem. The risk constraint aims to bound the time-averaged cost of state variability and thus can improve the worst-case performance for large disturbances. The resulting risk-constrained LQR problem is solved through the dual reformulation to a minimax problem, by using a reinforcement learning (RL) method termed as stochastic gradient-descent with max-oracle (SGDmax). In particular, the zero-order policy gradient (ZOPG) approach is used to simplify the gradient estimation using simulated system trajectories. Numerical tests conducted on the IEEE 68-bus system have validated the convergence of our proposed SGDmax for GFM model and corroborate the effectiveness of the risk constraint in improving the worst-case performance while reducing the variability of the overall control cost.

Frequency control↗

DOE OSTI2024

Cyberwheel

Cyberwheel is a high fidelity training environment for autonomous cyber defense agents that includes both simulation and emulation capabilities. Cyberwheel simplifies customization of the training network and easily allows redefining the agent’s reward function, observation space, and action space to support rapid experimentation of novel approaches to agent design. It also provides visibility into agent behaviors necessary for agent evaluation and sufficient documentation and examples to lower the barrier to entry.

Oesch, TimothySean [Oak Ridge National Laboratory ↗

DOE OSTI2021

Hybrid (PDE+ML) models in the context of land ice modeling

Focal Area(s): Primary focal areas are: predictive modeling through the use of AI-derived model components; advanced methods including network design/optimization/deep learning. Science Challenge: The atmospheric, ocean and ice dynamics components of the Energy Exascale Earth System Model (E3SM) are governed by Partial Differential Equations (PDEs) and significant efforts have been made during the last decades to develop such computational models. Here we propose to fundamentally improve these PDE-based codes by enhancing them with Machine Learning (ML) sub-models for complex, poorly understood physical processes in the context of ice sheet modeling. We propose to train these models with a novel approach that allows the assimilation of the different sources of data available (direct/indirect observations and possibly simulation data), improving on existing simplified models. We also highlight computational challenges originating from the coexistence of PDE-based and ML-based models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Improving Common PV Module Temperature Models by Incorporating Radiative Losses to the Sky

PV module operating temperature is the second-most important factor influencing PV system yield–after irradiance–and a substantial contributor to uncertainty in energy system yield predictions. Models commonly used to predict operating temperature in system simulations are based on a simplified energy balance that lumps together different heat loss mechanisms–including radiation–and assumes an overall linear behavior. Radiative heat loss to the sky is usually substantial, but modeling it accurately requires additional information about down-welling long-wave radiation or sky temperature and increases the complexity of temperature model equations. In this work we show how radiative losses to the sky can be separated into two parts to improve the accuracy of modeling without additional complexity. We also predict and demonstrate the variation of these losses at different tilt angles and show that the effective view factor is reduced by the non- isotropic distribution of down-welling long-wave radiation. Finally, we demonstrate substantial reduction in bias (MBE) and scatter (RMSE) when the new radiative loss term is added to the Faiman model using one year of measurements at Sandia National Labs.

14 SOLAR ENERGY↗

NASA NTRS1990

Structure and equilibrium of coronal magnetic fields

In 'closed' magnetic structures (i.e., coronal loops) the random shuffling of magnetic footpoints in the photosphere causes twisting and braiding of field lines in the corona. If the motions are sufficiently slow, the coronal field evolves through a sequence of force-free equilibrium states. Numerical simulations are presented for a simplified model in which the overall curvature of the coronal loop is neglected. It is shown that magnetic fine structures develop on spatial scales significantly smaller than those of the imposed 'photospheric' velocity field.

Van Ballegooijen, A. A.↗

NASA NTRS1993

Approximate optimal guidance for the advanced launch system

A real-time guidance scheme for the problem of maximizing the payload into orbit subject to the equations of motion for a rocket over a spherical, non-rotating earth is presented. An approximate optimal launch guidance law is developed based upon an asymptotic expansion of the Hamilton - Jacobi - Bellman or dynamic programming equation. The expansion is performed in terms of a small parameter, which is used to separate the dynamics of the problem into primary and perturbation dynamics. For the zeroth-order problem the small parameter is set to zero and a closed-form solution to the zeroth-order expansion term of Hamilton - Jacobi - Bellman equation is obtained. Higher-order terms of the expansion include the effects of the neglected perturbation dynamics. These higher-order terms are determined from the solution of first-order linear partial differential equations requiring only the evaluation of quadratures. This technique is preferred as a real-time, on-line guidance scheme to alternative numerical iterative optimization schemes because of the unreliable convergence properties of these iterative guidance schemes and because the quadratures needed for the approximate optimal guidance law can be performed rapidly and by parallel processing. Even if the approximate solution is not nearly optimal, when using this technique the zeroth-order solution always provides a path which satisfies the terminal constraints. Results for two-degree-of-freedom simulations are presented for the simplified problem of flight in the equatorial plane and compared to the guidance scheme generated by the shooting method which is an iterative second-order technique.

Feeley, T. S.↗

NASA NTRS1995

Modeling the decay of energy containing eddies: A source of solar wind heating

To understand the solar wind heating and acceleration mechanisms one needs to understand the decay of energy containing eddies. With this goal in mind, attempts have been made to extend the fluid dynamic phenomenology of large scale quasi-equilibrium to the case of magnetohydrodynamics. Matthaeus et al. have proposed a model for the inhomogeneous transport and decay of five mean variables, namely, two mean square Elsasser variables z(exp 2) (sub +/-) their correlation lengths, and the difference between the kinetic and magnetic energies. We test the validity of this model in the simplified case of homogeneous turbulence simulated in a periodic box. We propose a class of models and show that they may fit the simulation satisfactorily. Analytic solutions of this class of model reveal their inherent properties and demonstrate the difficulties associated with finite cross helicity. It is noted that adjustments are required to make the simplest models, which are based upon isotropic turbulence, scale properly with respect to the strength of the mean magnetic field. This can be interpreted as due to anisotropic turbulence, which can be modelled by simple parameterization in the phenomenology.

Hossain, M.↗

NASA NTRS1996

LES on unstructured deforming meshes: Towards reciprocating IC engines

A variable explicit/implicit characteristics-based advection scheme that is second-order accurate in space and time has been developed recently for unstructured deforming meshes (O'Rourke & Sahota 1996a). To explore the suitability of this methodology for Large-Eddy Simulation (LES), three subgrid-scale turbulence models have been implemented in the CHAD CFD code (O'Rourke & Sahota 1996b): a constant-coefficient Smagorinsky model, a dynamic Smagorinsky model for flows having one or more directions of statistical homogeneity, and a Lagrangian dynamic Smagorinsky model for flows having no spatial or temporal homogeneity (Meneveau et al. 1996). Computations have been made for three canonical flows, progressing towards the intended application of in-cylinder flow in a reciprocating engine. Grid sizes were selected to be comparable to the coarsest meshes used in earlier spectral LES studies. Quantitative results are reported for decaying homogeneous isotropic turbulence, and for a planar channel flow. Computations are compared to experimental measurements, to Direct-Numerical Simulation (DNS) data, and to Rapid-Distortion Theory (RDT) where appropriate. Generally satisfactory evolution of first and second moments is found on these coarse meshes; deviations are attributed to insufficient mesh resolution. Issues include mesh resolution and computational requirements for a specified level of accuracy, analytic characterization of the filtering implied by the numerical method, wall treatment, and inflow boundary conditions. To resolve these issues, finer-mesh simulations and computations of a simplified axisymmetric reciprocating piston-cylinder assembly are in progress.

Haworth, D. C.↗