Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Porous 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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Scaling of two-phase water-steam relative permeability and thermal fluxes in porous media

Two-phase water-steam flow conditions are frequently encountered in many engineering applications, including geothermal reservoirs. Although routine calculations are based on the multiphase Darcy’s law, the role of the topology of the flowing phases at the pore-scale is usually neglected in the estimation of relative permeabilities. Instead, the latter are frequently computed using empirical models like the Corey correlation. In this work, we first apply the model for relative permeabilities based on pore-scale flow regimes developed by Picchi and Battiato (2019), Relative permeability scaling from pore-scale flow regimes, Water Resour. Res. 55, 3215–3233, to scenarios typical of geothermal reservoirs and then extend it by deriving the scaling laws for the transmissibilities and the thermal properties as a function of temperature. First, we discuss the scaling behavior of normalized relative permeabilities in terms of viscosity ratio and capillary number of water-steam systems and, then, we provide a validation of the model against experimental data available in the literature. The model captures well the data trends collected in real 3D porous media. These results suggest that water-steam relative permeabilities follow the same scaling behavior of gas-liquid systems where the non-wetting phase is much less viscous than the wetting phase. Finally, we investigate the impact that relative permeabilities have on heat transfer rates at two-phase flow conditions and the scaling of mass and energy transmissibility and thermal properties of the mixture. An estimation of the exergy carried by the two-phase water-steam mixture is also included.

42 ENGINEERING↗
DOE OSTI2023

Heat Transfer Characteristics of Particle and Air Flow Through Additively Manufactured Lattice Frame Material Based on Octet-Shape Topology

Particle-to-supercritical carbon dioxide (sCO 2 ) heat exchanger is a critical component in next-generation concentrating solar power (CSP) plants. The inherently low heat transfer between falling particles and sCO 2 imposes a challenge toward economic justification of levelized cost of electricity produced through solar energy. Introduction of integrated porous media with the walls bounding particle flow has the potential to enhance the overall particle-to-sCO 2 heat exchanger performance. Here, this paper presents an experimental study on heat transfer characterization of additively manufactured lattice frame material based on Octet-shaped unit cell with particles and air as working fluids. The lattice structures were additively manufactured in stainless steel (SS) 316L and SS420 (with 40% bronze infiltration) via Binder jetting process, where the lattice porosities were varied between 0.75 and 0.9. The mean particle diameters were varied from 266 μm to 966 μm. The effective thermal conductivity and averaged heat transfer coefficient were determined through steady-state experiments. It was found that the presence of lattice enhances the effective thermal conductivity by 2–4 times when compared to packed bed of particles alone. Furthermore, for gravity-assisted particle flow through lattice panel, significantly high convective heat transfer coefficients ranging from 200 W/m2K to 400 W/m2K were obtained for the range of particle diameters tested. The superior thermal transport properties of Octet-shape-based lattice frame for particle flow makes it a very promising candidate for particle-to-sCO 2 heat exchanger for CSP application.

14 SOLAR ENERGY↗
DOE OSTI2024

Porous Iron Electrodes Reduce Energy Consumption During Electrocoagulation of a Virus Surrogate: Insights into Performance Enhancements Using Three-Dimensional Neutron Computed Tomography

Electrocoagulation has attracted significant attention as an alternative to conventional chemical coagulation because it is capable of removing a wide range of contaminants and has several potential advantages. In contrast to most electrocoagulation research that has been performed with nonporous electrodes, in this study, we demonstrate energy-efficient iron electrocoagulation using porous electrodes. In batch operation, investigation of the external pore structures through optical microscopy suggested that a low porosity electrode with sparse connection between pores may lead to mechanical failure of the pore network during electrolysis, whereas a high porosity electrode is vulnerable to pore clogging. Electrodes with intermediate porosity, instead, only suffered a moderate surface deposition, leading to electrical energy savings of 21% and 36% in terms of electrocoagulant delivery and unit log virus reduction, respectively. Neutron computed tomography revealed the critical role of electrode porosity in utilizing the electrode’s internal surface for electrodissolution and effective delivery of electrocoagulant to the bulk. Energy savings of up to 88% in short-term operation were obtained with porous electrodes in a continuous flow-through system. Further investigation on the impact of current density and porosity in long-term operation is desired as well as the capital cost of porous electrodes.

42 ENGINEERING↗
NASA NTRS1995

HPTAM, a two-dimensional Heat Pipe Transient Analysis Model, including the startup from a frozen state

A two-dimensional Heat Pipe Transient Analysis Model, 'HPTAM,' was developed to simulate the transient operation of fully-thawed heat pipes and the startup of heat pipes from a frozen state. The model incorporates: (a) sublimation and resolidification of working fluid; (b) melting and freezing of the working fluid in the porous wick; (c) evaporation of thawed working fluid and condensation as a thin liquid film on a frozen substrate; (d) free-molecule, transition, and continuum vapor flow regimes, using the Dusty Gas Model; (e) liquid flow and heat transfer in the porous wick; and (f) thermal and hydrodynamic couplings of phases at their respective interfaces. HPTAM predicts the radius of curvature of the liquid meniscus at the liquid-vapor interface and the radial location of the working fluid level (liquid or solid) in the wick. It also includes the transverse momentum jump condition (capillary relationship of Pascal) at the liquid-vapor interface and geometrically relates the radius of curvature of the liquid meniscus to the volume fraction of vapor in the wick. The present model predicts the capillary limit and partial liquid recess (dryout) in the evaporator wick, and incorporates a liquid pooling submodel, which simulates accumulation of the excess liquid in the vapor core at the condenser end.

Tournier, Jean-Michel↗
NASA NTRS1991

Surface reaction and pore diffusion in flow-tube reactors

The interaction of gas diffusion with surface reaction in porous solids is discussed and applied specifically to heterogeneous rate measurements in flow-tube reactors. External diffusion to the outer surface of a reactive solid, internal diffusion within the pores, surface reaction, and laminar flow are considered. A procedure is developed to correct observed surface rate constants for the interaction of these processes. Measured surface areas and bulk densities are used to construct a semiempirical model for porous diffusion in vapor-formed HNO3-H2O ices which are used to simulate polar stratospheric cloud surfaces. The model is tested experimentally by varying the thickness of these ices from about 15 to 120 microns. The results are consistent with the model predictions and show that the HNO3-H2O ices used are highly porous, and the internal surface must be considered in calculating kinetic parameters from observed loss rates. The best fit of the data yields a tortuosity factor of 3.3 +/-1.1 for the ice substrates.

Keyser, Leon F.↗
DOE OSTI2024

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗
NASA NTRS1975

Development, fabrication and calibration of a porous surface microphone in an aerofoil

The development of a porous surface microphone in an airfoil intended to measure acoustic signals in a turbulent airflow and to minimize the flow noise is described. The sensor because of its airfoil operates over a wide range of yaw angles and flow velocities without excessive flow noise. The acoustic properties of the porous materials used in the airfoil sensor and their effects on the frequency response of the sensor were analyzed and tested. An accurate airfoil was selected, having a smaller thickness-to-chord ratio and an airfoil sensor was designed. The sensor was calibrated acoustically and its flow noise evaluated in the quiet BBN wind tunnel at flow velocities up to 70 m/sec. Results are presented.

Noiseux, D. U.↗
DOE OSTI2025

An experimental study on gas-liquid phase fluid migration in hydrate-bearing sediments during hydrate dissociation

Natural gas hydrate production tests face problems such as severe sand blockage, poor gas-liquid phase separation, and significant land subsidence. This is because of the insufficient understanding of the complex phase transition and gas-liquid multi-phase fluid migration during hydrate dissociation. In hydrate-bearing sediment systems, hydrate phase transition couples with gas-liquid fluid migration. The phase transition causes changes in pore structure, which in turn modifies porous infiltration parameters and fluid flow capacity. Meanwhile, alterations in phase interfaces affect key parameters like surface tension and wettability. Gas-liquid fluid migration influences heat and mass transfer, thus affecting phase equilibrium and dissociation rates. To bridge the gap in describing gas-liquid fluid migration during hydrate dissociation in experiments, this research innovatively integrated an unsteady-state gas displacement by water and a quantitative hydrate dissociation process, independently developed a multi-phase seepage experimental system suitable for hydrate dissociation and determined the relationship between seepage parameters and hydrate saturation under different porosity. The results are as follows: a) Core samples with higher initial porosity show a greater recovery rate of fluid flow capacity. b) The retarding effect of multi-phase fluid has a more significant impact on the migration of the wetting phase fluid (water) than that of the non-wetting phase fluid (methane). c) During hydrate dissociation, the evolution of absolute permeability shows an “S-shaped” pattern, and the evolution of relative permeability shows a “wiring-harness” pattern. In conclusion, the findings can provide a theoretical basis for preventing geological disasters and for geotechnical engineering design during hydrate production.

58 GEOSCIENCES↗
DOE OSTI2024

Macroscopic modeling of gas permeability in hierarchical micro/nanoporous media: A unified characterization of rarefaction using Klinkenberg theory and equivalent diameter

Estimating gas transport through a hierarchical micro/nanoporous system is challenging due to non-equilibrium gas dynamics. The primary difficulty lies in determining the rarefaction level, because identifying a representative flow dimension in a complex porous system with multiple pore scales is not straightforward. Our study performed a pore-level analysis for gas permeability in dual-scale porous media with varying porosity, throat size, and secondary pore size under different rarefaction conditions. We found that secondary porosity negatively affects permeability due to increased friction forces, with this influence growing as the secondary pore size and porosity increase until the secondary pore becomes comparable to the throat. However, rarefaction reduces the effects of secondary pores due to boundary slip. Traditional Knudsen number (Kn) calculations based on Darcy-defined height failed to accurately describe the rarefaction effects on gas permeability. Instead, we introduced an equivalent diameter to calculate the Kn, which provided an accurate normalization of apparent gas permeability independent of pore geometry. Furthermore, the extended Kozeny–Carman–Klinkenberg model developed in our previous study successfully yielded a macroscopic model for apparent gas permeability in hierarchical micro/nanoporous systems as a function of the traditional Darcy height and porosity.

04 OIL SHALES AND TAR SANDS↗
NASA NTRS1968

An economical method for the continuous production of iodine-123

Simple and inexpensive method produces iodine 123, in a conventional cyclotron. Tellurium 122, a stable isotope available in enrichments exceeding 95 percent, is held on a porous metal plate by a flowing stream of helium and bombarded with either alpha particles or helium 3.

Blue, J. W.↗
NASA NTRS1973

Nutation damper

A nutation damper for use on a spinning body is disclosed. The damper is positioned parallel to the spin axis of the body and radially displaced therefrom. The damper is partially filled with a fluid and contains a porous media to impede the flow of the fluid induced by nutation.

Evans, J.↗
NASA NTRS1978

Experimental investigation of contamination prevention techniques for a cryogenically cooled telescope in earth-orbit

A one-tenth scale model of the Shuttle Infrared Telescope Facility (SIRTF) was constructed to experimentally investigate the feasibility of a helium purge system to reduce contamination on cryogenic surfaces. The on-orbit ambient contaminant environment was simulated with ion beams of N2(+), O(+) and H2O(+). Helium was injected into the model telescope at the simulated plane of the primary optics using two techniques: a) porous plug, b) nozzle. Helium flow rates and telescope pressures were varied to determine the purge effectiveness against potential contaminants. Based on the subscale results, a 90% reduction in O, N2, and H2O at the primary optics of the SIRTF can be obtained using 20 K helium with a pressure of .00001 torr and a flow rate of less than 0.1 g/s.

Hetrick, M. A.↗
NASA NTRS1981

Icing tunnel tests of a glycol-exuding porous leading edge ice protection system on a general aviation airfoil

Tests were conducted in the Icing Research Tunnel at the NASA Lewis Research Center to determine the characteristics of an ice protection system that distributes a glycol solution onto the leading edge of an airfoil through a porous surface material. Minimum fluid flow rates required to achieve anti-icing (no ice formation) were determined for various flight conditions and angles of attack. The ability of the system to remove ice formed on the airfoil before system activation was also investigated.

Kohlman, D. L.↗
NASA NTRS1990

An adaptive grid algorithm for one-dimensional nonlinear equations

Richards' equation, which models the flow of liquid through unsaturated porous media, is highly nonlinear and difficult to solve. Step gradients in the field variables require the use of fine grids and small time step sizes. The numerical instabilities caused by the nonlinearities often require the use of iterative methods such as Picard or Newton interation. These difficulties result in large CPU requirements in solving Richards equation. With this in mind, adaptive and multigrid methods are investigated for use with nonlinear equations such as Richards' equation. Attention is focused on one-dimensional transient problems. To investigate the use of multigrid and adaptive grid methods, a series of problems are studied. First, a multigrid program is developed and used to solve an ordinary differential equation, demonstrating the efficiency with which low and high frequency errors are smoothed out. The multigrid algorithm and an adaptive grid algorithm is used to solve one-dimensional transient partial differential equations, such as the diffusive and convective-diffusion equations. The performance of these programs are compared to that of the Gauss-Seidel and tridiagonal methods. The adaptive and multigrid schemes outperformed the Gauss-Seidel algorithm, but were not as fast as the tridiagonal method. The adaptive grid scheme solved the problems slightly faster than the multigrid method. To solve nonlinear problems, Picard iterations are introduced into the adaptive grid and tridiagonal methods. Burgers' equation is used as a test problem for the two algorithms. Both methods obtain solutions of comparable accuracy for similar time increments. For the Burgers' equation, the adaptive grid method finds the solution approximately three times faster than the tridiagonal method. Finally, both schemes are used to solve the water content formulation of the Richards' equation. For this problem, the adaptive grid method obtains a more accurate solution in fewer work units and less computation time than required by the tridiagonal method. The performance of the adaptive grid method tends to degrade as the solution process proceeds in time, but still remains faster than the tridiagonal scheme.

Gutierrez, William E.↗
NASA NTRS1991

Progress and recent developments in sodium-metal chloride batteries

Significant strides have been made in the development of high-temperature rechargeable sodium batteries utilizing transition metal chloride cathodes in the last decade, mainly due to the expertise available on Na/S batteries. These systems have already performed attractively in the various feasibility studies and have an excellent safety record. Despite the encouraging figures obtained for specific energies, certain design changes such as modifying the geometry of the beta alumina electrolyte and optimization of the porous cathodes for enhanced electrolyte flow need to be made to achieve high power densities required in applications such as electric vehicles and space. The chemistry of MCl2 cathodes, electrode fabrication, and design options are discussed, and performance data are examined.

Ratnakumar, B. V.↗
NASA NTRS1994

Capillary Movement of Liquid in Granular Beds in Microgravity

A more complete understanding of the dynamics of capillary flow through an g unsaturated porous medium would be useful for a number of space and terrestrial applications. One such application involves the development of an effective water and nutrient delivery system for the growth of plants in space. An experiment was conducted on the Mir Space Station that used an experimental cuvette called "Capillary Test Bed". An objective of this experiment was to compare fluid migration in the "Capillary Test Bed" under terrestrial laboratory simulated microgravity conditions by positioning the cuvette such that the hydrostatic force is negated and on Mir at microgravity conditions. Difference in fluid migration in the cuvette were observed. Based on these observations, improvements in the design of the cuvette are under-way that should allow for collection of more precise data in future flight experiments. The results of the MIR experiment and modifications being made to the cuvette for improved data collection will be discussed.

Yendeler, Boris S.↗