Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfacial velocity”

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 73 records · Page 4

Experimental Investigation of Flow Orientation Effects on Cryogenic Flow Boiling

Cryogenic flow boiling physics is of paramount importance for future cryogenic space applications employing Cryogenic Fluid Management (CFM) technologies. This experimental study investigates the two-phase flow and heat transfer performance of LN 2 flow boiling across five distinct flow orientations: vertical upflow, vertical downflow, horizontal flow, 45° inclined upflow, and 45° inclined downflow. The study employed the steady- state heating method within a circular heated tube featuring an 8.5-mm inner diameter and a 680-mm heated length. High-speed video recordings were utilized to capture two-phase flow patterns and interfacial behaviors across various flow orientations. The experiments covered a wide range of operating conditions, including mass velocities ranging from 351.80 to 1572.77 kg/m 2 s and inlet pressures from 297.14 to 1032.97 kPa, primarily with near-saturated inlet subcooling. Distinct two-phase flow patterns and regime transitions were identified for each flow orientation. Symmetrical flow patterns were evident in vertical orientations, whereas non-vertical orientations exhibited asymmetric flow stratifications, primarily influenced by the buoyancy force in a terrestrial gravity environment. Bubble dynamic parameters were quantified, and bubble collision and dispersion phenomena were visualized. Analyzing heat transfer performance based on local flow boiling curves and variations in heat transfer coefficients (HTCs), it was observed that vertical upflow demonstrated the most enhanced heat transfer performance, while vertical downflow exhibited the lowest. As mass velocity exceeded 830 kg/ m 2 s, the differences in heat transfer among orientations became less distinct, emphasizing the role of flow inertia in mitigating the influence of flow orientation. A direct comparison of microgravity HTC data from a recent study by the authors against the 1-𝘨 e HTC data indicated enhanced heat transfer performance in microgravity for flow orientation configurations in terrestrial gravity. However, this enhancement progressively diminished as the heat flux increased. Distinct HTC trends were observed as mass velocity increased for different flow orientations and gravity levels.

Cryogens↗

Sensitivity of cellular automata grain structure predictions for high solidification rates

Understanding as-solidified grain structure development during constrained (G 0) alloy solidification conditions is a necessary component of understanding the larger process-microstructure relationships during additive manufacturing (AM) processing. This study applies a cellular automata (CA) model, for static solidification velocities and thermal gradients on the order of those expected in AM melt pools, to investigate grain shape development across ranges of input parameters governing temperature field evolution, nucleation, and growth. Variation in nucleated grain shape from columnar to equiaxed as a function of thermal gradient and nucleation density is modeled with CA and verified against the established transition model of Hunt and Gaumann. Furthermore, while thermal conditions and nucleation density are shown to have the largest roles on microstructure development, the sensitivity of grain structure was very unevenly distributed across parameter space. In particular, it was found that microstructure is very sensitive to mean nucleation undercooling and the interfacial response function in process space regimes near grain structure transitions, but entirely controlled by a combination of thermal gradient, solidification velocity, and nucleation density away from such transitions. The understanding of grain shape sensitivity obtained through CA modeling, particularly at large thermal gradient and nucleation density, will be necessary for accurate microstructure prediction during the extreme processing conditions and large input parameter uncertainty common to AM.

36 MATERIALS SCIENCE↗

Film condensation in a horizontal rectangular duct

Condensation heat transfer in an annular flow regime with and without interfacial waves was experimentally investigated. The study included measurements of heat transfer rate with condensation of vapor flowing inside a horizontal rectangular duct and experiments on the initiation of interfacial waves in condensation, and adiabatic air-liquid flow. An analytical model for the condensation was developed to predict condensate film thickness and heat transfer coefficients. Some conclusions drawn from the study are that the condensate film thickness was very thin (less than 0.6 mm). The average heat transfer coefficient increased with increasing the inlet vapor velocity. The local heat transfer coefficient decreased with the axial distance of the condensing surface, with the largest change at the leading edge of the test section. The interfacial shear stress, which consisted of the momentum shear stress and the adiabatic shear stress, appeared to have a significant effect on the heat transfer coefficients. In the experiment, the condensate flow along the condensing surface experienced a smooth flow, a two-dimensional wavy flow, and a three-dimensional wavy flow. In the condensation experiment, the local wave length decreased with the axial distance of the condensing surface and the average wave length decreased with increasing inlet vapor velocity, while the wave speed increased with increasing vapor velocity. The heat transfer measurements are reliable. And, the ultrasonic technique was effective for measuring the condensate film thickness when the surface was smooth or had waves of small amplitude.

Lu, Qing↗

Combustion Of Moving Droplets And Of Droplets Suspended Within A Convective Environment: Transient Numerical Results

The problem considered is that of a single-component liquid fuel (n-heptane) droplet undergoing evaporation and combustion in a hot, convective, low pressure, zero-gravity environment of infinite expanse. For a moving droplet, the relative velocity (U(sub infinity)) between the droplet and freestream is subject to change due to the influence of the drag force on the droplet. For a suspended droplet, the relative velocity is kept constant. The governing equations for the gas-phase and the liquid-phase consist of the unsteady, axisymmetric equations of mass, momentum, species (gas-phase only) and energy conservation. Interfacial conservation equations are employed to couple the two phases. Variable properties are used in the gas- and liquid-phase. Multicomponent diffusion in the gas-phase is accounted for by solving the Stefan-Maxwell equations for the species diffusion velocities. A one-step overall reaction is used to model the combustion. The governing equations are discretized using the finite volume and SIMPLEC methods. A colocated grid is adopted. Hyperbolic tangent stretching functions are used to concentrate grid points near the fore and aft lines of symmetry and at the droplet surface in both the gas- and liquid-phase. The discretization equations are solved using the ADI method with the TDMA used on each line of the two alternating directions. Iterations are performed within each time-step until convergence is achieved. The grid spacing, size of the computational domain and time-step were tested to ensure that all solutions are independent of these parameters. A detailed discussion of the numerical model is given.

Gogos, George↗

Mechanistic mass transfer in hollow fiber membrane solvent extraction for bio-based isobutanol

Membrane solvent extraction (MSE) has emerged as a promising method for selectively recovering bioproducts from complex aqueous streams. Bio-isobutanol, a next-generation feedstock for biofuel, remains challenging to recover because of its low concentration and the presence of inhibitory substances. This study explores the potential of hollow fiber (HF) MSE for bio-isobutanol recovery and systematically examines the coupled effects of fiber packing, shell-side flow dynamics, and aqueous chemistry on performance. A resistance-in-series model is applied to understand mass transfer in the HF MSE modules, quantify local resistances, and validate overall performance. The results show that increasing the fiber packing provides a larger interfacial area but induces poor flow distribution and channeling, hindering effective isobutanol transport. Meanwhile, increasing the shell-side velocity improves isobutanol recovery due to reductions in the boundary layer thickness. The presence of salts, added to mimic fermentation broth, increases the partition coefficient through salting-out effects, further improving isobutanol flux. A modified correlation for the shell-side mass transfer coefficient (k s,ϕ+v ), integrating geometric and hydrodynamic effects, was developed and validated. The proposed model achieves highly predictive accuracy (r 2 = 0.9808) across a wide range of conditions, outperforming previous models. The findings provide mechanistic insight into the interaction of geometric packing, hydrodynamics, and chemistry in governing mass transfer in HF MSE. Overall, this work demonstrates the potential of HF MSE for efficient bio-isobutanol recovery and also provides practical guidelines on critical factors (packing fraction, partition coefficient, and shell-side velocity), aiding in the design and scaling of MSE systems for resource recovery.

Aqueous chemistry↗

A phase-field method for boiling heat transfer

Here we present a phase field method for heat transfer in two-phase flow with boiling. The vapor/liquid interface evolution is modeled by the Cahn-Hilliard equation. The phase change rate is determined by accounting for the heat conduction balance on either vapor or liquid side of the interfacial area, depending on which side the temperature is assumed to be maintained at the saturation temperature during boiling. The velocity correction scheme proposed by Dong & Shen [27] is extended to solve the Navier-Stokes equations for a non-solenoidal velocity field, and the entropy viscosity method is employed for stabilization. The phase change model is verified by two-dimensional simulations of a vapor bubble growing in super-heated liquid and in film boiling. In both cases, mesh independence of the results is systematically performed. Subsequently, the method is applied to predict the growth of three-dimensional vapor bubble in a rectangular microchannel with boiling flow, achieving good agreement with experimental measurements and available simulation results using the level-set method. The numerical experiments demonstrate that the required mesh resolution for the phase field method is comparable with that of volume of fluid (VoF) and level-set methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Motion driven by the interface

Due to the reduction in buoyant forces aboard orbiting spacecraft such as the Space Shuttle, fluid motion driven by gradients in interfacial tension will be important in the processing of materials in space. In this paper, preliminary results from a study of surface tension driven flow in a pendant drop are reported. The drop is heated from above, and the resulting temperature gradients on the drop surface give rise to interfacial tension gradients. These, in turn, drive a circulation in the drop which is made visible by suitable tracers. The velocities are measured using a video technique, and the data on core velocities are found to agree well with results from a predictive theoretical model.

Jayaraj, K.↗

Implementation of a Drift-Flux Model in SAM for Modeling of Passively Transported GAS in Molten Salt Reactors

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is further developing the Multiphysics Object Oriented Simulation Environment (MOOSE) system thermal hydraulics (T/H) code, System Analysis Module (SAM), to include modeling of molten salt reactor (MSR) designs. MSR designs in which the fuel is dissolved in the salt coolant poses unique modeling and simulation (M&S) challenges because fission products—some of which will be noncondensable gasses—travel with the coolant, thus allowing for deposition and interaction throughout the entire system flow loop. This behavior must be modeled to successfully understand MSR system behavior during normal and off-normal operation events. In past work, SAM has been extended to include a passive species transport capability for modeling of MSR fission products that travel with the salt coolant. This work extends the model to include a drift-flux model to capture the tendency of the gas phase to move at a different velocity than the bulk liquid. Additional constitutive models have also been added for calculation of gas volume fraction, bubble radius, and interfacial area, which will be important for calculation of species phase migration in future coupling activities. Testing has been performed to ensure that the new model agrees with expected values for several simple models, but additional future work will be required to test the model for more complex geometries that will be needed for full-scale MSR simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Velocity and Attenuation of Acoustic Emission Waves in SiC/SiC Composites Loaded in Tension

The behavior of acoustic waves produced by microfracture events and from pencil lead breaks was studied for two different silicon carbide fiber-reinforced silicon carbide matrix composites. The two composite systems both consisted of Hi-Nicalon (trademark) fibers and carbon interfaces but had different matrix compositions that led to considerable differences in damage accumulation and acoustic response. This behavior was primarily due to an order of magnitude difference in the interfacial shear stress for the two composite systems. Load/unload/reload tensile tests were performed and measurements were made over the entire stress range in order to determine the stress-dependence of acoustic activity for increasing damage states. It was found that using the extensional wave velocities from acoustic emission (AE) events produced from pencil lead breaks performed outside of the transducers enabled accurate measurements of the stiffness of the composite. The extensional wave velocities changed as a function of the damage state and the stress where the measurement was taken. Attenuation for AE waveforms from the pencil lead breaks occurred only for the composite possessing the lower interfacial shear stress and only at significantly high stresses. At zero stress after unloading from a peak stress, no attenuation occurred for this composite because of crack closure. For the high interfacial stress composite no attenuation was discernable at peak or zero stress over the entire stress-range of the composite. From these observations, it is believed that attenuation of AE waveforms is dependent on the magnitude of matrix crack opening.

Morscher, Gregory N.↗

Modeling of vapor-liquid interactions in condensing ejectors

Ejectors are compact mechanical devices that utilize the expansion of a high-pressure primary fluid to entrain and compress a low-pressure secondary fluid by means of momentum transfer between the two streams of fluid. Condensing ejectors feature both momentum and heat transfer through the interaction between the vapor and liquid streams. The goal of this paper is to develop and validate a one-dimensional slug flow model to simulate the vapor–liquid interactions in two types of condensing ejectors – one with primary liquid and secondary vapor flows (Type I), and the other with primary vapor and secondary liquid flows (Type II). Control volume analysis of the mass, momentum, and energy balance in each phase and across the interface was conducted for the ejectors. The slug flow models for both types of ejectors are validated with published results in the literature. The friction coefficient on the inner wall of the ejectors and the interfacial heat transfer coefficient are identified as controlling parameters for the ejector performance in terms of pressure and temperature distributions along the axis. The detailed parametric study shows that the liquid inlet velocity and mixing tube diameter have a significant impact on the performance of the Type I ejectors, and the performance of the Type II ejectors is mainly controlled by the secondary liquid flow inlet temperature.

42 ENGINEERING↗

Experimental Heat Transfer Results and Flow Visualization of Vertical Upflow Boiling in Earth Gravity with Subcooled Inlet Conditions – In Preparation for Experiments Onboard the International Space Station

Since 2012, researchers at the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and NASA Glenn Research Center have been collaborating on a long-term effort to study flow boiling and condensation in microgravity. The ultimate goal has been to develop the Flow Boiling and Condensation Experiment (FBCE) for the International Space Station (ISS). Based on the findings from prior flow boiling experiments both at different orientations in Earth gravity and onboard parabolic flights simulating short durations of microgravity, a final refined experiment design, construction, and operating procedure have been arrived at for long-duration microgravity flow boiling experiments onboard the ISS. This study investigates flow boiling of n-Perfluorohexane with subcooled inlet in a rectangular channel of dimensions 114.6 mm heated length, 2.5 mm width, and 5 mm height. These pre-launch experiments (Mission Sequence Testing) were conducted in vertical upflow orientation in Earth gravity using the same experimental rig that was launched to the ISS in August 2021. The various operating parameters varied are heating configuration (single- and double-sided), mass velocity (180 – 3200 kg/m2s), inlet subcooling (+0 – 32°C, encompassing both highly subcooled and near-saturated inlet conditions), and inlet pressure (119 – 191 kPa). High-speed video flow visualization images are presented to explain the two-phase interfacial physics within the channel’s heated section. Heat transfer results in terms of flow boiling curves, streamwise profiles of wall temperature and heat transfer coefficient, and averaged heat transfer coefficients are analyzed and parametric effects elucidated. Severe temporal thermodynamic equilibrium is observed for near-saturated inlet at very low velocities. Nucleate boiling degradation starts at larger heat fluxes for single-sided heating than double sided at low mass velocities with highly subcooled inlet, and conversely at high mass velocities with near-saturated inlet. Nucleate boiling degradation can be delayed to higher heat fluxes by highly subcooling the inlet and increasing mass velocity. The entire local heat transfer coefficient profiles are degraded at higher heat fluxes for near-saturated inlet, but only the downstream part for highly subcooled inlet. This study also confirmed reliability of the upcoming ISS experimental data for subcooled inlet conditions and the collected Earth-gravity data will be used for comparison against the ISS data.

flow boiling↗

Transport Processes in Dendritic Crystallization

Free dentritic growth refers to the unconstrained development of crystals within a supercooled melt, which is the classical dendrite problem. The development of theoretical understanding of dendritic growth and its experimental status is sketched showing that transport theory and interfacial thermodynamics (capillarity theory) are insufficient ingredients to develop a truly predictive model of dendrite formation. The convenient, but incorrect, notion of maximum velocity was used for many years to estimate the behavior of dendritic transformations until supplanted by modern dynamic stability theory. The proper combinations of transport theory and morphological stability seem to be able to predict the salient aspects of dendritic growth, especially in the neighborhood of the tip.

Glicksman, M. E.↗

Free dendritic growth

Free dendritic growth refers to the unconstrained development of crystals within a supercooled melt, which is the classical 'dendrite problem'. Great strides have been taken in recent years in both the theoretical understanding of dendritic growth and its experimental status. The development of this field will be sketched, showing that transport theory and interfacial thermodynamics (capillarity theory) were sufficient ingredients to develop a truly predictive model of dendrite formation. The convenient, but incorrect, notion of 'maximum velocity' was used for many years to estimate the behavior of dendritic transformations until supplanted by modern dynamic stability theory. The proper combinations of transport theory and morphological stability seem to able to predict the salient aspects of dendritic growth, especially in the neighborhood of the tip. The overall development of cast microstructures, such as equiaxed zone formation, rapidly solidified microstructures, etc., also seems to contain additional non-deterministic features which lie outside the current theories discussed here.

Glicksman, M. E.↗

Heat Transfer and Interfacial Flow Physics of Microgravity Flow Boiling in Single-Side-Heated Rectangular Channel with Subcooled Inlet Conditions – Experiments Onboard the International Space Station

This study is the culmination of a long-term collaborative effort between researchers from the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and the NASA Glenn Research Center to investigate gravitational effects on flow boiling and flow condensation. The science and design concepts for this large-scale effort were initiated in 2011 and included several studies detailing various aspects of two-phase fluid physics in both Earth gravity and microgravity, culminating in construction of the large-scale experimental facility named “Flow Boiling and Condensation Experiment (FBCE)”. The experiment was launched to the International Space Station (ISS) in August 2021. Following the successful installation of FBCE, equipped with the Flow Boiling Module (FBM), onboard the ISS and completion of several safety checks, flow boiling experiments were performed for five months from February 2022 until July 2022. This resulted in a large flow boiling database covering broad ranges of operating parameters and heating configurations spanning several research objectives. This study investigates microgravity flow boiling of n-perfluorohexane with subcooled inlet in a single-side-heated rectangular channel of dimensions 114.6-mm heated length, 2.5-mm heated width, and 5.0-mm height. Key operating parameters investigated are mass velocity (199.90 – 3200.13 kg/m 2 s), inlet subcooling (0.10 – 45.76°C), and inlet pressure (113.30 – 164.29 kPa). Images and image sequences acquired via high-speed-video are presented to elucidate the interfacial flow physics. To analyze and explain the effects of various parameters in microgravity, heat transfer results are presented as flow boiling curves, streamwise profiles of wall temperature and heat transfer coefficient, and parametric trends of local and averaged heat transfer coefficient. Mass velocity and inlet subcooling significantly influenced most of the aforementioned aspects of flow boiling, whereas effects of inlet pressure were comparatively insignificant. Although the data and observed flow physics might be different, the parametric effects and trends in microgravity are similar to vertical upflow in Earth gravity. Some cases, especially low mass velocities, high heat fluxes, and large degrees of inlet subcooling, experienced temporally anomalous flow behaviors caused by two-phase flow instabilities manifesting as flow reversals and resulted in deviations in overall trends. Severe thermodynamic non-equilibrium is observed throughout the channel. Overall, FBCE’s ISS experiments were successful for subcooled inlet with single-sided heating of rectangular channel, and the collected data well established the various effects on flow boiling physics in highly controlled long-term microgravity conditions.

flow boiling↗

Flow Visualization, Heat Transfer, and Critical Heat Flux of Flow Boiling in Earth Gravity with Saturated Liquid-Vapor Mixture Inlet Conditions – In Preparation for Experiments Onboard the International Space Station

This study investigates flow boiling of n-Perfluorohexane with saturated two-phase mixture inlet in a rectangular channel of dimensions 114.6-mm heated length, 2.5-mm width, and 5 mm-height. The experiments were performed as part of the Mission Sequence Testing of the Flow Boiling and Condensation Experiment’s (FBCE) Flow Boiling Module (FBM) in the vertical upflow configuration in Earth gravity using the same experimental system that was launched to the International Space Station (ISS) in August 2021. The operating parameters varied are heating configuration (single- and double-sided), mass velocity (380 – 2400 kg/m2s), inlet quality (0.011 – 0.519), and inlet pressure (120 – 179 kPa). High-speed video photographs are presented to explain the two-phase flow patterns within the channel’s heated length. Flow patterns are constituted by low-density and high-density fronts moving along the channel, with the high-density fronts gradually reducing in length due to evaporation. Heat transfer results in terms of flow boiling curves, streamwise wall temperature profiles, streamwise heat transfer coefficient profiles, and average heat transfer coefficients are presented and trends discussed. CHF data from the present experiments are combined with prior databases to compile a consolidated FBCE-CHF database for saturated inlet to expand the ranges of operating conditions and include other flow orientations in Earth gravity. Experimental CHF trends are also discussed. The interfacial lift-off model shows a good CHF predictive accuracy evidenced by a mean absolute error of 11.97% for this consolidated database after constraining it to mass velocities greater than or equal to 500 kg/m2s. Finally, this study confirmed reliability of the upcoming ISS experiments for saturated inlet conditions and the collected Earth-gravity data will be compared to ISS microgravity data.

Flow boiling↗

Continuum shock mixture models for Ni+Al multilayers: Inert mesoscale simulations

Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2 × 6 μm 2 ) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ = 15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.

Kittell, David E. [Sandia National Laboratories (S↗

Optical Studies of model binary miscibility gap system

In order to develop a better understanding of separation processes in binary miscibility gap metal alloys, model transparent fluid systems were studied. The system selected was diethylene glycol-ethyl salicylate which has convenient working temperatures (288 to 350 K), low toxicity, and is relatively easy to purify. The system is well characterized with respect to its phase diagram, density, surface and interfacial tensions, viscosity and other pertinent physical properties. Studies of migration of the dispersed phase in a thermal gradient were performed using conventional photomicroscopy. Velocities of the droplets of the dispersed phase were measured and compared to calculated rates which included both Stokes and thermal components. A holographic microscopy system was used to study growth, coalescence, and particle motions. Sequential holograms allowed determination of particle size distribution changes with respect to time and temperature. Holographic microscopy is capable of recording particle densities up to 10 to the 7th power particles/cu cm and is able to resolve particles of the order of 2 to 3 microns in diameter throughout the entire volume of the test cell. The reconstructed hologram produces a wavefront that is identical to the original wavefront as it existed when the hologram was made. The reconstructed wavefront is analyzed using a variety of conventional optical methods.

Lacy, L. L.↗

Study of the film boiling heat transfer and two-phase flow interface behavior using image processing

Here, in this study, based on a small-scale quench test facility, the two-phase flow interface behavior during quench transients is visualized and analyzed utilizing an image processing framework. The high-fidelity experimental results obtained for two-phase flow in the current framework can support various studies both in the time domain and in the frequency domain. In particular, visualization of the data obtained from different heating surfaces under different test conditions are used to perform a full-scale transient 2-D vapor film reconstruction. The liquid-vapor interface variations in various heat transfer regimes as well as at the initial film breakup point can be directly obtained through the processed data. Moreover, the temporal variation of the interfacial wave frequency approaching quench is investigated in detail. Based on the high-resolution data obtained for the liquid-vapor interface, the detailed phase velocity and temperature profiles are obtained through theoretic analysis, based on which the film boiling heat transfer coefficient (HTC) can be determined. In addition, an improved film boiling HTC model is developed considering the effects of wall superheat, liquid subcooling temperature, vapor film thickness as well as fluid properties. The model is found to predict film boiling HTC well within 15% error.

42 ENGINEERING↗