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.

A note on the efficient evaluation of a modified Hilbert transformation

2020 ◽  
Vol 0 (0) ◽  
Author(s):  
Olaf Steinbach ◽  
Marco Zank
Keyword(s):  
Finite Element ◽  
Evolution Equations ◽  
Time Derivative ◽  
Space Time ◽  
Boundary Integral ◽  
Boundary Element Methods ◽  
Hilbert Transformation ◽  
Element Discretization ◽  
First Order ◽  
Parabolic Evolution Equations

AbstractIn this note we consider an efficient data–sparse approximation of a modified Hilbert type transformation as it is used for the space–time finite element discretization of parabolic evolution equations in the anisotropic Sobolev space H1,1/2(Q). The resulting bilinear form of the first order time derivative is symmetric and positive definite, and similar as the integration by parts formula for the Laplace hypersingular boundary integral operator in 2D. Hence we can apply hierarchical matrices for data–sparse representations and for acceleration of the computations. Numerical results show the efficiency in the approximation of the first order time derivative. An efficient realisation of the modified Hilbert transformation is a basic ingredient when considering general space–time finite element methods for parabolic evolution equations, and for the stable coupling of finite and boundary element methods in anisotropic Sobolev trace spaces.

scholarly journals A priori error estimates for the space-time finite element approximation of a quasilinear gradient enhanced damage model

2021 ◽  
Author(s):  
Marita Holtmannspötter
Keyword(s):  
Finite Element ◽  
Error Estimates ◽  
A Priori ◽  
Damage Model ◽  
Space Time ◽  
A Priori Error Estimates ◽  
Element Approximation ◽  
Element Discretization ◽  
Priori Error Estimates ◽  
Almost All

In this paper we investigate a priori error estimates for the space-time Galerkin finite element discretization of a quasilinear gradient enhanced damage model. The model equations are of a special structure as the state equation consists of two quasilinear elliptic PDEs which have to be fulfilled at almost all times coupled with a nonsmooth, semilinear ODE that has to hold true in almost all points in space. The system is discretized by a constant discontinuous Galerkin method in time and usual conforming linear finite elements in space. Numerical experiments are added to illustrate the proven rates of convergence.

Numerical Analysis of a Finite Element Approximation to a Level Set Model for Free-Surface Flows

2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Tomás Chacón Rebollo ◽  
Macarena Gómez Mármol ◽  
Isabel Sánchez Muñoz
Keyword(s):  
Finite Element ◽  
Free Surface ◽  
Level Set Method ◽  
Level Set ◽  
Free Surface Flows ◽  
Surface Flow ◽  
Optimal Error Estimates ◽  
Element Approximation ◽  
Element Discretization ◽  
Level Set Function

Abstract In this paper, we study a finite element discretization of a Level Set Method formulation of free-surface flow. We consider an Euler semi-implicit discretization in time and a Galerkin discretization of the level set function. We regularize the density and viscosity of the flow across the interface, following the Level Set Method. We prove stability in natural norms when the viscosity and density vary from one to the other layer and optimal error estimates for smooth solutions when the layers have the same density. We present some numerical tests for academic flows.

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.

A Multi-Field Space-Time Finite Element Method for Structural Acoustics

1995 ◽  
Author(s):  
Lonny L. Thompson
Keyword(s):  
Finite Element ◽  
Time Domain ◽  
Acoustic Waves ◽  
Space Time ◽  
High Order ◽  
Elastic Structure ◽  
Structural Acoustics ◽  
Field Representation ◽  
Element Discretization ◽  
Acoustic Fluid

Abstract A Computational Structural Acoustics (CSA) capability for solving scattering, radiation, and other problems related to the acoustics of submerged structures has been developed by employing some of the recent algorithmic trends in Computational Fluid Dynamics (CFD), namely time-discontinuous Galerkin Least-Squares finite element methods. Traditional computational methods toward simulation of acoustic radiation and scattering from submerged elastic bodies have been primarily based on frequency domain formulations. These classical time-harmonic approaches (including boundary element, finite element, and finite difference methods) have been successful for problems involving a limited range of frequencies (narrow band response) and scales (wavelengths) that are large compared to the characteristic dimensions of the elastic structure. Attempts at solving large-scale structural acoustic systems with dimensions that are much larger than the operating wavelengths and which are complex, consisting of many different components with different scales and broadband frequencies, has revealed limitations of many of the classical methods. As a result, there has been renewed interest in new innovative approaches, including time-domain approaches. This paper describes recent advances in the development of a new class of high-order accurate and unconditionally stable space-time methods for structural acoustics which employ finite element discretization of the time domain as well as the usual discretization of the spatial domain. The formulation is based on a space-time variational equation for both the acoustic fluid and elastic structure together with their interaction. Topics to be discussed include the development and implementation of higher-order accurate non-reflecting boundary conditions based on the exact impedance relation through the. Dirichlet-to-Neumann (DtN) map, and a multi-field representation for the acoustic fluid based on independent pressure and velocity potential variables. Numerical examples involving radiation and scattering of acoustic waves are presented to illustrate the high-order accuracy achieved by the new methodology for CSA.

Sign in / Sign up

Export Citation Format

Share Document