SEARCH · Search NASA

Results for “interfacial heat transfer”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Numerical analysis of bipropellant combustion in orbit maneuvering vehicle thrust chamber

The combustion and aerodynamic processes of liquid-liquid bipropellant; monomethylhydrazine and nitrogen tetroxide in the combustion chamber of a variable thrust engine (VTE) have been examined to assess the flow structure, combustion characteristics, and wall heat transfer under the simulated engine operations at two power levels. The Generalized Elliptic Multi-phase flow with Chemical reaction, Heat transfer and Interfacial Processes (GEM-CHIP) code and the Bipropellant Combustion (BICOMB) code, which is the GEMCHIP algorithm reformulated in the frame work of a body fitted coordinate system, have been used for the numerical study. Combustion flow structure with non-premixed type injection that simulates the impingement type injector is predicted and the results are compared with that of the premixed injector.

Chiu, H. H.↗

NASA NTRS2014

CVB: the Constrained Vapor Bubble Capillary Experiment on the International Space Station MARANGONI FLOW REGION

The Constrained Vapor Bubble (CVB) is a wickless, grooved heat pipe and we report on a full- scale fluids experiment flown on the International Space Station (ISS). The CVB system consists of a relatively simple setup a quartz cuvette with sharp corners partially filled with either pentane or an ideal mixture of pentane and isohexane as the working fluids. Along with temperature and pressure measurements, the two-dimensional thickness profile of the menisci formed at the corners of the quartz cuvette was determined using the Light Microscopy Module (LMM). Even with the large, millimeter dimensions of the CVB, interfacial forces dominate in these exceedingly small Bond Number systems. The experiments were carried out at various power inputs. Although conceptually simple, the transport processes were found to be very complex with many different regions. At the heated end of the CVB, due to a high temperature gradient, we observed Marangoni flow at some power inputs. This region from the heated end to the central drop region is defined as a Marangoni dominated region. We present a simple analysis based on interfacial phenomena using only measurements from the ISS experiments that lead to a predictive equation for the thickness of the film near the heated end of the CVB. The average pressure gradient for flow in the film is assumed due to the measured capillary pressure at the two ends of the liquid film and that the pressure stress gradient due to cohesion self adjusts to a constant value over a distance L. The boundary conditions are the no slip condition at the wall interface and an interfacial shear stress at the liquid- vapor interface due to the Marangoni stress, which is due to the high temperature gradient. Although the heated end is extremely complex, since it includes three- dimensional variations in radiation, conduction, evaporation, condensation, fluid flow and interfacial forces, we find that using the above simplifying assumptions, a simple successful model can be developed.

heat pipes↗

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 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 VOF and Sharp Interface multiphase 2Daxisymmetric 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. These parameters include: the effect of turbulence modeling and turbulence damping at the interface; the effect of constant vs. temperature dependent fluid properties; the effect of accommodation coefficient; the effect of modeling phase change at the interface and the effect of the fill level.

Computational Fluid Dynamics↗

NASA NTRS2015

The Constrained Vapor Bubble Experiment - Interfacial Flow Region

Internal heat transfer coefficient of the CVB correlated to the presence of the interfacial flow region. Competition between capillary and Marangoni flow caused Flooding and not a Dry-out region. Interfacial flow region growth is arrested at higher power inputs. 1D heat model confirms the presence of interfacial flow region. 1D heat model confirms the arresting phenomena of interfacial flow region Visual observations are essential to understanding.

microgravity↗

NASA NTRS2004

Boiling Heat Transfer Mechanisms in Earth and Low Gravity: Boundary Condition and Heater Aspect Ratio Effects

Kim, Jungho↗

NASA NTRS2002

Length Scale and Gravity Effects on Microgravity Boiling Heat Transfer

Boiling is a complex phenomenon where hydrodynamics, heat transfer, mass transfer, and interfacial phenomena are tightly interwoven. An understanding of boiling and critical heat flux in microgravity environments is of importance to space based hardware and processes such as heat exchange, cryogenic fuel storage and transportation, electronic cooling, and material processing due to the large amounts of heat that can be removed with relatively little increase in temperature. Although research in this area has been performed in the past four decades, the mechanisms by which heat is removed from surfaces in microgravity are still unclear. In earth gravity, buoyancy is an important parameter that affects boiling heat transfer through the rate at which bubbles are removed from the surface. A simple model describing the bubble departure size based on a quasistatic force balance between buoyancy and surface tension is given by the Fritz [I] relation: Bo(exp 1/2) = 0.0208 theta where Bo is the ratio between buoyancy and surface tension forces. For small, rapidly growing bubbles, inertia associated with the induced liquid motion can also cause bubble departure. In microgravity, the magnitude of effects related to natural convection and buoyancy are small and physical mechanisms normally masked by natural convection in earth gravity such as Marangoni convection can substantially influence the boiling and vapor bubble dynamics. CHF (critical heat transfer) is also substantially affected by microgravity. In 1 g environments, Bo has been used as a correlating parameter for CHF. Zuber's CHF model for an infinite horizontal surface assumes that vapor columns formed by the merger of bubbles become unstable due to a Helmholtz instability blocking the supply of liquid to the surface. The jets are spaced lambda(sub D) apart, where lambda(sub D) = 2pi square root of 3[(sigma)/(g(rho(sub l) - rho(sub v)](exp 1/2) = 2pi square root of 3 L Bo(exp -1/2) = square root of 3 lambda(sub c) and is the wavelength that amplifies most rapidly. The critical wavelength, lambda(sub c), is the wavelength below which a vapor layer underneath a liquid layer is stable. For heaters with Bo smaller than about 3 (heaters smaller than lambda(sub D)), the above model is not applicable, and surface tension effects dominate. Bubble coalescence is thought to be the mechanism for CHF under these conditions. Small Bo can result by decreasing the size of a heater in earth gravity, or by operating a large heater in a lower gravity environment. In the microgravity of space, even large heaters can have low Bo, and models based on Helmholtz instability should not be applicable. The macrolayer model of Haramura and Katto is dimensionally equivalent to Zuber's model and has the same dependence on gravity, so it should not be applicable as well. The goal of this work is to determine how boiling heat transfer mechanisms in a low-g environment are altered from those at higher gravity levels. Boiling data using a microheater array was obtained under gravity environments ranging from 1.8 g to 0.02 g with heater sizes ranging from 2.7 mm to 1 mm. The boiling behavior for 2.7 mm at 0.02 g looked quite similar to boiling on the 1 mm heater at 1 g-the formation of a large primary bubble surrounded by smaller satellite bubbles was observed under both conditions. The similarity suggests that for heaters smaller than some fraction of I(sub c), coalescence and surface tension dominate boiling heat transfer. It also suggests that microgravity boiling can be studied by studying boiling on very small heaters.

Kim, Jungho↗

DOE OSTI2021

Enhancing Steam-Side Heat Transfer via Microdroplet Ejection using Inorganic Coatings

During condensation of water on a superhydrophobic nanostructured surfaces, droplets coalesce and jump (droplet ejection) due to the release of excess surface energy. Meticulously designed nanostructured surfaces or coatings capable of droplet ejection can enhance heat transfer performance by easy removal of droplets during condensation. In the present work, the capabilities of the Nelumbo droplet ejection coatings technology currently used in commercial heat exchangers are explored and optimized for materials and conditions relevant to steam condensers. Specifically, this includes (1) adapting and optimizing the procedures and chemistries to deposit droplet ejecting coatings on materials used in the steam condenser, (2) testing of the heat transfer, durability, and fouling of the fabricated coatings under a variety of steam condenser relevant conditions, and (3) modeling of the impact of droplet dynamics and interfacial properties on heat transfer in steam condensers. We experimentally demonstrated that at low supersaturation conditions (S~1.02), developed superhydrophobic nanostructured surface gives a 40% higher peak heat transfer flux and a 17% higher peak overall heat transfer coefficient (U) with droplet ejection mechanism compared to the dropwise condensation mechanism on the same nanostructured surface. Though the experimentally observed droplet ejection mechanism was short-lived (~3 min) due to the limitation of the chiller to maintain a constant supersaturation condition in fabricated steam condenser. And large variations in water inlet temperature (+ to - 2 degrees C) caused flooding of nanostructure (Wenzel state) at higher supersaturation (S > or = 1.12). This is an important finding because typical supersaturation levels in field operation are over 1.8 and as such, either the flooding potential for these materials should be improved or operation at lower supersaturation may be required to maximize heat transfer efficiency. If supersaturation conditions can be maintained reliably in the steam condenser in addition to droplet ejection mechanism during condensation, this could provide a potential 2% reduction in water flow rate which equates to a savings of over 3900 gallons of water per minute for a 500 MW steam turbine and over $0.3MM savings annually given a 0.02 cent marginal water withdrawal cost. It is also worth noting that these savings are based on improvement of droplet ejection coatings over coatings capable of performing consistent dropwise mechanism in in-field surface condensers. Savings due to the droplet ejections coatings when compared to bare substrates could be much higher. Additionally, the present work provides the importance of steam condenser operating conditions and insights into the challenges in modeling the condensing droplet dynamics on surfaces.

20 FOSSIL-FUELED POWER PLANTS↗

NASA NTRS2015

TankSIM: A Cryogenic Tank Performance Prediction Program

Accurate prediction of the thermodynamic state of the cryogenic propellants in launch vehicle tanks is necessary for mission planning and successful execution. Cryogenic propellant storage and transfer in space environments requires that tank pressure be controlled. The pressure rise rate is determined by the complex interaction of external heat leak, fluid temperature stratification, and interfacial heat and mass transfer. If the required storage duration of a space mission is longer than the period in which the tank pressure reaches its allowable maximum, an appropriate pressure control method must be applied. Therefore, predictions of the pressurization rate and performance of pressure control techniques in cryogenic tanks are required for development of cryogenic fluid long-duration storage technology and planning of future space exploration missions. This paper describes an analytical tool, Tank System Integrated Model (TankSIM), which can be used for modeling pressure control and predicting the behavior of cryogenic propellant for long-term storage for future space missions. It is written in the FORTRAN 90 language and can be compiled with any Visual FORTRAN compiler. A thermodynamic vent system (TVS) is used to achieve tank pressure control. Utilizing TankSIM, the following processes can be modeled: tank self-pressurization, boiloff, ullage venting, and mixing. Details of the TankSIM program and comparisons of its predictions with test data for liquid hydrogen and liquid methane will be presented in the final paper.

Bolshinskiy, L. G.↗

DOE OSTI2021

Flow boiling in geometrically modified microchannels

Microchannels are a promising solution for high-heat-flux thermal management scenarios, including high-performance microelectronics cooling and power electronics cooling. However, thermohydraulic instabilities result from the rapid vapor bubble formation. The prior literature has examined several methods, including constricted inlet microchannels, expanding microchannels, and auxiliary jetting microchannels, to mitigate the effect of these instabilities. Computational fluid dynamics and heat transfer (CFD/HT) modeling of the flow boiling phenomena in these microchannel configurations has seen limited examination, and one-to-one numerical comparisons of the different mitigation strategies have not been performed. In the present investigation, CFD/HT analyses using a three-dimensional (3D) volume of fluid model coupled with a phase-change model for the interfacial heat and mass transfer were performed for multiple microchannel configurations (constricted inlet, expanding, and auxiliary jetting microchannels). A benchmark case of a rectangular microchannel was examined to quantify baseline thermohydraulic performance. Results demonstrated slight to significant thermal performance improvements for all cases, and significant pressure benefits for the expanding and jetting cases, consistent with experimental results in the literature. Bubble dynamics and visualization for the baseline and alternative configurations are provided to give insight into their underlying physics, and the differences in performance were investigated and compared with available literature.

Broughton, Justin (ORCID:000000026848457X)↗

NASA NTRSDate not supplied

Mastering the Cryogenic Frontier: Predicting Sloshing in Cryogenic Propellant Tanks

In recent years, NASA has undertaken the development and validation of computation fluid dynamic (CFD) models for predicting cryogenic two-phase sloshing within propellant tanks. These models have been rigorously anchored against cryogenic ground tests that induce sloshing, as well as limited microgravity tests. The primary focus of these models is to accurately predict the two-phase condensation and, at times, evaporation within the cryogenic propellant. This approach is important for estimating the required autogenous and/or helium pressurant for completing propulsion system burns or ensuring adequate pressure for in-space propellant transfers. Unlike in hypergolic propulsion systems, where the propellant is stored at room temperature, cryogenic propellant is highly volatile. Therefore, ensuring the accurate estimation of pressurant necessitates the precise modeling of two-phase interfacial heat and mass transfer processes. This paper aims to provide an overview of NASA’s continuing efforts in developing and validating these CFD models, with particular emphasis on benchmark cases with sloshing that enhances phase change, tracing their historical development, and demonstrating their predictive accuracy when compared to test results.

low-gravity fluid modeling↗

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↗

NASA NTRS2021

Experimental and Computational Investigation of Flow Boiling in Microgravity

This study explores use of Computational Fluid Dynamics (CFD) to predict near-saturated flow boiling of FC-72 in microgravity. The computational method employs transient analysis to predict detailed interfacial behavior and heat transfer characteristics along a rectangular channel heated along two opposite walls. Predicted results are validated against experimental temperature measurements and high-speed video images captured during a series of parabolic aircraft maneuvers for three sets of operating conditions which include variations of both mass velocity and wall heat flux. The computational method is based on the multi-phase volume of fluid (VOF) model, which is combined with appropriate phase change and turbulence models, and accounts for both shear-lift force on bubbles and conjugate heat transfer along the heating walls. A key advantage of the CFD method is ability to capture details that are very difficult to measure experimentally, including detailed spatial variations of bubble shape, void fraction, mixture fluid temperature, liquid velocity, and vapor velocity, results for which are presented for each of the three test cases. Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with experiment. Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations. Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method. Predicted wall temperature is fairly uniform in the middle of the heated length but increases in the entrance region, due to sensible heat transfer in the subcooled liquid, and decreases toward the exit, mostly because of flow acceleration resulting from increased void fraction.

two-phase cooling; flow boiling; CFD; microgravity↗

DOE OSTI2022

Numerical simulation of evaporating wavy falling liquid films in laminar gas streams

Numerical simulations are performed to investigate the interfacial heat and mass transfer from evaporating wavy falling liquid films in interaction with laminar gas streams. The OpenFOAM solver has been used to conduct the simulations where the liquid-gas interface is resolved using the Volume of Fluid method of solving for three phases, i.e., liquid, vapor, and air. The configuration considered is a falling liquid (water) film on a heated vertical plate with a confined laminar moist air (gas) flow that is either (a) co-current or (b) counter-current to the downward liquid flow. The evaporation at the liquid-gas interface is driven by the interfacial gradient of the vapor mass fraction. Interfacial waves are triggered using a monochromatic forcing disturbance that leads to sinusoidal or solitary waves forming at the liquid-gas interface under respective forcing frequencies. Further, the numerical model is validated with the available experimental data. The results show nearly a 15% enhancement in time-averaged Sherwood number (Sh) due to film waviness (sinusoidal or solitary) at the lower volumetric gas flow rate, Q g = +50 (co-current) and Q g = -50 (counter-current). This enhancement in the Sh for both the waves further increases by 11% with Q g = +800 and 196% with Q g = -800. A closer examination of the mass transfer process over a wave demonstrates that with Q g = +50, the concentration of the gas side streamlines at the trough locations of the wave leads to higher values of Sh at these locations. However, with Q g = +800, although the overall Sh increases, vortices appear at the wave trough locations, leading to a corresponding decrease in the local Sh values. Correlations are proposed for predicting Sh under co-current and counter-current gas flow effects.

36 MATERIALS SCIENCE↗

DOE OSTI2021

A two-phase three-field modeling framework for heat pipe application in nuclear reactors

Heat pipes and two-phase thermosyphons are highly efficient heat transfer devices utilizing continuous evaporation and condensation of working fluid for two-phase heat transport in closed systems. Because of the nearly isothermal and fully passive phase-change heat transfer mechanism, heat pipes and thermosyphons have found many applications in nuclear engineering, space technologies, and other energy systems. High-temperature heat pipes are used in nuclear microreactors to remove fission power from the primary system and are coupled with power conversion systems or process heat applications. Modeling of the two-phase flow phenomena inside a heat pipe is essential to its design and safety analysis. In this study, a comprehensive one-dimensional two-phase three-field flow model has been developed for the analysis of heat pipes in normal operation conditions and transients. The conservation or field equations of mass, momentum, and energy were developed for the liquid film, vapor, and droplet. In addition, constitutive models or correlations were reviewed thoroughly and provided for the closure of the three-field equations. Specific constitutive equations regarding interfacial mass and heat transfer at two interfaces, namely film-gas interface and gas-droplet interface, were reviewed for droplet entrainment and deposition rates as well as film and droplet evaporation rates. Furthermore, mechanistic correlations of annular flow film thickness were recommended for the modeling of the thermosyphons without a wick as a critical constitutive correlation. Furthermore, experimental data needs from new experiments using a prototype working fluid or surrogate fluids for the model validation of high-temperature heat pipes in microreactors were recommended for future research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1984

Morphological stability and fluid dynamics of vapor crystal growth

Research on morphological stability and fluid dynamics of crystal growth is discussed. Interfacial heat and mass transfer research is discussed. The finding of surface roughening is a precursor to a solid-solid phase transition was further quantified. Progress was obtained with the mass spectroscopic characterization of GeSe-Ge I sub 4.

Rosenberger, F. E.↗

NASA NTRS1993

The roles of fluid motion and other transport phenomena in the morphology of materials

Two crystallization problems were studied: the growth of protein crystals, in particular the influence of colloidal forces and convection, and the influence of interface resistance on the growth of dendritic crystals. The protein study involved both experimental and theoretical work; the work of dendrites was entirely theoretical. In the study of protein crystallization, experiments were carried out where crystals were grown in the presence and absence of natural convection. No evidence was found that convection retards crystal growth. The theoretical study focused on the influence of colloidal forces (electrostatic and London-van der Waals) on the interaction between a protein molecule and a flat crystal surface. It was shown that the interaction is extremely sensitive to colloidal forces and that electrostatic interactions play a strong role in deciding whether or not a molecule will find a favorable site for adsorption. In the study of dendritic growth, the role of an interfacial resistance on the selection processes was examined. Using a computational scheme, it was found that the selected velocity is strongly dependent on the magnitude of the interfacial resistance to heat transfer. This is a possible explanation for discrepancies between the theoretical and experimental results on succinonitrile.

Saville, D. A.↗

NASA NTRS1996

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗