SEARCH · Search NASA

Results for “engine 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 1,459 records · Page 81

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Impact of mechanical tolerances on partial turn skipping in helical flux compression generator

Maintaining tight mechanical tolerances of the components used in helical flux compression generators is crucial for optimizing device performance. Moving beyond current literature, which provides various approximations for acceptable tolerances, a 3D simulation of the armature stator interaction during generator operation was developed, revealing that existing equations frequently misestimate the tolerance limits. Without restricting the generality of the presented approach, the focus was primarily on the armature’s concentricity and roundness, comparing the acceptable tolerance limits between generators of different sizes. An analytical approach for select geometries was developed, confirming the 3D results. Quantitative results reveal that a steeper Gurney angle and increased winding pitch allow for less restrictive tolerances for eccentricity and ellipticity. Even small deviations from the ideal armature positioning and shape will result in fluctuations of the output current time-derivative, observed experimentally and in simulation. Larger deviations then cause partial turn skipping, associated with loss of magnetic flux. The simulation elucidates the impact of tolerance deviations on generator performance by offering precise tolerance ranges for various generator sizes, thereby facilitating informed design decisions.

42 ENGINEERING↗

DOE OSTI2019

Development, Demonstration and Validation of Data-Driven Compact Diode Models for Circuit Simulation and Analysis

Compact semiconductor device models are essential for efficiently designing and analyzing large circuits. However, traditional compact model development requires a large amount of manual effort and can span many years. Moreover, inclusion of new physics (e.g., radiation effects) into an existing model is not trivial and may require redevelopment from scratch. Machine Learning (ML) techniques have the potential to automate and significantly speed up the development of compact models. In addition, ML provides a range of modeling options that can be used to develop hierarchies of compact models tailored to specific circuit design stages. In this paper, we explore three such options: (1) table-based interpolation, (2) Generalized Moving Least-Squares, and (3) feedforward Deep Neural Networks, to develop compact models for a p-n junction diode. We evaluate the performance of these "data-driven" compact models by (1) comparing their voltage-current characteristics against laboratory data, and (2) building a bridge rectifier circuit using these devices, predicting the circuit's behavior using SPICE-like circuit simulations, and then comparing these predictions against laboratory measurements of the same circuit.

1N4148 switching diode↗

DOE OSTI2016

Benchmarking Initiatives

Benchmarking has been recognized as an efficient tool for model validation as well as scientific collaboration. Several workshop series have been established to foster the benchmarking idea: setting up test cases with increasing complexity for method development and code comparison. In addition to representing model complexity, one of the key efforts is to develop codes that are suitable for modern HPC platforms such as PetaFlop supercomputers. It has been realized that those challenges are beyond single team capabilities. The Subsurface Environmental Simulation Benchmarking (SSBench) working group was initiated to develop, publish, and make publicly available well-described benchmark problems that could be used to demonstrate simulator conformance with norms established by the subsurface science and engineering community. SSBench supports the development of environmentally-relevant benchmark problem sets that provide rigorous tests for process model conceptualization, numerical implementation, process coupling, and accuracy. The initial SSBench focus on multicomponent reactive transport modeling reflects the expanding range of mechanisms and increasing complexity needed to address issues relevant to climate change, nuclear waste management, contaminant remediation, and pollution prevention.

benchmarking↗

DOE OSTI2025

Toward engineering lattice structures with the material point method (MPM)

This study examines the potential of two variants of the material point method—the generalized interpolation material point (GIMP) and dual domain material point (DDMP) methods—in developing a robust computational framework for engineering lattice structures under different loading conditions. The study begins with assessing the ability of the two methods in predicting elastic buckling phenomena using column geometries with and without initial geometric imperfections. The results indicate that both methods effectively capture buckling phenomena when initial geometric imperfections are introduced. After this verification step, we create several models of tetrahedral lattice structures with varying strut diameter and orientation and subject them to quasi-static loading. We then validate the numerical results using laboratory test results. The results show that, while both methods accurately predict load–displacement curves in the pre-buckling regime, their predictive capabilities diminish in the post-buckling regime. Through visual comparison between the numerical and experimental deformed shapes, it appears that the discrepancies between model and experimental results are attributed to initial geometric imperfections in the lattices that occurred during 3D printing. We then establish a second set of lattice models where different types of initial geometric imperfections are considered. The results from these models show that imperfections have a negligible influence in the pre-buckling regime but affect the behavior considerably in the post-buckling regime. As a final step in this work, we subject the lattice models to impact loading and employ hypothetical soft and stiff materials. These results show that the lattice stiffness, which depends on material stiffness, strut diameter, and orientation, significantly influences the ability of a lattice structure to resist impact. In particular, we find that a stiffer lattice (i.e., one made with a stiff material and thicker struts) is capable of absorbing more energy than a softer one during impact. Although material nonlinearities, inelasticity, and detailed contact formulations are not considered in this study, the findings obtained herein lay the groundwork for engineering lattice structures under extreme loading conditions through a simulation-driven framework based on particle-based methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Influence of Fiber Orientation on Deformation of Additive Manufactured Composites

Anisotropy caused by flow-induced fiber orientation of discontinuous fibers within the extrudate in the Extrusion Deposition Additive Manufacturing (EDAM) process gives rise to anisotropic shrinkage in printed parts and thereby, final part deformations. Three fiber orientation states, described by the second-order orientation tensor, were investigated to determine their influence upon final geometry. Two material systems were investigated including a short glass fiber-reinforced polyamide and a short carbon fiber-reinforced polyamide. A curvilinear geometry was modeled and printed virtually in the thermo-mechanical simulation, ADDITIVE3D©. The scale of the printed geometry was in the range of 300–400 mm and contained three semicircular geometry regions connected to linear regions, thereby inducing significant magnification of the spring-in deformation. The predicted deformation of the printed geometry for the three fiber orientation states and two material systems were compared. The deformation predicted for glass fiber-filled polyamide and carbon fiber-filled polyamide of the same fiber orientation states showed comparable deformations. In conclusion, the largest residual deformation was shown to correspond to the orientation tensor with the greatest degree of anisotropy.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Ten questions on future and extreme weather data for building simulation and analysis in a changing climate

Weather plays a significant role in building operations as it directly influences HVAC loads and in turn the building energy and thermal performance. In a changing climate, future trends and extreme weather events become critical concerns in the global building decarbonization and clean energy transition. This paper aims to address ten key questions concerning extreme and future weather data for building applications, and more importantly to identify research gaps and guide the curation and selection of future and extreme weather data for use in building performance simulation and assessment. The paper intends to inform architects and engineers, operators, owners, policy makers, and other stakeholders on considering the impacts of future and extreme weather data and adopting strategies for selecting and applying this data in various use cases related to building design, operation, and retrofit for energy efficiency, electrification, and climate resilience.

Yan, Da↗

DOE OSTI2022

Forward and inverse modeling of fault transmissibility in subsurface flows

Characterizing physical properties of faults, such as their transmissibility, is crucial for performing predictive numerical simulation of subsurface flows, such as those encountered in petroleum engineering and remediation of subsurface contamination. Here, this paper provides a complete investigation of the inverse problem for fault transmissibility in subsurface flow models, under appropriate assumptions on fault structure. In particular, the following aspects are considered: 1) fault modeling and well-posedness of the forward problem; 2) finite element (FEM) discretizations of the forward problem and their rigorous a priori convergence analysis; 3) Well-posedness of the Bayesian inverse problem, FEM discretization of the infinite dimensional Bayesian inverse formulation, and its rigorous a priori analysis. Moreover, computation of the maximum a posteriori (MAP) point via fast inexact Newton-conjugate gradient optimization and a Laplace approximation of the Bayesian posterior are also presented. Numerical results illustrate the use of the proposed fault model in forward and inverse problems for subsurface flows in two dimensional domains with multiple faults.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Algorithms and analyses for stochastic optimization for turbofan noise reduction using parallel reduced-order modeling

Simulation-based optimization of acoustic liner design in a turbofan engine nacelle for noise reduction purposes can dramatically reduce the cost and time needed for experimental designs. Because uncertainties are inevitable in the design process, a stochastic optimization algorithm is posed based on the conditional value-at-risk measure so that an ideal acoustic liner impedance is determined that is robust in the presence of uncertainties. A parallel reduced-order modeling framework is developed that dramatically improves the computational efficiency of the stochastic optimization solver for a realistic nacelle geometry. The reduced stochastic optimization solver takes less than 500 s to execute. Additionally, well-posedness and finite element error analyses of the state system and optimization problem are provided.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Bayesian calibration of strength model parameters from Taylor impact data

Materials strength plays a key role in determining the mechanical response of engineered structures. As such, accurate strength models are crucial in simulations involving complex loading conditions, particularly when deformation in the plastic regime is deemed important. In this work, a Gaussian process based surrogate for a finite element simulation of a Taylor impact test is developed and used for Bayesian calibration of the Preston–Tonks–Wallace strength model. Here, the surrogate model is shown to closely approximate the salient features of the Taylor cylinder deformation and is validated against simulation output before being used in the strength model calibration routine. The results show that Taylor impact test data can be used in the calibration of constitutive equations through the use of a combination of data science techniques, namely Gaussian processes and Bayesian inference.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A parallel and performance portable implementation of a full-field crystal plasticity model

We have developed a parallel implementation of an Elasto-Viscoplastic Fast Fourier Transform-based (EVPFFT) micromechanical solver to enable computationally efficient crystal plasticity modeling for polycrystalline materials. Our primary focus lies in achieving performance portability, allowing a single EVPFFT implementation to run optimally on various homogeneous architectures, including multi-core Central Processing Units (CPUs), as well as on heterogeneous computer architectures comprising multi-core CPUs and Graphics Processing Units (GPUs) from different vendors. To accomplish this goal, we have leveraged MATAR, a C++ software library that simplifies the creation and utilization of multidimensional dense or sparse matrix and array data structures. These data structures are designed to be portable across diverse architectures through the use of Kokkos, a performance-portable library. Additionally, we have employed the Message Passing Interface (MPI) to efficiently distribute the computational workload among processors. The heFFTe (Highly Efficient FFT for Exascale) library is used to facilitate the performance portability of the fast Fourier transforms (FFTs) computation. The computational performance of EVPFFT is evaluated and presented in terms of parallel scalability and simulation runtime on different high-performance computing (HPC) architectures. As a result, the utility of the developed framework to efficiently simulate the micro-mechanical fields in polycrystalline microstructures in engineering applications is discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Compressive sensing adaptation for polynomial chaos expansions

Basis adaptation in Homogeneous Chaos spaces rely on a suitable rotation of the underlying Gaussian germ. Several rotations have been proposed in the literature resulting in adaptations with different convergence properties. Here, we present a new adaptation mechanism that builds on compressive sensing algorithms, resulting in a reduced polynomial chaos approximation with optimal sparsity. The developed adaptation algorithm consists of a two-step optimization procedure that computes the optimal coefficients and the input projection matrix of a low dimensional chaos expansion with respect to an optimally rotated basis. Finally, we demonstrate the attractive features of our algorithm through several numerical examples including the application on Large-Eddy Simulation (LES) calculations of turbulent combustion in a HIFiRE scramjet engine.

42 ENGINEERING↗

DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗

DOE OSTI2023

On the evolution of scalar iso-surface area density in a turbulent mixing layer

In this paper, the transport equation for scalar iso-surface area density (Σ) in a turbulent, temporally developing mixing layer is examined. Exploring the spatial and temporal evolution of the terms in the Σ transport equation is vital to improving our understanding of turbulent flows characterized by distinct interfaces, e.g. the flame surface or the turbulent/non-turbulent interface. Previous work reported by the authors identified that Σ exhibits self-similar behaviour consistent with the development of the temporal mixing layer. Accordingly, each of the terms in the Σ transport equation is found to behave in a self-similar manner, though there are notable differences in the self-similar behaviours for each term. Based on the results presented herein, it is suggested that the rate of change of Σ (∂Σ/∂t) and the advection term scale with hλ Φ /ΔU, where h is the width of the mixing layer, λ Φ is the scalar Taylor length scale and ΔU is the velocity difference. The production and destruction terms are found to scale with an additional factor (ReSc) 1/2 . In contrast, the molecular diffusion term is found to scale with a factor (ReSc) –1/2 compared to ∂Σ/∂t. Importantly, it is found that the difference between the production and destruction terms, or net surface ‘stretch’, scales with the same factor as ∂Σ/∂t and the advection term, which may have a significant impact on how the evolution of Σ is understood and modelled in turbulent flows.

42 ENGINEERING↗

DOE OSTI2019

Assessing the Quality of Molecular Simulations for Vapor–Liquid Equilibria: An Analysis of the TraPPE Database

As molecular modeling and simulation techniques become increasingly important sources of thermophysical property and phase equilibrium data, the ability to assess the robustness of that data becomes more critical. Recently, the use of the compressibility factor (Z) has been suggested as a metric for testing the quality of simulation data for vapor–liquid equilibria (VLE). Int his paper, we analyze predicted VLE data from the transferable potentials for phase equilibria (TraPPE) database and show that, apart from data entry or typographical errors, Z will always be well-behaved in Gibbs ensemble Monte Carlo (GEMC) simulations even when the simulations are not sufficiently equilibrated. However, this is not true for grand canonical Monte Carlo simulations. When the pressure is calculated from the internal forces, then pressure and density are strongly correlated for the vapor phase and, for GEMC simulations, it is recommended to treat Z as an instantaneous mechanical property. From analysis of the TraPPE VLE data, we propose a complementary metric based on the predicted vapor pressures at three neighboring temperatures and their deviation from a local Clausius–Clapeyron fit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Solar-driven selective conversion of millimolar dissolved carbon to fuels with molecular flux generation

Abstract The direct utilization of dissolved inorganic carbon in seawater for CO 2 conversion promises chemical production on-demand and with zero carbon footprint. Photoelectrochemical (PEC) CO 2 reduction (CO 2 R) devices promise the sustainable conversion of dissolved carbon in seawater to carbon products using sunlight as the only energy input. However, the diffusion-dominant transport mechanism and the near-zero concentration of CO 2 (aq) (CO 2 dissolved in aqueous solution) in static seawater has made it extremely challenging to achieve high solar-to-fuel (STF) efficiency and high carbon-product selectivity. Here, where CO 2 (aq) as a reactant generated in situ by acidification of HCO 3 - flows continuously from BiVO 4 photoanodes to Si photocathodes, enabling a single-step conversion of dissolved carbon into products. Our PEC device significantly increases the CO selectivity from 3% to 21%, which approaches the 30% theoretical limit according to multi-physics modeling. Meanwhile, the Si/BiVO 4 PEC CO 2 R device achieved a STF efficiency of 0.71%. Such flow engineering achieves flow-dependent selectivity, rate, and stability in simulated seawater, thus promising practical solar fuel production at scale.

Liu, Bin↗

DOE OSTI2026

Burn mode signatures enabled by high dynamic range fusion reaction history

The evolution of the fusion burn of a compressed inertial confinement fusion (ICF) implosion gives information about the evolution of the temperature, mass, and volume of the hot spot. Currently, the fusion reaction history has been measured with about a decade of dynamic range, giving information about just the peak of the fusion burn. There are proposals for extending the dynamic range to 1000×, measuring the rising edge of the burn earlier in time. Using fusion hot spot theory and a set of ICF simulations, we identify and categorize different stages of fusion burn and what signatures could be measured. For ice layered implosions, the details, conditions, and dynamics of the fusion burn propagation could be observed. For double shell and OMEGA scale implosions, the evolution of the final compression and rebounding shock before ignition could be observed. Higher dynamic range expands investigation into the logarithmic derivative, α, as a signature of various mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Master equation emulation and coherence preservation with classical control of a superconducting qubit

Open quantum systems are a topic of intense theoretical research. The use of master equations to model a system's evolution subject to an interaction with an external environment is one of the most successful theoretical paradigms. General experimental tools to study different open system realizations have been limited, and so it is highly desirable to develop experimental tools which emulate diverse master equation dynamics and give a way to test open systems theories. In this paper we demonstrate a systematic method for engineering specific system-environment interactions and emulating master equations of a particular form using classical stochastic noise in a superconducting transmon qubit. We also demonstrate that non-Markovian noise can be used as a resource to extend the coherence of a quantum system and counteract the adversarial effects of Markovian environments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Logical and Physical Reversibility of Conservative Skyrmion Logic

Magnetic skyrmions are nanoscale whirls of magnetism that can be propagated with electrical currents. The repulsion between skyrmions inspires their use for reversible computing based on the elastic billiard ball collisions proposed for conservative logic in 1982. In this letter, we evaluate the logical and physical reversibility of this skyrmion logic paradigm, as well as the limitations that must be addressed before dissipation-free computation can be realized.

42 ENGINEERING↗