Search NASA⌕ Search

SEARCH · Search NASA

Results for “Interfacial flows”

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

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

At least 271 records · Page 15

Two-layer viscous instability in a rotating couette device

A novel experiment to study the interfacial shear instability between two liquids is described. Density-matched immiscible liquids are confined between concentric cylinders such that the interface is parallel to the cylinder walls. Interfacial waves that develop because of viscosity differences between the shearing fluids are studied as a function of rotation rate and depth ratio using optical techniques. Conditions neutral stability and the most unstable wavenumber agree reasonably well with predictions from linear stability analysis of the Navier-Stokes equations. Illumination using laser sheets allows precise measurement of the interface shape. Future experiments will verify the correctness of weakly nonlinear theories that describe energy transfer and saturation of wave growth by nonlinear effects. Measurements of solitary wave shapes, that occur far above neutral stability, will be compared to similar measurements for systems that have gravity as an important force to determine how gravity effects large disturbances. These results will be used to interpret slug and annular flow data that have been obtained in other mu g studies.

Gallagher, Christopher↗

Droplet breakup by multimodal nonlinear Rayleigh Taylor instability

A droplet impacted by a shock wave will undergo a process of fragmentation due to the development of interfacial hydrodynamic instabilities. The interface experiences variable acceleration and shear that result in the development of both inertial (Rayleigh–Taylor) and shear (Kelvin–Helmholtz) instabilities. These perturbations grow in time and drive the fragmentation and breakup of the deformed droplet. Experiments are performed on nominally 0.86 mm water droplet subjected to a Mach 7.6 detonation wave, resulting in a high Weber number (~36,000) breakup event. Perturbation growth is measured from a series of high-speed (> MHz) shadowgraph images. It is proposed that, given the size of the large-scale perturbations observed in experiments, these instabilities are growing in the non-linear regime and can be described by bubble-merger models for nonlinear mixing. Calculations are performed for the growth rates and size of these instabilities using deformation and external flow models to establish the time-dependent boundary conditions. The concurrence of the measured perturbation widths and the predictions of the simple model lend credence to the theory. This novel approach serves to open a new avenue in the characterization of droplet breakup via hydrodynamic instabilities.

Droplet breakup↗

Development of an Atomization Methodology for Spray Combustion

In liquid rocket propulsion, the knowledge and the understanding of liquid-gas interfacial phenomena are very important. This is important for predicting the onset of cavitation occurring in swirl injection elements used in STME, as well as atomization processes in shear-induced injectors (co-axial) and impinging injector elements. From the fact that all the physical processes including droplet size distribution, droplet dispersion, mixing and combustion are controlled by atomization processes, it is expected that the successful incorporation of the volume of fraction (VOF) will greatly enhance the analytical capability of predicting spray combustion processes in liquid-fueled engines. In this paper, a methodology is developed to define and track interfaces between two fluids in non-orthogonal, body-fitted grids using a single fractional volume of fluid (VOF) variable to describe the distribution of the liquid phase in a gas-liquid flow field. This method was implemented in a mature CFD code MAST (Multiphase All-Speed Transient) utilizing the general PISO-C algorithm. For the preliminary study on the analysis of spray combustion and tracking of the interface between two phases, we will report on the progress of the simulation of the instability on the liquid column; the surface wave instability and the droplet breakup from the liquid surface.

Seung, S. P.↗

Unusual Contact-Line Dynamics of Thick Films and Drops

We report several novel phenomena In contact-line and fingering dynamics of macroscopic spinning drops and gravity-driven films with dimensions larger than the capillary length. It is shown through experimental and theoretical analysis that such macroscopic films can exhibit various interfacial shapes, including multi valued ones, near the contact line due to a balance between the external body forces with capillarity. This rich variety of front shapes couples with the usual capillary, viscous, and intermolecular forces at the contact line to produce a rich and unexpected spectrum of contact-line dynamics. A single finger develops when part of the front becomes multivalued on a partially wetting macroscopic spinning drop in contrast to a different mechanism for microscopic drops of completely wetting fluids. Contrary to general expectation, we observe that, at high viscosity and low frequencies of rotation, the speed of a glycerine finger increases with increasing viscosity. Completely wetting Dow Corning 200 Fluid spreads faster over a dry inclined plane than a prewetted one. The presence of a thin prewetted film suppresses fingering both for gravity-driven flow and for spin coating. We analyze some of these unique phenomena in detail and offer qualitative physical explanations for the others.

Veretennikov, Igor↗

Wave effects on ocean-ice interaction in the marginal ice zone

The effects of wave train on ice-ocean interaction in the marginal ice zone are studied through numerical modeling. A coupled two-dimensional ice-ocean model has been developed to include wave effects and wind stress for the predictions of ice edge dynamics. The sea ice model is coupled to the reduced-gravity ocean model through interfacial stresses. The main dynamic balance in the ice momentum is between water-ice stress, wind stress, and wave radiation stresses. By considering the exchange of momentum between waves and ice pack through radiation stress for decaying waves, a parametric study of the effects of wave stress and wind stress on ice edge dynamics has been performed. The numerical results show significant effects from wave action. The ice edge is sharper, and ice edge meanders form in the marginal ice zone owing to forcing by wave action and refraction of swell system after a couple of days. Upwelling at the ice edge and eddy formation can be enhanced by the nonlinear effects of wave action; wave action sharpens the ice edge and can produce ice meandering, which enhances local Ekman pumping and pycnocline anomalies. The resulting ice concentration, pycnocline changes, and flow velocity field are shown to be consistent with previous observations.

Liu, Antony K.↗

Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

Flexible inhibitors are generally used in solid rocket motors (SRMs) as a means to control the burning of propellant. Vortices generated by the flow of propellant around the flexible inhibitors have been identified as a driving source of instabilities that can lead to thrust oscillations in launch vehicles. Potential coupling between the SRM thrust oscillations and structural vibration modes is an important risk factor in launch vehicle design. As a means to predict and better understand these phenomena, a multidisciplinary simulation capability that couples the NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This capability is crucial to the development of NASA's new space launch system (SLS). This paper summarizes the efforts in applying the coupled software to demonstrate and investigate fluid-structure interaction (FSI) phenomena between pressure waves and flexible inhibitors inside reusable solid rocket motors (RSRMs). The features of the fluid and structural solvers are described in detail, and the coupling methodology and interfacial continuity requirements are then presented in a general Eulerian-Lagrangian framework. The simulations presented herein utilize production level CFD with hybrid RANS/LES turbulence modeling and grid resolution in excess of 80 million cells. The fluid domain in the SRM is discretized using a general mixed polyhedral unstructured mesh, while full 3D shell elements are utilized in the structural domain for the flexible inhibitors. Verifications against analytical solutions for a structural model under a steady uniform pressure condition and under dynamic modal analysis show excellent agreement in terms of displacement distribution and eigenmode frequencies. The preliminary coupled results indicate that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor. This insight could have profound implications for SRM and flexible inhibitor designs for current and future launch vehicles including SLS.

Yang, H. Q.↗

Validation of a Two-Phase CFD Model for Predicting Tank Self-Pressurization in the Ground-Based K-Site Experiment

A two-phase CFD model for self-pressurization of a cryogenic storage tank partially filled with liquid hydrogen is presented using the Volume-Of-Fluid approach for representing the phase boundary and the associated interfacial heat, mass and momentum transfer between the liquid and the vapor regions. The CFD model is validated against self-pressurization experiment performed in the K-site flightweight hydrogen storage tank at NASA Glenn Research Center. Laminar and turbulent simulations together with conjugated heat transfer analysis are performed. Effects of turbulence, as well as tank wall conduction are presented and discussed. Predicted tank pressures and fluid temperatures are compared with the experimental data at two different heat loads for validating the CFD model.

phase change↗

A well-posed numerical method to track isolated conformal map singularities in Hele-Shaw flow

We present a new numerical method for calculating an evolving 2D Hele-Shaw interface when surface tension effects are neglected. In the case where the flow is directed from the less viscous fluid into the more viscous fluid, the motion of the interface is ill-posed; small deviations in the initial condition will produce significant changes in the ensuing motion. This situation is disastrous for numerical computation, as small round-off errors can quickly lead to large inaccuracies in the computed solution. Our method of computation is most easily formulated using a conformal map from the fluid domain into a unit disk. The method relies on analytically continuing the initial data and equations of motion into the region exterior to the disk, where the evolution problem becomes well-posed. The equations are then numerically solved in the extended domain. The presence of singularities in the conformal map outside of the disk introduces specific structures along the fluid interface. Our method can explicitly track the location of isolated pole and branch point singularities, allowing us to draw connections between the development of interfacial patterns and the motion of singularities as they approach the unit disk. In particular, we are able to relate physical features such as finger shape, side-branch formation, and competition between fingers to the nature and location of the singularities. The usefulness of this method in studying the formation of topological singularities (self-intersections of the interface) is also pointed out.

Baker, Gregory↗

Influence of Nanoconfinement and Elevated Temperatures on Geocolloid-Facilitated Transport of Energy-related Contaminants

The main objective of this project is to obtain a better understanding of the transport behavior and interactions of geocolloids in the presence of energy- related contaminants under bulk and nanoconfined conditions. In that, we first established solid protocols to synthesize and fabricate numerous types of geocolloids with a. various sizes ranging from 40 – 800 nm and b. different hydrophile-lipophile balance (HLB) ratios in the range of 25:75 to 75:25. Direct force measurements with geocolloids having different degrees of surface coverages (i.e. HLB to realistically mimic adsorption of energy related contaminants) were conducted using the Surface Forces Apparatus (SFA) over a distance regime starting from 8μm all the way down to molecular contact. Repulsive forces were observed on approach starting from > 3μm, followed by an exponential increase of which magnitude appears to be larger than a decay length obtained from Derjaguin–Landau–Verwey–Overbeek (DLVO) theory in pure water. When the geocolloids were confined in salted water, the magnitude of onset of repulsion was varied as a function of salinity in solution, which can significantly alter a purely repulsive screened electrostatic (coulombic) interaction arising from, among geocolloids as well as between geocolloids and geosurfaces. The viscosity and flow characteristics of geocolloidal suspensions at different degrees of confinement were also investigated where we identified highly discontinuous rheological behaviors below a critical nanoconfinement level. We anticipate that the knowledge gained through this study will enable the scientists and researchers to better assess transport and fate behaviors of geocolloidal dispersions that can carry energy-related contaminants under realistically emulated geosystem.

03 NATURAL GAS↗

Flame Synthesis of Single- and Multi-Walled Carbon Nanotubes and Nanofibers

Metal-catalyzed carbon nanotubes are highly sought for a diverse range of applications that include nanoelectronics, battery electrode material, catalysis, hydrogen storage media and reinforcing agents in polymer composites. These latter applications will require vast quantities of nanotubes at competitive prices to be economically feasible. Moreover, reinforcing applications may not require ultrahigh purity nanotubes. Indeed, functionalization of nanotubes to facilitate interfacial bonding within composites will naturally introduce defects into the tube walls, lessening their tensile strength. Current methods of aerosol synthesis of carbon nanotubes include laser ablation of composite targets of carbon and catalyst metal within high temperature furnaces and decomposition of a organometallics in hydrocarbons mixtures within a tube furnace. Common to each approach is the generation of particles in the presence of the reactive hydrocarbon species at elevated temperatures. In the laser-ablation approach, the situation is even more dynamic in that particles and nanotubes are borne during the transient cooling phase of the laser-induced plasma for which the temperature far exceeds that of the surrounding hot gases within the furnace process tube. A shared limitation is that more efficient methods of nanoparticle synthesis are not readily incorporated into these approaches. In contrast, combustion can quite naturally create nanomaterials such as carbon black. Flame synthesis is well known for its commercial scalability and energy efficiency. However, flames do present a complex chemical environment with steep gradients in temperature and species concentrations. Moreover, reaction times are limited within buoyant driven flows to tens of milliseconds. Therein microgravity can greatly lessen temperature and spatial gradients while allowing independent control of flame residence times. In preparation for defining the microgravity experiments, the work presented here focuses on the effect of catalyst particle size and reactant gas in 1g.

VanderWal, R. L.↗

Computer Code for Nanostructure Simulation

Due to their small size, nanostructures can have stress and thermal gradients that are larger than any macroscopic analogue. These gradients can lead to specific regions that are susceptible to failure via processes such as plastic deformation by dislocation emission, chemical debonding, and interfacial alloying. A program has been developed that rigorously simulates and predicts optoelectronic properties of nanostructures of virtually any geometrical complexity and material composition. It can be used in simulations of energy level structure, wave functions, density of states of spatially configured phonon-coupled electrons, excitons in quantum dots, quantum rings, quantum ring complexes, and more. The code can be used to calculate stress distributions and thermal transport properties for a variety of nanostructures and interfaces, transport and scattering at nanoscale interfaces and surfaces under various stress states, and alloy compositional gradients. The code allows users to perform modeling of charge transport processes through quantum-dot (QD) arrays as functions of inter-dot distance, array order versus disorder, QD orientation, shape, size, and chemical composition for applications in photovoltaics and physical properties of QD-based biochemical sensors. The code can be used to study the hot exciton formation/relation dynamics in arrays of QDs of different shapes and sizes at different temperatures. It also can be used to understand the relation among the deposition parameters and inherent stresses, strain deformation, heat flow, and failure of nanostructures.

Filikhin, Igor↗

Numerical Modeling of Crystal of ZnSe by Physical Vapor Transport - Towards a more Comprehensive Formulations

Crystal growth from the vapor phase has various advantages over melt growth. The main advantage is from a lower processing temperature which makes the process more amenable in instances where the melting temperature of the crystal is high. Other benefits stem from the inherent purification mechanism in the process due to differences in the vapor pressures of the native elements and impurities, and the enhanced interfacial morphological stability during the growth process. Further, the implementation of PVT growth in closed ampoules affords experimental simplicity with minimal needs for complex process control which makes it an ideal candidate for space investigations in systems where gravity tends to have undesirable effects on the growth process. Bulk growth of wide band gap II-VI semiconductors by physical vapor transport has been developed and refined over the past several years at NASA MSFC. Results from a modeling study of PVT crystal growth of ZnSe are reported in this paper. The PVT process is numerically investigated using both two-dimensional and fully three-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and the compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles will be provided in the final manuscript along with computed heat and mass transfer rates. Comparisons with the 1-D model will also be provided. The effect of gravity on the process from numerical computations shows subtle effects although experimental evidence from vertically and horizontally grown samples show dramatic evidence of gravitational effects. The shortcomings of the problem formulation will be discussed and a framework will be provided leading up towards a more comprehensive model of PVT systems.

Ramachandran, N.↗

Crystal Growth of ZnSe by Physical Vapor Transport: A Modeling Study

Crystal growth from the vapor phase has various advantages over melt growth. The main advantage is from a lower processing temperature which makes the process more amenable in instances where the melting temperature of the crystal is high. Other benefits stem from the inherent purification mechanism in the process due to differences in the vapor pressures of the native elements and impurities, and the enhanced interfacial morphological stability during the growth process. Further, the implementation of Physical Vapor Transport (PVT) growth in closed ampoules affords experimental simplicity with minimal needs for complex process control which makes it an ideal candidate for space investigations in systems where gravity tends to have undesirable effects on the growth process. Bulk growth of wide band gap II-VI semiconductors by physical vapor transport has been developed and refined over the past several years at NASA MSFC. Results from a modeling study of PVT crystal growth of ZnSe arc reported in this paper. The PVI process is numerically investigated using both two-dimensional and fully three-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and the compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The preliminary results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles will be provided in the final manuscript along with computed heat and mass transfer rates. Comparisons with the 1-D model will also be provided.

Ramachandran, Narayanan↗

CFD Modeling of Tank Pressurization and Axial Jet Mixing Experiments with and without Non-Condensable Gas in the Ullage

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.

Computational Fluid Dynamics↗

CFD Modeling of Tank Pressurization and Axial Jet Mixing Experiments with and without Non-Condensable Gas in the Ullage

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 Bullard. 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.

Computational Fluid Dynamics↗

Scaling microstructural processes in the sintering of ionic ceramics

A multi-scale framework, combining a multiphase field formulation and large deformation mechanics, was developed as a stepping stone to perform the data analytics of the microstructural level kinetics of a sintering solid. Relevant microstructural information from this framework, such as grain, stress, and porosity statistics, was scaled up to describe the macroscopic level sintering kinetics. Here, the developed formulation was applied to describe the electric field assisted sintering of Y 2 O 3 . Microstructural inhomogeneities in a multi-granular solid result in the formation of a field of compressive stress networks, which interleave with low compression and weakly tensile regions, defining a scaffolding for sintering concentration regions to develop. A Poisson effect-induced lateral stress network is also naturally self-induced as a result of the mechanical constraints imposed by the sintering apparatus. For long sintering times, localized shear stresses enhancing mass flow along grain boundaries and internal surfaces develop. Three-sided pores are removed by either vacancy transport to the surrounding pores, or move towards the external surfaces through grain boundary diffusion. Four- and higher order-sided pores stabilize because an equal amount of vacancies are gained and lost through the connecting grain boundaries. Grain dewetting contributes to pore coalescence, suggesting that pore kinetics and grain growth are coupled and should be analyzed in concert. The combined sintering and grain growth kinetics define six regimes of sintering behavior: (1) T, the transient regime; (2) E$_Υ$, the surface energy dominated, early sintering regime, where the grain growth exponent, p = 1, and the stress concentration factor, $f$ ~ $1/\hat{ρ}^{4.6}$; (3) E S , the stress dominated, early sintering regime, where p = 1 and $f$ ~ $1/\hat{ρ}^{4}$; (4) I$_Υ$, the surface energy dominated, intermediate sintering regime, where p = 2 and $f$ ~ $1/\hat{ρ}^{4.6}$; (5) I S , the stress dominated, intermediate sintering regime, where p = 2 and $f$ ~ $1/\hat{ρ}^{4}$; and (6) L, the late sintering regime, where p = 3 and $f$ ~ 1. At the macroscopic level, the rapid densification and suppression of grain growth observed in the electric field assisted sintering process is a consequence of the compounding effects of the underlying stress-, transport-, and interfacial-energy-induced energy minimization kinetics, as predicted by the multi-scale framework.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Bonding by Hydroxide-Catalyzed Hydration and Dehydration

A simple, inexpensive method for bonding solid objects exploits hydroxide-catalyzed hydration and dehydration to form silicate-like networks in thin surface and interfacial layers between the objects. The method can be practiced at room temperature or over a wide range of temperatures. The method was developed especially to enable the formation of precise, reliable bonds between precise optical components. The bonds thus formed exhibit the precision and transparency of bonds formed by the conventional optical-contact method and the strength and reliability of high-temperature frit bonds. The method also lends itself to numerous non-optical applications in which there are requirements for precise bonds and/or requirements for bonds, whether precise or imprecise, that can reliably withstand severe environmental conditions. Categories of such non-optical applications include forming composite materials, coating substrates, forming laminate structures, and preparing objects of defined geometry and composition. The method is applicable to materials that either (1) can form silicate-like networks in the sense that they have silicate-like molecular structures that are extensible into silicate-like networks or (2) can be chemically linked to silicate-like networks by means of hydroxide-catalyzed hydration and dehydration. When hydrated, a material of either type features surface hydroxyl (-OH) groups. In this method, a silicate-like network that bonds two substrates can be formed either by a bonding material alone or by the bonding material together with material from either or both of the substrates. Typically, an aqueous hydroxide bonding solution is dispensed and allowed to flow between the mating surfaces by capillary action. If the surface figures of the substrates do not match precisely, bonding could be improved by including a filling material in the bonding solution. Preferably, the filling material should include at least one ingredient that can be hydrated to have exposed hydroxyl groups and that can be chemically linked, by hydroxide catalysis, to a silicate-like network. The silicate-like network could be generated in situ from the filling material and/or substrate material, or could be originally present in the bonding material.

Gwo, Dz-Hung↗

Hydrodynamic characterization of the redox chemistry of crown-encapsulated uranyl complexes

The redox properties of actinide-containing species strongly influence their reactivity, speciation, and interfacial behavior, but the experimental quantification of the electrochemical characteristics of molecular actinide complexes in nonaqueous media has not received the attention it deserves. Here, results from hydrodynamic methods and electrochemical simulations of U(VI)/U(V) redox are reported, including quantification of heterogeneous electron-transfer kinetics and estimation of chemical reversibility of U(VI)/U(V) interconversion at electrodes in acetonitrile-based electrolyte. The complexes investigated are recently reported U(VI) and U(V) complexes in which the uranyl ion (UO 2 n+ ) is encapsulated in a macrocyclic 18-crown-6-like moiety templated by a Pt(II) center. These complexes feature the most positive value U VI /U V reduction potential yet reported and are thus particularly relevant to study of facile U(V) generation from U(VI) precursors as well as uranium electroanalysis. Rotating disk electrode (RDE) studies have been used to quantify the diffusion coefficients of the U(VI) and U(V) complexes, and standard heterogeneous electron transfer rate constants ( k 0 ) for the redox have been determined using a conventional Koutecký-Levich analysis. Rotating ring-disk electrode (RRDE) studies have been used to directly interrogate the chemical reversibility of U(VI)-U(V) interconversion, confirming that reduction of the U(VI) complex at an Au disk is associated with formation of the U(V) analogue that can be readily re-oxidized at a Pt ring under hydrodynamic (rotating) conditions. Because measurements of the type reported here are generally associated with current flows that are larger than those found in corresponding quiescent (unstirred) conditions, our findings suggest that hydrodynamic methods could be advantageous for design of electroanalytical approaches to detection of actinide species and study of their redox properties.

actinides↗