Basic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
|
|
- Julian Dawson
- 6 years ago
- Views:
Transcription
1 Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg Summer School Challenges in Applied Control and Optimal Control Basque Center for Applied Mathematics Derio, July 4-8, 2011
2 The Loop in Adaptive Finite lement Methods (AFM) Adaptive Finite lement Methods (AFM) consist of successive loops of the cycle SOLV: SOLV = STIMAT = MARK = RFIN Numerical solution of the F discretized problem STIMAT: Residual and hierarchical a posteriori error estimators rror estimators based on local averaging Goal oriented weighted dual approach Functional type a posteriori error bounds MARK: RFIN: Strategies based on the max. error or the averaged error Bulk criterion for AFMs Bisection or red/green refinement or combinations thereof
3 A Posteriori rror stimation I For a closed subspace V H 1 (Ω) we assume a(, ) : V V R to be a bounded, V-elliptic bilinear form, i.e., a(v,w) C v k,ω w k,ω, v,w V, a(v,v) γ v 2 k,ω, v V, for some constants C > 0 and γ > 0. We further assume l V where V denotes the algebraic and topological dual of V and consider the variational equation: Find u V such that a(u,v) = l(v), v V. It is well-known by the Lax-Milgram Lemma that under the above assumptions the variational problem admits a unique solution.
4 A Posteriori rror stimation II xample. The standard example is Poisson s equation: Let Γ = Γ D Γ N,Γ D Γ N = and V := {v H 1 (Ω) v ΓD = 0} and a(v,w) := ( v, w) 0,Ω, v,w V, l(v) := (f,v) 0,Ω + (g,v) 0,ΓN, v V, where f L 2 (Ω) and g L 2 (Γ N ). Then, the variational equation represents the weak form of Poisson s problem u = f in Ω, u = 0 on Γ D, n u = g on Γ N, where n stands for the unit outward normal on Γ N.
5 A Posteriori rror stimation III xample cont d. We remark that in case Γ N = we have V = H 1 0(Ω) (pure homogeneous Dirichlet boundary conditions), whereas in case Γ D = the appropriate function space is V = {v H 1 (Ω) (v,χ Ω ) 0,Ω = 0}, where χ Ω stands for the characteristic function of Ω. In these cases, the V-ellipticity of the bilinear form a(, ) follows from the Poincaré-Friedrichs inequalities ( ) v 0,Ω C v 1,Ω + vds, ( v 0,Ω C v 1,Ω + Γ ) vdx. Ω
6 A Posteriori rror stimation IV An important issue in the theory of partial differential equations is the regularity of a solution. For instance, considering Poisson s problem, it is well-known that for a convex domain Ω R 2 with piecewise smooth boundary Γ and f L 2 (Ω) the weak solution u V satisfies u V H 2 (Ω) and u 2,Ω C f 0,Ω. However, there is less regularity, if the domain Ω is no longer convex. A classical example is the so-called L-shaped domain: Consider the Poisson equation in Ω := ( 1, +1) (0,1) ( 1,0) ( 1,0], Γ D := {0} [ 1,0] [0,1] {0}, Γ N := Γ \ Γ D, and assume f 0 and g such that u(r,ϕ) = r 2/3 sin( 2 3 ϕ) is the exact solution of the problem. The solution belongs to V H 5/3 ε (Ω) for any ε > 0 (but not to H 5/3 (Ω)!) and has a singularity in the origin.
7 A Posteriori rror stimation V Finite element approximations are based on the Ritz-Galerkin approach: Given a finite dimensional subspace V h V of test/trial functions, find u h V h such that a(u h,v h ) = l(v h ), v h V h. Since V h V, the existence and uniqueness of a discrete solution u h V h follows readily from the Lax-Milgram Lemma. Moreover, we deduce that the error e u := u u h satisfies the Galerkin orthogonality a(u u h,v h ) = 0, v h V h, i.e., the approximate solution u h V h is the projection of the solution u V onto V h with respect to the inner product a(, ) on V (elliptic projection). Using the Galerkin orthogonality, it is easy to derive the a priori error estimate u u h 1,Ω M inf v h V h u v h 1,Ω, where M := C/γ. This result tells us that the error is of the same order as the best approximation of the solution u V by functions from the finite dimensional subspace V h. It is known as Céa s Lemma.
8 A Posteriori rror stimation VI The Ritz-Galerkin method also gives rise to an a posteriori error estimate in terms of the residual r : V R In fact, it follows that for any v V r(v) := l(v) a(u h,v), v V. γ u u h 2 1,Ω a(u u h,u u h ) = r(u u h ) r 1,Ω u u h 1,Ω, whence u u h 1,Ω 1 γ sup r(v). v V v 1,Ω
9 Department of Mathematics, University of Houston Reliability and fficiency
10 Reliability and fficiency of rror stimators I Definition. An error estimator η h is called reliable, if it provides an upper bound for the error up to data oscillations osc rel h, i.e., if there exists a constant C rel > 0, independent of the mesh size h of the underlying triangulation, such that e u a C rel η h + osc rel h. On the other hand, an estimator η h is said to be efficient, if up to data oscillations osc eff h it gives rise to a lower bound for the error, i.e., if there exists a constant C eff > 0, independent of the mesh size h of the underlying triangulation, such that η h C eff e u a + osc eff h. Finally, an estimator η h is called asymptotically exact, if it is both reliable and efficient with C rel = C 1 eff.
11 Reliability and fficiency of rror stimators II Remark. The notion reliability is motivated by the use of the error estimator in error control. Given a tolerance tol, an idealized termination criterion would be e u a tol. Since the error e u a is unknown, we replace it with the upper bound, i.e., C rel η h + osc rel h tol. We note that the termination criterion both requires the knowledge of C rel and the incorporation of the data oscillation term osc rel and osc rel h 0, it reduces to η h tol. h. In the special case C rel = 1 An alternative, but less used termination criterion is based on the lower bound, i.e., we require 1 ( ) η h osc eff h tol. C eff Typically, this criterion leads to less refinement and thus requires less computational time which motivates to call the estimator efficient.
12 Department of Mathematics, University of Houston Residual-Type A Posteriori rror stimation
13 The error estimate The Role of the Residual u u h 1,Ω 1 γ sup r(v) v V v 1,Ω shows that in order to assess the error e u a we are supposed to evaluate the norm of the residual with respect to the dual space V, i.e., r V := r(v) sup. v V\{0} v a In particular, we have the equality r V = e u a, whereas for the relative error of r(v),v V, as an approximation of e u a we obtain ( e u a r(v)) = 1 e u a 2 v e u 2 e u a, v V with v a = 1. a The goal is to obtain lower and upper bounds for r V at relatively low computational expense.
14 Model problem: Let Ω be a bounded simply-connected polygonal domain in uclidean space lr 2 with boundary Γ = Γ D Γ N, Γ D Γ N = and consider the elliptic boundary value problem Lu := (a u) = f in Ω, u = 0 on Γ D, n a u = g on Γ N, where f L 2 (Ω), g L 2 (Γ N ) and a = (a ij ) 2 i,j=1 is supposed to be a matrix-valued function with entries a ij L (Ω), that is symmetric and uniformly positive definite. The vector n denotes the exterior unit normal vector on Γ N. Setting H 1 0,Γ D (Ω) := { v H 1 (Ω) v ΓD = 0 }, the weak formulation is as follows: Find u H 1 0,Γ D (Ω) such that where a(v,w) := Ω a(u,v) = l(v), v H 1 0,Γ D (Ω), a v w dx, l(v) := f v dx + Ω Γ N g v dσ, v H 1 0,Γ D (Ω).
15 F Approximation: Given a geometrically conforming simplicial triangulation T h of Ω, we denote by S 1,ΓD (Ω; T h ) := { v h H 1 0,Γ D (Ω) v h T P 1 (K), T T h } the trial space of continuous, piecewise linear finite elements with respect to T h. Note that P k (T), k 0, denotes the linear space of polynomials of degree k on T. In the sequel we will refer to N h (D) and h (D), D Ω as the sets of vertices and edges of T h on D. We further denote by T the area, by h T the diameter of an element T T h, and by h = the length of an edge h (Ω Γ N ). We refer to f T := T 1 T fdx the integral mean of f with respect to an element T T h and to g := 1 gds the mean of g with respect to the edge h(γ N ). The conforming P1 approximation reads as follows: Find u h S 1,ΓD (Ω; T h ) such that a(u h,v h ) = l(v h ), v h S 1,ΓD (Ω; T h ).
16 The residual r is given by r(v) := f v dx + Representation of the Residual I Ω Γ N g vds a(u h,v), v V. Applying Green s formula elementwise yields a(u h,v) = a u h v dx = [n a u h ] v ds + n a u h v ds, T T h T h (Ω) h (Γ N ) where [n a u h ] denotes the jump of the normal derivative of u h across h (Ω) and where we have used that u h 0 on T T h, since u h T P 1 (T). We thus obtain r(v) := r T (v) + r (v). T T h h (Ω Γ N )
17 Department of Mathematics, University of Houston Representation of the Residual II Here, the local residuals r T (v),t T h, are given by r T (v) := (f Lu h )v dx, whereas for r (v) we have r (v) := r (v) := T [n a u h ]v ds, h (Ω), (g n a u h ) v ds, h (Γ N ).
18 A Posteriori rror stimator and Data Oscillations The error estimator η h consists of element residuals η T,T T h, and edge residuals η, H (Ω Γ N ), according ( to 1/2, η h := ηt 2 + η) 2 T T h H (Ω Γ N ) where η T and η are given by η T := h T f T Lu h 0,T, T T h, { h 1/2 η := [n a u h] 0,, h (Ω), h 1/2 g n a u h 0,, h (Γ N ). The a posteriori error analysis ( further invokes the data oscillations 1/2, osc h := osc 2 T(f) + osc(g)) 2 T T h h (Γ N ) where osc T (f) and osc (g) are given by osc T (f) := h T f f T 0,T, osc (g) := h 1/2 g g 0,.
19 Clément s Quasi-Interpolation Operator I For p N h (Ω) N h (Γ N ) we denote by ϕ p the basis function in S 1,ΓD (Ω; T h ) with supporting point p, and we refer to D p as the set D p := { T T h p N h (T) }. We refer to π p as the L 2 -projection onto P 1 (D p ), i.e., (π p (v),w) 0,Dp = (v,w) 0,Dp, w P 1 (D p ), where (, ) 0,Dp stands for the L 2 -inner product on L 2 (D p ) L 2 (D p ). Then, Clément s interpolation operator P C is defined as follows P C : L 2 (Ω) S 1,ΓD (Ω, T h ), P C v := π P (v) ϕ P. p N h (Ω) N h (Γ N )
20 Clément s Quasi-Interpolation Operator II Theorem. Let v H 1 0,Γ D (Ω). Then, for Clément s interpolation operator there holds P C v 0,T C v 0,D (1) T v P C v 0,T C h T v 1,D (1) T Further, we have ( ( h (Ω) h (Γ N ), P C v 0, C v (1) 0,D, v P C v 0, C h 1/2 v 2 µ,d (1) K T T h v 2 µ,d (1), P C v 0,T C v (1) 0,D, T. v 1,D (1) ) 1/2 C v µ,ω, 0 µ 1, ) 1/2 C v µ,ω, 0 µ 1. where D (1) T := { T T h N h (T ) N h (T) }, D (1) := { T T h N h () N h (T ) }.
21 Reliability of the A posteriori rror stimator I Theorem. There exist constants Γ R and Γ osc > 0 depending only on the shape regularity of T h such that e u a Γ R η h + Γ osc osc h. Proof. Setting v = e u, we have e u a = a(e u,e u ) = r(e u ) = r(p C e u ) + r(e u P C e u ). For the first term on the right-hand side, by Galerkin orthogonality we obtain r(p C e u ) = f P C e u dx + g P C e u ds a(u h,p C e u ) = 0. Γ N Ω
22 Proof cont d. On the other hand, for the second term on the right-hand side of Green s formula yields r(e u P C e u ) = f (e u P C e u ) dx + g (e u P C e u ) ds + } a {{ u } h (e u P C e u ) dx Ω Γ T T h N T = Lu h n T a u h (e u P C e u ) ds = (f T Lu h ) (e u P C e u ) dx T T h T T h + T h (Ω) + T T h T [n a u h ] (e u P C e u ) ds + (f f T ) (e u P C e u ) dx + h (Γ N ) T h (Γ N ) (g n a u h ) (e u P C e u ) ds (g g ) (e u P C e u ) ds.
23 Proof cont d. In view of the local approximation properties of Clément s quasiinterpolation operator, it follows that ( r(e u P C e u ) C h 2 T f T Lu h 2 0,T) 1/2 e u 1,Ω T T h + h [n a u h ] J 2 0,) 1/2 e u 1,Ω + + h (Ω) h (Γ N ) h (Γ N ) h g n a u h 2 0,) 1/2 e u 1,Ω + ( T T h h 2 T f f T 2 0,T) 1/2 e u 1,Ω h g g 2 0,) 1/2 e u 1,Ω ), from which we may conclude.
24 lement and dge Bubble Functions I The element bubble function ψ T is defined by means of the barycentric coordinates λ T i,1 i 3, according to ψ T := 27 λ T 1 λ T 2 λ T 3. Note that supp ψ T = T int, i.e., ψ T T = 0, T T h. On the other hand, for h (Ω) h (Γ N ) and T T h such that T and p i N h (), 1 i 2, we introduce the edge-bubble functions ψ ψ := 4 λ T 1 λ T 2. Note that ψ = 0 for h (T),.
25 lement and dge Bubble Functions II The bubble functions ψ T and ψ have the following important properties that can be easily verified taking advantage of the affine equivalence of the finite elements: Lemma. There holds p h 2 0,T C p h 2 0, C T p 2 h ψ T dx, p 2 h ψ dσ, p h P 1 (T), p h P 1 (), p h ψ T 1,T C h 1 T p h 0,T, p h P 1 (T), p h ψ T 0,T C p h 0,T, p h P 1 (T), p h ψ 0, C p h 0,, p h P 1 ().
26 lement and dge Bubble Functions III For functions p h P 1 (), h (Ω) h (Γ N ) we further need an extension p h L2 (T) where T T h such that T. For this purpose we fix some T,, and for x T denote by x that point on such that (x x ). For p h P 1 () we then set p h := p h (x ). as the union of elements T T h con- Further, for h (Ω) h (Γ N ) we define D (2) taining as a common edge D (2) := { K T h h (T) }.
27 Department of Mathematics, University of Houston Lemma. There holds lement and dge Bubble Functions IV p h ψ (2) 1,D p h ψ (2) 0,D C h 1/2 p h 0,e, p h P 1 (), C h 1/2 p h 0,, p h P 1 (). Further, for all v V and µ = 0,1 there holds ( h 1 µ ) 1/2 C ( h 1 µ T v 2 µ,t) 1/2. T T h h (Ω) h (Γ N ) v 2 µ,d (2)
28 fficiency of the A posteriori rror stimator I Theorem. There exist constants γ R,γ > 0, depending only on the shape regularity of T h such that γ R η h γ osc h e 1,Ω. Proof. The proof will be given by a series of lemmas establishing the local efficiency of the estimator.
29 Lemma. Let T T h. Then there holds: h T f T Lu h 0,T C e 1,T + h T f f T 0,T. Proof. We set p h := f T. Observing ψ T T = 0, by Green s formula a T (u h,p h ψ T ) = (a u h ) p h ψ T dx + n } T a u {{ h p h ψ T} = 0 T Denoting by π h ψ T the L 2 -projection onto the linear space of elementwise constants and taking advantage of the properties of the bubble functions, it follows that ( f T Lu h 2 0,T C (f T Lu h ) π h ψ T dx = f π h ψ T dx a T (u h, π h ψ T ) + T ) (π h f f) π h ψ T dx = ( a T (e u,π h ψ T ) + T T ) (f T f) π h ψ T dx ds. T C h 1 T e u 1,T p h 0,T + f T f 0,T p h 0,T, T
30 Department of Mathematics, University of Houston Lemma. Let h (Ω). Then there holds: h 1/2 [n a u h ] 0, C e (2) 1,D + h f f T (2) 0,D D (2) D (2) + h f T Lu h (2) 0,D. Proof. We set p h := [n a u h ]. In view of ψ = 0,,, Green s formula gives n (2) D a u h p h ψ ds = a (2) D (u h,p h ψ ) + a u h }{{} p h ψ dx, = Lu h where p h ψ is the L 2 -projection onto the edgewise constants.
31 Proof cont d. If we use the properties of the bubble functions, it follows that [n (2) D a u h ] 2 0, C a u h ] p h ψ ds = a u h ] p h ψ ds = ( a D (2) (u h,p h ψ ) = a (2) D (e,p h ψ ) + C h 1/2 e u 1,D (2) D (2) ( D (2) [n D (2) f p h ψ dx + D (2) (f f T ) p h ψ dx + p h 0, + h 1/2 f f T 0,D (2) D (2) [n D (2) (f f T ) p h ψ dx + D (2) from which the assertion can be easily deduced. D (2) ) (f T Lu h ) p h ψ dx ) (f T Lu h ) p h ψ dx p h 0, + h 1/2 f T Lu h 0,D (2) p h 0,.
32 Lemma. Let h (Γ N ). Then there holds: h 1/2 g n u h 0, C e (2) 1,D + h 1/2 g g 0, + h f f T 0,D (2) + h f T Lu h (2) 0,D. Proof. We set p h := g n a u h. Observing ψ = 0,, by Green s formula we obtain n a u h p h ψ ds = a u h p h ψ ds = a (2) D (u h,p h ψ ) + D (2) D (2) n D (2) a u h }{{} = Lu h p h ψ dx.
33 Department of Mathematics, University of Houston Goal-Oriented Dual Weighted Approach
34 Goal-Oriented Dual Weighted Approach I The goal oriented dual weighted approach allows to control the error e u := u u h with respect to a rather general error functional or output functional J : V H 1 (Ω) lr. The goal oriented dual weighted approach strongly uses the solution z V of the associated dual problem a(v,z) = J(v), v V, and its finite element approximation z h V h, i.e., a(v h,z h ) = J(v h ), v h V h. Using Galerkin orthogonality, we readily deduce that J(e u ) = a(e u,z) = a(e u,z v h ) = r(z v h ), v h V h, where r( ) stands for the residual with respect to the computed finite element approximation u h.
35 Goal-Oriented Dual Weighted Approach II Theorem. Let u h V h := S 1,Γ (Ω; T h (Ω)) be the conforming P1 approximation of the solution u H 1 0(Ω) of Poisson s equation with f L 2 (Ω) and homogeneous Dirichlet boundary data. Then, the following error representation holds true J(e u ) = ) ((r T,z v h ) 0,T + (r T,z v h ) 0, T, v h V h, T T h (Ω) where the element residuals r T and the edges residuals r T are given by { 1 r T := f, T T h (Ω), r T := 2 ν [ u h ], h ( T Ω) 0, h ( T Γ) Moreover, we have the error estimate J(e u ) η DW := ω T ρ T, T T h (Ω) where for v h V h the element residuals ρ T and the weights ω T read ( 1/2, ( 1/2. ρ T := r T 2 0,T + h 1 T r T 0, T) 2 ωt := z v h 2 0,T + h T z v h 0, T) 2
36 Department of Mathematics, University of Houston Goal-Oriented Dual Weighted Approach III We remark that the previous result is not really a posteriori, since the solution z V of the dual solution is not known. Therefore, information about the weights ω T,T T h (Ω) has to be provided either by means of an a priori analysis or by the numerical solution of the dual problem. Theorem. Under the assumptions of the previous theorem let the error functional be given by Then, there holds J(v) := ( v, e u) 0,Ω e u 0,Ω, v V. e u 0,Ω C ( T T h (Ω) h 2 T ρ 2 T) 1/2.
37 Proof. The dual solution z V satisfies a(v,z) = ( v, e u) 0,Ω e u 0,Ω, v V, from which we readily deduce the a priori bound z 0,Ω 1. In view of the basic error estimate it follows that ( ) 1/2 ( J(e u ) = e u 0,Ω h 2 T ρ 2 T T T h (Ω) T T h (Ω) h 2 T ω2 T) 1/2. Choosing v h = P C z, where P C is Clément s quasi-interpolation operator, we find ( 1/2 inf (h 2 T z v h 2 0,T + h 1 T z v h 0, T) 2 C z 0,Ω. v h V h T T h (Ω) Using the last inequality in the previous one and observing the error representation gives the assertion.
38 Goal-Oriented Dual Weighted Approach IV Theorem. Consider the conforming P1 approximation of Poisson s equation under homogeneous Dirichlet boundary conditions and assume that the solution u V := H 1 0(Ω) is 2-regular. Using the the error functional J(v) := (v,e u) 0,Ω e u 0,Ω, v V, gives rise to the a posteriori error estimate ( e u 0,Ω C T T h (Ω) h 4 T ρ 2 T) 1/2.
39 Goal-Oriented Dual Weighted Approach V Finally, we apply the goal-oriented dual weighted approach to the pointwise estimation of the error at some point a Ω. Given some tolerance ε > 0, we consider the ball K ε (a) := {x Ω x a < ε} around the point a and define the regularized error functional J(v) := Kε(a) 1 v dx. Kε (a) The dual solution z of a(v,z) = J(v) behaves like a regularized Green s function With the residual ρ T we obtain z(x) log(r(x)), r(x) := x a 2 + ε 2. (u u h )(a) T T h (Ω) h 3 T r 2 T ρ T, r T := max x T r(x).
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Department of Mathematics, University of Houston Basic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationFEniCS Course. Lecture 8: A posteriori error estimates and adaptivity. Contributors André Massing Marie Rognes
FEniCS Course Lecture 8: A posteriori error estimates and adaptivity Contributors André Massing Marie Rognes 1 / 24 A priori estimates If u H k+1 (Ω) and V h = P k (T h ) then u u h Ch k u Ω,k+1 u u h
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
More informationConvergence and optimality of an adaptive FEM for controlling L 2 errors
Convergence and optimality of an adaptive FEM for controlling L 2 errors Alan Demlow (University of Kentucky) joint work with Rob Stevenson (University of Amsterdam) Partially supported by NSF DMS-0713770.
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationR T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T
2 R.H. NOCHETTO 2. Lecture 2. Adaptivity I: Design and Convergence of AFEM tarting with a conforming mesh T H, the adaptive procedure AFEM consists of loops of the form OLVE ETIMATE MARK REFINE to produce
More informationAxioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates)
Axioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates) Carsten Carstensen Humboldt-Universität zu Berlin 2018 International Graduate Summer School on Frontiers of Applied and
More informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Dietrich Braess, Carsten Carstensen, and Ronald H.W. Hoppe Convergence Analysis of a Conforming Adaptive Finite Element Method for an Obstacle
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationPDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM.
PDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM. Sergey Korotov Basque Center for Applied Mathematics / IKERBASQUE http://www.bcamath.org & http://www.ikerbasque.net 1 Introduction The analytical solution
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differential Equations Eric de Sturler University of Illinois at Urbana-Champaign Read section 8. to see where equations of type (au x ) x = f show up and their (exact) solution
More informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. NUMR. ANAL. Vol. 45, No. 1, pp. 68 82 c 2007 Society for Industrial and Applied Mathematics FRAMWORK FOR TH A POSTRIORI RROR ANALYSIS OF NONCONFORMING FINIT LMNTS CARSTN CARSTNSN, JUN HU, AND ANTONIO
More information2 Elliptic Differential Equations
2 Elliptic Differential Equations 2.1 Classical solutions As far as existence and uniqueness results for classical solutions are concerned, we restrict ourselves to linear elliptic second order elliptic
More informationA robust a-posteriori error estimate for hp-adaptive DG methods for convection-diffusion equations
IMA Journal of Numerical Analysis Page 1 of 3 doi:10.1093/imanum/drnxxx A robust a-posteriori error estimate for hp-adaptive DG methods for convection-diffusion equations LIANG ZHU AND DOMINIK SCHÖTZAU
More informationOverview. A Posteriori Error Estimates for the Biharmonic Equation. Variational Formulation and Discretization. The Biharmonic Equation
Overview A Posteriori rror stimates for the Biharmonic quation R Verfürth Fakultät für Mathematik Ruhr-Universität Bochum wwwruhr-uni-bochumde/num1 Milan / February 11th, 013 The Biharmonic quation Summary
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
More informationQUASI-OPTIMAL CONVERGENCE RATE OF AN ADAPTIVE DISCONTINUOUS GALERKIN METHOD
QUASI-OPIMAL CONVERGENCE RAE OF AN ADAPIVE DISCONINUOUS GALERKIN MEHOD ANDREA BONIO AND RICARDO H. NOCHEO Abstract. We analyze an adaptive discontinuous finite element method (ADFEM) for symmetric second
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationA WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS
A WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS LIN MU, JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This article introduces and analyzes a weak Galerkin mixed finite element method
More informationWeak Formulation of Elliptic BVP s
Weak Formulation of Elliptic BVP s There are a large number of problems of physical interest that can be formulated in the abstract setting in which the Lax-Milgram lemma is applied to an equation expressed
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationNumerical Methods for Large-Scale Nonlinear Systems
Numerical Methods for Large-Scale Nonlinear Systems Handouts by Ronald H.W. Hoppe following the monograph P. Deuflhard Newton Methods for Nonlinear Problems Springer, Berlin-Heidelberg-New York, 2004 Num.
More informationApproximation in Banach Spaces by Galerkin Methods
2 Approximation in Banach Spaces by Galerkin Methods In this chapter, we consider an abstract linear problem which serves as a generic model for engineering applications. Our first goal is to specify the
More informationAnalysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient
Analysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient Yifeng Xu 1 Jun Zou 2 Abstract Based on a new a posteriori error estimator, an adaptive finite element method is proposed
More informationThe Dirichlet-to-Neumann operator
Lecture 8 The Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator plays an important role in the theory of inverse problems. In fact, from measurements of electrical currents at the surface
More informationXingye Yue. Soochow University, Suzhou, China.
Relations between the multiscale methods p. 1/30 Relationships beteween multiscale methods for elliptic homogenization problems Xingye Yue Soochow University, Suzhou, China xyyue@suda.edu.cn Relations
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
More informationA posteriori error estimation of approximate boundary fluxes
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2008; 24:421 434 Published online 24 May 2007 in Wiley InterScience (www.interscience.wiley.com)..1014 A posteriori error estimation
More informationCONVERGENCE ANALYSIS OF AN ADAPTIVE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD FOR THE HELMHOLTZ EQUATION
CONVERGENCE ANALYSIS OF AN ADAPTIVE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD FOR THE HELMHOLTZ EQUATION A Dissertation Presented to the Faculty of the Department of Mathematics University of Houston
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationA brief introduction to finite element methods
CHAPTER A brief introduction to finite element methods 1. Two-point boundary value problem and the variational formulation 1.1. The model problem. Consider the two-point boundary value problem: Given a
More informationAn a posteriori error indicator for discontinuous Galerkin discretizations of H(curl) elliptic partial differential equations
IMA Journal of Numerical Analysis 005) Page of 7 doi: 0.093/imanum/ An a posteriori error indicator for discontinuous Galerkin discretizations of Hcurl) elliptic partial differential equations PAUL HOUSTON
More informationDownloaded 07/25/13 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMR. ANAL. Vol. 39, No. 6, pp. 2034 2044 c 2002 Society for Industrial and Applied Mathematics RSIDUAL-BASD A POSRIORI RROR SIMA FOR A NONCONFORMING RISSNR MINDLIN PLA FINI LMN CARSN CARSNSN Abstract.
More informationDomain Perturbation for Linear and Semi-Linear Boundary Value Problems
CHAPTER 1 Domain Perturbation for Linear and Semi-Linear Boundary Value Problems Daniel Daners School of Mathematics and Statistics, The University of Sydney, NSW 2006, Australia E-mail: D.Daners@maths.usyd.edu.au
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationA posteriori error estimation in the FEM
A posteriori error estimation in the FEM Plan 1. Introduction 2. Goal-oriented error estimates 3. Residual error estimates 3.1 Explicit 3.2 Subdomain error estimate 3.3 Self-equilibrated residuals 3.4
More informationFINITE ELEMENT METHODS
FINITE ELEMENT METHODS Lecture notes arxiv:1709.08618v1 [math.na] 25 Sep 2017 Christian Clason September 25, 2017 christian.clason@uni-due.de https://udue.de/clason CONTENTS I BACKGROUND 1 overview of
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More information1. Introduction. We consider the following convection-diffusion-reaction problem: Find u = u(x) such that
SIAM J. UMR. AAL. Vol. 42, o. 6, pp. 2496 2521 c 2005 Society for Industrial and Applied Mathematics XPLICIT AD AVRAGIG A POSTRIORI RROR STIMATS FOR ADAPTIV FIIT VOLUM MTHODS C. CARSTS, R. LAZAROV, AD
More informationMULTIGRID METHODS FOR MAXWELL S EQUATIONS
MULTIGRID METHODS FOR MAXWELL S EQUATIONS A Dissertation Submitted to the Graduate Faculty of the Louisiana State University and Agricultural and Mechanical College in partial fulfillment of the requirements
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationTechnische Universität Graz
Technische Universität Graz Error Estimates for Neumann Boundary Control Problems with Energy Regularization T. Apel, O. Steinbach, M. Winkler Berichte aus dem Institut für Numerische Mathematik Bericht
More informationFinite Element Methods for Fourth Order Variational Inequalities
Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2013 Finite Element Methods for Fourth Order Variational Inequalities Yi Zhang Louisiana State University and Agricultural
More informationKey words. a posteriori error estimators, quasi-orthogonality, adaptive mesh refinement, error and oscillation reduction estimates, optimal meshes.
CONVERGENCE OF ADAPTIVE FINITE ELEMENT METHOD FOR GENERAL ECOND ORDER LINEAR ELLIPTIC PDE KHAMRON MEKCHAY AND RICARDO H. NOCHETTO Abstract. We prove convergence of adaptive finite element methods (AFEM)
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationDiscontinuous Galerkin Methods: An Overview and Some Applications. Daya Reddy UNIVERSITY OF CAPE TOWN
Discontinuous Galerkin Methods: An Overview and Some Applications Daya Reddy UNIVRSITY OF CAP TOWN Joint work with Beverley Grieshaber and Andrew McBride SANUM, Stellenbosch, 22 24 March 2016 Structure
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationAn A Posteriori Error Estimate for Discontinuous Galerkin Methods
An A Posteriori Error Estimate for Discontinuous Galerkin Methods Mats G Larson mgl@math.chalmers.se Chalmers Finite Element Center Mats G Larson Chalmers Finite Element Center p.1 Outline We present an
More informationUniversité de Metz. Master 2 Recherche de Mathématiques 2ème semestre. par Ralph Chill Laboratoire de Mathématiques et Applications de Metz
Université de Metz Master 2 Recherche de Mathématiques 2ème semestre Systèmes gradients par Ralph Chill Laboratoire de Mathématiques et Applications de Metz Année 26/7 1 Contents Chapter 1. Introduction
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationAN EQUILIBRATED A POSTERIORI ERROR ESTIMATOR FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD
AN EQUILIBRATED A POSTERIORI ERROR ESTIMATOR FOR THE INTERIOR PENALTY DISCONTINUOUS GALERIN METHOD D. BRAESS, T. FRAUNHOLZ, AND R. H. W. HOPPE Abstract. Interior Penalty Discontinuous Galerkin (IPDG) methods
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationOn Dual-Weighted Residual Error Estimates for p-dependent Discretizations
On Dual-Weighted Residual Error Estimates for p-dependent Discretizations Aerospace Computational Design Laboratory Report TR-11-1 Masayuki Yano and David L. Darmofal September 21, 2011 Abstract This report
More informationTheory of PDE Homework 2
Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an
More informationWeak Solutions, Elliptic Problems and Sobolev Spaces
3 Weak Solutions, Elliptic Problems and Sobolev Spaces 3.1 Introduction In Chapter 2 we discussed difference methods for the numerical treatment of partial differential equations. The basic idea of these
More informationSTABILIZED GALERKIN FINITE ELEMENT METHODS FOR CONVECTION DOMINATED AND INCOMPRESSIBLE FLOW PROBLEMS
NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING BANACH CENTER PUBLICATIONS, VOLUME 29 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 1994 STABILIZED GALERIN FINITE ELEMENT METHODS FOR CONVECTION
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationarxiv: v1 [math.na] 19 Nov 2018
QUASINORMS IN SEMILINEAR ELLIPTIC PROBLEMS JAMES JACKAMAN AND TRISTAN PRYER arxiv:1811.07816v1 [math.na] 19 Nov 2018 (1.1) Abstract. In this note we examine the a priori and a posteriori analysis of discontinuous
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More information