Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Liquid flow”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Experimentally validated high-fidelity simulations of a liquid jet in supersonic crossflow

Here, utilizing recent advancements in computational schemes for compressible, multiphase flows, this work features a parametric study of a pure liquid jet in supersonic crossflow that involves simulating the atomization process for four values of momentum-flux ratio. These simulations are validated against experimental results measured with high-speed X-ray imaging, which confirm the accuracy of the numerical approach. Also, the effect of numerical resolution on some flow behavior is investigated, revealing convergence of the jet shape and surface instability wavelength. Analysis of the resulting sprays includes statistical descriptions of the liquid distribution, liquid structures created through breakup, interfacial instabilities, and dominant flow features. As the flowrate increases, the spray penetrates further, it becomes more disperse, and less liquid impacts the wall, but the droplet size distribution changes little. The wavelength of instabilities on the windward side of the jet diverges from measured trends in subsonic crossflows. In a visualization of the time-averaged flow, counter-rotating vortices are observed along the jet core and in the wake, affecting the process of primary atomization and early droplet trajectories.

42 ENGINEERING↗
DOE OSTI2020

Deep learning interfacial momentum closures in coarse-mesh CFD two-phase flow simulation using validation data

Multiphase flow phenomena have been widely observed in the industrial applications while it remains a challenging yet unsolved problems. Three-dimensional computational fluid dynamics (CFD) approaches resolve the flow fields on a finer special and temporal scales which can complement the dedicated experimental study. However, closures have to be introduced to reflect the underlying physics in multiphase flow. Among them, the interfacial forces, including drag, lift, turbulent dispersion and wall lubrication forces, play in important role on the bubble’s distribution and migration in liquid-vapor two-phase flow. Development of those closures traditionally rely on the experimental data and analytical derivation with simplified assumptions which usually cannot deliver a universal solution across wide range of flow conditions. In this paper, a data-driven approach, named as Feature Similarity Measurement (FSM), is developed and applied to improve the simulation capability of two-phase flow with coarse-mesh CFD approach. Interfacial momentum transfer in adiabatic bubbly flow serves as the focus of the present study. Both a mature and a simplified set of interfacial closures are taken as the low fidelity data. Experimental data and fine mesh CFD simulations results are adopted as high-fidelity data. Qualitative and quantitative analysis are performed in this paper which reveals that FSM can substantially improve the prediction of coarse mesh CFD model regardless of the choice of interfacial closures and it provides scalability and consistency across discontinuous flow regimes. Furthermore, it demonstrates that data-driven method can aid the multiphase flow modeling by exploring the connections between local physical features and simulation errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Simulation of Boiling Two-Phase Flow in a Helical Coil Steam Generator Using the Spectral Element Code Nek-2P

To enable the design of a light water small modular reactor, the boiling flow inside a helical coil steam generator has been simulated with the two-fluid model in Nek-2P. Nek-2P is the multiphase branch of the spectral element code Nek5000. Details of the implementation of the two-fluid model and the included closure models are discussed. The presented closure models include interactions for momentum, heat, and mass transfer between phases. Models for the drag, lift, and turbulent dispersion forces are included. The complete model is fully consistent in the limits of both phasic volume fractions approaching zero and is able to simulate flows of dispersed vapor, continuous liquid, dispersed liquid, continuous vapor, or any combination thereof. The closure models and their implementation in Nek-2P have been validated by comparing to experimental data for a boiling flow, demonstrating excellent agreement. Results from the simulation of the helical coil indicate strong phasic separation driven by the effects of buoyancy and inertia. Finally, significant differences were observed in the results compared to simulations performed using Star-CCM+, although these differences were somewhat expected.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

System and method for rapid, high throughput, high pressure synthesis of materials from a liquid precursor

The present disclosure relates to a system and method for synthesis of condensed nano-materials to at least one of create nanoparticles or modify existing nanoparticles. In one embodiment the system may have a source of liquid precursor, with the liquid precursor including a compound therein. A flow control element and a compression wave generating subsystem are also included. The flow control element is in communication with the source of the liquid precursor and creates a jet of liquid precursor. The compression wave generating subsystem drives a compression wave through at least a substantial portion of a thickness of the jet of liquid precursor to sufficiently compress the jet of liquid precursor, and to increase pressure and temperature of the jet of liquid precursor, to at least one of create nanoparticles or modify existing nanoparticles.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Predicted viscosity of liquid HMX up to 40 GPa

Viscous flow serves as a significant heating mechanism during the formation of hot spots, but the shear viscosity which determines this response is poorly characterized for most high explosives. Most recently, a model was proposed for the shear viscosity of liquid HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane) that was fit to pressures reaching 5 GPa, but this work uncovered uncertainties in the viscosity at 0 GPa and remains untested at higher pressures. We use molecular dynamics (MD) simulations and the Green-Kubo formalism to predict the temperature- and pressure-dependent shear viscosity of HMX over the pressure interval 0 GPa≤P≤40 GPa. Reassessment of the viscosity at 0 GPa rules out several potential explanations for discrepancies between earlier reports; we attribute these differences to details of MD trajectory integration impacting molecular flexibility. The shear viscosity of HMX exhibits an Arrhenius temperature dependence at each pressure considered, with exponential prefactor and activation energy terms that are also strong functions of pressure. An analytic form for the viscosity is developed based on an extension of the well-known Andrade equation that simultaneously captures the temperature and pressure dependencies in the MD data up to 40 GPa. Comparison against a recently developed model for the viscosity of liquid RDX (1,3,5-trinitro-1,3,5-triazinane) shows that both materials exhibit similar functional dependencies with the viscosity of HMX being higher by roughly an order of magnitude at a given temperature-pressure state.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Are Contact Angle Measurements Useful for Oxide-Coated Liquid Metals?

This work establishes that static contact angles for gallium-based liquid metals have no utility despite the continued and common use of such angles in the literature. In the presence of oxygen, these metals rapidly form a thin (~1–3 nm) surface oxide “skin” that adheres to many surfaces and mechanically impedes its flow. This property is problematic for contact angle measurements, which presume the ability of liquids to flow freely to adopt shapes that minimize the interfacial energy. We show here that advancing angles for a metal are always high (>140°)—even on substrates to which it adheres—because the solid native oxide must rupture in tension to advance the contact line. The advancing angle for the metal depends subtly on the substrate surface chemistry but does not vary strongly with hydrophobicity of the substrate. During receding measurements, the metal droplet initially sags as the liquid withdraws from the “sac” formed by the skin and thus the contact area with the substrate initially increases despite its volumetric recession. Furthermore, the oxide pins at the perimeter of the deflated “sac” on all the surfaces are tested, except for certain rough surfaces. With additional withdrawal of the liquid metal, the pinned angle gets smaller until eventually the oxide “sac” collapses. Thus, static contact angles can be manipulated mechanically from 0° to >140° due to hysteresis and are therefore uninformative. We also provide recommendations and best practices for wetting experiments, which may find use in applications that use these alloys such as soft electronics, composites, and microfluidics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Polymer Stabilized Paraboloid Liquid Crystal Microlenses with Integrated Pancharatnam–Berry Phase

Aspheric lenses reduce aberration and provide sharper images with improved spot size compared to spherical lenses. Here this paper demonstrates that applying shear flow can produce plano-concave liquid crystal (LC) lens arrays with paraboloid aspheric profiles. The focal length of individual lenses, with a 0.2 mm aperture, decreases from 0.67 to 0.45 mm as the chiral dopant increases from 0 to 6 wt%. The focal length is also sensitive to the polarization state of the incoming light. The lenses are stabilized by photopolymerizing with 6 wt% of reactive monomer added to the LC. A qualitative explanation for the flow-induced lens formation and the optical properties of the lenses is provided. The potential tunability of the lenses in various fields and their use as paraboloid reflectors are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Large Eddy Simulation of Low-Reynolds-Number Turbulent Flow of Low-Prandtl-Number Fluid in a Tight Lattice Bundle for Assessment of Reynolds-Averaged Navier-Stokes Turbulence Model

The MARVEL (Microreactor Applications Research Validation and Evaluation) microreactor utilizes natural circulation as core cooling mechanism and liquid metal as a primary coolant. Moreover, the reactor core has a pitch-to-diameter ratio of 1.056, which is considered a tight lattice configuration. Numerous studies have widely reported that Reynolds-Averaged Navier-Stokes (RANS) turbulence models inaccurately predict heat transfer in liquid metals and fail to capture flow pulsations that can occur within tight lattices, leading to further inaccuracies in simulation results. Therefore, evaluating the accuracy of RANS turbulence models in the thermal-hydraulic analysis of the MARVEL microreactor core is crucial for assessing reactor safety. In this study, a Large Eddy Simulation (LES) of the MARVEL microreactor core subchannel was conducted and compared with a RANS simulation to evaluate the accuracies and conservatism of the RANS model. The flow pulsation in a tight lattice predicted by LES enhanced the heat transfer, whereas the RANS model underpredicted it. Consequently, the RANS model predicted the peak cladding temperature higher than the LES model, but the discrepancy between the two approaches was not significant due to the good thermal characteristics of the liquid metal. It can be concluded that the steady-state RANS model is effective for the thermal analysis of liquid-metal-cooled MARVEL microreactor core and can provide conservative predictions from a safety analysis perspective.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2022

Supercritical, liquid, and gas CO 2 reactive transport and carbonate formation in portland cement mortar

In this paper, we investigate carbonate formation and reactive transport rate in variably saturated portland cement mortars when high concentrations of gas, liquid, or supercritical CO 2 flow through their pore network. Xray computed tomography completed during CO 2 flow is used to quantify the microstructural evolution as the mortar carbonates. After in situ tests, higher resolution scans, thermogravimetric analysis, and desorption isotherm analysis are performed to further quantify microstructural changes. We found that at dry conditions supercritical CO 2 moves more rapidly through the pore space and precipitates more carbonates than liquid or gas CO 2 . However, at 50% degree of saturation (DOS) the CO 2 state did not affect the rate of transport in that each specimen exposed to a different CO 2 state carbonated within the first hour of CO 2 exposure. When the pore space is at 50 or 100% DOS, supercritical CO 2 did not react with hydration products more rapidly nor did it result in more carbonate formation during exposure compared to gas or liquid CO 2 . The amount of Ca(OH) 2 that contributes to CaCO 3 formation is correlated to the DOS. For the mortar composition analyzed, Ca(OH) 2 contributes to approximately 40% of the carbonates formed in the 50% DOS specimens and 15% in the 100% DOS specimens. In other words, as the amount of moisture in the pore space increases, phases other than Ca(OH) 2 contribute to more than 50% of the total CaCO 3 formed.

42 ENGINEERING↗
DOE OSTI2021

The effect of transient needle lift on the internal flow and near-nozzle spray characteristics for modern GDI systems investigated by high-speed X-ray imaging

The development of the injector nozzle is a dynamic area in regard of several technical aspects. At first, the internal flow influences the near-field spray characteristics via various phenomena such as cavitation and turbulence. However, these phenomena are not fully understood due to their extremely fast, complex and multiscale nature. Furthermore, it governs the spray targeting inside the combustion chamber. High-speed X-ray imaging of GDI injector nozzles is performed in this study. The experimental results presented are related to the internal flow and primary breakup of discharged liquid jets. The injectors used are equipped with nozzles made of aluminum which have been specially developed for these investigations to enhance optical accessibility. The visualization of the needle motion, in-nozzle flow and the primary breakup region provides several exciting observations. First, the needle lift tracking exhibits short overshooting right before the steady-state of the injection phase. Additionally, this event leads to a short-term, however, significant change in the associated performance of the breakup. This phenomenon is found to be a consequence of the transient behavior of the in-nozzle flow. It is shown that under some circumstances hydraulic flip may occur during this overshooting period. The primary jet breakup region is visualized and evaluated by means of image processing. Thus, the transient behavior of liquid jet expansion is quantified in the vicinity of the nozzle. It is observed that the liquid jet direction deviates from the hole axis already at the nozzle outlet, which is caused by internal flow characteristics.

42 ENGINEERING↗
DOE OSTI2024

Phase field-volumetric lattice Boltzmann model of ion uptake in porous nuclear waste form materials under continuous flow

The flow field within the mesopores of sorbent particles plays a crucial role in radionuclide diffusion and ion uptake kinetics, thus, impacting the overall performance of porous nuclear waste form materials. To fundamentally understand the influence of microstructures and material properties on the radionuclide absorption and retention processes requires a coupled multi-physics model that considers the advection and diffusion within the flow field, the reaction at liquid-solid interfaces, and finally, the solid-state diffusion within a complex nanoporous medium. Here, this study employs the volumetric lattice Boltzmann method (VLBM) to accurately and efficiently calculate the steady state velocity field inside the mesopores of sorbent particles. The obtained velocity field is then utilized to calculate the advection of ions in the steady flow. A phase field (PF) model of ion uptake is used to describe the reaction occurring at the solid-liquid interface and diffusion inside the porous medium. The integrated PF-VLBM model is verified in terms of the mass conservation and numerical efficiency and validated qualitatively with experimental observation data. Then, it is applied to study the influence of thermodynamic and kinetic properties, as well as flow field conditions on the ion uptake kinetics. The numerical results demonstrate that the ion uptake kinetics in porous particles has three distinct stages, which is in agreement with the observations in continuous flow experiments. In the first stage, the kinetics is predominantly controlled by the flow field and ion diffusivity in the liquid phase. The kinetics in the second stage is primarily governed by ion diffusivity in the solid phase. In the third stage the system reaches a dynamic equilibrium with a net zero uptake flux at the interface. It is also found that porous structures significantly affect the efficiency and capacity of ion uptake. The simulation results can help to understand the physics behind the observed ion uptake kinetics in experiments and to facilitate the development of constitutive equations that can account for heterogeneous microstructures in engineering performance codes.

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

Methods and systems for redox flow battery electrolyte hydration

Methods and systems are provided for transporting and hydrating a redox flow battery system with a portable field hydration system. In one example, the redox flow battery system may be hydrated with the portable field hydration system in a dry state, in the absence of liquids. In this way, a redox flow battery system may be assembled and transported from a battery manufacturing facility to an end-use location off-site while the redox flow battery system is in the dry state, thereby reducing shipping costs, design complexities, as well as logistical and environmental concerns.

Song, Yang↗
DOE OSTI2021

Modeling impurity concentrations in liquid argon detectors

Impurities in noble liquid detectors used for neutrino and dark matter experiments can significantly impact the quality of data. We present an experimentally verified model for describing the dynamics of impurity distributions in liquid argon (LAr) detectors. The model considers sources, sinks, and transport of impurities within and between the gas and liquid argon phases. Measurements of oxygen concentrations in a 20-L LAr multi-purpose test stand are compared to calculations made with this model to show that an accurate description of the concentrations under various operational conditions can be obtained. A result of this analysis is a determination of Henry’s coefficient for oxygen in LAr. Additionally, these calculations also show that some processes have small effects on the impurity dynamics and excluding them yields a solution as a sum of two exponential terms. This solution provides a simple way to extract Henry’s coefficient with negligible approximation error. It is applied to the data and the Henry’s coefficient for oxygen in LAr is obtained as 0.84 +0.09 -0.05 , consistent with literature results. Finally, based on the analysis of the data with the model, we further suggest that, for a large liquid argon detector, barriers to flow (“baffles”) installed in the gas phase to restrict flow can help reduce the ultimate impurity concentration in the LAr.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2022

On the effect of mixing-driven vaporization in a homogeneous relaxation modeling framework

The homogeneous relaxation model (HRM) is one of the most widely used models to describe the liquid–gas phase transition in multiphase flows due to the occurrence of cavitation. However, in its original formulation, the HRM does not account for the presence of ambient gas species, which generally limits its applicability to the injector's internal flow where ambient gases are negligible. In this work, a mixing-driven vaporization (MDV) model was developed to extend the capability of the HRM in handling the mixing effect in the regions external to the nozzle, where vapor–liquid equilibrium for multi-species mixtures of fuel and ambient gas is considered. Herein, to assess the model performance, simulations of the Engine Combustion Network's Spray G injector were performed with the HRM and the MDV model under both flash-boiling and evaporating conditions. It was found that the MDV model led to a better match against x-ray measurements of fuel density in the near-nozzle region. In contrast to the HRM, the MDV model was able to reproduce the vaporization process in the mixing zone at the edge of the fuel jet, which aligns with the expected physics. This resulted in substantial differences in the prediction of other flow characteristics such as mixture temperature and pressure. Furthermore, this work demonstrates that evaporation timescales have a considerable effect on the MDV model's predictions, as shown by a parametric study in which a time factor was introduced to mimic the effect of different timescales due to different phase change mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Alkaline manganese redox flow battery with inhibitor

A redox flow battery includes a redox flow cell and a supply and storage system external of the redox flow cell. The supply and storage system includes first and second electrolytes for circulation through the redox flow cell. The first electrolyte is a liquid electrolyte having electrochemically active manganese species with multiple, reversible oxidation states in the redox flow cell. The electrochemically active manganese species may undergo reactions that cause precipitation of manganese oxide solids. The first electrolyte includes an inhibitor that limits the self-discharge reactions. The inhibitor includes an oxoanion compound.

Davenport, Timothy↗
DOE OSTI2020

Development of Methods for In-line Monitoring of Yield Stress During the Transfer of Radioactive Waste - 20435

Radioactive waste slurries at the Savannah River Site (SRS) and Hanford often behave as non-Newtonian fluids with a yield stress. The yield stress increases the size of equipment needed to transport or mix these slurries and can increase the risk of process upsets. This property is often measured by collection of samples and transportation to a laboratory, which can increase the risk of personnel exposure and even change the rheology. In addition, slurry rheology has been observed to change with time. Real-time in-line monitoring (RTIM) of rheological properties of slurries can have substantial benefits to critical processes which are currently in operation or will be part of future plans in Department of Energy (DOE) complex processes. The need for an appropriate RTIM technology that can potentially minimize or eliminate numerous sampling operations was recently identified. This research at Florida International University focuses on addressing the need for a reliable RTIM technology at US Department of Energy sites. The objective of this work is to investigate and develop appropriate and reliable technologies that can perform accurate and easy measurements of yield stress with minimal changes to the flow conditions at US Department of Energy sites. Use of ultrasound in longitudinal and shear modes was considered as a method for monitoring of yield stress variations with no disturbance to the flow. Measurement of shear wave velocity (in shear mode) and monitoring of sound velocity variations with start of flow from rest (longitudinal mode) are the underlying principles. Measurement of pressure loss and liquid rise (due to static differential pressure) during flow of slurries were two other methods that were considered and investigated. In addition, analysis of stability depth of a penetrometer in slurry medium (penetrometer method) revealed a potential suitability. This paper presents efforts associated with investigation of different methods mentioned earlier for kaolin-water mixtures at different concentrations. Theoretical analysis showed that all methods under investigation were capable of monitoring yield stress in a desirable range of 0 to 70 pascal. A bench scale test setup was constructed for investigation of pressure loss and liquid rise methods. Simulants at different concentrations were created for initial testing campaigns. For ultrasonic testing, effects of yield stress on sound propagation speed in longitudinal and shear modes were investigated for simulants placed inside PVC pipes. Parameters such as simulant concentration and pipe diameter were varied to find optimal values. It was found that spacing (or medium thickness) of 1/8 inch or less was necessary to observe echoes with sufficient strength. In addition, a significant effect of yield stress on sound propagation velocity was observed. For the liquid rise method, kaolin-water simulant at 1.17 specific gravity was pumped through a 1/4 inch pipe and liquid height was measured in a vertical riser branching off a three-way connection. Static gauge pressure was varied within [0.3 to 1.1] psi range by increase of pumping pressure and flow. Results showed consistent yield stress data within 10% deviation from values reported by a rheometer. (authors)

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

Development of Explainable Data-Driven Turbulence Models with Application to Liquid Fuel Nuclear Reactors

Liquid fuel nuclear reactors offer innovative possibilities in terms of nuclear reactor designs and passive safety systems. Molten Salts Reactors (MSRs) with a fast spectrum are a particular type of these reactors using liquid fuel. MSFRs often involve large open cavities in their core in which the liquid fuel circulates at a high speed to transport the heat generated by the nuclear reactions into the heat exchangers. This high-speed flow yields a turbulent field with large Reynolds numbers in the reactor core. Since the nuclear power, the neutron precursor’s transport and the thermal exchanges are strongly coupled in the MSFR’s core cavity, having accurate turbulent models for the liquid fuel flow is necessary to avoid introducing significant errors in the numerical simulations of these reactors. Nonetheless, high-accuracy simulations of the turbulent flow field in the reactor cavity of these reactors are usually prohibitively expensive in terms of computational resources, especially when performing multiphysics numerical calculations. Therefore, in this work, we propose a novel method using a modified genetic algorithm to optimize the calculation of the Reynolds Shear Stress Tensor (RST) used for turbulence modeling. The proposed optimization methodology is particularly suitable for advanced liquid fuel reactors such as the MSFRs since it allows the development of high-accuracy but still low-computational-cost turbulence models for the liquid fuel. We demonstrate the applicability of this approach by developing high accuracy Reynolds-Averaged Navier–Stokes (RANS) models (averaged flow error less than 5%) for a low and a large aspect ratio in a Backward-Facing Step (BFS) section particularly challenging for RANS models. The newly developed turbulence models better capture the flow field after the boundary layer tipping, over the extent of the recirculation bubble, and near the boundary layer reattachment region in both BFS configurations. The main reason for these improvements is that the developed models better capture the flow field turbulent anisotropy in the bulk region of the BFS. Then, we illustrate the interest in using this turbulence modeling approach for the case of an MSFR by quantifying the impact of the turbulence modeling on the reactor key parameters.

42 ENGINEERING↗