Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interfacial flows”

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 OSTI2022

Non-equilibrium molecular dynamics simulation to evaluate the effect of confinement on fluid flow in silica nanopores

The flow of fluids in nano-confinement has applications in separations, water purification, medical systems and in the recovery of fluids in petroleum systems. It is believed that fluids do not obey continuum laws when flowing in nanopores (or in confinement). This paper attempts to shed light on various nano confinement effects such as fluid-wall interactions, pore size and molecular geometry, using molecular dynamic simulation. Here, water, hexane, and methanol flow behaviors were simulated for pore diameters ranging from 1 to 8 nm. In addition to the density analysis of the confined fluids, the fluid flows in saturated nanopores were simulated by the sectional flow method. Water and methanol molecules were completely stabilized in the 1 nm pore. Hexane molecules could not enter the 1 nm pore, due to geometric considerations. For the 2 – 8 nm pores, all fluids showed reduced flow rate compared to the Hagen–Poiseuille flow, due to an interfacial molecular layer stabilized at the pore surface. Flow reduction observed in these studies is contrary to significant flow enhancements observed for fluid flow in carbon nanotubes of similar dimensions. Water flow showed almost constant stick length for all the pore sizes. Methanol flow had the largest stick length. Hexane flow was reduced because of overcrowding of molecules at the pore surface. The effect of confinement diminished with an increase in pore diameter.

42 ENGINEERING↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2016

Newton-Krylov Method in Applications of Solving Two-phase Problems using Drift Flux Model

In this paper, Jacobian-free Newton-Krylov (JFNK) method is investigated in an application to implicitly solve two-phase flow problems using four-equation drift flux model. The closure models include the EPRI drift flux correlations and additional constitutive models to determine flow regimes, wall boiling and interfacial heat/mass transfer, two-phase flow wall friction, etc. Different from many traditional computer codes, fully implicit methods are used for the time integration. The resulted nonlinear discretized equation system is solved using the JFNK method. Expensive and error prone derivation and implementation of the analytical Jacobian matrix are avoided. Numerical results are successfully validated using existing experimental data on flow boiling under forced convection conditions in both pipe and rod bundle geometries.

42 ENGINEERING↗
NASA NTRS2021

Experimental Results and Interfacial Lift-off Model Predictions of Critical Heat Flux for Flow Boiling with Subcooled Inlet Conditions – In Preparation for Experiments Onboard the International Space Station

This study investigates critical heat flux (CHF) for subcooled flow boiling of n-Perfluorohexane based on results of pre-launch Earth-gravity Mission Sequence Tests (MSTs) of the Flow Boiling and Condensation Experiment (FBCE), which was launched to the International Space Station (ISS) in August 2021. CHF measurements were made in a rectangular channel having a 2.5 mm by 5 mm cross-section and a 114.6-mm long heated segment. Both single-sided and double-sided heating were tested in vertical upflow in Earth gravity for a variety of inlet conditions. The inlet subcooling was varied in the range of 0.4 – 32.0°C and encompassed both near-saturated and highly subcooled conditions. Experimental trends and high-speed video records were investigated to better understand the mechanism of CHF. Overall trends show CHF increases as flow rate and/or inlet subcooling are increased. Flow features from the events around CHF justify the applicability of the Interfacial Lift-off Model and the determination of limiting criteria for its application. The present experimental data are combined with prior databases for various flow orientations with respect to Earth gravity and microgravity data collected on parabolic flights. Predictions are made using the Interfacial Lift-off Model for this consolidated subcooled-inlet FBCE-CHF database. A heat utility ratio was included in the model to capture the effects of subcooling and corresponding thermodynamic non-equilibrium. An overall mean absolute error of 19.04% indicates good predictive capability of the model for both heating configurations, different gravity environments, and a wide range of inlet subcooling.

flow boiling↗
DOE OSTI2021

Transport Modeling of Kinetically Limited Microscale Extraction Systems: Droplet and Supported Liquid Membrane Separations

Microfluidic droplet extraction devices can be used to efficiently measure solute phase transfer kinetics in liquid–liquid extraction systems. In this work, a droplet extraction apparatus was constructed using off-the-shelf and custom 3D printed parts for less than 110 USD. The interfacial kinetics of the transfer of Pu(IV) between an aqueous nitric acid phase and a 5 v/v% tridodecylamine (TDA) in toluene phase were characterized by determining the extent of Pu(IV) transfer between phases at different flow rates. Forward interfacial rate constants calculated from this data using an analytical model are 5.6 × 10 –3 mm s –1 under extraction (aqueous-to-organic solute transfer) conditions and 3.3 × 10 –3 mm s –1 under stripping (organic-to-aqueous solute transfer) conditions. The analytical model is valid when the rate of solute transfer across the interface between the two liquid phases has pseudo-first-order reversible reaction kinetics and is slow relative to the rate of solute transport to the interface. The change in the outlet concentration with residence time determined using the analytical model was compared with the results of a numerical model accounting for advection/diffusion to the interface in addition to kinetically limited interfacial solute transfer. These calculations revealed that for systems with forward interfacial rate constants greater than approximately 10 –2 mm s –1 , the analytical model must be supplemented with a numerical transport model accounting for an additional advection/diffusion mass transfer resistance. This threshold corresponds to a Damköhler number of one. The measured interfacial kinetic rate constants for Pu(IV)/TDA were used in a supported liquid membrane (SLM) numerical model, which showed that the Pu(IV) phase transfer kinetics are fast enough for use in a practical separation.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Asymmetric fluctuations and self-folding of active interfaces

We study the structure and dynamics of the interface separating a passive fluid from a microtubule-based active fluid. Turbulent-like active flows power giant interfacial fluctuations, which exhibit pronounced asymmetry between regions of positive and negative curvature. Experiments, numerical simulations, and theoretical arguments reveal how the interface breaks up the spatial symmetry of the fundamental bend instability to generate local vortical flows that lead to asymmetric interface fluctuations. The magnitude of interface deformations increases with activity: In the high activity limit, the interface self-folds invaginating passive droplets and generating a foam-like phase, where active fluid is perforated with passive droplets. These results demonstrate how active stresses control the structure, dynamics, and break-up of soft, deformable, and reconfigurable liquid–liquid interfaces.

active fluid↗
DOE OSTI2016

Modeling Electric Double Layer Effects on Charge Transfer at Flow Battery Electrode/Electrolyte Interfaces

The proposal aims to model interfacial processes associated with redox flow batteries (RFB) for grid/stationary storage, with the ultimate goal of helping design new electrolytes, electrodes, redox species, and interfaces. It addresses the linked problems of electric double layer (EDL) structure and electron transfer at model electrode/electrolyte interfaces. Fast electron transfer partly determines which redox species are viable for flow batteries (along with solubility, viscosity, etc.). Commercial systems apply aqueous, vanadium-based complexes, but other choices and organic electrolytes with larger voltage windows are the subject of active research at Sandia and elsewhere. The high salt concentration present in flow battery electrolytes yields non-trivial EDL consisting of solvents, counter-ions, co-ions, and redox species (at different charge states at the charging voltage threshold). EDL strongly influences electron transfer, especially when redox “mediators” are used. Modeling this key missing information is the main challenge. Our proposal seeks to apply Sandia’s LAMMPS molecular dynamics code. As a proof of principle, we study minimal model systems (graphite electrode, with ferrocene and fluornone as redox-active species). These model "catholyte" and "anolyte" molecules exhibit low reorganization energies (in the sense of Marcus theory) and are most amenable to EDL simulations which involve switching of redox states to mimic electron transfer. This study will enable future modeling and design of redox species relevant to flow batteries and redox mediators used at Sandia for other purposes.

25 ENERGY STORAGE↗
NASA NTRS1991

Electrical impedance imaging in two-phase, gas-liquid flows: 1. Initial investigation

The determination of interfacial area density in two-phase, gas-liquid flows is one of the major elements impeding significant development of predictive tools based on the two-fluid model. Currently, these models require coupling of liquid and vapor at interfaces using constitutive equations which do not exist in any but the most rudimentary form. Work described herein represents the first step towards the development of Electrical Impedance Computed Tomography (EICT) for nonintrusive determination of interfacial structure and evolution in such flows.

Lin, J. T.↗
NASA NTRS2016

Flow Boiling and Condensation Experiment (FBCE) for the International Space Station

An effective means to reducing the size and weight of future space vehicles is to replace present mostly single-phase thermal management systems with two-phase counterparts. By capitalizing upon both latent and sensible heat of the coolant rather than sensible heat alone, two-phase thermal management systems can yield orders of magnitude enhancement in flow boiling and condensation heat transfer coefficients. Because the understanding of the influence of microgravity on two-phase flow and heat transfer is quite limited, there is an urgent need for a new experimental microgravity facility to enable investigators to perform long-duration flow boiling and condensation experiments in pursuit of reliable databases, correlations and models. This presentation will discuss recent progress in the development of the Flow Boiling and Condensation Experiment (FBCE) for the International Space Station (ISS) in collaboration between Purdue University and NASA Glenn Research Center. Emphasis will be placed on the design of the flow boiling module and on new flow boiling data that were measured in parabolic flight, along with extensive flow visualization of interfacial features at heat fluxes up to critical heat flux (CHF). Also discussed a theoretical model that will be shown to predict CHF with high accuracy.

flow boiling↗
NASA NTRS1989

Effect of particle velocity fluctuations on the inertia coupling in two-phase flow

Consistent forms for the interfacial force, the interfacial pressure, the Reynolds stresses and the particle stress have been derived for the inviscid, irrotational incompressible flow of fluid in a dilute suspension of spheres. The particles are assumed to have a velocity distribution, giving rise to an effective pressure and stress in the particle phase. The velocity fluctuations also contribute in the fluid Reynolds stress and in the (elastic) stress field inside the spheres. The relation of these constitutive equations to the force on an individual sphere is discussed.

Drew, Donald A.↗
DOE OSTI2018

Experimental Investigation of the System Nonmetallic Inclusion-Molten Steel-Refractory System at High Temperatures

In continuous casting, the buildup of particles in the flow control system is called clogging. Clogging may be caused by preexisting nonmetallic inclusions attached to the refractory wall, interrupting or disturbing the fluid flow conditions. The interfacial properties of the nonmetallic inclusion-molten steel-refractory system are believed to play a key role in the clogging tendency of selected steel grades in contact with ceramics. The combination of a high-temperature laser scanning confocal microscope (HT-LSCM), which is employed for in situ observation of inclusions’ behavior in the liquid steel or at the steel/ceramic interface, and a high-temperature drop shape analyzer (HT-DSA), which is utilized for contact angle measurements, enabled the investigation of reactions and interactions in the described system at experimental temperatures of approximately 1600 °C. The objective of this study was to develop a method to obtain qualitative information on the deposition of nonmetallic inclusions (NMIs) in the submerged entry nozzle (SEN) during casting. The behavior of solid NMIs from a calcium-treated steel at the interface between liquid steel and Al{sub 2}O{sub 3} or MgO was investigated by means of HT-LSCM. During the test evaluation, inclusions were found to separate at the steel/ceramic interface due to an attractive force. This force increased with the NMIs’ size as well as with the ceramic particle change. The HT-DSA provided contact angle data between molten steel and the ceramic used for implementation into the theoretical capillary force model. A clear correspondence between the estimated and calculated capillary forces was observed. The change in the wettability of the system, as induced by the change in the ceramic from MgO to Al{sub 2}O{sub 3}, showed a reduction of the attractive forces, indicating a weaker separation toward the SEN for the same inclusion type. This study confirmed the strong potential of wettability investigations to predict the clogging tendency of certain steel grades.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Petascale Application of a Coupled CPU-GPU Algorithm for Simulation and Analysis of Multiphase Flow Solutions in Porous Medium Systems, In: 2014 IEEE 28th International Parallel and Distributed Processing Symposium

Large-scale simulation can provide a wide range of information needed to develop and validate theoretical models for multiphase flow in porous medium systems. In this paper, we consider a coupled solution in which a multiphase flow simulator is coupled to an analysis approach used to extract the interfacial geometries as the flow evolves. This has been implemented using MPI to target heterogeneous nodes equipped with GPUs. The GPUs evolve the multiphase flow solution using the lattice Boltzmann method while the CPUs compute up scaled measures of the morphology and topology of the phase distributions and their rate of evolution. Our approach is demonstrated to scale to 4,096 GPUs and 65,536 CPU cores to achieve a maximum performance of 244,754 million-lattice-node updates per second (MLUPS) in double precision execution on Titan. In turn, this approach increases the size of systems that can be considered by an order of magnitude compared with previous work and enables detailed in situ tracking of averaged flow quantities at temporal resolutions that were previously impossible. Furthermore, it virtually eliminates the need for post-processing and intensive I/O and mitigates the potential loss of data associated with node failures.

Computer Science↗
NASA NTRS1999

Surface Tension and Fingering of Miscible Interfaces

Experiments on miscible, buoyantly unstable reaction-diffusion fronts and non-reacting displacement fronts in Hele-Shaw cells show a fingering-type instability whose wavelengths (lambda*) are consistent with an interfacial tension (sigma) at the front caused by the change in chemical composition, even though the solutions are miscible in all proportions. In conjunction with the Saffman-Taylor model, the relation sigma = K/tau, where tau is the interface thickness and K approximately equal 4 +/- 2 x 10(exp -6) dyne, enables prediction of our measured values of lambda* as well as results from prior experiments on miscible interfaces. These results indicate that even for miscible fluids, surface tension is generally a more significant factor than diffusion in interfacial stability and flow characteristics.

Abib, Mohammed↗
DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part I

In this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In conclusion, these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II

Here, in this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

Multiphase flow and phase change in microgravity: Fundamental research and strategic research for exploration of space

NASA is preparing to undertake science-driven exploration missions. The NASA Exploration Team's vision is a cascade of stepping stones. The stepping-stone will build the technical capabilities needed for each step with multi-use technologies and capabilities. An Agency-wide technology investment and development program is necessary to implement the vision. The NASA Exploration Team has identified a number of areas where significant advances are needed to overcome all engineering and medical barriers to the expansion of human space exploration beyond low-Earth orbit. Closed-loop life support systems and advanced propulsion and power technologies are among the areas requiring significant advances from the current state-of-the-art. Studies conducted by the National Academy of Science's National Research Council and Workshops organized by NASA have shown that multiphase flow and phase change play a crucial role in many of these advanced technology concepts. Lack of understanding of multiphase flow, phase change, and interfacial phenomena in the microgravity environment has been a major hurdle. An understanding of multiphase flow and phase change in microgravity is, therefore, critical to advancing many technologies needed. Recognizing this, the Office of Biological and Physical Research (OBPR) has initiated a strategic research thrust to augment the ongoing fundamental research in fluid physics and transport phenomena discipline with research especially aimed at understanding key multiphase flow related issues in propulsion, power, thermal control, and closed-loop advanced life support systems. A plan for integrated theoretical and experimental research that has the highest probability of providing data, predictive tools, and models needed by the systems developers to incorporate highly promising multiphase-based technologies is currently in preparation. This plan is being developed with inputs from scientific community, NASA mission planners and industry personnel. The fundamental research in multiphase flow and phase change in microgravity is aimed at developing better mechanistic understanding of pool boiling and ascertaining the effects of gravity on heat transfer and the critical heat flux. Space flight experiments conducted in space have shown that nucleate pool boiling can be sustained under certain conditions in the microgravity environment. New space flight experiments are being developed to provide more quantitative information on pool boiling in microgravity. Ground-based investigations are also being conducted to develop mechanistic models for flow and pool boiling. An overview of the research plan and roadmap for the strategic research in multiphase flow and phase change as well as research findings from the ongoing program will be presented.

Singh, Bhim S.↗
DOE OSTI2016

An Alternative Frictional Boundary Condition for Computational Fluid Dynamics Simulation of Friction Stir Welding

For better application of numerical simulation in optimization and design of friction stir welding (FSW), this paper presents a new frictional boundary condition at the tool/workpiece interface for computational fluid dynamics (CFD) modeling of FSW. The proposed boundary condition is based on an implementation of the Coulomb friction model. Using the new boundary condition, the CFD simulation yields non-uniform distribution of contact state over the tool/workpiece interface, as validated by the experimental weld macrostructure. It is found that interfacial sticking state is present over large area at the tool-workpiece interface, while significant interfacial sliding occurs at the shoulder periphery, the lower part of pin side, and the periphery of pin bottom. Due to the interfacial sticking, a rotating flow zone is found under the shoulder, in which fast circular motion occurs. The diameter of the rotating flow zone is smaller than the shoulder diameter, which is attributed to the presence of the interfacial sliding at the shoulder periphery. For the simulated welding condition, the heat generation due to friction and plastic deformation makes up 54.4 and 45.6% of the total heat generation rate, respectively. In conclusion, the simulated temperature field is validated by the good agreement to the experimental measurements.

36 MATERIALS SCIENCE↗