Two level homogenization of flows in deforming double porosity media: biot-darcy-brinkman model

We present the two-level homogenization of the flow in a deformable double-porous structure described at two characteristic scales: the higher level porosity associated with the mesoscopic structure is constituted by channels in an elastic skeleton which is made of a microporous material. The macroscopic model is derived by the asymptotic analysis of the viscous flow in the heterogeneous structure characterized by two small parameters. The first level upscaling yields a Biot continuum model coupled with the Stokes flow. The second step of the homogenization leads to a macroscopic flow model which attains the form of the Darcy-Brinkman flow model coupled with the deformation of the poroelastic continuum involving the effective parameters given by the microscopic and the mesoscopic porosity features

Homogenization of the fluid-saturated piezoelectric porous metamaterials

The paper is devoted to the homogenization approach in modelling of peri- odic porous media constituted by piezoelectric porous skeleton with pores saturated by viscous fluid. The representative volume element contains the piezoelectric solid part (the matrix) and the fluid saturated pores (the channels). Both the matrix and the channels form connected subdomains. The mathematical model describing the material behaviour at the microscopic scale involves the quasi-static equilibrium equation governing the solid piezoelectric skeleton, the Stokes model of the viscous fluid flow in the channels and the coupling interface conditions on the transmission interface. The macroscopic model is derived using the unfolding method of homogenization. The effective material coefficients are computed using characteristic responses of the porous microstructure. The consti- tutive law for the upscaled piezo-poroelastic material involves a coefficient coupling the electric field and the pore pressure. A numerical example illustrates different responses of the porous medium subject to the drained and undrained loading

Flow rate--pressure drop relation for deformable shallow microfluidic channels

Laminar flow in devices fabricated from soft materials causes deformation of the passage geometry, which affects the flow rate--pressure drop relation. For a given pressure drop, in channels with narrow rectangular cross-section, the flow rate varies as the cube of the channel height, so deformation can produce significant quantitative effects, including nonlinear dependence on the pressure drop [{Gervais, T., El-Ali, J., G\"unther, A. \& Jensen, K.\ F.}\ 2006 Flow-induced deformation of shallow microfluidic channels.\ \textit{Lab Chip} \textbf{6}, 500--507]. Gervais et. al. proposed a successful model of the deformation-induced change in the flow rate by heuristically coupling a Hookean elastic response with the lubrication approximation for Stokes flow. However, their model contains a fitting parameter that must be found for each channel shape by performing an experiment. We present a perturbation approach for the flow rate--pressure drop relation in a shallow deformable microchannel using the theory of isotropic quasi-static plate bending and the Stokes equations under a lubrication approximation (specifically, the ratio of the channel's height to its width and of the channel's height to its length are both assumed small). Our result contains no free parameters and confirms Gervais et. al.'s observation that the flow rate is a quartic polynomial of the pressure drop. The derived flow rate--pressure drop relation compares favorably with experimental measurements.Comment: 20 pages, 6 figures; v2 minor revisions, accepted for publication in the Journal of Fluid Mechanic

Homogenization of coupled flow and deformation in a porous material

In this paper we present a framework for computational homogenization of the fluid-solid interaction that pertains to the coupled deformation and flow of pore fluid in a fluid-saturated porous material. Large deformations are considered and the resulting problem is established in the material setting. In order to ensure a proper FE-mesh in the fluid domain of the RVE, we introduce a fictitious elastic solid in the pores; however, the adopted variational setting ensures that the fictitious material does not obscure the motion of the (physical) solid skeleton. For the subsequent numerical evaluation of the RVE-response, hyperelastic properties are assigned to the solid material, whereas the fluid motion is modeled as incompressible Stokes\u27 flow. Variationally consistent homogenization of the standard first order is adopted. The homogenization is selective in the sense that the resulting macroscale (upscaled) porous media model reminds about the classical one for a quasi-static problem with displacements and pore pressure as the unknown macroscale fields. Hence, the (relative) fluid velocity, i.e. seepage, "lives" only on the subscale and is part of the set of unknown fields in the RVE-problem. As to boundary conditions on the RVE, a mixture of Dirichlet and weakly periodic conditions is adopted. In the numerical examples, special attention is given to an evaluation of the Biot coefficient that occurs in classical phenomenological models for porous media

Upscaling of the Coupling of Hydromechanical and Thermal Processes in a Quasi-static Poroelastic Medium

We undertake a formal derivation of a linear poro-thermo-elastic system within the framework of quasi-static deformation. This work is based upon the well-known derivation of the quasi-static poroelastic equations (also known as the Biot consolidation model) by homogenization of the fluid-structure interaction at the microscale. We now include energy, which is coupled to the fluid-structure model by using linear thermoelasticity, with the full system transformed to a Lagrangian coordinate system. The resulting upscaled system is similar to the linear poroelastic equations, but with an added conservation of energy equation, fully coupled to the momentum and mass conservation equations. In the end, we obtain a system of equations on the macroscale accounting for the effects of mechanical deformation, heat transfer, and fluid flow within a fully saturated porous material, wherein the coefficients can be explicitly defined in terms of the microstructure of the material. For the heat transfer we consider two different scaling regimes, one where the Péclet number is small, and another where it is unity.We also establish the symmetry and positivity for the homogenized coefficients.publishedVersio

Upscaling of the Coupling of Hydromechanical and Thermal Processes in a Quasi-static Poroelastic Medium

We undertake a formal derivation of a linear poro-thermo-elastic system within the framework of quasi-static deformation. This work is based upon the well-known derivation of the quasi-static poroelastic equations (also known as the Biot consolidation model) by homogenization of the fluid-structure interaction at the microscale. We now include energy, which is coupled to the fluid-structure model by using linear thermoelasticity, with the full system transformed to a Lagrangian coordinate system. The resulting upscaled system is similar to the linear poroelastic equations, but with an added conservation of energy equation, fully coupled to the momentum and mass conservation equations. In the end, we obtain a system of equations on the macroscale accounting for the effects of mechanical deformation, heat transfer, and fluid flow within a fully saturated porous material, wherein the coefficients can be explicitly defined in terms of the microstructure of the material. For the heat transfer we consider two different scaling regimes, one where the Péclet number is small, and another where it is unity.We also establish the symmetry and positivity for the homogenized coefficients.publishedVersio

Upscaling of the Coupling of Hydromechanical and Thermal Processes in a Quasi-static Poroelastic Medium

We undertake a formal derivation of a linear poro-thermo-elastic system within the framework of quasi-static deformation. This work is based upon the well-known derivation of the quasi-static poroelastic equations (also known as the Biot consolidation model) by homogenization of the fluid-structure interaction at the microscale. We now include energy, which is coupled to the fluid-structure model by using linear thermoelasticity, with the full system transformed to a Lagrangian coordinate system. The resulting upscaled system is similar to the linear poroelastic equations, but with an added conservation of energy equation, fully coupled to the momentum and mass conservation equations. In the end, we obtain a system of equations on the macroscale accounting for the effects of mechanical deformation, heat transfer, and fluid flow within a fully saturated porous material, wherein the coefficients can be explicitly defined in terms of the microstructure of the material. For the heat transfer we consider two different scaling regimes, one where the Péclet number is small, and another where it is unity.We also establish the symmetry and positivity for the homogenized coefficients.publishedVersio

Ion transport through deformable porous media: derivation of the macroscopic equations using upscaling

We study the upscaling or homogenization of the transport of a multicomponentelectrolyte in a dilute Newtonian solvent through a deformable porous medium.The pore scale interaction between the flow and the structure deformation is taken into account.After a careful adimensionalization process, we first consider so-called equilibrium solutions,in the absence of external forces, for which the velocity and diffusive fluxes vanish andthe electrostatic potential is the solution of a Poisson-Boltzmann equation.When the motion is governed by a small static electric field and small hydrodynamic and elastic forces,we use O'Brien's argument to deduce a linearized model. Then we perform the homogenizationof these linearized equations for a suitable choice of time scale. It turns out thatthe deformation of the porous medium is weakly coupled to the electrokinetics systemin the sense that it does not influence electrokinetics although the latter one yieldsan osmotic pressure term in the mechanical equations. As a byproduct we find that theeffective tensor satisfies Onsager properties, namely is symmetric positive definite

Solvability of a fluid-structure interaction problem with semigroup theory

Continuous semigroup theory is applied to proof the existence and uniqueness of a solution to a fluid-structure interaction (FSI) problem of non-stationary Stokes flow in two bulk domains, separated by a 2D elastic, permeable plate. The plate's curvature is proportional to the jump of fluid stresses across the plate and the flow resistance is modeled by Darcy's law. In the weak formulation of the considered physical problem, a linear operator in space is associated with a sum of two bilinear forms on the fluid and the interface domains, respectively. One attains a system of equations in operator form, corresponding to the weak problem formulation. Utilizing the sufficient conditions in the Lumer-Phillips theorem, we show that the linear operator is a generator of a contraction semigroup, and give the existence proof to the FSI problem