Engineering Papers⌕ Search

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.

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

At least 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
NASA NTRS1944

The "limiting line" in mixed subsonic and supersonic flow of compressible fluids

It is well known that the vorticity for any fluid element is constant if the fluid is non-viscous and the change of state of the fluid is isentropic. When a solid body is placed in a uniform stream, the flow far ahead of the body is irrotational. Then if the flow is further assumed to be isentropic, the vorticity will be zero over the whole filed of flow. In other words, the flow is irrotational. For such flow over a solid body, it is shown by Theodorsen that the solid body experiences no resistance. If the fluid has a small viscosity, its effect will be limited in the boundary layer over the solid body and the body will have a drag due to the skin friction. This type of essentially isentropic irrotational flow is generally observed for a streamlined body placed in a uniform stream, if the velocity of the stream is kept below the so-called "critical speed." At the critical speed or rather at a certain value of the ratio of the velocity of the undisturbed flow and the corresponding velocity of sound, shock waves appear. This phenomenon is called the "compressibility bubble." Along a shock wave, the change of state of the fluid is no longer isentropic, although still adiabatic. This results in an increase in entropy of the fluid and generally introduces vorticity in an originally irrotational flow. The increase in entropy of the fluid is, of course, the consequence of changing part of the mechanical energy into heat energy. In other words, the part of fluid affected by the shock wave has a reduced mechanical energy. Therefore, with the appearance of shock waves, the wake of the streamline body is very much widened, and the drag increases drastically. Furthermore, the accompanying change in the pressure distribution over the body changes the aerodynamic moment acting on it and in the case of an airfoil decreases the lift force. All these consequences of the breakdown of isentropic irrotational flow are generally undesirable in applied aerodynamics. Its occurrence should be delayed as much as possible by modifying the shape or contour of the body. However, such endeavor will be very much facilitated if the cause or the criterion for the breakdown can be found first.

MACH NUMBER↗
NASA NTRS2003

Technical Evaluation Report, Part A - Vortex Flow and High Angle of Attack

A symposium entitled Vortex Flow and High Angle of Attack was held in Loen, Norway, from May 7 through May 11, 2001. The Applied Vehicle Technology (AVT) panel, under the auspices of the Research and Technology Organization (RTO), sponsored this symposium. Forty-eight papers, organized into nine sessions, addressed computational and experimental studies of vortex flows pertinent to both aircraft and maritime applications. The studies also ranged from fundamental fluids investigations to flight test results, and significant results were contributed from a broad range of countries. The principal emphasis of this symposium was on "the understanding and prediction of separation-induced vortex flows and their effects on military vehicle performance, stability, control, and structural design loads." It was further observed by the program committee that "separation- induced vortex flows are an important part of the design and off-design performance of conventional fighter aircraft and new conventional or unconventional manned or unmanned advanced vehicle designs (UAVs, manned aircraft, missiles, space planes, ground-based vehicles, and ships)." The nine sessions addressed the following topics: vortical flows on wings and bodies, experimental techniques for vortical flows, numerical simulations of vortical flows, vortex stability and breakdown, vortex flows in maritime applications, vortex interactions and control, vortex dynamics, flight testing, and vehicle design. The purpose of this paper is to provide brief reviews of these papers along with some synthesizing perspectives toward future vortex flow research opportunities. The paper includes the symposium program. (15 refs.)

VORTEX BREAKDOWN↗
DOE OSTI2021

High-speed characterization of two-phase flow and bubble dynamics in titanium felt porous media for hydrogen production

Bubble dynamics and two-phase flow phenomena are closely related to the performance of proton exchange membrane electrolyzer cells (PEMECs). This paper reports an in-situ study of the oxygen bubble behavior and associated multiphase evolutions in the anode side of PEMECs with titanium (Ti) felt liquid gas diffusion layers (LGDLs) by a high-speed visualization system. The micro oxygen bubble dynamics was captured and analyzed at different locations and virous operating conditions. The results show that the bubble detachment frequency and detachment diameter greatly increase with the operating current density. Additionally, they are significantly impacted by the local pore structure and morphology of Ti felt LGDLs. In the flow channels, there exist only several discrete micro bubbles at a low current density (0.04 A/cm 2 ) and a large flow velocity (133 mm/s). At a current density (0.2 A/cm 2 ) and a flow velocity (67 mm/s), a number of gas slugs are formed in the follow channels, in addition to discrete micro bubbles. At a high current density (1 A/cm 2 ) and a flow velocity (67 mm/s), more bubbles appear in the channel, and the flow field is dominated by slug or annular flows. These investigations can help to better understand the two-phase flow and bubble detachment mechanism, and provide a foundation for electrochemical reaction, multiphase flow studies and optimize the design of gas diffusion layers and flow fields for PEMECs in the future.

08 HYDROGEN↗
DOE OSTI2020

Visco-inertial gas flow through wellbore cement fractures

Understanding the nature of fluid flow through fractured wellbore cement is fundamental for evaluating the leakage potential and risk assessments of leaky wellbores. In this study, the conditions that require considering visco-inertial flow for describing the gas flow through wellbore cement fractures were investigated. Nitrogen gas flow tests were conducted on fractured cement samples under varying pressure conditions and flow rates, covering both viscous and visco-inertial flow regimes. The data substantially deviated from Darcy's law at higher flowrates and were well-fit to Forchheimer's equation for visco-inertial flow. The inertial coefficient and critical Reynolds number were expressed as a function of the hydraulic aperture. Finally, the empirical function obtained from the experiments was used as an input to numerical simulations which showed the significant role of visco-inertial flow in wellhead pressure build-up and leakage rates, and demonstrated the importance of visco-inertial flow when modeling gas flow through wellbore cement fractures.

03 NATURAL GAS↗
DOE OSTI2020

Viscometric flow of dense granular materials under controlled pressure and shear stress

This study examines the flow of dense granular materials under external shear stress and pressure using discrete element method simulations. In this method, the material is allowed to strain along all periodic directions and adapt its solid volume fraction in response to an imbalance between the internal state of stress and the external applied stress. By systematically varying the external shear stress and pressure, the steady rheological response is simulated for: (1) rate-independent quasi-static flow; and (2) rate-dependent inertial flow. The simulated flow is viscometric with non-negligible first and second normal stress differences. While both normal stress differences are negative in inertial flows, the first normal stress difference switches from negative to slightly positive, and second normal stress difference tends to zero in quasi-static flows. The first normal stress difference emerges from a lack of coaxiality between a second-rank contact fabric tensor and strain rate tensor in the flow plane, while the second normal stress difference is linked to an excess of contacts in the shear plane compared with the vorticity direction. A general rheological model of second order (in terms of strain rate tensor) is proposed to describe the two types of flow, and the model is calibrated for various values of interparticle friction from simulations on nearly monodisperse spheres. The model incorporates normal stress differences in both regimes of flow and provides a complete viscometric description of steady dense granular flows.

42 ENGINEERING↗
DOE OSTI2023

Finite element modeling of biomass hopper flow

Hoppers are widely used biomass handling devices that channel bulk biomass from storage to subsequent handling equipment. Jenike’s longstanding approach, based on the Mohr-Coulomb model, has been successfully used to design hoppers handling cohesionless granular materials such as grains and other agricultural produces. However, designing a hopper to ensure reliable biomass flow is found to be challenging due to cohesion, irregular particle shape, and bulk material elastoplasticity. This study aims to address the biomass handling engineering challenge with alternative constitutive material models concerning the flow behavior of bulk solids. Finite element modeling is an approach that allows for implementation of different material models, whose underlying constitutive theories assist in investigating the origin and manifestation of bulk mechanical behavior of granular materials. This study focuses on the incipient gravity hopper flow of two types of biomass feedstocks, i.e., ground corn stover and Douglas fir wood. Three widely used constitutive material models, i.e., Mohr-Coulomb model, modified Cam-Clay model, and Drucker-Prager/Cap model, are implemented. Using the flow pattern represented by the volume of biomass exhibiting more than 7% of axial strain (Kamath and Puri, 1999), the finite element model predicts that the bulk corn stover particulate material forms an arch, which represents a hampered transition from the static state to the dynamic flow-state out of the hopper, whereas bulk Douglas fir wood particulate material develops a reliable mass flow pattern. A laboratory scale hopper was used to experimentally determine the biomass flow conditions, which were subsequently compared with the predicted onset of flow by a finite element model (FEM). The developed FEM was found to correctly predict the initiation of mass flow for the milled Douglas fir wood, whereas corn stover was predicted to establish a strong core flow suggesting an unreliable handling characteristic. This observation aligns with the reported poor handling of milled corn stover.

09 BIOMASS FUELS↗
DOE OSTI2020

Validation of Streaklines as Recorders of Synoptic Flow Direction in a Deltaic Setting

Knowledge of the flow patterns within distributary systems is key for understanding deltaic hydro- and morpho-dynamics, yet synoptic measurements of flow fields remain virtually nonexistent. As a means of overcoming this problem, a small number of studies have used biogenic surface films as synoptic flow tracers, under the assumption that biofilm streaklines are tangent to the local flow direction. Here we rigorously test this assumption and show that, despite flow patterns that change severely in space and time (over a range >270°), streaklines are relatively accurate synoptic flow tracers for the Wax Lake Delta, in Louisiana. When the incoming discharge was greater than 2400 m3/s with stable or falling tides, the streakline-derived flow direction departed from near bed flow direction measurements of 22.8° (root mean square). When the discharge was greater than 2400 m3/s and the tides were rising greater than 0.03 m/hr, they were accurate within 28.0°. Under conditions of discharge less than 2400 m3/s and tidal change less than a positive 0.03 m/hr, they were accurate within 33.3°, while during low discharge and rising tides they were accurate within 58.9°. Accuracy varied with distance from the delta, with proximal sites having greater precision. Our results demonstrate that a streakline-derived flow direction can characterize the spatiotemporal variability in the flow directions, but that the accuracy is significantly influenced by the hydrodynamic conditions and location within the network.

58 GEOSCIENCES↗
NASA NTRS1999

The value of assessing pulmonary venous flow velocity for predicting severity of mitral regurgitation: A quantitative assessment integrating left ventricular function

Although alteration in pulmonary venous flow has been reported to relate to mitral regurgitant severity, it is also known to vary with left ventricular (LV) systolic and diastolic dysfunction. There are few data relating pulmonary venous flow to quantitative indexes of mitral regurgitation (MR). The object of this study was to assess quantitatively the accuracy of pulmonary venous flow for predicting MR severity by using transesophageal echocardiographic measurement in patients with variable LV dysfunction. This study consisted of 73 patients undergoing heart surgery with mild to severe MR. Regurgitant orifice area (ROA), regurgitant stroke volume (RSV), and regurgitant fraction (RF) were obtained by quantitative transesophageal echocardiography and proximal isovelocity surface area. Both left and right upper pulmonary venous flow velocities were recorded and their patterns classified by the ratio of systolic to diastolic velocity: normal (>/=1), blunted (<1), and systolic reversal (<0). Twenty-three percent of patients had discordant patterns between the left and right veins. When the most abnormal patterns either in the left or right vein were used for analysis, the ratio of peak systolic to diastolic flow velocity was negatively correlated with ROA (r = -0.74, P <.001), RSV (r = -0.70, P <.001), and RF (r = -0.66, P <.001) calculated by the Doppler thermodilution method; values were r = -0.70, r = -0.67, and r = -0.57, respectively (all P <.001), for indexes calculated by the proximal isovelocity surface area method. The sensitivity, specificity, and predictive values of the reversed pulmonary venous flow pattern for detecting a large ROA (>0.3 cm(2)) were 69%, 98%, and 97%, respectively. The sensitivity, specificity, and predictive values of the normal pulmonary venous flow pattern for detecting a small ROA (<0.3 cm(2)) were 60%, 96%, and 94%, respectively. However, the blunted pattern had low sensitivity (22%), specificity (61%), and predictive values (30%) for detecting ROA of greater than 0.3 cm(2) with significant overlap with the reversed and normal patterns. Among patients with the blunted pattern, the correlation between the systolic to diastolic velocity ratio was worse in those with LV dysfunction (ejection fraction <50%, r = 0.23, P >.05) than in those with normal LV function (r = -0.57, P <.05). Stepwise linear regression analysis showed that the peak systolic to diastolic velocity ratio was independently correlated with RF (P <.001) and effective stroke volume (P <.01), with a multiple correlation coefficient of 0.71 (P <.001). In conclusion, reversed pulmonary venous flow in systole is a highly specific and reliable marker of moderately severe or severe MR with an ROA greater than 0.3 cm(2), whereas the normal pattern accurately predicts mild to moderate MR. Blunted pulmonary venous flow can be seen in all grades of MR with low predictive value for severity of MR, especially in the presence of LV dysfunction. The blunted pulmonary venous flow pattern must therefore be interpreted cautiously in clinical practice as a marker for severity of MR.

Non-NASA Center↗
DOE OSTI2017

Flow Rate Dependency of Permeability Changes with the Deposition of CSH under the Highly Alkaline Ca-rich Condition - 17175

The formation of calcium silicate hydrate (CSH) around the geological repository is one of the notable factors for the assessment of the radionuclide migration. This study discussed the clogging effect of CSH in flow path with the change in the flow rate based on the flow experiments using micro flow cell. In the experiment, the change of the permeability with the alteration of granite under the highly alkaline and Ca-rich condition was examined. As a result, the permeability in flow path gradually decreased with the periodical increase and decrease of the permeability. In addition, under the condition of the lower flow rate (0.06 ml/h), the permeability became much smaller than those for other conditions of the higher flow rates. These suggest that the clogging of flow path will progress with the local deposition/moving of CSH, in particular, under the condition of the lower flow rate. Retardation effect of radionuclide migration will be expected because it is prospected that flow rate of groundwater is lower than that of this study in the actual condition of repository. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Improved Safety Injection Flow Map Associated with Target RCS Depressurization to Maintain Core Coolability of OPR1000

In Korea, Severe Accident Management Guidance (SAMG) for Optimized Power Reactor (OPR1000) has been developed through the research results related to the WOG SAMG, in which the operation phase shifts from the Emergency Operation Procedure (EOP) to the SAMG when the Core Exit Temperature (CET) reaches 923 K. Among 7 mitigation strategies envisioned, the third management strategy, injecting coolant into Reactor Coolant System (RCS) by depressurizing the RCS is one of the key strategies to fortify the in-vessel retention capability. As an aid to such mitigation strategy, a safety injection (SI) flow map has been developed by Lee et al. for predicting required flow rate associated with target RCS depressurization to recover in-vessel coolability. The mechanistic SI flow map used the CET increase rate and core water level decrease rate obtained from MELCOR code simulation. It should be noted that during the severe accident, the total heat generated is sum of the decay heat and the oxidation heat. Thus, it is necessary to consider additional heat sources by oxidation reaction between cladding, structure, control rod, and hot steam prior to the injection of coolant. It was believed that the SI model estimated required flow rate rather conservatively by considering the expected amount decay heat and the accumulated heat in the core based on the RCS pressure. However, it was confirmed that the SI model developed by Lee et al. could not represent the total heat removal especially since start of oxidation. Under transient conditions, the heat sources accumulated in the heat structures (fuel, cladding, supporting structures) should be embarked in the model for accurate contribution of the total heat source to the required flow rate. Therefore, the objective of this study is to improve the previous SI flow map predicting the injection flow rate associated with target RCS depressurization to satisfy the in-vessel coolability during severe accidents. Additional terms were developed to reflect accurate heat accumulating and included in the existing model. Resulting SI flow map includes the target RCS depressurization, which enables discharge of the required flow rate by utilizing the flow rate-pressure curve with operation of two high pressure safety injections (HPSIs). Also recent MELCOR simulation result of a postulated steam generator tube rupture (SGTR) accident verified the accuracy of the improved map.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Heat Transfer Efficiency in Various Prandtl Number Turbulent Flows using DNS approach

Convective heat transfer in turbulent flows utilized in many nuclear reactor designs with various coolants. While the traditional light water reactors (LWR) rely on water coolant, some advanced next-generation designs propose to use gas, molten salt, liquid metal coolants. Due to large presence of LWRs, most of the convective heat transfer correlations are focused on fluids with properties similar to water. However, some of the next-generation reactor designs cannot directly apply those correlations to advanced coolants, such as liquid metal, which has very low Prandtl number. Direct Numerical Simulation (DNS) is one of the promising approaches to study the details of fluid flow and heat transfer using first principles. The advances of high performance computing (HPC) in recent years made it possible to apply DNS to a wide variety of flow conditions. Kim and Moin studied heat transfer in fully developed turbulent channel flow using DNS. Kawamura et al. investigated the turbulent heat transfer in channel flow with low to medium-high Prandtl number fluid (ranging between 0.025 and 5.0). They take Reynolds number effects into consideration as well. However, most of these studies focus on medium and high Prandtl number. Low Prandtl number flow (like liquid metal) DNS studies are very limited. In the presented research the effect of low Prandtl number is studied in the turbulent channel flow. The second-order accurate spatial and time discretization is used for the flow solver in the presented work. The instantaneous velocity and temperature fields are visualized in the simulation. The effect of Prandtl number on the heat transfer is investigated by comparing different turbulent flow cases. In the same hydrodynamic conditions, the heat conduction enhances as the Prandtl number decreases in the low Prandtl number range. The behavior of the thermal boundary layer follows the empirical correlation. As the Prandtl number decreases, the conductive sublayer region becomes wider, where is out of the range of the empirical correlation. In the future work, we plan to estimate the Nusselt number as well to compare the total heat transfer rate with respect to different Prandtl number flow, especially the very low Prandtl number cases.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Onset of instabilities in rotating flows by direct numerical simulation

A rotating disk is the canonical experiment for measuring surface reaction rates in geochemical and electrochemical systems. Using the similarity solution for laminar flow around an infinite disk, the mass transfer coefficient can be simply related to the intrinsic reaction rate at the surface. However, measurements of mass transfer rates use a finite-size disk within a larger container of solution; here the flow is no longer strictly laminar, but there must always be some recirculation. Our interest was initially in the assumption of a uniform radial concentration field, how this breaks down near the perimeter of the disk, and what effect that might have on the measured mass transfer rates. However, our numerical simulations suggest that the flow around a finite-size disk becomes time dependent at Reynolds number ($\textit{Re}$) below 1000, which is much smaller than the typical values in mass-transfer measurements ($\textit{Re}$~10 4 ). In this work, we observe the formation of coherent structures in the flow, which suggest a non-uniform mass transfer at the disk surface. The rotating-disk flow follows a similar sequence of instabilities to the Taylor–Couette flow: a centrifugal instability leading an axisymmetric, time-invariant flow, followed by a Hopf bifurcation to a time-periodic flow. To minimise the possibility that our results are a numerical artefact, we have also simulated the instability in the stationary boundary layer of a rotor–stator flow, comparing with self-similar solutions at low $\textit{Re}$ and with spectral methods near the critical Reynolds number.

42 ENGINEERING↗
DOE OSTI2020

3D modelling of MHD mixed convection flow in a vertical duct with transverse magnetic field and volumetric or surface heating

Understanding phenomena associated with the multiple effects/interactions of the fusion nuclear environment on liquid metal flow is required to correctly design liquid metal (LM) blankets for fusion facilities. These effects are investigated in the present work by numerically simulating 3D LM MHD flow. The simulated geometry consists of a straight, vertical duct which runs perpendicular to a strong, fringing applied magnetic field. There is also a region of applied heating as the primary goal is to explore buoyancy effects in MHD duct flows. Results are presented for both buoyancy assisted (upwards) and buoyancy opposed (downwards) flows in conducting and insulating ducts for a range of Hartmann numbers (Ha) up to 100, Reynolds numbers (Re) from 10 3 to 10 4 and Grashof (Gr) numbers from 10 7 to 10 8 . While increasing Gr or decreasing Re increases buoyancy effects, increasing Ha was shown to increase maximum temperature through turbulence reduction. The extent to which the MHD mixed convection flows are quasi-2D is analyzed and buoyant effects, in competition with electromagnetic forces, are shown to bring about 3D flow features not seen in purely MHD flows. Volumetric nuclear heating with steep gradients is applied to the vertical MHD flows for comparison to flows with surface heating only. Finally, surface heating generates stronger buoyancy effects than volumetric heating of the same total power; however, many of the same phenomena occur. Therefore, surface heating, the only option for lab experiments, can provide indication of the effects of volumetric heating in MHD flows.

42 ENGINEERING↗
DOE OSTI2019

A massively-parallel, unstructured overset method to simulate moving bodies in turbulent flows

An unstructured overset method capable of performing direct numerical simulation (DNS) and large eddy simulation (LES) of many (O(10 5 )) moving bodies, utilizing many computational cores (O(10 5 )), in turbulent, incompressible fluid flow is presented. Unstructured meshes are attached to bodies and placed within a fixed background domain. Body meshes are allowed to arbitrarily overlap and move throughout the domain. Within each mesh a high resolution, unstructured, non-dissipative finite volume method is used to solve for the flow field. Boundary conditions for each mesh are provided by interpolation from flow solutions on overlapping meshes. When many unstructured meshes of different resolution overlap, care is required in the connection between the different flow solutions. An interpolant is created which seeks to preserve volume conservation of flow quantities between meshes regardless of overlapping mesh differences. An implicit fractional step method is used for time advancement, requiring the calculation of a predicted fluid velocity and corrector pressure field. For the predictor step, the resulting interpolation is directly introduced into the implicit equations for the predicted flow field. For the corrected pressure field, the continuity between meshes is weakly enforced using a penalty formulation. The pressure formulation is symmetric, positive-definite and non-singular resulting in a formulation which is readily solvable using traditional iterative matrix inversion techniques. An Arbitrary Euler-Lagrangian (ALE) method coupled to a 6 degrees of freedom rigid body equation system (6-DOF) is used for body motion. For rotation, a quaternion representation is used to solve Euler's equations of rigid body motion. A linear spring damper model, which uses geometry information readily available from the overset assembly process, is used for collisions. Validation of the method for canonical flow fields is presented including assessment of order of accuracy and kinetic energy conservation properties. Particle-resolved direct numerical simulation (PR-DNS) of single particles in various flow fields are presented for validation. PR-DNS results of 50,000 spherical particles freely moving within turbulent channel flow are shown as a demonstration of the method at full scale. LES results of a marine propeller under crashback conditions are shown to demonstrate the ability to simulate highly unsteady turbulent flows over complex, moving geometries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗
DOE OSTI2020

The Complexity of Nonlinear Flow and non-Fickian Transport in Fractures Driven by Three-Dimensional Recirculation Zones

Subsurface fluid flow and solute transport are typically described by first-order or linear rate laws. However, deviations from the first-order rate laws (or anomalous behavior) are typical and lead to nonlinear flow and non-Fickian transport phenomena. The shared underlying mechanisms for these “anomalous” flow and transport phenomena, which are actually typical of fractures, are seldom analyzed. We investigated the co-occurrence of nonlinear flow and non-Fickian transport through numerical simulation experiments of pore-scale flow and solute transport processes through three-dimensional natural fractures. In these fractures, 3-D recirculation zones (RZs) developed and grew with increasing flow rates. Flow nonlinearity resulted from the RZ development and growth, which shrank the main flow channel. Moreover, the same RZs also captured solutes from the main flow channel and released them back later. This retardation resulted in non-Fickian transport. Thus, nonlinear flow and non-Fickian transport are intertwined via their shared dependence on RZs.

58 GEOSCIENCES↗
DOE OSTI2018

Generation of parasitic axial flow by drift wave turbulence with broken symmetry: Theory and experiment

Detailed measurements of intrinsic axial flow generation parallel to the magnetic field in the controlled shear decorrelation experiment linear plasma device with no axial momentum input are presented and compared to theory. The results show a causal link from the density gradient to drift-wave turbulence with broken spectral symmetry and development of the axial mean parallel flow. As the density gradient steepens, the axial and azimuthal Reynolds stresses increase and radially sheared azimuthal and axial mean flows develop. A turbulent axial momentum balance analysis shows that the axial Reynolds stress drives the radially sheared axial mean flow. The turbulent drive (Reynolds power) for the azimuthal flow is an order of magnitude greater than that for axial flow, suggesting that the turbulence fluctuation levels are set by azimuthal flow shear regulation. The direct energy exchange between axial and azimuthal mean flows is shown to be insignificant. Therefore, the axial flow is parasitic to the turbulence-zonal flow system and is driven primarily by the axial turbulent stress generated by that system. The non-diffusive, residual part of the axial Reynolds stress is found to be proportional to the density gradient and is formed due to dynamical asymmetry in the drift-wave turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Oil and gas flow through fractures and along interfaces in well cement

A number of wells at the Strategic Petroleum Reserve (SPR) have shown sustained, positive pressure (referred to as sustained casing pressure or SCP) in the cemented annulus behind the production casing. To better understand how SCP may develop for SPR conditions, we conducted gas and oil flow tests on cement specimens with flaws including cement fractures and discrete interfaces along a cement-steel contact. Many specimens were tested initially with gas, followed by oil, and finally with gas again to identify how the fluid type may affect the flow through flaws in the well cement. Nitrogen was used as the gas, and silicone oil with properties similar to typical crude oil was used in most tests. One set of measurements were made with crude oil. Composite steel-cement specimen with corroded steel were also tested. For both gas and oil tests, the measured flow test data were used to interpret permeability and, assuming the cubic law for flow between parallel plates, the corresponding hydraulic aperture of the flaw. The hydraulic apertures for the flawed specimens ranged from about 20 to >100 μm, which corresponds to permeabilities of about 10 -14 to 10 -12 m 2 , respectively; these hydraulic apertures are consistent with the range of values interpreted from field measurements on leaky wells. Observed differences between the flow of gas and oil were attributed to a number of factors, including non-linear flow of gas, possible blocking of flow paths by solids within the crude oil, two-phase gas and oil flow and the presence of residual oil in the cement flaw. Furthermore, we determined that corroded steel itself is permeable. Using input values consistent with the gas and oil flow measurements, we conducted one-dimensional simulations of gas and oil flow through cemented annulus systems to investigate the role of flaw size and fluid type on the expected response (i.e., pressure build-up) in the annular cement.

02 PETROLEUM↗