Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Properties”

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 451 records · Page 25

Computational Fluid Dynamics Simulation of Compressible Non-Newtonian Biomass in a Compression-Screw Feeder

Compression-screw feeders play a critical role in biorefineries to transport lignocellulosic biomass feedstocks from the feeding hoppers to biomass-conversion reactors in order to pretreat and convert biomass to hydrocarbon liquid biofuels and other power and energy resources. One of the main challenges in the operation of screw feeder is plugging and jamming of compressed biomass with high concentration of insoluble solids. The focus of this paper is to numerically investigate the screw feeder at these challenging operating conditions and help with the optimization of the screw feeder design to avoid operation failure. In this work a customized CFD model based on open-source OpenFOAM package [1] was developed to simulate the concentrated biomass as a highly viscous non-Newtonian fluid in the screw feeder. The biomass is modeled as a single-phase compressible Bingham fluid with a plastic viscosity as well as a density-dependent yield stress. The compressibility formulation (pressure-dependent density) and the density-dependent yield stress formulation in the governing equations follow the suggestions from a recent study by Duncan et al. [2]. A pilot-scale hopper/screw feeding system at NREL [3] is used to compare the experimental observations with our simulation results. The auger is 280 mm long and tapered with outer diameter changing from 80 mm to 35 mm. The auger rotates from 10 to 60 rpm in a conical throat which contains anti-rotational bars. The simulations predicted the required torque for the screw feeder and the pressure increase at the exit for biomass feedstocks with various fluid viscosity properties and auger rotating speeds. The analysis of the stress forces helped to identify the critical conditions were the screw feeder excessive wear or jamming could occur.

biofuels↗

Geothermal well testing pressure prediction by using a hybrid transformer model system: FORGE well use case

Geothermal has huge potential to become an indispensable component in achieving the goal of sustainable energy economy, given its capability to provide consistent baseload power to the electric grid. Injection tests are crucial in geothermal energy system as they naturally help to evaluate reservoir properties, understand fluid flow and even enhance reservoir performance. In this research, we developed a hybrid model system that integrates machine learning (ML) regression, a physics-based mathematical model, and transformer deep learning. Trained and validated using FORGE injection test dataset, this system can forecast the pressure variations both upward and downward over time. The pressure prediction achieved prediction accuracy within 3-6% variance of true pressure values. The system can significantly save time and reduce costs by testing only a few cycles and then using model predictions for further analysis, instead of conducting additional real injection cycle tests. The developed model system also holds promise for designing injection test processes and maintaining well production in geothermal energy. Presented at the IMAGE ‘25 Conference led by Shell.

FORGE↗

Thermodynamic and transport properties of cryogenic propellants and related fluids

Significant advances have been made in recent years in the quality and range of thermophysical data for the cryogenic propellants, pressurants, and inertants. A review of recently completed and current data compilation projects for helium, hydrogen, argon, nitrogen, oxygen, fluorine, and methane is given together with recommended references for thermodynamic and transport property data tables for these fluids. Modern techniques in the plotting of thermodynamic charts from tabular data (or from functions such as the equation of state) have greatly improved their precision and value. A list of such charts is included.

Johnson, V. J.↗

Water Hammer Simulations of Monomethylhydrazine Propellant

Fluid Transient analysis is important for the design of spacecraft propulsion system to ensure structural stability of the system in the event of sudden closing or opening of the valve. Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code developed at NASA/MSFC is capable of simulating pressure surge due to sudden opening or closing of valve when thermodynamic properties of real fluid are available for the entire range of simulation. Specifically GFSSP needs an accurate representation of pressure density relationship in order to predict pressure surge during a fluid transient. Unfortunately, the available thermodynamic property programs such as REFPROP, GASP or GASPAK do not provide the thermodynamic properties of Monomethylhydrazine(MMH). This work illustrates the process used for building a customized table of properties of state variables from available properties and speed of sound that is required by GFSSP for simulation. Good agreement was found between the simulations and measured data. This method can be adopted for modeling flow networks and systems with other fluids whose properties are not known in detail in order to obtain general technical insight.

Burkhardt, Zachary↗

Water Hammer Simulations of MMH Propellant - New Capability Demonstration of the Generalized Fluid Flow Simulation Program

Fluid Transient analysis is important for the design of spacecraft propulsion system to ensure structural stability of the system in the event of sudden closing or opening of the valve. Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code developed at NASA/MSFC is capable of simulating pressure surge due to sudden opening or closing of valve when thermodynamic properties of real fluid are available for the entire range of simulation. Specifically GFSSP needs an accurate representation of pressure-density relationship in order to predict pressure surge during a fluid transient. Unfortunately, the available thermodynamic property programs such as REFPROP, GASP or GASPAK does not provide the thermodynamic properties of Monomethylhydrazine (MMH). This paper will illustrate the process used for building a customized table of properties of state variables from available properties and speed of sound that is required by GFSSP for simulation. Good agreement was found between the simulations and measured data. This method can be adopted for modeling flow networks and systems with other fluids whose properties are not known in detail in order to obtain general technical insight. Rigorous code validation of this approach will be done and reported at a future date.

Burkhardt, Z.↗

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, Joule-Thompson device, control valve, heat exchanger core, parallel tube, and compressible orifice. The program has the provision of including additional resistance options through User Subroutines. GFSSP employs a finite volume formulation of mass, momentum, and energy conservation equations in conjunction with the thermodynamic equations of state for real fluids as well as energy conservation equations for the solid. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. The application and verification of the code has been demonstrated through 30 example problems.

Majumdar, A. K.↗

Critical Transport Phenomena in Fluid Helium Under Low Gravity

The feasibility of carrying out measurements of certain critical transport properties of pure fluid under conditions of low gravity was studied. These properties are the thermal conductivity, kappa, the shear viscosity zeta and the diffusive relaxation time tau, which are predicted to diverge (tend to infinity) as the liquid-vapor critical point is approached. However, in this critical region, the Earth's gravity effect becomes very important. As the critical point is approached, the gravity effects increasingly distort the results. The reason for this is that the compressibility of the fluid also diverges and under the influence of gravity causes a vertical density gradient in the fluid, which is significant even when very thin fluid layers (typically 1 mm high) are being used. The result is that the temperature dependence of kappa, zeta, and tau tends to flatten off as T sub c is approached instead of continuing to increase, and therefore the predictions from the renormalization group and mode coupling theories cannot be subjected to a satisfactory test.

Meyer, H.↗

Dehydrating and Sterilizing Wastes Using Supercritical CO2

A relatively low-temperature process for dehydrating and sterilizing biohazardous wastes in an enclosed life-support system exploits (1) the superior mass-transport properties of supercritical fluids in general and (2) the demonstrated sterilizing property of supercritical CO2 in particular. The wastes to be treated are placed in a chamber. Liquid CO2, drawn from storage at a pressure of 850 psi (approx.=5.9 MPa) and temperature of 0 C, is compressed to pressure of 2 kpsi (approx.=14 MPa) and made to flow into the chamber. The compression raises the temperature to 10 C. The chamber and its contents are then further heated to 40 C, putting the CO2 into a supercritical state, in which it kills microorganisms in the chamber. Carrying dissolved water, the CO2 leaves the chamber through a back-pressure regulator, through which it is expanded back to the storage pressure. The expanded CO2 is refrigerated to extract the dissolved water as ice, and is then returned to the storage tank at 0 C

Brown, Ian J.↗

Heterogeneous Cluster Energetics and Nonlinear Thermodynamic Response in Supercritical Fluids

Microstructural heterogeneities arising from molecular clusters directly affect the nonlinear thermodynamic properties of supercritical fluids. We present a physical model to elucidate the relation between energy exchange and heterogeneous cluster dynamics during the transition from liquidlike to gaslike conditions. By analyzing molecular-dynamics data and employing physical principles, the model considers contributions from three key processes, namely, changing cluster density, cluster separation, and transfer of molecules between clusters. We show that the proposed model is consistent with the energetics at subcritical conditions and can be used to explain the nonlinear behavior of thermodynamic response functions, including the peak in the isobaric heat capacity.

42 ENGINEERING↗

Using onset times from frequent geophysical surveys to understand reservoir fluid flow and to estimate flow properties

Seismic characterization of flow properties is a difficult prospect, due to the indirect relationship between permeability and seismic velocity and attenuation. At best, seismic time-lapse changes can detect the effects of saturation and pressure changes in a reservoir due to fluid flow. Even in this case, the interpretation of the seismic observation in terms of the state of the reservoir depends intimately on the properties of the rock physics model, which are usually poorly known. The onset time, the calendar time of geophysical changes, provides an alternative datum for characterizing properties such as reservoir permeability. The main advantage of an onset time is that it is sensitive to the flow properties of the reservoir yet insensitive to the details of the rock physics model. Two examples of the utility of onset times are discussed: The use of travel time shifts induced by the injection of carbon dioxide between two wells, and the time-lapse time shifts for elastic waves propagating through a reservoir undergoing enhanced oil recovery. In both examples, the onset times are mapped into permeability estimates using a trajectory-based approach akin to seismic tomography.

58 GEOSCIENCES↗

Blending and cavern formation within non-Newtonian fluids in stirred tanks: Application to nuclear waste fluid processing

A numerical approach for predicting the time-accurate fluid flow and mixing properties of non-Newtonian fluids is presented. This approach, which is based on the lattice Boltzmann method, is used to characterize the blending in shear-thinning fluids, as well as cavern formation in yield-stress fluids. Predictions compare favorably to measured data and expectations from first-principles theory. Importantly, across the range of fluids and systems here, the simulations require no reparameterization or retuning between scenarios. This generality is exploited to model a nuclear waste processing unit operation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Phase change in liquid face seals. II - Isothermal and adiabatic bounds with real fluids

Analytical studies of phase change effects in parallel and tapered liquid face seals are presented. An isothermal and adiabatic model of low Reynolds number flow are considered by numerical integration of the descriptive equations for a real fluid, and its thermodynamic properties are calculated for each step, using a computer program for the steam tables or fluid thermodynamic properties. It was shown that for low leakage rate the isothermal model is more accurate and for high leakage rates the adiabatic model is more accurate; that both models yield the same conclusions regarding stability; and that the transient of collapse is described by the adiabatic model which predicts a catastrophic collapse and then either failure or explosive return to a larger film thickness value. Finally, it is shown that converging seals may become unstable and the mass leakage rate is reduced significantly below the all liquid value when boiling occurs.

Hughes, W. F.↗

Internal Flow Thermal/Fluid Modeling of STS-107 Port Wing in Support of the Columbia Accident Investigation Board

As part of the aero-thermodynamics team supporting the Columbia Accident Investigation Board (CAB), the Marshall Space Flight Center was asked to perform engineering analyses of internal flows in the port wing. The aero-thermodynamics team was split into internal flow and external flow teams with the support being divided between shorter timeframe engineering methods and more complex computational fluid dynamics. In order to gain a rough order of magnitude type of knowledge of the internal flow in the port wing for various breach locations and sizes (as theorized by the CAB to have caused the Columbia re-entry failure), a bulk venting model was required to input boundary flow rates and pressures to the computational fluid dynamics (CFD) analyses. This paper summarizes the modeling that was done by MSFC in Thermal Desktop. A venting model of the entire Orbiter was constructed in FloCAD based on Rockwell International s flight substantiation analyses and the STS-107 reentry trajectory. Chemical equilibrium air thermodynamic properties were generated for SINDA/FLUINT s fluid property routines from a code provided by Langley Research Center. In parallel, a simplified thermal mathematical model of the port wing, including the Thermal Protection System (TPS), was based on more detailed Shuttle re-entry modeling previously done by the Dryden Flight Research Center. Once the venting model was coupled with the thermal model of the wing structure with chemical equilibrium air properties, various breach scenarios were assessed in support of the aero-thermodynamics team. The construction of the coupled model and results are presented herein.

Sharp, John R.↗

CFD-DEM Modeling of Autothermal Pyrolysis of Corn Stover with a Coupled Particle- and Reactor-Scale Framework

Autothermal operation of fast pyrolysis is an efficient process-intensification technique wherein exothermic oxidation reactions are used to overcome the heat-transfer bottleneck of conventional pyrolysis. The development of accurate, reliable modeling toolsets is imperative to generating a deeper understanding of biomass autothermal pyrolysis systems to support scale-up and industrial deployment. This modeling effort describes the development of single-particle and reactor models which incorporate detailed reaction schemes and simultaneous exothermic oxidation reactions. The particle-scale model was parameterized for corn stover feedstock with particle morphology, density, ash content, and biopolymer composition, all of which impact the emergent conversion characteristics during pyrolysis. Results were then used to parameterize a reactor-scale autothermal pyrolysis model, which was developed using a coarse-grained computational fluid dynamic-discrete element method. The simulation results compared well with experimental results, with the predicted bio-oil, light gas, and biochar yield within 3.0 wt% of the experimental yields. Further analyses were performed to test the influence of equivalence ratio, biomass injection position, and particle size distribution on autothermal pyrolysis. The analysis of the physio-chemical properties of the fluid and solid phase inside the reactor and at the reactor outlet help reveal important process interactions of autothermal pyrolysis.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

Rheology of synovial fluid and its role in joint lubrication

A review of studies by Barnett, Davies, MacConaill, Charnley, Ogston, Sundblad, and others indicates that there exists a tendency to explain the lubricating mechanism existing in synovial joints by either the lubricating properties of synovial fluid or the lubricating properties of the articular surfaces. The proponents of the respective hypothesis attempted to apply either the theory of hydrodynamic lubrication or the theory of boundary lubrication. Neither of these theories, however, can represent the lubrication mechanism in joints; both theories are too crude for this purpose. The nearest approach to reality appears to be represented by more recent elastohydrodynamic theories. It is the aim of this paper to show that both the synovial fluid and the articular cartilage play an important role in the joint lubrication, and that thixotropy of synovial fluid is the key parameter in the lubricating performance of the synovial fluid itself, and of the cartilage whose pores it fills. It will be suggested also that any pathological changes in either the synovial fluid or the articular cartilage, lead to a breakdown of joint lubrication.

LUBRICATION↗

Rotating-Pump Design Code

Pump Design (PUMPDES) is a computer program for designing a rotating pump for liquid hydrogen, liquid oxygen, liquid nitrogen, water, methane, or ethane. Using realistic properties of these fluids provided by another program called GASPAK, this code performs a station-by-station, mean-line analysis along the pump flow path, obtaining thermodynamic properties of the pumped fluid at each station and evaluating hydraulic losses along the flow path. The variables at each station are obtained under constraints that are consistent with the underlying physical principles. The code evaluates the performance of each stage and the overall pump. In addition, by judiciously choosing the givens and the unknowns, the code can perform a geometric inverse design function: that is, it can compute a pump geometry that yields a closest approximation of given design point. The code contains two major parts: one for an axial-rotor/inducer and one for a multistage centrifugal pump. The inducer and the centrifugal pump are functionally integrated. The code can be used in designing and/or evaluating the inducer/centrifugal-pump combination or the centrifugal pump alone. The code is written in standard Fortran 77.

Walker, James F.↗

Convective Instability in Ice I with Non-Newtonian Rheology: Application to the Galilean Satellites

At the temperatures and stresses associated with the onset of convection in an ice I shell of the Galilean satellites, ice behaves as a non-Newtonian fluid with a viscosity that depends on both temperature and strain rate. The convective stability of a non-Newtonian ice shell can be judged by comparing the Rayleigh number of the shell to a critical value. Previous studies suggest that the critical Rayleigh number for a non-Newtonian fluid depends on the initial conditions in the fluid layer, in addition to the thermal, rheological, and physical properties of the fluid. We seek to extend the existing definition of the critical Rayleigh number for a non-Newtonian, basally heated fluid by quantifying the conditions required to initiate convection in an ice I layer initially in conductive equilibrium. We find that the critical Rayleigh number for the onset of convection in ice I varies as a power (-0.6 to -0.5) of the amplitude of the initial temperature perturbation issued to the layer, when the amplitude of perturbation is less than the rheological temperature scale. For larger-amplitude perturbations, the critical Rayleigh number achieves a constant value. We characterize the critical Rayleigh number as a function of surface temperature of the satellite, melting temperature of ice, and rheological parameters so that our results may be extrapolated for use with other rheologies and for a generic large icy satellite. The values of critical Rayleigh number imply that triggering convection from a conductive equilibrium in a pure ice shell less than 100 km thick in Europa, Ganymede, or Callisto requires a large, localized temperature perturbation of a few kelvins to tens of kelvins to soften the ice and therefore may require tidal dissipation in the ice shell.

Barr, A. C.↗

Time-Lapse Integration at FWU: Fluids, Rock Physics, Numerical Model Integration, and Field Data Comparison

We present the current status of time-lapse seismic integration at the Farnsworth (FWU) CO2 WAG (water-alternating-gas) EOR (Enhanced Oil Recovery) project at Ochiltree County, northwest Texas. As a potential carbon sequestration mechanism, CO2 WAG projects will be subject to some degree of monitoring and verification, either as a regulatory requirement or to qualify for economic incentives. In order to evaluate the viability of time-lapse seismic as a monitoring method the Southwest Partnership (SWP) has conducted time-lapse seismic monitoring at FWU using the 3D Vertical Seismic Profiling (VSP) method. The efficacy of seismic time-lapse depends on a number of key factors, which vary widely from one application to another. Most important among these are the thermophysical properties of the original fluid in place and the displacing fluid, followed by the petrophysical properties of the rock matrix, which together determine the effective elastic properties of the rock fluid system. We present systematic analysis of fluid thermodynamics and resulting thermophysical properties, petrophysics and rock frame elastic properties, and elastic property modeling through fluid substitution using data collected at FWU. These analyses will be framed in realistic scenarios presented by the FWU CO2 WAG development. The resulting fluid/rock physics models will be applied to output from the calibrated FWU compositional reservoir simulation model to forward model the time-lapse seismic response. Modeled results are compared with field time-lapse seismic measurements and strategies for numerical model feedback/update are discussed. While mechanical effects are neglected in the work presented here, complementary parallel studies are underway in which laboratory measurements are introduced to introduce stress dependence of matrix elastic moduli.

4D↗