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     

Solvability of the clamped Euler–Bernoulli beam equation

In this study, solvability of the initial boundary value problem for general form Euler–Bernoulli beam equation which includes also moving point-loads is investigated. The complete proof of an existence and uniqueness properties of the weak solution of the considered equation with Dirichlet type boundary conditions is derived. The method used here is based on Galerkin approximation which is the main tool for the weak solution theory of linear evolution equations as well as in derivation of a priori estimate for the approximate solutions. All steps of the proposed technique are explained in detail.     

Boundary element methods for potential problems with nonlinear boundary conditions

Galerkin boundary element methods for the solution of novel first kind Steklov-Poincaré and hypersingular operator boundary integral equations with nonlinear perturbations are investigated to solve potential type problems in two- and three-dimensional Lipschitz domains with nonlinear boundary conditions. For the numerical solution of the resulting Newton iterate linear boundary integral equations, we propose practical variants of the Galerkin scheme and give corresponding error estimates. We also discuss the actual implementation process with suitable preconditioners and propose an optimal hybrid solution strategy.

     

On spline collocation methods for boundary integral equations in the plane

Nowadays boundary elemen; methods belong to the most popular numerical methods for solving elliptic boundary value problems. They consist in the reduction of the problem to equivalent integral equations (or certain generalizations) on the boundary Γ of the given domain and the approximate solution of these boundary equations. For the numerical treatment the boundary surface is decomposed into a finite number of segments and the unknown functions are approximated by corresponding finite elements and usually determined by collocation and Galerkin procedures. One finds the least difficulties in the theoretical foundation of the convergence of Galerkin methods for certain classes of equations, whereas the convergence of collocation methods, which are mostly used in numerical computations, has yet been proved only for special equations and methods. In the present paper we analyse spline collocation methods on uniform meshes with variable collocation points for one-dimensional pseudodifferential equations on a closed curve with convolutional principal parts, which encompass many classes of boundary integral equations in the plane. We give necessary and sufficient conditions for convergence and prove asymptotic error estimates. In particular we generalize some results on nodal and midpoint collocation obtained in [2], [7] and [8]. The paper is organized as follows. In Section 1 we formulate the problems and the results, Section 2 deals with spline interpolation in periodic Sobolev spaces, and in Section 3 we prove the convergence theorems for the considered collocation methods.     

A boundary integral equation method for three-dimensional crack problems in elasticity

This paper presents a solution procedure for three-dimensional crack problems via first kind boundary integral equations on the crack surface. The Dirichlet (Neumann) problem is reduced to a system of integral equations for the jump of the traction (of the field) across the crack surface. The calculus of pseudodifferential operators is used to derive existence and regularity of the solutions of the integral equations. With the concept of the principal symbol and the Wiener-Hopf technique we derive the explicit behavior of the densities of the integral equations near the edge of the crack surface. Based on the detailed regularity results we show how to improve the boundary element Galerkin method for our integral equations. Quasi-optimal asymptotic estimates for the Galerkin error are given.     

二维N－S方程的Fourier非线性Galerkin方法

本文对周期边界条件Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，证明了其Ｆｏｕｒｉｅｒ非线性Ｇａｌｅｒｋｉｎ逼近解的存在唯一性，同时给出了逼近解的误差估计．     

A stabilized semi-implicit Galerkin scheme for Navier-Stokes equations

By introducing a time relaxation term for the time derivative of higher frequency components, we proposed a stabilized semi-implicit Galerkin scheme for evolutionary Navier-Stokes equations in this paper. Analysis shows that such a scheme has weaker stability conditions than that of a classical semi-implicit Galerkin scheme and, when a suitable relaxation parameter σ is chosen, it generates an approximate solution with the same accuracy as the classical one. That means the proposed scheme might use a larger time step to generate a bounded approximate solution. Thus it is more suitable for long time simulations.     

Application of spectral filtering to discontinuous Galerkin methods on triangulations

We adapt the spectral viscosity (SV) formulation implemented as a modal filter to a discontinuous Galerkin (DG) method solving hyperbolic conservation laws on triangular grids. The connection between SV and spectral filtering, which is undertaken for the first time in the context of DG methods on unstructured grids, allows to specify conditions on the filter strength regarding time step choice and mesh refinement. A crucial advantage of this novel damping strategy is its low computational cost. We furthermore obtain new error bounds for filtered Dubiner expansions of smooth functions. While high order accuracy with respect to the polynomial degree N is proven for the filtering procedure in this case, an adaptive application is proposed to retain the high spatial approximation order. Although spectral filtering stabilizes the scheme, it leaves weaker oscillations. Therefore, as a postprocessing step, we apply the image processing technique of digital total variation (DTV) filtering in the new context of DG solutions and prove conservativity in the limit for this filtering procedure. Numerical experiments for scalar conservation laws confirm the designed order of accuracy of the DG scheme with adaptive modal filtering for polynomial degrees up to 8 and the viability of spectral and DTV filtering in case of shocks. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011     

Efficient methods using high accuracy approximate inertial manifolds

Summary. We extend the idea of the post-processing Galerkin method, in the context of dissipative evolution equations, to the nonlinear Galerkin, the filtered Galerkin, and the filtered nonlinear Galerkin methods. In general, the post-processing algorithm takes advantage of the fact that the error committed in the lower modes of the nonlinear Galerkin method (and Galerkin method), for approximating smooth, bounded solutions, is much smaller than the total error of the method. In each case, an improvement in accuracy is obtained by post-processing these more accurate lower modes with an appropriately chosen, highly accurate, approximate inertial manifold (AIM). We present numerical experiments that support the theoretical improvements in accuracy. Both the theory and computations are presented in the framework of a two dimensional reaction-diffusion system with polynomial nonlinearity. However, the algorithm is very general and can be implemented for other dissipative evolution systems. The computations clearly show the post-processed filtered Galerkin method to be the most efficient method. Received September 10, 1998 / Revised version received April 26, 1999 / Published online July 12, 2000     

An Element-Free Galerkin-Based Multi-objective Optimization of Laminated Composite Plates

An element-free Galerkin method is presented to analyze isotropic and laminated composite plates. This method employs the moving least square technique to approximate functions. In the analysis procedure, a collocation method is used to enforce boundary conditions. A consistent multi-objective optimization procedure is also applied. The function introduced here consists of minimizing the weight and cost, as well as of maximizing the load. A genetic algorithm is used for the optimization process.     

A new symmetric interior penalty discontinuous Galerkin formulation for the Serre–Green–Naghdi equations

In this work, we investigate the construction of a new discontinuous Galerkin discrete formulation to approximate the solution of Serre–Green–Naghdi (SGN) equations in the one-dimensional horizontal framework. Such equations describe the time evolution of shallow water free surface flows in the fully nonlinear and weakly dispersive asymptotic approximation regime. A new non-conforming discrete formulation belonging to the family of symmetric interior penalty discontinuous Galerkin methods is introduced to accurately approximate the solutions of the second order elliptic operator occurring in the SGN equations. We show that the corresponding discrete bilinear form enjoys some consistency and coercivity properties, thus ensuring that the corresponding discrete problem is well-posed. The resulting global discrete formulation is then validated through an extended set of benchmarks, including convergence studies and comparisons with data taken from experiments.     

Imposing boundary conditions in Sinc method using highest derivative approximation

Sinc approximate methods are often used to solve complex boundary value problems such as problems on unbounded domains or problems with endpoint singularities. A recent implementation of the Sinc method [Li, C. and Wu, X., Numerical solution of differential equations using Sinc method based on the interpolation of the highest derivatives, Applied Mathematical Modeling 31 (1) 2007 1–9] in which Sinc basis functions are used to approximate the highest derivative in the governing equation of the boundary value problem is evaluated for structural mechanics applications in which interlaminar stresses are desired. We suggest an alternative approach for specifying the boundary conditions, and we compare the numerical results for analysis of a laminated composite Timoshenko beam, implementing both Li and Wu’s approach and our alternative approach for applying the boundary conditions. For the Timoshenko beam problem, we obtain accurate results using both approaches, including transverse shear stress by integration of the 3D equilibrium equations of elasticity. The beam results indicate our approach is less dependent on the selection of the Sinc mesh size than Li and Wu’s SIHD. We also apply SIHD to analyze a classical laminated composite plate. For the plate example, we experience difficulty in obtaining a complete system of equations using Li and Wu’s approach. For our approach, we suggest that additional necessary information may be obtained by applying the derivatives of the boundary conditions on each edge. Using this technique, we obtain accurate results for deflection and stresses, including interlaminar stresses by integration of the 3D equilibrium equations of elasticity. Our results for both the beam and the plate problems indicate that this approach is easily implemented, has a high level of accuracy, and good convergence properties.     

Bernstein Ritz‐Galerkin method for solving an initial‐boundary value problem that combines Neumann and integral condition for the wave equation

In this article, the Ritz‐Galerkin method in Bernstein polynomial basis is implemented to give an approximate solution of a hyperbolic partial differential equation with an integral condition. We will deal here with a type of nonlocal boundary value problem, that is, the solution of a hyperbolic partial differential equation with a nonlocal boundary specification. The nonlocal conditions arise mainly when the data on the boundary cannot be measured directly. The properties of Bernstein polynomial and Ritz‐Galerkin method are first presented, then Ritz‐Galerkin method is used to reduce the given hyperbolic partial differential equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique presented in this article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Galerkin方法的一种简单后处理过程

基于近似惯性流形思想,以流函数形式定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程为例,给出了一种简单的后处理Ｇａｌｅｒｋｉｎ方法。其主要思想是利用近似惯性流形概念和对真解的一种新的分解,构造高低频分量间的近似作用规律。文中证明了这种简单的后处理Ｇａｌｅｒｋｉｎ方法可以较小的代价获得较经典Ｇａｌｅｒｋｉｎ方法高得多的精度。     

The h-p version of the coupling of finite element and boundary element methods for transmission problems in polyhedral domains

Summary. This paper analyzes the rate of convergence of the h-p version of the coupling of the finite element and boundary element method for transmission problems with a linear differential operator with variable coefficients in a bounded polyhedral domain and with constant coefficients in the exterior domain . This procedure uses the variational formulation of the differential equation in and involves integral operators on the interface between and . The finite elements are used to obtain approximate solutions of the differential equation in and the boundary elements are used to obtain approximate solutions of the integral equations. For given piecewise analytic data we show that the Galerkin solution of this coupling procedure converges exponentially fast in the energy norm if the h-p version is used both for finite elements and boundary elements. Received February 10, 1996 / Revised version received April 4, 1997     

Boundary integral approximation of a heat-diffusion problem in time-harmonic regime

In this paper we propose and analyse numerical methods for the approximation of the solution of Helmholtz transmission problems in the half plane. The problems we deal with arise from the study of some models in photothermal science. The solutions to the problem are represented as single layer potentials and an equivalent system of boundary integral equations is derived. We then give abstract necessary and sufficient conditions for convergence of Petrov–Galerkin discretizations of the boundary integral system and show for three different cases that these conditions are satisfied. We extend the results to other situations not related to thermal science and to non-smooth interfaces. Finally, we propose a simple full discretization of a Petrov–Galerkin scheme with periodic spline spaces and show some numerical experiments.     

High-order discontinuous element-based schemes for the inviscid shallow water equations: Spectral multidomain penalty and discontinuous Galerkin methods

Two commonly used types of high-order-accuracy element-based schemes, collocation-based spectral multidomain penalty methods (SMPM) and nodal discontinuous Galerkin methods (DGM), are compared in the framework of the inviscid shallow water equations. Differences and similarities in formulation are identified, with the primary difference being the dissipative term in the Rusanov form of the numerical flux for the DGM that provides additional numerical stability; however, it should be emphasized that to arrive at this equivalence between SMPM and DGM requires making specific choices in the construction of both methods; these choices are addressed. In general, both methods offer a multitude of choices in the penalty terms used to introduce boundary conditions and stabilize the numerical solution. The resulting specialized class of SMPM and DGM are then applied to a suite of six commonly considered geophysical flow test cases, three linear and three non-linear; we also include results for a classical continuous Galerkin (i.e., spectral element) method for comparison. Both the analysis and numerical experiments show that the SMPM and DGM are essentially identical; both methods can be shown to be equivalent for very special choices of quadrature rules and Riemann solvers in the DGM along with special choices in the type of penalty term in the SMPM. Although we only focus our studies on the inviscid shallow water equations the results presented should be applicable to other systems of nonlinear hyperbolic equations (such as the compressible Euler equations) and extendable to the compressible and incompressible Navier-Stokes equations, where viscous terms are included.     

Almost global existence for quasilinear wave equations in waveguides with Neumann boundary conditions

In this paper, we prove almost global existence of solutions to certain quasilinear wave equations with quadratic nonlinearities in infinite homogeneous waveguides with Neumann boundary conditions. We use a Galerkin method to expand the Laplacian of the compact base in terms of its eigenfunctions. For those terms corresponding to zero modes, we obtain decay using analogs of estimates of Klainerman and Sideris. For the nonzero modes, estimates for Klein-Gordon equations, which provide better decay, are available.

     

Galerkin approximations with embedded boundary conditions for retarded delay differential equations

Finite-dimensional approximations are developed for retarded delay differential equations (DDEs). The DDE system is equivalently posed as an initial-boundary value problem consisting of hyperbolic partial differential equations (PDEs). By exploiting the equivalence of partial derivatives in space and time, we develop a new PDE representation for the DDEs that is devoid of boundary conditions. The resulting boundary condition-free PDEs are discretized using the Galerkin method with Legendre polynomials as the basis functions, whereupon we obtain a system of ordinary differential equations (ODEs) that is a finite-dimensional approximation of the original DDE system. We present several numerical examples comparing the solution obtained using the approximate ODEs to the direct numerical simulation of the original non-linear DDEs. Stability charts developed using our method are compared to existing results for linear DDEs. The presented results clearly demonstrate that the equivalent boundary condition-free PDE formulation accurately captures the dynamic behaviour of the original DDE system and facilitates the application of control theory developed for systems governed by ODEs.     

A Comparison Study of Deep Galerkin Method and Deep Ritz Method for Elliptic Problems with Different Boundary Conditions

Recent years have witnessed growing interests in solving partial differential equations by deep neural networks, especially in the high-dimensional case. Unlike classical numerical methods, such as finite difference method and finite element method, the enforcement of boundary conditions in deep neural networks is highly nontrivial. One general strategy is to use the penalty method. In the work, we conduct a comparison study for elliptic problems with four different boundary conditions, i.e., Dirichlet, Neumann, Robin, and periodic boundary conditions, using two representative methods: deep Galerkin method and deep Ritz method. In the former, the PDE residual is minimized in the least-squares sense while the corresponding variational problem is minimized in the latter. Therefore, it is reasonably expected that deep Galerkin method works better for smooth solutions while deep Ritz method works better for low-regularity solutions. However, by a number of examples, we observe that deep Ritz method can outperform deep Galerkin method with a clear dependence of dimensionality even for smooth solutions and deep Galerkin method can also outperform deep Ritz method for low-regularity solutions.Besides, in some cases, when the boundary condition can be implemented in an exact manner, we find that such a strategy not only provides a better approximate solution but also facilitates the training process.