Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLOW - CHANNELS”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Screening of resonant magnetic perturbation fields assuming various plasma flow models

A recently updated version of the MARS-F code [Y. Q. Liu et al., Phys. Plasmas 7, 3681 (2000); L. Li et al., Phys. Plasmas 25, 082512 (2018); and G. L. Xia et al., Nucl. Fusion 59, 126035 (2019)] is utilized to numerically investigate the plasma screening effect on the applied resonant magnetic perturbation (RMP) field, assuming various equilibrium flow models, including the toroidal flow, the parallel flow and their combinations, and poloidal and toroidal projections of the parallel flow. A parallel equilibrium flow with a uniform radial profile is found to have no effect on plasma screening of the RMP field. A sheared parallel flow, however, does change plasma screening. The poloidal projection of the parallel flow weakens plasma screening in the resistive-inertial regime. The effect on the favorable average curvature regime is found, however, to be non-monotonic. With the increasing flow speed, the poloidal projection first weakens Glasser-Green-Johnson (GGJ)-screening. Further increase in the flow speed results in enhanced GGJ-screening again. This non-monotonic behavior is related to the perturbed parallel shielding current, which appears also off the mode rational surface at fast flow due to additional resonances between the RMP perturbation and the sound wave continuum. These results indicate that flow induced plasma screening to the RMP field can have complicated characteristics, which, in turn, can have implications on the RMP field penetration into the plasma in experiments for controlling the edge localized modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental investigation of abnormal transverse flow enhancement of α particles in heavy-ion collisions

The mass dependence of the transverse flow for Z = 1–5 fragments from the collisions of 40 Ar + 27 Al, 40 Ar + 48 Ti, and 40 Ar + 58 Ni at 47 MeV/nucleon is investigated experimentally in this article. The transverse flow values are determined using the in-plane components of the fragment transverse momenta, where three conventional methods, i.e., the kinetic flow tensor method, the transverse momentum analysis method, and the azimuthal correlation method, are applied to reconstruct the reaction plane in an event-by-event basis. It is demonstrated from the comparison of the present experimental mass dependent flow measurements and the model simulations using an improved antisymmetrized molecular dynamics model that the experimentally observed abnormal α transverse flow enhancement is closely related to the reaction plane reconstruction procedure in the flow extraction. We further investigate the physical existence of the abnormal α flow behavior using a two-particle azimuthal correlation method, which allows us to provide the relative flow magnitude information with an identification of fragment charge number without the knowledge of the reaction plane differing from the three conventional methods. It is found that the relative flow magnitudes deduced from the two-particle azimuthal correlation functions with an identification of Z, with the correction for the recoil effect imposed by the momentum conservation, show a monotonically increasing trend as a function of fragment charge number, with no exception of the α flow enhancement. These results, in addition to those from the improved antisymmetrized molecular dynamics model simulations, definitely provide experimental evidences for the inexistence of the abnormal α flow behavior in the heavy-ion collisions at the present incident energy region in nature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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↗
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 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

Validation of RELAP-7 for Forced Convective and Natural Circulation Reactor Flows

New reactor designs such as small modular reactors and newly proposed passive safety systems rely on natural circulation two-phase flow cooling of the reactor core under accident conditions. Without a pump to drive the flow, the circulation rate is a strong function of void fraction. Under these conditions where the flow rate, pressure drop, and heat transfer are all closely coupled through the mixture density, the accurate prediction of void fraction is pivotal. Furthermore, while these passive systems boast a more reliable cooling option, the sensitivity of the system to two-phase flow instabilities is more likely than in forced convective cooling and can result in dangerous density wave or flow excursion conditions. Improvement of one-dimensional system codes and further development of three-dimensional code capability are all important efforts in addressing these issues. The proposed work is a critical step towards providing confidence through validation of the capability of thermal-hydraulic tools to meet the emerging challenges in reactor simulation. Historically, system analysis codes develop a stiffness over time due to fine-tuning of adjustable code parameters that result in a reluctance or inability to adopt improved models. Therefore, a rigorous effort should be invested in preparation for RELAP-7 to identify required modeling improvements of constitutive relations and validation of the code applicability to the wide range of possible reactor flow conditions. Part of this validation involves the uncertainty quantification which is a critical step in the Code Scaling, Applicability, and Uncertainty (CSAU) methodology and Best-Estimate Plus Uncertainty (BEPU) licensing method. BEPU allows recovery of safety margins and potentially optimizing plant performance, resulting in economic benefit for the reactor operators. Beginning with the collection of past experimental data covering prototypic LWR conditions, this document reviews three critical datasets available for the validation of system codes. The experimental methods, procedures, facilities, instrumentation, and supplementary information are consolidated to aid the validation efforts of this work and future projects. These three datasets (824 conditions) are chosen due to their wide acceptance and demonstrated application for code validation. Two experimental datasets from the Nuclear Power Engineering Corporation (NUPEC) span prototypic conditions of Boiling Water Reactors (BWR) and Pressurized Water Reactors (PWR). These datasets are labeled as the BFBT (BWR Full-size Fine-mesh Bundle Test) and the PSBT (PWR Subchannel and Bundle Test) datasets respectively. The third dataset is taken in an internally heated annulus that is scaled based on a BWR subassembly and includes two-phase flows important to post-blowdown analysis. Together, the datasets consider the full range of pressure and mass flux expected in LWR forced convective flows from steady-state to long term cooling conditions. Recognizing the deficiency of available data under low-pressure, low-flow, a new natural circulation dataset consisting of 107 conditions is obtained. With system pressures ranging between 145 and 950 kPa, natural circulation is driven by a heat flux up to 275 kW/m2 in a 3 m long heated region of the annulus test section which is followed by a 2 m long unheated region before entering a condenser. Stable natural circulation flow rates in this facility were measured to be between 182 and 590 kg/m2-s. Pressure and temperature are measured at ten locations in the facility. In addition, measurements of void fraction, interfacial area concentration, and gas velocity are obtained at five axial locations in the annulus test section. The RELAP5-3D code demonstrates good accuracy for void fraction predictions in BFBT and PSBT benchmarks, but shows relatively large discrepancy in the annulus experiments at lower pressures. In two-phase natural circulation, large oscillations are observed in many system parameters during the simulation of 19 experiment cases which were conducted at low pressures, and the model shows good predictive performance with the remaining 88 experiment cases. However, some discrepancies in void fraction and temperature predictions are found which may require improved constitutive modeling under these low-pressure, low-flow condition. The quantified uncertainty in void fraction predictions in three forced convection benchmarks due to uncertainties in boundary conditions and geometries are relatively small, except for a few cases in the annulus experiment. The output uncertainty of the natural circulation model is quantified by considering the measurement error in boundary conditions, and it is shown that the output uncertainty ranges are relatively small for the range of conditions considered. Through this work, a clear pathway is provided for validation of RELAP7 once the code is available.

42 ENGINEERING↗
DOE OSTI2022

Assessment of Flow-Enhanced Sensors for Actinide Quantification in MSRs

This report serves as the deliverable for Milestone M2RS-22AN0401062 that is part of Work Package RS-22AN040106 (Flow Enhanced Sensors for MSRs – ANL). The goal of this milestone is to assess the capabilities of flow-enhanced electrochemical sensors (FEES) to quantify actinides in molten salt reactors. Flow enhanced electrochemical sensors are a type of electroanalytical sensor that has been developed at Argonne National Laboratory to be installed directly into MSR flow conduits in order to make measurements of the salt composition. These sensors represent a significant improvement in capabilities compared to earlier electroanalytical sensors that instead can only be operated in quiescent conditions. To support this work, Argonne constructed a molten salt flow system to serve as a testbed for the development of safeguards and process monitoring-relevant sensors in challenging molten salt flow conditions. Through extensive flow testing using UCl3-bearing salts, Argonne has been able to show that the FEES are able to make good measurements of salt composition. Concentration measurements with a mean absolute error of 0.09 wt% have been achieved for representative fuel salt mixtures across a range of UCl3 loadings. Repeatability was also good, with the relative standard deviation from repeated measurements being less than 1.0%. These sensors have also been shown to be able to measure the flow rate of the salt within the flow system’s transfer line. Although the accuracy of the sensors has been good, the sensors have thus far not been able to achieve performance good enough to satisfy the nuclear material accountancy requirements in 10 CFR 74. To do so, accuracy levels of 0.1% must be achieved; the FEES when operated in Argonne’s flow system have so far only been able to achieve relative errors of approximately 5.0%. Much of the error that has been observed at this stage is due to challenges in maintaining a stable salt composition in the molten salt flow system, and it is not believed that there are any inherent flaws in the sensor design. Nonetheless, the high degree of accuracy required by 10 CFR 74 represents a significant challenge, and further design evolution and integration of the sensors into multimodal sensing frameworks will likely be needed to push the measurement accuracy to the needed level.

Moore, Colin E.↗
DOE OSTI2024

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗
DOE OSTI2022

Turbulent flow characteristics in an 84-pin rod bundle for typical and damaged spacer grids

Hexagonal rod bundles arranged in a tightly packed triangular lattice are extensively used for heat transfer and energy generation applications. Staggered spacer grids are used to maintain the structural integrity of gas-cooled fast reactor (GFR) fuel assemblies, while inducing localized turbulence in flow. Damage to these spacer grids results in a disruption of flow fields within these hexagonal fuel bundles. Experimental flow visualizations are critical to identify the differences in local flow properties that the structural damage may cause. This experimental research investigates the flow-field characteristics at a near-wall and center plane in a prototypical 84-pin GFR fuel assembly. Newly installed typical spacers and spacers subject to naturally occurring damage due to material degradation over prolonged experimentation were investigated. Velocity fields were acquired by utilizing the matched-index-of-refraction method to obtain time-resolved particle image velocimetry measurements for a Reynolds number of 12 000. Reynolds decomposition statistical results divulged differences in the time-averaged velocity, velocity fluctuations, flow anisotropy, and Reynolds stress distributions. Galilean decomposition demarcated the influence of spacer grid damage on the velocity fields. To extract turbulent structures and elucidate mechanisms of flow instabilities, proper orthogonal decomposition analysis was employed. Reduced order flow reconstructions enabled the application of vortex identification algorithms to determine the spatial and statistical characteristics of vortices generated. This research work provides unique experimental data on the spacer grid condition-dependent flow. The results offer a deeper understanding of fluid dynamics behavior to support GFR rod bundle design efforts and computational fluid dynamics model validation.

36 MATERIALS SCIENCE↗
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 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↗
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 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↗