Generalized Finite Element Methods

Generalized Finite Element Methods Spring 2003 Properties of Galerkin approximations Suvranu De

Last class Equivalent representations of the mathematical model Strong formulation (BVP) Minimization statement Weak formulation (VBVP) Approximate solution techniques: Rayleigh Ritz Method Galerkin Method Other techniques Equivalence

This class Does the Galerkin method always work? (existence and uniqueness) How significantly does the approximate solution differ from the exact solution? (a priori error analysis) How fast does the approximate solution converge to the analytical solution? (rate of convergence)

VBVP The Galerkin Method • Takes the weak form (VBVP)as the starting point. • Very general. Works for problems that do not have a minimization principle. • For problems that do have a minimization principle, the Galerkin method and the the Rayleigh-Ritz methods produce exactly the same set of discrete equations.

VBVP The Galerkin Method Approximation Pose the weak form: Find such that in the subspace Xh(APPROXIMATION) Find such that

VBVP The Galerkin Method Using the property of bilinearity and linearityas before

VBVP The Galerkin Method Discrete equations The discrete problem Find such that Choose First row

VBVP The Galerkin Method Discrete equations The discrete problem Find such that Choose Second row and so on....

VBVP The Galerkin Method Discretized set of equations The discrete problem Find such that is equivalent to Find such that

A bilinear form is continuous(bounded) if  b > 0 such that A bilinear form is elliptic if  a > 0 such that Existence and uniqueness Definitions... Note that a elliptic bilinear form is always positive definite.

Existence and uniqueness Theorem We state without proof: Let Xh be a finite dimensional subspace of a Hilbert space X, a: XhXhR be a continuous ellipticbilinear form, and l: XhR be a bounded linear functional. Then there exists a unique function uh Xh that satisfies Henceforth we will assume the bilinear form to be continuous and elliptic

Goals Theory A priori... A priorierror estimates: bound various “measures” of u[ exact solution]-uh [approximate solution]; in terms of a constant C (function of W, problem parameters), h [some measure of discretization], and u • Useful for: • Comparison of different discretizations ( which converges faster?) • To understand the conditions that must be satisfied for rapid convergence • To understand if the method is properly implemented (to check a test problem for which the analytical solution is known)

Goals Theory A posteriori A posteriorierror estimates: bound various “measures” of u[ exact solution]-uh [approximate solution]; in terms of a constant C (function of W, problem parameters), h [some measure of discretization], and uh More useful than a priori error estimators during the solution of practical problems (solution unknown). But much more difficult to obtain the constant C (needs additional computational effort)

Error: Energy norm Theory Definition... Define the energy norm||| v ||| as: e.g., for the Dirichlet problem hence For an elastic bar, this is twice the strain energy stored in the bar Note: ||| ||| is problem dependent

Error: Energy norm Theory ...Definition Exact solution: u(x)  X Approximate solution: uh(x)  Xh Discretization error:

Error: Energy norm Theory Orthogonality... Property 1:The discretization error is orthogonal to the approximation space Xh in the energy norm Since then but Adding the last two equations

Error: Energy norm Theory ...Orthogonality Discretization error, eh Exact solution, u Approximate solution, uh Xh

Error: Energy norm Theory General bound Property 2: Example: in the case of linear elasticity, this property means that the strain energy of the discretized model is always bounded from above by the strain energy of the mathematical model. Strain energy of mathematical model Strain energy of discretized model Increased refinement

Error: Energy norm Theory General bound Property 2: 0 Bilinearity and symmetry Since a( , ) is SPD Therefore

Error: Energy norm Theory Best approximation Property 3:Best approximation property In words: even if you knew the exact solution (u) you could not find a wh in Xh more accurate than uh in the energy norm “inf” means infinum

eh Xh u uh wh Error: Energy norm Theory ...Orthogonality Combining properties 1 and 3: of all possible functions in Xh, the Galerkin process picks the one that makes the error orthogonal to Xh and this is the best possible choice in the energy norm.

Error: Energy norm Theory Best approximation Property 3:Best approximation property For any

Error: H1 norm Theory Reminders.. TheH1norm: measures e as well as ex

Error: H1 norm Theory Reminders.. Ellipticity of a( , ): Continuity of a( , ):

Error: H1 norm Theory Cea’s Lemma For a linear problem whose bilinear form is symmetric, continuous and elliptic the error e = u - uh satisfies where What happens if a0 ?

Error: H1 norm Theory Cea’s Lemma From ellipticity: Property 3 Continuity

Error: H1 norm Theory Link to interpolation error Cea’s Lemma From definition where is an interpolant of u.

Error: H1 norm Theory Link to interpolation error Example : if Then we may choose any as a function that interpolates uat 4 points. This function has nothing to do with the approximate solution uh Leads to interpolation error estimate The problem of trying to determine how good our approximate solution is boils down to finding out the interpolation error ( )

Error: H1 norm Theory Example Interpolation error estimate (from last slide) e.g., piecewise linear interpolation (finite element) rate of convergence

Convergence rate Theory Example slope = 1 same slope but lower C(more accurate)

Convergence rate Theory Example Problems in linear elasticity, easier to measure error in strain energy (E-Eh) slope = 2 Reason: from continuity and ellipticity

Summary • Of all possible functions in Xh, the Galerkin process picks the one that makes the error orthogonal to Xh and this is the best possible choice in the energy norm. • Cea’s lemma is a expression of the error estimate in H1 norm and allows us to link the error estimate to the interpolation error estimate. • If piece-wise Lagrange interpolation (of degree k) is used • then the numerical solution exhibits an order k convergence in the H1 norm (in the limit that ‘h’ is small) and an order 2k convergence in the energy norm.

Generalized Finite Element Methods