Comparison of Control Volume Analysis and Porous Media Averaging for Formulation of Porous Media Transport

Author(s):  
F. Civan
2021 ◽  
Author(s):  
S.R. Zhu ◽  
L.Z. Wu ◽  
T. Ma ◽  
S.H. Li

Abstract The numerical solution of various systems of linear equations describing fluid infiltration uses the Picard iteration (PI). However, because many such systems are ill-conditioned, the solution process often has a poor convergence rate, making it very time-consuming. In this study, a control volume method based on non-uniform nodes is used to discretize the Richards equation, and adaptive relaxation is combined with a multistep preconditioner to improve the convergence rate of PI. The resulting adaptive relaxed PI with multistep preconditioner (MP(m)-ARPI) is used to simulate unsaturated flow in porous media. Three examples are used to verify the proposed schemes. The results show that MP(m)-ARPI can effectively reduce the condition number of the coefficient matrix for the system of linear equations. Compared with conventional PI, MP(m)-ARPI achieves faster convergence, higher computational efficiency, and enhanced robustness. These results demonstrate that improved scheme is an excellent prospect for simulating unsaturated flow in porous media.


2021 ◽  
Vol 929 ◽  
Author(s):  
M. Ungarish

Previous studies have extended Benjamin's theory for an inertial steady-state gravity current of density $\rho _{c}$ in a homogeneous ambient fluid of density $\rho _{o} < \rho _{c}$ to the counterpart propagation in a linearly stratified (Boussinesq) ambient (density decreases from $\rho _b$ to $\rho _{o}$ ). The extension is typified by the parameter $S = (\rho _{b}-\rho _{o})/(\rho _{c}-\rho _{o}) \in (0,1]$ , uses Long's solution for the flow over a topography to model the flow of the ambient over the gravity current, and reduces well to the classical theory for small and moderate values of $S$ . However, for $S=1$ , i.e. $\rho _b = \rho _c$ , which corresponds to a symmetric intrusion, various idiosyncrasies appear. Here attention is focused on this case. The control-volume analysis (balance of volume, mass, momentum and vorticity) produces a fairly compact analytical formulation, pending a closure for the head loss, and subject to stability criteria (no inverse stratification downstream). However, we show that plausible closures that work well for the non-stratified current (like zero head loss on the stagnation line, or zero vorticity diffusion) do not produce satisfactory results for the intrusion (except for some small ranges of the height ratio of current to channel, $a = h/H$ ). The reasons and insights are discussed. Accurate data needed for comparison with the theoretical model are scarce, and a message of this paper is that dedicated experiments and simulations are needed for the clarification and improvement of the theory.


In this chapter, the non-Darcy model is employed for porous media filled with nanofluid. Both natural and forced convection heat transfer can be analyzed with this model. The governing equations in forms of vorticity stream function are derived and then they are solved via control volume-based finite element method (CVFEM). The effect of Darcy number on nanofluid flow and heat transfer is examined.


2019 ◽  
Vol 2019 ◽  
pp. 1-16 ◽  
Author(s):  
Hong Thai Vu ◽  
Evangelos Tsotsas

The modelling and numerical simulation of the drying process in porous media are discussed in this work with the objective of presenting the drying problem as the system of governing equations, which is ready to be solved by many of the now widely available control-volume-based numerical tools. By reviewing the connection between the transport equations at the pore level and their up-scaled ones at the continuum level and then by transforming these equations into a format that can be solved by the control volume method, we would like to present an easy-to-use framework for studying the drying process in porous media. In order to take into account the microstructure of porous media in the format of pore-size distribution, the concept of bundle of capillaries is used to derive the needed transport parameters. Some numerical examples are presented to demonstrate the use of the presented formulas.


Sign in / Sign up

Export Citation Format

Share Document