Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation optimization”

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 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

HPC for optimizing process parameters to control material evolution in seamless induction hardening of wind turbine main shaft bearings

Work proposed in this project focused on understanding the effect of martensitic transformation in the steel on the potential for cracking during seamless induction hardening (SIH) as a function of process conditions to allow the process to optimally scale up. Large-scale, three-dimensional phase-field simulations of martensitic transformation were performed using MEUMAPPS-SS (Microstructure Evolution Using Massively Parallel Phase-field Simulations – Solid State) code developed at Oak Ridge National Laboratory. The simulations were guided by location-specific thermal history generated by experimental measurements of time-temperature history generated at The Timken Company. The simulations were able to capture the morphological evolution of the martensite variants in an Fe-1.0C-1.5Cr steel based on the Nishiyama-Wasserman (NW) orientation relationship. The simulations were also able to quantify the stress-state at the interface between impinging martensite variants. The simulations indicated that the magnitude of the various stress and strain components were dependent on the sizes of the impinging plates with a reduction in these quantities with reduced plate size in agreement with experimental findings. The results obtained from the simulations will be used to guide the optimization of the alloy thermal conditions to eliminate quench cracking during SIH of bearing steels.

17 WIND ENERGY↗
NASA NTRS1974

Computer simulation of the NASA water vapor electrolysis reactor

The water vapor electrolysis (WVE) reactor is a spacecraft waste reclamation system for extended-mission manned spacecraft. The WVE reactor's raw material is water, its product oxygen. A computer simulation of the WVE operational processes provided the data required for an optimal design of the WVE unit. The simulation process was implemented with the aid of a FORTRAN IV routine.

Bloom, A. M.↗
DOE OSTI2026

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

Self-Tuning of Design Variables for Generalized Predictive Control

Three techniques are introduced to determine the order and control weighting for the design of a generalized predictive controller. These techniques are based on the application of fuzzy logic, genetic algorithms, and simulated annealing to conduct an optimal search on specific performance indexes or objective functions. Fuzzy logic is found to be feasible for real-time and on-line implementation due to its smooth and quick convergence. On the other hand, genetic algorithms and simulated annealing are applicable for initial estimation of the model order and control weighting, and final fine-tuning within a small region of the solution space, Several numerical simulations for a multiple-input and multiple-output system are given to illustrate the techniques developed in this paper.

Lin, Chaung↗
DOE OSTI2025

Optimizing Control for Efficient Load Shifting with Thermal Energy Storage in Existing HVAC Systems

This project developed the integration of a direct-contact heat exchanger (DCHX) based thermal energy storage (TES) system with a chiller–air handling unit (AHU) plant to evaluate its potential for reducing building energy costs. Detailed physical models of the TES unit, building envelope, and HVAC components were developed alongside simplified control-oriented models to support both high-fidelity simulation and real-time optimization. Two control strategies were implemented and compared: a rule-based control (RBC) aligned with utility time-of-use (ToU) rates, and a model predictive control (MPC) framework leveraging forecasts of building load, weather, and internal gains. Simulation results show that the RBC strategy reduced daily electricity costs by around 30% by shifting cooling production from on-peak to off-peak hours. In contrast, the MPC strategy achieved significantly greater performance, reducing daily operating costs by up to 44% and peak-hour costs by more than 60%. Both strategies maintained indoor thermal comfort within acceptable limits, with MPC further improving load distribution and reducing equipment cycling. The outcomes confirm that TES integration, particularly when coordinated with advanced predictive control, can provide substantial cost savings and on-peak demand reduction. . These findings directly support the U.S. Department of Energy’s goals for grid-interactive efficient buildings and demonstrate the potential of TES-enabled HVAC systems for scalable deployment across the commercial building.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1981

Solar Power Satellite antenna phase control system hardware simulation, phase 4. Volume 2: Analytical simulation of SPS system performance

The pilot signal parameter optimization and power transponder analyses are presented. The SPS antenna phase control system is modeled and the hardware simulation study described. Ionospheric and system phase error effects and the effects of high power amplifier phase and amplitude jitters are considered. Parameter optimization of the spread spectrum receiver, consisting of the carrier tracking loop and the code tracking loop, is described.

Lindsey, W. C.↗
DOE OSTI2017

Dual-Polarity Ion Trap Mass Spectrometry: Dynamic Monitoring and Controlling Gas-phase Ion–Ion Reactions

A dual-polarity linear ion trap (LIT) mass spectrometer was developed in this study, and the method for simultaneously controlling and detecting cations and anions was proposed and realized in the LIT. With the application of an additional dipolar DC field on the ejection electrodes of an LIT, dual-polarity mass spectra could be obtained, which include both the mass-to-charge (m/z) ratio and charge polarity information of an ion. Compared with conventional method, the ion ejection and detection efficiency could also be improved by about one-fold. Furthermore, ion–ion reactions within the LIT could be dynamically controlled and monitored by manipulating the distributions of ions with opposite charge polarities. This method was then used to control and study the reaction kinetics of ion–ion reactions, including electron transfer dissociation (ETD) and charge inversion reactions. A dual-polarity collision-induced dissociation (CID) experiment was proposed and performed to enhance the sequence coverage of a peptide ion. Ion trajectory simulations were also carried out for concept validation and system optimization. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1990

Applications of artificial neural nets in structural mechanics

A brief introduction to the fundamental of Neural Nets is given, followed by two applications in structural optimization. In the first case, the feasibility of simulating with neural nets the many structural analyses performed during optimization iterations was studied. In the second case, the concept of using neural nets to capture design expertise was studied.

Berke, Laszlo↗
NASA NTRS1992

Applications of artificial neural nets in structural mechanics

A brief introduction to the fundamental of Neural Nets is given, followed by two applications in structural optimization. In the first case, the feasibility of simulating with neural nets the many structural analyses performed during optimization iterations was studied. In the second case, the concept of using neural nets to capture design expertise was studied.

Berke, L.↗
NASA NTRS2014

Frequency-Modulated, Continuous-Wave Laser Ranging Using Photon-Counting Detectors

Optical ranging is a problem of estimating the round-trip flight time of a phase- or amplitude-modulated optical beam that reflects off of a target. Frequency- modulated, continuous-wave (FMCW) ranging systems obtain this estimate by performing an interferometric measurement between a local frequency- modulated laser beam and a delayed copy returning from the target. The range estimate is formed by mixing the target-return field with the local reference field on a beamsplitter and detecting the resultant beat modulation. In conventional FMCW ranging, the source modulation is linear in instantaneous frequency, the reference-arm field has many more photons than the target-return field, and the time-of-flight estimate is generated by balanced difference- detection of the beamsplitter output, followed by a frequency-domain peak search. This work focused on determining the maximum-likelihood (ML) estimation algorithm when continuous-time photoncounting detectors are used. It is founded on a rigorous statistical characterization of the (random) photoelectron emission times as a function of the incident optical field, including the deleterious effects caused by dark current and dead time. These statistics enable derivation of the Cramér-Rao lower bound (CRB) on the accuracy of FMCW ranging, and derivation of the ML estimator, whose performance approaches this bound at high photon flux. The estimation algorithm was developed, and its optimality properties were shown in simulation. Experimental data show that it performs better than the conventional estimation algorithms used. The demonstrated improvement is a factor of 1.414 over frequency-domainbased estimation. If the target interrogating photons and the local reference field photons are costed equally, the optimal allocation of photons between these two arms is to have them equally distributed. This is different than the state of the art, in which the local field is stronger than the target return. The optimal processing of the photocurrent processes at the outputs of the two detectors is to perform log-matched filtering followed by a summation and peak detection. This implies that neither difference detection, nor Fourier-domain peak detection, which are the staples of the state-of-the-art systems, is optimal when a weak local oscillator is employed.

Erkmen, Baris I.↗
DOE OSTI2018

Exploring the posterior surface of the large scale structure reconstruction

The large scale structure (LSS) of the universe is generated by the linear density gaussian modes, which are evolved into the observed nonlinear LSS. Given a set of data the posterior surface of the modes is convex in the linear regime, leading to a unique global maximum (MAP), but this is no longer guaranteed in the nonlinear regime. In this paper we investigate the nature of posterior surface using the recently developed MAP reconstruction method, with a simplified but realistic N-body simulation as the forward model. The reconstruction method uses optimization with analytic gradients from back-propagation through the simulation. For low noise cases we recover the initial conditions well into the nonlinear regime (k~ 1 h/Mpc) nearly perfectly. We show that the large scale modes can be recovered more precisely than the linear expectation, which we argue is a consequence of nonlinear mode coupling. For noise levels achievable with current and planned LSS surveys the reconstruction cannot recover very small scales due to noise. We see some evidence of non-convexity, specially for smaller scales where the non-injective nature of the mappings: several very different initial conditions leading to the same near perfect final data reconstruction. We investigate the nature of these phenomena further using a 1-d toy gravity model, where many well separated local maximas are found to have identical data likelihood but differ in the prior. We also show that in 1-d the prior favors some solutions over the true solution, though no clear evidence of these in 3-d. Our main conclusion is that on very small scales and for a very low noise the posterior surface is multi-modal and the global maximum may be unreachable with standard methods, while for realistic noise levels in the context of the current and next generation LSS surveys MAP optimization method is likely to be nearly optimal.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2007

A Case Study Using Modeling and Simulation to Predict Logistics Supply Chain Issues

Optimization of critical supply chains to deliver thousands of parts, materials, sub-assemblies, and vehicle structures as needed is vital to the success of the Constellation Program. Thorough analysis needs to be performed on the integrated supply chain processes to plan, source, make, deliver, and return critical items efficiently. Process modeling provides simulation technology-based, predictive solutions for supply chain problems which enable decision makers to reduce costs, accelerate cycle time and improve business performance. For example, United Space Alliance, LLC utilized this approach in late 2006 to build simulation models that recreated shuttle orbiter thruster failures and predicted the potential impact of thruster removals on logistics spare assets. The main objective was the early identification of possible problems in providing thruster spares for the remainder of the Shuttle Flight Manifest. After extensive analysis the model results were used to quantify potential problems and led to improvement actions in the supply chain. Similarly the proper modeling and analysis of Constellation parts, materials, operations, and information flows will help ensure the efficiency of the critical logistics supply chains and the overall success of the program.

Tucker, David A.↗
DOE OSTI2025

Design Overview of a High-Pressure Helium Flow Visualization Apparatus for Blanket Cooling Studies

Cooling of the fusion blanket first wall remains a significant challenge given the adverse conditions of heat and particle flux encountered near the plasma. Helium emerges as an attractive cooling candidate because of its chemical and neutronic inertness and separability from hydrogenic species (e.g. tritium). Because of the low thermal mass of helium, optimization of these coolant channels is warranted to provide high heat transfer performance at low pumping costs. Increasingly, computational fluid dynamics (CFD) simulations are employed to model and optimize these flow channels, and accompanying experimental data are needed to validate the predictions of these models. To provide the aforementioned experimental data, a high-pressure helium flow visualization upgrade has been designed for the Helium Flow Loop Experiment facility. This apparatus was built to American Society of Mechanical Engineers boiler and pressure vessel standards to withstand operating pressure of 4 MPa and mated to high-pressure glass windows. Seedless flow visualization is performed via high-speed background oriented schlieren (BOS), with image correlation used for time-resolved two-dimensional velocimetry at frequencies in excess of 60 kHz. Rectangular flow channel test articles are additively manufactured via laser powder bed fusion and installed into this visualization apparatus, with one-sided heating supplied by resistive heaters. In conclusion, the chosen test geometries were informed by prior CFD simulations, and the helium flow structures observed via BOS (detachment, recirculation, etc.) will be used for the validation of these accompanying models, in support of the design and optimization of blanket cooling channel configurations.

Helium flow↗
DOE OSTI2017

Design and optimization of 3D printed air-cooled heat sinks based on genetic algorithms

Enhancing power density and reliability of powerelectronics is extremely important in power electronicsapplications. One of the key challenges in the design process isto design the optimum heat sink. In this paper, an algorithm isproposed to design air-cooled heat sinks using genetic algorithm(GA) and finite element analysis (FEA) simulations. While theGA generates a population of candidate heat sinks in eachiteration, FEA simulations are used to evaluate the fitnessfunction of each. The fitness function considered in this paper isthe maximum junction temperature of the semiconductordevices. With an approach that prefers “survival of the fittest”,a heat sink providing better performance than the conventionalheat sinks is obtained. The simulation and experimentalevaluations of the optimized air-cooled heat sink are alsoincluded in the paper.

Wu, Tong↗
DOE OSTI2024

Bayesian blacksmithing: discovering thermomechanical properties and deformation mechanisms in high-entropy refractory alloys

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Bayesian optimal experimental design for constitutive model calibration

Computational simulation is increasingly relied upon for high/consequence engineering decisions, which necessitates a high confidence in the calibration of and predictions from complex material models. However, the calibration and validation of material models is often a discrete, multi-stage process that is decoupled from material characterization activities, which means the data collected does not always align with the data that is needed. To address this issue, an integrated workflow for delivering an enhanced characterization and calibration procedure—Interlaced Characterization and Calibration (ICC)—is introduced and demonstrated. Further, this framework leverages Bayesian optimal experimental design (BOED), which creates a line of communication between model calibration needs and data collection capabilities in order to optimize the information content gathered from the experiments for model calibration. Eventually, the ICC framework will be used in quasi real-time to actively control experiments of complex specimens for the calibration of a high-fidelity material model. This work presents the critical first piece of algorithm development and a demonstration in determining the optimal load path of a cruciform specimen with simulated data. Calibration results, obtained via Bayesian inference, from the integrated ICC approach are compared to calibrations performed by choosing the load path a priori based on human intuition, as is traditionally done. The calibration results are communicated through parameter uncertainties which are propagated to the model output space (i.e. stress–strain). In these exemplar problems, data generated within the ICC framework resulted in calibrated model parameters with reduced measures of uncertainty compared to the traditional approaches.

42 ENGINEERING↗
DOE OSTI2020

Calorimetry for low-energy electrons using charge and light in liquid argon

Precise calorimetric reconstruction of 5–50 MeV electrons in liquid argon time projection chambers (LArTPCs) will enable the study of astrophysical neutrinos in DUNE and could enhance the physics reach of oscillation analyses. Liquid argon scintillation light has the potential to improve energy reconstruction for low-energy electrons over charge-based measurements alone. Here we demonstrate light-augmented calorimetry for low-energy electrons in a single-phase LArTPC using a sample of Michel electrons from decays of stopping cosmic muons in the LArIAT experiment at Fermilab. Michel electron energy spectra are reconstructed using both a traditional charge-based approach as well as a more holistic approach that incorporates both charge and light. A maximum-likelihood fitter, using LArIAT’s well-tuned simulation, is developed for combining these quantities to achieve optimal energy resolution. A sample of isolated electrons is simulated to better determine the energy resolution expected for astrophysical electron-neutrino charged-current interaction final states. In LArIAT, which has very low wire noise and an average light yield of 18 pe/MeV, an energy resolution of σ/E≃9.3%/E⊕1.3% is achieved. Samples are then generated with varying wire noise levels and light yields to gauge the impact of light-augmented calorimetry in larger LArTPCs. At a charge-readout signal-to-noise of S/N≃30, for example, the energy resolution for electrons below 40 MeV is improved by ≈10%, ≈20%, and ≈40% over charge-only calorimetry for average light yields of 10 pe/MeV, 20 pe/MeV, and 100 pe/MeV, respectively.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Medium- and Heavy-Duty Truck Duty Cycles

This dataset provides second-by-second duty cycle data for Class 6 and Class 8 diesel trucks in Texas, including key vehicle metrics, engine-related data, and GPS data (excluding GPS latitude and longitude to ensure confidentiality). The data were collected via tablets installed on the trucks and organized into daily datasets, each associated with a unique vehicle ID and date. There are 12 daily datasets for Class 6 diesel trucks (three unique vehicle IDs) and 43 daily datasets for Class 8 diesel trucks (six unique vehicle IDs). The units associated with each column are included in the name. The engine performance data include columns such as engine speed, engine percent torque, and engine fuel rate. Road grade (%/100) was estimated using the GPS altitude and wheel-based vehicle speed, which is used as an input for FASTSim. Cumulative distance was also calculated using the wheel-based vehicle speed. Additional columns include: - Engine Speed (RPM): Removed inaccurate readings and used to calculate angular velocity (radians/second). - Torque (N·m): Calculated using engine percent torque, nominal friction percent torque, and engine reference torque values (those columns were removed from dataset), then normalized to express as torque (%). - Flywheel Power (%): Calculated using the angular velocity and torque (in kW), then normalized as a percentage of the maximum value. - Engine Fuel Rate (%) and Torque (%): Both metrics were normalized by dividing by their respective maximum values within each dataset to express them as percentages. The datasets were analyzed to assess the energy impact of various driving behaviors, simulate energy efficiency, and recommend optimal routes for diesel trucks using NLR’s tool called RouteE. For driver coaching, factors like speed and acceleration limits were considered, and idle periods were reduced (assuming the engine was off during idling) to adjust each drive cycle. These adjusted drive cycles were then simulated in FASTSim to evaluate their effect on fleet energy consumption and estimate potential energy savings. The original cycles are available for download on this page ![image](CoVaR_Image_for_Data_Page_Kenworth_Truck.jpg)

1Hz↗