Search NASA⌕ Search

SEARCH · Search NASA

Results for “fluid mixing”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗
DOE OSTI2024

Simulations of three-layer Richtmyer–Meshkov mixing in a shock tube

Here, the Richtmyer–Meshkov instability causes perturbations to grow after a shock traverses a fluid density interface. This increases the mixing rate between fluid from either side of the interface. We use the Flash Eulerian hydrodynamic code to investigate alterations when a thin third layer of intermediate density is placed along the interface, effectively creating two adjacent unstable interfaces. This is a common occurrence in engineering applications where a thin barrier initially separates two materials. We find that the width of the mixing layer is similar or slightly reduced; however, the total mass of mixed material can actually increase. The mixing layer becomes more compact and efficient. However, the normalized mixed mass decreases, meaning that finger entrainment becomes more important than in the simple two-layer case. The effect of adding the central layer appears to decrease when the Atwood number is decreased. The Flash results are also benchmarked against two-layer experimental data from a shock tube at the University of Arizona.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impacts of Pore Network-Scale Wettability Heterogeneity on Immiscible Fluid Displacement: A Micromodel Study

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, here we show the classic fingering and transitions in flow regimes at $logCa$ (capillary number) from -7.2 to -3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO 2 (scCO 2 ) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO 2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

2.5-D micromodel↗
DOE OSTI2024

Numerical Simulation and Experimental Comparison of System Analysis Module 1D Mixing Model for Cold Shock Transients in the Gallium Thermal-Hydraulic Mixing Facility

Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.

stratification↗
DOE OSTI2024

High-Fidelity Simulation of the Light-to-Dense Stratification Transient in the HiRJET Facility

Density stratification in a large enclosure is a crucial phenomenon to heat transfer and sustainable passive heat removal of a sodium fast reactor during reactivity transients. However, engineering turbulence models were identified to have unsatisfactory performance in predicting propagation of a stratified front. Yet, the scarcity of high-resolution data for stratification hampers the development of models. To explor e applications of leveraging direct numerical simulation (DNS) data to support turbulence model development, this work conducted DNS using NekRS to study a long stratification transient in the High-Resolution Jet (HiRJET) experimental facility. Here, this work considers an experiment run where light fluid is injected into a tank containing a denser fluid with a relative density difference of 1.5%. Formation of the stratified layer is identified as impingement of the buoyant jet promoting mixing of the two fluids. Based on the transient statistics, transport of the concentration can be characterized by regions with dominating effects of turbulent mixing, buoyant dissipation, and molecular diffusion, respectively, as moving away from the elevation of jet impingement. Concentration near the stratified front also exhibits oscillation at Brunt-Väisälä frequency. Preliminary validation of the simulation showed encouraging agreement of the concentration distribution with the reference experiment.

buoyant mixing↗
DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗
DOE OSTI2021

Kinetic mix at gas-shell interface in inverted corona fusion targets

Gas-filled, laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. At the fill pressures under investigation, ejected particles from the shell can penetrate deeply into the gas before colliding, leading to significant mixing across the gas–shell interface. Here, we use kinetic-ion, fluid-electron hybrid particle-in-cell simulations to explore the nature of that mix. Simulations of the system demonstrate characteristics of a weakly collisional electrostatic shock, whereby a strong electric field accelerates shell ions into the rarefied gas and reflects upstream gas ions. This interpenetration is mediated by collisional processes: At higher initial gas pressure, fewer shell particles pass into the mix region and reach the hotspot. This effect is detectable through neutron yield scaling vs gas pressure. Predictions of neutron yield scaling show excellent agreement with experimental data recorded at the OMEGA laser facility, suggesting that 1D kinetic mechanisms are sufficient to capture the mix process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Feasibility of Lead Fast Reactor Heat Exchanger Tube Online Monitoring

Within virtually all reactor systems with heat exchangers, the thinnest barriers between primary and secondary fluids are the heat exchanger tubes. Consistent with light water reactor experience, the advanced reactor heat exchangers are the most susceptible regions for corrosion and leakage (NRC, 2019). Moreover, traditional inspections and repairs using robotics will be expensive or impractical. On-line structural health corrosion monitoring (OLSHM) is therefore desirable to detect tube degradation prior to leaks that may allow mixing of heat exchanger fluids. Rhetorically many molten salt reactor (MSR) designers plan for on-line monitoring of critical reactor components including heat exchangers but there are no proven OLSHM sensors or systems for MSR heat exchangers. Guided wave ultrasound from the tube-end or ID using conventional low-temperature piezoelectric sensors has been shown to be sensitive to pits and cracks less than 50% through-wall in long pipes and tubes (Narayanan M.M., 2019) . Mode and frequency selection are important for practical sensor design decisions and ultimate performance. The L-0 wave-mode is easily excited by L-mode transducers but they have surface normal wave motion that will be sensitive to and more attenuative because of fluid interaction with the tube surface (Rose J.L., 1994). Shear Horizontal (SH) wave modes have material motion parallel to the tube surface therefore are virtually unaffected by fluid interaction. For SH-0 to be nondispersive in steel, the frequency-thickness product (MHz-mm) must be less than approximately 1.5. SH mode transducers must be coupled by a material that can support shear stress; a thin layer of honey works well for temporary sensors and an adhesive or a brazed metallurgical bond for permanent sensors. If sensors are incorporated into the fabrication process, the same kinds of SH-0 waves can be generated from the stagnant flow area of the tube OD without compromising the heat-exchange tube function. Piezoelectric materials are also available that can withstand harsh temperatures of molten salt heat exchangers (Stevenson T., 2015) (Machura et al., 2008). Though they may not have as high-performance characteristics, their performance only needs to be good enough to detect flaws of interest. This project attempts to demonstrate feasibility to develop a sensor that can be sensitive to flaws of interest and perform long-term monitoring of a molten salt reactor heat exchanger. Benefits to the advanced reactor system include: • Minimizing the expense of a reactor shut down for periodic inspection—only shutdown for cause. • Reduce risk of unexpected tube rupture by detecting flaws before they grow to through-wall leaks. • Cost reduction by minimizing need to design for manual/robotic inspectability. Transducer design, characterization of the beam profile, wedges vs. direct application of piezoelectric sensors, and ultrasonic pulse-echo responses from a representative heat exchange tube with fabricated flaw types are discussed.

On-Line Tube Monitor, Guided Wave Ultrasound, SH-0↗
DOE OSTI2020

Progress Toward Molten Salt Reactor Heat Exchanger On-Line Monitoring

Within virtually all reactor systems with heat exchangers, the thinnest barriers between primary and secondary fluids are the heat exchanger tubes. Consistent with light water reactor experience, the advanced reactor heat exchangers are the most susceptible regions for corrosion and leakage (NRC, 2019). Moreover, traditional inspections and repairs using robotics will be expensive or impractical. On-line structural health corrosion monitoring (OLSHM) is therefore desirable to detect tube degradation prior to leaks that may allow mixing of heat exchanger fluids. Rhetorically many molten salt reactor (MSR) designers plan for on-line monitoring of critical reactor components including heat exchangers but there are no proven OLSHM sensors or systems for MSR heat exchangers. Guided wave ultrasound from the tube-end or ID using conventional low-temperature piezoelectric sensors has been shown to be sensitive to pits and cracks less than 50% through-wall in long pipes and tubes (Narayanan M.M., 2019) . Mode and frequency selection are important for practical sensor design decisions and ultimate performance. The L-0 wave-mode is easily excited by L-mode transducers but they have surface normal wave motion that will be sensitive to and more attenuative because of fluid interaction with the tube surface (Rose J.L., 1994). Shear Horizontal (SH) wave modes have material motion parallel to the tube surface therefore are virtually unaffected by fluid interaction. For SH-0 to be nondispersive in steel, the frequency-thickness product (MHz-mm) must be less than approximately 1.5. SH mode transducers must be coupled by a material that can support shear stress; a thin layer of honey works well for temporary sensors and an adhesive or a brazed metallurgical bond for permanent sensors. If sensors are incorporated into the fabrication process, the same kinds of SH-0 waves can be generated from the stagnant flow area of the tube OD without compromising the heat-exchange tube function. Piezoelectric materials are also available that can withstand harsh temperatures of molten salt heat exchangers (Stevenson T., 2015) (Machura et al., 2008). Though they may not have as high-performance characteristics, their performance only needs to be good enough to detect flaws of interest. This project attempts to demonstrate feasibility to develop a sensor that can be sensitive to flaws of interest and perform long-term monitoring of a molten salt reactor heat exchanger. Benefits to the advanced reactor system include: • Minimizing the expense of a reactor shut down for periodic inspection—only shutdown for cause. • Reduce risk of unexpected tube rupture by detecting flaws before they grow to through-wall leaks. • Cost reduction by minimizing need to design for manual/robotic inspectability. Transducer design, characterization of the beam profile, and ultrasonic pulse-echo responses from a representative heat exchange tube with fabricated flaw types are discussed.

On-Line Tube Monitor, Guided Wave Ultrasound, SH-0↗
DOE OSTI2020

Effect of the numerical discretization scheme in Shock-Driven turbulent mixing simulations

In this work, we evaluate the effects of distinct numerical strategies in simulations of shock-driven turbulent mixing. An air-SF 6 -air gas curtain subjected to Mach 1.2 shock-waves is computed with Implicit Large-Eddy Simulation based on three numerical schemes: directional-split with Harten-Lax-van Leer (HLL) solver; directional-unsplit using HLL-Contact (HLLC) solver; and directional-unsplit with HLLC solver and a Low Mach number Correction (LMC). The results illustrate the importance of the numerical strategy to the accuracy of the predictions. Whereas both split and unsplit schemes result in a similar spatial development of the initial shock-driven instability, only the unsplit schemes can predict the turbulent mixing transition after reshock observed by the laboratory experiments. Such feature increases the mixing rate of the two fluids, this being particularly pronounced when the LMC is active due to i) the reduced flow Mach number, and ii) the larger effective Reynolds number. Since the selected mixing problem is driven by the deposition of vorticity at the fluids’ interface, the resultant flow physics is analyzed by investigating the contribution of distinct inviscid mechanisms to the production of vorticity: baroclinicity, stretching, and dilatation. As expected, it is observed that the production of vorticity is initially dominated by the baroclinicity mechanism. Yet, the relevance of the remainder mechanisms is enhanced after the reshock and may even surpass the baroclinicity term.

42 ENGINEERING↗
DOE OSTI2021

Fully consistent Eulerian Monte Carlo fields method for solving probability density function transport equations in turbulence modeling

The probability density function (PDF) transport equation method is a sophisticated model for the closure of turbulent mixing and turbulent reactive flows. An efficient solution approach for solving the PDF transport equation has been vital for the method to be widely used in applications. The Eulerian Monte Carlo fields (EMCF) method has been developed to solve the PDF transport equation efficiently for decades. A recent work by Wang et al. revealed a serious issue of the EMCF method for not being fully consistent with the PDF transport equation for which the method is designed to solve. This work advances the state of the art by introducing fully consistent EMCF methods for solving the PDF transport equation. The fully consistent EMCF formulations are derived for two different PDF equation forms. The consistency of the EMCF formulations is mathematically confirmed by examining the derived moment transport equations from the EMCF formulations and from the PDF transport equation. The method of manufactured solutions is employed to further verify the consistency and convergence of the different EMCF formulations numerically. The newly introduced EMCF formulations bring the EMCF method to full consistency with the PDF transport equations for the first time.

42 ENGINEERING↗
DOE OSTI2021

Experimental study on rheological and settling properties of shape memory polymer for fracture sealing in geothermal formations

This article studies the rheology and annular flow of a smart lost circulation material (LCM) at 160 °C, using a high-temperature flow loop. The smart LCM is a shape memory polymer (SMP) that activates at high temperatures, and its particle size increases to seal fractures in geothermal formations. SMP was mixed with the base fluid in two different concentrations, 1.0 and 3.0 wt%, to study rheology, wellbore hydraulics, activation process, and settling behavior under different testing conditions. The results of this study showed that the SMP could be activated at high temperatures. An increase of 80–100% in the particle size was observed at 160 °C. The mud samples showed a high shear-thinning behavior at the two concentrations with a Power-law flow index (n) ranging between 0.025 and 0.101. No additional frictional pressure losses were observed when SMP was added to the base fluid. SMP particles showed an excellent suspension at 1.0 wt% while, at 3.0 wt%, a bed was formed at a low flow rate and without pipe rotation. Increasing the drill pipe rotational speed or flow rate effectively removed the bed and homogeneously dispersed the SMP particles, ensuring a better sealing efficiency. SMP particle dispersion in inclined wells was better than in horizontal wells. Moreover, the findings of this study help optimize the lost circulation treatment by considering a wide range of operating parameters that can further be extended to different systems and geometries.

02 PETROLEUM↗
DOE OSTI2024

Prandtl number effects on extreme mixing events in forced stratified turbulence

Relatively strongly stratified turbulent flows tend to self-organise into a ‘layered anisotropic stratified turbulence’ (LAST) regime, characterised by relatively deep and well-mixed density ‘layers’ separated by relatively thin ‘interfaces’ of enhanced density gradient. Understanding the associated mixing dynamics is a central problem in geophysical fluid dynamics. It is challenging to study LAST mixing, as it is associated with Reynolds numbers $Re := UL/\nu \gg 1$ and Froude numbers $Fr :=(2{\rm \pi} U)/(L N) \ll 1$ ( $U$ and $L$ being characteristic velocity and length scales, $\nu$ the kinematic viscosity and $N$ the buoyancy frequency). Since a sufficiently large dynamic range (largely) unaffected by stratification and viscosity is required, it is also necessary for the buoyancy Reynolds number $Re_{b} := \epsilon /(\nu N^{2}) \gg 1$ , where $\epsilon$ is the (appropriately volume-averaged) turbulent kinetic energy dissipation rate. This requirement is exacerbated for oceanically relevant flows, as the Prandtl number $Pr := \nu /\kappa = {O}(10)$ in thermally stratified water (where $\kappa$ is the thermal diffusivity), thus leading (potentially) to even finer density field structures. We report here on four forced fully resolved direct numerical simulations of stratified turbulence at various Froude ( $Fr=0.5, 2$ ) and Prandtl ( $Pr=1, 7$ ) numbers forced so that $Re_{b}=50$ , with resolutions up to $30\,240 \times 30\,240 \times 3780$ . We find that, as $Pr$ increases, emergent ‘interfaces’ become finer and their contribution to bulk mixing characteristics decreases at the expense of the small-scale density structures populating the well-mixed ‘layers’. However, extreme mixing events (as quantified by significantly elevated local destruction rates of buoyancy variance $\chi _0$ ) are always preferentially found in the (statically stable) interfaces, irrespective of the value of $Pr$ .

Mechanics↗
DOE OSTI2022

Synthesis and structure of carbon-doped H 3 S compounds at high pressure

Understanding of recently reported putative close-to-room-temperature superconductivity in C-S-H compounds at 267 GPa demands a reproducible synthesis protocol as well as knowledge of the compounds' structure and composition. Here, we synthesized C-S-H compounds with various carbon compositions at high pressures from elemental carbon C and methane CH 4 , sulfur S, and molecular hydrogen H 2 . Here, we focus on compounds synthesized using methane as these allow a straightforward determination of their structure and composition by combining single-crystal x-ray diffraction and Raman spectroscopy. We applied a two-stage synthesis of [(CH 4 ) x (H 2 S) (1-x) ] 2 H 2 compounds with various compositions by first reacting sulfur and mixed methane-hydrogen fluids and forming CH 4 -doped H 2 S crystals at 0.5-3 GPa and then by growing single crystals of the desired hydrogen-rich compound. Raman spectroscopy applied to this material shows the presence of CH 4 molecules incorporated into the lattice and allows the determination of the CH 4 content, while single-crystal x-ray diffraction results suggest that the methane molecules substitute H 2 S molecules. The structural behavior of these compounds is very similar to the previously investigated methane-free crystals demonstrating a transition from Al 2 Cu type I4/mcm molecular crystal to a modulated molecular structure at 20-30 GPa and back to the same basic I4/mcm structure in an extended modification with greatly modified Raman spectra. This latter phase demonstrates a distortion into a Pnma structure at 132-159 GPa and then transforms into a common Im$\bar{3}$m H 3 S phase at higher pressures; however, no structural anomaly is detected near 220 GPa, where a sharp upturn in T c has been reported.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Data assimilation empowered neural network parametrizations for subgrid processes in geophysical flows

In the past couple of years, there has been a proliferation in the use of machine learning approaches to represent subgrid-scale processes in geophysical flows with an aim to improve the forecasting capability and to accelerate numerical simulations of these flows. Despite its success for different types of flow, the online deployment of a data-driven closure model can cause instabilities and biases in modeling the overall effect of subgrid-scale processes, which in turn leads to inaccurate prediction. To tackle this issue, we exploit the data assimilation technique to correct the physics-based model coupled with the neural network as a surrogate for unresolved flow dynamics in multiscale systems. In particular, we use a set of neural network architectures to learn the correlation between resolved flow variables and the parametrizations of unresolved flow dynamics and formulate a data assimilation approach to correct the hybrid model during their online deployment. We illustrate our framework in a set of applications of the multiscale Lorenz 96 system for which the parametrization model for unresolved scales is exactly known, and the two-dimensional Kraichnan turbulence system for which the parametrization model for unresolved scales is not known a priori. Our analysis, therefore, comprises a predictive dynamical core empowered by (i) a data-driven closure model for subgrid-scale processes, (ii) a data assimilation approach for forecast error correction, and (iii) both data-driven closure and data assimilation procedures. We show significant improvement in the long-term prediction of the underlying chaotic dynamics with our framework compared to using only neural network parametrizations for future prediction. Moreover, we demonstrate that these data-driven parametrization models can handle the non-Gaussian statistics of subgrid-scale processes, and effectively improve the accuracy of outer data assimilation workflow loops in a modular nonintrusive way.

42 ENGINEERING↗
DOE OSTI2022

Online Lead/Water Heat Exchanger Sensor/System Feasibility-PNNL (Final Report)

On-line structural health corrosion monitoring in advanced lead fast reactor heat exchangers and molten salt reactor heat exchangers is desirable for detecting tube degradation prior to leaks that may allow mixing of heat exchanger fluids or release of radiological contamination beyond the design containment boundary. This program, On-Line Lead/Water Heat Exchanger Sensor/System Feasibility – PNNL 76092, demonstrates feasibility for an ultrasonic torsional wave mode sensor attached to the outside of a long (30.5-m) heat exchanger tube in the stagnant flow area where the tube joins a heat exchanger plenum and where it is possible to protect a sensor and cable from high-force flows. The sensor must be cable connected to a monitoring instrument near the heat exchanger. The sensor and cable management approach for periodic in-service inspection will be impractical to implement on existing heat exchangers; rather permanently mounted sensors must be installed in conjunction with heat exchanger fabrication. Previous work has shown low-temperature lead zirconate titanate (PZT) piezoceramic sensors are able to detect anomalies of interest in 3.0-m long tubes. These sensors have hereby been extended to a 30.5-m long tube more representative of commercial power heat exchanger designs. The program will continue to investigate higher temperature piezoelectric ceramics and long-term performance of high temperature adhesives and sealants and mechanical pressure coupling for 350-500 °C lead reactor environments and higher temperature (700 °C) molten salt environments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Algae cultivation systems and methods adapted for weather variations

Algae cultivation systems and methods account for weather variations that can affect algae cultivation. In one system, an open raceway algae cultivation system includes a channel having a high section and a low liquid collection section. The channel is sloped to allow substantially all of an algae cultivation fluid in the high section to flow downwardly into the low liquid collection section. A barrier is removably positioned in the high section and a drain is positioned in the high section such that, when substantially all of the algae cultivation fluid has collected in the low liquid collection section, any rainwater that falls in the high section flows into the drain, without the rainwater mixing with the algae cultivation fluid in the low liquid collection section.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Robotic pendant drop: containerless liquid for μs-resolved, AI-executable XPCS

Abstract The dynamics and structure of mixed phases in a complex fluid can significantly impact its material properties, such as viscoelasticity. Small-angle X-ray Photon Correlation Spectroscopy (SA-XPCS) can probe the spontaneous spatial fluctuations of the mixed phases under various in situ environments over wide spatiotemporal ranges (10 −6 –10 3 s /10 −10 –10 −6 m). Tailored material design, however, requires searching through a massive number of sample compositions and experimental parameters, which is beyond the bandwidth of the current coherent X-ray beamline. Using 3.7-μs-resolved XPCS synchronized with the clock frequency at the Advanced Photon Source, we demonstrated the consistency between the Brownian dynamics of ~100 nm diameter colloidal silica nanoparticles measured from an enclosed pendant drop and a sealed capillary. The electronic pipette can also be mounted on a robotic arm to access different stock solutions and create complex fluids with highly-repeatable and precisely controlled composition profiles. This closed-loop, AI-executable protocol is applicable to light scattering techniques regardless of the light wavelength and optical coherence, and is a first step towards high-throughput, autonomous material discovery.

36 MATERIALS SCIENCE↗