Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two Phase Flow”

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 1,675 records · Page 93

The effect of gravity induced free convection upon the melting phenomena of a finite paraffin slab for thermal control

Two theoretical models were developed to predict the thermal response of the phase change material to a given hot plate temperature. A two-dimensional pure conduction model was developed to predict the melting of the phase change material when heat transfer was a function of conduction. A combined conduction-convection model, also two-dimensional, was developed to predict the phase change phenomena when heat transfer was a function of conduction and gravity-induced free convection. Both models were solved using explicit finite difference approximations on a digital computer. The experimental equipment consisted of a rectangular cell utilizing a heat chamber, an expansion chamber, and a test chamber; a sixteen channel multipoint recorder, and a fluid flow system. The recorder monitered hot and cold plate temperatures and interior node temperatures at two second intervals. A comparison of theoretical temperature profiles and experimental temperature profiles is presented for six runs at various angles of inclination of the test cell with respect to the horizontal direction. A detailed discussion of results is presented.

Bain, R. L.↗

Development of PFLOTRAN Transport Capability for Use in the Waste Isolation Pilot Plant Performance Assessment - 20545

Waste Isolation Pilot Plant (WIPP) performance assessment (PA) calculations estimate the probability of radionuclide release from the repository to the land surface and across the land withdrawal boundary for a regulatory period of 10,000 years after facility closure. Simulations of flow and transport in the repository and the surrounding Salado Formation are foundational to the PA. Because proposed additional waste emplacement panels would result in an asymmetric repository layout, the US Department of Energy (DOE) is preparing to transition to use of a three-dimensional (3-D) model domain for simulation of flow and transport instead of the two-dimensional (2-D) flared grid domain currently used. DOE has charged Sandia National Laboratories with developing the capability necessary to simulate processes affecting flow and transport in the WIPP in PFLOTRAN, an open-source massively parallel multi-phase flow and reactive transport code. The new flow and transport capabilities developed in PFLOTRAN incorporate WIPP-specific process models and will replace the 2-D simulators (BRAGFLO and NUTS) that are currently utilized for Salado flow and transport calculations in WIPP PA. The focus of this paper is on the development of a new Nuclear Waste Transport (NWT) mode in PFLOTRAN that has all of the capabilities necessary for Salado transport simulations, including the ability to handle complete dry-out (100% gas saturation) of arbitrary cells in the model domain, radionuclide mass conservation at step changes in porosity associated with borehole intrusion, and the ability to calculate fluxes on a flared grid. The new PFLOTRAN transport capability and a suite of verification tests were designed around a list of functional requirements for WIPP PA calculations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Turbine blade and vane heat flux sensor development, phase 2

The development of heat flux sensors for gas turbine blades and vanes and the demonstration of heat transfer measurement methods are reported. The performance of the heat flux sensors was evaluated in a cylinder in cross flow experiment and compared with two other heat flux measurement methods, the slug calorimeter and a dynamic method based on fluctuating gas and surface temperature. Two cylinders, each instrumented with an embedded thermocouple sensor, a Gardon gauge, and a slug calorimeter, were fabricated. Each sensor type was calibrated using a quartz lamp bank facility. The instrumented cylinders were then tested in an atmospheric pressure combustor rig at conditions up to gas stream temperatures of 1700K and velocities to Mach 0.74. The test data are compared to other measurements and analytical prediction.

Atkinson, W. H.↗

Compressible 2D Flow Field Interaction of Two Contra-Rotating Blades

The blades of coaxial, contra-rotating rotor systems cross each other in close proximity and at high relative speeds. This crossing event is a potential source of noise and severe blade loads. Effects of compressibility can aggravate the interaction and significantly alter the pressure field signature and phase relationships. A 2-D analysis of this phenomenon is performed by simulating two airfoils passing each other at specified speeds and vertical separation distances. Several test cases spanning a relevant range of Reynolds numbers, angles of attack, and relative Mach number are considered. The Mach number is varied to simulate the radial variation of velocity from the root to tip of a rotor blade to capture the pressure signature, lift, and drag of the airfoils. The velocity and pressure distributions on the airfoils, and in the space between the airfoils are computed before, at, and after airfoil crossing. The variations of lift and drag coefficients through the interaction are captured. The upper airfoil experiences an increase in lift followed by a very sharp drop in lift during the interaction. When relative Mach numbers are transonic, the region of interaction is greatly extended, with shock interactions occurring. The results show the complex nature of the aerodynamic and fluid dynamic impulses generated by blade-blade interactions, with implications to aeroelastic loads and aeroacoustic sources.

Barbely, Natasha L.↗

Catalytic Microtube Rocket Igniter

Devices that generate both high energy and high temperature are required to ignite reliably the propellant mixtures in combustion chambers like those present in rockets and other combustion systems. This catalytic microtube rocket igniter generates these conditions with a small, catalysis-based torch. While traditional spark plug systems can require anywhere from 50 W to multiple kW of power in different applications, this system has demonstrated ignition at less than 25 W. Reactants are fed to the igniter from the same tanks that feed the reactants to the rest of the rocket or combustion system. While this specific igniter was originally designed for liquid methane and liquid oxygen rockets, it can be easily operated with gaseous propellants or modified for hydrogen use in commercial combustion devices. For the present cryogenic propellant rocket case, the main propellant tanks liquid oxygen and liquid methane, respectively are regulated and split into different systems for the individual stages of the rocket and igniter. As the catalyst requires a gas phase for reaction, either the stored boil-off of the tanks can be used directly or one stream each of fuel and oxidizer can go through a heat exchanger/vaporizer that turns the liquid propellants into a gaseous form. For commercial applications, where the reactants are stored as gases, the system is simplified. The resulting gas-phase streams of fuel and oxidizer are then further divided for the individual components of the igniter. One stream each of the fuel and oxidizer is introduced to a mixing bottle/apparatus where they are mixed to a fuel-rich composition with an O/F mass-based mixture ratio of under 1.0. This premixed flow then feeds into the catalytic microtube device. The total flow is on the order of 0.01 g/s. The microtube device is composed of a pair of sub-millimeter diameter platinum tubes connected only at the outlet so that the two outlet flows are parallel to each other. The tubes are each approximately 10 cm long and are heated via direct electric resistive heating. This heating brings the gasses to their minimum required ignition temperature, which is lower than the auto-thermal ignition temperature, and causes the onset of both surface and gas phase ignition producing hot temperatures and a highly reacting flame. The combustion products from the catalytic tubes, which are below the melting point of platinum, are injected into the center of another combustion stage, called the primary augmenter. The reactants for this combustion stage come from the same source but the flows of non-premixed methane and oxygen gas are split off to a secondary mixing apparatus and can be mixed in a near-stoichiometric to highly lean mixture ratio. The primary augmenter is a component that has channels venting this mixed gas to impinge on each other in the center of the augmenter, perpendicular to the flow from the catalyst. The total crosssectional area of these channels is on a similar order as that of the catalyst. The augmenter has internal channels that act as a manifold to distribute equally the gas to the inward-venting channels. This stage creates a stable flame kernel as its flows, which are on the order of 0.01 g/s, are ignited by the combustion products of the catalyst. This stage is designed to produce combustion products in the flame kernel that exceed the autothermal ignition temperature of oxygen and methane.

Schneider, Steven J.↗

Effects of upstream-biased third-order space correction terms on multidimensional Crowley advection schemes

The impact of upstream-biased corrections for third-order spatial truncation error on the stability and phase error of the two-dimensional Crowley combined advective scheme with the cross-space term included is analyzed, putting primary emphasis on phase error reduction. The various versions of the Crowley scheme are formally defined, and their stability and phase error characteristics are intercompared using a linear Fourier component analysis patterned after Fromm (1968, 1969). The performances of the schemes under prototype simulation conditions are tested using time-dependent numerical experiments which advect an initially cone-shaped passive scalar distribution in each of three steady nondivergent flows. One such flow is solid rotation, while the other two are diagonal uniform flow and a strongly deformational vortex.

Schlesinger, R. E.↗

Simulation of mineral dissolution at the pore scale with evolving fluid-solid interfaces: review of approaches and benchmark problem set

This manuscript presents a benchmark problem for the simulation of single-phase flow, reactive transport, and solid geometry evolution at the pore scale. The problem is organized in three parts that focus on specific aspects: flow and reactive transport (part I), dissolution-driven geometry evolution in two dimensions (part II), and an experimental validation of three-dimensional dissolution-driven geometry evolution (part III). Five codes are used to obtain the solution to this benchmark problem, including Chombo-Crunch, OpenFOAM-DBS, a lattice Boltzman code, Vortex, and dissolFoam. These codes cover a good portion of the wide range of approaches typically employed for solving pore-scale problems in the literature, including discretization methods, characterization of the fluid-solid interfaces, and methods to move these interfaces as a result of fluid-solid reactions. A short review of these approaches is given in relation to selected published studies. Results from the simulations performed by the five codes show remarkable agreement both quantitatively--based on upscaled parameters such as surface area, solid volume, and effective reaction rate--and qualitatively--based on comparisons of shape evolution. This outcome is especially notable given the disparity of approaches used by the codes. Therefore, these results establish a strong benchmark for the validation and testing of pore-scale codes developed for the simulation of flow and reactive transport with evolving geometries. They also underscore the significant advances seen in the last decade in tools and approaches for simulating this type of problem.

58 GEOSCIENCES↗

Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport

In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transport barriers to self-propelled particles in fluid flows

Here, we present theory and experiments demonstrating the existence of invariant manifolds that impede the motion of microswimmers in two-dimensional fluid flows. One-way barriers are apparent in a hyperbolic fluid flow that block the swimming of both smooth-swimming and run-and-tumble Bacillus subtilis bacteria. We identify key phase-space structures, called swimming invariant manifolds (SwIMs), that serve as separatrices between different regions of long-time swimmer behavior. When projected into xy space, the edges of the SwIMs act as one-way barriers, consistent with the experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phase separation, edge currents, and Hall effect for active matter with Magnus dynamics

Here we examine run-and-tumble disks in two-dimensional systems where the particles also have a Magnus component to their dynamics. For increased activity, we find that the system forms a motility-induced phase-separated (MIPS) state with chiral edge flow around the clusters, where the direction of the current is correlated with the sign of the Magnus term. The stability of the MIPS state is non-monotonic as a function of increasing Magnus term amplitude, with the MIPS region first extending down to lower activities followed by a break up of MIPS at large Magnus amplitudes into a gel-like state. We examine the dynamics in the presence of quenched disorder and a uniform drive and find that the bulk flow exhibits a drive-dependent Hall angle. This is a result of the side jump effect produced by scattering from the pinning sites and is similar to the behavior found for skyrmions in chiral magnets with quenched disorder.

97 MATHEMATICS AND COMPUTING↗

Surface-Tension Driven Convection

The present experiments are directed to examine surface-tension gradients. The surface-tension induced convection under reduced gravity it's comparison with normal gravity convection with all other identical conditions are studied. The intensities of such flows and their penetration depth below the free surface are also determined. In this way, the dominant role of a free surface in a reduced-gravity environment is indicated. The interface between two fluid phases can influence the motion of fluids when either the interface has finite curvature or when the interfacial tension varies from point to point. In both cases, forces appear in the interfacial region that can affect or generate fluid motions.

Ostrach, S.↗

Thermal management of high power space based systems

Conventional techniques of using a portion of the spacecraft skin for radiation of waste heat will be inadequate for high powered payloads (50 to 100 kWe) due to the lack of sufficient area. A Shuttle type system using a pumped single phase fluid loop could be scaled up to higher power but this type of system would require excessive pump power and weight. A pumped two-phase heat transfer loop has a much lower pumping requirement due to the higher latent heat of vaporization of the fluid in comparison to the sensible heat it can absorb through a temperature change. Concepts for an evaporator and a condenser for a pumped two-phase system are described. The condenser uses capillary grooves and a separate pumped condensate return line to achieve high heat transfer coefficients and stable operation due to the separation of the vapor and liquid flows. The cold plate evaporator uses wicks to contain the liquid and transport it to the heated surface. It can also function as a condenser for warming components. Control concepts for the cold plate are discussed. Concepts for deployment or erection of large space radiators are also considered.

Hwangbo, H.↗

Ignition delay of a gas mixture above a liquid fuel pool

A computational study has been made of transient heat transfer and fluid flow of an axisymmetric two-layer gas-and-liquid system heated from above by a small hot spot (e.g., a hot wire or pilot flame) located close to the liquid surface. The gas phase is unconfined above the liquid pool. The effects of varying gravity level (0.0001 to 1 gn), liquid pool height, and heater height are investigated. Thermocapillary convection induced by the nonuniform heating of the liquid surface combines with buoyancy forces to affect the heat transfer and the transport of fuel vapor toward the heat source. At reduced gravity, gas-phase conduction is comparable to the incident radiation at the liquid surface, whereas at 1 gn, buoyant convection carries the heat upward from the heat source and, therefore, the liquid is heated primarily by radiation.

Schiller, D. N.↗

Observation of anomalous slow-mode shock and reconnection layer in the dayside magnetospause

Plasma and magnetic field data from the International Sun-Earth Explorer ISEE 2 spacecraft recorded during an outbound crossing of the dayside, northern hemisphere magnetopause in October 29, 1979, provide evidence for a slow shock (SS) in the observed reconnection layer. This layer is found to be bounded on the magnetosheath side by the SS; near the magnetospheric side of the layer, a second current sheet is found that may have been rotational discontinuity (RD). The direction of the accelerated plasma flow, the earthward sense of the normal magnetic filed across the SS and RD, and the relative orientation of the SS and the RD all indicate that the reconection site was located south of the spacecraft. Quantitative tests show that, allowing for experiemental uncertanties, data taken upstream and downstream of the SS are consistent with coplanarity and other Rankine-Hugoniot (RH) conditions. Examiniations of the flow parameters indicates two anomalous properties of the SS: the upstream flow, viewed in the deHoffmann-Teller frame, is superalfvenic and the downstream plasma is firehose unstable. In comparison to the long-wavelength slow-mode phase speed, however, the flow in the upstream region is super slow, while in the downstream region it is subslow, as requeired for a slow-mode shock. Further properties of the shock include a large decrease in total enthalpy across it, indicating the escape of a sizable heat flux from the shock structure, and the occurrence of a polarization reversal of the tangential magnetic field within the shock layer, a feature that is predicted by linear double-polytropic Hall-MHD and results from a large increase in pressure anistropy from the upstream region, where p(sub parallel) approximately equal to p(sub perp), to the downstream region, where p(sub parrallel) greater than p(sub perp). Quantitative tests of the RD-like discontinuity show that it satisfies the necessary RH conditions within experimental uncertainties and indicate that the flow across it is earthward, that is, it is the same as for the SS. Thus the RD is propagating behind the SS in the direction away from the Earth. The reversal in the order of the discontinuities is attributed to the reversal of the slow and intermediate phase speeds produced by the pressure anistropy in the region between the RD and the SS. The linear dispersion relation and polarization properties of dispersive MHD waves in an anisotropic plasma are examined in the appendix.

Walthour, D. W↗

Effect of Directional Array Size on the Measurement of Airframe Noise Components

A study was conducted to examine the effects of overall size of directional (or phased) arrays on the measurement of aeroacoustic components. An airframe model was mounted in the potential core of an open-jet windtunnel, with the directional arrays located outside the flow in an anechoic environment. Two array systems were used; one with a solid measurement angle that encompasses 31.6 deg.of source directivity and a smaller one that encompasses 7.2 deg. The arrays, and sub-arrays of various sizes, measured noise from a calibrator source and flap edge model setups. In these cases, noise was emitted from relatively small, but finite size source regions, with intense levels compared to other sources. Although the larger arrays revealed much more source region detail, the measured source levels were substantially reduced due to finer resolution compared to that of the smaller arrays. To better understand the measurements quantitatively, an analytical model was used to define the basic relationships between array to source region sizes and measured output level. Also, the effect of noise scattering by shear layer turbulence was examined using the present data and those of previous studies. Taken together, the two effects were sufficient to explain spectral level differences between arrays of different sizes. An important result of this study is that total (integrated) noise source levels are retrievable and the levels are independent of the array size as long as certain experimental and processing criteria are met. The criteria for both open and closed tunnels are discussed. The success of special purpose diagonal-removal processing in obtaining integrated results is apparently dependent in part on source distribution. Also discussed is the fact that extended sources are subject to substantial measurement error, especially for large arrays.

Brooks, Thomas F.↗

Evaporating Spray in Supersonic Streams Including Turbulence Effects

Evaporating spray plays an important role in spray combustion processes. This paper describes the development of a new finite-conductivity evaporation model, based on the two-temperature film theory, for two-phase numerical simulation using Eulerian-Lagrangian method. The model is a natural extension of the T-blob/T-TAB atomization/spray model which supplies the turbulence characteristics for estimating effective thermal diffusivity within the droplet phase. Both one-way and two-way coupled calculations were performed to investigate the performance of this model. Validation results indicate the superiority of the finite-conductivity model in low speed parallel flow evaporating sprays. High speed cross flow spray results indicate the effectiveness of the T-blob/T-TAB model and point to the needed improvements in high speed evaporating spray modeling.

Balasubramanyam, M. S.↗

Observations of Traveling Crossflow Resonant Triad Interactions on a Swept Wing

Experimental evidence indicates the presence of a triad resonance interaction between traveling crossflow modes in a swept wing flow. Results indicate that this interaction occurs when the stationary and traveling crossflow modes have similar and relatively low amplitudes (approx.1% to 6% of the total freestream velocity). The resonant interaction occurs at instability amplitudes well below those typically known to cause transition, yet transition is observed to occur just downstream of the resonance. In each case, two primary linearly unstable traveling crossflow modes are nonlinearly coupled to a higher frequency linearly stable mode at the sum of their frequencies. The higher-frequency mode is linearly stable and presumed to exist as a consequence of the interaction of the two primary modes. Autoand cross-bicoherence are used to determine the extent of phase-matching between the modes, and wavenumber matching confirms the triad resonant nature of the interaction. The bicoherence results indicate a spectral broadening mechanism and the potential path to early transition. The implications for laminar flow control in swept wing flows are significant. Even if stationary crossflow modes remain subcritical, traveling crossflow interactions can lead to early transition.

Eppink, Jenna L.↗

Fluid flow and heat convection studies for actively cooled airframes

This report details progress made on the jet impingement - liquid crystal - digital imaging experiment. With the design phase complete, the experiment is currently in the construction phase. In order to reach this phase two design related issues were resolved. The first issue was to determine NASP leading edge active cooling design parameters. Meetings were arranged with personnel at SAIC International, Torrance, CA in order to obtain recent publications that characterized expected leading edge heat fluxes as well as other details of NASP operating conditions. The information in these publications was used to estimate minimum and maximum jet Reynolds numbers needed to accomplish the required leading edge cooling, and to determine the parameters of the experiment. The details of this analysis are shown in Appendix A. One of the concerns for the NASP design is that of thermal stress due to large surface temperature gradients. Using a series of circular jets to cool the leading edge will cause a non-uniform temperature distribution and potentially large thermal stresses. Therefore it was decided to explore the feasibility of using a slot jet to cool the leading edge. The literature contains many investigations into circular jet heat transfer but few investigations of slot jet heat transfer. The first experiments will be done on circular jets impinging on a fiat plate and results compared to previously published data to establish the accuracy of the method. Subsequent experiments will be slot jets impinging on full scale models of the NASP leading edge. Table 1 shows the range of parameters to be explored. Next a preliminary design of the experiment was done. Previous papers which used a similar experimental technique were studied and elements of those experiments adapted to the jet impingement study. Trade-off studies were conducted to determine which design was the least expensive, easy to construct, and easy to use. Once the final design was settled, vendors were contacted to verify that equipment could be obtained to meet our specifications. Much of the equipment required to complete the construction of the experiment has been ordered or received. The material status list is shown in Appendix B.

Mills, A. F.↗