Search NASA⌕ Search

SEARCH · Search NASA

Results for “viscous media”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
NASA NTRS1988

A survey of oscillating flow in Stirling engine heat exchangers

Similarity parameters for characterizing the effect of flow oscillation on wall shear stress, viscous dissipation, pressure drop and heat transfer rates are proposed. They are based on physical agruments and are derived by normalizing the governing equations. The literature on oscillating duct flows, regenerator and porous media flows is surveyed. The operating characteristics of the heat exchanger of eleven Stirling engines are discribed in terms of the similarity parameters. Previous experimental and analytical results are discussed in terms of these parameters and used to estimate the nature of the oscillating flow under engine operating conditions. The operating points for many of the modern Stirling engines are in or near the laminar to turbulent transition region. In several engines, working fluid does not pass entirely through heat exchangers during a cycle. Questions that need to be addressed by further research are identified.

Simon, Terrence W.↗
NASA NTRS1985

Autoclave processing for composite material fabrication. 1: An analysis of resin flows and fiber compactions for thin laminate

High quality long fiber reinforced composites, such as those used in aerospace and industrial applications, are commonly processed in autoclaves. An adequate resin flow model for the entire system (laminate/bleeder/breather), which provides a description of the time-dependent laminate consolidation process, is useful in predicting the loss of resin, heat transfer characteristics, fiber volume fraction and part dimension, etc., under a specified set of processing conditions. This could be accomplished by properly analyzing the flow patterns and pressure profiles inside the laminate during processing. A newly formulated resin flow model for composite prepreg lamination process is reported. This model considers viscous resin flows in both directions perpendicular and parallel to the composite plane. In the horizontal direction, a squeezing flow between two nonporous parallel plates is analyzed, while in the vertical direction, a poiseuille type pressure flow through porous media is assumed. Proper force and mass balances have been made and solved for the whole system. The effects of fiber-fiber interactions during lamination are included as well. The unique features of this analysis are: (1) the pressure gradient inside the laminate is assumed to be generated from squeezing action between two adjacent approaching fiber layers, and (2) the behavior of fiber bundles is simulated by a Finitely Extendable Nonlinear Elastic (FENE) spring.

Hou, T. H.↗
DOE OSTI2025

A rheological model for loose sands with insights from DEM

A rheological model for loose granular media is developed to capture both solid-like and fluid-like responses during shearing. The proposed model is built by following the mathematical structure of an extended Kelvin–Voigt model, where an elastic spring and plastic slider act in parallel to a viscous damper. This arrangement requires the partition of the total stress into rate-independent and rate-dependent stress components. To model the solid-like behavior, a simple frictional plasticity model is adopted without modifications, thus contributing to the rate-independent stress. Instead, the fluid-like or rate-dependent stress is further decomposed into deviatoric and volumetric parts, by proposing a new formulation based on a combination of the μ(I) relation, originally developed under pressure-controlled shear, with a pressure-shear rate relation derived under volume-controlled shear. The proposed formulation allows the model to capture both the increase in the friction coefficient and the enhanced dilation at high shear rates. High-fidelity simulation data, obtained from discrete element method and multiscale modelling, are used to evaluate the performance of the proposed constitutive model. The model provides accurate results under both drained and undrained simple shear paths across a wide range of shear rates. Furthermore, it successfully reproduces at much lower computational cost the flowslide mobility computed through multiscale simulations, which is primarily regulated by the shear rate dependence of the material properties during the dynamic runout stage.

Elasticity↗
NASA NTRS2020

Gas-Liquid Flows Through Porous Media in Microgravity: The International Space Station Packed Bed Reactor Experiment

Experimental results on pressure drop and flow patterns for gas-liquid flow through packed beds obtained in the International Space Station with two types of packing are presented and analyzed. It is found that the pressure drop depends on the packing wettability in the viscous-capillary (V-C) regime and this dependence is compared with previously published results developed using short duration low-gravity aircraft tests. Within the V-C regime, the capillary contribution is the dominant force contributing to the pressure drop for the wetting case (glass) versus the viscous contribution dominating for the non-wetting case (Teflon). Outside of the V-C regime, it is also found that hysteresis effects that are often strong in normal gravity gas-liquid flows are greatly diminished in microgravity and pressure drop is nearly independent of packing wettability. A flow pattern transition map from bubble to pulse flow is also compared with the earlier aircraft data.

Viscous-capiliary regime↗
NASA NTRS1975

Propagation of acoustic modes in the transitional ionosphere

Continuum theory is unsatisfactory for describing wave propagation in the transitional ionosphere, since the mean free path is of the same order as the characteristic length (i.e., Knudsen number in this region is of the order of unity). Therefore a transition model in which isotropic electrons behave as a fluid while the ions are governed by the kinetic equation modified by ion-electron collision effects is proposed to study acoustic wave propagation in the upper atmosphere in the altitude region from 500 to 2000 km. The results show that the dissipation of acoustic waves by electrons is through the viscous and thermal conduction effects, and the dissipation of acoustic waves by ions is through Landau damping. Comparisons of the characteristics of acoustic mode propagation and wave dissipation mechanisms in the collision-dominated, transitional from collisional to collisionless, and collisionless media are also discussed.

Hung, R. J.↗
DOE OSTI2020

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1981

A title-gap flow model for use in aerodynamic loads assessment of space shuttle thermal protection system: Parallel gap faces

The problem of predicting aerodynamic loads on the insulating tiles of the space shuttle thermal protection system (TPS) is discussed and seen to require a method for predicting pressure and mass flux in the gaps between tiles. A mathematical model of the tile-gap flow is developed, based upon a slow viscous (Stokes) flow analysis, and is verified against experimental data. The tile-gap pressure field is derived from a solution of the two-dimensional Laplace equation; the mass-flux vector is then calculated from the pressure gradient. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is given. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is given. The flow model shows tile-gap mass flux to be very sensitive to the gap width indicating a need for coupling the TPS flow and tile displacement calculation. Analytical and experimental work to improve TPS flow predictions and a possible shuttle TPS hardware modification are recommended.

Dwoyer, D. L.↗
NASA NTRS1982

Flow through the tile gaps in the Space Shuttle Thermal Protection System

The problem of predicting aerodynamic loads on the insulating tiles of the Space Shuttle Thermal Protection System (TPS) is discussed and seen to require a method for predicting pressure and mass flux in the gaps between tiles. A mathematical model of the tile-gap flow is developed based upon a slow viscous (Stokes) flow analysis and is verified against available experimental data. This model derives the tile-gap pressure field from a solution of the two-dimensional Laplace equation; the mass flux vector is then calculated from the pressure gradient. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is also given. The flow model shows tile-gap mass flux to be very sensitive to the gap width indicating a need for coupling the TPS flow and tile displacement calculations. Finally recommendations are made concerning additional analytical and experimental work to improve TPS flow predictions.

Dwoyer, D. L.↗
NASA NTRSDate not supplied

The Plant Water Management Experiments: Soil

A simple means of watering plants in the low-g environment aboard orbiting spacecraft is not obvious. Since the beginning of spaceflight, numerous approaches have been pursued to water plants that seek to maximize plant viability and system reliability, while minimizing crew time and system complexity. We are not there yet. The Plant Water Management (PWM) Soil experiments seek to apply recent advances in low-g capillary fluidics phenomena to the challenges faced by plant growth operations aboard spacecraft. The primary challenge is to establish earth-like flows minimizing low-g specific adaptations required of the plants. This is difficult due to the ever-present fluid physics challenges of poorly-wetting multiphase inertial-visco-capillary flows in geometrically complex conduits and containers. In this paper, we present recent flight results for the PWM Soil experiments where arcillite ‘soil reservoirs’ are arranged in a non-wetting host soil that serves as an O2-breathing wetting barrier. In this way, a largely terrestrial water-soil environment is mimicked where, as liquid is evapo-transpired through the growing plant foliage, the effective water table passively ‘falls’ reducing viscous lengths and increasing water uptake for the plant. We present data from 6 days of 24-7 experiments on the ISS testing 3 different plant root models. We also present and correlate a capillary flow model which captures the primary features of the flow. Our summary is valued for the assessment of current and future low-g plant watering systems employing soil media.

microgravity↗
NASA NTRS1972

Mathematical characterization of mechanical behavior of porous frictional granular media

A new definition of loading and unloading along the yield surface of Roscoe and Burland is introduced. This is achieved by noting that the strain-hardening parameter in the plastic potential function is deduced from the yield locus equation of Roscoe and Burland. The analytical results are compared with the experimental results for plate-bearing and cone-penetrometer problems and close agreements are demonstrated. The wheel-soil interaction is studied under dynamic loading. The rate-dependent plasticity or viscoelastoplastic behavior is considered. This is accomplished by the internal (hidden) variables associated with time-dependent viscous properties directly superimposed with inelastic behavior governed by the yield criteria of Roscoe and Burland. Effects of inertia and energy dissipation are properly accounted for. Example problems are presented.

Chung, T. J.↗
NASA NTRS2003

Suppression of Buoyancy in Gaseous Media at High Temperatures

Consider a rectangular box filled with a fluid having a heated bottom and a cold top surface, and insulated side-walls (Benard problem). As the temperature difference between the horizontal top and bottom surfaces increases, a critical condition, defined quantitatively by the Rayleigh number, is reached beyond which density stratification can no longer be sustained by conduction and the fluid disrupts from its stable, quiescent state into an unstable, convective mode in which lighter and heavier gas mix. This paper suggests that such a statement is not necessarily true for gaseous media under normalized temperature differences that are much larger than justifiable for the Boussinesq approximation! In fact, there may be situations where a system cannot ever be made unstable with respect to the onset on buoyant convection no matter how large the temperature (density) difference becomes at a given pressure even under normal gravity! This unexpected behavior is primarily attributed to highly temperature-sensitive kinematic viscosity which counteracts the tendency toward instability and dampens convection by making the gas more viscous at higher temperatures. This compensation of the buoyant force by the viscous force exhibits itself by the formation of a peak hot-surface temperature beyond which a system will tend to be more stable as the hot-surface temperature increases.

Gokoglu, Suleyman A.↗
DOE OSTI2023

Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate

Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

42 ENGINEERING↗
NASA NTRS2010

Multiphysics Modeling for Dimensional Analysis of a Self-Heated Molten Regolith Electrolysis Reactor for Oxygen and Metals Production on the Moon and Mars

The technology of direct electrolysis of molten lunar regolith to produce oxygen and molten metal alloys has progressed greatly in the last few years. The development of long-lasting inert anodes and cathode designs as well as techniques for the removal of molten products from the reactor has been demonstrated. The containment of chemically aggressive oxide and metal melts is very difficult at the operating temperatures ca 1600 C. Containing the molten oxides in a regolith shell can solve this technical issue and can be achieved by designing a self-heating reactor in which the electrolytic currents generate enough Joule heat to create a molten bath. In a first phase, a thermal analysis model was built to study the formation of a melt of lunar basaltic regolith irradiated by a focused solar beam This mode of heating was selected because it relies on radiative heat transfer, which is the dominant mode of transfer of energy in melts at 1600 C. Knowing and setting the Gaussian-type heat flux from the concentrated solar beam and the phase and temperature dependent thermal properties, the model predicts the dimensions and temperature profile of the melt. A validation of the model is presented in this paper through the experimental formation of a spherical cap melt realized by others. The Orbitec/PSI experimental setup uses an 3.6-cm diameter concentrated solar beam to create a hemispheric melt in a bed of lunar regolith simulant contained in a large pot. Upon cooling, the dimensions of the vitrified melt are measured to validate the thermal model. In a second phase, the model is augmented by multiphysics components to compute the passage of electrical currents between electrodes inserted in the molten regolith. The current through the melt generates Joule heating due to the high resistivity of the medium and this energy is transferred into the melt by conduction, convection and primarily by radiation. The model faces challenges in two major areas, the change of phase as temperature increases, and the dominance of radiative heat flux as heat transfer mechanism within the melt the change of phase concerns the regolith itself which is present in states ranging from a fine grain regolith with low thermal conductivity and low density to a vitrified melt with much higher thermal conductivity, and higher density. As the regolith is heated, it starts to soften around 1300 C the melt iS very viscous and evolving gas bubbles out in thick, lava-like fashion. By 1600 C the regolith is completely melted and the viscosity is low The second challenge resides in the proper modeling of the radiative heat flux requiring the addition of the computing-demanding radiative-heat-transfer function to the general heat transfer equation. The model Includes temperature-dependent properties (density, thermal conductivity, heat capacity, and viscosity, and absorption coefficients) and solves the radiative heat flux equation assuming gray (fine grains) and semi-transparent (melt) media and using an absorption coefficient spectral found in the literature for terrestrial minerals similar in composition to those of lunar regolith simulant

Dominguez, Jesus A.↗
DOE OSTI2021

Optimal Wetting Angles in Lattice Boltzmann Simulations of Viscous Fingering

We conduct pore-scale simulations of two-phase flow using the 2D Rothman–Keller colour gradient lattice Boltzmann method to study the effect of wettability on saturation at breakthrough (sweep) when the injected fluid first passes through the right boundary of the model. We performed a suite of 189 simulations in which a “red” fluid is injected at the left side of a 2D porous model that is initially saturated with a “blue” fluid spanning viscosity ratios M=ν r /ν b ∈[0.001,100] and wetting angles θ w ∈[0°,180°]. As expected, at low-viscosity ratios M=ν r /ν b $\ll$1 we observe viscous fingering in which narrow tendrils of the red fluid span the model, and for high-viscosity ratios M$\gg$1, we observe stable displacement. The viscous finger morphology is affected by the wetting angle with a tendency for more rounded fingers when the injected fluid is wetting. However, rather than the expected result of increased saturation with increasing wettability, we observe a complex saturation landscape at breakthrough as a function of viscosity ratio and wetting angle that contains hills and valleys with specific wetting angles at given viscosity ratios that maximize sweep. This unexpected result that sweep does not necessarily increase with wettability has major implications to enhanced oil recovery and suggests that the dynamics of multiphase flow in porous media has a complex relationship with the geometry of the medium and the hydrodynamical parameters.

Engineering↗
NASA NTRS2002

RWPV bioreactor mass transport: earth-based and in microgravity

Mass transport and mixing of perfused scalar quantities in the NASA Rotating Wall Perfused Vessel bioreactor are studied using numerical models of the flow field and scalar concentration field. Operating conditions typical of both microgravity and ground-based cell cultures are studied to determine the expected vessel performance for both flight and ground-based control experiments. Results are presented for the transport of oxygen with cell densities and consumption rates typical of colon cancer cells cultured in the RWPV. The transport and mixing characteristics are first investigated with a step change in the perfusion inlet concentration by computing the time histories of the time to exceed 10% inlet concentration. The effects of a uniform cell utilization rate are then investigated with time histories of the outlet concentration, volume average concentration, and volume fraction starved. It is found that the operating conditions used in microgravity produce results that are quite different then those for ground-based conditions. Mixing times for microgravity conditions are significantly shorter than those for ground-based operation. Increasing the differential rotation rates (microgravity) increases the mixing and transport, while increasing the mean rotation rate (ground-based) suppresses both. Increasing perfusion rates enhances mass transport for both microgravity and ground-based cases, however, for the present range of operating conditions, above 5-10 cc/min there are diminishing returns as much of the inlet fluid is transported directly to the perfusion exit. The results show that exit concentration is not a good indicator of the concentration distributions in the vessel. In microgravity conditions, the NASA RWPV bioreactor with the viscous pump has been shown to provide an environment that is well mixed. Even when operated near the theoretical minimum perfusion rates, only a small fraction of the volume provides less than the required oxygen levels. 2002 Wiley Periodicals, Inc.

Models, Biological↗
DOE OSTI2024

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗
NASA NTRS2006

Pattern Formation in Diffusion Flames Embedded in von Karman Swirling Flows

Pattern formation is observed in nature in many so-called excitable systems that can support wave propagation. It is well-known in the field of combustion that premixed flames can exhibit patterns through differential diffusion mechanism between heat and mass. However, in the case of diffusion flames where fuel and oxidizer are separated initially there have been only a few observations of pattern formation. It is generally perceived that since diffusion flames do not possess an inherent propagation speed they are static and do not form patterns. But in diffusion flames close to their extinction local quenching can occur and produce flame edges which can propagate along stoichiometric surfaces. Recently, we reported experimental observations of rotating spiral flame edges during near-limit combustion of a downward-facing polymethylmethacrylate disk spinning in quiescent air. These spiral flames, though short-lived, exhibited many similarities to patterns commonly found in quiescent excitable media including compound tip meandering motion. Flame disks that grow or shrink with time depending on the rotational speed and in-depth heat loss history of the fuel disk have also been reported. One of the limitations of studying flame patterns with solid fuels is that steady-state conditions cannot be achieved in air at normal atmospheric pressure for experimentally reasonable fuel thickness. As a means to reproduce the flame patterns observed earlier with solid fuels, but under steady-state conditions, we have designed and built a rotating, porous-disk burner through which gaseous fuels can be injected and burned as diffusion flames. The rotating porous disk generates a flow of air toward the disk by a viscous pumping action, generating what is called the von K rm n boundary layer which is of constant thickness over the entire burner disk. In this note we present a map of the various dynamic flame patterns observed during the combustion of methane in air as a function of fuel flow rate and the burner rotational speed.

Nayagam, Vedha↗