Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow fields”

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

Quadrupolar density structures in driven magnetic reconnection experiments with a guide field

Magnetic reconnection is a ubiquitous process in plasma physics, driving rapid and energetic events such as coronal mass ejections. Reconnection between magnetic fields with arbitrary shear can be decomposed into an anti-parallel reconnecting component and a non-reconnecting guide-field component, which is parallel to the reconnecting electric field. This guide field modifies the structure of the reconnection layer and the reconnection rate. We present results from experiments on the MAIZE pulsed-power generator (500 kA peak current, 200 ns rise time), which use two exploding wire arrays, tilted in opposite directions, to embed a guide field in the plasma flows with a relative strength b≡B g /B rec =0, 0.4, or 1. The reconnection layers in these experiments have widths that are less than the ion skin depth, d i =c/ω pi , indicating the importance of the Hall term, which generates a distinctive quadrupolar magnetic field structure along the separatrices of the reconnection layer. Using laser imaging interferometry, we observe quadrupolar structures in the line-integrated electron density, consistent with the interaction of the embedded guide field with the quadrupolar Hall field. Our measurements extend over much larger length scales (40d i ) at higher β (∼1) than previous experiments, providing an insight into the global structure of the reconnection layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Localized Evaluation for Constructing Discrete Vector Fields

Topological abstractions offer a method to summarize the behavior of vector fields, but computing them robustly can be challenging due to numerical precision issues. One alternative is to represent the vector field using a discrete approach, which constructs a collection of pairs of simplices in the input mesh that satisfies criteria introduced by Forman's discrete Morse theory. While numerous approaches exist to compute pairs in the restricted case of the gradient of a scalar field, state-of-the-art algorithms for the general case of vector fields require expensive optimization procedures. This paper introduces a fast, novel approach for pairing simplices of two-dimensional, triangulated vector fields that do not vary in time. The key insight of our approach is that we can employ a local evaluation, inspired by the approach used to construct a discrete gradient field, where every simplex in a mesh is considered by no more than one of its vertices. Specifically, we observe that for any edge in the input mesh, we can uniquely assign an outward direction of flow. We can further expand this consistent notion of outward flow at each vertex, which corresponds to the concept of a downhill flow in the case of scalar fields. Working with outward flow enables a linear-time algorithm that processes the (outward) neighborhoods of each vertex one-by-one, similar to the approach used for scalar fields. Here, we couple our approach to constructing discrete vector fields with a method to extract, simplify, and visualize topological features. Empirical results on analytic and simulation data demonstrate drastic improvements in running time, produce features similar to the current state-of-the-art, and show the application of simplification to large, complex flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Emergent electric field induced by dissipative sliding dynamics of domain walls in a Weyl magnet

The dynamic motion of topological defects in magnets induces an emergent electric field, as exemplified by the continuous flow of skyrmion vortices. However, the electrodynamics underlying this emergent field remains poorly understood. In this context, magnetic domain walls—one-dimensional topological defects with two collective modes, sliding and spin-tilt—offer a promising platform for exploration. Here we demonstrate that the dissipative motion of domain walls under oscillatory current excitation generates an emergent electric field. We image domain patterns and quantify the domain-wall length under applied magnetic fields in mesoscopic devices based on the magnetic Weyl semimetal NdAlSi. These devices exhibit exceptionally strong domain-wall scattering and a pronounced emergent electric field, as observed in the imaginary component of the complex impedance. Spin dynamics simulations reveal that domain-wall sliding dominates over spin-tilting, in which the phase delay of the domain-wall motion with respect to the driving force impacts the emergent electric field. Our findings establish domain-wall dynamics as a platform for studying emergent electromagnetic fields and motivate further investigations of the coupled motion of magnetic solitons and conduction electrons.

Yamada, Rinsuke [The University of Tokyo, Japan]↗
DOE OSTI2026

Efficient Probabilistic Visualization of Local Divergence of 2D Vector Fields with Independent Gaussian Uncertainty

This work focuses on visualizing uncertainty of local divergence of two-dimensional vector fields. Divergence is one of the fundamental attributes of fluid flows, as it can help domain scientists analyze potential positions of sources (positive divergence) and sinks (negative divergence) in the flow. However, uncertainty inherent in vector field data can lead to erroneous divergence computations, adversely impacting downstream analysis. While Monte Carlo (MC) sampling is a classical approach for estimating divergence uncertainty, it suffers from slow convergence and poor scalability with increasing data size and sample counts. Thus, we present a two-fold contribution that tackles the challenges of slow convergence and limited scalability of the MC approach. (1) We derive a closed-form approach for highly efficient and accurate uncertainty visualization of local divergence, assuming independently Gaussian-distributed vector uncertainties. (2) We further integrate our approach into Viskores, a platform-portable parallel library, to accelerate uncertainty visualization. In our results, we demonstrate significantly enhanced efficiency and accuracy of our serial analytical (speed-up up to 1946×) and parallel Viskores (speed-up up to 19698×) algorithms over the classical serial MC approach. We also demonstrate qualitative improvements of our probabilistic divergence visualizations over traditional mean-field visualization, which disregards uncertainty. We validate the accuracy and efficiency of our methods on wind forecast and ocean simulation datasets.

Ouermi, Timbwaoga [University of Utah]↗
DOE OSTI2024

Understanding the impact of an applied axial magnetic field on efficient current coupling on the Z machine

Magnetized liner inertial fusion (MagLIF) is an attractive concept for producing thermonuclear fusion reactions. The MagLIF platform involves the operation of Helmholtz coils to apply a 15 Tesla axial magnetic field to the load region, where a cylindrical, fuel-filled metal liner is imploded by a ∼ 2 0 MA current pulse. The fringe field from these coils extends into the transmission line that delivers the current to the target. We investigated the extent to which this applied field disturbs the nominal power flow within that transmission line. A simplified model of the geometry shows that adding the applied magnetic field results in magnetic field lines that connect the cathode to the anode, suggesting electrons may not be magnetically insulated in this region. Particle-in-cell simulations indicated the addition of the applied magnetic field would not significantly impact the current delivery to the load. Velocimetry was used to experimentally assess the current delivery with and without the applied magnetic field. We find no measurable effects of the applied field on current delivery in the configuration investigated in this study. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Jet cone radius dependence of R AA and v 2 at PbPb 5.02 TeV from JEWEL+T R ENTo+v-USPhydro

We combine, for the first time, event-by-event T R ENTo initial conditions with the relativistic viscous hydrodynamic model v-USPhydro and the Monte Carlo event generator JEWEL to make predictions for the nuclear modification factor R AA and jet azimuthal anisotropies v n { 2 } in $\sqrt{s_{NN}}$ = 5.02 TeV PbPb collisions for multiple centralities and values of the jet cone radius R. The R-dependence of R AA and v 2 { 2 } strongly depends on the presence of recoiling scattering centers. We find a small jet v 3 { 2 } in mid-central collisions and consistent results in wide jet p T regions and centralities with ATLAS data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗
DOE OSTI2025

Hydrologic Model Data for the East Fork Poplar Creek Watershed Simulated with the Advanced Terrestrial Simulator (ATS): Streamflow and Network Expansion–Contraction Dynamics

This dataset supports hydrologic modeling and stream network expansion–contraction analysis for the East Fork Poplar Creek (EFPC) Watershed in Tennessee. It includes a Jupyter notebook for model setup, model configuration files, simulation outputs, and derived products used to evaluate model performance and investigate stream dynamics under varying hydrologic conditions. The dataset was generated using the Watershed Workflow Python package and the Advanced Terrestrial Simulator (ATS), enabling integrated surface–subsurface hydrologic simulations using a stream-aligned mesh. Outputs include high-resolution time series of streamflow, active network length, water table depth, and related hydrologic variables. Also included are spatially explicit stream persistency indices and classifications of reaches as perennial or non-perennial. These data facilitate reproducibility and support further research on stream intermittency and variability in network extent.The model data archive is organized in following directories:1) model_setup_inputsContains the Watershed Workflow Jupyter notebooks (accessed through any open source code editor), selected input datasets, and resulting ATS input files, including XML files (access through any open source code editor), computational mesh (.exo files can be viewed using Paraview), and meteorological forcing files (.h5 files can be accessed through h5py python package and HDFView open source software). 2) model_outputsIncludes ATS simulation outputs relevant to this study. Time series of spatially integrated or averaged variables (e.g., streamflow, water table depth) are provided as CSV files. Select spatial fields (e.g., ponded depth and water table depth) are saved as pickled Python objects to reduce file size, and can be accessed through pickle package in Python. Key geometry objects from Watershed Workflow—such as the surface mesh and river tree—are also included to support analysis of streamflow persistency and expansion–contraction dynamics. These files can also be accessed through Watershed Workflow Python package.3) model_evaluationProvides observed streamflow time series and field survey-based flow regime classifications used to evaluate model performance. Jupyter notebooks for processing ATS outputs and comparing model predictions with observations to build confidence in the model prior to scientific analysis are also included.4) Q_L_relationshipsContains workflows for generating time series of discharge, active network length, and related hydrologic variables used in the stream network expansion–contraction analysis. Includes routines for delineating baseflow-dominated periods. For each catchment, notebooks and processed data (as pickled DataFrames accessed through Pandas Python package) are provided. 5) figure_scriptsProvides the Jupyter notebooks used to generate the figures presented in the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Nonadiabatic Force Matching for Alchemical Free-Energy Estimation

We propose a method to compute free-energy differences from nonadiabatic alchemical transformations by using flow-based generative models. The method, nonadiabatic force matching, hinges on estimating the dissipation along an alchemical switching process in terms of a nonadiabatic force field that can be learned through stochastic flow matching. The learned field can be used in conjunction with short-time trajectory data to evaluate upper and lower bounds on the alchemical free energy that variationally converge to the exact value if the field is optimal. Applying the method to evaluate the alchemical free energy of atomistic models shows that it can substantially reduce the simulation cost of a free-energy estimate at a negligible loss of accuracy when compared with thermodynamic integration.

Computational chemistry↗
DOE OSTI2025

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Bubble Transport through a Porous Lattice with an Applied Inlet Flow

Within gas-evolving electrochemical systems, bubbles negatively impact performance by covering electrode active sites for reactions, blocking electric field lines, and obstructing liquid electrolyte flow causing pressure buildup. Recent additive manufacturing advances have enabled tuned porous electrode microstructures to be created, but producing systems that maximize electrochemical throughput and minimize bubble impact remains challenging. Thus, improved physical understanding of and modeling capabilities for bubble behavior are critical to improve electrolyzer design. To address this need, this study examines rising stage bubbles within a lattice with an applied liquid flow—an underexplored regime that strongly influences an electrochemical bubble’s fate. Notably, theoretical predictions and resolved bubble simulations are complemented by experiments from a 3D-printed visualization cell that matches the simulation geometry. The minimum threshold flow rate to achieve bubble breakthrough is found to be larger for higher porosities and for smaller bubbles. Different-sized bubbles decrease expected electrochemical performance in different ways; smaller bubbles tend to stay stuck but cover less solid surface, while larger bubbles more readily break through but cover more surface while in the lattice. The bubble trajectory, deformation, and contact area provide insight into these different behaviors. These findings provide design guidelines toward creating more effective electrolyzers.

Guo, Jack [Lawrence Livermore National Laboratory ↗
DOE OSTI2025

Competition for Management and Operation of Los Alamos National Laboratory

Reconnection is one of the key nonlinear processes in a plasma, in which great amounts of magnetic energy can be rapidly transformed and released via a change in magnetic topology. Field lines of opposite polarity flow together, annihilate, and release pent-up magnetic stress. Snapshots from one of our 2D test simulations (Fig 1) demonstrate key stages in magnetic reconnection. First, antiparallel magnetic field lines are brought together by an inflow in the z direction (the “upstream”). Second, the field lines begin reconnecting and an outflow forms in the x direction (the “downstream”). Finally, in the fully nonlinear stage, structures known as plasmoids are formed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Analyticity and the Unruh effect: a study of local modular flow

The Unruh effect can be formulated as the statement that the Minkowski vacuum in a Rindler wedge has a boost as its modular flow. In recent years, other examples of states with geometrically local modular flow have played important roles in understanding energy and entropy in quantum field theory and quantum gravity. Here I initiate a general study of the settings in which geometric modular flow can arise, showing (i) that any geometric modular flow must be a conformal symmetry of the background spacetime, and (ii) that in a well behaved class of “weakly analytic” states, geometric modular flow must be future-directed. I further argue that if a geometric transformation is conformal but not isometric, then it can only be realized as modular flow in a conformal field theory. Finally, I discuss a few settings in which converse results can be shown — i.e., settings in which a state can be constructed whose modular flow reproduces a given vector field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Effect of insulator end cap thickness on time-dependent Hartmann flow in a rotating mirror

We present a framework for analysing plasma flow in a rotating mirror. By making a series of physical assumptions, we reduce the magnetohydrodynamic (MHD) equations in a three-dimensional cylindrical system to a one-dimensional system in a shallow, cuboidal channel within a transverse magnetic field, similar to the Hartmann flow in ducts. We then solve the system both numerically and analytically for a range of values of the Hartmann number and calculate the dependence of the plasma flow speed on the thickness of the insulating end cap. We observe that the mean flow overshoots and decelerates before achieving a steady-state value, a phenomenon that the analytical model cannot capture. This overshoot is directly proportional to the thickness of the insulating end cap and the external electric field, with a weak dependence on the external magnetic field. Our simplified model can act as a benchmark for future simulations of the supersonic mirror device CMFX (centrifugal magnetic fusion experiment), which will employ more sophisticated physics and realistic magnetic field geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evidence of a toroidal magnetic field in the core of 3C 84

The spatial scales of relativistic radio jets, probed by relativistic magneto-hydrodynamic (RMHD) jet launching simulations and by most very long baseline interferometry (VLBI) observations differ by an order of magnitude. Bridging the gap between these RMHD simulations and VLBI observations requires selecting nearby active galactic nuclei (AGN), the parsec-scale region of which can be resolved. The radio source 3C 84 is a nearby bright AGN fulfilling the necessary requirements: it is launching a powerful, relativistic jet powered by a central supermassive black hole, while also being very bright. Using 22 GHz globe-spanning VLBI measurements of 3C 84 we studied its sub-parsec region in both total intensity and linear polarisation to explore the properties of this jet, with a linear resolution of ~0.1 parsec. We tested different simulation set-ups by altering the bulk Lorentz factor Γ of the jet, as well as the magnetic field configuration (toroidal, poloidal, helical). We confirm the persistence of a limb brightened structure, which reaches deep into the sub-parsec region. The corresponding electric vector position angles (EVPAs) follow the bulk jet flow inside but tend to be orthogonal to it near the edges. Our state-of-the-art RMHD simulations show that this geometry is consistent with a spine-sheath model, associated with a mildly relativistic flow and a toroidal magnetic field configuration.

3C 84 (NGC 1275)↗
DOE OSTI2024

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗