scholarly journals Numerical Solutions to Nonsmooth Dirichlet Problems Based on Lumped Mass Finite Element Discretization

2014 ◽  
Vol 2014 ◽  
pp. 1-9
Author(s):  
Haixiong Yu ◽  
Jinping Zeng
Keyword(s):  
Finite Element ◽  
Numerical Solutions ◽  
Piecewise Linear ◽  
Approximation Error ◽  
Uniform Mesh ◽  
Element Approximation ◽  
Dirichlet Problems ◽  
Lumped Mass ◽  
Element Discretization ◽  
Maximum Angle

We apply a lumped mass finite element to approximate Dirichlet problems for nonsmooth elliptic equations. It is proved that the lumped mass FEM approximation error in energy norm is the same as that of standard piecewise linear finite element approximation. Under the quasi-uniform mesh condition and the maximum angle condition, we show that the operator in the finite element problem is diagonally isotone and off-diagonally antitone. Therefore, some monotone convergent algorithms can be used. As an example, we prove that the nonsmooth Newton-like algorithm is convergent monotonically if Gauss-Seidel iteration is used to solve the Newton's equations iteratively. Some numerical experiments are presented.

An Angular Finite Element Method for the Solution of the Even-Parity Equation

2009 ◽  
Author(s):  
R. Becker ◽  
R. Koch ◽  
M. F. Modest ◽  
H.-J. Bauer
Keyword(s):  
Finite Element ◽  
Finite Element Approximation ◽  
Basis Functions ◽  
New Method ◽  
Discrete Ordinates ◽  
Test Case ◽  
Element Approximation ◽  
Promising Alternative ◽  
Element Discretization ◽  
Spatial Coordinates

The present article introduces a new method to solve the radiative transfer equation (RTE). First, a finite element discretization of the solid angle dependence is derived, wherein the coefficients of the finite element approximation are functions of the spatial coordinates. The angular basis functions are defined according to finite element principles on subdivisions of the octahedron. In a second step, these spatially dependent coefficients are discretized by spatial finite elements. This approach is very attractive, since it provides a concise derivation for approximations of the angular dependence with an arbitrary number of angular nodes. In addition, the usage of high-order angular basis functions is straightforward. In the current paper the governing equations are first derived independently of the actual angular approximation. Then, the design principles for the angular mesh are discussed and the parameterization of the piecewise angular basis functions is derived. In the following, the method is applied to two-dimensional test cases which are commonly used for the validation of approximation methods of the RTE. The results reveal that the proposed method is a promising alternative to the well-established practices like the Discrete Ordinates Method (DOM) and provides highly accurate approximations. A test case known to exhibit the ray effect in the DOM verifies the ability of the new method to avoid ray effects.

Composite Finite Element Approximation for Parabolic Problems in Nonconvex Polygonal Domains

2020 ◽  
Vol 20 (2) ◽  
pp. 361-378
Author(s):  
Tamal Pramanick ◽  
Rajen Kumar Sinha
Keyword(s):  
Finite Element ◽  
Finite Elements ◽  
Error Estimates ◽  
Finite Element Approximation ◽  
Parabolic Problems ◽  
Element Approximation ◽  
Dirichlet Problems ◽  
Element Mesh ◽  
Polygonal Domains ◽  
Numer Math

AbstractThe purpose of this paper is to generalize known a priori error estimates of the composite finite element (CFE) approximations of elliptic problems in nonconvex polygonal domains to the time dependent parabolic problems. This is a new class of finite elements which was introduced by [W. Hackbusch and S. A. Sauter, Composite finite elements for the approximation of PDEs on domains with complicated micro-structures, Numer. Math. 75 1997, 4, 447–472] and subsequently modified by [M. Rech, S. A. Sauter and A. Smolianski, Two-scale composite finite element method for Dirichlet problems on complicated domains, Numer. Math. 102 2006, 4, 681–708] for the approximations of stationery problems on complicated domains. The basic idea of the CFE procedure is to work with fewer degrees of freedom by allowing finite element mesh to resolve the domain boundaries and to preserve the asymptotic order convergence on coarse-scale mesh. We analyze both semidiscrete and fully discrete CFE methods for parabolic problems in two-dimensional nonconvex polygonal domains and derive error estimates of order {\mathcal{O}(H^{s}\widehat{\mathrm{Log}}{}^{\frac{s}{2}}(\frac{H}{h}))} and {\mathcal{O}(H^{2s}\widehat{\mathrm{Log}}{}^{s}(\frac{H}{h}))} in the {L^{\infty}(H^{1})}-norm and {L^{\infty}(L^{2})}-norm, respectively. Moreover, for homogeneous equations, error estimates are derived for nonsmooth initial data. Numerical results are presented to support the theoretical rates of convergence.

scholarly journals Fully discrete finite element data assimilation method for the heat equation

2018 ◽  
Vol 52 (5) ◽  
pp. 2065-2082 ◽  
Author(s):  
Erik Burman ◽  
Jonathan Ish-Horowicz ◽  
Lauri Oksanen
Keyword(s):  
Finite Element ◽  
Initial Data ◽  
Heat Equation ◽  
Time Derivative ◽  
Classical Problem ◽  
Classical Case ◽  
Space Time ◽  
Optimal Error Estimates ◽  
Element Approximation ◽  
Element Discretization

We consider a finite element discretization for the reconstruction of the final state of the heat equation, when the initial data is unknown, but additional data is given in a sub domain in the space time. For the discretization in space we consider standard continuous affine finite element approximation, and the time derivative is discretized using a backward differentiation. We regularize the discrete system by adding a penalty on the H2-semi-norm of the initial data, scaled with the mesh-parameter. The analysis of the method uses techniques developed in E. Burman and L. Oksanen [Numer. Math. 139 (2018) 505–528], combining discrete stability of the numerical method with sharp Carleman estimates for the physical problem, to derive optimal error estimates for the approximate solution. For the natural space time energy norm, away from t = 0, the convergence is the same as for the classical problem with known initial data, but contrary to the classical case, we do not obtain faster convergence for the L2-norm at the final time.

scholarly journals Adaptive Semidiscrete Finite Element Methods for Semilinear Parabolic Integrodifferential Optimal Control Problem with Control Constraint

2013 ◽  
Vol 2013 ◽  
pp. 1-11
Author(s):  
Zuliang Lu
Keyword(s):  
Optimal Control ◽  
Finite Element ◽  
Optimal Control Problem ◽  
Control Problem ◽  
A Posteriori Error Estimates ◽  
Finite Element Approximation ◽  
Element Approximation ◽  
Adaptive Finite Element ◽  
Element Discretization ◽  
Semilinear Parabolic

The aim of this work is to study the semidiscrete finite element discretization for a class of semilinear parabolic integrodifferential optimal control problems. We derive a posteriori error estimates inL2(J;L2(Ω))-norm andL2(J;H1(Ω))-norm for both the control and coupled state approximations. Such estimates can be used to construct reliable adaptive finite element approximation for semilinear parabolic integrodifferential optimal control problem. Furthermore, we introduce an adaptive algorithm to guide the mesh refinement. Finally, a numerical example is given to demonstrate the theoretical results.

scholarly journals AN IMPROVED FINITE ELEMENT APPROXIMATION AND SUPERCONVERGENCE FOR TEMPERATURE CONTROL PROBLEMS

2013 ◽  
Vol 18 (5) ◽  
pp. 631-640 ◽  
Author(s):  
Yuelong Tang
Keyword(s):  
Finite Element ◽  
Temperature Control ◽  
Piecewise Linear ◽  
A Priori ◽  
Finite Element Approximation ◽  
Linear Functions ◽  
Control Problems ◽  
Element Approximation ◽  
Adjoint State ◽  
Admissible Controls

In this paper, we consider an improved finite element approximation for temperature control problems, where the state and the adjoint state are discretized by piecewise linear functions while the control is not discretized directly. The numerical solution of the control is obtained by a projection of the adjoint state to the set of admissible controls. We derive a priori error estimates and superconvergence of second-order. Moreover, we present some numerical examples to illustrate our theoretical results.

scholarly journals A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations

2012 ◽  
Vol 2012 ◽  
pp. 1-11 ◽  
Author(s):  
Chuanjun Chen ◽  
Wei Liu
Keyword(s):  
Finite Element ◽  
Piecewise Linear ◽  
A Priori ◽  
Grid Method ◽  
Finite Element Approximation ◽  
Nonlinear Parabolic Equations ◽  
Coarse Grid ◽  
Element Approximation ◽  
Fine Grid ◽  
Nonlinear Parabolic

A two-grid method is presented and discussed for a finite element approximation to a nonlinear parabolic equation in two space dimensions. Piecewise linear trial functions are used. In this two-grid scheme, the full nonlinear problem is solved only on a coarse grid with grid sizeH. The nonlinearities are expanded about the coarse grid solution on a fine gird of sizeh, and the resulting linear system is solved on the fine grid. A priori error estimates are derived with theH1-normO(h+H2)which shows that the two-grid method achieves asymptotically optimal approximation as long as the mesh sizes satisfyh=O(H2). An example is also given to illustrate the theoretical results.

scholarly journals A survey on numerical methods for spectral Space-Fractional diffusion problems

2020 ◽  
Vol 23 (6) ◽  
pp. 1605-1646
Author(s):  
Stanislav Harizanov ◽  
Raytcho Lazarov ◽  
Svetozar Margenov
Keyword(s):  
Finite Element ◽  
Numerical Solution ◽  
Rational Approximation ◽  
Positive Definite ◽  
Second Order ◽  
Local Operator ◽  
Uniform Mesh ◽  
Element Approximation ◽  
Definite Operator ◽  
Positive Definite Operator

AbstractThe survey is devoted to numerical solution of the equation $ {\mathcal A}^\alpha u=f $, 0 < α<1, where $ {\mathcal A} $ is a symmetric positive definite operator corresponding to a second order elliptic boundary value problem in a bounded domain Ω in ℝd. The fractional power $ {\mathcal A}^\alpha $ is a non-local operator and is defined though the spectrum of $ {\mathcal A} $. Due to growing interest and demand in applications of sub-diffusion models to physics and engineering, in the last decade, several numerical approaches have been proposed, studied, and tested. We consider discretizations of the elliptic operator $ {\mathcal A} $ by using an N-dimensional finite element space Vh or finite differences over a uniform mesh with N points. In the case of finite element approximation we get a symmetric and positive definite operator $ {\mathcal A}_h: V_h \to V_h $, which results in an operator equation $ {\mathcal A}_h^{\alpha} u_h = f_h $ for uh ∈ Vh.The numerical solution of this equation is based on the following three equivalent representations of the solution: (1) Dunford-Taylor integral formula (or its equivalent Balakrishnan formula, (2.5), (2) extension of the a second order elliptic problem in Ω  ×  (0, ∞)⊂ ℝd+1 [17,55] (with a local operator) or as a pseudo-parabolic equation in the cylinder (x, t) ∈ Ω  ×  (0, 1), [70, 29], (3) spectral representation (2.6) and the best uniform rational approximation (BURA) of zα on [0, 1], [37,40].Though substantially different in origin and their analysis, these methods can be interpreted as some rational approximation of $ {\mathcal A}_h^{-\alpha} $. In this paper we present the main ideas of these methods and the corresponding algorithms, discuss their accuracy, computational complexity and compare their efficiency and robustness.

Sign in / Sign up

Export Citation Format

Share Document