Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow instability”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A subgrid model for electron-scale turbulent transport in global ion-scale gyrokinetic simulations of tokamak plasmas

A subgrid electron–temperature–gradient (ETG) model is demonstrated here which averages local electron-scale turbulence over intermediate scales in space and time to include in global ion–temperature–gradient (ITG) simulations. This approach results in ion-scale equations, which incorporate the electron heat transport from ETG turbulence and the effects of electron-scale turbulence on the ion scale. Flux-tube ETG Cyclone Base Case simulations are performed at different radial locations and a kinetic form of the flux is added to global ion-scale (IS) simulations as a source term. Analytic radial profiles of ETG heat flux are constructed and compared to flux-tube results at multiple radial locations. Different ratios of ITG to ETG heat flux levels are considered and the results of capturing ETG heat transport in global ion-scale simulations are discussed. The potential coupling of the ETG streamer potential and intermediate-scale zonal flows to the ion scale is further addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Time-resolved measurement of neutron energy isotropy in a sheared-flow-stabilized Z pinch

Previous measurements of neutron energy using fast plastic scintillators while operating the Fusion Z Pinch Experiment (FuZE) constrained the energy of any yield-producing deuteron beams to less than 4.65 keV. FuZE has since been operated at increasingly higher input power, resulting in increased plasma current and larger fusion neutron yields. A detailed experimental study of the neutron energy isotropy in these regimes applies more stringent limits to possible contributions from beam-target fusion. The FuZE device operated at -25 kV charge voltage has resulted in average plasma currents of 370 kA and D–D fusion neutron yields of $4\times10^7 \pm 4\times10^6$ neutrons per discharge. Measurements of the neutron energy isotropy under these operating conditions demonstrates the energy of deuteron beams is less than $7.4 \pm 5.6^\mathrm{(stat)} \pm 3.7^\mathrm{(syst)}$ keV. Characterization of the detector response has reduced the number of free parameters in the fit of the neutron energy distribution, improving the confidence in the forward-fit method. Gamma backgrounds have been measured and the impact of these contributions on the isotropy results have been studied. Additionally, a time dependent measurement of the isotropy has been resolved for the first time, indicating increases to possible deuteron beam energies at late times. This suggests the possible growth of m = 0 instabilities at the end of the main radiation event but confirms that the majority of the neutron production exhibits isotropy consistent with thermonuclear origin.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Importance of High Valence Element Nb in Ni-Rich Layered Cathodes for High-Voltage Lithium-Metal Batteries

Ni-rich layered cathode materials have attracted extensive attention due to their higher energy density and technological maturity in commercialization. As the nickel content is raised, especially surpassing 80%, the increased energy density comes with the tradeoff of diminished thermal stability and increased electrochemical structural instability of the cathode. Compared with Co, Al, B, and Ta, the introduction of high valence element Nb significantly improved the electrochemical cycling, delivering a capacity of 202 mAh/g, corresponding to a capacity retention of 92% after 200 cycles tested at 45 °C. Further, the ex situ differential scanning calorimetry and in situ isothermal microcalorimetry demonstrate that the Nb-modified cathode has the potential to enhance the safety of ultrahigh nickel (Ni) NMCs and displays remarkable resilience to extensive cycling by inhibiting high-temperature decomposition reactions and exhibiting a lower heat flow during electrochemical cycling.

25 ENERGY STORAGE↗
DOE OSTI2024

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

MARS-Q modeling of low- n resistive kink-peeling modes and edge harmonic oscillations in DIII-D

Linear and quasilinear magnetohydrodynamic (MHD) modeling is carried out for two DIII-D discharges that both featured a transition from quiescent H-mode (QH) to wide-pedestal QH (WPQH). The MHD perturbations, associated with the edge harmonic oscillations (EHOs) observed during the QH-phase in both discharges, are identified as low-n (n is the toroidal mode number) resistive kink-peeling instabilities, with the ideal MHD counterpart remaining stable. The quasilinear model successfully simulates EHO-like perturbations during the QH phase in both discharges, confirming the experimental observations. A less intuitive finding is the EHO-like behavior involving the n = 2 perturbation, simulated for one of the discharges during the WPQH-phase. This result, while consistent with experimental observations, is obtained despite the fact that the initial perturbation is linearly stable. The eventual growth of the perturbation is solely due to nonlinear interaction between the MHD perturbation and the plasma toroidal flow. The occurrence of the EHO-like perturbation during the WPQH-phase is found to be sensitive to the initial profile of the plasma edge rotation. For DIII-D plasmas considered, the neoclassical toroidal viscosity, generated by three-dimensional low-n perturbations, is found to play a dominant role in modifying the edge flow during EHOs.

EHO↗
DOE OSTI2023

A sui generis whipping-instability-based self-sequencing multi-monodisperse 2D spray from an anisotropic microfluidic liquid jet device

Well-defined aerosols pave the way for versatile basic and applied research. Here, we demonstrate a unique whipping instability that generates from a high-aspect-ratio microfluidic device resulting in a unique steady-state gas-focused whipping jet (WJ) without any need for electrification. This WJ device emanates a multi-monodisperse whipping spray jet with a two-dimensional (2D) profile. We demonstrate this phenomenon based on various fluidic parameters theoretically and experimentally. The 2D WJ’s unique behavior is derived using analytical fluid dynamics to explain jet diameter, whipping regime, and spreading angle. The phenomenon is further characterized experimentally by measuring the angle with respect to the flow rate, the distances between droplets, the droplet shapes, and the reproducibility of these parameters. We also explain the precise fabrication of such inexpensive devices. Lastly, we highlight these devices’ potential use as sample environments in versatile applications ranging from cryoelectron microscopy over mass spectrometry to drug formulation and structural studies at X-ray free-electron lasers.

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

A careful examination of closure models in Euler–Lagrange Simulations of compressible multiphase flow in a planar shock particle curtain problem

In this work we present a comprehensive investigation of state-of-the-art closure models employed to represent interphase momentum, thermal, and work exchange between the gas and particulate phases for Euler–Lagrange (EL) simulations in shock-driven flows. A complete list of closures for the force, torque, heat transfer, and work exchange models is provided. In particular, the present work includes a stochastic closure for the particle-to-particle variation in the quasi-steady force and a deterministic closure for particle-to-particle variation in the added mass force in an EL framework. These variations arise due to the presence of neighboring particles and particle–particle interactions. To investigate the importance of each closure term, we carry out fully three-dimensional simulations for a planar shock propagating over a random bed of inert particles. The primary goal is to evaluate the role of each closure term on the gas dynamic features (such as transmitted and reflected shock locations) and particle curtain features (such as upstream and downstream curtain locations). To this end, thirteen cases are considered, with each case progressively including a closure model with the goal to identify and quantify its contribution to the simulated dynamics. We show that the volume fraction dependence of the mean force models plays an important role in generating wave-like instabilities that lead to concentration bands. In addition, fluctuations in quasi-steady and added mass forces primarily decrease the internal instabilities that tend to enhance local volume fraction variations. Particle rotation is primarily due to inter-particle collisions, is generally weak, and does not play an important role in the translational dynamics for the present configuration. Inter-phase heat transfer has a strong effect on gas phase temperature, slows down the transmitted and reflected shocks, and decreases the width of the curtain. Furthermore, the absence of a work-coupling model fails to conserve the total energy, greatly under-predicts the gas temperature which in turn affects the particle dynamics.

Compressible flow↗
DOE OSTI2021

Symmetry-breaking-induced rare fluctuations in a time-delay dynamic system

Inspired by the experimental and numerical findings, we study the dynamic instabilities of two coupled nonlinear delay differential equations that are used to describe the coherent oscillations between the top and bottom boundary layers in turbulent Rayleigh–Bénard convection. By introducing two sensitivity parameters for the instabilities of the top and bottom boundary layers, we find three different types of solutions, namely in-phase single-period oscillations, multi-period oscillations and chaos. The chaos solution contains rare but large amplitude fluctuations. The statistical properties of these fluctuations are consistent with those observed in the experiment for the massive eruption of thermal plumes, which causes random reversals of the large-scale circulation in turbulent Rayleigh–Bénard convection. Finally, our study thus provides new insights into the origin of rare massive eruptions and sudden changes of large-scale flow pattern that are often observed in convection systems of geophysical and astrophysical scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Analytical model for balling defects in laser melting using rivulet theory and solidification

In laser welding and additive manufacturing communities, the balling (humping) defect is primarily attributed to the Plateau–Rayleigh fluid instability (PRI) with a few authors suggesting fluid jetting and volume conservation as alternative mechanisms. As analytical descriptions of these mechanisms are unavailable, combining them into a single formalism is unfeasible. We present a new model of PRI with higher accuracy, accounting for competition with solidification, to compare the expected behavior with known experimental trends when fluid jetting is neglected. We adapt a rivulet instability model from fluid physics to account for the stabilizing effects of the substrate, which the traditional cylindrical-jet geometry does not account for and estimate the instability growth rate. Our model yields a continuous transition from non-balling to balling hitherto lacking in current literature and predicts instability growth at higher wavelengths with strong sensitivity to the solidification front curvature. While the fluid surface is most unstable for shallow melt pools, the absolute magnitude of balling relevant to printing defects scales with melt pool depth and has a maximum for a given melt pool geometry. Synchrotron-based x-ray radiography of thin samples indicates that PRI growth rates and solidification can be comparable in magnitude and thus compete, as we find in this work. We predict that deviations between model predictions and our experimental results demonstrate the importance of fluid flows and heat transport in the balling process. Our experiments further demonstrate at least one mechanism by which the melt pool length and the balling wavelength are not equivalent, as commonly claimed.

Taylor, Z. [Stanford Univ., CA (United States); SL↗
DOE OSTI2024

Setups for eliminating static charge of the ATLAS18 strip sensors

Abstract Construction of the new all-silicon Inner Tracker (ITk), developed by the ATLAS collaboration to beable to track charged particles produced at the High-Luminosity LHC, started in 2020 and is expected tocontinue till 2028. The ITk detector will include 18,000 highly segmented and radiation hard n+-in-psilicon strip sensors (ATLAS18), which are being manufactured by Hamamatsu Photonics. Mechanical and electrical characteristics of produced sensors are measured upon their delivery at several institutes participating in a complex Quality Control (QC) program. The QC tests performed on each individual sensor check the overall integrity and quality of the sensor.During the QC testing of ATLAS18 strip sensors, an increased number of sensors thatfailed the electrical tests was observed. In particular, IV measurements indicated an early breakdown,while large areas containing several tens or hundreds of neighbouring strips with low interstrip isolationwere identified by the Full strip tests, and leakage current instabilities were measured in a long-term leakage current stability setup. Moreover, a high surface electrostatic charge reaching a level of several hundreds of volts per inch was measured on a large number of sensors and on the plastic sheets, whichmechanically protect these sensors in their paper envelopes. Accumulated data indicates a clearcorrelation between observed electrical failures and the sensor charge-up.To mitigate the above-described issues, the QC testing sites significantly modified the sensor handlingprocedures and introduced sensor recovery techniques based on irradiation of the sensor surface with UVlight or application of intensive flows of ionized gas. In this presentation, we will describe the setupsimplemented by the QC testing sites to treat silicon strip sensors affected by static charge andevaluate the effectiveness of these setups in terms of improvement of the sensor performance.

Instruments & Instrumentation↗
DOE OSTI2023

An Experimental Comparison of Cyclic Variations in Diesel–Natural Gas and POMDME–Natural Gas Dual Fuel Combustion

Abstract Cyclic variations in internal combustion engines are caused by various factors, including combustion mixture stratification, in-cylinder flows, local fluctuations in air-fuel ratio, etc. Cyclic variations have a profound impact on engine performance and emissions. In this study, cyclic variations in dual fuel combustion are analyzed, comparing diesel–natural gas (NG) and polyoxymethylene dimethyl ether (POMDME)-NG dual fuel combustion. Cyclic variability was initially quantified using the coefficient of variation of gross indicated mean effective pressure (IMEPg) computed from experimental cylinder pressure data. The cases analyzed in this study had a coefficient of variation (COV) of IMEPg greater than or around 5%, which was the lower limit of onset of instability for this engine. Experiments were performed at two fixed start of injection (SOI) of high-cetane fuel: 310 CAD and 350 CAD. For all experiments, a constant load of 5 bar IMEPg was maintained, and the intake boost pressure and rail pressure were fixed at 1.5 bar and 500 bar, respectively. For each case, 1000 cycles of cylinder pressure data were recorded, filtered, and processed using an in-house heat release analysis code for each cycle. A comparison between individual cycles and the “ensemble averaged cycle” was made for both diesel–NG and POMDME–NG combustion. For the early SOI of 310 CAD, the peak cylinder pressure fluctuations of individual cycle were found to be ± 15 bar for both fuel combinations, compared to the ensemble averaged cycle, and < 1/10th of the cycles had an IMEPg lower than 0.05 bar of the ensemble averaged cycle. However, the peak pressure fluctuations were found to be lower for POMDME–NG (±3 bar) than diesel–NG dual fuel combustion at 350 CAD SOI, indicating lower cyclic variations. The higher reactivity of POMDME helped reduce fluctuations in combustion phasing at the retarded SOI. The presence of cycles of deterioration and cycles of recovery were also observed with diesel–NG combustion for 310 CAD SOI, and the scatter in the IMEPg return map was similar for both fuel combinations. The IMEPg return map for POMDME–NG combustion was less scattered at the 350 CAD SOI.

Engineering↗
DOE OSTI2022

Non-disruptive error field identification based on magnetic island healing

Here a technique to identify intrinsic error fields (EFs) in tokamaks with minimized risk of disruption is demonstrated on the DIII-D tokamak. The method extends the conventional driven magnetic island ‘compass scan’ approach by modifying asynchronous control waveforms to enable prompt healing of the island instability. Healing of the island is achieved by reducing the imposed non-axisymmetric coil current and raising the density (here via gas fueling). The method is also shown to support multiple island threshold measurements per pulse, thus reducing the number of dedicated pulses necessary to conduct an EF identification. Non-linear modeling with the TM1 code reproduces the experimental results and approximately recovers the critical density required for island healing. Island healing is explained in the non-linear modeling by an increase in the viscous coupling between the static island and the nearby flowing plasma, thus healing the island as it accelerates into the plasma frame. Due to both simplicity and risk minimization, this technique is suitable for plasma-based EF identification in the early commissioning stages of future disruption-averse tokamaks such as ITER and SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Thermal confinement and transport in spherical tokamaks: a review

Here, we review the thermal plasma confinement and transport properties observed and predicted in low aspect ratio tokamaks, or spherical tokamaks (STs), which can depart significantly from those observed at higher aspect ratio. In particular, thermal energy confinement scalings show a strong, near linear dependence of energy confinement time on toroidal magnetic field, while the dependence on plasma current is more modest, the opposite of what is seen at higher aspect ratio. STs have revealed a very strong improvement in normalized confinement with decreasing collisionality, much stronger than at higher aspect ratio, which bodes well for an ST-based fusion pilot plant should this trend continue at an even lower collisionality than has already been accessed. These differences arise because of fundamental differences in transport in STs due to the more extreme toroidicity (i.e. reduced region of bad curvature), and to the relatively larger $E_r \times B$ shearing rates, both of which can suppress electrostatic drift wave instabilities at both ion and electron gyroradius scales. In addition, electromagnetic effects are much stronger in STs because they operate at high β T . Gyrokinetic (GK) studies, coupled with low- and high-k turbulence measurements, have shed light on the underlying physics controlling transport. At lower β T , both ion- and electron-scale electrostatic drift turbulence may be responsible for transport. At higher β T , microtearing, kinetic ballooning, and hybrid trapped electron/kinetic ballooning modes increasingly play a role, and they have a much stronger impact in the core of ST plasmas than at higher aspect ratio. Flow shear affects the balance between ion- and electron-scale modes. Non-linear GK simulations find regimes where the electron heat flux decreases with decreasing collisionality, consistent with the experimental global normalized confinement scaling. The ST is unique in that the relatively low toroidal magnetic field allows for localized measurements of electron-scale turbulence, and this coupled with turbulence measurements at ion-scales has facilitated detailed comparisons with GK simulations. These data have provided compelling evidence for the presence of ion temperature gradient and electron temperature gradient turbulence in some plasmas, and direct experimental support for the impact of experimental actuators like rotation shear, density gradient and magnetic shear on turbulence and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗
DOE OSTI2020

Ring patterns generated by an expanding colloidal meniscus

The drop-and-dry is a common technique allowing for creation of periodic nanoparticle (NP) structures for sensing, photonics, catalysis, etc. However, the reproducibility and scalability of this approach for fabrication of NP-based structures faces serious challenges due to the complexity of the simple, at first glance, evaporation process. In this work we study the effect of the spatial confinement on the NP self-assembly under slow solvent evaporation, when the air-liquid-substrate contact line (CL) expands from the center towards the walls of a cylindrical cell, forming a toroid. Using in situ video monitoring of the stick-slip CL motion, we find regular hydrodynamic perturbations in the meniscus, and reveal fine details of the formation of quasiperiodic rings of close packed NP layers. We report that drying of the toroidal NP droplet has a number of important differences from drying of the classical hemispherical colloidal drops. In toroidal drops we observe linear-in-time average meniscus motion, in contrast to the hemispherical drops where the meniscus moves as a square root of time. While both droplet geometries produce NP ring patterns, the ring width for the toroidal drop decreases with increasing ring radius, while it decreases with decreasing the radius of the hemispherical drop. We suggest that free ligands are the main cause of the Marangoni instabilities driving the periodic vorticity in the meniscus. Finally, we show that the usually ignored contact line tension may yield a considerable contribution to the CL pinning causing the CL slip-stick motion and the ring formation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comprehensive Evaluation for Protective Coatings: Optical, Electrical, Photoelectrochemical, and Spectroscopic Characterizations

Numerous efficient semiconductors suffer from instability in aqueous electrolytes. Strategies utilizing protective coatings have thus been developed to protect these photoabsorbers against corrosion while synergistically improving charge separation and reaction kinetics. Recently, various photoelectrochemical (PEC) protective coatings have been reported with suitable electronic properties to ensure low charge transport loss and reveal the fundamental photoabsorber efficiency. However, protocols for studying the critical figures of merit for protective coatings have yet to be established. For this reason, we propose four criteria for evaluating the performance of a protective coating for PEC water-splitting: stability, conductivity, optical transparency, and energetic matching. We then propose a flow chart that summarizes the recommended testing protocols for quantifying these four performance metrics. In particular, we lay out the stepwise testing protocols to evaluate the energetics matching at a semiconductor/coating/(catalyst)/liquid interface. Finally, we provide an outlook for the future benchmarking needs for coatings.

36 MATERIALS SCIENCE↗
DOE OSTI2020

AC loss and contact resistance in REBCO CORC ® , Roebel, and stacked tape cables

Abstract Many high-temperature superconductor (HTS) applications require superconducting cables with high currents while operating in an alternating magnetic field. HTS cables should be composed of numerous superconducting tapes to achieve the required current capacity. Alternating current and magnetic fields cause AC losses in such cables and can provoke conductor instability. AC losses and contact resistances were measured of several cable designs based on commercially available REBCO tapes at the University of Twente. The AC loss was measured under identical conditions for eight REBCO conductors manufactured according to three types of cabling methods—CORC ® (Conductor on Round Core), Roebel, and stacked tape, including a full-size REBCO CICC (cable in conduit conductor). The measurements were done at T = 4.2 K without transport current in a sinusoidal AC magnetic field of 0.4 T amplitude and frequencies from 5 to 55 mHz. The AC loss was measured simultaneously by calibrated gas flow calorimeter utilizing the helium boil-off method and by the magnetization method using pick-up coils. Also, the AC loss of two CORC® conductors and a Roebel cable was measured at 77 K. Each conductor was measured with and without background field of 1 T. The measured AC coupling loss in the CORC ® and Roebel conductors is negligible at 4.2 K for the applied conditions while at 77 K coupling loss was observed for all conductors. The absence of coupling loss at 4.2 K can be explained by shielding of the conductor interior; this is confirmed with measurement and calculation of the penetration field of CORC ® and Roebel cables. The inter-tape contact resistance was measured for CORC ® and stacked tape samples at 4.2 and 77 K. It was demonstrated that a short heat treatment of CORC ® conductor with solder-coated tapes activates tape-to-tape soldering and decreases the contact resistance. The reduction of contact resistance by two orders in magnitude to tens of nΩm is comparable with the interstrand contact resistance in ITER Nb 3 Sn type conductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Transition Zone Morphology Dynamics of Dissolved Oxygen (DO) in a Salinity-Impacted Hyporheic Zone

Freshwater salinization is a common yet underappreciated environmental problem in rivers, yet how this process impacts the transport of dissolved oxygen (DO) in the hyporheic zone remains unknown. Here, using flume experiments and numerical simulations, we have demonstrated that the morphological dynamics of the DO transition zone are controlled by the dimensionless number Da·Ra –2 , a measure of the rate of aerobic respiration to free convection. At low Da·Ra –2 , the oxic front breaks up into fingers which subsequently grow, and leads to accelerating of DO transport and shrinking of the anaerobic zone in the underlying sediment. As Da·Ra –2 increases, the growth rate of the instability decreases, and the DO plume fingers are suppressed and delayed compared with the saltwater-freshwater interface, the mixing area between saltwater and freshwater. These results indicate that the freshwater salinization syndrome can have significant impacts on the functioning of aquifer ecosystems unless regulated and managed effectively.

54 ENVIRONMENTAL SCIENCES↗