SEARCH · Engineering Papers

Results for “simulation analytics”

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 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Analysis of Automated Transit Network Systems with Battery-Electric Vehicles in Automated Mobility Districts: Preprint

The paper provides an overview of the insights and findings of the National Renewable Energy Laboratory's (NREL) ongoing research on the implementation of automated mobility districts (AMDs). AMD is a term coined by NREL to describe a geographically defined district or major activity center located in a dense urban setting with mobility applications provided by automated/autonomous vehicle (AV) systems spanning internal circulation and first-mile/last-mile connections to regional transportation hubs. Research over the past five years has focused on understanding the evolution of AMDs, beginning with demonstrations of automated shuttles prior to the pandemic, to more integrated on-demand mobility systems currently in initial stages of deployment. Initial insights and findings from earlier studies include the need for designation of a "jurisdiction having authority," a clear vision of a complete system to provide end-to-end mobility services, and the requisite intelligent infrastructure to complement AV technology. NREL's most recent Phase III research investigates station boarding/alighting (curb) issues, the full electrification of fleets, and the need for a systems engineering methodology (SEM) to properly analyze the complexities resulting from the convergence of automation, on-demand mobility, and electrification of the transit systems within the AMD. The paper reviews the findings of AMD research conducted during Phase I, II, and III, with special emphasis on Phase III results with respect to a descriptive example of the proposed SEM when a "digital twin" analytical model is used to simulate the transport fleet's battery-electric vehicle miles of travel and associated duty cycles through a rigorous analytical assessment with a comprehensive modeling process.

ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2026

Monochromation of pulsed electron beams with terahertz radiation at a planar mirror

Exquisite control of electron beam energy is required for many electron spectroscopy and imaging applications. For both continuous and pulsed beams, the beam energy spread is fundamentally limited by the electron source, and is typically a sizable fraction of an electron-volt. In this paper, we present a means to reduce electron beam energy spread after emission to the level of a few 10s of meV rms using femtosecond photoemission and an interaction with laser-derived single- to few-cycle terahertz (THz) radiation. We show analytically and in particle tracking simulations that this interaction can remove energy spread stored in both the transverse and longitudinal degrees of freedom. We analytically formulate the limit of energy spread that this technique can achieve, and map the non-ideal affects arising at high frequencies. The interaction is mediated by the beam's passage through a mirror which is reflective to terahertz radiation but allows transmission of the majority of the electron beam (e.g. a wire mesh). This method then only requires beam current losses of a few tens of percent, far smaller than what is achieved in prism and slit-based electron monochromators.

Accelerator Physics (physics.acc-ph)↗

NASA NTRS1992

Comparison of simulated and actual wind shear radar data products

Prior to the development of the NASA experimental wind shear radar system, extensive computer simulations were conducted to determine the performance of the radar in combined weather and ground clutter environments. The simulation of the radar used analytical microburst models to determine weather returns and synthetic aperture radar (SAR) maps to determine ground clutter returns. These simulations were used to guide the development of hazard detection algorithms and to predict their performance. The structure of the radar simulation is reviewed. Actual flight data results from the Orlando and Denver tests are compared with simulated results. Areas of agreement and disagreement of actual and simulated results are shown.

Britt, Charles L.↗

NASA NTRS2010

Multi-Evaporator Miniature Loop Heat Pipe for Small Spacecraft Thermal Control: New Technologies and Validation Approach - Part 1

Under NASA s New Millennium Program Space Technology 8 (ST 8) Project, four experiments Thermal Loop, Dependable Microprocessor, SAILMAST, and UltraFlex - were conducted to advance the maturity of individual technologies from proof of concept to prototype demonstration in a relevant environment , i.e. from a technology readiness level (TRL) of 3 to a level of 6. This paper presents the new technologies and validation approach of the Thermal Loop experiment. The Thermal Loop is an advanced thermal control system consisting of a miniature loop heat pipe (MLHP) with multiple evaporators and multiple condensers designed for future small system applications requiring low mass, low power, and compactness. The MLHP retains all features of state-of-the-art loop heat pipes (LHPs) and offers additional advantages to enhance the functionality, performance, versatility, and reliability of the system. Details of the thermal loop concept, technical advances, benefits, objectives, level 1 requirements, and performance characteristics are described. Also included in the paper are descriptions of the test articles and mathematical modeling used for the technology validation. An MLHP breadboard was built and tested in the laboratory and thermal vacuum environments for TRL 4 and TRL 5 validations, and an MLHP proto-flight unit was built and tested in a thermal vacuum chamber for the TRL 6 validation. In addition, an analytical model was developed to simulate the steady state and transient behaviors of the MLHP during various validation tests. Capabilities and limitations of the analytical model are also addressed.

Ku, Jentung↗

NASA NTRS1996

Electromagnetic Compatibility Testing Studies

This report discusses the results on analytical models and measurement and simulation of statistical properties from a study of microwave reverberation (mode-stirred) chambers performed at Texas Tech University. Two analytical models of power transfer vs. frequency in a chamber, one for antenna-to-antenna transfer and the other for antenna to D-dot sensor, were experimentally validated in our chamber. Two examples are presented of the measurement and calculation of chamber Q, one for each of the models. Measurements of EM power density validate a theoretical probability distribution on and away from the chamber walls and also yield a distribution with larger standard deviation at frequencies below the range of validity of the theory. Measurements of EM power density at pairs of points which validate a theoretical spatial correlation function on the chamber walls and also yield a correlation function with larger correlation length, R(sub corr), at frequencies below the range of validity of the theory. A numerical simulation, employing a rectangular cavity with a moving wall shows agreement with the measurements. The determination that the lowest frequency at which the theoretical spatial correlation function is valid in our chamber is considerably higher than the lowest frequency recommended by current guidelines for utilizing reverberation chambers in EMC testing. Two suggestions have been made for future studies related to EMC testing.

Trost, Thomas F.↗

NASA NTRS1981

A computer program for the simulation of heat and moisture flow in soils

A computer program that simulates the flow of heat and moisture in soils is described. The space-time dependence of temperature and moisture content is described by a set of diffusion-type partial differential equations. The simulator uses a predictor/corrector to numerically integrate them, giving wetness and temperature profiles as a function of time. The simulator was used to generate solutions to diffusion-type partial differential equations for which analytical solutions are known. These equations include both constant and variable diffusivities, and both flux and constant concentration boundary conditions. In all cases, the simulated and analytic solutions agreed to within the error bounds which were imposed on the integrator. Simulations of heat and moisture flow under actual field conditions were also performed. Ground truth data were used for the boundary conditions and soil transport properties. The qualitative agreement between simulated and measured profiles is an indication that the model equations are reasonably accurate representations of the physical processes involved.

Camillo, P.↗

DOE OSTI2020

Nuclear modification factors for jet fragmentation

Using a recently-developed perturbative-QCD approach for jet evolution in a dense quark-gluon plasma, we study the nuclear modification factor for the jet fragmentation function. The qualitative behaviour that we find is in agreement with the respective experimental observations in Pb+Pb collisions at the LHC: a pronounced nuclear enhancement at both ends of the spectrum. Our Monte Carlo simulations are supplemented with analytic estimates which clarify the physical interpretation of the results. The main source of theoretical uncertainty is the sensitivity of our calculations to a low-momentum cutoff which mimics confinement. To reduce this sensitivity, we propose a new observable, which describes the jet fragmentation into subjets and is infrared-and-collinear safe by construction. We present Monte Carlo predictions for the associated nuclear modification factor together with their physical interpretation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Additive friction stir deposition induced stress ripples in aluminum alloy

In additive manufacturing, modification of the viscosity of softened or molten metal plays a significant role in the development of mechanical properties distribution. Similar to soft polymers, molten or softened metal usually exhibits vastly different viscoelastic properties. During scratching tests, viscoelastic polymers undergo spatiotemporal stress development and relaxation, causing stress ripples to form on their surfaces. Here, it is observed that friction stir-based solid-state manufacturing can yield similar stress ripples on softened aluminum alloys due to the linear motion of the deposition tool. Softened aluminum alloys exhibit viscoelastic and viscoplastic responses similar to polymeric material during the solid-state deposition process. In the final product, residual stress ripples appear on the surface of the metal alloy due to periodic residual stress relaxation, leading to oscillatory behavior from the bulk mechanical properties of the material. The orientation of the in-built stress pattern formed in the bulk of the additively manufactured metal alloy underneath the visually evident onion ring-like surface pattern formation is diametrically in the opposite direction. The experimentally observed results are consistent with analytical calculation and numerical simulations in terms of behaviors.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Tuning crack-inclusion interaction with an applied T -stress

The interaction between cracks and inclusions plays an important role in the fracture behavior of particulate composites. It is commonly recognized that an inclusion stiffer than the matrix tends to deflect an approaching crack away while a softer inclusion attracts the crack. Here, we demonstrate by analytical modeling and numerical simulations that the crack-inclusion interaction can be tuned by an applied T-stress. Under a sufficiently large compressive applied T-stress, cracks can be attracted to stiffer inclusions while repelled by softer ones, therefore reversing the conventional trend. Potential applications of this work include composite electrodes in lithium-ion batteries and hydraulic fracturing.

36 MATERIALS SCIENCE↗

DOE OSTI2023

The Princeton Field-Reversed Configuration for Compact Nuclear Fusion Power Plants

The Princeton Field-Reversed Configuration (PFRC) nuclear fusion reactor concept is an innovative approach to fusion power generation prioritizing low neutron production and small size. Here, a combination of analytical modeling and numerical simulation shows that the novel heating approach generates an FRC with closed field lines. Simulation data from a single-particle Hamiltonian code predicts ms-scale plasma heating in reactor-scale conditions while PIC codes predict formation of warm FRC plasmas from initial mirror fields. The PFRC-1 and PFRC-2 experiments have heated electrons to energies well in excess of 100 eV and plasma durations to 300 ms, more than 10 4 times longer than the predicted tilt instability growth time. From these data, we have created a development plan and anticipated performance metrics for a fusion reactor based on the PFRC concept. The resulting 1–10 MW PFRC reactors would be suitable for diverse applications, from submarines to urban environments to space propulsion. PFRC is a steady-state, driven magnetic confinement device. Plasma, inside a cylindrical array of coils, is confined and heated by external RF antennae. PFRC would be ultra-low radiation due to both its fuel and small size. The choice of advanced fuels, deuterium and helium-3 (D– 3 He), may be enabled by the high-β FRC configuration. The small size of the reactor would enable rapid exhaust of the dangerous tritium ash. Low radiation would make the reactor safer to operate and, in combination with simple geometry and small size, dramatically lowers development and maintenance costs. This review paper gives an introduction to the physics of the PFRC and a summary of the PFRC-2 experiment results to date. It then discusses the future program plan and how PFRC reactors would be commercialized.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Non-dimensional performance and safety parameters for heat pipes

The use of heat pipes in safety-critical systems such as nuclear microreactors dictates the development of generalized, practical, scalable performance and safety parameters. Traditional dimensional metrics, while informative, lack the universality required for comparative analysis across varying designs and operating regimes. Here, this work introduces a comprehensive set of non-dimensional parameters to characterize heat pipe performance and safety, including capillary performance, effective thermal conductivity, response time, exergetic efficiency, allowable temperature gradients, allowable rate of temperature change, priming coefficients, and factor of safety. A reference heat pipe design representative of microreactor applications was analyzed via the developed parameters using both traditional analytical models and Sockeye simulations under transient and steady-state conditions. Sodium, potassium, and water were evaluated as working fluids to demonstrate the applicability of the framework across a broad temperature range. The proposed non-dimensional parameters effectively captured key thermal-hydraulic behaviors and safety concerns, as was demonstrated via Sockeye simulations. This framework supports the development of design optimization strategies, operational protocols, and safety assurance practices for advanced reactor systems and other high-reliability applications.

42 - ENGINEERING↗

DOE OSTI2020

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Investigation of nonlinear mechanical behavior of two superfine-grained graphites with DIC assisted disc splitting test

Using standardized uniaxial tensile test specimens is not practical due to the limited volumes in irradiation capsules, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed to provide a convenient way for estimating tensile strength. Unlike the traditional uniaxial tensile test in ASTM C749 (International Standard Test Method for Tensile Stress–Strain of Carbon and Graphite), which uses dog-bone shaped specimens larger than 12.95 mm × 120.65 mm, the ASTM D8289 standard uses smaller discs with diameters of 6–12.7 mm. A digital image correlation (DIC) system, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 disc splitting test on IG-110 samples and compared with previous measurements from Mersen 2114 samples. Results confirmed that the DIC technique can measure the surface displacement/strain on these small (Ø6 mm × 3 mm) graphite specimens with good repeatability. However, Mersen 2114 and IG-110 samples exhibited strain discrepancies when DIC measurements were compared with analytical and finite element simulation values. The loading history and strain results also indicated different mechanical behaviors between Mersen 2114 and IG-110, particularly the nonlinear behavior of the IG-110 samples. Good agreement was observed by comparing the splitting tensile strengths of two superfine-grained grades with results from other work. The specimen size effect is discussed when comparing the splitting tensile strength with corresponding uniaxial tensile strength of these two graphite grades.

Lin, Lianshan [ORNL] (ORCID:0000000203399219)↗

DOE OSTI2022

$\mathrm{SIM}$-$\mathrm{STEM}$ $\mathrm{Lab}$: Incorporating Compressed Sensing Theory for Fast STEM Simulation

We report recently it has been shown that precise dose control and an increase in the overall acquisition speed of atomic resolution scanning transmission electron microscope (STEM) images can be achieved by acquiring only a small fraction of the pixels in the image experimentally and then reconstructing the full image using an inpainting algorithm. In this paper, we apply the same inpainting approach (a form of compressed sensing) to simulated, sub-sampled atomic resolution STEM images. We find that it is possible to significantly sub-sample the area that is simulated, the number of g-vectors contributing the image, and the number of frozen phonon configurations contributing to the final image while still producing an acceptable fit to a fully sampled simulation. Here we discuss the parameters that we use and how the resulting simulations can be quantifiably compared to the full simulations. As with any Compressed Sensing methodology, care must be taken to ensure that isolated events are not excluded from the process, but the observed increase in simulation speed provides significant opportunities for real time simulations, image classification and analytics to be performed as a supplement to experiments on a microscope to be developed in the future

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Real-time insight into the multistage mechanism of nanoparticle exsolution from a perovskite host surface

In exsolution, nanoparticles form by emerging from oxide hosts by application of redox driving forces, leading to transformative advances in stability, activity, and efficiency over deposition techniques, and resulting in a wide range of new opportunities for catalytic, energy and net-zero-related technologies. However, the mechanism of exsolved nanoparticle nucleation and perovskite structural evolution, has, to date, remained unclear. Herein, we shed light on this elusive process by following in real time Ir nanoparticle emergence from a SrTiO 3 host oxide lattice, using in situ high-resolution electron microscopy in combination with computational simulations and machine learning analytics. We show that nucleation occurs via atom clustering, in tandem with host evolution, revealing the participation of surface defects and host lattice restructuring in trapping Ir atoms to initiate nanoparticle formation and growth. These insights provide a theoretical platform and practical recommendations to further the development of highly functional and broadly applicable exsolvable materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Generation of megatesla magnetic fields by intense-laser-driven microtube implosions

A microtube implosion driven by ultraintense laser pulses is used to produce ultrahigh magnetic fields. Due to the laser-produced hot electrons with energies of mega-electron volts, cold ions in the inner wall surface implode towards the central axis. By pre-seeding uniform magnetic fields on the kilotesla order, the Lorenz force induces the Larmor gyromotion of the imploding ions and electrons. Due to the resultant collective motion of relativistic charged particles around the central axis, strong spin current densities of ~ peta-ampere/cm 2 are produced with a few tens of nm size, generating megatesla-order magnetic fields. The underlying physics and important scaling are revealed by particle simulations and a simple analytical model. The concept holds promise to open new frontiers in many branches of fundamental physics and applications in terms of ultrahigh magnetic fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Anisotropic MagnetoMemristance

Abstract In the last decade, nanoscale resistive devices with memory have been the subject of intense study because of their possible use in brain-inspired computing. However, operational endurance is one of the limiting factors in the adoption of such technology. For this reason, we discuss the emergence of current-induced memristance in magnetic materials, known for their durability. We show analytically and numerically that a single ferromagnetic layer can possess GHz memristance, due to a combination of two factors: a current-induced transfer of angular momentum (Zhang-Li torque) and the anisotropic magnetoresistance (AMR). We term the resulting effect the anisotropic magneto-memristance (AMM). We connect the AMM to the topology of the magnetization state, within a simple model of a one-dimensional annulus-shaped magnetic layer, confirming the analytical results with micromagnetic simulations for permalloy. Our results open a new path towards the realization of single-layer magnetic memristive devices operating at GHz frequencies.

97 MATHEMATICS AND COMPUTING↗