Construction of locally conservative fluxes for the SUPG method

We consider the construction of locally conservative fluxes by means of a simple postprocessing technique obtained from the finite element solutions of advection diffusion equations. It is known that a naive calculation of fluxes from these solutions yields nonconservative fluxes. We consider two finite element methods: the usual continuous Galerkin finite element method for solving nondominating advection diffusion equations and the streamline upwind/Petrov‐Galerkin method for solving advection dominated problems. We then describe the postprocessing technique for constructing conservative fluxes from the numerical solutions of the general variational formulation. The postprocessing technique requires solving an auxiliary Neumann boundary value problem on each element independently and it produces a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. We provide a convergence analysis for the postprocessing technique. Performance of the technique and the convergence behavior are demonstrated through numerical examples including a set of test problems for advection diffusion equations, advection dominated equations, and drift‐diffusion equations. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1971–1994, 2015     

Finite dimensional approximation of nonlinear problems

Summary We continue here the study of a general method of approximation of nonlinear equations in a Banach space yet considered in [2]. In this paper, we give fairly general approximation results for the solutions in a neighborhood of a simple limit point. We the apply the previous analysis to the study of Galerkin approximations for a class of variationally posed nonlinear problems and to a mixed finite element method for the NavierStokes equations.This work has been completed during a visit at the Université Pierre et Marie Curic and at the Ecole PolytechniqueSupported by the Fonds National Suisse de la Recherche Scientifique     

A study of isogeometric analysis for scalar convection–diffusion equations

Isogeometric analysis (IGA), in combination with the streamline upwind Petrov–Galerkin (SUPG) stabilization, is studied for the discretization of steady-state convection–diffusion equations. Numerical results obtained for the Hemker problem are compared with results computed with the SUPG finite element method of the same order. Using an appropriate parameterization for IGA, the computed solutions are much more accurate than those obtained with the finite element method, both in terms of the size of spurious oscillations and of the sharpness of layers.     

A Petrov–Galerkin spectral method for fourth‐order problems

In this paper, we consider the Petrov–Galerkin spectral method for fourth‐order elliptic problems on rectangular domains subject to non‐homogeneous Dirichlet boundary conditions. We derive some sharp results on the orthogonal approximations in one and two dimensions, which play important roles in numerical solutions of higher‐order problems. By applying these results to a fourth‐order problem, we establish the H²‐error and L²‐error bounds of the Petrov–Galerkin spectral method. Numerical experiments are provided to illustrate the high accuracy of the proposed method and coincide well with the theoretical analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

The Crank–Nicolson–Galerkin finite element method for a nonlocal parabolic equation with moving boundaries

The aim of this article is to establish the convergence and error bounds for the fully discrete solutions of a class of nonlinear equations of reaction–diffusion nonlocal type with moving boundaries, using a linearized Crank–Nicolson–Galerkin finite element method with polynomial approximations of any degree. A coordinate transformation which fixes the boundaries is used. Some numerical tests to compare our Matlab code with some existing moving finite element methods are investigated. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1515–1533, 2015     

A posteriori error control for finite element approximations of elliptic eigenvalue problems

We develop a new approach to a posteriori error estimation for Galerkin finite element approximations of symmetric and nonsymmetric elliptic eigenvalue problems. The idea is to embed the eigenvalue approximation into the general framework of Galerkin methods for nonlinear variational equations. In this context residual-based a posteriori error representations are available with explicitly given remainder terms. The careful evaluation of these error representations for the concrete situation of an eigenvalue problem results in a posteriori error estimates for the approximations of eigenvalues as well as eigenfunctions. These suggest local error indicators that are used in the mesh refinement process.     

Well‐posedness of a generalized time‐harmonic transport equation for acoustics in flow

We study a generalized time‐harmonic transport equation, which appears in the Goldstein equations and allows us to model the acoustic radiation in a flow. We investigate the well‐posedness of this transport problem. The result will be established under the assumption of a Ω‐filling flow, which, in 2D, is simply equivalent to a flow that does not vanish. The approach relies on the method of characteristics, which leads to the resolution of the transport equation along the streamlines, and on general results of functional analysis. The theoretical results are illustrated with numerical results obtained with a Streamline Upwind Petrov‐Galerkin finite element scheme.     

Parameter identification in nonlocal nonlinear evolution equations

We develop an approximation framework for identifying parameters in a general class of nonautonomous, nonlocal and nonlinear evolution equations. After establishing existence and uniqueness of solutions, we present a convergence theory for Galerkin approximations to inverse problems involving these equations. Our approach relies on the theory of maximal monotone operators in Banach spaces. An application to a nonautonomous nonlinear integral equation arising in heat flow is also discussed.     

On a residual local projection method for the Darcy equation

A new symmetric local projection method built on residual bases (RELP) makes linear equal-order finite element pairs stable for the Darcy problem. The derivation is performed inside a Petrov–Galerkin enriching space approach (PGEM) which indicates parameter-free terms to be added to the Galerkin method without compromising consistency. Velocity and pressure spaces are augmented using solutions of residual dependent local Darcy problems obtained after a static condensation procedure. We prove the method achieves error optimality and indicates a way to recover a locally mass conservative velocity field. Numerical experiments validate theory. To cite this article: L.P. Franca et al., C. R. Acad. Sci. Paris, Ser. I 347 (2009).     

A moving Kriging‐based MLPG method for nonlinear Klein–Gordon equation

In this paper, the meshless local Petrov–Galerkin approximation is proposed to solve the 2‐D nonlinear Klein–Gordon equation. We used the moving Kriging interpolation instead of the MLS approximation to construct the meshless local Petrov–Galerkin shape functions. These shape functions possess the Kronecker delta function property. The Heaviside step function is used as a test function over the local sub‐domains. Here, no mesh is needed neither for integration of the local weak form nor for construction of the shape functions. So the present method is a truly meshless method. We employ a time‐stepping method to deal with the time derivative and a predictor–corrector scheme to eliminate the nonlinearity. Several examples are performed and compared with analytical solutions and with the results reported in the extant literature to illustrate the accuracy and efficiency of the presented method. Copyright © 2016 John Wiley & Sons, Ltd.     

热传导对流问题的自适应最小二乘Galerkin/Petrov混合有限元法

对热传导对流问题提出了自适应Galerkin/Petrov最小二乘混合有限元法．该算法对任何速度和压力有限元空间的组合是相容和稳定的(不需要满足Babuka-Brezzi稳定性条件)．利用Verfürth的一般理论,得到了热传导对流问题的残量型的后验误差估计．最后通过几个数值算例验证了方法的有效性．     

Hybridized schemes of the discontinuous Galerkin method for stationary convection–diffusion problems

For stationary linear convection–diffusion problems, we construct and study a new hybridized scheme of the discontinuous Galerkin method on the basis of an extended mixed statement of the problem. Discrete schemes can be used for the solution of equations degenerating in the leading part and are stated via approximations to the solution of the problem, its gradient, the flow, and the restriction of the solution to the boundaries of elements. For the spaces of finite elements, we represent minimal conditions responsible for the solvability, stability, accuracy, and superconvergence of the schemes. A new procedure for the post-processing of solutions of HDG-schemes is suggested.     

Finite difference reaction diffusion equations with nonlinear boundary conditions

This article is concerned with numerical solutions of finite difference systems of reaction diffusion equations with nonlinear internal and boundary reaction functions. The nonlinear reaction functions are of general form and the finite difference systems are for both time-dependent and steady-state problems. For each problem a unified system of nonlinear equations is treated by the method of upper and lower solutions and its associated monotone iterations. This method leads to a monotone iterative scheme for the computation of numerical solutions as well as an existence-comparison theorem for the corresponding finite difference system. Special attention is given to the dynamical property of the time-dependent solution in relation to the steady-state solutions. Application is given to a heat-conduction problem where a nonlinear radiation boundary condition obeying the Boltzmann law of cooling is considered. This application demonstrates a bifurcation property of two steady-state solutions, and determines the dynamic behavior of the time-dependent solution. Numerical results for the heat-conduction problem, including a test problem with known analytical solution, are presented to illustrate the various theoretical conclusions. © 1995 John Wiley & Sons, Inc.     

An initial boundary value problem for a pseudoparabolic equation with a nonlinear boundary condition

An initial boundary value problem for a quasilinear equation of pseudoparabolic type with a nonlinear boundary condition of the Neumann–Dirichlet type is investigated in this work. From a physical point of view, the initial boundary value problem considered here is a mathematical model of quasistationary processes in semiconductors and magnets, which takes into account a wide variety of physical factors. Many approximate methods are suitable for finding eigenvalues and eigenfunctions in problems where the boundary conditions are linear with respect to the desired function and its derivatives. Among these methods, the Galerkin method leads to the simplest calculations. On the basis of a priori estimates, we prove a local existence theorem and uniqueness for a weak generalized solution of the initial boundary value problem for the quasilinear pseudoparabolic equation. A special place in the theory of nonlinear equations is occupied by the study of unbounded solutions, or, as they are called in another way, blow-up regimes. Nonlinear evolutionary problems admitting unbounded solutions are globally unsolvable. In the article, sufficient conditions for the blow-up of a solution in a finite time in a limited area with a nonlinear Neumann–Dirichlet boundary condition are obtained.     

Numerical analysis of a Picard multilevel stabilization of mixed finite volume method for the 2D/3D incompressible flow with large data

In this article, we develop a branch of nonsingular solutions of a Picard multilevel stabilization of mixed finite volume method for the 2D/3D stationary Navier‐Stokes equations without relying on the unique solution condition. The method presented consists of capturing almost all information of initial problem (the nonlinear problems) on the coarsest mesh and then performs one Picard defect correction (the linear problems) on each subsequent mesh based on previous information thus only solving one large linear systems. What is more, the method presented can results in a better coefficient matrix in the model presented with small viscosity. Theoretical results show that the method presented is derived with the convergence rate of the same order as the corresponding finite volume method/finite element method solving the stationary Navier‐Stokes equations on a fine mesh. Therefore, the method presented is definitely more efficient than the standard finite volume method/finite element method. Finally, numerical experiments clearly show the efficiency of the method presented for solving the stationary Navier‐Stokes equations.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 34: 30–50, 2018     

非线性算子方程的泰勒展式算法

本文的目的是给出一种解Ｈｉｌｂｅｒｔ空间中非线性方程的ｋ阶泰勒展式算法（ｋ１）．标准Ｇａｌｅｒｋｉｎ方法可以看作１阶泰勒展式算法，而最优非线性Ｇａｌｅｒｋｉｎ方法可视为２阶泰勒展式算法．我们应用这种算法于定常的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的数值逼近．在一定情景下，最优非线性Ｇａｌｅｒｋｉｎ方法提供比标准Ｇａｌｅｒｋｉｎ方法和非线性Ｇａｌｅｒｋｉｎ方法更高阶的收敛速度．     

Asymptotically compatible schemes for the approximation of fractional Laplacian and related nonlocal diffusion problems on bounded domains

Approximations of solutions of fractional Laplacian equations on bounded domains are considered. Such equations allow global interactions between points separated by arbitrarily large distances. Two approximations are introduced. First, interactions are localized so that only points less than some specified distance, referred to as the interaction radius, are allowed to interact. The resulting truncated problem is a special case of a more general nonlocal diffusion problem. The second approximation is the spatial discretization of the related nonlocal diffusion problem. A recently developed abstract framework for asymptotically compatible schemes is applied to prove convergence results for solutions of the truncated and discretized problem to the solutions of the fractional Laplacian problems. Intermediate results also provide new convergence results for the nonlocal diffusion problem. Special attention is paid to limiting behaviors as the interaction radius increases and the spatial grid size decreases, regardless of how these parameters may or may not be dependent. In particular, we show that conforming Galerkin finite element approximations of the nonlocal diffusion equation are always asymptotically compatible schemes for the corresponding fractional Laplacian model as the interaction radius increases and the grid size decreases. The results are developed with minimal regularity assumptions on the solution and are applicable to general domains and general geometric meshes with no restriction on the space dimension and with data that are only required to be square integrable. Furthermore, our results also solve an open conjecture given in the literature about the convergence of numerical solutions on a fixed mesh as the interaction radius increases.     

Superconvergence of nonconforming finite element method for the Stokes equations

This article derives a general superconvergence result for nonconforming finite element approximations of the Stokes problem by using a least‐squares surface fitting method proposed and analyzed recently by Wang for the standard Galerkin method. The superconvergence result is based on some regularity assumption for the Stokes problem and is applicable to any nonconforming stable finite elements with regular but nonuniform partitions. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 143–154, 2002; DOI 10.1002/num.1036     

Post-Buckling Range of Plates in Axial Compression with Uncertain Initial Geometric Imperfections

The method of reliable solutions alias the worst scenario method is applied to the problem of von Karman equations with uncertain initial deflection. Assuming two-mode initial and total deflections and using Galerkin approximations, the analysis leads to a system of two nonlinear algebraic equations with one or two uncertain parameters-amplitudes of initial deflections. Numerical examples involve (i) minimization of lower buckling loads and (ii) maximization of the maximal mean reduced stress.     

Enriched finite element spaces for transient conduction heat transfer

This paper presents an alternative approach via finite elements to treat numerically the thermal shocks in heat transfer finite element analysis. The method consists in using the standard enriched finite element approaches with time-interpolation. It will be applied here to the transient conduction heat equation where the classical Galerkin method is shown to be unstable. The proposed method consists in adding and eliminating bubbles to the finite element space and then to interpolate the solution to the real time step. This modification is equivalent to the addition of a stabilizing term tuned by a local time-dependent stability parameter, which ensures an oscillating-free solution. To validate this approach, the numerical results obtained in classical 2D and 3D benchmark problems are compared with the Galerkin and the analytical solutions.