Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “physical processes”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Alpha

What is this parameter called alpha, and why is it so important? To answer these questions we consider the density of fission neutrons at any point $_r^→$ = (x, y, z) inside a volume of fissionable material. Ignoring neutron absorption, this time-varying density of neutrons is determined by two physical processes --- generation of more neutrons by fission chain reactions, and the loss of neutrons by diffusion to and through the surface. A simple, approximate mathematical description of these processes is the diffusion equation for the average neutron density n($_r^→$, t) including a fission source term.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Enabling Efficient Surveillance, Control, and Automation of Geothermal Operations with Advanced Predictive Analytics

Automation and control of geothermal energy production and operations require reliable and efficient predictive tools. While physics-based simulation offers a comprehensive tool for predicting energy production performance in geothermal systems, predicting the behavior of geothermal reservoirs involves complex multi-physics processes with coupling effects, highly uncertain input parameters and subsurface descriptions. Moreover, building, running, and integrating simulation models into standard model calibration and optimization workflows entail significant technical and computational efforts. An emerging alternative to physics-based simulation is data-driven predictive analytics models that have gained popularity in energy industry. In this report, we develop novel predictive models for integration into real-time fault diagnosis and model predictive control algorithms to improve the efficiency of energy production operations in geothermal reservoirs. The report includes two major research Thrust Areas, that is, the surface power plant and the subsurface reservoir.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

A Generalized Grain-Scale Model for the Non-Plasma and Plasma-Assisted Hydrogen Direct Reduction of Iron Ore

Direct Reduction of Iron ore using hydrogen (H-DRI) is a promising pathway towards efficient steelmaking and accurate predictive models are a necessity for scale-up and optimization of this technology. However, accurate models of this process remain limited because existing models oversimplify grain-scale phenomena, such as nonlinearity inside grain, self-sufficient porosity, surface reactions, and the role of plasma species. These phenomena are important for flash steelmaking and plasma-assisted H-DRI processes. To address this need, we present a phenomenological model for simulating H-DRI at the scale of a single micron-sized grain of the iron ore. We call this the Transient Reactive Grain Model (TRGM). TRGM incorporates key physical process: gas species transport, a chemical kinetics of material conversion, nanopore structural evolution and, adsorption-desorption surface kinetics at the reactive nanopore surface. The important contribution of this work is that the model provides a dependence on different reductant species, specifically hydrogen atoms versus molecules, so that role of hydrogen plasma reduction can be clarified compared to the use of pure hydrogen gas reduction. TRGM predictions agree well with experimental data for both molecular H2 reduction of Fe2O3 and plasma hydrogen reduction of Fe3O4. Results reveal species concentration gradients with a diffuse reaction zone, and enhanced hydrogen diffusion at the grain outer surface due to evolving porosity. These findings challenge common assumptions in existing models, including sharp reaction fronts, quasi-steady diffusion and kinetics, and the neglect of surface chemistry. As a generalized grain-scale model for H-DRI processes, TRGM has practical applications in flash steelmaking and in-flight reduction using both molecular and plasma hydrogen.

08 HYDROGEN↗
DOE OSTI2025

Advancing simulations of coupled electron and phonon nonequilibrium dynamics using adaptive and multirate time integration

Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.

Yao, Jia [California Institute of Technology (CalT↗
DOE OSTI2023

Exploring the parameter space of MagLIF implosions using similarity scaling. II. Current scaling

Magnetized liner inertial fusion (MagLIF) is a magneto-inertial-fusion (MIF) concept, which is presently being studied on the Z pulsed power facility. The MagLIF platform has achieved interesting plasma conditions at stagnation and produced significant fusion yields in the laboratory. Given the relative success of MagLIF, there is a strong interest to scale the platform to higher peak currents. However, scaling MagLIF is not entirely straightforward due to the large dimensionality of the experimental input parameter space and the numerous physical processes involved in MIF implosions. Here, in this work, we propose a novel method to scale MagLIF loads to higher currents. Our method is based on similarity (or similitude) scaling and attempts to preserve much of the physics regimes already known or being studied on today's Z pulsed-power driver. By avoiding significant deviations into unexplored and/or less well-understood regimes, the risk of unexpected outcomes on future scaled-up experiments is reduced. Using arguments based on similarity scaling, we derive the scaling rules for the experimental input parameters characterizing a MagLIF load (as functions of the characteristic current driving the implosion). We then test the estimated scaling laws for various metrics measuring performance against results of 2D radiation–magneto-hydrodynamic hydra simulations. Agreement is found between the scaling theory and the simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The Role of In-Cloud Wet Removal in Simulating Aerosol Vertical Profiles and Cloud Radiative Forcing

Among the physical processes controlling aerosol vertical profiles, in-cloud wet removal is of utmost importance while its representation in global climate models (GCMs) is crude. In this study, we implement into the Community Atmosphere Model version 6 (CAM6) a physically based aerosol wet removal parameterization scheme that explicitly treats aerosol activation, removal and resuspension. Evaluation against in-situ observations shows that the default scheme substantially overestimates the upper tropospheric black carbon (BC) and sea salt mass. Our physically based scheme reduces BC and sea salt mass by a factor of 10 and 1,000, respectively, in better agreement with observations. Also, the new scheme slightly increases number of aerosol particles between 12 nm and 4.8 μm in diameter, thereby mitigating the aerosol number underestimation in the default scheme. Our new scheme reduces the overestimation of coarse-mode aerosol (0.5–4.8 μm) number. Overall, the aerosol property changes (mass decrease and number increase) reduce the cloud condensation nuclei (CCN) concentration at low supersaturation (i.e., 0.02% and 0.1%), and increase CCN at high supersaturations (i.e., 0.5% and 1%). Consequently, the global annual mean cloud liquid water path increases by 1.89 g m -2 and the ice water path increases by 0.51 g m -2 . The global annual mean shortwave, longwave, and net cloud radiative forcing change by -1.06, 0.57, and -0.48 W m -2 , respectively. Further improvement is needed to reflect the real physics that the removal efficiencies for aerosol mass and number are disproportionate and to advect cloud-borne (activated) aerosols for a complete aerosol lifecycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Broadband multiwavelength properties of the archetypal blazar 3C 279 during the 2017 Event Horizon Telescope campaign

The archetypal blazar 3C 279 hosts a prominent relativistic jet and exhibits strong broadband variability across the electromagnetic spectrum. In April 2017, the Event Horizon Telescope (EHT) observed 3C 279, alongside one of the most extensive quasi-simultaneous multiwavelength (MWL) campaigns ever conducted. With the aim of investigating the physical processes governing 3C 279, we analyzed individual observations and multiband light curves, and constructed a new quasi-simultaneous MWL spectrum. We also performed phenomenological modeling using the turbulent extreme multi-zone (TEMZ) model to constrain the fundamental physical properties of the source. The EHT observations reveal a clear flux increase in the innermost core between April 5 and 11, 2017. Over a broader timescale, radio measurements at longer wavelengths show concurrent enhancements in core flux and polarization around mid-April, coinciding with the ejection of a superluminal knot. Record UV-optical flares with strong polarization variability occurred in late March, followed by gamma-ray activity that declined before the end of the EHT observing period. During this interval, the source remained in a low X-ray state and showed no detectable VHE emission. The TEMZ modeling suggests that the broadband spectrum and variability of 3C 279 can be explained within a jet scenario in which turbulent plasma cells are compressed by a stationary conical shock. However, alternative interpretations, such as magnetic reconnection or a moving shock-in-jet event, remain plausible. This coordinated MWL campaign advances our understanding of the origin of jet and gamma-ray emission in 3C 279, while also providing a comprehensive publicly available dataset that will serve as a valuable reference for future studies.

Principe, G. [Trieste U.; INFN, Trieste; Bologna, ↗
DOE OSTI2022

Correlative image learning of chemo-mechanics in phase-transforming solids

We report constitutive laws underlie most physical processes in nature. However, learning such equations in heterogeneous solids (for example, due to phase separation) is challenging. One such relationship is between composition and eigenstrain, which governs the chemo-mechanical expansion in solids. Here we developed a generalizable, physically constrained image-learning framework to algorithmically learn the chemo-mechanical constitutive law at the nanoscale from correlative four-dimensional scanning transmission electron microscopy and X-ray spectro-ptychography images. We demonstrated this approach on Li X FePO 4 , a technologically relevant battery positive electrode material. We uncovered the functional form of the composition–eigenstrain relation in this two-phase binary solid across the entire composition range (0 ≤ X ≤ 1), including inside the thermodynamically unstable miscibility gap. The learned relation directly validates Vegard’s law of linear response at the nanoscale. Our physics-constrained data-driven approach directly visualizes the residual strain field (by removing the compositional and coherency strain), which is otherwise impossible to quantify. Heterogeneities in the residual strain arise from misfit dislocations and were independently verified by X-ray diffraction line profile analysis. Our work provides the means to simultaneously quantify chemical expansion, coherency strain and dislocations in battery electrodes, which has implications on rate capabilities and lifetime. Broadly, this work also highlights the potential of integrating correlative microscopy and image learning for extracting material properties and physics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Polar primary aerosols across the ocean-sea ice-snow-atmosphere interface: From sources to impacts

Primary aerosols play a critical role in polar climate systems, influencing cloud formation, precipitation, radiative balance, and surface energy budgets. This paper provides a comprehensive synthesis of primary aerosol sources, transformation and removal processes, and broader atmospheric impacts in polar regions, emphasizing their links to ocean and sea ice biogeochemistry. These aerosols (including sea salt, primary organic aerosol, and primary biological aerosol particles) originate from marine and cryospheric environments and are emitted through physical processes, such as wave breaking, bubble bursting, and blowing snow. Emission sources include seawater, sea ice, snow, and freshwater from river discharge and glacial runoff. Once airborne, these particles can serve as a chemical reservoir, influencing atmospheric composition and reactivity, and as seeds for cloud droplet and ice crystal formation, influencing cloud microphysics and polar climate. Despite their importance, many of the processes governing primary aerosol emissions and transformations remain poorly constrained. The most pressing knowledge gaps pertain to emission processes, limited spatiotemporal observational coverage, instrumentation constraints, parameterization development, and the integration of interdisciplinary expertise. To improve our understanding of primary aerosol drivers and their response to climate, future research efforts should prioritize strategically coordinated and cross-disciplinary process studies, advancements in measurement technologies and coverage, and close collaboration between modelers and observational scientists to inform and refine model parameterizations. As polar regions continue to undergo profound changes marked by increased precipitation, reduced sea and land ice, freshening oceans, and shifting ecosystem dynamics, characterizing present-day primary aerosol populations is vital. Improved understanding will be essential for anticipating future changes in aerosol-radiation and aerosol-cloud interactions and their implications for polar and global climate systems.

Aerosol-cloud↗
DOE OSTI2024

A comprehensive numerical investigation on spray models for Direct-Injection Spark-Ignition engines

Gasoline direct-injection spark-ignition (DISI) engines generate a large portion of their unburned hydrocarbon (UHC) and soot emissions during the cold-start phase. A predictive computational fluid dynamics (CFD) modeling framework can be used to understand the physical processes that characterize fuel spray evolution and fuel-film formation at cold start conditions, which can help to reduce engine-out particulate emissions. This study systematically evaluated spray submodels and developed a set of simulation best practices for physical-numerical submodels with the goal of enabling accurate simulations of liquid spray behavior in a DISI engine. Three comprehensive experimental datasets containing free-spray projected liquid volume (PLV), liquid volume fraction (LVF), and near-field X-ray radiography data were used to validate the simulation results and evaluate the spray submodels. Systematic analysis delved into injected parcel distribution, droplet collision, spray breakup, and evaporation via a detailed assessment of the relevant spray submodels. Moreover, the effects of turbulence models and the initial turbulent flow properties on the liquid spray evolution were examined. Based on extensive calibration efforts, a set of simulation best practices for the free spray was developed and validated against the PLV/LVF data. Simulation results indicated that the uniform distribution for parcel initialization, coupled with appropriate droplet collision submodels, provides an improved spray morphology compared to the cluster distribution. The findings also underscored the importance of calibrating the Kelvin-Helmholtz Rayleigh-Taylor (KH-RT) breakup model constants and droplet heat transfer coefficient scaling factor to achieve favorable agreement regarding measured liquid penetration and spray widths. In conclusion, this study marks a substantial stride towards accurately predicting fuel film evolution and soot formation within DISI engine performance.

ECN Spray G↗
DOE OSTI2022

Overview of the COMPASS results *

COMPASS addressed several physical processes that may explain the behaviour of important phenomena. This paper presents results related to the main fields of COMPASS research obtained in the recent two years, including studies of turbulence, L–H transition, plasma material interaction, runaway electron, and disruption physics: Tomographic reconstruction of the edge/SOL turbulence observed by a fast visible camera allowed to visualize turbulent structures without perturbing the plasma. Dependence of the power threshold on the X-point height was studied and related role of radial electric field in the edge/SOL plasma was identified. The effect of high-field-side error fields on the L–H transition was investigated in order to assess the influence of the central solenoid misalignment and the possibility to compensate these error fields by low-field-side coils. Results of fast measurements of electron temperature during ELMs show the ELM peak values at the divertor are around 80% of the initial temperature at the pedestal. Liquid metals were used for the first time as plasma facing material in ELMy H-mode in the tokamak divertor. Good power handling capability was observed for heat fluxes up to 12 MW m –2 and no direct droplet ejection was observed. Partial detachment regime was achieved by impurity seeding in the divertor. The evolution of the heat flux footprint at the outer target was studied. Runaway electrons were studied using new unique systems—impact calorimetry, carbon pellet injection technique, wide variety of magnetic perturbations. Radial feedback control was imposed on the beam. Forces during plasma disruptions were monitored by a number of new diagnostics for vacuum vessel (VV) motion in order to contribute to the scaling laws of sideways disruption forces for ITER. Current flows towards the divertor tiles, incl. possible short-circuiting through PFCs, were investigated during the VDE experiments. The results support ATEC model and improve understanding of disruption loads.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

1D model of tokamak scrape-off layer and divertor plasmas: Dimensionless parameters and the results of numerical simulations

Solutions of a 1D model of the scrape-off layer and a divertor plasma in a tokamak are analyzed on the basis of dimensionless parameters governing the physics of both plasma and neutral transport and atomic physics processes. Numerical simulations based on fluid plasma and Monte Carlo neutral modeling confirm earlier results obtained with simplified analytic considerations and numerical simulations performed with fluid neutral description and demonstrate bifurcation of plasma parameters. Further, these simulations also confirm that the major physical reasons for the rollover of the plasma flux to the target, which is the key feature of the phenomenon of divertor plasma detachment, are the impurity radiation loss and plasma recombination. It is suggested that simulations of the edge plasma based on dimensionless parameters provide a powerful tool for checking the consistency of the physics incorporated in existing complex 2D codes used for edge plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Unsupervised quantum circuit learning in high energy physics

Unsupervised training of generative models is a machine learning task that has many applications in scientific computing. Here, in this work, we evaluate the efficacy of using quantum circuit-based generative models to generate synthetic data of high energy physics processes. We use nonadversarial, gradient-based training of quantum circuit Born machines to generate joint distributions over two and three variables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Understanding Processes Controlling the Temporal and Spatial Variations of PBL Structures Over the ARM SGP Site

The surface heat, moisture, and momentum fluxes are transferred to the atmosphere above through the planetary boundary layer (PBL), where vertical mixing due to turbulent eddies of different sizes plays critical roles. Therefore, reliably representing PBL processes in numerical models is critical for weather, climate, and air quality prediction. Currently, there are over ten PBL schemes that are selectable within the advanced research version of the Weather Research and Forecasting (WRF) model, indicative of the challenges in capturing the impacts of turbulence within the PBL in models. Further improvements in PBL parameterizations are needed for both weather and climate models, as emphasized in many recent national reports, but require an advanced understanding of the underlying boundary layer processes from observations. This project takes advantage of Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) investments in the atmospheric boundary layer observations and Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) simulations to characterize PBL structures, understand key physical processes controlling the mixed layer development, and to evaluate PBL parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Improving Gas Reactor Design with Complex Non-Standard Reaction Mechanisms in a Reactive Flow Model

Fluorinated gases are critical to a number of high-technology industries including semiconductor manufacturing and related processes such as flat panel display, LED, and photovoltaic manufacturing. These fluorinated gases have very high global warming potentials (GWPs) and, unfortunately are chemically very stable, in some cases lasting for thousands of years in the atmosphere. For environmental reasons, it is imperative that the release of these gases be minimized, but existing methods to destroy them requires a significant amount of energy. The goal of this project was to develop a more energy efficient method of destroying these high GWP gases using computational fluid dynamics, but the computational requirements are significant. In this project, we bring together resources available at the national labs, including high-performance computing hardware and modern computational methods, in order to explore an industrially important problem with key environmental and economic impact. In this project, we examined the state-of-the-art models for the physical processes involved, and made subtle changes to arrive at alternative formulations considerably better suited for computation and study. We examined several representative scenarios and developed a set of analyses techniques and tools that could be used to explore the systems in much greater detail. This exploratory project has demonstrated to ALZETA that HPC can indeed provide a unique view into their complex engineering problem - one that is simply inaccessible by any other means. The tools allow us to explore the details of the chemical process, and allow one to construct "what-if" scenarios to address shortcomings of the physical device in its current operating modes. However, as with many sources of completely new information, many new questions are illuminated with this capability. It will take considerable time, effort and experience before we learn to probe the solutions, pose operating scenarios and confirm modifications to the device that lead to real engineering advances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Low-dose high-resolution $\mathrm{TOF}$-$\mathrm{PET}$ using ionization-activated multi-state low-$\mathrm{Z}$ detector media

We propose PET scanners using low atomic number media that undergo a persistent local change of state along the paths of the Compton recoil electrons. Measurement of the individual scattering locations and angles, deposited energies, and recoil electron directions allows using the kinematical constraints of the 2-body Compton scattering process to perform a statistical time-ordering of the scatterings, with a high probability of precisely identifying where the gamma first interacted in the detector. In these cases the Line-of-Response is measured with high resolution, determined by the underlying physics processes and not the detector segmentation. There are multiple such media that act through different mechanisms. As an example in which the change of state is quantum-mechanical through a change in molecular configuration, rather than thermodynamic, as in a bubble chamber, we present simulations of a two-state photoswitchable organic dye, a ‘Switchillator’, that is activated to a fluorescent-capable state by the ionization of the recoil electrons. The activated state is persistent, and can be optically excited multiple times to image individual activated molecules. Energy resolution is provided by counting the activated molecules. Location along the LOR is implemented by large-area time-of-flight MCP-PMT photodetectors with single photon time resolution in the tens of ps and sub-mm spatial resolution. Simulations indicate a large reduction of dose.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗