Computational Models for Fluid-to-Solid Transitions in Yield Stress Fluids.
Abstract not provided.
McConnell, Joshua↗
SEARCH · Engineering Papers
Results for “fluid”
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.
Generating Inviscid and Viscous Fluid Flow Simulations over a Surface Using a Quasi-simultaneous Technique
A fluid-flow simulation over a computer-generated surface is generated using a quasi-simultaneous technique. The simulation includes a fluid-flow mesh of inviscid and boundary-layer fluid cells. An initial fluid property for an inviscid fluid cell is determined using an inviscid fluid simulation that does not simulate fluid viscous effects. An initial boundary-layer fluid property a boundary-layer fluid cell is determined using the initial fluid property and a viscous fluid simulation that simulates fluid viscous effects. An updated boundary-layer fluid property is determined for the boundary-layer fluid cell using the initial fluid property, initial boundary-layer fluid property, and an interaction law. The interaction law approximates the inviscid fluid simulation using a matrix of aerodynamic influence coefficients computed using a two-dimensional surface panel technique and a fluid-property vector. An updated fluid property is determined for the inviscid fluid cell using the updated boundary-layer fluid property.
Sturdza, Peter↗
Ultrasound Thermal Field Imaging of Opaque Fluids
We have initiated an experimental program to develop an ultrasound system for non-intrusively imaging the thermal field in opaque fluids under an externally imposed temperature gradient. Many industrial processes involve opaque fluids, such as molten metals, semiconductors, and polymers, often in situations in which thermal gradients are important. For example, one may wish to understand semiconductor crystal growth dynamics in a Bridgman apparatus. Destructive testing of the crystal after the process is completed gives only indirect information about the fluid dynamics of the formation process. Knowledge of the coupled thermal and velocity fields during the growth process is then essential. Most techniques for non-intrusive velocity and temperature measurement in fluids are optical in nature, and hence the fluids studied must be transparent. In some cases (for example, LDV (laser Doppler velocimetry) and PIV (particle imaging velocimetry)) the velocities of small neutrally buoyant seed particles suspended in the fluid, are measured. Without particle seeding one can use the variation of the index of refraction of the fluid with temperature to visualize, through interferometric, Schlieren or shadowgraph techniques, the thermal field. The thermal field in turn gives a picture of the pattern existing in the fluid. If the object of study is opaque, non-optical techniques must be used. In this project we focus on the use of ultrasound, which propagates easily through opaque liquids and solids. To date ultrasound measurements have almost exclusively relied on the detection of sound scattered from density discontinuities inside the opaque material of interest. In most cases it has been used to visualize structural properties, but more recently the ultrasound Doppler velocimeter has become available. As in the optical case, it relies on seed particles that scatter Doppler shifted sound back to the detector. Doppler ultrasound techniques are, however, not useful for studying convective fluid flow in crystal growth, because particle seeding is unacceptable and flow velocities are typically too low to be resolved, and may be even lower in microgravity conditions where buoyancy forces are negligible. We will investigate a different use of ultrasound to probe the flows of opaque fluids non-intrusively and without the use of seed particles: our goal is to ultrasonically visualize the thermal field of opaque fluids with relatively high spatial resolution. The proposed technique relies upon the variation of sound speed with temperature of the fluid. A high frequency ultra-sound pulse passing through a fluid-filled chamber will traverse the chamber in a time determined by the relevant chamber dimension and the temperature of the fluid through which the pulse passes. With high time-resolution instrumentation that compares the excitation signal with the received pulse we can detect the influence of the fluid temperature on the pulse travel time. This is effectively an interferometric system, which in its optical form is an extremely sensitive approach to measuring thermal fields in fluids. Moreover, the temperature dependence of sound velocity in liquid metals is comparable to the temperature dependence of the speed of light required for accurate interferometric thermal images in transparent fluids. With an array of transducers scanned electronically a map of the thermal field over the chamber could be produced. An alternative approach would be to use the ultrasound analog of the shadowgraph technique. In the optical version, collimated light passes through the fluid, where it is focused or defocused locally by temperature field induced variations of the index of refraction. The resulting image reveals the thermal field through the spatial pattern of light intensity variations. By analogy, an ultrasound plane wave traversing an opaque fluid sample would be also locally focused or defocused depending on the speed of sound variations, giving rise to spatial variations in sound intensity that will reveal the thermal field pattern. These approaches could be applied to any situation in which temperature differences are expected to occur, and will rapidly provide information about the flow that simply cannot be obtained by any current intrusive or non-intrusive diagnostic technique. As materials processing in microgravity matures it will become increasingly important to have available simple and versatile diagnostic tools, such as we will develop, for studying the flows of opaque fluids under thermal forcing.
Andereck, C. David↗
Evaluation of Propylene Glycol-Based Fluids for Constellation Habitats and Vehicles
Two fluid life tests have been conducted to evaluate propylene glycol-based fluids for use in Constellation habitats and vehicles. The first test was conducted from November 2008 to January 2009 to help determine the compatibility of the propylene glycol-based fluid selected for Orion at the time. When the first test uncovered problems with the fluid selection, an investigation and selection of a new fluid were conducted. A second test was started in March 2010 to evaluate the new selection. For the first test, the fluid was subjected to a thermal fluid loop that had flight-like properties, as compared to Orion. The fluid loop had similar wetted materials, temperatures, flow rates, and aluminum wetted surface area to fluid volume ratio. The test was designed to last for 10 years, the life expectancy of the lunar habitat. However, the test lasted less than two months. System filters became clogged with precipitate, rendering the fluid system inoperable. Upon examination of the precipitate, it was determined that the precipitate composition contained aluminum, which could have only come from materials in the test stand, as aluminum is not part of the original fluid composition. Also, the fluid pH was determined to have increased from 10.1, at the first test sample, to 12.2, at the completion of the test. This high of a pH is corrosive to aluminum and was certainly a contributing factor to the development of precipitate. Due to the problems encountered during this test, the fluid was rejected as a coolant candidate for Orion. A new propylene glycol-based fluid was selected by the Orion project for use in the Orion vehicle. The Orion project has conducted a series of screening tests to help verify that there will be no problems with the new fluid selection. To compliment testing performed by the Orion project team, a new life test was developed to test the new fluid. The new test bed was similar to the original test bed, but with some improvements based on experience gained from the earlier test bed. The surface area of both aluminum and nickel in the test bed were designed to be similar to that of the Orion fluid loop, since the Orion fluid loop was expected to have high concentrations of both metals in the system. Also, additional sample materials were added to the test bed to match recent updates to materials selections for Orion. At the time of this paper publication, approximately five months of testing will have been completed. This paper gives a status of the testing completed to date.
Lee, Steve↗
Analysis of Direct Samples of Early Solar System Aqueous Fluids
Over the past three decades we have become increasingly aware of the fundamental importance of water, and aqueous alteration, on primitive solar-system bodies. Some carbonaceous and ordinary chondrites have been altered by interactions with liquid water within the first 10 million years after formation of their parent asteroids. Millimeter to centimeter-sized aggregates of purple halite containing aqueous fluid inclusions were found in the matrix of two freshly-fallen brecciated H chondrite falls, Monahans (1998, hereafter simply "Monahans") (H5) and Zag (H3-6) (Zolensky et al., 1999; Whitby et al., 2000; Bogard et al., 2001) In order to understand origin and evolution of the aqueous fluids inside these inclusions we much measure the actual fluid composition, and also learn the O and H isotopic composition of the water. It has taken a decade for laboratory analytical techniques to catch up to these particular nanomole-sized aqueous samples. We have recently been successful in (1) measuring the isotopic composition of H and O in the water in a few fluid inclusions from the Zag and Monahans halite, (2) mineralogical characterization of the solid mineral phases associated with the aqueous fluids within the halite, and (3) the first minor element analyses of the fluid itself. A Cameca ims-1270 equipped with a cryo-sample-stage of Hokkaido University was specially prepared for the O and H isotopic measurements. The cryo-sample-stage (Techno. I. S. Corp.) was cooled down to c.a. -190 C using liquid nitrogen at which the aqueous fluid in inclusions was frozen. We excavated the salt crystal surfaces to expose the frozen fluids using a 15 keV Cs+ beam and measured negative secondary ions. The secondary ions from deep craters of approximately 10 m in depth emitted stably but the intensities changed gradually during measurement cycles because of shifting states of charge compensation, resulting in rather poor reproducibility of multiple measurements of standard fluid inclusions of +/- 90 0/00(2 sigma) for delta D, and +/- 29 0/00 (2 sigma) for delta O-18. On the other hand, the reproducibility of Delta O-17 is plus or minus 8 /00 (2 sigma ) because the observed variations of isotope ratios follow a mass dependent fractionation law. Variations of delta D of the aqueous fluids range over sog,a 330(90; 2 sigma ) to +1200(90) 0/00 for Monahans and delta 300(96) 0/00 to +90(98)0/00 for Zag. Delta O-17 of aqueous fluids range over delta 16(22) 0/00 to +18(10) 0/00 for Monahans and +3(10) 0/00 to +27(11) 0/00 for Zag. These variations are larger than the reproducibility of standard analyses and suggest that isotope equilibria were under way in the fluids before trapping into halite. The mean values of delta D and Delta O-17 are +290 0/00 and +9 0/00, respectively. The mean values and the variations of these fluids are different from the representative values of ordinary chondrites, verifying our working hypothesis that the fluid inclusion-bearing halites were not indigenous to the H chondrite parent-asteroid but rather represent exogenous material delivered onto the asteroid from a separate cryovolcanically-active body. This initial isotopic work has demonstrated the feasibility of the measurements, but also revealed sample processing and analytical shortcomings that are now being addressed. Examination of solid mineral inclusions within Monahans and Zag halite grains by confocal Raman spectroscopy at the Carnegie Geophysical Laboratory has revealed them to be metal, magnetite, forsteritic olivine (Fo.98), macromolecular carbon (MMC), pyroxenes, feldspar with Raman spectral affinity to anorthoclase and, probably, fine-grained lepidocrocite (FeO(OH)). In addition, one inclusion features aliphatic material with Raman spectral features consistent with a mixture of short-chain aliphatic compounds. We have initiated analyses of the bulk composition of the fluids within the inclusions in Zag and Monahans halites at Virginia Tech by LA ICPMS using angilent 7500ce quadrupole ICPMS and a Lambda Physik GeoLas 193 nm excimer laser ablation system. Preliminary results reveal that the inclusion aqueous fluids contain highly charged cations of Ca, Mg and Fe. The minerals and compounds discovered thus far within Monahans/Zag halites are indicative of an originating body at least partly composed of unequilibrated anhydrous materials (high Fo olivine, pyroxenes, feldspars, possibly the metal) which were subjected to aqueous alteration (the halite parent brine) and containing a light organic component (the short-chain aliphatic compounds). This material was ejected from the originating body with little or no disruption, as evidenced with the presence of fluid inclusions. An actively geysering body similar to modern Enceladus (Postberg et al., 2011) may be a reasonable analogue in this respect. Also, the originating body should have been within close proximity to the H chondrite parent in order to generate the number of halite grains seen in Monahans and Zag. Other candidates for Monahans/Zag halite parent bodie(s) may include a young Ceres with its possible liquid ocean, or Main Belt comets.
Zolensky, Michael E.↗
Reaction path model of the formation of abiotic immiscible hydrocarbon fluids in subducted carbonated serpentinites, Lanzo Massif (Western Italian Alps)
Fluids generated from subducted slabs participate in the cycling of deep carbon in the crust and upper mantle. In these fluids, aqueous carbon species vary in oxidation state between +IV and -IV depending on whether the fluids are oxidizing or reducing, respectively. Most studies of subduction-zone fluids have focused on oxidized carbon species. However, recent studies of natural samples have demonstrated the occurrence of deep, reducing fluids, generated in both subducted oceanic upper mantle and crustal rocks. CH 4 -H 2 -rich fluid inclusions in subducted carbonated serpentinites have demonstrated the existence of abiotic, immiscible, hydrocarbon fluids at upper mantle conditions. To investigate the formation of such immiscible hydrocarbon fluids during the evolution of subducted carbonated serpentinites, we used equilibrium constants from the Deep Earth Water model to carry out predictive chemical mass transfer modeling to simulate the alteration reactions. A novel feature of the models was the inclusion of an immiscible hydrocarbon fluid containing six components (CH 4,f , C 2 H 6,f , C 3 H 8,f , isoC 4 H 10,f , CO 2,f , H 2,f ). This feature enabled prediction of the formation of a separate immiscible fluid in equilibrium with aqueous species and minerals. We developed a predictive reaction path model of invasive H 2,f reacting with carbonated serpentinites and interstitial aqueous fluids for comparison with the natural samples from the Lanzo Massif, western Italian Alps. Over a range of temperatures and pressures, immiscible hydrocarbon fluids formed in association with altered mineral assemblages. Reaction progress caused the transformation of carbonated serpentinites and the formation of clinopyroxene, brucite, graphite, and hydrocarbon fluids, along with changes of pH, logfO 2 , and aqueous species. CH 4,f was the most abundant hydrocarbon species in all the models. The overall results at 2.0 GPa and 400 to 450 °C were consistent with the natural samples from the Lanzo Massif. Interestingly, large amounts of H 2 O formed due to oxidation of H 2 . More hydrocarbons and H 2 O formed in models with lower fluid/rock mass ratios or with more reactant H 2 . Models at different pressure and temperature conditions showed similar results with some variation in the relative stabilities of aragonite, graphite and olivine solid solution, and associated differences in mineral sequences, hydrocarbon fluids, values of aqueous species, and the final log fO 2 and pH. As a result, our models strongly support the laboratory and field evidence that reduction of carbonated serpentinites by infiltrating H 2 fluids can cause the formation of immiscible, abiotic hydrocarbon fluids in subduction zones.
58 GEOSCIENCES↗
Fluid Film Bearing Code Development
The next generation of rocket engine turbopumps is being developed by industry through Government-directed contracts. These turbopumps will use fluid film bearings because they eliminate the life and shaft-speed limitations of rolling-element bearings, increase turbopump design flexibility, and reduce the need for turbopump overhauls and maintenance. The design of the fluid film bearings for these turbopumps, however, requires sophisticated analysis tools to model the complex physical behavior characteristic of fluid film bearings operating at high speeds with low viscosity fluids. State-of-the-art analysis and design tools are being developed at the Texas A&M University under a grant guided by the NASA Lewis Research Center. The latest version of the code, HYDROFLEXT, is a thermohydrodynamic bulk flow analysis with fluid compressibility, full inertia, and fully developed turbulence models. It can predict the static and dynamic force response of rigid and flexible pad hydrodynamic bearings and of rigid and tilting pad hydrostatic bearings. The Texas A&M code is a comprehensive analysis tool, incorporating key fluid phenomenon pertinent to bearings that operate at high speeds with low-viscosity fluids typical of those used in rocket engine turbopumps. Specifically, the energy equation was implemented into the code to enable fluid properties to vary with temperature and pressure. This is particularly important for cryogenic fluids because their properties are sensitive to temperature as well as pressure. As shown in the figure, predicted bearing mass flow rates vary significantly depending on the fluid model used. Because cryogens are semicompressible fluids and the bearing dynamic characteristics are highly sensitive to fluid compressibility, fluid compressibility effects are also modeled. The code contains fluid properties for liquid hydrogen, liquid oxygen, and liquid nitrogen as well as for water and air. Other fluids can be handled by the code provided that the user inputs information that relates the fluid transport properties to the temperature.
Source record↗
Hydrophobicity of Cryogenic Fluids for Fuel Transfer in Space Applicaitons
Introduction Hydrophobicity is the tendency of a fluid to repel another material or fluid. Hydrophobiciy can be measured through the wetting angle of the fluid on the surface of the other material. The greater the wetting angle, the more hydrophobic the surface is against the fluid. Hydrophobicity properties can be caused by one of two means, either through physical properties or chemical properties [\citenum{NARBUTT2020121}]. Physical hydrophobicity due to the surface of a solid creating surface roughness/patterning, minimizing the contact area between the fluid and the surface. This can be observed throughout nature, such as the surface of lotus leaves or butterfly wings. This physical property can be induced on a variety of surfaces, namely through laser etching, allowing a surface to be finely lineated to imitate these natural surfaces while controlling quantity, depth, and patterns of the etching on the surface [\citenum{10.1063/1.4905616}]. Physical hydrophobicity is therefore dependent on several variables, including the surface that the fluid is on and properties of the fluid itself, such as density and surface tension, both of which are dependent on temperature and/or pressure. Chemical hydrophobicity [\citenum{MadeiraHydrophobocicicici}] is due to the inherent chemical properties of the materials being used. This is most commonly due to the molecular structures changing the polarity of the materials. Depending on a nonpolar material will repel a polar material, proving to be hydrophobic, conversely if the polarity of the materials is the same (polar – polar, nonpolar – nonpolar), they will attract each other. Cryogenics refers to the behavior of materials at very low temperatures (<\ang{-100} C) [\citenum{ZOHURI20181}]. In space applications many things are inherently cryogenic, therefore this is an important field. There has been little research in how hydrophobicity changes at cryogenic temperatures. Cryogenic fluids are commonly used as fuels for spacecraft, therefore, integrating a hydrophobic surface can increase the transfer rate of the fuel. To test this, several experiments were set up to determine the hydrophobic properties of cryogenic fluids, including nitrogen (LN2), argon (LAr), oxygen (LOx), hydrogen (LH2), and methane (CH4), at cryogenic temperatures. Etched Wafer Testing A silicon wafer cut from a crystal of silicon [100] was used to model the potential hydrophobicity of various cryogenic fluids. Silicon [100] references to the crystallographic orientation of the silicon crystals in the wafer. These wafers were then laser etched to create a surface that is more likely to be hydrophobic. To test the hydrophobicity of the wafers at cryogenic temperatures, the temperatures of the wafers must be reduced to the same temperature as the cryogenic fluid being used to prevent rapid boil-off. To achieve this a double-walled vacuum insulated glass chamber was utilized. The chamber is open, with a double-walled glass that can be placed in a vacuum to remove condensation from the outside to allow easier viewing of the experiment. Furthermore, the silicon wafer's temperature must be lowered to the temperature of LN2, as well as maintain the temperature throughout the experiment. To achieve this a piece of 6061 aluminum was used, creating a stand-off for the wafer and it would allow for the insulation of the temperature of the wafers. The container was then filled with LN2 and once the LN2 stabilized the vaporization and the levels of LN2 dropped under the height of the wafer, the wafer was then allowed to air dry. Once the wafer was dry from the LN2, drops of LN2 were placed on the surface of the wafer for observation. During the first trial, the LN2 that was dropped on the surface displayed nonhydrophobic behaviors, spreading out along the surface, with a minimal wetting angle (too small to be measured). This process was repeated for liquid argon. Further testing with other cryogenic liquids will require different testing apparatuses due to being more volatile. Furthermore, several papers [\citenum{voltvolt7}] have suggested that running a voltage can induce a hydrophobic effect throughout a surface, further testing will include a voltage (constant and oscillating) to determine if voltage influences inducing hydrophobicity at cryogenic temperatures. Chemical Testing Coating the interior walls of the fuel tanks and fuel lines can successfully create a hydrophobic surface. The inherent problem is finding a material that can be used to coat the surface, furthermore, at these temperatures, the coating will remain solid, which could pose issues in maintaining the hydrophobic properties. To test this water (polar) is hydrophobic against oils and fats (nonpolar), therefore the same experiment as the etched wafer testing was conducted to determine if the hydrophobic properties will persist as the nonpolar material remains solid and the polar material remains a fluid. The water remained hydrophobic, allowing testing can expand to the cryogenic fluids. However, another hurdle is faced in finding a chemically opposite material to the cryogenic fluids being tested. As many of the cryogenic fluids being tested are diatomic, they are inherently nonpolar, therefore, the material used for the hydrophobic coating must be polar. Polar greases and lubricants are difficult to come by, however, lithium stearate, appears to be a potential candidate for creating a coated hydrophobic surface for cryogenic fluids at cryogenic temperatures.
Cryogenics↗
Two-way coupled particle-turbulence interaction: Effect of numerics and resolution on fluid and particle statistics
Euler-Lagrange point-particle simulation has emerged as a premier methodology for studying dispersed particle-laden flows. This method's popularity stems from its ability to resolve fine-scale fluid structures while also tracking individual particles at reduced cost using an appropriate particle acceleration model. However, the point-particle model has known convergence issues in that refinement of the fluid grid can lead to changes in the predicted statistics. The reasons for non-convergence are twofold: the point-particle two-way coupling force in the Navier-Stokes equations requires a numerical regularization and, without careful implementation, yields a singular force on the fluid with grid refinement. The second factor that yields grid-dependent statistics is that the point-particle force model in general depends on the undisturbed fluid velocity. When the undisturbed fluid velocity is not robustly modeled in a grid-insensitive way, the calculated force for both particles and fluid will be grid-dependent, contaminating their respective statistics. While the first issue regarding regularizing the point-particle source term has received attention in the literature, the consequences of robustly modeling the undisturbed velocity in the context of grid refinement of turbulence has received little attention. Here, we consider decaying homogeneous isotropic turbulence laden with particles at different Stokes numbers. For a given Stokes number, we systematically refine the grid and demonstrate that explicitly modeling the undisturbed fluid velocity yields relative grid insensitivity for the energy of the particle and fluid phases, as well as acceleration of the particles. We also demonstrate that an appropriately defined dissipation rate is also grid-insensitive when an undisturbed fluid velocity correction is used. In contrast, when the undisturbed fluid velocity is modeled using the conventional approach of interpolating the local fluid velocity to the particle location, we show this procedure yields divergent statistics with grid refinement. In particular, we show that higher-order interpolation of the fluid velocity in two-way coupled problems is worse than lower-order interpolation, in the absence of a correction procedure to estimate the undisturbed fluid velocity. We also examine velocity derivative statistics of the fluid phase and demonstrate that these statistics are not in general convergent even when the undisturbed fluid velocity is explicitly modeled. Collectively, the observations in this paper are used to present a philosophy on the types of questions which are answerable with point-particle methods.
42 ENGINEERING↗
Experimental Simulation of Shock Reequilibration of Fluid Inclusions During Meteorite Impact
Fluid inclusions are microscopic volumes of fluid trapped within minerals as they precipitate. Fluid inclusions are common in terrestrial minerals formed under a wide array of geological settings from surface evaporite deposits to kimberlite pipes. While fluid inclusions in terrestrial rocks are the rule rather than the exception, only few fluid inclusion-bearing meteorites have been documented. The rarity of fluid inclusions in meteoritic material may be explained in two ways. First, it may reflect the absence of fluids (water?) on meteorite parent bodies. Alternatively, fluids may have been present when the rock formed, but any fluid inclusions originally trapped on the parent body were destroyed by the extreme P-T conditions meteorites often experience during impact events. Distinguishing between these two possibilities can provide significant constraints on the likelihood of life on the parent body. Just as textures, structures, and compositions of mineral phases can be significantly altered by shock metamorphism upon hypervelocity impact, fluid inclusions contained within component minerals may be altered or destroyed due to the high pressures, temperatures, and strain rates associated with impact events. Reequilibration may occur when external pressure-temperature conditions differ significantly from internal fluid isochoric conditions, and result in changes in fluid inclusion properties and/or textures. Shock metamorphism and fluid inclusion reequilibration can affect both the impacted target material and the meteoritic projectile. By examining the effects of shock deformation on fluid inclusion properties and textures we may be able to better constrain the pressure-temperature path experienced by shocked materials and also gain a clearer understanding of why fluid inclusions are rarely found in meteoritic samples.
Madden, M. E. Elwood↗
Fluid channeling system
A fluid channeling system includes a fluid ejector, a heat exchanger, and a fluid pump disposed in series flow communication The ejector includes a primary inlet for receiving a primary fluid, and a secondary inlet for receiving a secondary fluid which is mixed with the primary fluid and discharged therefrom as ejector discharge. Heat is removed from the ejector discharge in the heat exchanger, and the heat exchanger discharge is compressed in the fluid pump and channeled to the ejector secondary inlet as the secondary fluid In an exemplary embodiment, the temperature of the primary fluid is greater than the maximum operating temperature of a fluid motor powering the fluid pump using a portion of the ejector discharge, with the secondary fluid being mixed with the primary fluid so that the ejector discharge temperature is equal to about the maximum operating temperature of the fluid motor.
Davis, Donald Y.↗
Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions
This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.
Majumdar, A. K.↗
Axisymmetric flows from fluid injection into a confined porous medium
In this work, we study the axisymmetric flows generated from fluid injection into a horizontal confined porous medium that is originally saturated with another fluid of different density and viscosity. Neglecting the effects of surface tension and fluid mixing, we use the lubrication approximation to obtain a nonlinear advection-diffusion equation that describes the time evolution of the sharp fluid-fluid interface. The flow behaviors are controlled by two dimensionless groups: M , the viscosity ratio of displaced fluid relative to injected fluid, and Γ , which measures the relative importance of buoyancy and fluid injection. For this axisymmetric geometry, the similarity solution involving R 2 /T (where R is the dimensionless radial coordinate and T is the dimensionless time) is an exact solution to the nonlinear governing equation for all times. Four analytical expressions are identified as asymptotic approximations (two of which are new solutions): (i) injection-driven flow with the injected fluid being more viscous than the displaced fluid ( Γ $\ll$ 1 and M < 1) where we identify a self-similar solution that indicates a parabolic interface shape; (ii) injection-driven flow with injected and displaced fluids of equal viscosity ( Γ $\ll$ 1 and M = 1), where we find a self-similar solution that predicts a distinct parabolic interface shape; (iii) injection-driven flow with a less viscous injected fluid ( Γ $\ll$ 1 and M > 1) for which there is a rarefaction wave solution, assuming that the Saffman-Taylor instability does not occur at the reservoir scale; and (iv) buoyancy-driven flow ( Γ $\gg$ 1) for which there is a well-known self-similar solution corresponding to gravity currents in an unconfined porous medium [S. Lyle et al. “Axisymmetric gravity currents in a porous medium,” J. Fluid Mech. 543, 293–302 (2005)]. The various axisymmetric flows are summarized in a Γ - M regime diagram with five distinct dynamic behaviors including the four asymptotic regimes and an intermediate regime. The implications of the regime diagram are discussed using practical engineering projects of geological CO 2 sequestration, enhanced oil recovery, and underground waste disposal.
42 ENGINEERING↗
Methodologies for Reservoir Characterization Using Fluid Inclusion Gas Chemistry
The purpose of this project was to: 1) evaluate the relationship between geothermal fluid processes and the compositions of the fluid inclusion gases trapped in the reservoir rocks; and 2) develop methodologies for interpreting fluid inclusion gas data in terms of the chemical, thermal and hydrological properties of geothermal reservoirs. Phase 1 of this project was designed to conduct the following: 1) model the effects of boiling, condensation, conductive cooling and mixing on selected gaseous species; using fluid compositions obtained from geothermal wells, 2) evaluate, using quantitative analyses provided by New Mexico Tech (NMT), how these processes are recorded by fluid inclusions trapped in individual crystals; and 3) determine if the results obtained on individual crystals can be applied to the bulk fluid inclusion analyses determined by Fluid Inclusion Technology (FIT). Our initial studies however, suggested that numerical modeling of the data would be premature. We observed that the gas compositions, determined on bulk and individual samples were not the same as those discharged by the geothermal wells. Gases discharged from geothermal wells are CO 2 -rich and contain low concentrations of light gases (i.e. H 2 , He, N, Ar, CH4). In contrast many of our samples displayed enrichments in these light gases. Efforts were initiated to evaluate the reasons for the observed gas distributions. As a first step, we examined the potential importance of different reservoir processes using a variety of commonly employed gas ratios (e.g. Giggenbach plots). The second technical target was the development of interpretational methodologies. We have develop methodologies for the interpretation of fluid inclusion gas data, based on the results of Phase 1, geologic interpretation of fluid inclusion data, and integration of the data. These methodologies can be used in conjunction with the relevant geological and hydrological information on the system to create fluid models for the system. The hope is that the methodologies developed will allow bulk fluid inclusion gas analysis to be a useful tool for estimating relative temperatures, identifying the sources and origins of the geothermal fluids, and developing conceptual models that can be used to help target areas of enhanced permeability.
15 GEOTHERMAL ENERGY↗
Chemical and isotopic evolution of flowback fluids from the Utica Gas Shale Play, Eastern Ohio USA
For this work, hydraulic fracturing flowback fluids were collected from two Utica/Point Pleasant well pads in eastern Ohio. One site was in the wet gas zone (UPPW), while the site ~50 km south consisted of four wells on the same pad in the dry gas zone (UPPS). Samples of input fluids also were collected before and during hydraulic fracturing. Flowback fluids are Na-Ca-Cl brines with total dissolved salt (TDS) concentrations that increase over several months from ~100 to 200 g/L. The slightly higher TDS of the dry gas fluids are in part due to recycled flowback used as input fluids, and in part due to lower volume of water used in hydraulic fracturing. Concentrations of most major ions (Ca 2+ , Mg 2+ , Na + , Sr2 + , Fe, Mn, Cl – , Br – ) are similar for the five wells sampled, although small but systematic changes occur in the major element ratios over time. Most notably, an increase in the Sr/Cl ratio corresponds to an increase in the 87 Sr/ 86 Sr ratio in the fluids, suggesting interactions with a more radiogenic Sr source in the subsurface. Dissolved Ba concentrations and Ra activities were different between the two sites, reflecting a high SO 4 2– fluid used at the UPPW site, and water-rock reactions occurring during hydraulic fracturing at the UPPS site. Water oxygen (δ 18 O) and hydrogen (δD) isotopes for input fluids used in the UPPW4 well and fresh water used for the UPPS wells fall on the Global Meteoric waterline (GMWL). Flowback fluids from both sites are relatively enriched in δ 18 O and δD compared to the input, but do not appear to follow a simple mixing trend, suggesting reaction and isotopic exchange with carbonates and fractionation due to imbibition in the rock. Cl isotopes, δ 37 Cl in the FP fluids varied from ~ –0.43 to +0.13‰, the largest variation was observed in the earlier stages of flowback, while in the later stages δ 37 Cl exhibited a small but systematic increase over time, suggesting diffusion control of isotopic composition. Some of the trace species measured (dissolved Fe, Mn, thiosulfate and organic acids) do not follow the same trends as the major ions, suggesting contributions from input fluids, or microbially mediated reactions are exerting control on their concentrations. Input water chemistry exerts an important control on the Sr and Ba concentrations in flowback water. High SO 4 2– in the input fluids used for hydraulically fracturing the wet gas well leads to precipitation of barite-celestite in flowback fluids.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A study of fluid-structure interaction on overlapping grids with a focus on flexible beams
Fluid-structure interaction problems appear in many areas of science and engineering. Adequately understanding these problems is crucial to the study of such disciplines as aeroelastic flutter, wind turbine performance, implosions and explosions, and hemodynamics, to name a few. One technique to explore these problems is the use of numerical simulations. Although methods for the accurate and efficient numerical simulation of many problems involving only flexible structures or those involving only fluids have existed for a long time, the simulation of coupled problems, in which a flexible structure (or structures) interact with fluids remains challenging. The source of the difficulty with these problems is the construction of schemes for the coupling between the fluid and the solid that are both accurate and stable. There are two general approaches to designing a coupling procedure. In the first, the so-called monolithic approach, a single solver is implemented that handles both the fluid equations and the solid equations in a single framework. In this scheme, the solver handles the fluid and solid equations simultaneously, resulting in a strong coupling between the two. The second approach, named the partitioned approach, handles the fluid and structure domains separately. Whereas in the monolithic scheme there is a single solver for both the fluid and the structure, in the partitioned scheme there are separate solvers for the fluid and the structure. This separation heavily restricts the communication and coordination that can occur between the two solvers. Typically the only communication done in these schemes is through the application of boundary conditions (where the fluid applies a boundary condition on the fluid, and the solid one (or more) on the fluid). This makes the implementation of partitioned schemes easy, especially because they can take advantage of existing algorithms and software developed for the simulation of fluid flow or computational mechanics. However, they have a weakness – accuracy and stability are much harder to achieve. Stability is a particular problem for these schemes, as they are prone to the so-called added-mass effect.
97 MATHEMATICS AND COMPUTING↗
Method and Apparatus for Measuring Fluid Flow
The invention is a method and apparatus for monitoring the presence, concentration, and the movement of fluids. It is based on utilizing electromagnetic measurements of the complex permittivity of the fluids for detecting and monitoring the fluid. More particularly the apparatus uses one or more microwave probes which are placed at the locations where the measurements are to be made. A radio frequency signal is transmitted to the probe and the reflected signal is phase and amplitude detected at a rapid rate for the purpose of identifying the fluids, based on their dielectric constant at the probe. The apparatus can be used for multiple purposes including measures of flow rates, turbulence, dispersion, fluid identification, and changes in flow conditions of multiple fluids or multiple states of a single fluid in a flowline or a holding container. The apparatus includes a probe consisting of two electrical conductors separated by an insulator. A radio frequency signal is communicated to the probe and is reflected back from the portion of the probe exposed to the fluid. The radio frequency signal also provides a reference signal. An oscillator generates a second signal which combined with each of the reference signal and the reflected signal to produce signals of lower frequencies to facilitate filtering and amplifying those signals. The two signals are then mixed in a detector to produce an output signal that is representative of the phase and amplitude change caused by the reflection of the signal at the probe exposed to the fluid. The detector may be a dual phase detector that provides two such output signals that are in phase quadrature. A phase shifter may be provided for selectively changing the phase of the reference signal to improve the sensitivity of at least one of the output signals for more accurate readings and/or for calibration purposes. The two outputs that are in quadrature with respect to each other may be simultaneously monitored to account for drift errors. The output signals are digitized and provided to a computer at a sample rate which may be very high. The computer is operable to identify the fluid based on its complex permittivity as may be useful for identifying the flow rates, determining the fluid mixture ratio, detecting impurities in the fluid, and so forth. Novelty is believed to reside in the use of the real part of complex permittivity to measure small difference in permittivity of the fluid.
Arndt, G. Dickey↗
A Geophysical Flow Experiment in a Compressible Critical Fluid
The first objective of this experiment is to build an experimental system in which, in analogy to a geophysical system, a compressible fluid in a spherical annulus becomes radially stratified in density through an A.C. electric field. When this density gradient is demonstrated, the system will be augmented so that the fluid can be driven by heating and rotation and tested in preparation for a microgravity experiment. This apparatus consists of a spherical capacitor filled with critical fluid in a temperature controlled environment. To make the fluid critical, the apparatus will be operated near the critical pressure, critical density, and critical temperature of the fluid. This will result in a highly compressible fluid because of the properties of the fluid near its critical point. A high voltage A.C. source applied across the capacitor will create a spherically symmetric central force because of the dielectric properties of the fluid in an electric field gradient. This central force will induce a spherically symmetric density gradient that is analogous to a geophysical fluid system. To generate such a density gradient the system must be small (approx. 1 inch diameter). This small cell will also be capable of driving the critical fluid by heating and rotation. Since a spherically symmetric density gradient can only be made in microgravity, another small cell, of the same geometry, will be built that uses incompressible fluid. The driving of the fluid by rotation and heating in these small cells will be developed. The resulting instabilities from the driving in these two systems will then be studied. The second objective is to study the pattern forming instabilities (bifurcations) resulting from the well controlled experimental conditions in the critical fluid cell. This experiment will come close to producing conditions that are geophysically similar and will be studied as the driving parameters are changed.
Hegseth, John↗