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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Microblowing Technique Demonstrated to Reduce Skin Friction

One of the most challenging areas of research in aerodynamics is the reduction of skin friction, especially for turbulent flow. Reduced skin friction means less drag. For aircraft, less drag can lead to less fuel burned or to a greater flight range for a fixed amount of fuel. Many techniques and methods have been tried; however, none of them has significantly reduced skin friction in the flight environment. An innovative skin-friction reduction technique, the Microblowing Technique (MBT), was invented in 1993. This is a unique concept in which an extremely small amount of air is blown vertically at a surface through very small holes. It can be used for aircraft or marine vehicles, such as submarines (where water is blown through the holes instead of air). As shown in the figure, the outer layer, which controls vertical flow, is a plate with high-vertical holes. The inner layer, which produces evenly distributed flow, is a low-permeability porous plate. Microblowing reduces the surface roughness and changes the flow velocity profile on the surface, thereby reducing skin friction.

Hwang, Danny P.↗
NASA NTRS1997

The Effect of Gravity on the Combustion Synthesis of Porous Ceramics and Metal Matrix Composites

Combustion synthesis (self propagating, high temperature synthesis-SHS) is a novel technique that is capable of producing many advanced materials. The ignition temperature (Tig) of such combustion synthesis reactions is often coincident with that of the lowest melting point reactant. The resultant liquid metal wets and spreads around the other solid reactant particles of higher melting points, thereby improving the reactant contact and kinetics, followed by formation of the required compounds. This ignition initiates a combustion propagating wave whose narrow reaction front rapidly travels through the reactants. Since this process is highly exothermic, the heat released by combustion often melts the reactant particles ahead of the combustion front and ignites the adjacent reactant layer, resulting in a self-sustaining reaction. Whenever a fluid phase (liquid or gas) is generated by the reaction system, gravity-driven phenomena can occur. Such phenomena include convective flows of fluid by conventional or unstable convection and settling of the higher density phases. A combustion process is often associated with various kinds of fluid flow. For instance, if the SHS reaction is carried out under inert or reactive gas atmospheres, or a volatile, e.g., B2O3, is deliberately introduced as a reactant, convective flows of the gas will occur due to a temperature gradient existing in the atmosphere when a combustion wave is initiated. The increased gas flow will produce a porous (or expanded) SHS product. Owing to the highly exothermic nature of many SHS reactions, liquid phase(s) can also form before, at, or after the combustion front. The huge temperature gradient at the combustion front can induce convective flows (conventional or unstable) of the liquid phase. Each of these types of convective fluid flow can change the combustion behavior of the synthesizing reaction, and, therefore, the resultant product microstructure. In addition, when two or more phases of different density are produced at or ahead of the propagating combustion front settling of the higher density phase will occur resulting in a non-uniform product microstructure and properties.

Moore, J. J.↗
DOE OSTI2023

Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations

In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Electrochemical Residence Time Distribution as a Diagnostic Tool for Redox Flow Batteries

The fluid dynamic and electrochemical performance of redox flow batteries (RFBs) stems from the relationship between the flow field and the porous electrode, whose interplay determines how active species move and react during device operation. While characterization techniques, such as residence time distribution, offer insights into species mobility within a reactive volume for a traditional chemical reactor, electrochemical reactors also enable simultaneous measurement of the redox reactions, unlocking another dimension of analysis. Herein, we demonstrate how potentiodynamic measurements, using injections of electrolyte examined through moment analysis, can provide electrode-specific performance scaling relationships across a matrix of carbon paper and cloth electrodes with flow through and interdigitated flow fields. We further combine experimental campaigns with multiphysics simulations to demonstrate how electrode surface area can be estimated with this technique, which we then validate with activated and unactivated commercial carbon cloth electrodes. These studies reveal the multiscale observations that potentiodynamic measurements afford, augmenting existing electrochemical techniques for holistic electrochemical reactor diagnostics.

Electrochemistry↗
DOE OSTI2022

June 2022 Permeability Measurements in P-Tunnel, Nevada National Security Site

This document describes gas permeability measurements made in June 2022 in P-Tunnel at the Nevada National Security Site (NNSS). Four boreholes (AC-1, GI-5, GI-6, and HF-1) were tested with a dual packer system. The pressure versus flow rate data that were collected are analyzed using a numerical simulation method that captures compressible gas flow in partially saturated porous media using a close approximation to the geometry of the test setup. The first part of this document begins with data collection and a description of the packer assembly followed by details on data acquisition, measurement instrumentation, and equipment followed by a description of pre-test system verification. Next, the implementation of testing in the field is described, including how specific locations were chosen for testing. The second part of this document describes the numerical simulations used to back-out permeability from the pressure versus flow rate data. First, the numerical model is described, including assumptions. Next the details of the mesh that represents the experiment are described. A method for generating type curves in pressure versus flow rate space is presented and finally the data are mapped to permeability. We conclude with a discussion of the results.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

Development of a hydraulic model of the human systemic circulation

Physical and numeric models of the human circulation are constructed for a number of objectives, including studies and training in physiologic control, interpretation of clinical observations, and testing of prosthetic cardiovascular devices. For many of these purposes it is important to quantitatively validate the dynamic response of the models in terms of the input impedance (Z = oscillatory pressure/oscillatory flow). To address this need, the authors developed an improved physical model. Using a computer study, the authors first identified the configuration of lumped parameter elements in a model of the systemic circulation; the result was a good match with human aortic input impedance with a minimum number of elements. Design, construction, and testing of a hydraulic model analogous to the computer model followed. Numeric results showed that a three element model with two resistors and one compliance produced reasonable matching without undue complication. The subsequent analogous hydraulic model included adjustable resistors incorporating a sliding plate to vary the flow area through a porous material and an adjustable compliance consisting of a variable-volume air chamber. The response of the hydraulic model compared favorably with other circulation models.

NASA Program Biomedical Research and Countermeasur↗
DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗
NASA NTRS1990

The liquid/vapor boundary in a porous plug

Results are presented from measurements of pressure and temperature differences as a function of the flow rate and hydrostatic head for two types of porous plugs for phase separation of superfluid He and He vapor: (1) a plug made of sintered stainless steel, similar to that in the IRAS, IRT, and COBE dewars except for larger pores, and (2) a glass plug with pores of uniform cross section. Results show that a liquid/vapor phase boundary exists at the outer surface of the porous plug for low flows or large hydrostatic heads; as the flow increases or the hydrostatic head decreases, the boundary retreats into the plug. Some hysteresis near the phase boundary retreat was observed in the stainless steel plug.

Dipirro, M. J.↗
NASA NTRS1993

Evaluation of a distributed catchment scale water balance model

The validity of some of the simplifying assumptions in a conceptual water balance model is investigated by comparing simulation results from the conceptual model with simulation results from a three-dimensional physically based numerical model and with field observations. We examine, in particular, assumptions and simplifications related to water table dynamics, vertical soil moisture and pressure head distributions, and subsurface flow contributions to stream discharge. The conceptual model relies on a topographic index to predict saturation excess runoff and on Philip's infiltration equation to predict infiltration excess runoff. The numerical model solves the three-dimensional Richards equation describing flow in variably saturated porous media, and handles seepage face boundaries, infiltration excess and saturation excess runoff production, and soil driven and atmosphere driven surface fluxes. The study catchments (a 7.2 sq km catchment and a 0.64 sq km subcatchment) are located in the North Appalachian ridge and valley region of eastern Pennsylvania. Hydrologic data collected during the MACHYDRO 90 field experiment are used to calibrate the models and to evaluate simulation results. It is found that water table dynamics as predicted by the conceptual model are close to the observations in a shallow water well and therefore, that a linear relationship between a topographic index and the local water table depth is found to be a reasonable assumption for catchment scale modeling. However, the hydraulic equilibrium assumption is not valid for the upper 100 cm layer of the unsaturated zone and a conceptual model that incorporates a root zone is suggested. Furthermore, theoretical subsurface flow characteristics from the conceptual model are found to be different from field observations, numerical simulation results, and theoretical baseflow recession characteristics based on Boussinesq's groundwater equation.

Troch, Peter A.↗
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↗
DOE OSTI2024

Quantifying local and global mass balance errors in physics-informed neural networks

Physics-informed neural networks (PINN) have recently become attractive for solving partial differential equations (PDEs) that describe physics laws. By including PDE-based loss functions, physics laws such as mass balance are enforced softly in PINN. This paper investigates how mass balance constraints are satisfied when PINN is used to solve the resulting PDEs. We investigate PINN’s ability to solve the 1D saturated groundwater flow equations (diffusion equations) for homogeneous and heterogeneous media and evaluate the local and global mass balance errors. We compare the obtained PINN’s solution and associated mass balance errors against a two-point finite volume numerical method and the corresponding analytical solution. We also evaluate the accuracy of PINN in solving the 1D saturated groundwater flow equation with and without incorporating hydraulic heads as training data. We demonstrate that PINN’s local and global mass balance errors are significant compared to the finite volume approach. Tuning the PINN’s hyperparameters, such as the number of collocation points, training data, hidden layers, nodes, epochs, and learning rate, did not improve the solution accuracy or the mass balance errors compared to the finite volume solution. Mass balance errors could considerably challenge the utility of PINN in applications where ensuring compliance with physical and mathematical properties is crucial.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

PFLOTRAN-SIP: A PFLOTRAN Module for Simulating Spectral-Induced Polarization of Electrical Impedance Data

Spectral induced polarization (SIP) is a non-intrusive geophysical method that collects chargeability information (the ability of a material to retain charge) in the time domain or its phase shift in the frequency domain. Although SIP is a temporal method, it cannot measure the dynamics of flow and solute/species transport in the subsurface over long times (i.e., 10–100 s of years). Data collected with the SIP technique need to be coupled with fluid flow and reactive-transport models in order to capture long-term dynamics. To address this challenge, PFLOTRAN-SIP was built to couple SIP data to fluid flow and solute transport processes. Specifically, this framework couples the subsurface flow and transport simulator PFLOTRAN and geoelectrical simulator E4D without sacrificing computational performance. PFLOTRAN solves the coupled flow and solute-transport process models in order to estimate solute concentrations, which were used in Archie’s model to compute bulk electrical conductivities at near-zero frequency. These bulk electrical conductivities were modified while using the Cole–Cole model to account for frequency dependence. Using the estimated frequency-dependent bulk conductivities, E4D simulated the real and complex electrical potential signals for selected frequencies for SIP. These frequency-dependent bulk conductivities contain information that is relevant to geochemical changes in the system. This study demonstrated that the PFLOTRAN-SIP framework is able to detect the presence of a tracer in the subsurface. SIP offers a significant benefit over ERT in the form of greater information content. It provided multiple datasets at different frequencies that better constrained the tracer distribution in the subsurface. Consequently, this framework allows for practitioners of environmental hydrogeophysics and biogeophysics to monitor the subsurface with improved resolution.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

A Novel Approach for Reducing Rotor Tip-Clearance Induced Noise in Turbofan Engines

Rotor tip-clearance induced noise, both in the form of rotor self noise and rotor-stator interaction noise , constitutes a significant component of total fan noise. Innovative yet cost effective techniques to suppress rotor-generated noise are, therefore, of foremost importance for improving the noise signature of turbofan engines. To that end, the feasibility of a passive porous treatment strategy to positively modify the tip-clearance flow field is addressed. The present study is focused on accurate viscous flow calculations of the baseline and the treated rotor flow fields. Detailed comparison between the computed baseline solution and experimental measurements shows excellent agreement. Tip-vortex structure, trajectory, strength, and other relevant aerodynamic quantities are extracted from the computed database. Extensive comparison between the untreated and treated tip-clearance flow fields is performed. The effectiveness of the porous treatment for altering the rotor-tip vortex flow field in general and reducing the intensity of the tip vortex, in particular, is demonstrated. In addition, the simulated flow field for the treated tip clearly shows that substantial reduction in the intensity of both the shear layer roll-up and boundary layer separation on the wall is achieved.

Khorrami, Mehdi R.↗
DOE OSTI2021

Geometrical-Based Generative Adversarial Network to Enhance Digital Rock Image Quality

X-ray microcomputed tomography (micro-CT) is a common tool for the study of porous media structures and properties. High-quality micro-CT data are required to accurately capture pore structures. Acquiring high-quality micro-CT data, however, is not always possible, owing to application limitations and experimental constraints. Therefore, we propose a geometrical-based generative adversarial network (GAN) to rapidly restore noisy micro-CT images to their clean counterparts. The training data and related ground-truth (GT) data are scanned for 7 min and 9.5 h, respectively. To evaluate the performance of the geometrical-based GAN, a 6003 voxel image that has never been used for training is reconstructed and compared with the corresponding GT image. Histogram matching and linear normalization are implemented to adjust the histogram of the reconstructed image to that of the GT image. A watershed-based segmentation method is then applied to delineate pore and solid phases. Lastly, we measure the Minkowski functionals and petrophysical properties, including absolute permeability, pore size distribution, drainage capillary pressure-saturation curve, and imbibition relative permeability, to estimate the physical accuracy of the denoised image. The results show that the proposed geometrical-based GAN can accurately restore noisy micro-CT data. By reducing the scanning time from 9.5 h to 7 min, the expenditure of collecting micro-CT can be decreased significantly. This is particularly important for applications where time-lapse images of a dynamic process are required, high-throughput imaging is necessary for real-time data analysis, or where the quantification of large sample volumes is required.

58 GEOSCIENCES↗
NASA NTRS1992

An Experimental Study of the Fluid Mechanics Associated With Porous Walls

The fluid mechanics of air exiting from a porous material is investigated. The experiments are filter rating dependent, as porous walls with filter ratings differing by about three orders of magnitude are studied. The flow behavior is investigated for its spatial and temporal stability. The results from the investigation are related to jet behavior in at least one of the following categories: (1) jet coalescence effects with increasing flow rate; (2) jet field decay with increasing distance from the porous wall; (3) jet field temporal turbulence characteristics; and (4) single jet turbulence characteristics. The measurements show that coalescence effects cause jet development, and this development stage can be traced by measuring the pseudoturbulence (spatial velocity variations) at any flow rate. The pseudoturbulence variation with increasing mass flow reveals an initial increasing trend followed by a leveling trend, both of which are directly proportional to the filter rating. A critical velocity begins this leveling trend and represents the onset of fully developed jetting action in the flow field. A correlation is developed to predict the onset of fully developed jets in the flow emerging from a porous wall. The data further show that the fully developed jet dimensions are independent of the filter rating, thus providing a length scale for this type of flow field (1 mm). Individual jet characteristics provide another unifying trend with similar velocity decay behavior with distance; however, the respective turbulence magnitudes show vast differences between jets from the same sample. Measurements of the flow decay with distance from the porous wall show that the higher spatial frequency components of the jet field dissipate faster than the lower frequency components. Flow turbulence intensity measurements show an out of phase behavior with the velocity field and are generally found to increase as the distance from the wall is increased.

N Ramachandran↗
DOE OSTI2020

Modeling CO 2 plume migration using an invasion-percolation approach that includes dissolution

Geological carbon dioxide (CO 2 ) sequestration has been proposed as a viable technique to decrease effective emissions of CO 2 into the atmosphere. However, the security of this sequestration is tied to our understanding of the long-term migration of CO 2 in subsurface. The dissolution of CO 2 in the reservoir brine is one of the main long-term trapping mechanisms. However, the assumptions used in large-scale reservoir simulations usually lead to an overestimation of the dissolution volume. We propose a modified approach based on the macroscopic invasion-percolation (MIP) theory that allows the dissolution of CO 2 into brine. We used a high-resolution geological model to compare the Darcy-, modified MIP-, and classic MIP-based simulation results. We observed a significant shrinkage in the nonaqueous plume volume when dissolution is considered during the MIP simulation. In the case of Darcy-based simulation, the plume was completely trapped inside the reservoir with limited migration even after a thousand-year simulation. On the other hand, the majority of the plume migrated out of the simulated reservoir in the case of MIP. Our approach provides more realistic estimation of the dissolution volume and nonaqueous plume extent while leveraging the computational efficiency enjoyed by MIP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An experimental investigation on the CO 2 storage capacity of the composite confining system

Assuring secure containment of stored CO 2 is of paramount importance—for climate change mitigation, for permitting, and for reassuring the public. Regional seals, such as those sealing petroleum accumulations, have proven to be effective for securing CO 2 . However, the goal of permanent sequestration can also be satisfied with “composite confining systems” consisting of multiple, possibly discontinuous flow barriers that, in aggregate, create a system with very high permeability anisotropy and effectively retard the vertical migration of CO 2 . This study focuses on investigating the barrier characteristics necessary for effective containment of CO 2 , using both physical flow experiments and modified invasion percolation simulations. The simulations are calibrated to the physical experiments and are used to further extend the analysis. Results show that for a composite confining system, a) even barriers with low capillary entry pressure contrast to the underlying flow unit can divert rising CO 2 ; b) curved or anticlinal barrier topography can enhance CO 2 trapping; c) fining-upward gradations make little difference to barrier effectiveness or CO 2 retention; d) longer barriers retain more CO 2 regardless of barrier topography. Finally, field-scale simulations have demonstrated the importance of barrier length for increasing CO 2 storage capacity. Furthermore, the results presented can be used to inform the development of new screening criteria for characterization and effectiveness of composite confining systems.

58 GEOSCIENCES↗