Superconvergence for a mixed finite elmenent method for elastic wave propagation in a plane domain

Summary A higher order mixed finite element method is introduced to approximate the solution of wave propagation in a plane elastic medium. A quasi-projection analysis is given to obtain error estimates in Sobolev spaces of nonpositive index. Estimates are given for difference quotients for a spatially periodic problem and superconvergence results of the same type as those of Bramble and Schatz for Galerkin methods are derived.     

A nonoverlapping domain decomposition method for stabilized finite element approximations of the Oseen equations

A nonoverlapping domain decomposition algorithm of Robin–Robin type is applied to the discretized Oseen equations using stabilized finite element approximations of velocity and pressure thus allowing in particular equal-order interpolation. As a crucial result we have to inspect the proof of a modified inf–sup condition, in particular, the dependence of the stability constant with respect to the Reynolds number (cf. appendix). After proving coercivity and strong convergence of the method, we derive an a posteriori estimate which controls convergence of the discrete subdomain solutions to the global discrete solution provided that jumps of the discrete solution converge at the interface. Furthermore, we obtain information on the design of some free parameters within the Robin-type interface condition which essentially influence the convergence speed. Some numerical results confirm the theoretical ones.     

A symmetric mixed finite element method for nearly incompressible elasticity based on biorthogonal systems

We present a symmetric version of the nonsymmetric mixed finite element method presented in (Lamichhane, ANZIAM J 50 (2008), C324–C338) for nearly incompressible elasticity. The displacement–pressure formulation of linear elasticity is discretized using a Petrov–Galerkin discretization for the pressure equation in (Lamichhane, ANZIAM J 50 (2008), C324–C338) leading to a non‐symmetric saddle point problem. A new three‐field formulation is introduced to obtain a symmetric saddle point problem which allows us to use a biorthogonal system. Working with a biorthogonal system, we can statically condense out all auxiliary variables from the saddle point problem arriving at a symmetric and positive‐definite system based only on the displacement. We also derive a residual based error estimator for the mixed formulation of the problem. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

Augmented hybrid finite element method for domain decomposition

We present an Augmented Hybrid Finite Element Method for Domain Decompositon. In this method, finite element approximations are defined independently on each subdomain and do not match at interface. This dows the user to mda local change of design, of meshes on one aubdomain without modifying other subdomains. Optimal reaults are obtained for a second-order model problem.     

Domain decomposition and splitting methods for Mortar mixed finite element approximations to parabolic equations

Summary We introduce in this article a new domain decomposition algorithm for parabolic problems that combines Mortar Mixed Finite Element methods for the space discretization with operator splitting schemes for the time discretization. The main advantage of this method is to be fully parallel. The algorithm is proven to be unconditionally stable and a convergence result in (Δt/h ^1/2) is presented.     

A mixed finite element method for the Stokes equations

A new mixed finite element for the Stokes equations is considered. This new finite element is based on a mixed formulation of the Stokes problem in which the gradient of the velocity is introduced and the velocity is approximated by the Raviart-Thomas element [1]. Optimal error estimates are derived. The number of degrees of freedom, for this element, is the lowest possible, and the local conservation of the mass is assured. A hybrid version of the mixed method is also considered. Finally, some numerical results for the incompressible Navier-Stokes equations are presented. © 1994 John Wiley & Sons, Inc.     

An analysis of a mixed finite element method for the Navier-Stokes equations

Summary A simple mixed finite element method is developed to solve the steady state, incompressible Navier-Stokes equations in a neighborhood of an isolated—but not necessarily unique—solution. Convergence is established under very mild restrictions on the triangulation, and, when the solution is sufficiently smooth, optimal error bounds are obtained.     

A mixed virtual element method for nearly incompressible linear elasticity equations

Over the last few years deep artificial neural networks (ANNs) have very successfully been used in numerical simulations for a wide variety of computational problems including computer vision, image classification, speech recognition, natural language processing, as well as computational advertisement. In addition, it has recently been proposed to approximate solutions of high-dimensional partial differential equations (PDEs) by means of stochastic learning problems involving deep ANNs. There are now also a few rigorous mathematical results in the scientific literature which provide error estimates for such deep learning based approximation methods for PDEs. All of these articles provide spatial error estimates for ANN approximations for PDEs but do not provide error estimates for the entire space-time error for the considered ANN approximations. It is the subject of the main result of this article to provide space-time error estimates for deep ANN approximations of Euler approximations of certain perturbed differential equations. Our proof of this result is based (i) on a certain ANN calculus and (ii) on ANN approximation results for products of the form \([0,T]\times \mathbb {R}^{d}\ni (t,x){\kern -.5pt}\mapsto {\kern -.5pt} tx{\kern -.5pt}\in {\kern -.5pt} \mathbb {R}^{d}\) where \(T{\kern -.5pt}\in {\kern -.5pt} (0,\infty )\), \(d{\kern -.5pt}\in {\kern -.5pt} \mathbb {N}\), which we both develop within this article.

     

B-spline based finite element method in one-dimensional discontinuous elastic wave propagation

The B-spline variant of the finite element method (FEM) is tested in one-dimensional discontinuous elastic wave propagation. The B-spline based FEM (called Isogeometric analysis IGA) uses spline functions as testing and shape functions in the Galerkin continuous content. Here, the accuracy of stress distribution and spurious oscillations of the B-spline based FEM are studied in numerical modeling of one-dimensional propagation of stress discontinuities in a bar, where the analytical solution is known. For time integration, the Newmark method, implicit form of the generalized-α method, the central difference method and the predictor/multi-corrector method are tested and compared. The use of lumped and consistent mass matrices in explicit time integration is discussed. Due to accuracy, the consistent mass matrix is preferred in explicit time integration in IGA.     

Limit analysis decomposition and finite element mixed method

This paper proposes an original decomposition approach to the upper bound method of limit analysis. It is based on a mixed finite element approach and on a convex interior point solver, using linear or quadratic discontinuous velocity fields. Presented in plane strain, this method appears to be rapidly convergent, as verified in the Tresca compressed bar problem in the linear velocity case. Then, using discontinuous quadratic velocity fields, the method is applied to the celebrated problem of the stability factor of a Tresca vertical slope: the upper bound is lowered to 3.7776-value to be compared to the best published lower bound 3.7752-by succeeding in solving a nonlinear optimization problem with millions of variables and constraints.     

Solution domain decomposition with finite difference methods for partial differential equations

The Solution Domain Decomposition method of Nguyen and Paik [J. Sci. Comput. 4 , 357 (1993)] originally developed in a pseudospectral context is extended for use with finite difference techniques to solve partial differential equations. The essential idea behind this method lies in an application of the superposition principle, which allows interactions between adjacent subdomains to be decoupled and the resulting equations to be solved in parallel. Several tests are performed to assess its accuracy and efficiency based on a model problem arising from thermal convection inside a fluid-saturated porous cavity with heating from a side. The results reveal a well adaptation of the methodology into the framework of finite differences. © 1995 John Wiley & Sons, Inc.     

Refined mixed finite element method for the elasticity problem in a polygonal domain

The purpose of this article is to study a mixed formulation of the elasticity problem in plane polygonal domains and its numerical approximation. In this mixed formulation the strain tensor is introduced as a new unknown and its symmetry is relaxed by a Lagrange multiplier, which is nothing else than the rotation. Because of the corner points, the displacement field is not regular in general in the vicinity of the vertices but belongs to some weighted Sobolev space. Using this information, appropriate refinement rules are imposed on the family of triangulations in order to recapture optimal error estimates. Moreover, uniform error estimates in the Lamé coefficient λ are obtained for λ large. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 323–339, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10009     

Least‐squares mixed finite element methods for the RLW equations

A least‐squares mixed finite element (LSMFE) schemes are formulated to solve the 1D regularized long wave (RLW) equations and the convergence is discussed. The L² error estimates of LSMFE methods for RLW equations under the standard regularity assumption on the finite element partition are given.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

The modified method of characteristics with mixed finite element domain decomposition procedures for the transient behavior of a semiconductor device

For the transient behavior of a semiconductor device, the modified method of characteristics with mixed finite element domain decomposition procedures applicable to parallel arithmetic is put forward. The electric potential equation is described by the mixed finite element method, and the electric, hole concentration and heat conduction equations are treated by the modified method of characteristics finite element domain decomposition methods. Some techniques, such as calculus of variations, domain decomposition, characteristic method, energy method, negative norm estimate and prior estimates and techniques are employed. Optimal order estimates in L² norm are derived for the error in the approximation solution. Thus the well‐known theoretical problem has been thoroughly and completely solved.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 353–368 2012     

Two‐Grid method for nonlinear parabolic equations by expanded mixed finite element methods

In this article, we develop a two‐grid algorithm for nonlinear reaction diffusion equation (with nonlinear compressibility coefficient) discretized by expanded mixed finite element method. The key point is to use two‐grid scheme to linearize the nonlinear term in the equations. The main procedure of the algorithm is solving a small‐scaled nonlinear equations on the coarse grid and dealing with a linearized system on the fine space using the Newton iteration with the coarse grid solution. Error estimation to the expanded mixed finite element solution is analyzed in detail. We also show that two‐grid solution achieves the same accuracy as long as the mesh sizes satisfy H = O(h^1/2). Two numerical experiments are given to verify the effectiveness of the algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A parallel finite element variational multiscale method based on fully overlapping domain decomposition for incompressible flows

Based on fully overlapping domain decomposition and a recent variational multiscale method, a parallel finite element variational multiscale method for convection dominated incompressible flows is proposed and analyzed. In this method, each processor computes a local finite element solution in its own subdomain using a global mesh that is locally refined around its own subdomain, where a stabilization term based on two local Gauss integrations is adopted to stabilize the numerical form of the Navier–Stokes equations. Using the technical tool of local a priori estimate for the finite element solution, error bounds of the discrete solution are estimated. Algorithmic parameter scalings are derived. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the method. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 856–875, 2015     

An iterative and parallel solver based on domain decomposition for the h-p version of the finite element method

In this paper, we propose a new interative and parallel solver, based on domain decomposition, for the h-p version of the finite element method in two dimensions. It improves our previous work in two aspects: (1) A subdomain may contain several super-elements of the coarse mesh, thus can be of arbitrary shape and size. This makes the solver more efficient and more flexible in computational practice. (2) The p-version components (i.e., the high order side and internal modes) in every element are treated separately, which results in better parallelism.     

Convergence of an adaptive mixed finite element method for convection-diffusion-reaction equations

We prove the convergence of an adaptive mixed finite element method (AMFEM) for (nonsymmetric) convection-diffusion-reaction equations. The convergence result holds for the cases where convection or reaction is not present in convection- or reaction-dominated problems. A novel technique of analysis is developed by using the superconvergence of the scalar displacement variable instead of the quasi-orthogonality for the stress and displacement variables, and without marking the oscillation dependent on discrete solutions and data. We show that AMFEM is a contraction of the error of the stress and displacement variables plus some quantity. Numerical experiments confirm the theoretical results.     

Two-grid stabilized mixed finite element method for fully discrete reaction-diffusion equations

Two-grid mixed finite element method is proposed based on backward Euler schemes for the unsteady reaction-diffusion equations. The scheme combines with the stabilized mixed finite element scheme by using the lowest equal-order pairs for the velocity and pressure. The space two-grid method is also used to reduce the time consuming. The benefits of this approach are to avoid the higher derivative, but to have more favorable stability, and to get the numerical solution of the two unknown variables simultaneously. Stability analysis and error estimates are given in this work. Finally, the theoretical results are verified by the numerical examples.