SEARCH · Engineering Papers

Results for “Fluid interfaces”

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

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Recent advances in computational mathematics and applications

We are honored to bring you this special issue dedicated to recent advances in computational mathematics and applications in science and engineering. The papers in this special issue were presented in ”Conference on Computational Mathematics and Applications (CCMA)” held at the University of Nevada Las Vegas (UNLV) during October 25–27, 2019. The conference has attracted about 90 attendees from six countries. The 15 papers in the special issue comprise a diverse collection of theoretical numerical analysis as well as applications in various subjects. Topics cover optimal control of PDEs, discontinuous Galerkin methods, fluid flow in porous media, turbulence flow, static and moving interface problems, complex fluids composed by the mixture of Newtonian fluid and nematic (liquid crystal) flows, modeling of ocean–atmosphere system, image segmentation algorithm, uncertainty quantification problems for turbulence model and Maxwell’s equations, iterative solver for systems resulting from mixed methods, and novel methods for solving systems of nonlinear equations. As the Guest Editors we would like to thank Dr. Leland Jameson at NSF (National Science Foundation) for kindly supporting our conference proposal which made this conference possible. We would also like to thank Darren Sugrue from Elsevier who provided partial support for our conference, and Thennarasu Gunasekaran and his production team from Elsevier who produced this nice special issue. We are grateful to many people (Dr. Zhijian Wu, Lori Ornelas, Elsa Juarez, and many Ph.D. students) at the Department of Mathematical Sciences of UNLV who provided tremendous support for the conference. Finally, we very much appreciate all authors who contributed their precious results to this special issue.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Influence of Cross Perturbations on Turbulent Kelvin–Helmholtz Instability

Kelvin–Helmholtz instability has been studied extensively in 2D. This study attempts to address the influence of turbulent flow and cross perturbation on the growth rate of the instability and the development of mixing layers in 3D by means of direct numerical simulation. Two perfect gases are considered to be working fluids moving as opposite streams, inducing shear instability at the interface between the fluids and resulting in Kelvin–Helmholtz instability. The results show that cross perturbation affects the instability by increasing the amplitude growth while adding turbulence has almost no effect on the amplitude growth. Furthermore, by increasing the turbulence intensity, a more distinct presence of the inner flow can be seen in the mixing layer of the two phases, and the presence of turbulence expands the range of high-frequency motion significantly due to turbulence structures. The results give a basis for which 3D Kelvin–Helmholtz phenomena should be further investigated using numerical simulation for predictive modeling, beyond the use of simplified 2D theoretical models.

Sementilli, Mae (ORCID:0000000207125928)↗

NASA NTRS1989

A free surface vortex method with weak viscous effects

The advantage of vortex methods for inviscid flows of several fluids separated by interfaces is that problems of this kind may be posed and solved numerically entirely on the interface surface. For free surface problems, where one of the fluids has negligible inertia, a small amount of viscosity produces a thin, weak rotational layer adjacent to the potential flow vortex sheet. Here, a method is described for including this weak viscous effect as a modification of the boundary conditions for the potential part of the flow, thus allowing this slightly viscous flow to be formulated as a vortex method on the free surface. This is important for unsupported free drops where the viscosity at the free surface is the only source of damping for oscillations, and for gas-filled bubbles, where it provides both drag and damping. Small applied shear stresses from aerodynamic, electrostatic, or other sources may be included provided they are of the same order as the potential flow shear stress or smaller.

Lundgren, T. S.↗

NASA NTRS2011

CFD Study of the Hydrocarbon Boost Low-Pressure Inducer and Kicker in the Presence of a Circumferential Groove

Results are presented of a computational fluid dynamics (CFD) study done in support of Marshall Space Flight Center's (MSFC) sub-scale water flow experiments of the Hydrocarbon Boost (HCB) Oxidizer Turbopump (OTP) being developed by the Air Force Research Laboratory (AFRL) and Aerojet. A circumferential groove may be added to the pump to reduce synchronous cavitation and subsequent bearing loads at a minimal performance cost. However, the energy may reappear as high order cavitation (HOC) that spans a relatively large frequency range. Thus, HOC may have implications for the full-scale OTP inducer in terms of reduced structural margin at higher mode frequencies. Simulations using the LOCI/Stream CFD program were conducted in order to explore the fluid dynamical impact of the groove on the low-pressure inducer and kicker. It was found that the circumferential groove has minimal head performance impact, but causes back-flowing high-swirl fluid to interact with the nearly-axial incoming fluid just above the inducer blades. The high-shear interface between the fluids is Kelvin-Helmholtz unstable, resulting in trains of low pressure regions or 'pearls' forming near the upstream edge of the groove. When the static pressure in these regions becomes low enough and they get cut by the blade leading edge, HOC is thought to occur. Although further work is required, the numerical models indicate that HOC will occur in the runbox of the AFRL/Aerojet HCB OTP. Comparisons to the ongoing water flow experiments will be discussed, as well as possible designs that may mitigate HOC while continuing to reduce synchronous cavitation. December 2011 MSS/LPS/SPS Joint Subcommittee Meeting ABSTRACT SUBMITTAL FORM

Coker, Robert↗

NASA NTRS2011

CFD Study of the J2-X LOX Two-Bladed Inducer in the Presence of a Circumferential Groove

Results are presented of a computational fluid dynamics (CFD) study done in support of water flow experiments of the J-2X Oxidizer Turbopump (OTP) 2-bladed inducer with a circumferential groove that were conducted at Marshall Space Flight Center (MSFC). Sub-scale water flow testing results indicate that the circumferential groove greatly reduces synchronous cavitation and subsequent bearing loads at a minimal performance cost, but the energy reappears as high order cavitation (HOC) that spans a relatively large frequency range. Thus, HOC may have implications for the full-scale OTP inducer in terms of reduced structural margin at higher mode frequencies. Simulations using the LOCI-Stream CFD program were conducted in order to explore the root physical cause of the HOC. It was found that the axial recirculation pattern in the circumferential groove causes high-swirl fluid to interact with the nearly-axial incoming fluid just above the inducer blades. The high-shear interface between the fluids is Kelvin-Helmholtz unstable, resulting in trains of low pressure regions or 'pearls' forming near the upstream edge of the groove. When the pressure in these regions becomes low enough and they get cut by the blade leading edge, HOC is thought to occur. Although further work is required, the numerical models indicate that the root cause of HOC is hydrodynamic. That is, the pearls are always present, even when cavitation is not. Comparisons to ongoing water flow experiments will be discussed, as well as predictions for the full-scale OTP inducer.

Coker, Robert↗

DOE OSTI2024

Heat exchanger assemblies having embedded sensors

Heat exchanger assemblies are provided that can include: a heat exchanger housing; at least one primary conduit operably coupled to the heat exchanger housing and configured to convey a primary heat exchange fluid; at least one secondary conduit operably coupled to the heat exchanger housing and configured to convey a secondary heat exchange fluid; at least one thermal interface between the primary and secondary fluids; and at least one sensor operably engaged with the thermal interface. Heat exchanger assemblies including molten salt, liquid metal, or water/steam as part of the heat exchange fluids of the heat exchanger assembly are provided. The heat exchanger assemblies can include: at least one thermal interface between primary and secondary heat exchange fluids of the heat exchanger assembly; and a sensor operably engaged with the at least one interface.

Glass, III, Samuel W.↗

DOE OSTI2020

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Rayleigh–Taylor and Richtmyer–Meshkov instabilities: A journey through scales

Hydrodynamic instabilities such as Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities usually appear in conjunction with the Kelvin–Helmholtz (KH) instability and are found in many natural phenomena and engineering applications. They frequently result in turbulent mixing, which has a major impact on the overall flow development and other effective material properties. This can either be a desired outcome, an unwelcome side effect, or just an unavoidable consequence, but must in all cases be characterized in any model. The RT instability occurs at an interface between different fluids, when the light fluid is accelerated into the heavy. The RM instability may be considered a special case of the RT instability, when the acceleration provided is impulsive in nature such as that resulting from a shock wave. Here in this pedagogical review, we provide an extensive survey of the applications and examples where such instabilities play a central role. First, fundamental aspects of the instabilities are reviewed including the underlying flow physics at different stages of development, followed by an overview of analytical models describing the linear, nonlinear and fully turbulent stages. RT and RM instabilities pose special challenges to numerical modeling, due to the requirement that the sharp interface separating the fluids be captured with fidelity. These challenges are discussed at length here, followed by a summary of the significant progress in recent years in addressing them. Examples of the pivotal roles played by the instabilities in applications are given in the context of solar prominences, ionospheric flows in space, supernovae, inertial fusion and pulsed-power experiments, pulsed detonation engines and Scramjets. Progress in our understanding of special cases of RT/RM instabilities is reviewed, including the effects of material strength, chemical reactions, magnetic fields, as well as the roles the instabilities play in ejecta formation and transport, and explosively expanding flows. The article is addressed to a broad audience, but with particular attention to graduate students and researchers who are interested in the state-of-the-art in our understanding of the instabilities and the unique issues they present in the applications in which they are prominent.

42 ENGINEERING↗

NASA NTRSDate not supplied

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the Volume of Fluid (VOF) and Sharp Interface multiphase 2D axisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. Turbulence modeling; turbulence damping at the interface; and using constant vs. temperature dependent fluid properties were shown to have the most profound influence on predicted tank pressures and temperatures. Current study indicates that including phase change at the interface into the computational model is crucial for accurate prediction of the tank self-pressurization process. The effect of accommodation coefficient was also studied. Tank pressure values predicted by the VOF, and Sharp Interface models are within 4% of the experimental ones.

Computational Fluid Dynamics↗

NASA NTRS2004

Thermocapillary migration of a drop: an exact solution with Newtonian interfacial rheology and stretching/shrinkage of interfacial area elements for small Marangoni numbers

In this paper we analyze the effects of the following phenomena associated with the thermocapillary migration of a drop. The first is the influence of Newtonian surface rheology of the interface and the second is that of the energy changes associated with stretching and shrinkage of the interfacial area elements, when the drop is in motion. The former occurs because of dissipative processes in the interfacial region, such as when surfactant molecules are adsorbed at the interface in sufficient concentration. The interface is typically modeled in this instance by ascribing to it a surface viscosity. This is a different effect from that of interfacial tension gradients arising from surfactant concentration gradients. The stretching and shrinkage of interfacial area elements leads to changes in the internal energy of these elements that affects the transport of energy in the fluids adjoining the interface. When an element on the interface is stretched, its internal energy increases because of the increase in its area. This energy is supplied by the neighboring fluids that are cooled as a consequence. Conversely, when an element on the interface shrinks, the adjoining fluids are warmed. In the case of a moving drop, elements of interfacial area are stretched in the forward half of the drop, and are shrunk in the rear half. Consequently, the temperature variation on the surface of the drop and its migration speed are modified. The analysis of the motion of a drop including these effects was first performed by LeVan in 1981, in the limit when convective transport of momentum and energy are negligible. We extend the analysis of LeVan to include the convective transport of momentum by demonstrating that an exact solution of the momentum equation is obtained for an arbitrary value of the Reynolds number. This solution is then used to calculate the slightly deformed shape of the drop from a sphere.

Rheology/methods↗

NASA NTRS1996

Interface Morphology During Crystal Growth: Effects of Anisotropy and Fluid Flow

The effect of a parallel shear flow and anisotropic interface kinetics on the onset of instability during growth from a supersaturated solution is analyzed. The model used for anisotropy is based on the microscopic picture of step motion. A shear flow (linear Couette flow or asymptotic suction profile) parallel to the crystal solution interface in the same direction as the step motion decreases interface stability. A shear flow counter to the step motion enhances stability and for sufficiently large shear rates the interface is absolutely morphologically stable. For large wave numbers, the perturbed flow field can be neglected and a simple analytic approximation for the stability-instability demarcation is found.

Coriell, S. R.↗

DOE OSTI2023

A Conservative Finite-Volume Based Interface-Tracking Algorithm Using the Signed Distance Function

Methods for tracking an interface between two fluid phases are developed to ensure desired fluid properties, conservation, and stability are preserved in a finitevolume (FV) discretization. Here, the interface is tracked using a level set method where the signed distance function implicitly defines the interface. Marching methods are used to evaluate the value of the signed distance function, including a novel initialization method to initialize any implicit function to the signed distance function around sharp corners in the level set. Global conservation and consistency with a set of governing equations is enforced by a compression coefficient that measures the volumetric compression or expansion due to inaccuracies in the level set evaluation. A redistribution method is integrated into the volume correction to eliminate the small-cell instability while maintaining global conservation. This suite of methods is implemented and tested using static uniform velocity, and potential flow cases with multiple interface geometries. Results show these methods achieve up to second order accuracy, and are conservative. The application for these methods is intended to track the interface of a 3D printing filament in a finite-volume discretization of the all-speed Navier-Stokes equations.

42 ENGINEERING↗

NASA NTRS1989

Numerical studies of the effects of jet-induced mixing on liquid-vapor interface condensation

Numerical solutions of jet-induced mixing in a partially full cryogenic tank are presented. An axisymmetric laminar jet is discharged from the central part of the tank bottom toward the liquid-vapor interface. Liquid is withdrawn at the same volume flow rate from the outer part of the tank. The jet is at a temperature lower than the interface, which is maintained at a certain saturation temperature. The interface is assumed to be flat and shear-free and the condensation-induced velocity is assumed to be negligibly small compared with radial interface velocity. Finite-difference method is used to solve the nondimensional form of steady state continuity, momentum, and energy equations. Calculations are conducted for jet Reynolds numbers ranging from 150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of above stated parameters on the condensation Nusselt and Stanton numbers which characterize the steady-state interface condensation process are investigated. Detailed analysis to gain a better understanding of the fundamentals of fluid mixing and interface condensation is performed.

Lin, Chin-Shun↗

NASA NTRS1989

Numerical studies of the effects of jet-induced mixing on liquid-vapor interface condensation

Lin, Chin-Shun↗

NASA NTRS2005

Concepts to Automate Fluid Transfer Capability of Low Impact Docking System (LIDS)

The capability to transfer mass between spacecraft is necessary for many mission scenarios. Docking and berthing operations have enabled fluid, electrical, crew and equipment transfers to some degree on all manned space operations since the Gemini program. The Apollo program performed some sophisticated docking maneuvers to land men on the moon and return them safely to Earth. These programs primarily transferred crew, equipment, and pressurized atmosphere between docked spacecraft. The International Space Station (ISS) U.S. modules are connected by Common Berthing Mechanism (CBM) portals. They provide many feed-through ports for electrical, and fluid transfer between modules, as well as a large diameter crew and equipment tunnel. Fluid and electrical jumpers are manually installed after the CBM sealing surfaces have been securely mated to maintain the pressurized cabin environment. CBM berthing and subsequent fluid transfer capability requires a lengthy manual process involving an active interface that mates with a passive half. The Androgynous Peripheral Attach System (MAS) a Russian technology that docked the Russian Zarya module to Unity, or Node 1, is a more complex system that also is capable of fuel transfer, enabling refueling of the Russian re-boost engines on ISS. For several years, a Low Impact Docking System (LIDS) has been under development at Johnson Space Center (JSC). This docking technology has a requirement to be androgynous in order to allow the fabrication of a single configuration that can dock with all other LIDS units. It is desired to make electrical and fluid coupling mating an automated process to enable routine docking and undocking operations to support future exploration missions. It is envisioned that modular design and vehicle assembly will require an efficient LIDS for fuel, electrical, crew, and equipment transfer. Marshall Space Flight Center (MSFC) has joined the LIDS development effort and plans to employ fluid transfer concepts and technology from previous development programs, such as the Automated Fluid Interface System (AFIS) and the Ohital Fluid Transfer System (OFTS) that were developed and tested by MSFC. Orbital Expressderived robotic satellite servicing technology is scheduled to fly soon to demonstrate fluid transfer technology developed for the Air Force. This paper will compare known technology against possible fluid transfer requirements for LIDS to develop design strategies for automated fluid transfer.

Miernik, Janie H.↗

NASA NTRS1991

Vapor condensation at the free surface of an axisymmetric liquid mixed by a laminar jet

This paper presents numerical solutions of jet-induced mixing in a partially full cryogenic tank. An axisymmetric laminar jet is discharged from the central part of the tank bottom toward the liquid-vapor interface. Liquid is withdrawn at the same volume flow rate from the outer part of the tank. The jet is at a temperature lower than the interface, which is maintained at a certain saturation temperature. The interface is assumed to be flat and shear free and the condensation-induced velocity is assumed to be negligibly small compared with radial interface velocity. Finite-difference method is used to slove the nondimensional form of steady-state continuity, momentum, and energy equations. Calculations are conducted for jet Reynolds numbers ranging from 150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of previously stated parameters on the condensation Nusselt and Stanton numbers that characterize the steady-state interface condensation process are investigated. Detailed analysis is performed to gain a better understanding of the fundamentals of fluid mixing and interface condensation.

Lin, Chin-Shun↗

NASA NTRS2007

Translator for Optimizing Fluid-Handling Components

A software interface has been devised to facilitate optimization of the shapes of valves, elbows, fittings, and other components used to handle fluids under extreme conditions. This software interface translates data files generated by PLOT3D (a NASA grid-based plotting-and- data-display program) and by computational fluid dynamics (CFD) software into a format in which the files can be read by Sculptor, which is a shape-deformation- and-optimization program. Sculptor enables the user to interactively, smoothly, and arbitrarily deform the surfaces and volumes in two- and three-dimensional CFD models. Sculptor also includes design-optimization algorithms that can be used in conjunction with the arbitrary-shape-deformation components to perform automatic shape optimization. In the optimization process, the output of the CFD software is used as feedback while the optimizer strives to satisfy design criteria that could include, for example, improved values of pressure loss, velocity, flow quality, mass flow, etc.

Landon, Mark↗

DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗