Stable cell-centered finite volume discretization for Biot equations

In this paper we discuss a new discretization for the Biot equations. The discretization treats the coupled system of deformation and flow directly, as opposed to combining discretizations for the two separate sub-problems. The coupled discretization has the following key properties, the combination of which is novel: 1) The variables for the pressure and displacement are co-located, and are as sparse as possible (e.g. one displacement vector and one scalar pressure per cell center). 2) With locally computable restrictions on grid types, the discretization is stable with respect to the limits of incompressible fluid and small time-steps. 3) No artificial stabilization term has been introduced. Furthermore, due to the finite volume structure embedded in the discretization, explicit local expressions for both momentum-balancing forces as well as mass-conservative fluid fluxes are available. We prove stability of the proposed method with respect to all relevant limits. Together with consistency, this proves convergence of the method. Finally, we give numerical examples verifying both the analysis and convergence of the method

Review of Summation-by-parts schemes for initial-boundary-value problems

High-order finite difference methods are efficient, easy to program, scales well in multiple dimensions and can be modified locally for various reasons (such as shock treatment for example). The main drawback have been the complicated and sometimes even mysterious stability treatment at boundaries and interfaces required for a stable scheme. The research on summation-by-parts operators and weak boundary conditions during the last 20 years have removed this drawback and now reached a mature state. It is now possible to construct stable and high order accurate multi-block finite difference schemes in a systematic building-block-like manner. In this paper we will review this development, point out the main contributions and speculate about the next lines of research in this area

A fast immersed boundary method for external incompressible viscous flows using lattice Green's functions

A new parallel, computationally efficient immersed boundary method for solving three-dimensional, viscous, incompressible flows on unbounded domains is presented. Immersed surfaces with prescribed motions are generated using the interpolation and regularization operators obtained from the discrete delta function approach of the original (Peskin's) immersed boundary method. Unlike Peskin's method, boundary forces are regarded as Lagrange multipliers that are used to satisfy the no-slip condition. The incompressible Navier-Stokes equations are discretized on an unbounded staggered Cartesian grid and are solved in a finite number of operations using lattice Green's function techniques. These techniques are used to automatically enforce the natural free-space boundary conditions and to implement a novel block-wise adaptive grid that significantly reduces the run-time cost of solutions by limiting operations to grid cells in the immediate vicinity and near-wake region of the immersed surface. These techniques also enable the construction of practical discrete viscous integrating factors that are used in combination with specialized half-explicit Runge-Kutta schemes to accurately and efficiently solve the differential algebraic equations describing the discrete momentum equation, incompressibility constraint, and no-slip constraint. Linear systems of equations resulting from the time integration scheme are efficiently solved using an approximation-free nested projection technique. The algebraic properties of the discrete operators are used to reduce projection steps to simple discrete elliptic problems, e.g. discrete Poisson problems, that are compatible with recent parallel fast multipole methods for difference equations. Numerical experiments on low-aspect-ratio flat plates and spheres at Reynolds numbers up to 3,700 are used to verify the accuracy and physical fidelity of the formulation.Comment: 32 pages, 9 figures; preprint submitted to Journal of Computational Physic