Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userViscous Transport in Eroding Porous Media
106
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Transport of viscous fluid through porous media is a direct consequence of the pore structure. Here we investigate transport through a specific class of two-dimensional porous geometries, namely those formed by fluid-mechanical erosion. We investigate the tortuosity and dispersion by analyzing the first two statistical moments of tracer trajectories. For most initial configurations, tortuosity decreases in time as a result of erosion increasing the porosity. However, we find that tortuosity can also increase transiently in certain cases. The porosity-tortuosity relationships that result from our simulations are compared with models available in the literature. Asymptotic dispersion rates are also strongly affected by the erosion process, as well as by the number and distribution of the eroding bodies. Finally, we analyze the pore size distribution of an eroding geometry. The simulations are performed by combining a high-fidelity boundary integral equation solver for the fluid equations, a second-order stable time stepping method to simulate erosion, and new numerical methods to stably and accurately resolve nearly-touching eroded bodies and particle trajectories near the eroding bodies.
rate research
Read More
We extend unsteady thin aerofoil theory to model aerofoils with generalised chordwise porosity distributions. The analysis considers a linearised porosity boundary condition where the seepage velocity through the aerofoil is related to the local pressure jump across the aerofoil surface and to the unsteady characteristics of the porous medium. Application of the Plemelj formulae to the resulting boundary value problem yields a singular Fredholm--Volterra integral equation which does not admit an analytic solution. Accordingly, we develop a numerical solution scheme by expanding the bound vorticity distribution in terms of appropriate basis functions. Asymptotic analysis at the leading- and trailing-edges reveals that the appropriate basis functions are weighted Jacobi polynomials whose parameters are related to the porosity distribution. The Jacobi polynomial basis enables the construction of a numerical scheme that is accurate and rapid, in contrast to the standard choice of Chebyshev basis functions that are shown to {be unsuitable} for porous aerofoils. Applications of the numerical solution scheme to discontinuous porosity profiles, quasi-static problems, and the separation of circulatory and non-circulatory contributions are presented. Analogues to the classical Theodorsen and Sears functions are computed numerically, which show that an effect of trailing-edge porosity is to reduce the amount of vorticity shed into the wake, thereby reducing the magnitude of the circulatory lift. {Fourier transform inversion of these frequency-domain functions produces porous extensions to the Wagner and K{u}ssner functions for transient aerofoil motions or gust encounters, respectively.}
We present a novel computational modeling framework to numerically investigate fluid-structure interaction in viscous fluids using the phase field embedding method. Each rigid body or elastic structure immersed in the incompressible viscous fluid matrix, grossly referred to as the particle in this paper, is identified by a volume preserving phase field. The motion of the particle is driven by the fluid velocity in the matrix for passive particles or combined with its self-propelling velocity for active particles. The excluded volume effect between a pair of particles or between a particle and the boundary is modeled by a repulsive potential force. The drag exerted to the fluid by a particle is assumed proportional to its velocity. When the particle is rigid, its state is described by a zero velocity gradient tensor within the nonzero phase field that defines its profile and a constraining stress exists therein. While the particle is elastic, a linear constitutive equation for the elastic stress is provided within the particle domain. A hybrid, thermodynamically consistent hydrodynamic model valid in the entire computational domain is then derived for the fluid-particle ensemble using the generalized Onsager principle accounting for both rigid and elastic particles. Structure-preserving numerical algorithms are subsequently developed for the thermodynamically consistent model. Numerical tests in 2D and 3D space are carried out to verify the rate of convergence and numerical examples are given to demonstrate the usefulness of the computational framework for simulating fluid-structure interactions for passive as well as self-propelling active particles in a viscous fluid matrix.
We investigate the elastoviscoplastic flow through porous media by numerical simulations. We solve the Navier-Stokes equations combined with the elastoviscoplastic model proposed by Saramito for the stress tensor evolution. In this model, the material behaves as a viscoelastic solid when unyielded, and as a viscoelastic Oldroyd-B fluid for stresses higher than the yield stress. The porous media is made of a symmetric array of cylinders, and we solve the flow in one periodic cell. We find that the solution is time-dependent even at low Reynolds numbers as we observe oscillations in time of the unyielded region especially at high Bingham numbers. The volume of the unyielded region slightly decreases with the Reynolds number and strongly increases with the Bingham number; up to 70% of the total volume is unyielded for the highest Bingham numbers considered here. The flow is mainly shear dominated in the yielded region, while shear and elongational flow are equally distributed in the unyielded region. We compute the relation between the pressure drop and the flow rate in the porous medium and present an empirical closure as function of the Bingham and Reynolds numbers. The apparent permeability, normalized with the case of Newtonian fluids, is shown to be greater than 1 at low Bingham numbers, corresponding to lower pressure drops due to the flow elasticity, and smaller than 1 for high Bingham numbers, indicating larger dissipation in the flow owing to the presence of the yielded regions. Finally we investigate the effect of the Weissenberg number on the distribution of the unyielded regions and on the pressure gradient.
This paper presents a time-space Hausdorff derivative model for depicting solute transport in aquifers or water flow in heterogeneous porous media. In this model, the time and space Hausdorff derivatives are defined on non-Euclidean fractal metrics with power law scaling transform which, respectively, connect the temporal and spatial complexity during transport. The Hausdorff derivative model can be transformed to an advection-dispersion equation with time- and space-dependent dispersion and convection coefficients. This model is a fractal partial differential equation (PDE) defined on a fractal space and differs from the fractional PDE which is derived for non-local transport of particles on a non-fractal Euclidean space. As an example of applications of this model, an explicit solution with a constant diffusion coefficient and flow velocity subject to an instantaneous source is derived and fitted to the breakthrough curves of tritium as a tracer in porous media. These results are compared with those of a scale-dependent dispersion model and a time-scale dependent dispersion model. Overall, it is found that the fractal PDE based on the Hausdorff derivatives better captures the early arrival and heavy tail in the scaled breakthrough curves for variable transport distances. The estimated parameters in the fractal Hausrdorff model represent clear physical mechanisms such as linear relationships between the orders of Hausdorff derivatives and the transport distance. The mathematical formulation is applicable to both solute transport and water flow in porous media.
Imbibition plays a central role in diverse energy, environmental, and industrial processes. In many cases, the medium has multiple parallel strata of different permeabilities; however, how this stratification impacts imbibition is poorly understood. We address this gap in knowledge by directly visualizing forced imbibition in three-dimensional (3D) porous media with two parallel strata. We find that imbibition is spatially heterogeneous: for small capillary number Ca, the wetting fluid preferentially invades the fine stratum, while for Ca above a threshold value, the fluid instead preferentially invades the coarse stratum. This threshold value depends on the medium geometry, the fluid properties, and the presence of residual wetting films in the pore space. These findings are well described by a linear stability analysis that incorporates crossflow between the strata. Thus, our work provides quantitative guidelines for predicting and controlling flow in stratified porous media.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
