SEARCH · Engineering Papers

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Digital Valve for High Pressure High Flow Applications

To address the challenges, which are involved with the development of flow control valves that can meet the requirements of deep oil wells such as high pressure, high flow rate, limited power and limited space, the authors have conceived a novel design configuration [1]. This design consists of a digitalized flow control valve with multipath and multistage pressure reduction structures. Specifically, the valve is configured as a set of parallel flow paths from the inlet to the outlet. A choke valve controls the total flow rate by digitally opening different paths or different combination of the paths. Each path is controlled by a poppet cap valve basically operated in on-off states. The number of flow states is 2N where N is the number of flow paths. To avoid erosion from sand in the oil and high speed flow, the seal area of the poppet cap valve is located at a distance from the flow inlet away from the high speed flow and the speed is controlled to stay below a predefined erosion safe limit. The path is a multistage structure composed of a set of serial nozzles-expansion chambers that equally distribute the total pressure drop to each stage. The pressure drop of each stage and, therefore, the flow speed at the nozzles and expansion chambers is controlled by the number of stages. The paths have relatively small cross section and could be relatively long for large number of stages and still fit in the strict annular space limit of a typical downhole region of an oil well. The paper will present the design configuration, analysis and preliminary test results.

downhole↗

DOE OSTI2017

FLOW TESTING AND ANALYSIS OF THE FSP-1 EXPERIMENT

The U.S. High Performance Research Reactor Conversions fuel development team is focused on developing and qualifying the uranium-molybdenum (U-Mo) alloy monolithic fuel to support conversion of domestic research reactors to low enriched uranium. Several previous irradiations have demonstrated the favorable behavior of the monolithic fuel. The Full Scale Plate 1 (FSP-1) fuel plate experiment will be irradiated in the northeast (NE) flux trap of the Advanced Test Reactor (ATR). This fueled experiment contains six aluminum-clad fuel plates consisting of monolithic U-Mo fuel meat. Flow testing experimentation and hydraulic analysis have been performed on the FSP-1 experiment to be irradiated in the ATR at the Idaho National Laboratory (INL). A flow test experiment mockup of the FSP-1 experiment was completed at Oregon State University. Results of several flow test experiments are compared with analyses. This paper reports and shows hydraulic analyses are nearly identical to the flow test results. A water velocity of 14.0 meters per second is targeted between the fuel plates. Comparisons between FSP-1 measurements and this target will be discussed. This flow rate dominates the flow characteristics of the experiment and model. Separate branch flows have minimal effect on the overall experiment. A square flow orifice was placed to control the flowrate through the experiment. Four different orifices were tested. A flow versus delta P curve for each orifice is reported herein. Fuel plates with depleted uranium in the fuel meat zone were used in one of the flow tests. This test was performed to evaluate flow test vibration with actual fuel meat densities and reported herein. Fuel plate deformation tests were also performed and reported.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2016

Characterization of Impedance Void Meter Performance in Two-Phase Flow Measurement

Over the past sixty years, exclusive work has been done on two-phase flow measurement and plenty of measurement methods have been developed. Among them, impedance void meter is considered to be a simple instrument with a quick and reliable response. Impedance measurement technique in two-phase flow has existed for more than five decades. Impedance void meter consists of two electrodes flush mounted on the flow channel which allows measurements to be performed over the entire flow area. In principle, impedance void meter is based on the different electrical conductance of water and air. The two electrodes connected with cable and voltage supplies and the air-water mixture between them form a complete circuit. The change of gas fraction occupied in the flow area is reflected by the change of impedance in this part, and then the change of the voltage signal of the circuit. Using the relationship between void fraction and voltage signal, area averaged void fraction can be measured. In a previous study by Bolman et al., the flow pattern was found that has an effect on the performance of the impedance void meter. This is due to the difference of the flow structure that affects the impedance of the flow area. Olsen suggested that the impedance void meter should fulfill the following factors in order to get good measurement results: gas phase is finely dispersed in the measuring flow area; the channel cross section should be uninterrupted; there should be several bubble diameters between sources of current. For pipe measurement, it is impossible to meet all these requirement. In this case, it is necessary to use other independent void fraction measurement methods to calibrate impedance void meter. However, the calibration result varies for different flow patterns; thus for those situations in which identifying the flow pattern is difficult, the results may not be reliable. The possible solution should be found in an aspect of impedance void meter design itself. Much work has been done focusing on the geometric design of impedance void meters. This study found that the distance between the two electrodes of the impedance void meter has an effect on the measurement performance. The primary focus of this study is to characterize the performance and measurement range of an impedance void meter in regards to its geometry aspect. In order to achieve this goal, a sensitivity study of the electrode distance was performed by numerical simulation. Based on the simulation results, an experiment was performed using an impedance void meter with different electrode distances. Additionally, a gamma densitometer and a multi single-sensor conductivity probe were also applied for the measurement comparison. Drift flux model was also used for the experimental results comparison. (authors)

42 ENGINEERING↗

DOE OSTI2016

Characterizing dynamic hysteresis and fractal statistics of chaotic two-phase flow and application to fuel cells

In this study, we analyze the stability of two-phase flow regimes and their transitions using chaotic and fractal statistics, and we report new measurements of dynamic two-phase pressure drop hysteresis that is related to flow regime stability and channel water content. Two-phase flow dynamics are relevant to a variety of real-world systems, and quantifying transient two-phase flow phenomena is important for efficient design. We recorded two-phase (air and water) pressure drops and flow images in a microchannel under both steady and transient conditions. Using Lyapunov exponents and Hurst exponents to characterize the steady-state pressure fluctuations, we develop a new, measurable regime identification criteria based on the dynamic stability of the two-phase pressure signal. We also applied a new experimental technique by continuously cycling the air flow rate to study dynamic hysteresis in two-phase pressure drops, which is separate from steady-state hysteresis and can be used to understand two-phase flow development time scales. Using recorded images of the two-phase flow, we show that the capacitive dynamic hysteresis is related to channel water content and flow regime stability. The mixed-wettability microchannel and in-channel water introduction used in this study simulate a polymer electrolyte fuel cell cathode air flow channel.

42 ENGINEERING↗

DOE OSTI2021

Turbulence-driven flow dynamics in general axisymmetric toroidal geometry

The present work gives the equations governing the generation of toroidally axisymmetric flows by turbulent Reynolds and Maxwell stresses in finite aspect ratio, general cross section tokamak plasmas. Inclusion of the divergence-free flow constraint in lowest order changes the time scale for evolution of the poloidal flows driven by turbulence by substantial factors. In the pedestal region for present day machines, comparing to earlier cylindrical models, the time scale evaluated using a large aspect ratio circular cross section model can be two orders of magnitude longer while the present, general geometry result can be about one order of magnitude longer. Inclusion of gyroviscosity in the calculation shows that the only lowest order radial velocity fluctuations that enter the problem are those due to fluctuating E×B flows. Toroidal and poloidal flow effects on the toroidally axisymmetric flows are inextricably coupled due to the neoclassical poloidal viscosity. Accordingly, the physics is inherently three dimensional and measurements of all three velocity components are required to obtain the information needed to quantitatively test the theory. The parallel and angular momentum equations for the lowest order, toroidally axisymmetric flows look like radial transport equations when the turbulence is included. The turbulence terms provide the radial transport fluxes. In standard neoclassical theory, the parallel flow equation is local on each flux surface; there is no radial derivative term. However, adding turbulence gives a way, in principle, for radial transport to lead to poloidal flows that deviate from the neoclassical prediction. As a result, the inclusion of the Maxwell stress provides a mechanism for MHD fluctuations to alter the toroidally axisymmetric flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

On the statistical theory of self-gravitating collisionless dark matter flow

Dark matter, if it exists, accounts for five times as much as the ordinary baryonic matter. Compared to hydrodynamic turbulence, the flow of dark matter might possess the widest presence in our universe. This paper presents a statistical theory for the flow of dark matter that is compared with N-body simulations. By contrast to hydrodynamics of normal fluids, dark matter flow is self-gravitating, long-range, and collisionless with a scale-dependent flow behavior. The peculiar velocity field is of constant divergence nature on small scale and irrotational on large scale. The statistical measures, i.e., correlation, structure, dispersion, and spectrum functions, are modeled on both small and large scales, respectively. Kinematic relations between statistical measures are fully developed for incompressible, constant divergence, and irrotational flow. Incompressible and constant divergence flow share the same kinematic relations for even order correlations. The limiting correlation of velocity $\mathrm{ρ_{L}=1/2}$ on the smallest scale ( r = 0) is a unique feature of collisionless flow (⁠$\mathrm{ρ_{L}=1}$ for incompressible flow). On large scale, transverse velocity correlation has an exponential form $T_{2}∝e^{–r/r_2}$ with a constant comoving scale r 2 =21.3 Mpc/h that may be related to the horizon size at matter–radiation equality. All other correlation, structure, dispersion, and spectrum functions for velocity, density, and potential fields are derived analytically from kinematic relations for irrotational flow. On small scale, longitudinal structure function follows one-fourth law of ${S}_{2}^{1}∝r^{1/4}$. All other statistical measures can be obtained from kinematic relations for constant divergence flow. Vorticity is negatively correlated for scale r between 1 and 7 Mpc/h. Divergence is negatively correlated for r > 30 Mpc/h that leads to a negative density correlation.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1997

Experimental Investigation of the Near Wall Flow Structure of a Low Reynolds Number 3-D Turbulent Boundary Layer

Laser Doppler velocimetry (LDV) measurements and hydrogen bubble flow visualization techniques were used to examine the near-wall flow structure of 2D and 3D turbulent boundary layers (TBLs) over a range of low Reynolds numbers. The goals of this research were (1) an increased understanding of the flow physics in the near wall region of turbulent boundary layers,(2) to observe and quantify differences between 2D and 3D TBL flow structures, and (3) to document Reynolds number effects for 3D TBLs. The LDV data have provided results detailing the turbulence structure of the 2D and 3D TBLs. These results include mean Reynolds stress distributions, flow skewing results, and U and V spectra. Effects of Reynolds number for the 3D flow were also examined. Comparison to results with the same 3D flow geometry but at a significantly higher Reynolds number provided unique insight into the structure of 3D TBLs. While the 3D mean and fluctuating velocities were found to be highly dependent on Reynolds number, a previously defined shear stress parameter was discovered to be invariant with Reynolds number. The hydrogen bubble technique was used as a flow visualization tool to examine the near-wall flow structure of 2D and 3D TBLs. Both the quantitative and qualitative results displayed larger turbulent fluctuations with more highly concentrated vorticity regions for the 2D flow.

FLOW VISUALIZATION↗

NASA NTRS2022

Design of Hypervelocity Flow Generator (HFG) and Its Diagnostics

Ground facilities for hypersonic research are a key element for successful development of aerodynamically proven hypervelocity vehicles. Design concepts and diagnostics of a hypervelocity flow generator (HFG) were made as a test platform for hypersonic/hypervelocity spacecraft models at the NASA Langley Research Center. The HFG is a hypersonic flow field generator using optically heated gas which is blown into an 80 cu m vacuum chamber. The vacuum chamber is kept at a stable vacuum pressure with a combination of three large vacuum pumps, while the HFG is in the test mode. The HFG provides a relatively small test section with approximately a 20 cm window. This facility was designed to generate 2.45 km/sec of flow speed, and potentially generate a continuous flow with the nozzle and vacuum system. The window in test section provides a direct view of the shock wave around a model in order to measure temperature, pressure, and density profile within the shock layer. One of the key test goals under this project is to understand why the emission spectra from the standing shockwave plasma predicted by Lora-Loran codes are significantly different from the measured emission spectra from the Flight Investigation Reentry Environment (FIRE II) Flight. The correct estimate of the thermal loading on the leading edge of hypersonic vehicles greatly affects the aerodynamic design, the material selection for the vehicle, and the cooling requirement and can be obtained by the precise modeling of emission spectra from shockwave plasma. However, the estimation of thermal loading is not an easy task due to complex non-equilibrium radiative process within high temperature shock layers that still falls into a category of cold plasma. Direct flight experiments are the most desirable, but not a cost-effective approach. Analysis by computational fluid dynamics (CFD) offers many test flexibilities. However, the CFD codes must be fully tested and validated with experimental data before the codes are effectively used for practical design. Large discrepancies between experiments and codes appear in hypersonic/hypervelocity flow regimes at high altitudes of 60 km ~ 90 km. This HFG facility offers some important parameters for CFD code validation, such as collision cross-sections, relaxation times, reaction rate coefficients and transportation coefficients. The HFG test facility is based on the ejection flow of high temperature gas heated over to 3500 K through a nozzle. The tungsten gas chamber of the HFG is heated up to a desired temperature by a 60 kW optical power beam source. This system consists of an optical power source, a thermal chamber, an expansion nozzle, a test section, and an 80 cu m vacuum tank. 60 kW optical input power is obtained from the 150 kW Vortek arc lamp system (by Vortek Industries, Vancouver, Canada). This optical beam is focused to heat the gas chamber within which a flow media is heated. The maximum achievable temperature of the flow medium reached approximately 3500 K or even higher but is limited by the melting point of the chamber material used. The exhaust velocity through the nozzle was determined by the stagnation temperature and the molecular weight of the working medium at the test section. To provide design parameters, a NASA Chemical Equilibrium with Application (CEA) computer program is used for the simulation of aerothermal data. This CEA program can calculate chemical equilibrium and properties of complex mixtures using shock tube parameters. For nitrogen gas at 2666 K stagnation temperature, the maximum achievable velocity at the test section is approximately 2.45 km/sec which is within the range of the thermal velocity of 8000 K shock layer. Based on the calculation through the CEA program, the design parameters of a HFG were determined and implemented for the test section that includes an expansion nozzle. The installed test section of the experimental facility can sustain a condition of a re-entry vehicle from the Space at an altitude of 60-90 km. After installation of the HFG, the system was fully tested and its operational parameters were measured. An 80 cu m vacuum chamber of HFG was set at 1 torr level to keep a stable downstream condition. The pumping time to reach the minimum vacuum pressure (~ 1 torr) at the test chamber from atmospheric pressure was approximately an hour using Kinney (MBV-14000/MB -1600/KT-300) and Stokes (Model 1772 and 412) pumps. Such a setup condition allowed a continued stable operation of the HFG experiment with flowrates through 1-, 3-, or 5-mm diameters of nozzle throat. The flow characteristics of the HFG for various operating conditions were performed using a focused Schlieren method. At 0.8 torr chamber pressure, a barrel shock was observed at the test section. With a sphere obstacle of 12 mm diameter in the flow, a bow shock (~ 2-mm thickness) was observed by a focused Schlieren visualization method.

Hypervelocity flow↗

DOE OSTI2022

Modeling the role of lateral surface flow in low‐relief polygonal tundra

Abstract Ice‐wedge polygon troughs play an important role in controlling the hydrology of low‐relief polygonal tundra regions. Lateral surface flow is confined to troughs only, but it is often neglected in model projections of permafrost thermal hydrology. Recent field and modeling studies have shown that, after rain events, increases in trough water levels are significantly more than the observed precipitation, highlighting the role of lateral surface flow in the polygonal tundra. Therefore, understanding how trough lateral surface flow can influence polygonal tundra thermal hydrology is important, especially under projected changes in temperatures and rainfall in the Arctic regions. Using an integrated cryohydrology model, this study presents plot‐scale end‐of‐century projections of ice‐wedge polygon water budget components and active layer thickness with and without trough lateral surface flow under the Representative Concentration Pathway 8.5 scenario. Trough lateral surface flow is incorporated through a newly developed empirical model, evaluated against field measurements. The numerical scenario that includes trough lateral surface flow simulates discharge (outflow from a polygon) and recharge (rain‐induced inflow to a polygon trough from upslope areas), while the scenario that does not include trough lateral surface flow ignores recharge. The results show considerable reduction (about 100–150%) in evapotranspiration and discharge in rainy years in the scenarios ignoring trough lateral surface flow, but less effect on soil water storage, in comparison with the scenario with trough lateral surface flow. In addition, the results demonstrate long‐term changes (~10–15 cm increase) in active layer thickness when trough lateral surface flow is modeled. This study highlights the importance of including lateral surface flow processes to better understand the long‐term thermal and hydrological changes in low‐relief polygonal tundra regions under a changing climate.

Jan, Ahmad↗

DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part I

In this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In conclusion, these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II

Here, in this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Shale-fluid interactions during forced lmbibition and flow-back

One way to minimize the need for disposal of flow-back waters produced during the hydraulic fracturing of unconventional reservoirs is, potentially, to recycle these fluids by recovering and re-injecting them during hydraulic fracturing of nearby new wells. These flow-back waters contain, however, large concentrations of dissolved salts. During hydraulic fracturing of unconventional reservoirs, a well may undergo extended shut-in periods, following the initial hydraulic-fracture stimulation. The concern here, therefore, is that the salts dissolved in the re-injected flow-back waters may precipitate in the subsurface, thus promoting scale formation; or alternately, because of their high salinity, these same fluids may cause the dissolution of various minerals present in the subsurface. Large divalent cations (e.g., Ca²⁺, Ba²⁺ and Sr²⁺) and anions (e.g., CO3²⁻ and SO4²) may chemically interact with and cause the swelling of clay fractions, thus blocking pores in certain pore size ranges. All these phenomena (i.e., scale formation, clay swelling, minerals dissolution) can influence the petrophysical properties of the formation, including its permeability, porosity, contact angle, etc. One way to potentially overcome some of these technical challenges, is to pre-treat these flow-back waters, prior to re-injection (e.g., via chemical precipitation, followed by membrane filtration) in order to remove unwanted salts. However, such pre-treatment is costly, and the degree of salt removal needed to prevent formation damage is important in determining the economic viability of flow-back water recycle via re-injection. A key focus of this work, therefore, is to study the impact that the injection of flow-back waters, that have undergone varying degrees of pre-treatment and clean-up, will have on the flow-through characteristics of a tight-gas reservoir, particularly of the Marcellus Shale formation. To do so, we have carried-out forced-imbibition experiments with shale samples from various depths in the Marcellus formation, whose properties span high to low ranges of TOC, clay content and matrix permeabilities. We have studied the forced-imbibition characteristics, at an injection pressure of 24,233 kPa (3500 psig), of various flow-back fluids into shale cores saturated with methane at a pore pressure of 17,338 kPa (2500 psig) and a confining stress of 27,680 kPa (4000 psig). Following the forced imbibition experiments, we analyzed the cores for potential changes in permeability and porosity. Quantification of the ion content of the injected and recovered fluids was also performed using Ion Chromatography (IC). The IC analysis of the recovered fluids indicates substantial interactions and geochemical reactions between the invading fluids and the shale cores. The injection of flow-back fluids shows significant impact on the flow-through characteristics, with untreated flow-back fluids, which have the highest cation and anion content, also exhibiting the highest observable reductions in permeability (>90%).

defects, mechanical behavior, carbon sequestration↗

DOE OSTI2022

Understanding the role of flow dynamics in thermoacoustic combustion instability

Thermoacoustic combustion instability is one of the most challenging operational issues in several high-performance, low-emissions combustion technologies, including gas turbines, aircraft engines, rockets, and industrial boilers. Driven by the coupling between combustor acoustics and flame heat release rate fluctuations, thermoacoustic combustion instability can lead to reduced operability, increased emissions, and, in the most extreme cases, catastrophic failure of combustor components. The feedback loop between acoustics and combustion is often facilitated by fluid mechanic oscillations, referred to as “velocity coupling,” whereby acoustic oscillations drive flow fluctuations, which in turn create fluctuations in the flame. The character of these fluid mechanic oscillations is highly dependent on the structure of the flow field and the receptivity of the flow to external excitation. Combustor flow fields use features like fluid recirculation and shear to enhance flame holding and reduce emissions, but these are also the same features that can make the flow receptive to acoustic excitation or even drive self-excited oscillations. Here, in this paper, we discuss the basics of thermoacoustic instability with a focus on the role of hydrodynamic oscillations in typical combustor flows. To facilitate this discussion, we explore the hydrodynamic instability characteristics of several key combustor unit flows (wakes, swirling jets, etc.) and show how the hydrodynamic stability of a flow is an important consideration in determining a combustor’s propensity for thermoacoustic oscillations. Several examples of coupling between hydrodynamics and thermoacoustics are discussed to illustrate this important link. The paper concludes by discussing the potential for designing flow fields that are thermoacoustic instability resistant, either through a reduction in the receptivity of the flow or through nonlinear coupling mechanisms by which self-excited flow instabilities can suppress velocity-coupled combustion oscillations.

42 ENGINEERING↗

DOE OSTI2024

Heat transfer coefficients of moving particle beds from flow-dependent thermal conductivity and near-wall resistance

Accurate determination of heat transfer coefficients for flowing packed particle beds is essential to the design of particle heat exchangers and other thermal and thermochemical equipment. While such dense granular flows mostly fall into the well-known plug-flow regime, the discrete nature of granular materials alters the thermal transport processes in both the near-wall and bulk regions of flowing particle beds from their stationary counterparts. As a result, heat transfer correlations based on the stationary particle bed thermal conductivity could be inadequate for flowing particles in a heat exchanger. Most earlier works have achieved a reasonable agreement with experiments by treating granular heat transfer media as a plug-flow continuum with a near-wall thermal resistance in series. However, the thermal conductivity values of the continuum were often obtained from measurements on stationary beds owing to the difficulty of flowing bed measurements. In this work, it was found that the properties of a stationary bed are highly sensitive to the method of particle packing and there is a decrease in the particle bed thermal conductivity and increase in the near-wall thermal resistance, measured as an effective air gap thickness, on the onset of particle flow. These variations in thermal conductivity of stationary and flowing particle beds can lead to errors in heat transfer coefficient calculations. Therefore, the heat transfer coefficients for granular flows were calculated using experimentally determined flowing particle bed thermal conductivity and near-wall air gap for ceramic particles – CARBO CP 40/100 (mean diameter = 275 µm), HSP 40/70 (404 µm) and HSP 16/30 (956 µm); at velocities of 5–15 mm·s –1 ; and temperatures of 300–650 °C. The thermal conductivity and air gap values for CP 40/100 and HSP 40/70 were further used to calculate heat transfer coefficients across different particle bed temperatures and velocities for different parallel-plate heat exchanger dimensions. Furthermore, these calculations, which show good agreement with measured HTC values reported in literature, can be used as a guide for heat exchanger designs. Graphical abstract

14 SOLAR ENERGY↗

DOE OSTI2024

Integrating a Water Tracer Model Into WRF‐Hydro for Characterizing the Effect of Lateral Flow in Hydrologic Simulations

Abstract Most current land models approximate terrestrial hydrological processes as one‐dimensional vertical flow, neglecting lateral water movement from ridges to valleys. Such lateral flow is fundamental at catchment scales and becomes crucial for finer‐scale land models. To test the effect of incorporating lateral flow toward three‐dimensional representations of hydrological processes in the next generation land models, we integrate a water tracer model into the WRF‐Hydro framework to track water movement from precipitation to discharge and evapotranspiration. This hydrologic‐tracer integrated system allows us to identify the key mechanisms by which lateral flow affects the flow paths and transit times in WRF‐Hydro. By comparing modeling experiments with and without lateral routing in two contrasting catchments, we determine the impacts of lateral flow on the transit times of precipitation event‐water. Results show that with limited hydrologic connectivity, lateral flow extends the transit times by reducing (increasing) event‐water drainage loss (accumulation) in ridges (valleys) and allowing reinfiltration of infiltration‐excess flow, which is missing in most land models. On the contrary with high hydrologic connectivity, lateral flow can effectively accelerate the water release to streams and reduce the transit time. However, the transit times are substantially underestimated by the model compared with isotope‐derived estimates, indicating model limitations in representing flow paths and transit times. This study provides some insights on the fundamental differences in terrestrial hydrology simulated by land models with and without lateral flow representation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Preferential Paths of Air-water Two-phase Flow in Porous Structures with Special Consideration of Channel Thickness Effects

Accurate understanding and predicting the flow paths of immiscible two-phase flow in rocky porous structures are of critical importance for the evaluation of oil or gas recovery and prediction of rock slides caused by gas-liquid flow. A 2D phase field model was established for compressible air-water two-phase flow in heterogenous porous structures. The dynamic characteristics of air-water two-phase interface and preferential paths in porous structures were simulated. The factors affecting the path selection of two-phase flow in porous structures were analyzed. Transparent physical models of complex porous structures were prepared using 3D printing technology. Tracer dye was used to visually observe the flow characteristics and path selection in air-water two-phase displacement experiments. The experimental observations agree with the numerical results used to validate the accuracy of phase field model. The effects of channel thickness on the air-water two-phase flow behavior and paths in porous structures were also analyzed. The results indicate that thick channels can induce secondary air flow paths due to the increase in flow resistance; consequently, the flow distribution is different from that in narrow channels. This study provides a new reference for quantitatively analyzing multi-phase flow and predicting the preferential paths of immiscible fluids in porous structures.

42 ENGINEERING↗

DOE OSTI2016

Composition and Realization of Source-to-Sink High-Performance Flows: File Systems, Storage, Hosts, LAN and WAN

A number of Department of Energy (DOE) science applications, involving exascale computing systems and large experimental facilities, are expected to generate large volumes of data, in the range of petabytes to exabytes, which will be transported over wide-area networks for the purpose of storage, visualization, and analysis. To support such capabilities, significant progress has been made in various components including the deployment of 100 Gbps networks with future 1 Tbps bandwidth, increases in end-host capabilities with multiple cores and buses, capacity improvements in large disk arrays, and deployment of parallel file systems such as Lustre and GPFS. High-performance source-to-sink data flows must be composed of these component systems, which requires significant optimizations of the storage-to-host data and execution paths to match the edge and long-haul network connections. In particular, end systems are currently supported by 10-40 Gbps Network Interface Cards (NIC) and 8-32 Gbps storage Host Channel Adapters (HCAs), which carry the individual flows that collectively must reach network speeds of 100 Gbps and higher. Indeed, such data flows must be synthesized using multicore, multibus hosts connected to high-performance storage systems on one side and to the network on the other side. Current experimental results show that the constituent flows must be optimally composed and preserved from storage systems, across the hosts and the networks with minimal interference. Furthermore, such a capability must be made available transparently to the science users without placing undue demands on them to account for the details of underlying systems and networks. And, this task is expected to become even more complex in the future due to the increasing sophistication of hosts, storage systems, and networks that constitute the high-performance flows. The objectives of this proposal are to (1) develop and test the component technologies and their synthesis methods to achieve source-to-sink high-performance flows, and (2) develop tools that provide these capabilities through simple interfaces to users and applications. In terms of the former, we propose to develop (1) optimization methods that align and transition multiple storage flows to multiple network flows on multicore, multibus hosts; and (2) edge and long-haul network path realization and maintenance using advanced provisioning methods including OSCARS and OpenFlow. We also propose synthesis methods that combine these individual technologies to compose high-performance flows using a collection of constituent storage-network flows, and realize them across the storage and local network connections as well as long-haul connections. We propose to develop automated user tools that profile the hosts, storage systems, and network connections; compose the source-to-sink complex flows; and set up and maintain the needed network connections. These solutions will be tested using (1) 100 Gbps connection(s) between Oak Ridge National Laboratory (ORNL) and Argonne National Laboratory (ANL) with storage systems supported by Lustre and GPFS file systems with an asymmetric connection to University of Memphis (UM); (2) ORNL testbed with multicore and multibus hosts, switches with OpenFlow capabilities, and network emulators; and (3) 100 Gbps connections from ESnet and their Openflow testbed, and other experimental connections. This proposal brings together the expertise and facilities of the two national laboratories, ORNL and ANL, and UM. It also represents a collaboration between DOE and the Department of Defense (DOD) projects at ORNL by sharing technical expertise and personnel costs, and leveraging the existing DOD Extreme Scale Systems Center (ESSC) facilities at ORNL.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Final Report on the Viability of Acoustic Techniques for Density and Mass Flow in Enrichment Plants

The primary purpose of this research is to study the viability of acoustic signatures and sensors that could support accurate, noninvasive, and unattended measurement of uranium hexafluoride (UF6) gas density and mass flow rate in situations representative of gaseous centrifuge enrichment plants (GCEP) under safeguards by the International Atomic Energy Agency. The measurement method under development can be readily retrofitted to the exterior of piping in existing systems, potentially facilitating acceptance by facility operators. The findings from the demonstration of a preliminary prototype design and first-generation analysis algorithms in this project will inform the safeguards community as to whether continued acoustic instrument and methods development are warranted. This project conducted two types of measurements, depending on the flow conditions of the gas. For static measurements, the focus is on understanding the acoustic properties at low pressure for non-flowing UF6. For flow measurements, the focus is on understanding the performance of the acoustic system in a simulation of a GCEP pipe using a surrogate gas under fully developed flow conditions. Testing with flowing UF6 was not considered because the cost of building a flowing UF6 system with flow similar to a GCEP was considerably greater than allocated project resources. This document summarizes the results of research conducted over the course of the project. Based on the results to date, the approach is seen to be promising with non-invasive acoustic measurements shown to be possible at the desired pressures. The project has demonstrated • 5% uncertainty of the UF6 gas density for static measurements at pressures relevant to GCEP operations over a 5-minute measurement window with a clear path forward to achieving considerably lower uncertainties over a two-hour period, • Measurements of flowing surrogate gases at pressures relevant to GCEP operations. However, additional research remains to be done to further develop the measurement approach and transition the technology for field use. These development needs include: • Understand the impact of flowing gas on gas density measurements, • Adding sensors and improving the measurement procedure, • Augmenting the data analysis methods by leveraging machine learning approaches, • Reevaluating uncertainty requirements for density and mass flow given observed Online Enrichment Monitor (OLEM) measurement uncertainties, • Designing a robust instrumentation package for field-testing, and • Evaluating measurement results to understand the potential impact on mass balance calculations within GCEP facilities.

42 ENGINEERING↗