Search NASA⌕ Search

SEARCH · Search NASA

Results for “STEADY FLOW”

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

A Mixed Length Scale Model for Migrating Fluvial Bedforms

With the expansion of hydropower, in-stream converters, flood-protection infrastructures, and growing concerns on deltas fragile ecosystems, there is a pressing need to evaluate and monitor bedform sediment mass flux. It is critical to estimate real-time bedform size and migration velocity and provide a theoretical framework to convert easily accessible time histories of bed elevations into spatially evolving patterns. In this study, we collected spatiotemporally resolved bathymetries from laboratory flumes and the Colorado River in statistically steady, homogeneous, subcritical flow conditions. Wave number and frequency spectra of bed elevations show compelling evidence of scale-dependent velocity for the hierarchy of migrating bedforms observed in the laboratory and field. New scaling laws were applied to describe the full range of migration velocities as function of two dimensionless groups based on the bed shear velocity, sediment diameter, and water depth. Further simplification resulted in a mixed length scale model estimating scale-dependent migration velocities, without requiring bedform classification or identification.

58 GEOSCIENCES↗
DOE OSTI2022

Physics‐Informed Machine Learning Method for Large‐Scale Data Assimilation Problems

Abstract We develop a physics‐informed machine learning approach for large‐scale data assimilation and parameter estimation and apply it for estimating transmissivity and hydraulic head in the two‐dimensional steady‐state subsurface flow model of the Hanford Site given synthetic measurements of said variables. In our approach, we extend the physics‐informed conditional Karhunen‐Loéve expansion (PICKLE) method to modeling subsurface flow with unknown flux (Neumann) and varying head (time‐dependent Dirichlet) boundary conditions. We demonstrate that the PICKLE method is comparable in accuracy with the standard maximum a posteriori (MAP) method, but is significantly faster than MAP for large‐scale problems. Both methods use a mesh to discretize the computational domain. In MAP, the parameters and states are discretized on the mesh; therefore, the size of the MAP parameter estimation problem directly depends on the mesh size. In PICKLE, the mesh is used to evaluate the residuals of the governing equation, while the parameters and states are approximated by the truncated conditional Karhunen‐Loéve expansions with the number of parameters controlled by the smoothness of the parameter and state fields, and not by the mesh size. For a considered example, we demonstrate that the computational cost of PICKLE increases near linearly (as N 1.15 ) with the number of grid nodes N , while that of MAP increases much faster (as N 3.28 ). We also show that once trained for one set of Dirichlet boundary conditions (i.e., one river stage), the PICKLE method provides accurate estimates of the hydraulic head for any value of the Dirichlet boundary conditions (i.e., for any river stage).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Relating Hydro–Mechanical and Elastodynamic Properties of Dynamically Stressed Tensile–Fractured Rock in Relation to Applied Normal Stress, Fracture Aperture, and Contact Area

We exploit nonlinear elastodynamic properties of fractured rock to probe the micro-scale mechanics of fractures and understand the relation between fluid transport and fracture aperture under dynamic stressing. Experiments were conducted on rough, tensile-fractured Westerly granite subject to triaxial stresses. We measure fracture permeability for steady-state fluid flow with deionized water. Pore pressure oscillations are applied at amplitudes ranging from 0.2 to 1 MPa at 1 Hz frequency. During dynamic stressing we transmit ultrasonic signals through the fracture using an array of piezoelectric transducers (PZTs) to monitor evolution of interface properties. We examine the influence of fracture aperture and contact area by conducting measurements at effective normal stresses of 10–20 MPa. Additionally, the evolution of contact area with stress is characterized using pressure sensitive film. These experiments are conducted separately with the same fracture and map contact area at stresses from 9 to 21 MPa. The measurements are a proxy for “true” contact area for the fracture surface and we relate them to elastic properties using the calculated PZT sensor footprints via numerical modeling of Fresnel zones. We compare the elastodynamic response of the fracture using the stress-induced changes in ultrasonic wave velocities for transmitter-receiver pairs to image spatial variations in contact properties. We show that nonlinear elasticity and permeability enhancement decrease with increasing normal stress. Additionally, post-oscillation wave velocity and permeability exhibit quick recoveries toward pre-oscillation values. Estimates of fracture contact area (global and local) demonstrate that the elastodynamic and permeability responses are dominated by fracture topology.

58 GEOSCIENCES↗
DOE OSTI2021

The Effect of Compressibility and Outer Boundaries on Incipient Viscous Fingering

The knowledge of the effects of instability and heterogeneity on displacements, primarily enhanced oil recovery, and carbon dioxide storage are well known, although they remain difficult to predict. Here, the usual recourse to modeling these effects is through numerical simulation. Simulation remains the gold standard for prediction; however, its results lack generality, being case-specific. There are also several analytic models for displacements that are usually more informative than simulation results. However, these methods apply to steady-state, incompressible flow.

02 PETROLEUM↗
DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Salinity Contributions from Geothermal Waters to the Rio Grande and Shallow Aquifer System in the Transboundary Mesilla (United States)/Conejos-Médanos (Mexico) Basin

Freshwater scarcity has raised concerns about the long-term availability of the water supplies within the transboundary Mesilla (United States)/Conejos-Médanos (Mexico) Basin in Texas, New Mexico, and Chihuahua. Analysis of legacy temperature data and groundwater flux estimates indicates that the region’s known geothermal systems may contribute more than 45,000 tons of dissolved solids per year to the shallow aquifer system, with around 8500 tons of dissolved solids being delivered from localized groundwater upflow zones within those geothermal systems. If this salinity flux is steady and eventually flows into the Rio Grande, it could account for 22% of the typical average annual cumulative Rio Grande salinity that leaves the basin each year—this salinity proportion could be much greater in times of low streamflow. Regional water level mapping indicates upwelling brackish waters flow towards the Rio Grande and the southern part of the Mesilla portion of the basin with some water intercepted by wells in Las Cruces and northern Chihuahua. Upwelling waters ascend from depths greater than 1 km with focused flow along fault zones, uplifted bedrock, and/or fractured igneous intrusions. Overall, this work demonstrates the utility of using heat as a groundwater tracer to identify salinity sources and further informs stakeholders on the presence of several brackish upflow zones that could notably degrade the quality of international water supplies in this developed drought-stricken region.

Pepin, Jeff D. (ORCID:0000000274109979)↗
DOE OSTI2026

Theoretical background for a fast flow liquid metal divertor experiment

Tokamak Energy has developed a liquid metal experiment featuring a lithium loop that circulates lithium through the HIDRA stellarator at the University of Illinois at Urbana-Champaign. The loop includes a replaceable divertor module, designed to demonstrate fast, steady, open-surface liquid metal flow. The first-generation divertor module was deliberately kept simple—an open-surface chute—to enable modelling validation and establish a benchmark for future design iterations. This paper presents the theoretical foundations of the experimental design. Specifically, we identify module overflooding as the primary challenge and determine the limits of fluid velocity and magnetic field strength necessary to prevent it. The steady-state flow patterns are expected to exhibit relatively slow surface velocities, which do not fully align with the fast-flow concept. Nevertheless, the experiment is designed to achieve controllable, continuous open-surface flow within an operational fusion device, representing a significant step toward the development of liquid metal divertor technology.

Experiment↗
DOE OSTI2023

Advanced two-phase subchannel method via non-linear iteration

A fast-running, robust two-phase flow, sub-channel model is presented based on non-linear solution of the steady-state subchannel fluid flow equations. The drift-flux model solves for conservation of liquid and vapor mass, mixture energy, and axial and transverse mixture momentum as part of an efficient planar marching scheme and nonlinear, nested outer and inner iteration. Here, models based on mechanistic subcooled boiling, two-phase turbulent void mixing, and drift are included. Solution verification and mesh convergence studies were performed for modern GE 10 × 10 fuel geometry and are shown to have excellent convergence behavior. Run time performance for a 50 axial mesh model showed 2.2 seconds on a single CPU core to tightly converge all 3D distributions (flow, void, pressure) for the GE 10 × 10 fuel geometry, supporting its efficient use within the Virtual Environment for Reactor Applications boiling water reactor framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

TRACE Analysis of the STS Target Flow Loop Conceptual Design

This report documents the efforts to model the conceptual design of the STS natural circulation flow loop using TRACE. The results, summarized in Section 2.1, support the conclusion that the conceptual STS natural circulation flow loop design provides sufficient cooling of the tungsten target and exhibits predictable flow behavior, provided that certain design features and/or operating conditions are used. In particular: A sufficient initial liquid volume fraction is required to avoid dryout of the target flow channels; Inclusion of a pressurizer or bellows volume at the top of the loop is recommended and is capable of achieving steady, non-oscillatory flow and temperature conditions at sufficiently high loop pressure. Section 1 provides a detailed description of the geometry and assumptions used in the model, and Section 2 presents the results for the current flow loop design as well as some investigation of potential design modifications with the goal of enhancing the thermal hydraulic stability of the loop.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Micro-CT Imaging and Fluid Flow Simulations of Fractures in MSEEL Shale

micro-CT scans of a naturally fractured MSEEL shale sample as a shear fracture is generated and displaced in the center, intact region of the sample. Steady State laminar fluid flow simulations were then performed on each individual fracture during each experimental segment as effective stress and shear displacement changed. Released under LA-UR-21-24335.

Fractures,Fracturing Initiation,Marcellus Shale↗
DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental characterization of extreme temperature granular flows for solar thermal energy transport and storage

High-temperature, dense granular flows along an inclined plane were considered for solar thermal energy transport and storage with sintered bauxite particles. A series of experiments was performed for particle inlet temperatures of ~ 200, 400, 600, and 800 °C to understand the mechanisms of granular flows at extreme temperatures. Mass flow rates were measured using a load cell and free-surface velocities were measured and computed using particle image velocimetry. Surface temperatures were measured using infrared cameras. A significant decrease in steady-state particle mass flow rate was observed with increasing temperature due to changing flow properties. A decrease in bulk particle free-surface velocities was observed at higher temperatures. Free-surface velocity measurement error between experiments were within 20% of the average. The particle surface temperatures decreased from inlet to outlet with larger gradients at higher temperatures observed due to increasing convection and radiative heat losses. Here, a decrease in temperature was observed along the side walls due to a decrease in particle velocities.

14 SOLAR ENERGY↗
DOE OSTI2023

Geometric Optimization of an Electrochemical Purification Cell to Prevent Corrosion in CSP Plants During Operation

When exposed to moisture or oxygen, molten chloride salts produce corrosive impurities which degrade containment alloys. This can significantly decrease the lifetime and increase costs of molten-salt-based systems. To overcome this barrier, we designed and modeled an electrochemical purification cell to remove the corrosive impurity MgOH+. Various reactor architectures, including continuous stirred tank reactors (CSTRs) and plug flow reactors (PFRs) were investigated. Steady-state thermoelectric properties were evaluated using analytical methods, allowing assessment of the effects of structure and design parameters such as flow rate, cell length, and cross-sectional area of molten salt. The results suggest that our design could most effectively increase reliability and decrease costs of molten-chloride-salt-based systems by protecting them during continuous operation using an annular plug flow reactor.

analytical modeling↗
DOE OSTI2021

Time-periodic steady-state solution of fluid-structure interaction and cardiac flow problems through multigrid-reduction-in-time

In this study, a time-periodic MGRIT algorithm is proposed as a means to reduce the time-to-solution of numerical algorithms by exploiting the time periodicity inherent to many applications in science and engineering. The time-periodic MGRIT algorithm is applied to a variety of linear and nonlinear single- and multiphysics problems that are periodic-in-time. It is demonstrated that the proposed parallel-in-time algorithm can obtain the same time-periodic steady-state solution as sequential time-stepping. It is shown that the required number of MGRIT iterations can be estimated a priori and that the new MGRIT variant can significantly and consistently reduce the time-to-solution compared to sequential time-stepping, irrespective of the number of dimensions, linear or nonlinear PDE models, single-physics or coupled problems and the employed computing resources. The numerical experiments demonstrate that the time-periodic MGRIT algorithm enables a greater level of parallelism yielding faster turnaround, and thus, facilitating more complex and more realistic problems to be solved.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Behavior of an air-cooled thermosiphon during a non-air gas ingress event

Here, this study investigates the transient response of an air-cooled thermosiphon subjected to a sudden non-air gas ingress event. Results presented are from experimental testing conducted on a large-scale thermal hydraulic facility that examines heat removal performance of the Reactor Cavity Cooling System (RCCS) concept, a passive safety system for advanced nuclear reactors. Relying on buoyancy driven natural circulation flow, these systems feature a series of riser standpipes and coupled chimney network that provide a pathway for heat rejection to the atmosphere. In this work, the facility was modified to allow a transition of the draft intake from normal ambient air, to draft intake from an open-top volume containing high purity argon. Testing began by establishing normal operation and allowing the facility to reach steady-state thermal hydraulic flow conditions with natural air draft. The transitioned sequence was then initiated which simultaneously changed the inlet boundary from open air to a pathway from the quiescent argon volume. This event caused 1,200 cu-feet of the heavy gas, twice the internal volume of total facility flow path, to be drawn into the inlet plenum and ingress into the heated riser standpipes. With the upper chimney network still containing residual air, there was insufficient density difference to maintain buoyancy driven natural circulation. After a period of only 90 s past the ingress event, the facility experienced complete flow stagnation before entering an extended period of severely degraded system flow. Due to the cessation of bulk fluid movement and subsequent failure of its heat removal function, fluid and structural temperatures began to rise sharply. Re-circulation patterns developed within the multiple parallel riser standpipes, where hot gas near the outlets was observed to travel downward and re-enter the inlet of adjacent channels. After approximately 18-minutes, fluid temperatures and their associated density difference rose to a level sufficient to allow reestablishment of buoyancy driven system flow, and ultimately, recovery of facility operation to normal behavior.

42 ENGINEERING↗
DOE OSTI2025

Liquid crystal torons in Poiseuille-like flows

Three-dimensional (3D) simulations of the structure of liquid crystal (LC) torons, topologically protected distortions of the LC director field, under material flows are rare but essential in microfluidic applications. Here, we show that torons adopt a steady-state configuration at low flow velocity before disintegrating at higher velocities, in line with experimental results. Furthermore, we show that under partial slip conditions at the boundaries, the flow induces a reversible elongation of the torons, also consistent with the experimental observations. These results are in contrast with previous simulation results for 2D skyrmions under similar flow conditions, highlighting the need for a 3D description of this LC soliton in relation to its coupling to the material flow. These findings pave the way for future studies of other topological solitons, like hopfions and heliknotons, in flowing soft matter systems.

Science & Technology - Other Topics↗
DOE OSTI2023

Extending the single-fluid solvability conditions for more general plasma systems

We extend the single-fluid solvability conditions to include plasma systems at arbitrary $\beta$ with arbitrary flows and external forcing terms. This treatment includes both the isotropic and the anisotropic cases. The generalized conditions that result can be used to generate certain classes of single-fluid equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Evolution of eddy viscosity in the wake of a wind turbine

The eddy viscosity hypothesis is a popular method in wind turbine wake modeling for estimating turbulent Reynolds stresses. We document the downstream evolution of eddy viscosity in the wake of a wind turbine from experimental and large-eddy-simulation data. Wake eddy viscosity is isolated from its surroundings by subtracting the inflow profile, and the driving forces are identified in each wake region. Eddy viscosity varies in response to changes in turbine geometry and nacelle misalignment with larger turbines generating stronger velocity gradients and shear stresses. We propose a model for eddy viscosity based on a Rayleigh distribution. Model parameters are obtained from scaling the eddy viscosity hypothesis and demonstrate satisfactory agreement with the reference data. The model is implemented in the curled wake formulation in the FLOw Redirection and Induction in Steady State (FLORIS) framework and assessed through comparisons with the previous formulation. Our approach produced more accurate flow field estimates with lower total error for the majority of cases.

17 WIND ENERGY↗