In this paper, we show that the permeability of a porous material (Tartar (1980)) and that of a bubbly fluid (Lipton and Avellaneda. Proc. R. Soc. Edinburgh Sect. A: Math. 114 (1990), 71–79) are limiting cases of the complexified version of the two-fluid models posed in Lipton and Avellaneda (Proc. R. Soc. Edinburgh Sect. A: Math. 114 (1990), 71–79). We assume the viscosity of the inclusion fluid is $z\mu _1$ and the viscosity of the hosting fluid is $\mu _1\in \mathbb {R}^{+}$, $z\in \mathbb {C}$. The proof is carried out by the construction of solutions for large $|z|$ and small $|z|$ with an iteration process similar to the one used in Bruno and Leo (Arch. Ration. Mech. Anal. 121 (1993), 303–338) and Golden and Papanicolaou (Commun. Math. Phys. 90 (1983), 473–491) and the analytic continuation. Moreover, we also show that for a fixed microstructure, the permeabilities of these three cases share the same integral representation formula (3.17) with different values of contrast parameter $s:=1/(z-1)$, as long as $s$ is outside the interval $\left [-\frac {2E_2^{2}}{1+2E_2^{2}},-\frac {1}{1+2E_1^{2}}\right ]$, where the positive constants $E_1$ and $E_2$ are the extension constants that depend only on the geometry of the periodic pore space of the material.

In this paper, we consider the Stokes equations in a perforated domain. When the number of holes increases while their radius tends to 0, it is proven in Desvillettes et al. [J. Stat. Phys. 131 (2008) 941–967], under suitable dilution assumptions, that the solution is well approximated asymptotically by solving a Stokes–Brinkman equation. We provide here quantitative estimates in $L^{p}$-norms of this convergence.

In this paper, we study an exponential time differencing method for solving the gauge system of incompressible viscous flows governed by Stokes or Navier-Stokes equations. The momentum equation is decoupled from the kinematic equation at a discrete level and is then solved by exponential time stepping multistep schemes in our approach. We analyze the stability of the proposed method and rigorously prove that the first order exponential time differencing scheme is unconditionally stable for the Stokes problem. We also present a compact representation of the algorithm for problems on rectangular domains, which makes FFT-based solvers available for the resulting fully discretized system. Various numerical experiments in two and three dimensional spaces are carried out to demonstrate the accuracy and stability of the proposed method.

A full multigrid method with coarsening by a factor-of-three to distributed control problems constrained by Stokes equations is presented. An optimal control problem with cost functional of velocity and/or pressure tracking-type is considered with Dirichlet boundary conditions. The optimality system that results from a Lagrange multiplier framework, form a linear system connecting the state, adjoint, and control variables. We investigate multigrid methods with finite difference discretization on staggered grids. A coarsening by a factor-of-three is used on staggered grids that results nested hierarchy of staggered grids and simplified the inter-grid transfer operators. A distributive-Gauss-Seidel smoothing scheme is employed to update the state- and adjoint-variables and a gradient update step is used to update the control variables. Numerical experiments are presented to demonstrate the effectiveness and efficiency of the proposed multigrid framework to tracking-type optimal control problems.

In this paper, we consider an optimal control problem governed by Stokes equations with H1-norm state constraint. The control problem is approximated by spectral method, which provides very accurate approximation with a relatively small number of unknowns. Choosing appropriate basis functions leads to discrete system with sparse matrices. We first present the optimality conditions of the exact and the discrete optimal control systems, then derive both a priori and a posteriori error estimates. Finally, an illustrative numerical experiment indicates that the proposed method is competitive, and the estimator can indicate the errors very well.

A nonconforming rectangular finite element method is proposed to solve a fluid structure interaction problem characterized by the Darcy-Stokes-Brinkman Equation with discontinuous coefficients across the interface of different structures. A uniformly stable mixed finite element together with Nitsche-type matching conditions that automatically adapt to the coupling of different sub-problem combinations are utilized in the discrete algorithm. Compared with other finite element methods in the literature, the new method has some distinguished advantages and features. The Boland-Nicolaides trick is used in proving the inf-sup condition for the multidomain discrete problem. Optimal error estimates are derived for the coupled problem by analyzing the approximation errors and the consistency errors. Numerical examples are also provided to confirm the theoretical results.

In this work, we study numerical methods for a coupled fluid-porous media flow model. The model consists of Stokes equations and Darcy's equations in two neighboring subdomains, coupling together through certain interface conditions. The weak form for the coupled model is of saddle point type. A mortar finite element method is proposed to approximate the weak form of the coupled problem. In our method, nonconforming Crouzeix-Raviart elements are applied in the fluid subdomain and the lowest order Raviart-Thomas elements are applied in the porous media subdomain; Meshes in different subdomains are allowed to be nonmatching on the common interface; Interface conditions are weakly imposed via adding constraint in the definition of the finite element space. The well-posedness of the discrete problem and the optimal error estimate for the proposed method are established. Numerical experiments are also given to confirm the theoretical results.

The goal of this work is to construct and study hybrid and multiplicative two-level overlapping Schwarz algorithms with standard coarse spaces for the almost incompressible linear elasticity and Stokes systems, discretized by mixed finite and spectral element methods with discontinuous pressures. Two different approaches are considered to solve the resulting saddle point systems: a) a preconditioned conjugate gradient (PCG) method applied to the symmetric positive definite reformulation of the almost incompressible linear elasticity system obtained by eliminating the pressure unknowns; b) a GMRES method with indefinite overlapping Schwarz preconditioner applied directly to the saddle point formulation of both the elasticity and Stokes systems. Condition number estimates and convergence properties of the proposed hybrid and multiplicative overlapping Schwarz algorithms are proven for the positive definite reformulation of almost incompressible elasticity. These results are based on our previous study [8] where only additive Schwarz preconditioners were considered for almost incompressible elasticity. Extensive numerical experiments with both finite and spectral elements show that the proposed overlapping Schwarz preconditioners are scalable, quasi-optimal in the number of unknowns across individual subdomains and robust with respect to discontinuities of the material parameters across subdomains interfaces. The results indicate that the proposed preconditioners retain a good performance also when the quasi-monotonicity assumption, required by the available theory, does not hold.

Within the framework of variational modelling we derive a one-phase moving boundary problem describing the motion of a semipermeable membrane enclosing a viscous liquid, driven by osmotic pressure and surface tension of the membrane. For this problem we prove the existence of classical solutions for a short-time.

Two-phase fluid flows on substrates (i.e. wetting phenomena) are important in many industrial processes, such as micro-fluidics and coating flows. These flows include additional physical effects that occur near moving (three-phase) contact lines. We present a new 2-D variational (saddle-point) formulation of a Stokesian fluid with surface tension that interacts with a rigid substrate. The model is derived by an Onsager type principle using shape differential calculus (at the sharp-interface, front-tracking level) and allows for moving contact lines and contact angle hysteresis and pinning through a variational inequality. Moreover, the formulation can be extended to include non-linear contact line motion models. We prove the well-posedness of the time semi-discrete system and fully discrete method using appropriate choices of finite element spaces. A formal energy law is derived for the semi-discrete and fully discrete formulations and preliminary error estimates are also given. Simulation results are presented for a droplet in multiple configurations to illustrate the method.

A straightforward method is presented for computing three-dimensional Stokes flow, due to forces on a surface, with high accuracy at points near the surface. The flow quantities are written as boundary integrals using the free-space Green’s function. To evaluate the integrals near the boundary, the singular kernels are regularized and a simple quadrature is applied in coordinate charts. High order accuracy is obtained by adding special corrections for the regularization and discretization errors, derived here using local asymptotic analysis. Numerical tests demonstrate the uniform convergence rates of the method.

We present an unconditionally energy stable and uniquely solvable finite difference scheme for the Cahn-Hilliard-Brinkman (CHB) system, which is comprised of a Cahn-Hilliard-type diffusion equation and a generalized Brinkman equation mod-eling fluid flow. The CHB system is a generalization of the Cahn-Hilliard-Stokes model and describes two phase very viscous flows in porous media. The scheme is based on a convex splitting of the discrete CH energy and is semi-implicit. The equations at the implicit time level are nonlinear, but we prove that they represent the gradient of a strictly convex functional and are therefore uniquely solvable, regardless of time step size. Owing to energy stability, we show that the scheme is stable in the time and space discrete and norms. We also present an efficient, practical nonlinear multigrid method . comprised of a standard FAS method for the Cahn-Hilliard part, and a method based on the Vanka smoothing strategy for the Brinkman part . for solving these equations. In particular, we provide evidence that the solver has nearly optimal complexity in typical situations. The solver is applied to simulate spinodal decomposition of a viscous fluid in a porous medium, as well as to the more general problems of buoyancy- and boundary-driven flows.

This paper aims at a general guideline to obtain a posteriori error estimates for the finite element error control in computational partial differential equations. In the abstract setting of mixed formulations, a generalised formulation of the corresponding residuals is proposed which then allows for the unified estimation of the respective dual norms. Notably, this can be done with an approach which is applicable in the same way to conforming, nonconforming and mixed discretisations. Subsequently, the unified approach is applied to various model problems. In particular, we consider the Laplace, Stokes, Navier-Lamé, and the semi-discrete eddy current equations.

We consider the effect of surface roughness on solid-solid contact in a Stokes flow.Various models for the roughness are considered, and a unified methodology is given toderive the corresponding asymptotics of the drag force in the close-contact limit. In thisway, we recover and clarify the various expressions that can be found in previousstudies.

In this paper, we study the numerical solution of the Stokes system in deformed axisymmetric geometries. In the azimuthal direction the discretization is carried out by using truncated Fourier series, thus reducing the dimension of the problem. The resulting two-dimensional problems are discretized using the spectral element method which is based on the variational formulation in primitive variables. The meridian domain is subdivided into elements, in each of which the solution is approximated by truncated polynomial series. The results of numerical experiments for several geometries are presented.

We describe possible solutions for a stationary flow of two superposed fluids between two close surfaces in relative motion. Physically, this study is within the lubrication framework, in which it is of interest to predict the relative positions of the lubricant and the air in the device. Mathematically, we observe that this problem corresponds to finding the interface between the two fluids, and we prove that this interface can be viewed as a square root of a polynomial of degree at most 6. We solve this equation using an original method. First, we check that our results are consistent with previous work. Next, we use this solution to answer some physically relevant questions related to the lubrication setting. For instance, we obtain theoretical and numerical results, which can predict the occurrence of a full film with respect to physical parameters (fluxes, shear velocity, viscosities). In particular, we present a figure giving the number of stationary solutions depending on the physical parameters. Moreover, we give some indications for a better understanding of the multi-fluid case.

A triangular spectral method for the Stokes equations is developed in this paper. The main contributions are two-fold: First of all, a spectral method using the rational approximation is constructed and analyzed for the Stokes equations in a triangular domain. The existence and uniqueness of the solution, together with an error estimate for the velocity, are proved. Secondly, a nodal basis is constructed for the efficient implementation of the method. These new basis functions enjoy the fully tensorial product property as in a tensor-produce domain. The new triangular spectral method makes it easy to treat more complex geometries in the classical spectral-element framework, allowing us to use arbitrary triangular and tetrahedral elements.

In this paper, new finite difference methods based on the augmented immersed interface method (IIM) are proposed for simulating an inextensible moving interface in an incompressible two-dimensional flow. The mathematical models arise from studying the deformation of red blood cells in mathematical biology. The governing equations are incompressible Stokes or Navier-Stokes equations with an unknown surface tension, which should be determined in such a way that the surface divergence of the velocity is zero along the interface. Thus, the area enclosed by the interface and the total length of the interface should be conserved during the evolution process. Because of the nonlinear and coupling nature of the problem, direct discretization by applying the immersed boundary or immersed interface method yields complex nonlinear systems to be solved. In our new methods, we treat the unknown surface tension as an augmented variable so that the augmented IIM can be applied. Since finding the unknown surface tension is essentially an inverse problem that is sensitive to perturbations, our regularization strategy is to introduce a controlled tangential force along the interface, which leads to a least squares problem. For Stokes equations, the forward solver at one time level involves solving three Poisson equations with an interface. For Navier-Stokes equations, we propose a modified projection method that can enforce the pressure jump condition corresponding directly to the unknown surface tension. Several numerical experiments show good agreement with other results in the literature and reveal some interesting phenomena.

This study is mainly dedicated to the development and analysis ofnon-overlapping domain decomposition methods for solving continuous-pressurefinite element formulations of the Stokes problem. These methods have thefollowing special features. By keeping the equations and unknowns unchanged atthe cross points, that is, points shared by more than two subdomains, one caninterpret them as iterative solvers of the actual discrete problem directlyissued from the finite element scheme. In this way, the good stabilityproperties of continuous-pressure mixed finite element approximations of theStokes system are preserved. Estimates ensuring that each iteration can beperformed in a stable way as well as a proof of the convergence of theiterative process provide a theoretical background for the application of therelated solving procedure. Finally some numerical experiments are given todemonstrate the effectiveness of the approach, and particularly to compare itsefficiency with an adaptation to this framework of a standard FETI-DP method.

The gradient based topological optimization tools introduced during thelast ten years tend naturally to modify the topology of a domain by creating small holes inside the domain. Once these holes have been created, they usually remainunchanged, at least during the topological phase of the optimizationalgorithm. In this paper, a new asymptotic expansion is introduced which allows to decide whether an existing hole must be removed or not forimproving the cost function. Then, two numerical examples are presented:the first one compares topological optimization with standard shape optimization, and the second one, issued from a lake oxygenation problem, illustrates the use of the new asymptotic expansion.