An adaptive finite element method for a time‐dependent Stokes problem

In this article, we conduct an a posteriori error analysis of the two‐dimensional time‐dependent Stokes problem with homogeneous Dirichlet boundary conditions, which can be extended to mixed boundary conditions. We present a full time–space discretization using the discontinuous Galerkin method with polynomials of any degree in time and the ? ₂ ? ?₁ Taylor–Hood finite elements in space, and propose an a posteriori residual‐type error estimator. The upper bounds involve residuals, which are global in space and local in time, and an L ²‐error term evaluated on the left‐end point of time step. From the error estimate, we compute local error indicators to develop an adaptive space/time mesh refinement strategy. Numerical experiments verify our theoretical results and the proposed adaptive strategy.     

Galerkin method for a Stefan-type problem in one space dimension

Based on a Landau-type transformation, both continuous and discrete in time L²-Galerkin methods are applied to a single-phase Stefan-type problem in one space dimension. Optimal rates of convergence in L^ℵ, L^∞, and H¹-norms are derived and computational results are presented. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 393–416, 1997     

Superconvergence analysis and a posteriori error estimation of a finite element method for an optimal control problem governed by integral equations

In this paper, we discuss the numerical simulation for a class of constrained optimal control problems governed by integral equations. The Galerkin method is used for the approximation of the problem. A priori error estimates and a superconvergence analysis for the approximation scheme are presented. Based on the results of the superconvergence analysis, a recovery type a posteriori error estimator is provided, which can be used for adaptive mesh refinement. The research project is supported by the National Basic Research Program under the Grant 2005CB321701 and the National Natural Science Foundation of China under the Grant 10771211.     

An adaptive weak Galerkin finite element method with hierarchical bases for the elliptic problem

Based on a posteriori error estimator with hierarchical bases, an adaptive weak Galerkin finite element method (WGFEM) is proposed for the elliptic problem with mixed boundary conditions. For the posteriori error estimator, we are only required to solve a linear algebraic system with diagonal entries corresponding to the degree of freedoms, which significantly reduces the computational cost. The upper and lower bounds of the error estimator are shown to addresses the reliability and efficiency of the adaptive approach. Numerical simulations are provided to demonstrate the effectiveness and robustness of the proposed method.     

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.     

Galerkin methods for a semilinear parabolic problem with nonlocal boundary conditions

We formulate and analyze a Crank-Nicolson finite element Galerkin method and an algebraically-linear extrapolated Crank-Nicolson method for the numerical solution of a semilinear parabolic problem with nonlocal boundary conditions. For each method, optimal error estimates are derived in the maximum norm.Dedicated to Professor J. Crank on the occasion of his 80th birthdaySupported in part by the National Science Foundation grant CCR-9403461.Supported in part by project DGICYT PB95-0711.     

Legendre multiscaling functions for solving the one‐dimensional parabolic inverse problem

An inverse problem concerning diffusion equation with a source control parameter is investigated. The approximation of the problem is based on the Legendre multiscaling basis. The properties of Legendre multiscaling functions are first presented. These properties together with Galerkin method are then utilized to reduce the inverse problem to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the new technique. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Discontinuous Galerkin methods for solving double obstacle problem

In this article the ideas in Wang et al. [SIAM J Numec Anal 48 (2010), 708–73] are extended to solve the double obstacle problem using discontinuous Galerkin methods. A priori error estimates are established for these methods, which reach optimal order for linear elements. We present a test example, and the numerical results on the convergence order match the theoretical prediction. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Mixed FEM and BEM coupling for the three‐dimensional magnetostatic problem

We present two new mixed finite element methods coupled with a boundary method for the three dimensional magnetostatic problem. Such formulations are obtained by coupling a finite element method inside a bounded domain with a boundary integral method involving either the Calderon equations or the inverse of Dirichlet Neumann operator to treat the exterior domain. First, we present the formulations and then prove that our mixed formulations are well posed and that they lead to a convergent Galerkin method. Finally, we give numerical results for a sphere immersed in a homogeneous (source) field in the two formulations. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 443–462, 2003     

On an exact control problem for a semilinear wave equation with an unknown source

The exact controllability of a semilinear wave equation, with Dirichlet boundary control on a part of the boundary and an unknown source, is shown. The nonlinear term has at most a linear growth, the initial and target spaces are L²(Ω)×H⁻¹(Ω).     

Approximation of an optimal control problem governed by a differential parabolic inclusion

In this article, we investigate the optimal control problem governed by parabolic inclusion. We describe the Galerkin approximation and we demonstrate the existence of the strong condensation points of the set of solutions of approximate optimization problems. Each of these points is a solution of the initial optimization problem.     

An inverse problem for operators of a triangular structure

An inverse problem for operators of a triangular structure is studied. An algorithm for the solution and necessary and sufficient conditions for the solvability of this problem are obtained, moreover uniqueness is proved. Applications to difference and differential operators are considered.     

On the coupling of LDG-FEM and BEM methods for the three dimensional magnetostatic problem

We present two new coupling models for the three dimensional magnetostatic problem. In the first model, we propose a new coupled formulation, prove that it is well posed and solves Maxwell’s equations in the whole space. In the second, we propose a new coupled formulation for the Local Discontinuous Galerkin method, the finite element method and the boundary element method. This formulation is obtained by coupling the LDG method inside a bounded domain Ω₁ with the FEM method inside a layer where Ω is a bounded domain which is made up of material of permeability μ and such that , and with a boundary element method involving Calderon’s equations. We prove that our formulation is consistent and well posed and we present some a priori error estimates for the method.     

Pointwise error estimates of the local discontinuous Galerkin method for a second order elliptic problem

In this paper we derive some pointwise error estimates for the local discontinuous Galerkin (LDG) method for solving second-order elliptic problems in (). Our results show that the pointwise errors of both the vector and scalar approximations of the LDG method are of the same order as those obtained in the norm except for a logarithmic factor when the piecewise linear functions are used in the finite element spaces. Moreover, due to the weighted norms in the bounds, these pointwise error estimates indicate that when at least piecewise quadratic polynomials are used in the finite element spaces, the errors at any point depend very weakly on the true solution and its derivatives in the regions far away from . These localized error estimates are similar to those obtained for the standard conforming finite element method.

     

An $h$-Adaptive Runge-Kutta Discontinuous Galerkin Method for Hamilton-Jacobi Equations

In [35, 36], we presented an $h$-adaptive Runge-Kutta discontinuous Galerkin method using troubled-cell indicators for solving hyperbolic conservation laws. A tree data structure (binary tree in one dimension and quadtree in two dimensions) is used to aid storage and neighbor finding. Mesh adaptation is achieved by refining the troubled cells and coarsening the untroubled "children". Extensive numerical tests indicate that the proposed $h$-adaptive method is capable of saving the computational cost and enhancing the resolution near the discontinuities. In this paper, we apply this $h$-adaptive method to solve Hamilton-Jacobi equations, with an objective of enhancing the resolution near the discontinuities of the solution derivatives. One- and two-dimensional numerical examples are shown to illustrate the capability of the method.     

定常不可压Navier—Stokes方程迎风格式

本文利用齐次定解条件对定常不可压Navier—Stokes方程的非线性项进行处理,给出了相应的一种迎风Galerkin有限元算法;针对这种迎风Galerkin有限元算法,在迎风参数满足一定条件下,利用其三项式具有的一些很好性质,更简单地证明了该问题解的存在唯一性。     

A locking-free scheme for the LQR control of a Timoshenko beam

In this paper we analyze a locking-free numerical scheme for the LQR control of a Timoshenko beam. We consider a non-conforming finite element discretization of the system dynamics and a control law constant in the spatial dimension. To solve the LQR problem we seek a feedback control which depends on the solution of an algebraic Riccati equation. An optimal error estimate for the feedback operator is proved in the framework of the approximation theory for control of infinite dimensional systems. This estimate is valid with constants that do not depend on the thickness of the beam, which leads to the conclusion that the method is locking-free. In order to assess the performance of the method, numerical tests are reported and discussed.     

An inverse problem for scattering by a doubly periodic structure

Consider scattering of electromagnetic waves by a doubly periodic structure with for integers , . Above the structure, the medium is assumed to be homogeneous with a constant dielectric coefficient. The medium is a perfect conductor below the structure. An inverse problem arises and may be described as follows. For a given incident plane wave, the tangential electric field is measured away from the structure, say at for some large . To what extent can one determine the location of the periodic structure that separates the dielectric medium from the conductor? In this paper, results on uniqueness and stability are established for the inverse problem. A crucial step in our proof is to obtain a lower bound for the first eigenvalue of the following problem in a convex domain :

     

The element‐free Galerkin method for a quasistatic contact problem with the Tresca friction in elastic materials

This paper is proposed for the error estimates of the element‐free Galerkin method for a quasistatic contact problem with the Tresca friction. The penalty method is used to impose the clamped boundary conditions. The duality algorithm is also given to deal with the non‐differentiable term in the quasistatic contact problem with the Tresca friction. The error estimates indicate that the convergence order is dependent on the nodal spacing, the time step, the largest degree of basis functions in the moving least‐squares approximation, and the penalty factor. Numerical examples demonstrate the effectiveness of the element‐free Galerkin method and verify the theoretical analysis.