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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Engineering Redox Flow Battery Electrodes with Spatially Varying Porosity Using Non‐Solvent‐Induced Phase Separation

Redox flow batteries (RFBs) are a promising electrochemical platform for efficiently and reliably delivering electricity to the grid. Within the RFB, porous carbonaceous electrodes facilitate electrochemical reactions and distribute the flowing electrolyte. Tailoring electrode microstructure and surface area can improve RFB performance, lowering costs. Electrodes with spatially varying porosity may increase electrode utilization and provide surface area in reaction‐limited zones; however, the efficacy of such designs remains an open area of research. Herein, a non‐solvent‐induced phase‐separation (NIPS) technique that enables the reproducible synthesis of macrovoid‐free electrodes with well‐defined across‐thickness porosity gradients is described. The monotonically varying porosity profile is quantified and the physical properties and surface chemistries of porosity‐gradient electrodes are compared with macrovoid‐containing electrode, also synthesized by NIPS. Then, the electrochemical and fluid dynamic performance of the porosity‐gradient electrodes is evaluated, exploring the effect of changing the direction of the porosity gradient and benchmarking against the macrovoid‐containing electrode. Lastly, the performance is examined in a vanadium RFB, finding that the porosity‐gradient electrode outperforms the macrovoid electrode, is independent of gradient direction, and performs favorably compared to advanced electrodes in the contemporary literature. It is anticipated that the approach motivates further exploration of microstructurally tailored electrodes in electrochemical systems.

25 ENERGY STORAGE↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2023

Modeling Geologic Waste Repository Systems Below Residual Saturation

The heat generated by high-level radioactive waste can pose numerical and physical challenges to subsurface flow and transport simulators if the liquid water content in a region near the waste package approaches residual saturation due to evaporation. Here, residual saturation is the fraction of the pore space occupied by liquid water when the hydraulic connectivity through a porous medium is lost, preventing the flow of liquid water. While conventional capillary pressure models represent residual saturation using asymptotically large values of capillary pressure, here, residual saturation is effectively modeled as a tortuosity effect alone. Treating the residual fluid as primarily dead-end pores and adsorbed films, relative permeability is independent of capillary pressure below residual saturation. To test this approach, PFLOTRAN is then used to simulate thermal-hydrological conditions resulting from direct disposal of a dual-purpose canister in unsaturated alluvium using both conventional asymptotic and revised, smooth models. Importantly, while the two models have comparable results over 100 000 years, the number of flow steps required is reduced by approximately 94%.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Numerical modeling of ion transport and adsorption in porous media: A pore-scale study for capacitive deionization desalination

In this report a pore-scale model is presented to simulate the dynamic ion transport and adsorption processes in porous electrodes used for capacitive deionization (CDI). The Stokes equation governing water flow is solved using the lattice Boltzmann method and Nernst-Planck equation describing ion transport is solved using the finite volume method. The ion adsorption process is considered at the surface of carbon electrodes. After validation against analytical solutions and published results, the model is used to study the coupled water flow, ion transport and adsorption in both two-dimensional and three-dimensional porous CDI electrodes at the pore scale. The effect of electrode microstructure, electrical potential and flow velocity on the adsorption processes is quantitatively investigated, and the relative importance of various parameters is determined. The presented model can be a powerful numerical tool to quantitatively analyze ion transport and adsorption in porous electrodes, and may provide useful information for the design and optimization of CDI electrodes and operating conditions for desired desalination efficiency, water throughput and cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Pressure Drop Correlation Improvement for the Near-Wall Region of Pebble-Bed Reactors

Packed beds play an important role in several engineering fields, with their applications in nuclear energy being driven by the development of next-generation reactors utilizing pebble fuel. The random nature of a packed pebble bed creates a flow field that is complex and difficult to predict. Porous media models are an attractive option for modeling pebble-bed reactors (PBRs), as they provide intermediate fidelity results and are computationally efficient. Porous media models, however, rely on the use of correlations to estimate the effect of complicated flow features on the pressure drop and heat transfer in the system. Existing correlations were developed to predict the average behavior of the bed, but they are inaccurate in the near-wall region where the presence of the wall affects the pebble packing. This work aims to investigate the accuracy of a porous media model using the Kerntechnischer Ausschuss (KTA) correlation, the most common pressure drop correlation for PBRs compared to the high-fidelity large eddy simulation (LES). A bed of 1568 pebbles is investigated at Reynolds numbers from 625 to 10 000. The bed is divided into five concentric subdomains to compare the average velocity, friction losses, and form losses between the porous media and LES codes. The comparison between the LES simulation and the KTA correlation revealed that the KTA correlation largely underpredicts the form losses in the near-wall region, leading to an overprediction of the velocity near the wall by nearly 30%. An investigation of the form losses across the range of Reynolds numbers in the LES results provided additional insight into how the KTA correlation may be improved to better predict these spatial effects in a pebble bed. These data suggest that the form coefficient near the wall must be increased by 48% while decreasing the form coefficient of the inner bulk region of the bed by 15%. The implementation of these improvements to the KTA correlation in a porous media model produced a radial velocity profile that saw significantly improved agreement with the LES results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗
DOE OSTI2025

Applying a Compact Porous Media Model to Numerically Derive Resistance Coefficients for Lattice Structures

Additive Manufacturing allows for exploring various geometries to achieve specific engineering criteria. Lattices are one geometry with unique properties, including being periodically repeating structures which allow flow through them to be represented as a porous media according to Darcy-Forchheimer equations. These equation’s coefficients are generally experimentally derived, but this work demonstrates the ability to numerically derive them with CFD. Simulations were performed using three-dimensional stead state Reynolds-averaged Navier-Stokes with a k-ω Shear Stress Transport turbulence model using Ansys Fluent. Three lattice geometries were investigated and drag coefficients were derived. The method was validated against externally published data for similar geometries demonstrating strong agreement, and grid convergence for all simulations was calculated with a Grid Convergence Index method. Wall roughness is demonstrated to have a non-negligible impact on results and roughness values are considered for the primary focus Octahedral geometry where both smooth wall and rough wall coefficients were derived. The porosity coefficients for the Octahedral geometry at 1.0 [m/s] were found to be 2.89×10 6 and 2.90×10 6 [1/(Pa*m*s)] for the permeability coefficients, 6.37×10 1 and 5.44×10 1 [m 2 /kg] for the inertial resistance coefficients, and with a max pressure drop of 5116.7 [Pa] and 4429.5 [Pa] for the smooth walls and rough walls, respectively. The derived numerical method enables rapid exploration and optimization of new lattice designs for diverse engineering applications.

42 ENGINEERING↗
DOE OSTI2024

Evaporation from thin porous Coatings: Pore size effects and predictive equation for homogeneous coatings

Evaporation of small water droplets on solids is hindered because surface tension pulls the droplet into a spherical cap that has a small perimeter. Our solution is to coat a solid with a very thin, porous layer into which the droplet flows to create a large-area disk with concomitant high rate of evaporation. We investigate evaporation by varying factors that have not been previously considered: pore size and distribution, contact angle, temperature, and relative humidity (RH). A larger pore size resulted in faster evaporation, which we explain through faster transport within the coating. Even faster evaporation occurred for a bilayer structure with small particles on the air side and larger particles on the solid side. Further, the water advancing contact angle had an insignificant effect in the range from < 10° through to 60°. Our results for different pore sizes, temperature, humidity, and contact angle all collapse onto a single curve when appropriately normalized. This validates an equation that can be used for the evaporation from a homogeneous coating that depends only one empirical factor and the droplet volume. Since the volume is often user-controlled, we envisage that this equation can be used to predict evaporation and guide design of fast-drying coatings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Pore-scale simulation of multiphase flow and reactive transport processes involved in geologic carbon sequestration

Multiphase flow and reactive transport are two essential physicochemical processes that govern the effectiveness of geological carbon sequestration (GCS). The interaction and feedback among different phases and components during intricate physicochemical processes hold great significance in understanding CO 2 sequestration. Pore-scale simulations can account for multiphase flow and reactive transport processes in porous media and obtain spatial distributions of parameters (density, velocity, concentration, etc.) in the pore space as well as their temporal evolutions. This proves especially valuable considering that experiments can be hindered by constraints in spatial and temporal resolution. The comprehensive insights garnered from pore-scale research can be leveraged for continuum modeling using the representative elementary volume (REV) concept. In this contribution, four sequential mechanisms of CO 2 -brine-rock interaction in three zones delineated by CO 2 saturation are elaborated to elucidate complicated physicochemical processes involved in GCS, which are followed by general descriptions of mathematical equations and pore-scale numerical methods. In addition, as interested and commonly encountered processes, leakage risks during GCS and CO 2 -enhanced oil recovery (CO 2 -EOR) processes are presented. The existing challenges and future directions are discussed for both the performance of the pore-scale models and the current gaps in the field of GCS. Importantly, we expect that this review will prove beneficial for researchers interested in pore-scale simulations, GCS, and related disciplines.

58 GEOSCIENCES↗
DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗
DOE OSTI2021

Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks

Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).

42 ENGINEERING↗
DOE OSTI2020

Calculation of Hydraulic Conductivity for a Porous Medium with Staggered Array of Solid Grains - 20087

For a porous medium with periodic array of solid grains. the flow of an incompressible viscous fluid in the pore space is considered with an aim of calculating the permeability of the medium. The solid grains are of ellipsoidal shape. For this purpose, a unit-cell boundary-value problem (Stokes problem) is solved numerically. With the assumption of periodicity of the medium structure and the fluid properties (in fact constant), the flow field and pressure distributions are determined numerically with inertial effects ignored. By applying the theory of homogenization to the viscous flow through pore space, a microscale boundary-value problem is defined. After solving the problem, the permeability is determined by taking the micro-cell average of the velocity components, Various porous structures are chosen to allow a wide range of porosity values. The permeability increases mildly for porosity values in the lower range. But it increases sharply for porosity values in the upper range of the porosity. The permeability of the present study agrees with the Kozeny-Carman relation excellently for smaller porosity. Certain discrepancy between the calculated permeability in this study and the Kozeny-Carman relation appears and increases with porosity for larger porosity. It is believed that the discrepancy is due to the replacement of the pores by a bundle of capillaries, which is unrealistic, in the Kozeny-Carman relation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

T2DM v1.0

T2DM adds two-dimensional hydrodynamic dispersion to TOUGH2/EOS7C. T2DM can be used when modeling flow and transport of brine in porous media to account for hydrodynamic dispersion.

Oldenburg, CurtisM↗
DOE OSTI2023

Scaling hyporheic nitrogen cycling in large river alluvial aquifers

Research Element 1: Denitrification is a difficult flux to measure when fluxes are low and when the ecosystem is difficult to sample. These attributes describe Nyack floodplain, a large gravel-bedded aquifer adjacent to the Middle Fork Flathead River. Given that denitrification is simply respiration with nitrate as a terminal electron acceptor, we reasoned that N 2 production should negatively covary with oxygen respiration so long that oxygen is available. We sampled 16 wells 8 times throughout a year and analyzed these water samples on a membrane-inlet mass spectrometer for N 2 :Ar and O 2 . We fit multilevel models to N 2 :Ar as a function of dissolved oxygen. Within wells, there was no relationship between O 2 and N 2 :Ar. Among wells N 2 :Ar increased as O 2 declined. Based on the slope of this relationship, ~ -0.07, and estimated oxygen respiration in Nyack derived from prior studies, a possible range for denitrification is 0-8 mmol m -2 d -1 corresponding to a turnover rate for nitrate of 0-0.2/d. This estimated turnover rate is in the lower 20% of rates measured in the LINX 2 study, yet shows that even in a low nitrate (0.005 mmol/L) aquifer denitrification can be a substantial sink for dissolved N. Research Element 3 and 4: We developed a new theoretical approach to describe the relationship between water inflow, water outflow, and water age in porous media control volumes under dynamic flow conditions. Our approach had three distinct steps. First we described the application of well-established equations borrowed from Chemical Engineering to describe water age distributions in the hyporheic zone under steady state conditions. We published the associated results in the scientific journal PLOS ONE. Second, we developed a new theory ("flow spans") to describe how predicted steady-state water age distributions will respond to time varying imbalances in hyporheic discharge and recharge. Finally, we developed a proof-of-concept simulation model based on the new theory and used the model to examine the effects of flow regulation on a simulated river with characteristics representative of western U.S. montane rivers with expansive hyporheic zones.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Machine Learning Based Prediction of Airflow Maldistribution in A-Type Heat Exchangers

Airflow maldistribution is one of the primary causes of performance degradation in air-to refrigerant heat exchangers (HX) and has been shown to decrease heat transfer by as much as 35%. As a result, many units are oversized to meet the target capacity, resulting in increased system cost and refrigerant charge. Several studies have explored how characteristics like package type and HX geometry impact the flow profile, but results are restricted to a limited range of parameters and cannot be extrapolated to new designs. In this work, a machine learning (ML) model is trained to predict the inlet flow profile of dry air entering A-type HXs across a broad range of geometries and conditions. Flow profiles are generated using a porous media CFD model and used to train an Artificial Neural Network (ANN) which exhibits maximum and average relative L2 norm errors of 0.48 and 0.05. Additionally, these predictions take less than a second to generate resulting in a speed up factor of 2.42E5 compared to CFD. Component-level simulations are conducted to determine the performance degradation resulting from the predicted airflow maldistribution profiles. The new ML model will enable rapid and accurate prediction of performance degradation resulting from airflow maldistribution in A-type HXs, allowing for more accurate and cost-effective HX design.

42 ENGINEERING↗
DOE OSTI2024

Quartz Dissolution Effects on Flow Channelization and Transport Behavior in Three‐Dimensional Fracture Networks

We perform a set of reactive transport simulations in three-dimensional fracture networks to characterize the impact of geochemical reactions on flow channelization. Flow channelization, a frequently observed phenomenon in porous and fractured subsurface rock formations, results from the spatially variable hydraulic resistance offered by a geological structure. In addition to geo-structural features such as network connectivity, geometry, and hydraulic resistance, geochemical reactions, for example, dissolution and precipitation, can dynamically inhibit or enhance flow channelization. These geochemical processes can change the fracture permeability leading to increased flow channelization, which are localized connected regions of high volumetric flow rates that are seemingly ubiquitous in the subsurface. In our simulations, fractures partially filled with quartz are gradually dissolved until quasi-steady state conditions are obtained. We compare the flow field's initial unreacted and final dissolved states in terms of flow and transport observations. We observe that the dissolved fracture networks provide less resistance to flow and exhibit increased flow channelization when compared to their unreacted counterparts. However, there is substantial variability in the magnitude of these changes which implies that the channelization strongly depends on the network structure. In turn, we identify the interplay between the particular network structure and the impact of geochemical dissolution on flow channelization. The presented results indicate that geological systems that have been weathering or reactive for longer times in older landscapes are likely to have increased flow channelization compared to their equivalent but younger counterparts, which implies a time dependence on flow channelization in fractured media.

Hyman, Jeffrey D.↗
DOE OSTI2022

Less can be more: Insights on the role of electrode microstructure in redox flow batteries from two-dimensional direct numerical simulations

Understanding how to structure a porous electrode to facilitate fluid, mass, and charge transport is key to enhancing the performance of electrochemical devices, such as fuel cells, electrolyzers, and redox flow batteries (RFBs). Here, using a parallel computational framework, direct numerical simulations are carried out on idealized porous electrode microstructures for RFBs. Strategies to improve an electrode design starting from a regular lattice are explored. By introducing vacancies in the ordered arrangement, it is possible to achieve higher voltage efficiency at a given current density, thanks to improved mixing of reactive species, despite reducing the total reactive surface. Careful engineering of the location of vacancies, resulting in a density gradient, outperforms disordered configurations. Our simulation framework is a new tool to explore transport phenomena in RFBs, and our findings suggest new ways to design performant electrodes.

25 ENERGY STORAGE↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗