Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase 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,081 records · Page 60

Surface Structure of Reacting Solid Ramjet Fuels

Solid fuel ramjets are a developing technology that will enable extended range for supersonic vehicles. The simplicity of the system in terms of moving parts results in the performance being dictated by the solid fuel combustion processes. The physics near the surface, particularly in the fuel-rich recirculation zone, are challenging to resolve due to flow opacity. Here, regression rate measurements of hydroxyl-terminated polybutadiene from a miniature, portable slab burner are reported. Mass-averaged regression rates agree with correlations from the literature, suggesting that the proper physics in the burner are captured. Synchrotron-based phase-contrast imaging is used to quantify local surface features, including a multiphase froth layer, molten polymer transport, and local regression rate measurements in the recirculation zone. Regression rates near the inlet are less than the mass-averaged values, whereas near the reattachment point the local value exceeds the mass-averaged values. Further, the observed froth layer is as large as 2 mm near the inlet and decreases to sub-millimeter thicknesses near the reattachment point. An analytical heat transfer model is used to quantify the influence of the surface froth layer on the regression rates. Results indicate that the froth layer has an insulating effect and will reduce the local regression rate.

33 ADVANCED PROPULSION SYSTEMS↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

Final report of “A Detailed Study of the Physical Mechanisms Controlling CO2-Brine Capillary Trapping in the Subsurface” (University of Arizona, DE-SC0006696)

Carbon capture and storage (CCS) of carbon dioxide emissions generated by production or combustion of fossil fuels is a technologically viable means to reduce the build-up of CO2 in the atmosphere and oceans. Using advantages of scale and location, CCS is particularly suitable for large point sources near ubiquitous deep saline aquifers, depleted gas reservoirs, or at production reservoirs for enhanced oil recovery (EOR). In the BES-funded research project, Oregon State University (OSU) carried out capillary trapping experiments with proxy fluids that mimic the properties of the scCO2/brine system under ambient temperatures and pressures, and successfully developed a unique and novel x-ray compatible, high-pressure, elevated temperature setup to study the scCO2/brine system under challenging reservoir conditions. Both methodologies were applied to a variety of porous media, including synthetic (glass bead) and geologic (Bentheimer sandstone) materials. The University of Arizona (UA) developed pore-scale lattice Boltzmann (LB) models which are able to handle the experimental conditions for proxy fluids, as well as the scCO2/brine system, that are capable of simulating permeability in volumes of tens of millions of fluid elements. We reached the following summary findings (main institute indicated): 1. (OSU/UA) To understand capillary trapping in a multiphase fluid-porous medium system, the system must be analyzed from a pore-scale force balance perspective; trapping can be enhanced by manipulating wetting and nonwetting phase fluid properties. 2. (OSU) Pore-scale fluid connectivity and topology has a clear and direct effect on nonwetting phase capillary trapping efficiency. 3. (OSU) Rock type and flow regime also have a pronounced effects on capillary trapping. 4. (OSU/UA) There is a predictable relationship between NWP connectivity and NWP saturation, which allows for development of injection strategies that optimize trapping. The commonly used Land model (Land, 1968) does not predict amount of trapped NWP accurately. 5. (UA) There are ambiguities regarding the segmentation of large-volume gray-scale CT data into pore-volumes suitable for pore-scale modeling. Simulated permeabilities vary by three orders of magnitude and do not resemble observed values very well. Small-volume synchrotron-based CT data (such as produced by OSU) does not suffer significantly from segmentation ambiguities. 6. (UA) A standard properly parameterized Shan-Chen model LB model is useful for simulating porous media with proxy fluids as well as the scCO2/brine system and produces results that are consistent with tomographic observations. 7. (UA) A LB model with fluid-interactions defined by a (modified) Peng-Robinson Equation of State is able to handle the scCO2/brine system with variable solid phase wettability. This model is numerically stable at temperatures between 0 and 250 °C and pressures between 3 and 50 MPa, and produces appropriate densities above the critical point of CO2 and exhibits three-phase separation below. Based on above findings OSU and UA have proposed continued experimentation and pore-scale modeling of the scCO2/brine system. The reported research has extensively covered capillary trapping using proxy fluids, but due to limited beam-time availability we were unable to apply our high-pressure CO2 setup to sufficient variation in fluid properties, and initial scCO2 connectivity. New data will also allow us to test, calibrate and apply our LB models to reservoir conditions beyond those that are currently feasible experimentally. Such experiments and simulations will also allow us to provide information how suitable proxy fluids are for the scCO2/brine system. We believe it would be worthwhile to pursue the following new research questions: 1. What are the fundamental differences in the physics underlying capillary trapping at ambient vs. supercritical conditions? 2. Do newly developed pore-scale trapping interactions and relationships translate to continuum scales? A motivation for these questions was elaborated in “Capillary Trapping of Super-Critical CO2: Linking Pore and Continuum Scales to Verify new Relationships” that was submitted to DOE-BES in 2015.

20 FOSSIL-FUELED POWER PLANTS↗

A New Analytic-Adaptive Model for EGS Assessment, Development and Management Support

To increase understanding of the energy extraction capacity of Enhanced Geothermal System(s) (EGS), a numerical model development and application project is completed. The general objective of the project is to develop and apply a new, data-coupled Thermal-Hydrological-Mechanical-Chemical (T-H-M-C) model in which the four internal components can be freely selected from existing simulation software without merging and cross-combining a diverse set of computational codes. Eight tasks are completed during the project period. The results are reported in five publications, an MS thesis, twelve quarterly, and two annual reports to DOE. Two US patents have also been issued during the project period, with one patent application originated prior to the start of the project. The “Multiphase Physical Transport Modeling Method and Modeling System” (U.S. Patent 8,396,693 B2, 2013), a key element in the GHE sub-model solution, is successfully used for EGS studies. The “Geothermal Energy Extraction System and Method" invention (U.S. Patent 8,430,166 B2, 2013) originates from the time of project performance, describing a new fluid flow control solution. The new, coupled T-H-M-C numerical model will help analyzing and designing new, efficient EGS systems.

15 GEOTHERMAL ENERGY↗

A New Analytic-Adaptive Model for EGS Assessment, Development and Management Support

To increase understanding of the energy extraction capacity of Enhanced Geothermal System(s) (EGS), a numerical model development and application project is completed. The general objective of the project is to develop and apply a new, data-coupled Thermal-Hydrological-Mechanical-Chemical (T-H-M-C) model in which the four internal components can be freely selected from existing simulation software without merging and cross-combining a diverse set of computational codes. Eight tasks are completed during the project period. The results are reported in five publications, an MS thesis, twelve quarterly, and two annual reports to DOE. Two US patents have also been issued during the project period, with one patent application originated prior to the start of the project. The “Multiphase Physical Transport Modeling Method and Modeling System” (U.S. Patent 8,396,693 B2, 2013), a key element in the GHE sub-model solution, is successfully used for EGS studies. The “Geothermal Energy Extraction System and Method" invention (U.S. Patent 8,430,166 B2, 2013) originates from the time of project performance, describing a new fluid flow control solution. The new, coupled T-H-M-C numerical model will help analyzing and designing new, efficient EGS systems.

15 GEOTHERMAL ENERGY↗

Integrated Lift- and Wall-Force Closures in Computational Fluid Dynamics (CMFD)

In this report, we provide the necessary background to put this work in context and for a more detailed background, the reader is referred to CASL 2013 report. Central constraints on the design of LWRs are thermal-structural limits, e.g., maximum fuel pin temperature and maximum allowable working pressure in a pressure vessel. In both of these cases, the knowledge of the spatially and temporally varying two-phase flow is required to make accurate predictions as the distribution of both the phases have an impact on heat extraction as well as neutron moderation. Engineered two-phase flow systems have complex geometry and often high Reynolds numbers which increase the computational cost of an analysis, and generally make direct numerical simulations of the flow field impossible even on the fastest computers today. Consequently, a number of models which seek to reduce this computational cost have been created over the years. One of these called the two-fluid model has become the standard in three dimensional CFD (Computational Fluid Dynamics) codes and is and is also referred to as CMFD (Computational Multiphase Fluid Dynamics). The two-fluid model splits each conserved quantity into two fields, one for each phase that is treated as interspersed continua. Splitting one field into two requires additional interfacial closure relationships. In the two-fluid model, the momentum flux between the phases is governed by a number of these closure relationships. Generally, the net force on the dispersed phase is decomposed into inertial, added mass, buoyancy, drag, lift, and wall forces. In addition, there is a time-dependent Basset force, as well as turbulent dispersion effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Cloud Coalescence: A Dynamical Instability Affecting Multiphase Environments

Mass and size distributions are the key characteristics of any astrophysical object, including the densest clumps comprising the cold phase of multiphase environments. In our recent papers, we showed how individual clouds of various sizes form and evolve in active galactic nuclei. In particular, we showed that large clouds undergo damped oscillations as a response to their formation process. Here we follow up this investigation, addressing how different size clouds interact. We find that smaller clouds become trapped in the advective flows generated by larger clouds. The explanation for this behavior leads to a rather remarkable conclusion: even in the absence of gravity, complexes of clouds are dynamically unstable. In an idealized environment (e.g., one free of turbulence and magnetic fields) a perfectly symmetric arrangement of static clouds will remain static, but any small spatial perturbation will lead to all clouds coalescing into a single, large cloud, given enough time. Here, by using numerical simulations, we investigate the main factors that determine the rate of coalescence. In addition to the cloud separation distance, we find that the transient response of clouds to a disturbance is the primary factor. Turbulent motions in the flow can easily suppress this tendency for spatially well-separated clouds to coalesce, so it is as yet unclear if this phenomenon can occur in nature. Nevertheless, this Letter casts strong doubts on a recent hypothesis that large clouds are prone to fragmentation.

79 ASTRONOMY AND ASTROPHYSICS↗

Experiments and CFD simulation of mass transfer and hydrodynamics in a cylindrical bubble column

Flow visualization, Particle Image Velocimetry (PIV) and dissolved oxygen charging experiments are carried out in a small-scale cylindrical air-water bubble column at different superficial gas velocities. Void fraction profile, total gas holdup, liquid phase velocity field and volumetric mass transfer coefficient are obtained. Details of each experiment are given, with a new procedure for obtaining a 2D void fraction profile from flow visualization by image processing and time averaging which could provide useful information for CFD model validation. The experimental results are analyzed and compared with simulations where good agreement is found for both local and global flow patterns. The averaged governing equations and constitutive relations used in the simulation are discussed briefly, with emphasis on the proper comparison between simulations and experiments based on the derivation of the two-fluid model. Finally, by conducting three sets of experiments related to phase distribution, velocity field and mass transfer in the same apparatus, a validated CFD model for multiphase mass transfer is established and can be used in the future design of a multi-phase reactor.

42 ENGINEERING↗

Vadose Zone Transport of Tritium and Nitrate under Ponded Water Conditions

Vadose zone transport of tritium and nitrate can be important considerations at radioactive waste sites, landfills, or areas with industrial impacts. These contaminants are of particular concern because they typically have a relatively higher mobility in the subsurface compared to other compounds. Here, we describe a semiarid site with tritium and nitrate contamination involving a manmade ponded water source above a thick unsaturated zone at Los Alamos National Laboratory in New Mexico, USA. This study demonstrates the value of vadose zone flow and transport modeling for the development of field investigation plans (i.e., identifying optimal borehole locations and depths for contaminant characterization), and how a combination of modeling with isotope and geochemical measurements can provide insight into how tritium and nitrate transport in the vadose zone in semiarid environments. Modeling results suggest that at this location, tritium transport is well predicted by classical multiphase theory. Our work expands the demonstrated usefulness of a standard tritium conceptual model to sites with ponded surface conditions and agrees with previous results where a standard model was able to explain the evolution of a tritium plume at an arid waste disposal site. In addition, depth-based analyses of δ 18 O and δ 2 H of pore waters helped confirm the extent of pond infiltration into the vadose zone, and the δ 15 N of nitrate showed that the contaminant release history of the site was different than originally assumed.

58 GEOSCIENCES↗

Down-the-barrel and transverse observations of the Large Magellanic Cloud: Evidence for a symmetric galactic wind on the near and far sides of the galaxy

We compare the properties of gas flows on both the near and far side of the Large Magellanic Cloud (LMC) disk using Hubble Space Telescope UV absorption-line observations toward an active galactic nucleus behind (transverse) and a star within (down-the-barrel) the LMC disk at an impact parameter of 3.2 kpc. We find that even in this relatively quiescent region gas flows away from the disk at speeds up to ∼100 km s{sup −1} in broad and symmetrical absorption in the low and high ions. The symmetric absorption profiles combined with previous surveys showing little evidence that the ejected gas returns to the LMC and provide compelling evidence that the LMC drives a global, large-scale outflow across its disk, which is the likely result of a recent burst of star formation in the LMC. We find that the outflowing gas is multiphase, ionized by both photoionization (Si ii and Si iii) and collisional ionization (Si iv and C iv). We estimate a total mass and outflow rate to be ≳10{sup 7} M{sub ⊙} and ≳ 0.4 M{sub ⊙} yr{sup −1}. Since the velocity of this large-scale outflow does not reach the LMC escape velocity, the gas removal is likely aided by either ram-pressure stripping with the Milky Way (MW) halo or tidal interactions with the surrounding galaxies, implying that the environment of LMC-like or dwarf galaxies plays an important role in their ultimate gas starvation. Finally we reassess the mass and plausible origins of the high-velocity complex toward the LMC given its newly determined distance that places it in the lower MW halo and sky-coverage that shows it extends well beyond the LMC disk.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiphase eruptions associated with the craters Tycho and Aristarchus.

In 1967 part of Tycho and all of Aristarchus were photographed from Orbiter 5 with a ground resolution of about 5 m. Geologic units of Tycho are discussed together with crater counts in and near to Tycho, counts of small craters near to Surveyor 7, geologic units of Aristarchus, and crater counts in and near to Aristarchus. It is concluded that both Tycho and Aristarchus probably originated in the same manner, although there were differences in their later development. Extensive lava flows issued from the inner and outer rims of both craters. Flows from small craters in the vicinity of the Tycho and Aristarchus structures show that these small craters are volcanic. This establishes that one class of lunar volcano, at least, resembles an impact crater.

Strom, R. G.↗

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the Volume of Fluid (VOF) and Sharp Interface multiphase 2D axisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. Turbulence modeling; turbulence damping at the interface; and using constant vs. temperature dependent fluid properties were shown to have the most profound influence on predicted tank pressures and temperatures. Current study indicates that including phase change at the interface into the computational model is crucial for accurate prediction of the tank self-pressurization process. The effect of accommodation coefficient was also studied. Tank pressure values predicted by the VOF, and Sharp Interface models are within 4% of the experimental ones.

Computational Fluid Dynamics↗

Fluid-based gating mechanism with tunable multiphase selectivity and antifouling behavior

A gating mechanism that uses a capillary stabilized liquid as a reversible, reconfigurable gate that fills and seals pores in the closed state, and creates a non-fouling, liquid-lined pore in the open state is disclosed. Theoretical modeling and experiments demonstrate that for each transport substance, the gating threshold—the pressure needed to open the pores—can be rationally tuned over a wide pressure range. This enables realizing in one system differential response profiles for a variety of liquids and gases, even letting liquids flow through the pore while preventing gas from escaping. These capabilities allow dynamic modulation of gas-liquid sorting and to separate multi-phase mixtures, with the liquid lining ensuring sustained antifouling behavior. Because the liquid gating strategy enables efficient short-term and long-term operation and can be applied to a variety of pore structures and membrane materials, and to nano, micro as well as macroscale fluid systems, the gating systems is useful in a wide range of applications.

36 MATERIALS SCIENCE↗

Lattice Boltzmann modeling of self-propelled Leidenfrost droplets on ratchet surfaces

Here in this paper, the self-propelled motion of Leidenfrost droplets on ratchet surfaces is numerically investigated with a thermal multiphase lattice Boltzmann model with liquid-vapor phase change. The capability of the model for simulating evaporation is validated via the D 2 law. Using the model, we first study the performances of Leidenfrost droplets on horizontal ratchet surfaces. It is numerically shown that the motion of self-propelled Leidenfrost droplets on ratchet surfaces is owing to the asymmetry of the ratchets and the vapor flows beneath the droplets. It is found that the Leidenfrost droplets move in the direction toward the slowly inclined side from the ratchet peaks, which agrees with the direction of droplet motion in experiments [Linke et al., Phys. Rev. Lett., 2006, 96, 154502]. Moreover, the influences of the ratchet aspect ratio are investigated. For the considered ratchet surfaces, a critical value of the ratchet aspect ratio is approximately found, which corresponds to the maximum droplet moving velocity. Furthermore, the processes that the Leidenfrost droplets climb uphill on inclined ratchet surfaces are also studied. Lastly, numerical results show that the maximum inclination angle at which a Leidenfrost droplet can still climb uphill successfully is affected by the initial radius of the droplet.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

RealGasBrine v1.0 option of TOUGH+ v1.5

RealGasBrine v1.0 is a numerical code that for the simulation of the behavior of gas-bearing porous and/fractured geologic media. It is an option of TOUGH+ v1.5 [Moridis, 2014], a successor to the TOUGH2 [Pruess et al., 1999; 2012] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. RealGasBrine v1.0 needs the TOUGH+ v1.5 core code in order to compile and execute. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. RealGasBrine v1.0 describes the non-isothermal two- (for pure water) or three-phase (for brine) flow of an aqueous phase and a real gas mixture in a gas-bearing medium, with a particular focus in ultra-tight (such as tight-sand and shale gas) systems. Up to 12 individual real gases can be tracked, and salt can precipitate as solid halite. The capabilities of the code include coupled flow and thermal effects, real gas behavior, Darcy and non-Darcy flow, several isotherm options of gas sorption onto the grains of the porous media, complex fracture descriptions, gas solubility into water, and geomechanical effects on flow properties. RealGasBrine v1.0 allows the study of flow and transport of fluids and heat over a wide range of time frames and spatial scales not only in gas reservoirs, but also in any problem involving the flow of gases in geologic media, including the geologic storage of greenhouse gas mixtures, the behavior of geothermal reservoirs with multi-component condensable (H2O and CO2) and non-condensable gas mixtures, the transport of water and released H2 in nuclear waste storage applications, etc.

Moridis, GeorgeJ↗

Transport Phenomena During Equiaxed Solidification of Alloys

Recent progress in modeling of transport phenomena during dendritic alloy solidification is reviewed. Starting from the basic theorems of volume averaging, a general multiphase modeling framework is outlined. This framework allows for the incorporation of a variety of microscale phenomena in the macroscopic transport equations. For the case of diffusion dominated solidification, a simplified set of model equations is examined in detail and validated through comparisons with numerous experimental data for both columnar and equiaxed dendritic growth. This provides a critical assessment of the various model assumptions. Models that include melt flow and solid phase transport are also discussed, although their validation is still at an early stage. Several numerical results are presented that illustrate some of the profound effects of convective transport on the final compositional and structural characteristics of a solidified part. Important issues that deserve continuing attention are identified.

Beckermann, C.↗

Volatile Removal Assembly Flight Experiment and KC-135 Packed Bed Experiment: Results and Lessons Learned

The Volatile Removal Assembly (VRA) is a high temperature catalytic oxidation process that will be used as the final treatment for recycled water aboard the International Space Station (ISS). The multiphase nature of the process had raised concerns as to the performance of the VRA in a microgravity environment. To address these concerns, two experiments were designed. The VRA Flight Experiment (VRAFE) was designed to test a full size VRA under controlled conditions in microgravity aboard the SPACEHAB module and in a 1 -g environment and compare the performance results. The second experiment relied on visualization of two-phase flow through small column packed beds and was designed to fly aboard NASA's microgravity test bed plane (KC-135). The objective of the KC-135 experiment was to understand the two-phase fluid flow distribution in a packed bed in microgravity. On Space Transportation System (STS) flight 96 (May 1999), the VRA FE was successfully operated and in June 1999 the KC-135 packed bed testing was completed. This paper provides an overview of the experiments and a summary of the results and findings.

Holder, Donald W.↗