SEARCH · Search NASA

Results for “FLUID MECHANICS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Microscopic activated dynamics theory of the shear rheology and stress overshoot in ultradense glass-forming fluids and colloidal suspensions

We formulate a particle and force level, activated dynamics-based statistical mechanical theory for the continuous startup nonlinear shear rheology of ultradense glass-forming hard sphere fluids and colloidal suspensions in the context of the elastically collective nonlinear Langevin equation approach and a generalized Maxwell model constitutive equation. Activated structural relaxation is described as a coupled local-nonlocal event involving caging and longer range collective elasticity which controls the characteristic stress relaxation time. Theoretical predictions for the deformation-induced enhancement of mobility, the onset of relaxation acceleration at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady-state structural relaxation time and viscosity, a nonvanishing activation barrier in the shear thinning regime, an apparent Herschel–Buckley form of the shear rate dependence of the steady-state shear stress, exponential growth of different measures of a yield or flow stress with packing fraction, and reduced fragility and dynamic heterogeneity under deformation were previously shown to be in good agreement with experiments. The central new question we address here is the defining feature of the transient response—the stress overshoot. In contrast to the steady-state flow regime, understanding the transient response requires an explicit treatment of the coupled nonequilibrium evolution of structure, elastic modulus, and stress relaxation time. We formulate a new quantitative model for this aspect in a physically motivated and computationally tractable manner. Theoretical predictions for the stress overshoot are shown to be in good agreement with experimental observations in the metastable ultradense regime of hard sphere colloidal suspensions as a function of shear rate and packing fraction, and accounting for deformation-assisted activated motion appears to be crucial for both the transient and steady-state responses.

Mechanics↗

DOE OSTI2024

Environmental damping and vibrational coupling of confined fluids within isolated carbon nanotubes

Abstract Because of their large surface areas, nanotubes and nanowires demonstrate exquisite mechanical coupling to their surroundings, promising advanced sensors and nanomechanical devices. However, this environmental sensitivity has resulted in several ambiguous observations of vibrational coupling across various experiments. Herein, we demonstrate a temperature-dependent Radial Breathing Mode (RBM) frequency in free-standing, electron-diffraction-assigned Double-Walled Carbon Nanotubes (DWNTs) that shows an unexpected and thermally reversible frequency downshift of 10 to 15%, for systems isolated in vacuum. An analysis based on a harmonic oscillator model assigns the distinctive frequency cusp, produced over 93 scans of 3 distinct DWNTs, along with the hyperbolic trajectory, to a reversible increase in damping from graphitic ribbons on the exterior surface. Strain-dependent coupling from self-tensioned, suspended DWNTs maintains the ratio of spring-to-damping frequencies, producing a stable saturation of RBM in the low-tension limit. In contrast, when the interior of DWNTs is subjected to a water-filling process, the RBM thermal trajectory is altered to that of a Langmuir isobar and elliptical trajectories, allowing measurement of the enthalpy of confined fluid phase change. These mechanisms and quantitative theory provide new insights into the environmental coupling of nanomechanical systems and the implications for devices and nanofluidic conduits.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Hydrodynamic computations of high-power laser drives generating metal ejecta jets from surface grooves

Understanding dynamic fragmentation in shock-loaded metals and predicting properties of the resulting ejecta are of considerable importance for both basic and applied science. The nature of material ejection has been shown to change drastically when the free surface melts on compression or release. In this work, we present hydrodynamic simulations of laser-driven microjetting from micron-scale grooves on a tin surface. We study microjet formation across a range of shock strengths from drives that leave the target solid after release to drives that induce shock melting in the target. The shock-state particle velocity (Up) varies from 0.3 to 3 km/s and the shock breakout pressure is 3–120 GPa. The microjet tip velocity is 1–8 km/s and the free-surface velocity varies from 0.1 to 5 km/s. Two tin equations of state are examined: a “soft” model (LEOS 501) where the target melts for U p > 1 km/s and a more detailed multiphase model (SESAME 2161) that melts for U p > 1.4 km/s. We use these two models to examine the influence of phase change and the choice of the material model on microjet formation and evolution. We observe in our computational results that jet formation can be classified into three regimes: a low-energy regime where material strength affects jet formation, a moderate-energy regime dominated by the changing phase of tin material, and a high-energy regime where results are insensitive to the material model and jet formation is described by an idealized steady-jet theory. Using an ensemble of 2D simulations, we show that these trends hold across a wide range of drive energies and groove angles.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

FY23 Status Report on MOOSE-Based Approaches to Modeling Core Bowing in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY23, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21 and FY22. Several benchmark verification examples were selected from the IAEA verification and validation report with increased number of ducts and more complex contact interaction behavior. These examples involve a full symmetric sector with restraint rings at multiple load pad locations to simulate a limited free-bow restraint system concept, as well as irradiation induced swelling and creep effects in a sector. In addition, improvements to the contact module sideset assignment were assessed and compared with previous MOOSE results to verify the contact behavior. The MOOSE-based results were compared to IAEA benchmark participants’ results. In every case, the MOOSE results agreed with other simulations results for estimating bowing behavior, providing additional verification basis of the tools for this particular physics application. Estimation of contact forces were mostly in agreement, with a few outlier results. A plan was suggested for dealing with the discrepancies with estimating contact force values. In addition, a 1-way coupling demonstration was performed using subchannel analysis code Pronghorn-SC and MOOSE on an ABR-1000-design sodium-cooled fast reactor assembly to evaluate coolant and duct wall temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Numerical calculation of the particle–fluid–particle stress in random arrays of fixed particles

Based on the nearest particle statistics [Zhang, J. Fluid Mech. 910, A16 (2021)], the phase interaction force in a multiphase flow is decomposed into a particle–mean-field force and the divergence of the particle–fluid–particle (PFP) stress. The PFP stress is proportional to the correlation product of the distance from a particle to its nearest neighbor and the force on the particle conditionally averaged on the nearest-neighbor location. In this work, a functional form of the stress is obtained corrected to the first order of the ratio between the interparticle distance and the macroscopic length scale of the flow. Particle-resolved numerical simulations are used to calculate the PFP stress in random arrays of fixed particles and to explore the physics represented by the stress. The numerical results show that the PFP stress is attractive along the direction of the flow and is repulsive in the directions perpendicular to the flow. In the flow regime simulated, this PFP stress can be considered as a macroscopic representation of the drafting–kissing–tumbling mechanism. Finally, the Reynolds stress for the fluid phase is also calculated and compared with the PFP stress.

42 ENGINEERING↗

DOE OSTI2020

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Investigating the Effects of Chemical Mechanism on Soot Formation Under High-Pressure Fuel Pyrolysis

We performed Computational Fluid Dynamics (CFD) simulations using a Reynolds-Averaged Navier-Stokes (RANS) turbulence model of high-pressure spray pyrolysis with a detailed chemical kinetic mechanism encompassing pyrolysis of n-dodecane and formation of polycyclic aromatic hydrocarbons. We compare the results using the detailed mechanism and those found using several different reduced chemical mechanisms to experiments carried out in an optically accessible, high-pressure, constant-volume combustion chamber. Three different soot models implemented in the CONVERGE CFD software are used: an empirical soot model, a method of moments, and a discrete sectional method. There is a large variation in the prediction of the soot between different combinations of chemical mechanisms and soot model. Furthermore, the amount of soot produced from all models is substantially less than experimental measurements. All of this indicates that there is still substantial work that needs to be done to arrive at simulations that can be relied on to accurately predict soot formation.

02 PETROLEUM↗

DOE OSTI2020

Reaction heterogeneity in practical high-energy lithium–sulfur pouch cells

The lithium–sulfur (Li–S) battery is a promising next-generation energy storage technology because of its high theoretical energy and low cost. Extensive research efforts have been made on new materials and advanced characterization techniques for mechanistic studies. However, it is uncertain how discoveries made on the material level apply to realistic batteries due to limited analysis and characterization of real high-energy cells, such as pouch cells. Evaluation of pouch cells (>1 A h) (instead of coin cells) that are scalable to practical cells provides a critical understanding of current limitations which enables the proposal of strategies and solutions for further performance improvement. Herein, we design and fabricate pouch cells over 300 W h kg -1 , compare the cell parameters required for high-energy pouch cells, and investigate the reaction processes and their correlation to cell cycling behavior and failure mechanisms. Spatially resolved characterization techniques and fluid-flow simulation reveal the impacts of the liquid electrolyte diffusion within the pouch cells. We found that catastrophic failure of high-energy Li–S pouch cells results from uneven sulfur/polysulfide reactions and electrolyte depletion for the first tens of cycles, rather than sulfur dissolution as commonly reported in the literature. The uneven reaction stems from limited electrolyte diffusion through the porous channels into the central part of thick cathodes during cycling, which is amplified both across the sulfur electrodes and within the same electrode plane. A combination of strategies is suggested to increase sulfur utilization, improve nanoarchitectures for electrolyte diffusion and reduce consumption of the electrolytes and additives.

25 ENERGY STORAGE↗

DOE OSTI2025

Damage modeling of power tower receiver tubes using the SRLIFE tool

Concentrating Solar Power (CSP) molten-salt central receivers are subject to high, transient incident flux during daily operation. The resulting creep-fatigue damage impacts the receiver’s reliability and restricts the permissible incident flux distribution for a given receiver. This paper aims to reduce CSP plants’ levelized cost of electricity by developing a methodology to predict lifetime and identifies the primary damage mechanism (creep vs fatigue) for any given fluid temperature and temperature gradient. Results are presented in the form of a damage map that serves as a valuable operation guide and design tool. Damage maps can be used to reduce maintenance costs by improving reliability and reduce receiver capital costs by better utilizing the receiver area. FEA simulation and damage modeling of tubes subject to asymmetrical flux conditions is performed in the open-source receiver design tool srlife. Parametric studies are performed over a range of inner tube temperatures and thermal gradients for A230, 316H, 740H, A282, A617, and 800H high temperature alloys. Damage maps are presented for each alloy. A parametric, FEA-based methodology is presented for comparison of fatigue-creep ratios and prediction of tube lifetime based on the critical thermal operating conditions. Fatigue is found to be negligible compared to creep for almost every case. Here, this finding suggests that fatigue effects associated with cloud events are insignificant compared to creep at these high temperature operating conditions. Additionally, lifetime predictions identify thermal conditions where small changes in operating conditions can result in large changes in predicted lifetime.

14 SOLAR ENERGY↗

DOE OSTI2022

Iron corrosion in the “inert” supercritical CO 2 , ab initio dynamics insights: How impurities matter

Metal corrosion from exposure to supercritical fluids is a long-standing challenge lacking atomistic understanding for effective corrosion prevention in diverse important applications. Here, in this work, reactions relevant to corrosion in supercritical CO 2 (sCO 2 ), with H 2 O and NO 2 as impurities, are revealed using ab initio molecular dynamics and enhanced sampling methods. While little reactivity is observed in bulk sCO 2 , the iron surface is found to play a crucial almost-catalytic role in its own corrosion, activating the inert CO2 and hence facilitating $\ast$HCO 3 – and $\ast$CO 3 2– formation at the interface. Moreover, the hydrogen bond network of $\ast$H 2 O assists the H shuttling to the activated CO 2 and impurities, forming reactive species observed in experiments. Such atomistic insight elucidates the reaction mechanisms of metal corrosion in aqueous and supercritical fluids and can point to possible remedies.

36 MATERIALS SCIENCE↗

DOE OSTI2020

A viscoelastic model for seismic attenuation using fractal mechanical networks

SUMMARY Seismic attenuation (quantified by the quality factor Q) has a significant impact on the seismic waveforms, especially in the fluid-saturated rocks. This dissipative process can be phenomenologically represented by viscoelastic models. Previous seismological studies show that the Q value of Earth media exhibits a nearly frequency-independent behaviour (often referred to as constant-Q in literature) in the seismic frequency range. Such attenuation can be described by the mathematical Kjartansson constant-Q model, which lacks of a physical representation in the viscoelastic sense. Inspired by the fractal nature of the pore fluid distribution in patchy-saturated rocks, here we propose two fractal mechanical network (FMN) models, that is, a fractal tree model and a quasi-fractal ladder model, to phenomenologically represent the frequency-independent Q behaviour. As with the classic viscoelastic models, the FMN models are composed of mechanical elements (spring and dashpots) arranged in different hierarchical patterns. A particular parametrization of each model can produce the same complex modulus as in the Kjartansson model, which leads to the constant-Q. Applying the theory to several typical rock samples, we find that the seismic attenuation signature of these rocks can be accurately represented by either one of the FMN models. Besides, we demonstrate that the ladder model in particular exhibits the realistic multiscale fractal structure of the saturated rocks. Therefore, the FMN models as a proxy could provide a new way to estimate the microscopic rock structure property from macroscopic seismic attenuation observation.

Geochemistry & Geophysics↗

DOE OSTI2021

Spontaneous imbibition in a square tube with corner films: theoretical model and numerical simulation

Spontaneous imbibition in an angular tube with corner films is a fundamental problem in many scientific and engineering processes. In this paper, a modified interacting capillary bundle model is developed to describe the liquid imbibition dynamics in a square tube with corner films. The square tube is decomposed into several interacting sub-capillaries and the local capillary pressure in each sub-capillary is derived based on the specific shape of its meniscus. The conductance of each sub-capillary is calculated using single-phase lattice Boltzmann simulation. The modified interacting capillary bundle model and color-gradient lattice Boltzmann method are used to simulate the liquid imbibition dynamics in the square tube with different fluid properties. The predictions by the modified interacting capillary bundle model match well with the lattice Boltzmann simulation results for different conditions, demonstrating the accuracy and robustness of the interacting capillary bundle model to describe the imbibition dynamics with corner films. In addition, the interacting capillary bundle model is helpful to investigate the mechanisms during spontaneous imbibition and the influences of fluid viscosity, surface tension, wetting phase contact angle and gravity on imbibition dynamics. Finally, a universal scaling law of imbibition dynamics for the main meniscus is developed and the scaling law for arc meniscus is also analyzed.

42 ENGINEERING↗

DOE OSTI2024

Implementation of Detailed Polyethylene Pyrolysis Kinetics into CFD Simulations using Machine Learning

Municipal solid waste (MSW) and waste plastics have received significant attention due to the issues of waste generation and storage, as well as their potential as an energy resource. High-density polyethylene (HDPE) makes up a large portion of plastic waste and has been the subject of several conversion studies. However, the mechanisms associated with converting HDPE through pyrolysis and gasification are extensive and complex making them difficult to implement into high-fidelity computational fluid dynamic (CFD) simulations. For this project, a primary pyrolysis mechanism containing 42 unique species and 737 heterogeneous reactions was used to generate kinetic data over a range of operating conditions. A machine learning (ML) model was developed to replicate the results of the detailed pyrolysis mechanism while significantly increasing the computational efficiency. A deep operator network (DeepONet) architecture was adopted to train the model using time steps relevant to CFD simulations. The ML used physics-based loss functions to ensure mass conservation. The ML model has been deployed in simple MFiX CFD simulations, single particle, and an experimental drop tube reactor, and has shown promising performance compared to the original scheme.

Houston, Ross↗

DOE OSTI2022

Facile polymerization in a bicellar template to produce polymer nano-rings

We report a well-defined discoidal bicelle composed of three lipids, specifically zwitterionic long-chain 1,2-dipalmitoyl phosphocholine (DPPC) and short-chain 1,2-dihexanoyl phosphocholine (DHPC) doped with anionic 1,2-dipalmitoyl phosphoglycerol (DPPG) provides a generalized template for the synthesis of hydrophobic polymer nano-rings. The lipid molar ratio of DPPC/DHPC/DPPG is 0.71/0.25/0.04. The detailed investigation and discussion were based on styrene but tested on three other vinyl monomers. The structure of nano-rings is identified through the detailed analysis of small angle X-ray/neutron scattering (SAXS and SANS) data and transmission electron micrographs (TEM), supported by the differential scanning calorimetric (DSC) data before and after polymerization. The investigation covers samples with a styrene-to-lipid ratio ranged varied from 1:50 to 1:10.The styrene monomers are initially located at both the discoidal planar (long-chain lipid rich) and rim (short-chain lipid rich) regions. During polymerization, they migrate to the more fluid rim region-section. The formation mechanism involves the interplay of hydrophobic interaction, mismatched miscibility of polystyrene between the ordered and disordered phases, and crystallinity of the long lipid acyl chains. This facile synthesis is proven applicable for several hydrophobic monomers. The well-defined nano-rings greatly enhance the interfacial area and have the potential to be the building blocks for functional materials, if monomers are incorporated with desirable functions, for future applications. Crosslinked polymer nano-rings are produced in well-defined discoidal bicelles. The enhanced immiscibility in the planar region drives the hydrophobic polymers to the fluidic rim. The platform can be generalized for synthesis of other hydrophobic polymer nano-rings.

36 MATERIALS SCIENCE↗

DOE OSTI2022

The role of surface forces in environment-enhanced cracking of brittle solids

Fracture initiation and propagation in brittle materials is promoted in surface-reactive (sorptive) environments, a phenomenon known as subcritical crack growth (SCG). Laboratory measured crack-propagation velocity vs. stress intensity factor relationships typically exhibit highly nonlinear, multi-stage characteristics that are sensitive to environmental factors such as adsorbate concentration and temperature. For practical purposes, empirical relationships (e.g., a power law) have been used to describe this complex phenomenon. However, how the overall SCG behavior emerges from the underlying fundamental processes near the crack tip, such as the interaction of the crack surfaces separated by only a few nanometers and mass transport within the nano-confined space, is still not well understood. This paper develops a mechanistic, surface-force-based fracture theory (SFFT) which integrates surface force models, fluid transport models, and linear elastic fracture mechanics to quantitatively explain the multi-stage characteristics of SCG in brittle solids. A numerical model is developed based on SFFT and solved through an implicit partitioned scheme for efficiency and modularity. The results are validated by Wiederhorn's data on crack propagation in soda-lime glasses at a wide range of relative humidity levels. We show that, for the first time, the entire range of an SCG curve can be captured by a single physics-based model. The predicted SCG curves reveal that the development of repulsive disjoining pressure behind the crack tip can be responsible for the reduced apparent fracture toughness in a sorptive environment. In conclusion, the shape of the SCG curve, and its changes with respect to the environment, is found to critically depend on the assumed transport models.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Facile Metal Release from Pore-Lining Phases Enables Unique Carbonate Zonation in a Basalt Carbon Mineralization Demonstration

Carbon-negative strategies such as geologic carbon sequestration in continental flood basalts offers a promising route to the removal of greenhouse gases, such as CO 2 , via safe and permanent storage as stable carbonates. This potential has been successfully demonstrated at a field scale at the Wallula Basalt Carbon Storage Pilot Project where supercritical CO 2 was injected into the Columbia River Basalt Group (CRBG). Here, we analyze recovered post-injection sidewall core cross-sections containing carbonate nodules using μ-XRF chemical mapping techniques that revealed compositional zonation within the nodules. The unique nature of the subsurface anthropogenic carbonates is highlighted by the near absence of Mg in an ankerite-like composition. Furthermore, a comparison between pre- and post-injection sidewall cores along with an in-depth chemical mapping of basalt pore lining cements provides a better understanding into the source and fate of critical cationic species involved in the precipitation of carbon mineralization products. Collectively, these results provide crucial insights into carbonate growth mechanisms under a time-dependent pore fluid composition. Finally, as such, these findings will enable parameterization of predictive models for future CO 2 sequestration efforts in reactive reservoirs around the world.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Statistical properties of homogeneous and isotropic turbulence in He II measured via particle tracking velocimetry

Despite being a quantum two-fluid system, superfluid helium-4 (He II) is observed to behave similarly to classical fluids when a flow is generated by mechanical forcing. This similarity has brought up the feasibility of utilizing He II for high Reynolds number classical turbulence research, considering the small kinematic viscosity of He II. However, it has been suggested that the nonclassical dissipation mechanism in He II at small scales may alter its turbulent statistics and intermittency. In this work, we report our study of a nearly homogeneous and isotropic turbulence (HIT) generated by a towed grid in He II. We measure the velocity field using particle tracking velocimetry with solidified deuterium particles as the tracers. By correlating the velocities measured simultaneously on different particle trajectories or at different times along the same particle trajectory, we are able to conduct both Eulerian and Lagrangian flow analyses. Spatial velocity structure functions obtained through the Eulerian analysis show scaling behaviors in the inertial subrange similar to that for classical HIT but with enhanced intermittency. The Lagrangian analysis allows us to examine the flow statistics down to below the dissipation length scale. Interestingly, abnormal deviations from the classical scaling behaviors are observed in this regime. Lastly, we discuss how these deviations may relate to the motion of quantized vortices in the superfluid component in He II.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗