Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Illuminating dark matter halo density profiles without subhaloes

Cold dark matter haloes consist of a relatively smooth dark matter component as well as a system of bound subhaloes. It is the prevailing practice to include all mass, including mass in subhaloes, in studies of halo density profiles in simulations. However, often in observational studies satellites are treated as having their own distinct dark matter density profiles in addition to the profile of the host. This difference can make comparisons between theoretical and observed results difficult. In this work, we investigate density profiles of the smooth components of host haloes by excluding mass contained within subhaloes. We find that the density profiles of the smooth halo component (without subhaloes) differ substantially from the conventional halo density profile, declining more rapidly at large radii. We also find that concentrations derived from smooth density profiles exhibit less scatter at fixed mass and a weaker mass dependence than standard concentrations. Both smooth and standard halo profiles can be described by a generalized Einasto profile, an Einasto profile with a modified central slope, with smaller residuals than either a Navarro–Frenk–White or Einasto profile. Furthermore, these results hold for both Milky Way-mass and cluster-mass haloes. This new characterization of smooth halo profiles can be useful for many analyses, such as lensing and dark matter annihilation, in which the smooth and clumpy components of a halo should be accounted for separately.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Implementation of a Model Predictive Control Strategy to Regulate Temperature Inside Plug-Flow Solar Reactor With Countercurrent Flow

Abstract Solar-driven thermochemical energy storage systems are proven to be promising energy carriers (solar fuels) to utilize solar energy by using reactive solid-state pellets. However, the production of solar fuel requires a quasi-steady-state process temperature, which represents the main challenge due to the transient nature of solar power. In this work, an adaptive model predictive controller (MPC) is presented to regulate the temperature inside a tubular solar reactor to produce solid-state solar fuel for long-term thermal storage systems. The solar reactor system consists of a vertical tube heated circumferentially over a segment of its length by concentrated solar power, and the reactive pellets (MgMn2O4) are fed from the top end and flow downwards through the heated tube. A countercurrent flowing gas supplied from the lower end interacts with flowing pellets to reduce it thermochemically at a temperature range of 1000—1500 °C. A low-order physical model was developed to simulate the dynamics of the solar reactor including the reaction kinetics, and the proposed model was validated numerically by using a 7-kW electric furnace. The numerical model then was utilized to design the MPC controller, where the control system consists of an MPC code linked to an adaptive system identification code that updates system parameters online to ensure system robustness against external disturbances (sudden change in the flow inside the reactor), model mismatches, and uncertainty. The MPC controller parameters are tuned to enhance the system performance with minimum steady-state error and overshoot. The controller is tested to track different temperature ranges between 500 °C and 1400 °C with different particles/gas mass flowrates and ramping temperature profiles. Results show that the MPC controller successfully regulated the reactor temperature within ± 1 °C of its setpoint and maintained robust performance with minimum input effort when subjected to sudden changes in the amount of flowing media and the presence of chemical reaction.

Engineering↗
DOE OSTI2023

Interface models for three-dimensional Rayleigh–Taylor instability

We derive interface models for three-dimensional Rayleigh–Taylor instability (RTI), making use of a novel asymptotic expansion in the non-locality of the fluid flow. These interface models are derived for the purpose of studying universal features associated with RTI such as the Froude number in single-mode RTI, the predicted quadratic growth of the interface amplitude under multi-mode random perturbations, the optimal (viscous) mixing rates induced by the RTI and the self-similarity of horizontally averaged density profiles and the remarkable stabilization of the mixing layer growth rate which arises for the three-fluid two-interface heavy–light–heavy configuration, in which the addition of a third fluid bulk slows the growth of the mixing layer to a linear rate. Our interface models can capture the formation of small-scale structures induced by severe interface roll-up, reproduce experimental data in a number of different regimes and study the effects of multiple interface interactions even as the interface separation distance becomes exceedingly small. Compared with traditional numerical schemes used to study such phenomena, our models provide a computational speed-up of at least two orders of magnitude.

Mechanics↗
DOE OSTI2020

Correlations between triaxial shapes and formation history of dark matter haloes

ABSTRACT The shape of dark matter (DM) haloes plays a critical role in constraining cosmology with upcoming large-scale structure surveys. In this paper, we study the correlations between the triaxial shapes and formation histories in DM haloes in the MultiDark Planck 2 N-body cosmological simulation. We find that halo ellipticity is strongly correlated with halo properties that serve as proxies of halo formation history, such as halo concentration and the normalized peak-centroid offset. These correlations are nearly independent of the halo density peak height. We present a simple model for the correlation between halo ellipticity and concentration using conditional abundance matching, and provide fitting formulae for the multidimensional distributions of triaxial halo shape as a function of halo peak height. We apply our halo shape model to gauge the effects of halo ellipticity and orientation bias on the excess surface mass density profiles in cluster-sized haloes. Our model should be useful for exploring the impact of triaxial halo shape on cosmological constraints in upcoming weak lensing surveys of galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Three‐Dimensional Magnetic Reconnection Spreading in Current Sheets of Non‐Uniform Thickness

Abstract Magnetic reconnection in naturally occurring and laboratory settings often begins locally and elongates, or spreads, in the direction perpendicular to the plane of reconnection. Previous work has largely focused on current sheets with a uniform thickness, for which the predicted spreading speed for anti‐parallel reconnection is the local speed of the current carriers. We derive a scaling theory of three‐dimensional (3D) spreading of collisionless anti‐parallel reconnection in a current sheet with its thickness varying in the out‐of‐plane direction, both for spreading from a thinner to thicker region and a thicker to thinner region. We derive an expression for calculating the time it takes for spreading to occur for a current sheet with a given profile of its thickness. A key result is that when reconnection spreads from a thinner to a thicker region, the spreading speed in the thicker region is slower than both the Alfvén speed and the speed of the local current carriers by a factor of the ratio of thin to thick current sheet thicknesses. This is important because magnetospheric and solar observations have previously measured the spreading speed to be slower than previously predicted, so the present mechanism might explain this feature. We confirm the theory via a parametric study using 3D two‐fluid numerical simulations. We use the prediction to calculate the time scale for reconnection spreading in Earth's magnetotail during geomagnetic activity. The results are also potentially important for understanding reconnection spreading in solar flares and the dayside magnetopause of Earth and other planets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nonlinear dimension reduction for surrogate modeling using gradient information

We introduce a method for the nonlinear dimension reduction of a high-dimensional function $u:{\mathbb{R}}^d\rightarrow{\mathbb{R}}$, $d\gg 1$. Our objective is to identify a nonlinear feature map $g:{\mathbb{R}}^d\rightarrow{\mathbb{R}}^m$, with a prescribed intermediate dimension $m\ll d$, so that $u$ can be well approximated by $f\circ g$ for some profile function $f:{\mathbb{R}}^m\rightarrow{\mathbb{R}}$. We propose to build the feature map by aligning the Jacobian $\nabla g$ with the gradient $\nabla u$, and we theoretically analyze the properties of the resulting $g$. Once $g$ is built, we construct $f$ by solving a gradient-enhanced least squares problem. Our practical algorithm uses a sample $\{{\textbf{x}}^{(i)},u({\textbf{x}}^{(i)}),\nabla u({\textbf{x}}^{(i)})\}_{i=1}^N$ and builds both $g$ and $f$ on adaptive downward-closed polynomial spaces, using cross validation to avoid overfitting. We numerically evaluate the performance of our algorithm across different benchmarks, and explore the impact of the intermediate dimension $m$. We show that building a nonlinear feature map $g$ can permit more accurate approximation of $u$ than a linear $g$, for the same input data set.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Case For Intra-rack Resource Disaggregation in HPC

The expected halt of traditional technology scaling is motivating increased heterogeneity in high-performance computing (HPC) systems with the emergence of numerous specialized accelerators. As heterogeneity increases, so does the risk of underutilizing expensive hardware resources if we preserve today’s rigid node configuration and reservation strategies. This has sparked interest in resource disaggregation to enable finer-grain allocation of hardware resources to applications. However, there is currently no data-driven study of what range of disaggregation is appropriate in HPC. To that end, we perform a detailed analysis of key metrics sampled in NERSC’s Cori, a production HPC system that executes a diverse open-science HPC workload. In addition, we profile a variety of deep-learning applications to represent an emerging workload. We show that for a rack (cabinet) configuration and applications similar to Cori, a central processing unit with intra-rack disaggregation has a 99.5% probability to find all resources it requires inside its rack. In addition, ideal intra-rack resource disaggregation in Cori could reduce memory and NIC resources by 5.36% to 69.01% and still satisfy the worst-case average rack utilization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Numerical investigation of alpha particle confinement under the perturbation of neoclassical tearing modes and toroidal field ripple in CFETR

The confinement of alpha particles in burning plasma is a key issue in fusion reactor design, including particle interaction with instabilities. Here, we include two topics: the effect of neoclassical tearing modes (NTMs) and toroidal field ripple on alpha particle loss, and the assessment of particle redistribution under an NTM with a reduced model. We consider Chinese fusion engineering test reactor parameters, the alpha particle distribution given by TRANSP/NUBEAM and the NTM perturbation function given by the initial value code TM1. We show that the synergistic effect of the NTM and ripple is negligible; the particle loss fraction does not change with increasing NTM amplitude. Only passing particles are affected by the mode particle resonance, producing profile flattening but no increased loss because only trapped particles are influenced by ripple. To study alpha particle profile flattening, the work adopts an innovative method of phase vector rotation to determine regions of good and broken Kolmogorov–Arnold–Moser surfaces and equilibrates the particle density according to local stochasticity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Reducing Uncertainty in Offshore Wind Energy Yield Estimates via a Metocean Reference Site

The offshore wind industry is burgeoning in the coastal waters of the United States, specifically along the Atlantic. For wind energy to be successful, reliable observations and model simulations are needed for resource assessment and forecasting. While many of these activities have already begun, there is currently an absence of observations at hub-height in these waters, with the closest available hub-height measurements usually taken onshore. Deployment of floating lidars has occurred through various federally funded projects, but only encapsulates time periods of a couple of years at best. Private industry is also beginning to leverage floating lidars, but this data is often proprietary, and not shared with the general public. In this work, we make the case for a metocean reference site for long-term offshore wind energy. Specifically, we quantify the impact of having a metocean reference site compared to other methods of determining hub-height winds and energy production. We use an offshore floating lidar to directly measure the wind resource, and compare these measurements to predictions derived from other widely-available surface meteorological variables. These prediction methods (vertical extrapolation, machine learning, and NWP output) produce a variety of vertical wind speed profiles, of which produce different energy yield estimates for a reference offshore turbine (Figure 1). While some methods perform reasonably well against the lidar, the uncertainty in these energy yield estimates has financial implications, further illustrating the need for long-term measurements in coastal waters.

machine learning↗
DOE OSTI2023

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗
DOE OSTI2024

Impact of a Bluff Body on a Pressurized Oxy-fuel Flame

Designing an effective burner is vital for the development of pressurized oxy-fuel combustion technologies. In the present work, turbulent jet diffusion burners are adopted for a pressurized oxy-combustor, with a bluff-body employed to create a recirculation zone, thereby stabilizing the flame in such a combustor. The objective of this numerical study is to perform a systematic analysis of the characteristics of such a pressurized non-premixed flame. Specifically, a 15-bar pressurized oxy-fuel combustor of power 100 kWth is modeled by means of the Ansys FLUENT commercial platform, using the Reynolds-averaged Navier-Stokes (RANS) approach. The present work focuses on identifying the aerothermodynamic features of the pressurized oxy-fuel burner with a disk-shape bluff-body. It is shown that the fuel-to-oxidizer stream momentum ratio has a great impact on the temperature profile of the down-fired, co-axial, pressurized oxy-coal, diffusion flame. A parametric study of the blockage ratio of the burner identifies an optimal range for this pressurized burner.

Li, Lei↗
DOE OSTI2021

The scaling of kinetic and transport behaviors in the solution-phase chemistry of a plasma–liquid interface

The reactions at a plasma–liquid interface often involve species such as the solvated electron or the hydroxyl radical, which initiate the reduction or oxidation of solution-phase reactants (so-called scavengers) or are consumed by their own second-order recombination. Here, the mathematical scaling of the reaction–diffusion equations at the interface is used to obtain a characteristic time that can be used to determine the transition from highly efficient scavenger reduction or oxidation to lower efficiencies due to transport limitations. The characteristic time (tc) is validated using numerical solutions of the reaction–diffusion equations. When the scavenger kinetics are faster than second-order recombination, this characteristic transition time scales proportionally with the scavenger diffusivity (Ds) and the square of the scavenger bulk concentration (SB) and inversely proportional to the electron flux (J) squared; that is, t c = D s S B 2 F 2 /J 2 , where F is Faraday’s constant. However, when the scavenger kinetics are comparable or slower than second-order recombination, this scaling does not hold. Extending this analysis to three-dimensional systems shows that the profile of the electron flux on the surface affects the spatial location where reactions are most effective. Finally, the assessment of the implications of these behaviors for the reactor design highlights how effectively controlling the electron flux and solution transport may be necessary to improve the efficiency of scavenger reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Experimental measurements of gas-particle flows in large-scale strippers

Measurements of gas–solid particle flow in a large-scale stripper unit are reported. The experiments target industrial (pilot) scale measurements of gas-particle flow that can be used to validate numerical methods for modeling multiphase flows, such as computational fluid dynamics coupled to the discrete element method (CFD-DEM). Specifically, experiments were performed with Geldart Group B (533 µm) glass beads in a stripper unit. In general, previous CFD-DEM validation datasets were performed in benchtop experimental systems and limited to O(10 5 ) particles. Here, in this work, the number of particles in the isolated section of the stripper is O(10 7 ) for all tests. Measurements are reported for a ~2 m tall section of the ~1 m diameter stripper. Radial mass flux profiles measured at five axial positions and three axial pressure drop measurements are reported for six operating conditions.

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

Controlling pollutant emissions in a high-pressure combustor with fuel-diluent blending

This paper investigates the formation of nitric oxide (NO) and carbon monoxide (CO) emissions as a function of nitrogen (N 2 ) or carbon dioxide (CO 2 ) diluent content, added to a premixed reacting jet in crossflow. Here, reaction characteristics of a rich methane-air jet injected into a lean vitiated crossflow were analyzed at an elevated pressure of 5 atm. The jet was pre-heated and enriched with 0%, 15%, 30%, and 50% mass fraction to quantify the effect on pollutant emissions. Simulated results of the full chemistry Star-CCM+CFD model were verified with data taken in the experimental high-pressure combustion facility, which provides pressure, temperature and velocity profiles, as well as line-of-sight chemiluminescence and exit emission measurements. The analysis revealed the significant influence of the diluent to delay axial combustion; this mixing delay increased the axial flame lift-off and controlled any thermal ignition limitations. Increasing the diluent content increased the timescale for the flame to stabilize, which allowed for greater entrainment of crossflow oxidizer into the axial jet stream, and led to decreased pollutant emissions. Hence, crossflow entrainment is a critical driving force at high diluent content. Local nitric oxide (NO) emission formation in the axial stage was predicted numerically, showing the correlation between diluent addition, axial heat release, and the formation of nitric oxide pollution. The high diluent levels leaned out the jet mixture and delayed axial combustion, while reducing prompt and thermal NOx by mitigating flame hotspots and minimizing the timescale that the products remain at high temperatures.

03 NATURAL GAS↗
DOE OSTI2024

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗
DOE OSTI2022

Physics Guided Simulation of Electrostatic Discharge: Technical Report

Triboelectrically-charged objects may create threshold sparks, electrostatic discharge (ESD) events, to equilibrate charge between themselves and other relatively charged objects. ESD events exhibit many complex physical phenomena. They are a nexus of several fields of physics with disparate characteristic scales: plasma physics, chemical kinetics, hydrodynamics, circuit models, etc. These scales can span many orders of magnitude from the varied collisions thermalizing information within a plasma on the $\mathcal{O}(fs/ps)$ to the physical size of the plasma channel on the $\mathcal{O}(100µm)$, to the speed of a nonlinear hydrodynamic wave propagating at $\mathcal{O}(µm, ns)$. These threshold ESD events may occur in situations of programmatic importance, delivering energy and power profiles to a “victim load” generating deleterious consequences. To predict and mitigate these consequences we must answer questions about the spark’s energy budget: how much energy goes into producing the spark channel; how much gets radiated away; how much energy is advected away into the hydrodynamics; and how much energy is delivered to a victim load. An ESD simulation toolset has been created and evolved in order to answer these questions. An appropriate, physics guided implementation for simulation can be done by gaining insight into its constituent physics and leveraging that intuition to choose a suitable numerical operator. We examine in detail the chemical kinetics, circuit discharge, and hydrodynamics to deter mine dominant regimes, values, timescales, and interactions to uncover the underlying physical dynamics. We also examine and propose model reduction schemes for high-dimensional chemical kinetics. We use past and current work with experimentally validated and theoretically-verified hydrodynamics to calculate applicability limits of the non-ionizing strong shock limit. We quantify the energy budget from a hydrodynamic perspective and demonstrate that a significant fraction of the stored energy is “earmarked” for hydrodynamic advection as an energy terminus. Lastly, we combine the constituent physics of an ESD event (chemical kinetics, circuit model, and hydrodynamics) into a cohesive, actionable toolset and obtain promising results from an isothermal test case. We then propose a viable, modular evolution of the ESD toolset based upon the performed examination of the physics uncovering dominant physical scales and the stiffness of the compositional differential system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Acute Stent-Induced Endothelial Denudation: Biomechanical Predictors of Vascular Injury

Recent concern for local drug delivery and withdrawal of the first Food and Drug Administration-approved bioresorbable scaffold emphasizes the need to optimize the relationships between stent design and drug release with imposed arterial injury and observed pharmacodynamics. In this study, we examine the hypothesis that vascular injury is predictable from stent design and that the expanding force of stent deployment results in increased circumferential stress in the arterial tissue, which may explain acute injury poststent deployment. Using both numerical simulations and ex vivo experiments on three different stent designs (slotted tube, corrugated ring, and delta wing), arterial injury due to device deployment was examined. Furthermore, using numerical simulations, the consequence of changing stent strut radial thickness on arterial wall shear stress and arterial circumferential stress distributions was examined. Regions with predicted arterial circumferential stress exceeding a threshold of 49.5 kPa compared favorably with observed ex vivo endothelial denudation for the three considered stent designs. In addition, increasing strut thickness was predicted to result in more areas of denudation and larger areas exposed to low wall shear stress. We conclude that the acute arterial injury, observed immediately following stent expansion, is caused by high circumferential hoop stresses in the interstrut region, and denuded area profiles are dependent on unit cell geometric features. Such findings when coupled with where drugs move might explain the drug–device interactions.

Conway, Claire↗
DOE OSTI2021

Effect of thermal environments on fast charging Li-ion batteries

Battery thermal management systems (BTMSs) are expected to keep the battery temperature at a moderate level (~30 °C) to minimize the thermally exacerbated degradation. However, during fast charging, a strong cooling system is required to restrict the temperature rise of Li-ion batteries (LiBs), which significantly increases the cost and weight of battery packs, and induces a large temperature variation inside the battery. Here, in this work, we find that all these drawbacks could be relieved by allowing LiBs to charge at higher temperatures. Since the fast charging of a LiB only takes a tiny fraction of its lifetime, the aging rate is limited even at a charging temperature of 60 °C. Three types of thermal environments are proposed: kept constant at 30 °C, preheated to 60 °C, and adiabatic fast charging. With an experimentally validated electrochemical-thermal (ECT) coupled model, we explore the interplay between thermal management and the fast-charging performance. It is found that a gradually increasing temperature profile is the best option to balance the lithium plating and thermal management of the battery. Combining adiabatic fast charging with a preheating step, we can achieve minimal cooling need, perfect temperature uniformity within a battery, and fast-charging capability simultaneously.

25 ENERGY STORAGE↗