Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational modeling and 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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Ambient and Initial Temperature Effects on Energy Consumption Rate Modeled in FASTSim

Ambient and initial temperatures significantly impact the energy consumption rate (ECR) of battery electric vehicles (BEVs) due to auxiliary loads and the temperature dependence of battery efficiency. This study introduces a streamlined, physics-based thermal modeling approach within the FASTSim tool that bridges the gap between oversimplified constant-load models and computationally expensive high-fidelity simulations. By employing a lumped thermal mass framework, the model captures fundamental energy balances and critical non-linear energy penalties while maintaining the computational efficiency required for expansive sensitivity studies. The simulations evaluated a compact BEV hatchback with a resistive heater over city (UDDS) and highway (HWFET) test cycles. Compared to a 22 degrees Celsius initial and ambient temperature baseline, a -7 degrees Celsius initial/ambient temperature resulted in a 221% increase in the ECR for the city cycle and a 100% increase for the highway cycle. Conversely, a 45 degrees Celsius initial / 40 degrees Celsius ambient temperature resulted in a 40% increase for UDDS and an 18% increase for HWFET. These results demonstrate that while cold conditions impose the most severe energy penalties due to resistive heating, the impact is consistently more pronounced in city driving where auxiliary loads represent a larger proportion of total energy. This lightweight yet robust framework enables researchers to rapidly quantify BEV thermal sensitivity across diverse climates without the need for high-overhead simulation environments.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Ion trap long-range XY model for quantum state transfer and optimal spatial search

Abstract Linear ion trap chains are a promising platform for quantum computation and simulation. The XY model with long-range interactions can be implemented with a single side-band Mølmer–Sørensen scheme, giving interactions that decay as 1 / r α , where α parameterises the interaction range. Lower α leads to longer range interactions, allowing faster long-range gate operations for quantum computing. However, decreasing α causes an increased generation of coherent phonons and appears to dephase the effective XY interaction model. We characterise and show how to correct for this effect completely, allowing lower α interactions to be coherently implemented. Ion trap chains are thus shown to be a viable platform for spatial quantum search in optimal O ( N ) time, for N ions. Finally, we introduce a O ( N ) quantum state transfer protocol, with a qubit encoding that maintains a high fidelity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1996

Simulations of light effects on the human circadian pacemaker: implications for assessment of intrinsic period

The sensitivity of the human circadian system to light has been the subject of considerable debate. Using computer simulations of a recent quantitative model for the effects of light on the human circadian system, we investigated these effects of light during different experimental protocols. The results of the simulations indicate that the nonuniform distribution over the circadian cycle of exposure to ordinary room light seen in classical free-run studies, in which subjects select their exposure to light and darkness, can result in an observed period of approximately 25 h, even when the intrinsic period of the subject's endogenous circadian pacemaker is much closer to 24 h. Other simulation results suggest that accurate assessment of the true intrinsic period of the human circadian pacemaker requires low ambient light intensities (approximately 10-15 lx) during scheduled wake episodes, desynchrony of the imposed light-dark cycle from the endogenous circadian oscillator, and a study length of at least 20 days. Although these simulations await further experimental substantiation, they highlight the sensitivity to light of the human circadian system and the potential confounding influence of light on the assessment of the intrinsic period of the circadian pacemaker.

Non-NASA Center↗
NASA NTRSDate not supplied

Scalable Variable Charge Molecular Dynamics Simulations of Metal-Oxide Systems

Interfaces between metals and oxides are an important feature in many technologically relevant materials, e.g., oxidation of metal surfaces, oxide-dispersion strengthened (ODS) alloys, dielectric components, and thermal barrier coatings among others. Experimental studies of such interfaces are challenging since the majority are buried within the bulk, making computational modeling an attractive alternative. Molecular dynamics (MD) simulations operate at the length scales relevant to many interface-mediated mechanisms, but the requisite interatomic potentials for metal-oxide systems require computationally expensive variable charge schemes to account for the disparate bonding types, thus often limiting their effectiveness. Here we introduce several improvements to the charge transfer interatomic potential (CTIP) model which enable greater computational efficiency for large scale MD simulations. Then, using a new CTIP parametrization for the Ni-O system, we demonstrate its capabilities to capture critical atomic scale mechanisms associated with metal-oxide interfaces. Long time scale simulations (>10 ns) are used to investigate high temperature oxidation and oxide precipitation from the melt, and large length scale simulations (> 1 million atoms) are used to study the interaction of dislocations with oxide particles. We have implemented the new CTIP model in the widely used, open-source MD code LAMMPS.

Gabriel Plummer↗
DOE OSTI2021

Simulation and Modeling of Time-Resolved X-Ray Detector for the Saturn Accelerator

In this work, we present the technology-aided computer design (TCAD) device simulation and modeling of a silicon p-i-n diode for detecting time-dependent X-ray radiation. We show that the simulated forward and reverse breakdown current–voltage characteristics agree well with the measured data under nonradiation environment by only calibrating carrier lifetimes for the forward bias case and avalanche model critical fields for the reverse bias condition. Using the calibrated parameters and other nominal material properties, we simulated the radiation responses of the p-i-n diode and compared with experimental data when the diode was exposed to X-ray radiation at Sandia’s Saturn facility and the Idaho State University (ISU) TriMeV facility. For Saturn’s Gaussian dose-rate pulses, we show three findings from TCAD simulations. First, the simulated photocurrents are in excellent agreement with the measured data for two dose-rate pulses with peak values of 1.16×10 10 and 1.88×10 10 rad(Si)/s. Second, the simulation results of high dose-rate pulses predict increased delayed photocurrents with longer time tails in the diode electrical responses due to excess carrier generation. Third, simulated peak values of diode radiation responses versus peak dose rates at different bias conditions provide useful guidance to determine the dose-rate range that the p-i-n diode can reliably detect in experiment. For TriMeV’s non-Gaussian dose-rate pulse, our simulated diode response is in decent agreement with the measured data without further calibration. We also studied the effects of device geometry, recombination process, and dose-rate enhancement via TCAD simulations to understand the higher measured response in the time after the peak dose-rate radiation for the p-i-n diode exposed to TriMeV irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1977

Effect of altered gravity on temperature regulation in mammals: Investigation of gravity effect on temperature regulation in mammals

Male, Long-Evans hooded rats were instrumented for monitoring core and hypothalamic temperatures as well as shivering and nonshivering thermogenesis in response to decreased ambient temperature in order to characterize the nature of the neural controller of temperature in rats at 1G and evaluate chronic implantation techniques for the monitoring of appropriate parameters at hypergravic fields. The thermoregulatory responses of cold-exposed rats at 2G were compared to those at 1G. A computer model was developed to simulate the thermoregulatory system in the rat. Observations at 1 and 2G were extended to acceleration fields of 1.5, 3.0 and 4.0G and the computer model was modified for application to altered gravity conditions. Changes in the acceleration field resulted in inadequate heat generation rather than increased heat loss. Acceleration appears to impair the ability of the neurocontroller to appropriately integrate input signals for body temperature maintenance.

Horwitz, B. A.↗
DOE OSTI2022

Exploring the Use of Novel Spatial Accelerators in Scientific Applications

Driven by the need to find alternative accelerators which can viably replace GPUs in next-generation Supercomputing systems, this paper proposes a methodology to enable agile application/hardware co-design. The application-first methodology provides the ability to come up with design of accelerators while working with real-world workloads, available accelerators, and system software. The iterative design process targets a set of kernels in a workload for performance estimates that can prune the design space for later phases of detailed architectural evaluations. To this effect, in this paper, a novel data-parallel device model is introduced that simulates the latency of performance-sensitive operations in an accelerator including data transfers and kernel computation using multi-core CPUs. The use of off-the-shelf simulators, such as pre-RTL simulator Aladdin or multiple tools available for exploring the design of deep neural network accelerators (e.g., Timeloop) is demonstrated for evaluation of various accelerator designs using applications with realistic inputs. Examples of multiple device configurations that are instantiable in a system are explored to evaluate the performance benefit of deploying novel accelerators. The proposed device is integrated with a programming model and system software to potentially explore the impacts of high-level programming languages/compilers and low-level effects such as task scheduling on multiple accelerators. We analyze our methodology for a set of applications that represent high-performance computing (HPC) and graph analytics. The applications include a computational chemistry kernel realized using tensor contractions, triangle counting, GraphSAGE and Breadth-first Search. These applications include kernels such as dense matrix-dense matrix multiplication, sparse matrix-spare matrix multiplication, and sparse matrix-dense vector multiplication. Our results indicate potential performance benefits and insights for system design by including accelerators that realize these kernels along-side general purpose accelerators.

AI, codesign, Accelerated Computing, Modeling and ↗
NASA NTRS1975

Performance analysis of wideband data and television channels

Several aspects are discussed of space shuttle communications, including the return link (shuttle-to-ground) relayed through a satellite repeater (TDRS). The repeater exhibits nonlinear amplification and an amplitude-dependent phase shift. Models were developed for various link configurations, and computer simulation programs based on these models are described. Certain analytical results on system performance were also obtained. For the system parameters assumed, the results indicate approximately 1 db degradation relative to a link employing a linear repeater. While this degradation is dependent upon the repeater, filter bandwidths, and modulation parameters used, the programs can accommodate changes to any of these quantities. Thus the programs can be applied to determine the performance with any given set of parameters, or used as an aid in link design.

Geist, J. M.↗
NASA NTRS2012

A Summary of Data and Findings from the First Aeroelastic Prediction Workshop

This paper summarizes data and findings from the first Aeroelastic Prediction Workshop (AePW) held in April, 2012. The workshop has been designed as a series of technical interchange meetings to assess the state of the art of computational methods for predicting unsteady flowfields and static and dynamic aeroelastic response. The goals are to provide an impartial forum to evaluate the effectiveness of existing computer codes and modeling techniques to simulate aeroelastic problems, and to identify computational and experimental areas needing additional research and development. For this initial workshop, three subject configurations have been chosen from existing wind tunnel data sets where there is pertinent experimental data available for comparison. Participant researchers analyzed one or more of the subject configurations and results from all of these computations were compared at the workshop. Keywords: Unsteady Aerodynamics, Aeroelasticity, Computational Fluid Dynamics, Transonic Flow, Separated Flow.

Schuster, David M.↗
NASA NTRS2013

Analysis of Test Case Computations and Experiments for the First Aeroelastic Prediction Workshop

This paper compares computational and experimental data from the Aeroelastic Prediction Workshop (AePW) held in April 2012. This workshop was designed as a series of technical interchange meetings to assess the state of the art of computational methods for predicting unsteady flowfields and static and dynamic aeroelastic response. The goals are to provide an impartial forum to evaluate the effectiveness of existing computer codes and modeling techniques to simulate aeroelastic problems and to identify computational and experimental areas needing additional research and development. Three subject configurations were chosen from existing wind-tunnel data sets where there is pertinent experimental data available for comparison. Participant researchers analyzed one or more of the subject configurations, and results from all of these computations were compared at the workshop.

Schuster, David M.↗
DOE OSTI2020

A transported Livengood-Wu integral model for knock prediction in CFD simulation

This work describes the development of a transported Livengood-Wu (L-W) integral model for computational fluid dynamics (CFD) simulation to predict auto-ignition and engine knock tendency. The currently employed L-W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition combustion for three different fuels, showing good accuracy in prediction of auto-ignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L-W integral model is coupled with G-equation model to indicate end-gas auto-ignition and knock tendency in CFD simulations of a direct injection spark ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely Co-Optima Alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multi-variable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency indicated by the fuel energy contained in the auto-ignited region of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.

Yue, Zongyu↗
NASA NTRS1995

Spectral variation of high power microwave pulse propagating in a self-generated plasma

A systematic study to understand the spectral variation of a high power microwave pulse propagating in a self-generated plasma is carried out. It includes the theoretical formulation, experimental demonstration, and computer simulations and computer experiments. The experiment of pulse propagation is conducted in a vacuum chamber filled with dry air (approximately 0.2 torr); the chamber is made of a 2 ft. cube of Plexiglas. A rectangular microwave pulse (1 microsec pulse width and 3.27 GHz carrier frequency) is fed into the cube through an S band microwave horn placed at one side of the chamber. A second S-band horn placed at the opposite side of the chamber is used to receive the transmitted pulse. The spectra of the incident pulse and transmitted pulse are then compared. As the power of the incident pulse is only slightly (less than 15%) above the breakdown threshold power of the background air, the peak of the spectrum of the transmitted pulse is upshifted from the carrier frequency 3.27 GHz of the incident pulse. However, as the power of the incident pulse exceeds the breakdown threshold power of the background air by 30%, a different phenomenon appears. The spectrum of the transmitted pulse begins to have two peaks. One is upshifted and the other one downshifted from the single peak location of the incident pulse. The amount of frequency downshift is comparable to that of the upshifted frequency. A theoretical model describing the experiment of pulse propagation in a self-generated plasma is developed. There are excellent agreements between the experimental results and computer simulations based on this theoretical model, which is also used to further carry out computer experiments identifying the role of plasma introduced wave loss on the result of frequency downshift phenomenon.

Ren, A.↗
DOE OSTI2022

In Situ X-ray Radiography and Computational Modeling to Predict Grain Morphology in $\beta$-Titanium during Simulated Additive Manufacturing

The continued development of metal additive manufacturing (AM) has expanded the engineering metallic alloys for which these processes may be applied, including beta-titanium alloys with desirable strength-to-density ratios. To understand the response of beta-titanium alloys to AM processing, solidification and microstructure evolution needs to be investigated. In particular, thermal gradients (Gs) and solidification velocities (Vs) experienced during AM are needed to link processing to microstructure development, including the columnar-to-equiaxed transition (CET). In this work, in situ synchrotron X-ray radiography of the beta-titanium alloy Ti-10V-2Fe-3Al (wt.%) (Ti-1023) during simulated laser-powder bed fusion (L-PBF) was performed at the Advanced Photon Source at Argonne National Laboratory, allowing for direct determination of Vs. Two different computational modeling tools, SYSWELD and FLOW-3D, were utilized to investigate the solidification conditions of spot and raster melt scenarios. The predicted Vs obtained from both pieces of computational software exhibited good agreement with those obtained from in situ synchrotron X-ray radiography measurements. The model that accounted for fluid flow also showed the ability to predict trends unobservable in the in situ synchrotron X-ray radiography, but are known to occur during rapid solidification. A CET model for Ti-1023 was also developed using the Kurz–Giovanola–Trivedi model, which allowed modeled Gs and Vs to be compared in the context of predicted grain morphologies. Both pieces of software were in agreement for morphology predictions of spot-melts, but drastically differed for raster predictions. The discrepancy is attributable to the difference in accounting for fluid flow, resulting in magnitude-different values of Gs for similar Vs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Learning compact physics‐aware delayed photocurrent models using dynamic mode decomposition

Abstract Radiation‐induced photocurrent in semiconductor devices can be simulated using complex physics‐based models, which are accurate, but computationally expensive. This presents a challenge for implementing device characteristics in high‐level circuit simulations where it is computationally infeasible to evaluate detailed models for multiple individual circuit elements. In this work we demonstrate a procedure for learning compact delayed photocurrent models that are efficient enough to implement in large‐scale circuit simulations, but remain faithful to the underlying physics. Our approach utilizes dynamic mode decomposition (DMD), a system identification technique for learning reduced‐order discrete‐time dynamical systems from time series data based on singular value decomposition. To obtain physics‐aware device models, we simulate the excess carrier density induced by radiation pulses by solving numerically the ambipolar diffusion equation, then use the simulated internal state as training data for the DMD algorithm. Our results show that the significantly reduced‐order delayed photocurrent models obtained via this method accurately approximate the dynamics of the internal excess carrier density—which can be used to calculate the induced current at the device boundaries—while remaining compact enough to incorporate into larger circuit simulations.

Hanson, Joshua↗
DOE OSTI2026

Trigonometric continuous-variable gates and hybrid quantum simulations of the sine-Gordon model

Hybrid qubit-qumode quantum computing platforms provide a natural setting for simulating interacting bosonic quantum field theories. However, existing continuous-variable gate constructions rely predominantly on polynomial functions of canonical quadratures. In this work, we introduce a complementary universality paradigm based on trigonometric continuous-variable gates, which enable a Fourier-like representation of bosonic operators and are particularly well suited for periodic and non-perturbative interactions. We present an ancilla-based framework for implementing trigonometric gates with arguments given by arbitrary Hermitian functions of qumode quadratures. The protocol yields unitary gates deterministically, and non-unitary gates through probabilistic post-selection. As a concrete application, we develop a hybrid qubit-qumode quantum simulation of the lattice sine-Gordon model. Using these gates, we prepare ground states via quantum imaginary-time evolution, simulate real-time dynamics, compute time-dependent vertex two-point correlation functions, and extract quantum kink profiles under topological boundary conditions. Our results demonstrate that trigonometric continuous-variable gates provide a physically natural framework for simulating interacting field theories on near-term hybrid quantum hardware, while establishing a parallel route to universality beyond polynomial gate constructions. We expect that the trigonometric gates introduced here to find broader applications, including quantum simulations of condensed matter systems, quantum chemistry, and biological models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRSDate not supplied

Unsteady CFD Simulations of a Compression Corner Geometry Using Wall-Modeled LES Methods in Loci/CHEM

As computing technology continues to improve, simulating unsteady aerodynamic environments on launch vehicle geometries becomes more viable. For NASA’s Space Launch System, the flowfield is characterized by large regions of unsteady turbulent wall-bounded flow. Simulating these unsteady wall-bounded environments is of interest to the SLS Aerodynamics team. Wall-modeled LES methods are currently the best approach for simulating these environments in engineering applications. Less-expensive computationally than fully-resolved methods while maintaining many of the benefits in simulating unsteady environments. Before applying these methods blindly to large geometries, best practices can first be applied and understood on smaller canonical problems.

Bryson Frank↗
DOE OSTI2025

Steps toward quantum simulations of hadronization and energy loss in dense matter

A framework for simulating the real-time dynamics of composite particles in a simple model of dense matter that is amenable to quantum computers is developed. As a demonstration, we perform classical simulations of heavy-hadrons propagating through a dense medium in the Schwinger model. Measurements of the time-dependent energy and charge density are used to identify mechanisms responsible for energy loss and hadron production (hadronization). A study of entanglement dynamics highlights the importance of quantum coherence between the particles that make up the dense medium. Throughout this work, care is taken to isolate, and remove, phenomena that arise solely from a finite lattice spacing. It is found that signatures of entanglement are more sensitive to lattice artifacts than other observables. Toward quantum simulations, we present an efficient method and the corresponding quantum circuits for preparing ground states in the presence of heavy mesons. Finally, these circuits are used to estimate the resources required to simulate in-medium energy loss and hadronization in the Schwinger model using quantum computers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

Evolutionary Computational Methods for Identifying Emergent Behavior in Autonomous Systems

A technique based on Evolutionary Computational Methods (ECMs) was developed that allows for the automated optimization of complex computationally modeled systems, such as autonomous systems. The primary technology, which enables the ECM to find optimal solutions in complex search spaces, derives from evolutionary algorithms such as the genetic algorithm and differential evolution. These methods are based on biological processes, particularly genetics, and define an iterative process that evolves parameter sets into an optimum. Evolutionary computation is a method that operates on a population of existing computational-based engineering models (or simulators) and competes them using biologically inspired genetic operators on large parallel cluster computers. The result is the ability to automatically find design optimizations and trades, and thereby greatly amplify the role of the system engineer.

Terrile, Richard J.↗