A fully discrete stabilized finite element method for the time-dependent Navier-Stokes equations

In this article, we consider a fully discrete stabilized finite element method based on two local Gauss integrations for the two-dimensional time-dependent Navier-Stokes equations. It focuses on the lowest equal-order velocity-pressure pairs. Unlike the other stabilized method, the present approach does not require specification of a stabilization parameter or calculation of higher-order derivatives, and always leads to a symmetric linear system. The Euler semi-implicit scheme is used for the time discretization. It is shown that the proposed fully discrete stabilized finite element method results in the optimal order bounds for the velocity and pressure.     

An explicitly uncoupled VMS stabilization finite element method for the time-dependent Darcy-Brinkman equations in double-diffusive convection

In this article, a full explicitly uncoupled variational multiscale (VMS) stabilization finite element method for solving the Darcy-Brinkman equations in double-diffusive convection is proposed. This method introduces three uncoupled VMS treatments for the velocity, the temperature, and the concentration as the postprocessing steps at each time step, respectively. We only need first to solve three full decoupled linear problems and then to solve three full decoupled postprocessing problems. This method is easy to implement because the existing codes can be used. The unconditional stability is proved and the a priori error estimates are derived. A series of numerical experiments are also given to confirm the theoretical analysis and to demonstrate the efficiency of the new method.     

Fully discrete finite element approaches for time-dependent Maxwell's equations

A fully discrete finite element method is used to approximate the electric field equation derived from time-dependent Maxwell's equations in three dimensional polyhedral domains. Optimal energy-norm error estimates are achieved for general Lipschitz polyhedral domains. Optimal -norm error estimates are obtained for convex polyhedral domains. Received February 3, 1997 / Revised version received February 27, 1998     

A partitioned finite element scheme based on Gauge-Uzawa method for time-dependent MHD equations

In this paper, we mainly introduce a partitioned scheme based on Gauge-Uzawa finite element method for the 2D time-dependent incompressible magnetohydrodynamics (MHD) equations. It is a fully decoupled projection method which combines the Gauge and Uzawa methods within a variational formulation. Firstly, the temporal discretization is based on backward Euler technique for the linear term and semi-implicit scheme for the nonlinear term. Secondly, the spatial approximation of fluid velocity, hydrodynamic pressure, and magnetic field apply the mixed element method. Finally, the validity, reliability, and accuracy of the proposed algorithms are supported by numerical experiments.     

An adaptive finite element method for singular parabolic equations

Summary. In this paper we consider additive Schwarz-type iteration methods for saddle point problems as smoothers in a multigrid method. Each iteration step of the additive Schwarz method requires the solutions of several small local saddle point problems. This method can be viewed as an additive version of a (multiplicative) Vanka-type iteration, well-known as a smoother for multigrid methods in computational fluid dynamics. It is shown that, under suitable conditions, the iteration can be interpreted as a symmetric inexact Uzawa method. In the case of symmetric saddle point problems the smoothing property, an important part in a multigrid convergence proof, is analyzed for symmetric inexact Uzawa methods including the special case of the additive Schwarz-type iterations. As an example the theory is applied to the Crouzeix-Raviart mixed finite element for the Stokes equations and some numerical experiments are presented. Mathematics Subject Classification (1991):65N22, 65F10, 65N30Supported by the Austrian Science Foundation (FWF) under the grant SFB F013}\and Walter Zulehner     

Mixed finite element approximation of eddy current problems

Finite element approximations of eddy current problems thatare entirely based on the magnetic field H are haunted by theneed to enforce the algebraic constraint curl H=0 in non-conductingregions. As an alternative to techniques employing combinatorialSeifert (cutting) surfaces, in order to introduce a scalar magneticpotential we propose mixed multi-field formulations, which enforcethe constraint in the variational formulation. In light of thefact that the computation of cutting surfaces is expensive,the mixed finite element approximation is a viable option despitethe increased number of unknowns.     

Optimal error estimate of the penalty finite element method for the time-dependent Navier-Stokes equations

A fully discrete penalty finite element method is presented for the two-dimensional time-dependent Navier-Stokes equations. The time discretization of the penalty Navier-Stokes equations is based on the backward Euler scheme; the spatial discretization of the time discretized penalty Navier-Stokes equations is based on a finite element space pair which satisfies some approximate assumption. An optimal error estimate of the numerical velocity and pressure is provided for the fully discrete penalty finite element method when the parameters and are sufficiently small.

     

A splitting positive definite mixed finite element method for two classes of integro-differential equations

In this article, we establish a new mixed finite element procedure to solve the second-order hyperbolic and pseudo-hyperbolic integro-differential equations, in which the mixed element system is symmetric positive definite without requiring the LBB consistency condition. Convergence analysis shows that the method yields the approximate solutions with optimal accuracy in L ²(??) norm for u and in H(?div;??) norm for the flux???. Numerical experiments are given to verify the theoretical results.     

Superconvergence analysis of splitting positive definite nonconforming mixed finite element method for pseudo-hyperbolic equations

In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bilinear element is used for u. Superconvergence results in ｜｜·｜｜div,h norm for p and optimal error estimates in L2 norm for u are derived for both semi-discrete and fully discrete schemes under almost uniform meshes.     

Stabilized finite element method based on the Crank-Nicolson extrapolation scheme for the time-dependent Navier-Stokes equations

This paper provides an error analysis for the Crank-Nicolson extrapolation scheme of time discretization applied to the spatially discrete stabilized finite element approximation of the two-dimensional time-dependent Navier-Stokes problem, where the finite element space pair for the approximation of the velocity and the pressure is constructed by the low-order finite element: the quadrilateral element or the triangle element with mesh size . Error estimates of the numerical solution to the exact solution with are derived.

     

Equivalence between two finite element methods for the eddy current problem

The goal of this Note is to prove that two, in principle different, well-known finite element approximations of the eddy current model are equivalent. The first one concerns a formulation involving the magnetic field in the conductor and the magnetic scalar potential in the dielectric. The second one solves another formulation of the same problem involving the magnetic field in both, the conductor and the dielectric, and a Lagrange multiplier in the dielectric. The latter is also shown to be equivalent to a third formulation involving two Lagrange multipliers, which leads to a well posed linear system.     

Local and parallel finite element algorithms for the time-dependent Oseen equations

Numerical Algorithms - Based on two-grid discretizations, local and parallel finite element algorithms are proposed and analyzed for the time-dependent Oseen equations. Using conforming finite...     

The finite element method for parabolic equations

Summary In this first of two papers, computable a posteriori estimates of the space discretization error in the finite element method of lines solution of parabolic equations are analyzed for time-independent space meshes. The effectiveness of the error estimator is related to conditions on the solution regularity, mesh family type, and asymptotic range for the mesh size. For clarity the results are limited to a model problem in which piecewise linear elements in one space dimension are used. The results extend straight-forwardly to systems of equations and higher order elements in one space dimension, while the higher dimensional case requires additional considerations. The theory presented here provides the basis for the analysis and adaptive construction of time-dependent space meshes, which is the subject of the second paper. Computational results show that the approach is practically very effective and suggest that it can be used for solving more general problems.The work was partially supported by ONR Contract N00014-77-C-0623     

An h-p Petrov-Galerkin finite element method for linear Volterra integro-differential equations

We analyze an h-p version Petrov-Galerkin finite element method for linear Volterra integrodifferential equations. We prove optimal a priori error bounds in the L ²- and H ¹-norm that are explicit in the time steps, the approximation orders and in the regularity of the exact solution. Numerical experiments confirm the theoretical results. Moreover, we observe that the numerical scheme superconverges at the nodal points of the time partition.     

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.     

半线性抛物方程的时空有限元方法

讨论半线性抛物方程的连续Galerkin时空有限元方法,利用有限元方法和有限差分方法相结合的技巧,证明了弱解的存在唯一性,给出了时间最大模,空间L~2模,即L~∞(L~2)模误差估计.并给出数值算例证明了连续时空有限元方法对于半线性抛物方程的有效性.     

Nonconforming finite element method for stochastic Stokes equations

In this paper, we study nonconforming finite element method for stochastic Stokes equation driven by white noise. We apply “green function framework” and standard duality technique to study the error estimate for velocity in L²$L^2$-norm and for pressure in H^-1$H^{-1}$-norm. Finally, numerical experiment proves our theoretical results.     

An algorithm using the finite volume element method and its splitting extrapolation

This paper is to present a new efficient algorithm by using the finite volume element method and its splitting extrapolation. This method combines the local conservation property of the finite volume element method and the advantages of splitting extrapolation, such as a high order of accuracy, a high degree of parallelism, less computational complexity and more flexibility than a Richardson extrapolation. Because the splitting extrapolation formulas only require us to solve a set of smaller discrete subproblems on different coarser grids in parallel instead of on the globally fine grid, a large scale multidimensional problem is turned into a set of smaller discrete subproblems. Additionally, this method is efficient for solving interface problems if we regard the interfaces of the problems as the interfaces of the initial domain decomposition.     

An hp finite element method for singularly perturbed systems of reactiondiffusion equations

We consider the approximation of a coupled system of two singularly perturbed reaction-diffusion equations by the finite element method. The solution to such problems contains boundary layers which overlap and interact, and the numerical approximation must take this into account in order for the resulting scheme to converge uniformly with respect to the singular perturbation parameters. We present results on a high order hp finite element scheme which includes elements of size O (εp) and O (μp) near the boundary, where ε, μ are the singular perturbation parameters and p is the degree of the approximating polynomials. Under the assumption of analytic input data, the method yields exponential rates of convergence as p → ∞, independently of ε and μ. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)