Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID 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,567 records · Page 87

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗

Design of Fiber Optic Sensors for Measuring Hydrodynamic Parameters

The science of optical hydrodynamics involves relating the optical properties to the fluid dynamic properties of a hydrodynamic system. Fiber-optic sensors are being designed for measuring the hydrodynamic parameters of various systems. As a flowing fluid makes an encounter with a flat surface, it forms a boundary layer near this surface. The region between the boundary layer and the flat plate contains information about parameters such as viscosity, compressibility, pressure, density, and velocity. An analytical model has been developed for examining the hydrodynamic parameters near the surface of a fiber-optic sensor. An analysis of the conservation of momentum, the continuity equation and the Navier-Stokes equation for compressible flow were used to develop expressions for the velocity and the density as a function of the distance along the flow and above the surface. When examining the flow near the surface, these expressions are used to estimate the sensitivity required to perform direct optical measurements and to derive the shear force for indirect optical measurements. The derivation of this result permits the incorporation of better design parameters for other fiber-based sensors. Future work includes analyzing the optical parametric designs of fiber-optic sensors, modeling sensors to utilize the parameters for hydrodynamics and applying different mixtures of hydrodynamic flow. Finally, the fabrication of fiber-optic sensors for hydrodynamic flow applications of the type described in this presentation could enhance aerospace, submarine, and medical technology.

Lyons, Donald R.↗

Flow measurement. II.

Cryogenic fluid mass flow measurement methods, noting flow meter calibration

FLOW MEASUREMENT↗

Simulating transitional flow and heat transfer over the flat plate and circular cylinder using a K-epsilon turbulence model

A boundary layer computer code, called TEXSTAN, has been used to simulate external, transitional, fluid flow and heat transfer over the flat plate and circular cylinder. Turbulence modeling was accomplished by using a low-Reynolds number K-epsilon turbulence model. The study was made for free-stream turbulence intensities less than 15 percent. The comparison of the flat plate heat transfer predictions to experimental data showed good agreement, for the location of transition as well as in the heat transfer coefficient. The typical discrepancy in predicting the heat transfer coefficient was less than 5 percent of the measured value.

Sullivan, T. J.↗

Study of FES/CAST/HGS

The microgravity materials processing program has been instrumental in providing the crystal growth community with an experimental environment to better understand the phenomena associated with the growing of crystals. In many applications one may pursue the growth of large single crystals which cannot be grown on earth due to convective driven flows. A microgravity environment is characterized by neither convection of buoyancy. Consequently superior crystals are able to be grown in space. On the other hand, since neither convection nor buoyancy dominates the fluid flow in a microgravity environment, then lesser dominating phenomena can affect crystal growth, such as surface driven flows or diffusion limited solidification. In the case of experiments that are to be flown in space using the Fluid Experiments System (FES), diffusion limited growth should be the dominating phenomenon. The use of holographic and Schlieren optical techniques for studying the concentration gradients in solidification processes has been used by several investigators over the years. The Holographic Ground System (HGS) facility at MSFC has been a primary resource in researching this capability. Consequently scientific personnel have been able to utilize these techniques in both ground based research and in space experiments. An important event in the scientific utilization of the HGS facilities was the TGS (triglycine sulfate) Crystal Growth and the Casting and Solidification Technology (CAST) experiments that were flown on the International Microgravity Lab (IML) mission in March of this year. The preparation and processing of these space observations are the primary experiments reported in this work. This project provides some ground-based studies to optimize on the holographic techniques used to acquire information about the crystal growth processes flown on IML. Since the ground-based studies will be compared with the space-based experimental results, it is necessary to conduct sufficient ground based studies to best determine how the experiment in space worked. The current capabilities in computer based systems for image processing and numerical computation have certainly assisted in those efforts. As anticipated, this study has certainly shown that these advanced computing capabilities are helpful in the data analysis of such experiments.

Workman, Gary L.↗

Structure and Shape of Surface-Mediated Assembly of Surfactants

Achieving controls on the self-assembly of cationic and block copolymer micelles has significant implications for advancing novel material systems for energy, environmental, and biological applications and for tuning fluid flow behavior in subsurface geologic environments. While it is known that the aggregation of micelles is influenced by the composition of the surrounding fluid, the influence of a solid interface has been less explored. In this study, we probe the organization of cetyltrimethylammonium bromide (CTAB) micelles in the absence and presence of the Pluronic P123 block copolymer and quartz substrate using transmission and grazing-incidence small-angle X-ray scattering measurements and classical molecular dynamics (MD) simulations. In the absence of the quartz interface, CTAB with and without P123 molecules assemble as ellipsoid core–shell micelles. Densification of the core on adding P123 is noted from the decrease in the core radius from 15.9 to 14.6 Å. Here, the presence of a quartz substrate causes the micelles to elongate, which is noted by the emergence of a power-law slope in the low q region (<0.02 Å –1 ). Moreover, in the presence of both P123 and the quartz substrate, the micelle shape changes from ellipsoid core–shell to cylindrical core–shell, and a significantly higher number of aggregates are formed. The higher number of aggregates, and faster aggregation kinetics are linked to the organization of the solvent structures as noted from the MD simulations. These findings demonstrate the effects of adding a block copolymer and including a quartz substrate on the self-assembly of cationic (CTAB) micelles, which can be used to inform the design of functional material systems that harness surfactant chemistries.

42 ENGINEERING↗

Numerical analysis of downhole heat exchanger designed for geothermal energy production

A traditional way of geothermal energy production implies hot brine extraction from the reservoir to a surface facility. This method is economically profitable only for high enthalpy reservoirs associated with shallow depths and high temperatures. The low enthalpy geothermal (LEG) projects are out of consideration due to the high cost of produced electric power related to expensive drilling, cost of the binary cycle system, and long term installation cost. This paper presents a new zero mass extraction method utilizing downhole heat exchanger (DHE) with no geo-fluid production to the surface. The well design stays as a single horizontal well with coupled production and injection sections. The brine pump is located between the sections providing hot brine circulation through the DHE. The coupled fluid flow and heat transfer mathematical model was developed and simulated using nodal analysis method. The LEG reservoir prototype located in South Louisiana was used to study several cases of the DHE lengths, inclination angles, reservoir permeabilities, and flow rates optimization. According to the analysis, the power unit is able to provide 160 kW of net electric power with CO2 working fluid circulating inside a single lateral well. Increasing the reservoir temperature and the number of laterals, the available power production rises up to 600 kW with an attractive electricity cost of $21.84/MWh at 7000 m well depth. This methods opens a new prospective for the depleted petroleum wells converting them to the electricity production units.

15 GEOTHERMAL ENERGY↗

Geomechanical and Fluid Transport Properties

Complex changes to fluid and mechanical properties of a pyrolyzing formation could lead to significant subsidence. Capturing subsidence requires accounting for changes in porosity due to generation and drainage of fluid. Pores, permeability and pore pressure are generated during pyrolysis of kerogen into oil and gas. In a series reaction scheme, oil is an intermediate product of pyrolysis. A large portion of oil-range fluids may be lost in an in-situ operation if the oil produced cannot be drained quickly enough. Petroleum fluid quality and production are highly dependent on the relative rates of creation, destruction, and transport of oil and gas in a given volume. All of these are dependent on the heating rate of oil shale. If oil shale is heated quickly, fluid flow may increase, but at a cost to the amount of oil that is finally extracted. To account for the effect of pore fluid on mechanical deformation, an elastic model was combined with a Mohr-Coulomb failure criterion and Biot’s linear consolidation model. The overall model is poroelastic, because it accounts for rock deforming by considering the response of the bulk system and the pore fluids it contains. Biot’s linear model includes all ways that pore fluid volume can change in a system, including stress-strain interactions between pores and the bulk material. Additionally, the model captures the distinction between fluid and solid matrix compression and describes how bulk deformation is affected by fluid drainage. The fluid response to changing stress is affected by changes in pore pressure, while the bulk material responds to effective stress changes. To capture all of these effects, a number of theoretical and empirical models were employed in FLAC3D™, a rock-mechanics software package produced by Itasca Consulting Group. The basic capabilities of FLAC3D™ have been supplemented by customized algorithms for the simulation package. The most important contribution has been to develop a methodology for describing poro-thermo-mechanical properties of oil shale across the complete chronology of in situ pyrolysis. This report highlights those methods and results for inferring realistic oil shale properties.

04 OIL SHALES AND TAR SANDS↗

Spectral SP: A New Approach to Mapping Reservoir Flow and Permeability

Our objectives for the current project were to develop an innovative inversion and analysis procedure for magnetotelluric field data and time variable self-potentials that will enable us to map not only the subsurface resistivity structure of a geothermal prospect but to also delineate the permeability distribution within the field. Hence, the ultimate objective were to provide better targeting information for exploratory and development drilling of a geothermal prospect. Field data were collected and analyzed from the Kilauea Summit, Kilauea East Rift Zone, and the Humuula Saddle between Mauna Loa and Mauna Kea volcanoes. All of these areas were known or suspected to have geothermal activity of varying intensities. Our results provided evidence for significant long-term coordinated changes in spontaneous potential that could be associated with subsurface flows, significant interferences were encountered that arose from surface environmental changes (rainfall, temperature) that rendered it nearly impossible to unequivocally distinguish between deep fluid flow changes and environmental effects. Further, the analysis of the inferred spontaneous potential changes in the context of depth of the signals, and hence, permeability horizons, were unable to be completed in the time available.

15 GEOTHERMAL ENERGY↗

Turbulence and numerical wind tunnels

The feasibility and possible significance of aerodynamic simulations are considered. Despite certain limitations, there is reason to believe that many useful and important flows containing turbulence can be computed with an accuracy acceptable for use in aircraft designs. Attention is given to the equations for unsteady fluid flow, the direct calculation of unsteady viscous flows, the direct calculation of nonlinear equations, the computation of flows about bodies using turbulence modeling, future developments in numerical methods, and future developments in computers.

Lomax, H.↗

Some flow phenomena in a constant area duct with a Borda type inlet including the critical region

Mass limiting flow characteristics for a 55 L/D tube with a Borda type inlet have been assessed over large ranges of temperature and pressure using fluid nitrogen. Under certain conditions, separation and pressure drop at the inlet was sufficiently strong to permit partial vaporization and the remaining fluid flowed through the tube as if it were a free jet. An empirical relation was determined which defines conditions under which this type of flow can occur. A flow coefficient is presented which enables one to estimate flow rates over the experimental range. A flow rate stagnation pressure map for selected stagnation isotherms and pressure profiles document these flow phenomena.

Hendricks, R. C.↗

Computation of axisymmetric and ionized hypersonic flows using particle and continuum methods

Comparisons between particle and continuum simulations of hypersonic near-continuum flows are presented. The particle approach employs the direct simulation Monte Carlo (DSMC) method, and the continuum approach solves the appropriate equations of fluid flow. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code that is efficiently vectorized is developed for this study. In this DSMC code, particular attention is paid to matching the relaxation rates employed in the continuum approach. This investigation represents a continuum of a previous study that considered thermochemical relaxation in one-dimensional shock waves of nitrogen. Comparison of the particle and continuum methods is first made for an axisymmetric blunt-body flow of air at 7 km/s. Very good agreement is obtained for the two solutions. The two techniques also compare well for a one-dimensional shock wave in air at 10 km/s. In both applications, the results are found to be sensitive to various aspects of the chemistry model employed.

Boyd, Iain D.↗

Flow in Rotating Serpentine Coolant Passages With Skewed Trip Strips

Laser velocimetry was utilized to map the velocity field in serpentine turbine blade cooling passages with skewed trip strips. The measurements were obtained at Reynolds and Rotation numbers of 25,000 and 0.24 to assess the influence of trips, passage curvature and Coriolis force on the flow field. The interaction of the secondary flows induced by skewed trips with the passage rotation produces a swirling vortex and a corner recirculation zone. With trips skewed at +45 deg, the secondary flows remain unaltered as the cross-flow proceeds from the passage to the turn. However, the flow characteristics at these locations differ when trips are skewed at -45 deg. Changes in the flow structure are expected to augment heat transfer, in agreement with the heat transfer measurements of Johnson, et al. The present results show that trips are skewed at -45 deg in the outward flow passage and trips are skewed at +45 deg in the inward flow passage maximize heat transfer. Details of the present measurements were related to the heat transfer measurements of Johnson, et al. to relate fluid flow and heat transfer measurements.

Tse, David G.N.↗

Waves in Radial Gravity Using Magnetic Fluid

Terrestrial laboratory experiments studying various fluid dynamical processes are constrained, by being in an Earth laboratory, to have a gravitational body force which is uniform and unidirectional. Therefore fluid free-surfaces are horizontal and flat. Such free surfaces must have a vertical solid boundary to keep the fluid from spreading horizontally along a gravitational potential surface. In atmospheric, oceanic, or stellar fluid flows that have a horizontal scale of about one-tenth the body radius or larger, sphericity is important in the dynamics. Further, fluids in spherical geometry can cover an entire domain without any sidewall effects, i.e. have truly periodic boundary conditions. We describe spherical body-force laboratory experiments using ferrofluid. Ferrofluids are dilute suspensions of magnetic dipoles, for example magnetite particles of order 10 nm diameter, suspended in a carrier fluid. Ferrofluids are subject to an additional body force in the presence of an applied magnetic field gradient. We use this body force to conduct laboratory experiments in spherical geometry. The present study is a laboratory technique improvement. The apparatus is cylindrically axisymmetric. A cylindrical ceramic magnet is embedded in a smooth, solid, spherical PVC ball. The geopotential field and its gradient, the body force, were made nearly spherical by careful choice of magnet height-to-diameter ratio and magnet size relative to the PVC ball size. Terrestrial gravity is eliminated from the dynamics by immersing the "planet" and its ferrofluid "ocean" in an immiscible silicone oil/freon mixture of the same density. Thus the earth gravity is removed from the dynamics of the ferrofluid/oil interface and the only dynamically active force there is the radial magnetic gravity. The entire apparatus can rotate, and waves are forced on the ferrofluid surface by exterior magnets. The biggest improvement in technique is in the wave visualization. Fluorescing dye is added to the oil/freon mixture and an argon ion laser generates a horizontal light that can be scanned vertically. Viewed from above, the experiment is a black circle with wave deformations surrounded by a light background. A contour of the image intensity at any light sheet position gives the surface of the ferrofluid "ocean" at that "latitude". Radial displacements of the waves as a function of longitude are obtained by subtracting the contour line positions from a no-motion contour at that laser sheet latitude. The experiments are run by traversing the forcing magnet with the laser sheet height fixed and images are frame grabbed to obtain a time-series at one latitude. The experiment is then re-run with another laser-sheet height to generate a full picture of the three-dimensional wave structure in the upper hemisphere of the ball as a function of time. We concentrate here on results of laboratory studies of waves that are important in Earth's atmosphere and especially the ocean. To get oceanic scaling in the laboratory, the experiment must rotate rapidly (4-second rotation period) so that the wave speed is slow compared to the planetary rotation speed as in the ocean. In the Pacific Ocean, eastward propagating Kelvin waves eventually run into the South American coast. Theory predicts that some of the wave energy should scatter into coastal-trapped Kelvin waves that propagate north and south along the coast. Some of this coastal wave energy might then scatter into mid-latitude Rossby waves that propagate back westward. Satellite observations of the Pacific Ocean sea-surface temperature and height seem to show signatures of westward propagating mid-latitude Rossby waves, 5 to 10 years after the 1982-83 El Nino. The observational data is difficult to interpret unambiguously owing to the large range of motions that fill the ocean at shorter timescales. This series of reflections giving eastward, north- ward, and then westward traveling waves is observed cleanly in the laboratory experiments, confirming the theoretical expectations.

Ohlsen, D. R.↗

Adjustable primary and supplemental power units

A propulsion system for an aircraft includes a core flow path in communication with a compressor section, combustor section and a turbine section. A first bottoming cycle system includes a bottoming working fluid flow in thermal communication with a high energy exhaust gas flow that is generated by the core engine. The first bottoming cycle is configured to recover power from the high energy exhaust gas flow in a first engine operating condition and in a second engine operating condition. A second bottoming cycle system is configured to recover power from the high energy exhaust gas flow in the first engine operating condition and not to recover power in a second engine operating condition.

MacDonald, Malcolm P.↗

A high-pressure flow through test vessel for neutron imaging and neutron diffraction-based strain measurement of geological materials

Neutron scattering and neutron imaging have emerged as powerful methods for experimentally investigating material deformation and fluid flow in the interior of otherwise inaccessible or opaque structures. This paper describes the design and provides example uses of a pressure cell developed for investigating such behaviors within geological materials. The cell can accommodate cylindrical samples with diameters up to 38.1 mm and lengths up to 154 mm. Additionally, ports in the cell and a pressure isolating sleeve around the sample allow the independent application of confining pressure up to 69 MPa and axial pressure up to 34.5 MPa. Furthermore, two material versions of the cell have been manufactured and used to date. An aluminum version is typically used for temperatures below 40 °C, because of its relative transparency to neutrons, while a titanium version, which is comparatively more neutron attenuating, is used for experiments requiring triaxial pressurization under conditions up to 350 °C. The pressure cells were commissioned at the VULCAN engineering diffractometer at the Oak Ridge National Laboratory (ORNL), Spallation Neutron Source, and have since been used at the ORNL high flux isotope reactor CG1-D imaging beamline, National Institute of Standard and Technology (NIST) BT-2, and NIST NG6 imaging beamlines.

47 OTHER INSTRUMENTATION↗

Actuator disk theory for incompressible highly rotating flows.

A solution has been obtained for a stator-rotor pair operating in an annulus with constant hub and tip radii. The stator and rotor are represented as actuator discs, and perfect fluid flow is assumed. The solutions are exact within these limitations, no linearization being required. The forms of blade loadings considered allow the introduction of large vorticity by either the rotor or stator. As a result, the rotor may be a ?nonconstant-work' row. The solutions obtained are of summational form, but many of the summations are obtained in closed form, the resultant formulas allowing rapid calculation of desired examples. An example numerical result is included.

Oates, G. C.↗

Preconditioned conjugate-gradient methods for low-speed flow calculations

An investigation is conducted into the viability of using a generalized Conjugate Gradient-like method as an iterative solver to obtain steady-state solutions of very low-speed fluid flow problems. Low-speed flow at Mach 0.1 over a backward-facing step is chosen as a representative test problem. The unsteady form of the two dimensional, compressible Navier-Stokes equations is integrated in time using discrete time-steps. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. The new iterative solver is used to solve a linear system of equations at each step of the time-integration. Preconditioning techniques are used with the new solver to enhance the stability and convergence rate of the solver and are found to be critical to the overall success of the solver. A study of various preconditioners reveals that a preconditioner based on the Lower-Upper Successive Symmetric Over-Relaxation iterative scheme is more efficient than a preconditioner based on Incomplete L-U factorizations of the iteration matrix. The performance of the new preconditioned solver is compared with a conventional Line Gauss-Seidel Relaxation (LGSR) solver. Overall speed-up factors of 28 (in terms of global time-steps required to converge to a steady-state solution) and 20 (in terms of total CPU time on one processor of a CRAY-YMP) are found in favor of the new preconditioned solver, when compared with the LGSR solver.

Ajmani, Kumud↗