Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLOW VELOCITY, HIGH - SUPERSONIC”

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

MHD flow in liquid metal blankets: Major design issues, MHD guidelines and numerical analysis

The design of breeding blankets represents the major challenge for fusion reactor engineering because of performance requirements and severe operating conditions in terms of heat load and neutron flux. Liquid metal alloys such as lead-lithium, due to their lithium content, can be used to breed tritium, one of the plasma fuel components, and owing to their high thermal conductivity, they may serve as coolants. On the other hand, there are technical issues related to the fact that the liquid metals are electrically conducting and interact with the plasma-confining magnetic field. Induced electric currents and generated electromagnetic forces affect velocity and pressure distribution in the blankets. Magnetohydrodynamic (MHD) flows for fusion applications have been often investigated in simplified geometries, such as pipes, ducts, bends, with focus on their fundamental features. These analyses are essential, since results remain valid as background for the development of blanket designs, even when a concept is dismissed. However, the conceptual study of fusion blankets requires to take into account the global multiple effects, that arise when the full system is considered. Progress made in fusion-related MHD research results from combined numerical and experimental activities. In this paper we review and summarize features of 2D and 3D MHD flows that are typical in liquid metal blankets, together with available correlations for MHD pressure losses. This knowledge can provide simple design MHD guidelines that support a preliminary estimate of MHD effects in a blanket concept, in terms of pressure drop and flow distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

CFD simulations of Molten Salt Fast Reactor core cavity flows

Computational Fluid Dynamics (CFD) has become increasingly important in the research and development of advanced nuclear reactors. Here, in the current study, extensive CFD simulations were conducted for the coolant flow in Molten Salt Fast Reactor (MSFR) core models using the state-of-the-art spectral element flow solver Nek5000 and multiscale coarse-mesh thermal-hydraulic software Pronghorn. The underlying motivation is to seek an in-depth understanding of how the internal velocity distribution can be influenced by the MSFR core cavity shape, the Reynolds number, turbulence modeling options and the inlet boundary conditions. The CFD techniques involved in this investigation range from coarse-mesh CFD, RANS modeling, to the high-fidelity LES calculations. Specifically, a series of RANS simulations were performed for the 2-D axisymmetric core model and 3-D wedge domains to study the flow distribution inside the MSFR core. It is observed that a proper representation of the MSFR inlet channel duct is important for the prediction of internal flow distribution. It is also showcased here how researchers can leverage the Nek5000 CFD results to calibrate more efficient coarse-mesh CFD tools, like Pronghorn, for the actual MSFR design needs. Moreover, this paper highlights a 3-D LES model for an entire MSFR core using the spectral element method and demonstrates the feasibility of this modeling approach. The readiness and potential limitations of the RANS approach are examined with respect to the high-fidelity LES simulations. The present investigation lays a solid foundation as we are leveraging the high-fidelity CFD capabilities to inform MSFR design efforts.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Fabrication and Preliminary Demonstration of Microwave Resonant Cavity Transducer Performance

We are investigating a microwave cavity-based transducer for in-core high-temperature fluid flow sensing in molten salt cooled reactors (MSCR) and sodium fast reactors (SFR). This sensor is a hollow metallic cylindrical cavity, which can be fabricated from stainless steel, and as such is expected to be resilient to radiation, high temperature and corrosive environment of MSCR and SFR. The principle of sensing consists of making one wall of the cylindrical cavity flexible enough so that dynamic pressure, which is proportional to fluid velocity, will cause membrane deflection. Membrane deflection causes cavity volume change, which leads to a shift in the resonant frequency. We have developed an initial design for proof-of-principle testing of the flow sensor performance in microwave K-band. A cylindrical resonator prototype was fabricated from brass for initial tests in water. The cavity size is matched to the flange of a standard WR-42 waveguide. Microwave field is coupled into the resonant cavity through a subwavelength-size aperture. A test article was developed, consisting of a piping Tee with bulkhead WR-42 microwave waveguide installed in leak-proof design. In the test article, the cylindrical cavity is positioned in the center of the pipe. A microwave waveguide circulator was installed in the setup to suppress microwave reflections at the cavity entrance. Preliminary spectral characterization of cavity spectral response was performed with microwave VNA. Applying mechanical pressure to cavity membrane showed a measurable shift in the microwave resonant frequency. We also investigate mechanical integrity of the flowmeter’s membrane through computer simulations. By calculating the stress on the plate due to deflection and, comparing the stress to the material ultimate tensile strength and yield strength, it can be estimated if the plate will fail. The stress on the plate was calculated with an analytic closed form solution model, and with COMSOL Structural Mechanics Module which does not involve any approximations. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Flow through negatively charged, nanoporous membranes separates Li + and K + due to induced electromigration

Flow through negatively charged nanopores separates Li + and K + with selectivities of up to 70 and Li + passages from 20% to above 100%. Remarkably, both the Li + /K + selectivity and Li + passage initially increase with flow rate, breaking the permeability/selectivity trade-off. Modelling demonstrates that flow through the membranes creates electric fields that retard transport of cations. Selectivity increases with flow rate because the K + electromigration velocity exceeds its convective velocity, but for Li + electromigration is weaker than convection. Here, modelling also shows the importance of controlling concentration polarization. With further work, related separations might provide highly pure Li salts for battery manufacturing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Heat transfer from glass melt to cold cap: Computational fluid dynamics study of cavities beneath cold cap

Efficient glass production depends on the continuous supply of heat from the glass melt to the floating layer of batch, or cold cap. Computational fluid dynamics (CFD) are employed to investigate the formation and behavior of gas cavities that form beneath the batch by gases released from the collapsing primary foam bubbles, ascending secondary bubbles, and in the case of forced bubbling, from the rising bubbling gas. The gas phase fraction, temperature, and velocity distributions below the cold cap are used to calculate local and average heat transfer rates as a function of the bubbling rate. It is shown that the thickness of the cavities is nearly independent of the cold cap shape and the amount of foam evolved during batch conversion. Furthermore, it is ~7 mm and up to ~15 mm for the cases without and with forced bubbling used to promote circulation within the melt, respectively. Using computed velocity and temperature profiles, the melting rate of the simulated high-level nuclear waste glass batch was estimated to increase with the bubbling rate to the power of ~0.3 to 0.9, depending on the flow pattern. The simulation results are in good agreement with experimental data from laboratory- and pilot-scale melter tests.

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

Strongly Heated Turbulent Flow in a Channel with Pin Fins

Large-eddy simulations (LES) were performed to study the turbulent flow in a channel of height H with a staggered array of pin fins with diameter D = H/2 as a function of heating loads that are relevant to the cooling of turbine blades and vanes. The following three heating loads were investigated—wall-to-coolant temperatures of Tw/Tc = 1.01, 2.0, and 4.0 - where the Reynolds number at the channel inlet was 10,000 and the back pressure at the channel outlet was 1 bar. For the LES, two different subgrid-scale models—the dynamic kinetic energy model (DKEM) and the wall-adapting local eddy-viscosity model (WALE)—were examined and compared. This study was validated by comparing with data from direct numerical simulation and experimental measurements. The results obtained show high heating loads to create wall jets next to all heated surfaces that significantly alter the structure of the turbulent flow. Results generated on effects of heat loads on the mean and fluctuating components of velocity and temperature, turbulent kinetic energy, the anisotropy of the Reynolds stresses, and velocity-temperature correlations can be used to improve existing RANS models.

42 ENGINEERING↗
DOE OSTI2020

A Model of Ice Wedge Polygon Drainage in Changing Arctic Terrain

As ice wedge degradation and the inundation of polygonal troughs become increasingly common processes across the Arctic, lateral export of water from polygonal soils may represent an important mechanism for the mobilization of dissolved organic carbon and other solutes. However, drainage from ice wedge polygons is poorly understood. We constructed a model which uses cross-sectional flow nets to define flow paths of meltwater through the active layer of an inundated low-centered polygon towards the trough. The model includes the effects of evaporation and simulates the depletion of ponded water in the polygon center during the thaw season. In most simulations, we discovered a strong hydrodynamic edge effect: only a small fraction of the polygon volume near the rim area is flushed by the drainage at relatively high velocities, suggesting that nearly all advective transport of solutes, heat, and soil particles is confined to this zone. Estimates of characteristic drainage times from the polygon center are consistent with published field observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Modeling Flow in a Microreactor Core: From High-Fidelity CFD to Subchannel Analysis

The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2023

Demonstration of Microwave Cavity Transducer Performance in High Temperature Fluid Flow Sensing: Development of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with one of the flat walls of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the microwave resonant frequency, which can be detected with a spectrum analyzer. We have performed a preliminary proof-of-principle test of flow sensing in high temperature liquid sodium in environment, which is similar to that of an advanced reactor. The liquid sodium setup consists of a cylindrical vessel with a center feed line, where a transducer inserted through the top cap of the vessel measures velocity of the impinging liquid jet. For this test, we have developed a cylindrical resonator, which was machined from stainless steel 316 and electroplated with silver on the interior surfaces. Outer diameter of the cylindrical cavity is approximately 3cm, and the thin wall is approximately 200µm thick. We have also developed and insertion probe consisting of a 50cm WR-42 brass waveguide enclosed in a protective SS316 tube. The cavity was excited through a WR-42 waveguide through a subwavelength hole on the side of the wall of the cylinder in the TE 011 mode with resonant frequency f ≈ 17.8GHz. As a calibration experiment, we have assembled a water vessel with a center feed, and with dimensions similar to those of the liquid sodium setup. Frequency shift of the cavity spectral response was obtained at room temperature and ambient pressure by gradually increasing water flow rate from 0 to 25GPM (gallons per minute). Corresponding monotonic increase of resonant frequency by several hundred KHz was observed. Next, a sodium flow rate sensing study was performed in a liquid sodium vessel at 340° C temperature and ambient pressure. The flow rate was changed by varying sodium pump power, and a frequency shift of several hundred KHz was observed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Hydrodynamics of a multicomponent vesicle under strong confinement

Herein we numerically investigate the hydrodynamics and membrane dynamics of a multicomponent vesicle in two strongly confined geometries. This serves as a simplified model for red blood cells undergoing large deformations while traversing narrow constrictions. We propose a new parameterization for the bending modulus that remains positive for all lipid phase parameter values. For a multicomponent vesicle passing through a stenosis, we establish connections between various properties: lipid phase coarsening, size and flow profile of the lubrication layers, excess pressure, and the tank-treading velocity of the membrane. For a multicomponent vesicle passing through a contracting channel, we find that the lipid always phase separates so that the vesicle is stiffer in the front as it passes through the constriction. For both cases of confinement we find that lipid coarsening is arrested under strong confinement, and resumes at a high rate upon relief from extreme confinement. The results may be useful for efficient sorting lipid domains using microfluidic flows by controlled release of vesicles passing through strong confinement.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Evaluating Transfer and Pumping of Slurries from Pulsed Jet Mixed Vessels

Due to gravity, solids in slurries will settle if density differences between the solids and liquid are positive (i.e., particle has a negative buoyant force) unless rheological properties and flow conditions are adequate to overcome the gravitational effects. The rate of settling depends on the force balance of the particle, which includes the surface forces associated with fluid rheology. Given the same fluid and solid properties, particles that are larger and denser tend to settle faster. When pumping slurry into a vessel at concentrations precluding hindered settling with insufficient mixing, particle and density distributions can result in preferential settling, creating stratification in the solids concentration within the vessel. For vessels with transfer line inlets located in the lower portion of the tank, the stratified solids concentration may be detrimental to the transfer system performance. Elevated concentrations of solids in the slurry entrained at the inlet to the transfer line can result in the effective viscosity or slurry bulk density exceeding the design limits of the pump. These conditions could result in plugging of the transfer line or onset of cavitation of the pumps because of excessive pressure drop. These conditions can be exacerbated with periodic inlet conditions existing at the transfer line inlet. Periodic conditions can result when vessel mixing is intermittent such as with pulsed jet mixers (PJM). The transfer line inlet conditions are impacted by the periodic nature of the PJM operations with respect to suspension of solids and their transport to the inlet of the transfer line. A scaling approach is presented, and corresponding test requirements are developed for assessing the prevention of plugging the pipeline. Line plugging mechanisms are addressed that do not include plugging due to steady-state high-density slurry entering the transfer line and reducing the net positive suction head available (NPSHA) at the pump inlet to below that required for pump operation. Items considered include the transition to reduced relative flow velocities, such that the critical pipe velocity for solids deposition, Ucd, is not maintained, and segregation of heavy solids during the transport. The recommended requirements to prevent plugging include: • Limits for viscosity and density for entrained slurry to prevent the pressure drop in the pipeline from exceeding pump capacity. • Limits for viscosity and density for entrained slurry to prevent the net positive suction head available (NPSHA) from falling below the net positive suction head required (NPSHR) for operating the pump. • Transfer line velocity and flow rate requirements to maintain solids in suspension, while avoiding line plugging that results from deposition of solids within the transfer line. This paper describes the development of the scaling and testing requirements to verify that proposed approaches for transfer and pump out are appropriately developed

periodic flow, mixing, pipeline plugging, scaling,↗
DOE OSTI2026

Background-Oriented Schlieren Velocimetry of Helium Coolant Flow in Additively Manufactured Channels

High-pressure helium gas cooling is an attractive solution for thermal management of the fusion blanket first wall, as this coolant is chemically and neutronically inert and separable from hydrogenic species. However, due to the low thermal mass of helium, geometric optimization of these channels is required to provide sufficient cooling at manageable flow rates and pumping burdens. Increasingly, analysis and optimization of these coolant channels rely on computational fluid dynamics (CFD) simulations, and these require relevant experimental data for turbulence model validation. Toward this end, a high-pressure helium gas flow visualization system has been employed to image the flow of helium in flow channels with one-sided heating, mimicking the blanket first wall environment. Flow of helium at 4 MPa pressure and flow rates up to 68 g/s (Reynolds number 57 000) is supplied to rectangular channel test sections, with uniform heating applied to the bottom wall of the channel at heat fluxes varied between roughly 50 and 130 kW/m2. A high-speed camera is used to image index of refraction (IOR) gradients in the fluid via background-oriented schlieren (BOS), and temperature and pressure instrumentation are used to characterize thermal-hydraulic performance of each channel. Cross correlation of time-resolved BOS images is then used to calculate time-averaged 2-D helium velocity fields. Flow in additively manufactured (AM) channels is examined in this manner, including both featureless channels and those containing baffling as a heat transfer enhancement. The flow distribution seen in the featureless case differs significantly from that seen in prior simulations, whereas the flow in the baffled case shows the predicted behavior of flow forced along the heated wall. This augmented flow distribution is seen to increase the heat transfer coefficient in the baffled test section. Here, strategies are discussed for ongoing and future validation of these simulations, with the aim of model deployment for blanket cooling design and optimization.

Additive manufacturing↗
DOE OSTI2024

Electromagnetic Inhomogeneous-Energy-Density-Driven Instability in the UCLA Large Plasma Device (EM-IEDDI in LAPD) (Final Technical Report)

Using the Large Plasma Device (LAPD), the electromagnetic signatures associated with the inhomogeneous-energy-density-driven (IEDD) instability were targeted for documentation in a lab experiment. Exploitable features of LAPD are the extreme-length, high-beta, low-collisionality, and excellent reproducibility. The signatures of interest were mode amplitude, mode structure, frequency spectrum, phase velocity, and excitation threshold. The abrupt inhomogeneity at the interface between two coaxial plasma columns, that have disparate space potential values, were expected to generate a strong, localized, transverse, dc electric field and form a strong ExB-flow inhomogeneity that was expected to trigger broadband, intermediate-frequency waves that propagate many gyroradii beyond the interfacial source region. WVU was unable to destabilize the IEDDI mode in LAPD. Experimental limitation of three different parameters conspire to prevent IEDDI-favorable conditions even though Alfven-favorable conditions are met. Electron-neutral collision frequency is too high, Structuring of the radial profile of the radial electric field is negligible and the control of that structure is minimal, and the Magnetic field is too large for some purposes and too small for other purposes. WVU designed a new antenna and relied on externally-driving the IEDDI (instead of spontaneous destabilization). The antenna operation failed for unforeseen reasons and lack of electrician time and availability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Field validation of an actuated convergent-beam acoustic Doppler profiler for high resolution flow mapping

Abstract Energetic river and tidal flow environments feature complex hydrodynamic conditions. Conventional acoustic Doppler profiling instrumentation typically requires assumptions of flow homogeneity over the spatial scales of the divergent beam separation. This removes the ability to measure spatio-temporal variability within the flow. However, velocity variability within these spatial scales is often important in the design of structures and devices exposed to such flow, informing dynamic and peak hydrodynamic load predictions. The research presented outlines the development and testing of a flow measurement instrument consisting of multiple spatially-separated single beam acoustic Doppler profilers converging on a remote focal point, the location of which can be programmatically adjusted through actuation. This increases the spatial resolution at which remote field measurements can be made in energetic flow environments. Field testing of the instrument was conducted in a tidal channel at Sequim Bay Inlet, WA, USA. Results are compared with independent reference measurements made by a co-located, motion-corrected acoustic Doppler velocimeter. This comparison, across tidal velocities of 0.4–0.7 m s −1 , showed a mean velocity error of 0.5%–13.2% across nine sample locations within a 3.00 m × 2.25 m plane. To the authors’ knowledge, this is the first publication of flow velocity results from non-intrusive acoustic measurements at an off-axis, focal point location that features 3D positional control.

13 HYDRO ENERGY↗
DOE OSTI2023

Analytical Study and Experimental Verification of Shear-Thinning Ink Flow in Direct Ink Writing Process

Abstract Direct ink writing (DIW) process is a facile additive manufacturing technology to fabricate three-dimensional (3D) objects with various materials. Its versatility has attracted considerable interest in academia and industry in recent years. As such, upsurging endeavors are invested in advancing the ink flow behaviors in order to optimize the process resolution and the printing quality. However, so far, the physical phenomena during the DIW process are not revealed in detail, leaving a research gap between the physical experiments and its underlying theories. Here, we present a comprehensive analytical study of non-Newtonian ink flow behavior during the DIW process. Different syringe-nozzle geometries are modeled for the comparative case studies. By using the computational fluid dynamics (CFD) simulation method, we reveal the shear-thinning property during the ink extrusion process. Besides, we study the viscosity, shear stress, and velocity fields, and analyze the advantages and drawbacks of each syringe-nozzle model. On the basis of these investigations and analyses, we propose an improved syringe-nozzle geometry for stable extrusion and high printing quality. A set of DIW printing experiments and rheological characterizations are carried out to verify the simulation studies. The results developed in this work offer an in-depth understanding of the ink flow behavior in the DIW process, providing valuable guidelines for optimizing the physical DIW configuration toward high-resolution printing and, consequently, improving the performance of DIW-printed objects.

Engineering↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗
DOE OSTI2020

Identifying recharge processes into a vast “fossil” aquifer based on dynamic groundwater 81 Kr age evolution

Water in deep aquifers in arid regions is often considered to be "fossil" when modern recharge rates are negligible relative to the reservoir capacity. Over the past five decades, the Nubian Sandstone Aquifer (NSA) in the arid region of the Sinai Peninsula (Egypt) and the Negev Desert (Israel) has been considered to contain fossil water based on 14 C dating, which revealed 14 C ages of about 30 kyr over most of the aquifer. However, this relatively homogeneous age distribution contradicts the expected increase in groundwater age in the direction of decreasing piezometric head along the flow trajectories. Here, dating results with the longer-lived 81 Kr radioisotope (t 1/2 = 229 ± 11 kyr) are presented, highlighting a wide age range of 40 kyr to 630 kyr in the confined sections of the aquifer, all with very low 14 C activity (<1 pmC). Elevated 81Kr and 14 C activities were only observed within or close to the system's recharge areas. These findings support a new perception of groundwater replenishment during different epochs from the early mid- Pleistocene to the Holocene. By tracking the downstream age evolution, rejuvenation was identified in places where the confinement had been breached. At other locations, the existence of an older groundwater body contributing to the aquifer was detected by means of strongly depleted 81 Kr activity. High spatial heterogeneity in groundwater ages close to the discharge zone of the system is attributed to pronounced age stratification with depth. Calculated ages in the more isolated sections of the system were used to assess regional flow velocity, hydraulic conductivity, and their agreement with present recharge rates. We conclude that regional aquifers in which low 14 C activity suggests homogeneous ancient groundwater age should be reevaluated with 81Kr. With an effective age range beyond one million years, this may enable the restoration of recharge history and provide crucial information for the management of groundwater resources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Clogging, dynamics, and reentrant fluid for active matter on periodic substrates

Here, we examine the collective states of run-and-tumble active matter disks driven over a periodic obstacle array. When the drive is applied along a symmetry direction of the array, we find a clog-free uniform liquid state for low activity, while at higher activity, the density becomes increasingly heterogeneous and an active clogged state emerges in which the mobility is strongly reduced. For driving along nonsymmetry or incommensurate directions, there are two different clogging behaviors consisting of a drive-dependent clogged state in the low activity thermal limit and a drive-independent clogged state at high activity. These regimes are separated by a uniform flowing liquid at intermediate activity. There is a critical activity level above which the thermal clogged state does not occur, as well as an optimal activity level that maximizes the disk mobility. Thermal clogged states are dependent on the driving direction while active clogged states are not. In the low activity regime, diluting the obstacles produces a monotonic increase in the mobility; however, for large activities, the mobility is more robust against obstacle dilution. We also examine the velocity-force curves for driving along nonsymmetry directions and find that they are linear when the activity is low or intermediate but become nonlinear at high activity and show behavior similar to that found for the plastic depinning of solids. At higher drives, the active clustering is lost. For low activity, we also find a reentrant fluid phase, where the system transitions from a high mobility fluid at low drives to a clogged state at higher drives and then back into another fluid phase at very high drives. We map the regions in which the thermally clogged, partially clogged, active uniform fluid, clustered fluid, active clogged, and directionally locked states occur as a function of disk density, drift force, and activity.

36 MATERIALS SCIENCE↗