Film-stability in a vertical rotating tube with a core-gas flow
Linear hydrodynamic stability of interface between Newtonian liquid film and core fluid under influence of swirl, core flow, and gravity
SEARCH · Search NASA
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.
Linear hydrodynamic stability of interface between Newtonian liquid film and core fluid under influence of swirl, core flow, and gravity
During testing on the E-1 Test Stand, the High Pressure (HP) Liquid Hydrogen (LH) vent line was displaced off its supports, damaging the expansion joints (EJ), piping, and pipe supports of the line. The damage to the vent line is believed to have occurred during Cell 3 testing on Thursday, September 21, during the first test of the day. The fuel systems being used for Cell 3 testing are connected to the LP and HP LH Vent systems and are being used in Liquid Natural Gas (LNG) service. The Cell 3 interface bleed valve that discharges to the HP LH Vent system was inadvertently left at 100% open during the chill-in process for longer than the typical amount of time. Since this valve is usually at a lower set point or intermittently closed during the chill-in process, more than the typical amount of LNG flowed into the vent line as pre-test procedures continued as normal. The open bleed valve was noticed before test preparations were completed and prior to refilling of the HP LH Tank before test. The test sequence was started, run, and shutdown normally. Video of the flare stack showed an atypical discharge of liquid at the flare stack exit. Pressure downstream of the Cell 3 interface bleed valve (a valve connected to the southernmost part of the vent system) spiked to a higher-than-normal level. No evidence of damage was captured at the time of the event. The system was used in a typical fashion after the event and for 4 subsequent tests on Thursday and Friday, before the damage was observed. A drive-by inspection of the flare stack was performed on Friday (September 22nd) morning before sunrise and before test activities began, but no damage of the vent line was observed at that time. No abnormal system behavior was observed during testing on Friday. The damage was discovered and reported to NASA Test Operations on Monday, initiating the Incident Response Team. The Mishap appointment letter was issued September 27th to formally initiate the Mishap Investigation.
The trapping and behavior of gas bubbles were studied during low-gravity solidification of carbon tetrabromide, a transparent metal-model material. The experiment was performed during a NASA-sponsored sounding rocket flight and involved gradient freeze solidification of a gas-saturated melt. Gas bubbles were evolved at the solid-liquid interface during the low-gravity interval. No large-scale thermal migration of bubbles, bubble pushing by the solid-liquid interface, or bubble detachment from the interface were observed during the low-gravity experiment. A unique bubble motion-fluid flow event occurred in one specimen: a large bubble moved downward and caused some circulation of the melt. The gas bubbles that were trapped by the solid in commercial-purity material formed voids that had a cylindrical shape, in contrast to the spherical shape that had been observed in a prior low-gravity experiment. These shapes were not influenced by the gravity level (0.0001 g-0 vs g-0), but were dependent upon the initial temperature gradient. In higher purity material, however, the shape of the voids changed from cylindrical in 1g to spherical in low gravity.
The trapping and behavior of gas bubbles were studied during low gravity solidification of carbon tetrabromide. The flight experiments were performed during two sounding rocket flights (SPAR 1 and SPAR 3) and involved gradient freeze solidification of gas saturated melts. Gas bubbles were evolved at the solid-liquid interfaces during the low gravity intervals. No large-scale thermal migration of bubbles, bubble pushing by the solid-liquid interface, or bubble detachment from the interface were observed during the low gravity experiments. During the SPAR 3 experiment, a unique bubble motion-fluid flow event occurred in one specimen: a large bubble moved downward and caused some circulation of the melt. The gas bubbles that were trapped by the solid in commercial purity material formed voids that had a cyclindrical shape in SPAR 3, in contrast to the spherical shape that had been observed in SPAR 1. These shapes were not influenced by the gravity level, but were dependent upon the initial temperature gradient. In higher purity material the shape of the voids changed from cylindrical in one-g to spherical in low-g.
We report nonthermal plasmas in contact with liquids have been shown to generate a variety of reactive species capable of initiating reduction–oxidation (redox) reactions at the electrochemically active plasma–liquid interface. In conventional electrochemical cells, selective redox chemistry is achieved by controlling the reduction potential at the solid electrode–electrolyte interface by applying a bias via an external circuit. In the case of plasma–liquid systems, an analogous means of tuning the reduction potential near the interface has not clearly been identified. When treated as a floating surface, the liquid is expected to adopt a net negative charge to balance the flux of hot electrons and relatively cold positive ions. The reduction potential near the plasma–liquid interface is hypothesized to be proportional to the floating potential, which can be approximated using an analytical model provided the plasma parameters are known. Herein, we present a framework for correlating the electron density and electron temperature of a noble gas plasma jet to the reduction potential near the plasma–liquid interface. The plasma parameters were acquired for an argon atmospheric plasma jet in contact with an aqueous solution by means of laser Thomson scattering. The reduction potential was determined using identical reference electrodes to measure the potential difference between the plasma–liquid interface and bulk solution. Interestingly, the measured reduction potentials near the plasma–liquid interface were found to be in good agreement with the model-predicted values determined using the plasma parameters obtained from the Thomson scattering experiments.
A method was developed for determining the surface tension coefficient at the interface of two immiscible liquids or of a gas and a liquid. Basic in this method is the measurement of the capillary rise of a droplet of one of the phases placed on the surface of the other phase. It is shown that the capillary rise is a function of the surface energy of the droplet and the viscous friction of the phases. A mathematical model of the dynamics of a liquid droplet is supported by experiments.
This paper describes a novel method of coupling the atomization and spray combustion processes encountered in coaxial injection elements of liquid rocket engine thrust chambers. This method is based on the Jet-Embedding technique in which the liquid jet core equations and the gas phase equations are solved separately. The liquid and gas phase solutions, however, are coupled through the boundary conditions at the interface between the phases. The computational grid for the gas phase calculations are adapted to the shape of the liquid jet core. The axial variation of droplet sizes are calculated using a stability analysis appropriate for the atomization regime of liquid jet breakup. The predictions of this method have been validated with experimental data on low speed water jets. Using this method, calculations are performed for the SSME fuel preburner single injector flow field. The results obtained are in good agreement with the predictions of the volume-of-fluid method.
A major hurdle in utilizing carbon dioxide (CO 2 ) lies in separating it from industrial flue gas mixtures and finding suitable storage methods that enable its application in various industries. To address this issue, we utilized a combination of molecular dynamics simulations and experiments to investigate the behavior of CO 2 in common room-temperature ionic liquids (RTIL) when in contact with aqueous interfaces. Here, our investigation of RTILs, [EMIM][TFSI] and [OMIM][TFSI], and their interaction with a pure water layer mimics the environment of a previously developed ultrathin enzymatic liquid membrane for CO 2 separation. We analyzed diffusion constants and viscosity, which reveals that CO 2 molecules exhibit faster mobility within the selected ILs compared to what would be predicted solely based on the viscosity of the liquids using the standard Einstein–Stokes relation. Moreover, we calculated the free energy of translocation for various species across the aqueous–IL interface, including CO 2 and HCO 3 – . Free energy profiles demonstrate that CO 2 exhibits a more favorable partitioning behavior in the RTILs compared to that in pure water, while a significant barrier hinders the movement of HCO 3 – from the aqueous layer. Experimental measurement of the CO 2 transport in the RTILs corroborates the model. These findings strongly suggest that hydrophobic RTILs could serve as a promising option for selectively transporting CO 2 from aqueous media and concentrating it as a preliminary step toward storage.
Evaporation is a ubiquitous and complex phenomenon of importance to many natural and industrial systems. Evaporation occurs when molecules near the free interface overcome intermolecular attractions with the bulk liquid. As molecules escape the liquid phase, heat is removed, causing evaporative cooling. The influence of evaporative cooling on inducing a temperature difference with the surrounding atmosphere as well as within the liquid is poorly understood. Here, we develop a technique to overcome past difficulties encountered during the study of heterogeneous droplet evaporation by coupling a piezo-driven droplet generation mechanism to a controlled micro-thermocouple to probe microdroplet evaporation. The technique allowed us to probe the gas-phase temperature distribution using a micro-thermocouple (50 μm) in the vicinity of the liquid–vapor interface with high spatial (±10 μm) and temporal (±100 ms) resolution. We experimentally map the temperature gradient formed surrounding sessile water droplets having varying curvature dictated by the apparent advancing contact angle (100° ≲ θ ≲ 165°). The experiments were carried out at temperatures below and above ambient for a range of fixed droplet radii (130 μm ≲ R ≲ 330 μm). Our results provide a primary validation of the centuries-old theoretical framework underpinning heterogeneous droplet evaporation mediated by the working fluid, substrate, and gas thermophysical properties, droplet apparent contact angle, and droplet size. We show that microscale droplets residing on low-thermal-conductivity substrates such as glass absorb up to 8× more heat from the surrounding gas compared to droplets residing on high-thermal-conductivity substrates such as copper. Finally, our work not only develops an experimental understanding of the heat transfer mechanisms governing droplet evaporation but also presents a powerful platform for the study and characterization of liquid–vapor transport at curved interfaces wetting and nonwetting advanced functional surfaces.
Porous liquids (PLs) are an attractive material for gas separation and carbon sequestration due to their permanent internal porosity and high adsorption capacity. PLs that contain zeolitic imidazole frameworks (ZIFs), such as ZIF-8, form PLs through exclusion of aqueous solvents from the framework pore due to its hydrophobicity. The gas adsorption sites in ZIF-8 based PLs are historically unknown; gas molecules could be captured in the ZIF-8 pore or adsorb at the ZIF-8 interface. To address this question, ab initio molecular dynamics was used to predict CO 2 binding sites in a PL composed of a ZIF-8 particle solvated in a water, ethylene glycol, and 2-methylimidazole solvent system. Further, the results show that CO 2 energetically prefers to reside inside the ZIF-8 pore aperture due to strong van der Waals interactions with the terminal imidazoles. However, the CO 2 binding site can be blocked by larger solvent molecules that have greater adsorption interactions. CO 2 molecules were unable to diffuse into the ZIF-8 pore, with CO 2 adsorption occurring due to binding with the ZIF-8 surface. Therefore, future design of ZIF-based PLs for enhanced CO 2 adsorption should be based on the strength of gas binding at the solvated particle surface.
Developing efficient energy storage systems that capture and convert CO 2 is critical for mitigating carbon emissions. Here, we report a Na–CO 2 battery with ruthenium dioxide (RuO 2 ) cathode catalysts and propane-1,3-diamine (PDA) as an electrolyte additive to enhance CO 2 capture and conversion efficiency. The integration of CO 2 adsorption and electrochemical reduction facilitates activation of the inert CO 2 molecule and circumvents gas–solid–liquid ternary-phase reactions at the interface. We employed density functional theory (DFT) calculations to systematically unravel the reaction mechanisms and energetics governing CO 2 reduction, both with and without PDA. Our results reveal an energetically favorable pathway toward the formation of Na 2 CO 3 and C as final discharge products, rather than sodium oxalate (Na 2 C 2 O 4 ). The CO 2 –amine adduct facilitates charge transfer from PDA to CO 2 , which results in activation of CO 2 . The kinetics of CO 2 conversion and regeneration of PDA were found to be significantly enhanced on the RuO 2 surface compared to the bulk electrolyte. More importantly, pre-activation of CO 2 via the amine–CO 2 adduct lowers the total overpotential to 2.44 V, compared to 3.13 V without PDA. This study provides fundamental insights into CO 2 electroreduction in Na–CO 2 batteries and underscores the promise of electrolyte engineering for sustainable CO 2 utilization and high-performance energy storage.
The influence of gravity on the two-phase flow dynamics is obvious.As the gravity level is reduced,there is a new balance between inertial and interfacial forces, altering the behavior of the flow. In bubbly flow,the absence of drift velocity leads to spherical-shaped bubbles with a rectilinear trajectory.Slug flow is a succession of long bubbles and liquid slug carrying a few bubbles. There is no flow reversal in the thin liquid film as the long bubble and liquid slug pass over the film. Although the flow structure seems to be simpler than in normal gravity conditions,the models developed for the prediction of flow behavior in normal gravity and extended to reduced gravity flow are unable to predict the flow behavior correctly.An additional benefit of conducting studies in microgravity flows is that these studies aide the development of understanding for normal gravity flow behavior by removing the effects of buoyancy on the shape of the interface and density driven shear flows between the gas and the liquid phases. The proposal calls to study specifically the following: 1) The dynamics of isolated bubbles in microgravity liquid flows will be analyzed: Both the dynamics of spherical isolated bubbles and their dispersion by turbulence, their interaction with the pipe wall,the behavior of the bubbles in accelerated or decelerated flows,and the dynamics of isolated cylindrical bubbles, their deformation in accelerated/decelerated flows (in converging or diverging channels), and bubble/bubble interaction. Experiments will consist of the use of Particle Image Velocimetry (PIV) and Laser Doppler Velocimeters (LDV) to study single spherical bubble and single and two cylindrical bubble behavior with respect to their influence on the turbulence of the surrounding liquid and on the wall 2) The dynamics of bubbly and slug flow in microgravity will be analyzed especially for the role of the coalescence in the transition from bubbly to slug flow (effect of fluid properties and surfactant), to identify clusters that promote coalescence and transition the void fraction distribution in bubbly and slug flow,to measure the wall friction in bubbly flow. These experiments will consist of multiple bubbles type flows and will utilize hot wire and film anemometers to measure liquid velocity and wall shear stress respectively and double fiber optic probes to measure bubble size and velocity as a function of tube radius and axial location.
Experiments have been conducted to investigate the effect of fluid mixing on the depressurization of a large liquid hydrogen storage tank. The test tank is approximately ellipsoidal, having a volume of 4.89 m(exp 3) and an average wall heat flux of 4.2 W/m(exp 2) due to external heat input. A mixer unit was installed near the bottom of the tank to generate an upward directed axial jet flow normal to the liquid-vapor interface. Mixing tests were initiated after achieving thermally stratified conditions in the tank either by the introduction of hydrogen gas into the tank or by self-pressurization due to ambient heat leak through the tank wall. The subcooled liquid jet directed towards the liquid-vapor interface by the mixer induced vapor condensation and caused a reduction in tank pressure. Tests were conducted at two jet submergence depths for jet Reynolds numbers from 80,000 to 495,000 and Richardson numbers from 0.014 to 0.52. Results show that the rate of tank pressure change is controlled by the competing effects of subcooled jet flow and the free convection boundary layer flow due to external tank wall heating. It is shown that existing correlations for mixing time and vapor condensation rate based on small scale tanks may not be applicable to large scale liquid hydrogen systems.
Semiconductor photocatalyst approaches for solar CO 2 reduction are attractive due to their simplicity but have lagged in efficiency compared to less-integrated photoelectrochemical (PEC) approaches and to electrolysis reactors. We identify poor mass transport and catalyst deactivation as key constraints. To address them, we have developed a continuous-flow photocatalytic reactor system allowing us to control the triple-phase interface on the photocatalyst surface using the liquid and reactant gas flow rates. With the goal of selectively producing CO, the reactor is optimized by controlling the pressure and flow rates of the reactant gas and electrolyte in contact with both sides with the intermediately placed catalyst. In comparison to batch reactors with an immobile photocatalyst bed and gas phase CO 2 or CO 2 purged water, 10–24 times higher production rates are achieved for photocatalysts such as TiO 2 , ZnO, C 3 N 4 , and CdS by simply changing to the designed flow-type photoreactor without any catalyst modification. In addition, CO selectivity (93.2%) and long-term stability (>780 min) using the designed reactor are significantly enhanced compared to using the batch reactors (71.7%, <180 min for reduced 50% activity). Here, we propose that the enhanced mass transport on the photocatalyst surface accelerates the desorption of the initial photolysis product, CO, and prevents the poisoning effect from deactivating photocatalyst activity. This study has the potential to facilitate the utilization of semiconductor-based photocatalytic reactions for achieving superior performance wih gaseous reactants.
An analysis to predict the pressurant gas requirements for the discharge of cryogenic liquid propellants from storage tanks is presented, along with an algorithm and two computer programs. One program deals with the pressurization (ramp) phase of bringing the propellant tank up to its operating pressure. The method of analysis involves a numerical solution of the temperature and velocity functions for the tank ullage at a discrete set of points in time and space. The input requirements of the program are the initial ullage conditions, the initial temperature and pressure of the pressurant gas, and the time for the expulsion or the ramp. Computations are performed which determine the heat transfer between the ullage gas and the tank wall. Heat transfer to the liquid interface and to the hardware components may be included in the analysis. The program output includes predictions of mass of pressurant required, total energy transfer, and wall and ullage temperatures. The analysis, the algorithm, a complete description of input and output, and the FORTRAN 4 program listings are presented. Sample cases are included to illustrate use of the programs.
A lattice Boltzmann (LB) model with an efficient and accurate interface treatment for conjugate heat transfer across a thin wall between two different media is developed. The proposed interface treatment avoids fine meshing and computation within the thin layer; instead, the energy balance within the thin layer and the conjugate conditions on each interface are utilized to construct explicit updating schemes for the microscopic distribution functions of the LB model at the interior lattice nodes of the two media next to the thin layer. The proposed interface schemes reduce to the standard interface scheme for conjugate conditions in the literature in the limit of zero thickness of the thin layer, and thus it can be considered a more general interface treatment. A simplified version of the interface treatment is also proposed when the heat flux variation along the tangential direction of the thin layer is negligible. Three representative numerical tests are conducted to verify the applicability and accuracy of the proposed interface schemes. The results demonstrate that the intrinsic second-order accuracy of the LB model is preserved with the proposed interface schemes for thin layers with constant tangential fluxes, while for general situations with varying tangential fluxes, first-order accuracy is obtained. Finally, this interface treatment within the LB framework is attractive in conjugate heat transfer modeling involving thin layers for its simplicity, accuracy, and significant reduction in computational resources.
A steady shear flow is set up in a planar liquid layer bounded below by a rigid plate, and above by an interface with a passive gas, by imposing a temperature gradient along the layer and driving the motion by thermocapillarity. The thermal-convective instabilities to which this dynamic state is susceptible include stationary longitudinal rolls involving the classical Marangoni instability, and unsteady hydrothermal waves, whose new instability mechanism derives its energy from the horizontal temperature gradients. The thermal stability characteristics for liquid layers with and without return flow profiles are presently given as functions of the Prandtl number for the liquid and the Biot number for the interface, and comparisons are made with available experimental observations.
A two-phase CFD model for tank pressurization and following jet mixing of a cryogenic storage tank is presented using VOF approach for representing the phase boundary and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. The CFD model was validated against pressurization and liquid jet mixing data for a 110-inch diameter tank provided by Bullard1. Cases with like-gas and non-condensable gas pressurization are studied. The results of the cases with like-gas pressurization and following mixing are presented first, focusing on the effects of turbulence at the vapor-liquid interface on the tank pressure and phase change rates predictions. The second part of this paper is devoted to testing and validating the CFD model against non-condensable gas pressurization and following mixing. Tank pressures predicted in both cases are compared with each other and with the experimental data.