No Arabic abstract
Many applications of porous media research involves high pressures and, correspondingly, exchange of thermal energy between the fluid and the matrix. While the system is relatively well understood for the case of non-moving porous media, the situation when the elastic matrix can move and deform, is much more complex. In this paper we derive the equations of motion for the dynamics of a deformable porous media which includes the effects of friction forces, stresses, and heat exchanges between the media, by using the new methodology of variational approach to thermodynamics. This theory extends the recently developed variational derivation of the mechanics of deformable porous media to include thermodynamic processes and can easily include incompressibility constraints. The model for the combined fluid-matrix system, written in the spatial frame, is developed by introducing mechanical and additional variables describing the thermal energy part of the system, writing the action principle for the system, and using a nonlinear, nonholonomic constraint on the system deduced from the second law of thermodynamics. The resulting equations give us the general version of possible friction forces incorporating thermodynamics, Darcy-like forces and friction forces similar to those used in the Navier-Stokes equations. The equations of motion are valid for arbitrary dependence of the kinetic and potential energies on the state variables. The results of our work are relevant for geophysical applications, industrial applications involving high pressures and temperatures, food processing industry, and other situations when both thermodynamics and mechanical considerations are important.
We derive the equations of motion for the dynamics of a porous media filled with an incompressible fluid. We use a variational approach with a Lagrangian written as the sum of terms representing the kinetic and potential energy of the elastic matrix, and the kinetic energy of the fluid, coupled through the constraint of incompressibility. As an illustration of the method, the equations of motion for both the elastic matrix and the fluid are derived in the spatial (Eulerian) frame. Such an approach is of relevance e.g. for biological problems, such as sponges in water, where the elastic porous media is highly flexible and the motion of the fluid has a primary role in the motion of the whole system. We then analyze the linearized equations of motion describing the propagation of waves through the media. In particular, we derive the propagation of S-waves and P-waves in an isotropic media. We also analyze the stability criteria for the wave equations and show that they are equivalent to the physicality conditions of the elastic matrix. Finally, we show that the celebrated Biots equations for waves in porous media are obtained for certain values of parameters in our models.
Many parts of biological organisms are comprised of deformable porous media. The biological media is both pliable enough to deform in response to an outside force and can deform by itself using the work of an embedded muscle. For example, the recent work (Ludeman et al., 2014) has demonstrated interesting sneezing dynamics of a freshwater sponge, when the sponge contracts and expands to clear itself from surrounding polluted water. We derive the equations of motion for the dynamics of such an active porous media (i.e., a deformable porous media that is capable of applying a force to itself with internal muscles), filled with an incompressible fluid. These equations of motion extend the earlier derived equation for a passive porous media filled with an incompressible fluid. We use a variational approach with a Lagrangian written as the sum of terms representing the kinetic and potential energy of the elastic matrix, and the kinetic energy of the fluid, coupled through the constraint of incompressibility. We then proceed to extend this theory by computing the case when both the active porous media and the fluid are incompressible, with the porous media still being deformable, which is often the case for biological applications. For the particular case of a uniform initial state, we rewrite the equations of motion in terms of two coupled telegraph-like equations for the material (Lagrangian) particles expressed in the Eulerian frame of reference, particularly suitable for numerical simulations, formulated for both the compressible media/incompressible fluid case and the doubly incompressible case. We derive interesting conservation laws for the motion, perform numerical simulations in both cases and show the possibility of self-propulsion of a biological organism due to particular running wave-like application of the muscle stress.
The flow of viscoelastic fluids in porous media is encountered in many practical applications, such as in the enhanced oil recovery process or in the groundwater remediation. Once the flow rate exceeds a critical value in such flows, an elastic instability with fluctuating flow field is observed, which ultimately transits to a more chaotic and turbulence-like flow structure as the flow rate further increases. In a recent study, it has been experimentally shown that this chaotic flow behaviour of viscoelastic fluids can be suppressed by increasing the geometric disorder in a model porous media consisting of a microchannel with several micropillars placed in it. However, the present numerical study demonstrates that this is not always true. We show that it depends on the initial arrangement of the micropillars for mimicking the porous media. In particular, we find that for an initial ordered and aligned configuration of the micropillars, the introduction of geometric order actually increases the chaotic flow dynamics as opposed to that seen for an initial ordered and staggered configuration of the micropillars. We suggest that this chaotic flow behaviour actually depends on the number of the stagnation points revealed to the flow field where maximum stretching of the viscoelastic microstructure happens. Our findings and explanation are perfectly in line with that observed and provided in a more recent experimental study.
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.
The Brinkman equation has found great popularity in modelling the interfacial flow between free fluid and a porous medium. However, it is still unclear how to determine an appropriate effective Brinkman viscosity without resolving the flow at the pore scale. Here, we propose to determine the Brinkman viscosity for rough porous media from the interface slip length and the interior permeability. Both slip and permeability can be determined from unit-cell analysis, thus enabling an a priori estimate of the effective viscosity. By comparing the velocity distribution in the porous material predicted from the Brinkman equation with that obtained from pore-scale resolved simulations, we show that modelling errors are $sim 10%$ and not larger than $40%$. We highlight the physical origins of the obtained errors and demonstrate that the Brinkman model can be much more accurate for irregular porous structures.