A hybrid numerical model for multiphase fluid flow in a deformable porous medium

In this paper, a fully coupled finite volume-finite element model for a deforming porous medium interacting with the flow of two immiscible pore fluids is presented. The basic equations describing the system are derived based on the averaging theory. Applying the standard Galerkin finite element method to solve this system of partial differential equations does not conserve mass locally. A non-conservative method may cause some accuracy and stability problems. The control volume based finite element technique that satisfies local mass conservation of the flow equations can be an appropriate alternative. Full coupling of control volume based finite element and the standard finite element techniques to solve the multiphase flow and geomechanical equilibrium equations is the main goal of this paper. The accuracy and efficiency of the method are verified by studying several examples for which analytical or numerical solutions are available. The effect of mesh orientation is investigated by simulating a benchmark water-flooding problem. A representative example is also presented to demonstrate the capability of the model to simulate the behavior in heterogeneous porous media.     

MIXED FINITE ELEMENT METHODS FOR FRACTIONAL NAVIER-STOKES EQUATIONS

This paper gives the detailed numerical analysis of mixed finite element method for fractional Navier-Stokes equations.The proposed method is based on the mixed finite element method in space and a finite difference scheme in time.The stability analyses of semi-discretization scheme and fully discrete scheme are discussed in detail.Furthermore,We give the convergence analysis for both semidiscrete and flly discrete schemes and then prove that the numerical solution converges the exact one with order O(h2+k),where h and k:respectively denote the space step size and the time step size.Finally,numerical examples are presented to demonstrate the effectiveness of our numerical methods.     

WEB‐spline based mesh‐free finite element analysis for the heat equation and the time‐dependent Navier–Stokes equation: A survey

In this article, we give some numerical techniques and error estimates using web‐spline based mesh‐free finite element method for the heat equation and the time‐dependent Navier–Stokes equations on bounded domains. The web‐spline method uses weighted extended B‐splines on a regular grid as basis functions and does not require any grid generation. We demonstrate the method by providing numerical results for the Poisson's and stationary Stokes equation. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

粘性流体两相多孔介质非线性动力问题的罚有限元法

采用基于混合物理论的多孔介质模型,给出粘性流体饱和两相多孔介质非线性动力问题的控制场方程以及相应边值和初值问题的提法,用Galerkin加权残值法导出罚有限元公式,并给出该非线性方程组的迭代求解方法。考虑了体积分数和渗透率与变形相关的情况。用编制的有限元程序计算分析了一维多孔柱体在脉冲载荷作用下的瞬态响应,数值结果表明文中方法正确有效。     

MONOLITHIC MULTIGRID FOR REDUCED MAGNETOHYDRODYNAMIC EQUATIONS

In this paper,the monolithic multigrid method is investigated for reduced magnetohydrodynamic equations.We propose a diagonal Braess-Sarazin smoother for the finite element discrete system and prove the uniform convergence of the MMG method with respect to mesh sizes.A multigrid-preconditioned FGMRES method is proposed to solve the magnetohydrodynamic equations.It turns out to be robust for relatively large physical parameters.By extensive numerical experiments,we demonstrate the optimality of the monolithic multigrid method with respect to the number of degrees of freedom.     

A second-order finite element variational multiscale scheme for the fully discrete unsteady Navier–Stokes equations

In this report, we present and study a fully discrete finite element variational multiscale scheme for the unsteady incompressible Navier–Stokes equations where high Reynolds numbers are allowed. The scheme uses conforming finite element pairs for spatial discretization and a three-point difference formula for temporal discretization which is of second-order, where a stabilization term based on two local Gauss integrations is employed to stabilize the numerical scheme. We prove stability of the scheme, derive a priori error estimates for the fully discrete solution, and finally, give some numerical results to verify the theoretical predictions and demonstrate the effectiveness of the proposed numerical scheme.     

On finite element schemes of the dirichlet problem

In this paper, we consider finite element schemes applied to the Dirichlet problem for the system of nonlinear elliptic equations, based on piecewise linear polynomials, and present iterative methods for solving algebraic nonlinear equations, which construct monotone sequences. Furthermore, we derive error estimates which imply uniform convergence. Our results are based on the discrete maximum principle. Finally, some typical numerical examples are given to demonstrate the usefulness of convergence results.     

基于参数变分原理的非均质材料弹塑性有限元分析的Voronoi单元法

在非均质材料的有限元数值模拟中,采用了Voronoi单元(VCFEM)以克服经典位移元的局限性．基于参数变分原理和二次规划法进行了Voronoi单元的二维弹塑性分析A·D2推导了有限元列式并形成最终的二次规划求解模型．研究了非均质材料微观夹杂对整体力学性能的影响．数值算例证明了该方法的正确和可行性．     

入水冲击问题变分原理及其它

首先建立入水前后两个衔接阶段的较为严密的场方程.再得到与之对应的各类变分原理,界限定理,第二阶段问题的边界积分方程.证明了解的存在性并提供了求解实施方案.最后以船舶兴波阻力问题的算例,论证了第二阶段问题的一种特殊应用及其正确性.从而为求取较为精确的入水冲击问题基本方程的变分有限元及边界元方法奠定了严密的理论基础.     

速度追踪问题中鞍点系统的新分裂迭代

 有限元离散一类速度追踪问题后得到具有鞍点结构的线性系统,针对该鞍点系统,本文提出了一种新的分裂迭代技术.证明了新的分裂迭代方法的无条件收敛性,详细分析了新的分裂预条件子对应的预处理矩阵的谱性质.数值结果验证了对于大范围的网格参数和正则参数,新的分裂预条件子在求解有限元离散速度追踪问题得到的鞍点系统时的可行性和有效性.     

Competing through altering the environment: A cross-diffusion population model coupled to transport–Darcy flow equations

We consider an evolution model describing the spatial population distribution of two salt tolerant plant species, such as mangroves, which are affected by inter- and intra-specific competition (Lotka–Volterra), population pressure (cross-diffusion) and environmental heterogeneity (environmental potential). The environmental potential and the Lotka–Volterra terms are assumed to depend on the salt concentration in the root region, which may change as a result of mangroves’ ability to uptake fresh water and leave the salt of the solution behind, in the saturated porous medium. Consequently, partial differential equations modelling the population dynamics on the surface are coupled with Darcy–transport equations modelling the salt and pressure-velocity distribution in the subsurface. We prove the existence of weak solutions of the coupled problem and provide a numerical discretization based on a stabilized mixed finite element method, which we use to numerically demonstrate the behaviour of the system.     

A Three-Field Finite Element Method for Elliptic Partial Differential Equations Driven by Stochastic Loads

Abstract

This paper is concerned with the application of nonconforming finite element methods to stochastic partial differential equations. We present a mixed formulation of a three-field finite element method applied to an elliptic model problem involving stochastic loads. We then derive the exact form for the expected value and variance of the solution. Additionally, the rate of convergence for the stochastic error is presented. Finally, we demonstrate through numerical experiments that the method is robust and reliable.     

ERROR ESTIMATES FOR TWO-SCALE COMPOSITE FINITE ELEMENT APPROXIMATIONS OF NONLINEAR PARABOLIC EQUATIONS

We study spatially semidiscrete and fully discrete two-scale composite finite element method for approximations of the nonlinear parabolic equations with homogeneous Dirich-let boundary conditions in a convex polygonal domain in the plane.This new class of finite elements,which is called composite finite elements,was first introduced by Hackbusch and Sauter[Numer.Math.,75(1997),pp.447-472]for the approximation of partial differential equations on domains with complicated geometry.The aim of this paper is to introduce an efficient numerical method which gives a lower dimensional approach for solving par-tial differential equations by domain discretization method.The composite finite element method introduces two-scale grid for discretization of the domain,the coarse-scale and the fine-scale grid with the degrees of freedom lies on the coarse-scale grid only.While the fine-scale grid is used to resolve the Dirichlet boundary condition,the dimension of the finite element space depends only on the coarse-scale grid.As a consequence,the resulting linear system will have a fewer number of unknowns.A continuous,piecewise linear composite finite element space is employed for the space discretization whereas the time discretization is based on both the backward Euler and the Crank-Nicolson methods.We have derived the error estimates in the L∞(L2)-norm for both semidiscrete and fully discrete schemes.Moreover,numerical simulations show that the proposed method is an efficient method to provide a good approximate solution.     

A Second Order Unconditionally Convergent Finite Element Method for the Thermal Equation with Joule Heating Problem

In this paper, we study the finite element approximation for nonlinear thermal equation. Because the nonlinearity of the equation, our theoretical analysis is based on the error of temporal and spatial discretization. We consider a fully discrete second order backward difference formula based on a finite element method to approximate the temperature and electric potential, and establish optimal $L^2$error estimates for the fully discrete finite element solution without any restriction on the time-step size. The discrete solution is bounded in infinite norm. Finally, several numerical examples are presented to demonstrate the accuracy and efficiency of the proposed method.     

Preconditioners for higher order edge finite element discretizations of Maxwell's equations

In this paper,we are concerned with the fast solvers for higher order edge finite element discretizations of Maxwell's equations.We present the preconditioners for the first family and second family of higher order N′ed′elec element equations,respectively.By combining the stable decompositions of two kinds of edge finite element spaces with the abstract theory of auxiliary space preconditioning,we prove that the corresponding condition numbers of our preconditioners are uniformly bounded on quasi-uniform grids.We also present some numerical experiments to demonstrate the theoretical results.     

The coupling system of Navier–Stokes equations and elastic Navier–Lame equations in a blood vessel

In this paper, the blood flow problem is considered in a blood vessel, and a coupling system of Navier–Stokes equations and linear elastic equations, Navier–Lame equations, in a cylinder with cylindrical elastic shell is given as the governing equations of the problem. We provide two finite element models to simulating the three-dimensional Navier–Stokes equations in the cylinder while the asymptotic expansion method is used to solving the linearly elastic shell equations. Specifically, in order to discrete the Navier–Stokes equations, the dimensional splitting strategy is constructed under the cylinder coordinate system. The spectral method is adopted along the rotation direction while the finite element method is used along the other directions. By using the above strategy, we get a series of two-dimensional-three-components (2D-3C) fluid problems. By introduce the S-coordinate system in E³ and employ the thickness of blood vessel wall as the expanding parameter, the asymptotic expansion method can be established to approximate the solution of the 3D elastic problem. The interface contact conditions can be treated exactly based on the knowledge of tensor analysis. Finally, numerical test shows that our method is reasonable.     

Meshfree Finite Volume Element Method for Constrained Optimal Control Problem Governed by Random Convection Diffusion Equations

In this paper, we investigate a stochastic meshfree finite volume element method for an optimal control problem governed by the convection diffusion equations with random coefficients. There are two contributions of this paper. Firstly, we establish a scheme to approximate the optimality system by using the finite volume element method in the physical space and the meshfree method in the probability space, which is competitive for high-dimensional random inputs. Secondly, the a priori error estimates are derived for the state,the co-state and the control variables. Some numerical tests are carried out to confirm the theoretical results and demonstrate the efficiency of the proposed method.     

求解Brinkman方程的修正弱有限元方法

本文针对Brinkman方程引入了一种修正弱Galerkin(MWG)有限元方法.我们通过具有两个离散弱梯度算子的变分形式来逼近模型. 在MWG方法中, 分别用次数为$k$和$k-1$的不连续分段多项式来近似速度函数$u$和压力函数$p$. MWG方法的主要思想是用内部函数的平均值代替边界函数. 因此, 与WG方法相比, MWG方法在不降低准确性的同时, 具有更少的自由度, 对于任意次数不超过$k-1$ 的多项式,MWG方法均可以满足稳定性条件. MWG 方法具有高度的灵活性, 它允许在具有一定形状正则性的任意多边形或多面体上使用不连续函数. 针对$H^1$和$L^22$范数下的速度和压力近似解, 建立了最优阶误差估计. 数值算例表明了该方法的准确性, 收敛性和稳定性.     

A finite element method for granular flow through a frictional boundary

A finite element method for the flow of dry granular solids through a domain involving a frictional contact boundary is formulated. The granular material is assumed as a compressible viscous-elastic–plastic continuum. Based on the principles of continuum mechanics, a complete set of equations is developed. The resulting boundary value problem is solved by the finite element method in space and by the finite difference method in time. The derivation of the finite element equations and the mathematical framework of the numerical technique are presented, together with two illustrative examples to demonstrate the validity of the technique.     

Exponential decay rate of a nonlinear suspension bridge model by a local distributed and boundary dampings

In this paper, we investigate the decay properties of the unconstrained one dimensional suspension bridge model. With only partial damping acting on one or on both equations and with boundary dampings, we prove that the first order energy is decaying exponentially, our method of proof is based on the energy method to build the appropriate Lyapunov functional. Moreover, we develop a numerical algorithm which is based on the finite element method to approximate the spatial variable and the Crank–Nicolson type of symmetric difference scheme to discretize the time derivative, and also a Picard type iteration process for solving the system of nonlinear equations obtained by discretization. At the end, we present some numerical experiments to illustrate our theoretical results.