On the Cahn-Hilliard-Brinkman system

We consider a diffuse interface model for phase separation of an isothermal incompressible binary fluid in a Brinkman porous medium. The coupled system consists of a convective Cahn-Hilliard equation for the phase field $\phi$, i.e., the difference of the (relative) concentrations of the two phases, coupled with a modified Darcy equation proposed by H.C. Brinkman in 1947 for the fluid velocity $\mathbf{u}$. This equation incorporates a diffuse interface surface force proportional to $\phi \nabla \mu$, where $\mu$ is the so-called chemical potential. We analyze the well-posedness of the resulting Cahn-Hilliard-Brinkman (CHB) system for $(\phi,\mathbf{u})$. Then we establish the existence of a global attractor and the convergence of a given (weak) solution to a single equilibrium via {\L}ojasiewicz-Simon inequality. Furthermore, we study the behavior of the solutions as the viscosity goes to zero, that is, when the CHB system approaches the Cahn-Hilliard-Hele-Shaw (CHHS) system. We first prove the existence of a weak solution to the CHHS system as limit of CHB solutions. Then, in dimension two, we estimate the difference of the solutions to CHB and CHHS systems in terms of the viscosity constant appearing in CHB

Analysis of the Brinkman-Forchheimer equations with slip boundary conditions

In this work, we study the Brinkman-Forchheimer equations driven under slip boundary conditions of friction type. We prove the existence and uniqueness of weak solutions by means of regularization combined with the Faedo-Galerkin approach. Next we discuss the continuity of the solution with respect to Brinkman's and Forchheimer's coefficients. Finally, we show that the weak solution of the corresponding stationary problem is stable

Incorporating variable viscosity in vorticity-based formulations for Brinkman equations

In this brief note, we introduce a non-symmetric mixed finite element formulation for Brinkman equations written in terms of velocity, vorticity and pressure with non-constant viscosity. The analysis is performed by the classical Babu\v{s}ka-Brezzi theory, and we state that any inf-sup stable finite element pair for Stokes approximating velocity and pressure can be coupled with a generic discrete space of arbitrary order for the vorticity. We establish optimal a priori error estimates which are further confirmed through computational example

Analysis of the Brinkman-Forchheimer equations with slip boundary conditions

In this work, we study the Brinkman-Forchheimer equations driven under slip boundary conditions of friction type. We prove the existence and uniqueness of weak solutions by means of regularization combined with the Faedo-Galerkin approach. Next we discuss the continuity of the solution with respect to Brinkman's and Forchheimer's coefficients. Finally, we show that the weak solution of the corresponding stationary problem is stable

Boundary elements method for microfluidic two-phase flows in shallow channels

In the following work we apply the boundary element method to two-phase flows in shallow microchannels, where one phase is dispersed and does not wet the channel walls. These kinds of flows are often encountered in microfluidic Lab-on-a-Chip devices and characterized by low Reynolds and low capillary numbers. Assuming that these channels are homogeneous in height and have a large aspect ratio, we use depth-averaged equations to describe these two-phase flows using the Brinkman equation, which constitutes a refinement of Darcy's law. These partial differential equations are discretized and solved numerically using the boundary element method, where a stabilization scheme is applied to the surface tension terms, allowing for a less restrictive time step at low capillary numbers. The convergence of the numerical algorithm is checked against a static analytical solution and on a dynamic test case. Finally the algorithm is applied to the non-linear development of the Saffman-Taylor instability and compared to experimental studies of droplet deformation in expanding flows.Comment: accepted for publication, Computers and Fluids 201

On the nonlocal Cahn-Hilliard-Brinkman and Cahn-Hilliard-Hele-Shaw systems

The phase separation of an isothermal incompressible binary fluid in a porous medium can be described by the so-called Brinkman equation coupled with a convective Cahn-Hilliard (CH) equation. The former governs the average fluid velocity $\mathbf{u}$, while the latter rules evolution of $\varphi$, the difference of the (relative) concentrations of the two phases. The two equations are known as the Cahn-Hilliard-Brinkman (CHB) system. In particular, the Brinkman equation is a Stokes-like equation with a forcing term (Korteweg force) which is proportional to $\mu\nabla\varphi$, where $\mu$ is the chemical potential. When the viscosity vanishes, then the system becomes the Cahn-Hilliard-Hele-Shaw (CHHS) system. Both systems have been studied from the theoretical and the numerical viewpoints. However, theoretical results on the CHHS system are still rather incomplete. For instance, uniqueness of weak solutions is unknown even in 2D. Here we replace the usual CH equation with its physically more relevant nonlocal version. This choice allows us to prove more about the corresponding nonlocal CHHS system. More precisely, we first study well-posedness for the CHB system, endowed with no-slip and no-flux boundary conditions. Then, existence of a weak solution to the CHHS system is obtained as a limit of solutions to the CHB system. Stronger assumptions on the initial datum allow us to prove uniqueness for the CHHS system. Further regularity properties are obtained by assuming additional, though reasonable, assumptions on the interaction kernel. By exploiting these properties, we provide an estimate for the difference between the solution to the CHB system and the one to the CHHS system with respect to viscosity