Isoparametric finite element approximation of curvature on hypersurfaces
|
|
- Erick Lucas
- 6 years ago
- Views:
Transcription
1 Isoparametric finite element approximation of curvature on ypersurfaces C.-J. Heine Abstract Te discretisation of geometric problems often involves a point-wise approximation of curvature, altoug te discretised surfaces are globally only of class C 0,1. We define a discrete second fundamental form by means of an L -projection in te context of isoparametric Lagrange-elements of polynomial degree γ. We prove tat te order of convergence wit respect to te L -norm for tis discrete second fundamental form is γ 1. 1 Introduction In many interesting matematical, pysical and engineering problems te curvature of surfaces plays an important role. Often tese surfaces are free boundaries and are to be determined as one of te major unknowns; some examples are mean-curvature flow, dendrite growt or equilibrium figures of freely rotating drops. Te numerical treatment of tese problems can be difficult if te curvature enters te governing equations explicitly. Often it is possible to arrive at a weak formulation wic does no longer explicitly contain iger derivatives of te free boundary (Dziuk, 1991, for example). Wen suc a weak formulation is not known (Scmidt, 1996; Heine, 006, for example) one needs a pointwise definition of curvature for te discretised surface. Tis paper is concerned wit te approximation of te second fundamental form of a C k -ypersurface R n+1 in te context of an approximation of wit isoparametric Lagrange finite elements. Using te coordinates of te ambient space R n+1 we can write II = ν (1.1) for te Weingarten map, were ν denotes te outer unit normal field of and its tangential gradient. In Section 3 we use a weak formulation of (1.1) to define a discrete second fundamental form II for te C 0,1 -surface. Section 4 establises L -estimates for II for te case n 3. Finally Section 5 presents te results of experimental convergence tests wit respect to te L - and L -norm. Abteilung für Angewandte Matematik, Universität Freiburg, Hermann-Herder-Straße 10, Freiburg i.br., Germany. 1
2 c.-j. eine Notation and Fundamental Definitions.1 Differential geometry In te following R n+1 is an n-dimensional isometric C -immersion into R n+1. We identify te tangent space T P at a point P wit te affine yperplane of R n+1 wic is tangential to at P. We denote te (local) outer unit normal field by ν = (ν 1,..., ν n+1 ) : R n+1. For f C 1 () we define te tangential gradient of f on by f := f ( f ν ) ν R n+1, were f denotes te usual gradient of R n+1. Above definition makes sense only if f is extended to a neigbourood of. However, f does not depend on te extension but only on te values of f on. Likewise, for w = (w i ) n+1 i=1 ( C 1 () ) n+1 we define its tangential divergence by w := w ( w i ν ) ν i and we define te Laplace-Beltrami operator of g C () by g := g. Te Weingarten-map II = (II 1,..., II n+1 ) : R (n+1) (n+1) is defined by II i := ν i, 1 i n + 1. II is te usual second fundamental form but written in te coordinate frame of te ambient space R n+1, wit zero-extension to (T P ). For P te eigenvalues κ i, 1 i n, of te second fundamental form II(P ) TP T P are called principal curvatures of at P. Furter, we define te mean curvature H, te curvature vector Y and te Gauss curvature K at P by H := 1 n n κ i, Y := n H ν, K := i=1 n κ i. (.1) Elementary computation yields te formula id = Y for te curvature vector Y. Finally, for = te formula for integration by parts on is given by f do = f Y do, f C 1 (). (.). Function spaces For Ω R n open we denote by L p (Ω) and H m p (Ω) te usual Lebesgue- and Sobolevspaces wit te corresponding norms ( u Lp(Ω) = Ω u p dx ) 1 p i=1 ( m, u H m p (Ω) = D k u p L p(ω) wit te usual modification for p =. For p = we simply write H m (Ω) := H m (Ω). By C ( m,α Ω ) and C m,α (Ω) we denote te usual Hölder spaces and k=0 ) 1 p,
3 finite element approximation of curvature on ypersurfaces 3 teir associated norms. Using local coordinate-carts it is possible to carry over te definitions of Sobolev- and Hölder-spaces to a differential manifold. Te resulting spaces are independent from te special coice of te coordinate system up to norm-equivalence. In particular, Wloka (198, pp. 9) sows tat te spaces H m () are well defined for a C m 1,1 -manifold for any m N, wic covers especially te case of te C 0,1 -manifolds constructed by te finite element approac of Section.3 below. For an n-dimensional C 0,1 -immersion R n+1 we coose u H 1 () := u u do for te H 1 semi-norm..3 Isoparametric finite elements For simplicity we consider a compact closed C k, 1 -ypersurface R n+1 were k 1. We note tat te results could be extended to immersions but tis would involve some notational complications. Because te sectional curvature of is bounded tere is δ > 0 suc tat te decomposition x = π(x) + d(x) ν(x), wit π : U, is unique for x U := {x R n+1 dist(x, ) < δ}, were d(x) C k (U) is a signed distance function and ν(x) := d π(x). It follows tat ν(x) is normal to at π(x) wit ν(x) 1. Figure.1 Construction of γ π( Q i ) γ 1 Q i To construct an isoparametric finite element approximation γ of of given polynomial degree γ we first approximate by a simplicial polytope 1 wit associated triangulation T 1 = {T 1 T 1 is a face of 1 }. Te polytope 1 must be contained in te δ-strip U were te projection π is well-defined. Ten, for any simplex T 1 T 1, we take te images π(q i) of te canonical Lagrange-nodes Q i of T 1 and define an isoparametric simplex T γ by applying Lagrange-interpolation of degree γ to te coordinates of te projected Lagrange-nodes π(q i ), see Figure.1. Te result is a discrete surface γ of class C0,1 and associated isoparametric triangulation T γ were T γ := Φγ ( T ) T γ T 1 is an image of a reference simplex T R n under a polynomial parametrisation Φ γ : T T γ T 1 of degree γ. By construction 1 interpolates γ between its vertex Lagrange-nodes. Tis will be a elpful feature in Section 4 below.
4 4 c.-j. eine On γ we define parametric finite elements wic again are piece-wise polynomial of degree γ or over te reference simplex T : W γ (γ ) := {v C 0,1 ( γ ) Φγ T 1 v := (v Φ γ ) P T 1 γ ( T ) T γ T γ }. (.3) Te finite element space W γ (γ ) is virtually isoparametric as te local parametrisation of γ over te reference simplex is piecewise a polynomial of degree γ; so we ave in particular id γ (W γ (γ ))n+1. Because π γ is by construction a bijection we can lift v W γ (γ ) to v C 0,1 () by v (y) := (π 1 v )(y) = v (π 1 (y)), x γ, y = π(x). We call a family of isoparametric triangulations {T γ } >0 quasi-uniform if tere is c > 0 independent from wit c ρ(t ) T T γ wit = max (T ), T T γ were for T T γ we denote by (T ) te diameter of T and by ρ(t ) te radius of te largest ball contained in T. 3 Definition of te discrete curvature In te following always denotes a C 0,1 -approximation of a closed compact C immersion R n+1 of dimension n wit W H 1 ( ) suc tat id (W ) n+1. Definition 3.1 (Discrete Curvature Vector) Under te given assumptions te discrete curvature vector Y (W ) n+1 is defined as te solution of (Y, ψ ) L( ) = ( id, ψ ) L( ) ψ (W ) n+1. (3.1) Tis definition follows Scmidt (1993) and is motivated by te identity Y ψ do = id ψ do = id ψ do ψ (H 1 ()) n+1, wit Y := n H ν and R n+1 at least C 1, 1, see Section.1. Te experimental convergence tests for te case of a spere presented in Scmidt (1993) sow approximate linear convergence in te L and L norms for quadratic Lagrange elements and no convergence for linear elements. For te definition of te scalar mean and Gaussian curvature we ave a look at te Weingarten map II, see Section.1. To derive a weak formulation for II = (II 1,..., II n+1 ) R (n+1) (n+1) we multiply te identity II j = (ν ) j by a test-function ψ H 1 (), integrate over and take te integration by parts formula (.) on page into account. Carried over to and W tis yields:
5 finite element approximation of curvature on ypersurfaces 5 Definition 3. (Discrete Weingarten map) Te discrete Weingarten map II = (II 1,..., IIn+1 ) (W ) (n+1) (n+1) is defined as te solution of ( II j, ψ )L = ( ν j ( ), ψ )L ( ν j ( ), Y ) ψ ψ L ( ) (W ) n+1, (3.) were Y is te discrete curvature vector defined in (3.1) and ν = (ν 1,..., νn+1 ) is te (discontinuous) outer normal field of. We use II to define te discrete mean curvature H and discrete Gaussian curvature K in te obvious way, see Section.1. 4 Error estimates Teorem 4.1 Let R n+1 (n 3) be a closed compact embedded ypersurface of class C m, 1 wit m 1 and 1 γ m. Let {T γ } >0 be a quasi-uniform family of isoparametric triangulations of of polynomial degree γ. Ten it olds for te discrete curvature vector and te discrete Weingarten map Y, II W γ (γ ) tat were c = c(σ, γ, ) does not depend on. Y Y L() c γ 1, (4.1) II II L() c γ 1, (4.) Remark 4. Te regularity assumption for ensures tat Y, II H m 1, (). Remark 4.3 Our proof of Teorem 4.1 is based on Dziuk (1988) wo proves error estimates for te discretisation of te Poisson equation on arbitrary surfaces using linear finite elements. Te numerical experiments performed by Scmidt (1996) and our own numerical experiments (see Section 5) are in accordance wit tis proposition. We continue wit some intermediate results before we give a proof of Teorem 4.1. Lemma 4.4 For u W γ (γ ) it olds: Proof. See Dziuk (1988). 1 c u L( ) u L() c u L( ), 1 c u L( ) u L() c u L( ). (4.3) Lemma 4.5 (Inverse Inequality) Let T T γ. Ten it olds u L(T ) c σ(t )n/+1 u (T ) L(T γ ) u W γ (γ ). (4.4)
6 6 c.-j. eine Proof. Te proof is a one-to-one copy of te proof for non-parametric finite elements. Following Dziuk (1988) we introduce te following abbreviations: P u := u u ν u ν, P ik := δ ik ν i ν k, H ik := d xix k = ν i,xk = ν k,xi. We ave tat u (x) = (P (x) d(x) H(x)) u x d(x) ν(x) and P H = H P = H. Te defining equations (3.1) for Y can be rewritten using integrals over instead of : 1 Y ϕ do = (P (I d H) id ) : (P (I d H) ϕ ) 1 do, (4.5) wit := do do. Obviously we ave 1 c c <. By defining te abbreviation A := 1 P (I d H)P (I d H)P we arrive at 1 ( ) (4.5) Y ϕ do = A id : ϕ do. (4.6) Alternatively, we can use te following formulation: 1 1 Y ϕ do = (P (I d H) ) ϕ do = B : ϕ do, (4.7) wit B := 1 P (I d H)P. To cope wit te error introduced by lifting W γ (γ ) to we need te following estimates: Lemma 4.6 Let R n+1 be a closed compact embedded ypersurface of class C m, 1 wit m 1. Let T γ be a quasi-uniform isoparametric triangulation of of polynomial degree 1 γ m. Ten were c = c(σ, γ, ). d c γ+1, (4.8) 1 c γ+1, (4.9) (A I) P c γ+1, (4.10) (B I) P c γ, (4.11) (4.8) and (4.9) are already mentioned witout proof in Nédélec (1976) for surfaces in R 3. Dziuk (1988) proves (4.10) for te linear case (γ = 1) for surfaces in R 3. We neverteless give a detailed proof of (4.9) because we tink tat te result is not obvious, toug te calculations are elementary. Proof Lemma 4.6. For any isoparametric surface γ we also ave naturally a simplicial polytope 1 wic interpolates te Cm, 1 -surface as well as γ linearly in te vertex nodes of te triangulation T γ. Let T γ T γ be an isoparametric
7 finite element approximation of curvature on ypersurfaces 7 element, T = π(t γ ) its associated lifted element and T 1 te corresponding triangle wic interpolates T γ and T in te vertex Lagrange-nodes (Figure 4.1). Witout loss of generality we can assume tat T 1 Rn {0}. We use te following notations for te coordinates on R n+1, T 1, T γ and T : R n+1 : x = (x 1,..., x n+1 ). T 1 : η = (η 1,..., η n ). Witout loss of generality we can assume η i = x i, 1 i n. We parametrise T γ and T over T 1. T γ : ξ = (ξ 1(η),..., ξ n+1 (η)). Te coordinate functions ξ i (η) of T γ are polynomials of degree γ over T 1. T : π(ξ(η)) = (π 1 ξ(η),..., π n+1 ξ(η)) = ξ(η) d(ξ(η)) ν(ξ(η)). We remark tat te coordinate functions ξ i (η) of γ are ordinary (i.e. nonparametric) finite element interpolations of polynomial degree γ of te functions π i ξ(η) H m+1, (T 1 ). Tis means tat te usual estimates for Lagrange-interpolation apply wit respect to te interpolation of π i by ξ i. Equation (4.8) is ten clear because c γ+1 π i (ξ(η)) ξ i (η) = d(ξ(η)) ν i (ξ(η)), 1 i (n + 1). Te estimate (4.9) for te surface element is more complicated. For any (n + 1)- tuple (a 1,..., a n+1 ) we denote by a b i te n-tuple were te i-t component as been omitted (e.g. a b 1 = (a,..., a n+1 )). In tis notation te surface elements do on T γ and do on T are n do = π n+1 + π bk n dη, do = x n+1 + ξ bk dη. So we ave 1 do do k=1 x n+1 π n+1 + n ξ bk k=1 k=1 π bk. (4.1) Figure 4.1 Parametrisation of over T γ over T 1 π (ξ(η)) d( ξ ( η)) ν( ξ( η)) T ξ(η) T γ η T 1
8 8 c.-j. eine For estimating te individual determinants we employ Leibniz formula, using te notation k i := i for i < k, k i := i + 1 for i k: ξ bk = n ξ ki sgn(α), α S n i=1 α(i) π bk = n sgn(α) α S n i=1 π ki α(i). (4.13) First we note tat T 1 interpolates T γ as well as T linearly between teir vertices. So we ave te following interpolation estimates for te derivatives of te coordinate functions (1 j n): c ξ i i = ξ i δ ij, c π i δ ij (1 i n), (4.14) c ξ n+1, c π n+1. (4.15) It follows tat ξ b i c, π b i c (1 i n), (4.16) because eac of te summands in Leibniz formula (4.13) contains a factor ξn+1 respectively πn+1. Furter, te polynomials ξ i (η) interpolate π i (η), so we ave π k i ξ bk ξ k i c γ. Togeter wit (4.16) tis yields π bk It remains to estimate ξ n+1 = ξ bk π bk ξ bk + π bk c γ+1 (1 k n). }{{}}{{} =O( γ ) =O() π n+1 = ( + O()) (4.14) ξ n+1 (4.17) π n+1. (4.18) To find an estimate we ave to consider te differences ξi πi = d ν i d νi for 1 i, j n. From te approximation properties of ξ i (η) we first get c γ ξ k π k = d ν k + d ν k (1 k n + 1) = d c γ a.e., because d = O( γ+1 ) and ν k is bounded as is C 1, 1. (4.19) Furter, we ave ν i = O() for 1 i n. To see tis we apply te cain rule to d(ξ(η)) i = d x n+1 ξ(η) ξ n+1 i + n j=1 d x j ξ(η) ξ j i
9 finite element approximation of curvature on ypersurfaces 9 So we get ν i = d x i = O(), 1 i n, by estimating eac term according to (4.14), (4.15) and (4.19). Terefore it olds because νi yields ξ i π i = d ν i d ν i = O( γ+1 ), 1 i n, (4.0) is bounded as is at least C 1, 1. Togeter wit (4.13) and (4.18) tis ξ n+1 π n+1 = c x n+1 π n+1 c γ+1, wic completes te proof of (4.9). We quote te proof for (4.10) from Dziuk (1988) were te result is proved for te linear case (γ = 1). Te generalisation for γ > 1 is just a consequence of (4.9): (A I) P = 1 P ((I d H) P (I d H) P I) P (4.9) = P (I d H) P (I d H) P I) P + O( γ+1 ) (4.8) = P P P P + O( γ+1 ) = P ν ν t P + O( γ+1 ). It follows from (4.17) and (4.0) tat ν = ν + O( γ ). So we ave For (4.11) we similarly get (A I) P P ν + c γ+1 c γ+1, (B I) P = 1 P ((I d H) P I) P (4.9) = P (I d H) P I) P + O( γ+1 ) (4.8) = P P P + O( γ+1 ) = ν ν t P + O( γ+1 ) = O( γ ). Proof Teorem 4.1. Let ϕ W γ (γ ) be an arbitrary test-function. For brevity we set := L() and := L (). We ave Y Y (Y Y ) (ϕ Y ) do } {{ } :=I 1 + (Y Y ) (Y ϕ ) do. (4.1) } {{ } :=I
10 10 c.-j. eine To use te weak formulation (4.6) wit I 1 we ave to compensate for te factor 1 introduced by projecting te finite element functions to te smoot surface : I 1 = (4.7) (Y + Y Y Y ) (ϕ Y ) do (Y Y ) (ϕ Y ) do + 1 Y ϕ Y + Lem. 4.6 c γ+1 ϕ Y + c γ (ϕ Y ) 1 Lem. 4.5 c γ 1 ϕ Y. (Y 1 Y ) (ϕ Y ) do (I B ) P : (ϕ Y ) do So we finally get from (4.1): Y Y Y Y Y ϕ + c γ 1 ϕ Y ϕ W γ (γ ). = Y Y c γ 1, by coosing ϕ = I γ 1 Y, using tat Y H m 1, (). Tis proves equation (4.1) of Teorem 4.1 on page 5. Te proof of te estimate (4.) for te discrete second fundamental form II works very similar. As for te proof of (4.1) we start by transforming te integrals in te defining equation (3.) for II into surface integrals over ; we get j 1 II j ϕ do = ν ( A (id ) ) : ϕ do ν j 1 Y ϕ do. (4.) Terefore j II j II j j j (II j II ) (IIj ϕ ) do + (II j II ) (ϕ II ) do j j II j II II j ϕ + 1 II j ϕ II + (II j 1 j j II ) (ϕ II ). µ } {{} :=I 3 (4.3)
11 finite element approximation of curvature on ypersurfaces 11 We transform I 3 by employing te weak formulation (4.): I 3 (ν j id ν j j A id ) : ( ϕ II ) do + (ν j Y ν j j Y ) (ϕ II ) do ( (A I) P id + (id id ) ) (ϕ II j ) + ν j ν j id (ϕ II j ) + µ 1 Y Y II j ϕ + ( µ 1 ν j ν j Y + 1 µ 1 ) j Y II ϕ c γ 1 II j ϕ (see (4.1), (4.9), (4.10), Lemma 4.5). Te remaining terms in (4.3) are of iger order so tat we finally arrive at j II j II c γ 1, wic completes te proof of Teorem 4.1 on page 5. 5 Numerical experiments We ave implemented our definition of te discrete curvature (see Section 3) for isoparametric finite elements in te context of te finite element toolbox ALBERT (Scmidt and Siebert, 1998). Table 5.1 to 5.4 on pages sow te results of our experimental convergence tests wit piecewise polynomial discretisations of degree 1,, 3 and 4 for te following model problems (see Figure 5.1 on te next page): (x z ) + y + z 1 = 0, (5.1a) (x z ) + (y z ) + z 1 = 0, (5.1b) x + ( y) + ( 3 z) 1 = 0, x + y + z 1 = 0. (5.1c) (5.1d) Table 5. on page 14 corresponding to surface (5.1b) as already been publised in Heine (006). Te underlying macro triangulation for te discretisation of all model-problems is an icosaedron. Te discrete surface is just te interpolation of te exact sape wit Lagrange-polynomials. At first te mes was refined adaptively until te L error wit respect to te exact surface was below Afterwards te meses were refined globally by doubly bisecting te elements in eac refinement step. Te tables 5.1 to 5.4 on pages sow te maximal mes widt, te relative L and L errors and te experimental order of convergence (EOC) of te discrete mean and Gaussian curvatures H and K derived from te discrete second fundamental form as defined in Definition 3. on page 5.
12 1 c.-j. eine Figure 5.1 Surface meses for testing our definition of te discrete curvature. Te colour indicates te magnitude of te discrete Gaussian curvature. Te defining equations for te surfaces are (5.1b)-(5.1d). A picture for (5.1a) as been omitted as it is of te same type as te surface sown in figure (a). (a) Equation (5.1b) (b) Equation (5.1c) (c) Equation (5.1d) For te normalisation of te relative errors we use te curvature of te exact surfaces defined by (5.1a) to (5.1d): E C,r := C C Lr( ), were C {K, H}, r {, }. C Lr( ) In order to compare functions defined on wit tose defined on we were using an approximate ortogonal lift from onto wit respect to an averaged normal field on. For a sequence of isoparametric triangulations T j te experimental order of convergence (eoc) is defined by EOC j := (log Ej E j+1 )/(log j j+1 ), were E j := C C j denotes te error between te discrete solution C j and te exact solution C in te appropriate norm. Wile it is not surprising tat te discrete curvatures originating from Definition 3. on page 5 do not converge for linearly interpolated surfaces, te experimental order of convergence for iger order parametrisations is surprising. Te tables 5.1 to 5.4 on pages suggest tat for Lagrange-interpolates of polynomial degree k te experimental order of convergence for te discrete curvature is greater tan k 1, provided tat k > 1. We would ave expected values for te EOC of at most k 1, but not greater. However, one notes tat te results for te less complicated surfaces (5.1c) an ellipsoid and (5.1d) a spere are not as suggestive as te results for te more complicated surfaces: for tem table 5.3 on page 15 and 5.4 on page 16 sow a decreasing EOC as te number of unknowns increases. Tis migt indicate tat te values for te EOC as sown in te tables 5.1 on te next page and 5. on page 14 just did not stabilise yet.
13 finite element approximation of curvature on ypersurfaces 13 Table 5.1 Discrete curvature, surface (5.1a), relative errors Polynomial degree: Polynomial degree: Polynomial degree: Polynomial degree: e e e
14 14 c.-j. eine Table 5. Discrete curvature, surface (5.1b), relative errors Polynomial degree: Polynomial degree: Polynomial degree: Polynomial degree:
15 finite element approximation of curvature on ypersurfaces 15 Table 5.3 Discrete curvature, surface (5.1c), relative errors Polynomial degree: Polynomial degree: Polynomial degree: Polynomial degree: e e e e e e e 05 3.
16 16 c.-j. eine Table 5.4 Discrete curvature, surface (5.1d), relative errors Polynomial degree: Polynomial degree: e e e e e e e 05. Polynomial degree: e e e e Polynomial degree: e e e e e e e e e e e e e e
17 finite element approximation of curvature on ypersurfaces 17 References G. Dziuk. Finite elements for te Beltrami operator on arbitrary surfaces. In S. Hildebrandt and R. Leis, editors, Partial Differential Equations & Calculus of Variations, volume 1357 of Lecture Notes in Mat., pages Springer, Berlin, G. Dziuk. An algoritm for evolutionary surfaces. Numer. Mat., 58(6): , Te algoritm uses Dziuk (1988). Claus-Justus Heine. Computations of form and stability of rotating drops wit finite elements. IMA J. Numer. Anal., 6(4):73 751, 006. ISSN J.-C. Nédélec. Curved finite element metods for te solution of singular integral equations on surfaces in R 3. Comput. Metods Appl. Mec. Engrg., 8(1):61 80, A. Scmidt. Computation of tree dimensional dendrites wit finite elements. Journal of Computational Pysics, 15:93 31, A. Scmidt. Die Berecnung dreidimensionaler Dendriten mit Finiten Elementen. PD tesis, Albert-Ludwigs-Universität Freiburg i. Br., A. Scmidt and K.G. Siebert. ALBERT 0. Adaptive multilevel finite element toolbox using Bisectioning refinement and Error control by Residual Tecniques, August Scmidt and Siebert (1999). A. Scmidt and K.G. Siebert. Abstract data structures for a finite element package: Design principles of ALBERT. Zeitscrift für Angewandte Matematik und Mecanik, 79:49 5, J. Wloka. Partielle Differentialgleicungen. B.G. Teubner, Stuttgart, 198.
APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS. Introduction
Acta Mat. Univ. Comenianae Vol. LXVII, 1(1998), pp. 57 68 57 APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS A. SCHMIDT Abstract. Te pase transition between solid and liquid in an
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationTHE DISCRETE PLATEAU PROBLEM: CONVERGENCE RESULTS
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 0025-5718XX0000-0 THE DISCRETE PLATEAU PROBLEM: CONVERGENCE RESULTS GERHARD DZIUK AND JOHN E. HUTCHINSON Abstract. We solve te problem of
More informationError estimates for a semi-implicit fully discrete finite element scheme for the mean curvature flow of graphs
Interfaces and Free Boundaries 2, 2000 34 359 Error estimates for a semi-implicit fully discrete finite element sceme for te mean curvature flow of graps KLAUS DECKELNICK Scool of Matematical Sciences,
More informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationMA455 Manifolds Solutions 1 May 2008
MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationFINITE ELEMENT APPROXIMATIONS AND THE DIRICHLET PROBLEM FOR SURFACES OF PRESCRIBED MEAN CURVATURE
FINITE ELEMENT APPROXIMATIONS AND THE DIRICHLET PROBLEM FOR SURFACES OF PRESCRIBED MEAN CURVATURE GERHARD DZIUK AND JOHN E. HUTCHINSON Abstract. We give a finite element procedure for te Diriclet Problem
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationAN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS
Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationAnalysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem
Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationPriority Program 1253
Deutsce Forscungsgemeinscaft Priority Program 53 Optimization wit Partial Differential Equations K. Deckelnick, A. Günter and M. Hinze Finite element approximation of elliptic control problems wit constraints
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationBOUNDARY REGULARITY FOR SOLUTIONS TO THE LINEARIZED MONGE-AMPÈRE EQUATIONS
BOUNDARY REGULARITY FOR SOLUTIONS TO THE LINEARIZED MONGE-AMPÈRE EQUATIONS N. Q. LE AND O. SAVIN Abstract. We obtain boundary Hölder gradient estimates and regularity for solutions to te linearized Monge-Ampère
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationFinite element methods for surface PDEs
Acta Numerica (2013), pp. 289 396 c Cambridge University Press, 2013 doi:10.1017/s0962492913000056 Printed in te United Kingdom Finite element metods for surface PDEs Gerard Dziuk Abteilung für Angewandte
More informationarxiv: v2 [math.na] 5 Jul 2017
Trace Finite Element Metods for PDEs on Surfaces Maxim A. Olsanskii and Arnold Reusken arxiv:1612.00054v2 [mat.na] 5 Jul 2017 Abstract In tis paper we consider a class of unfitted finite element metods
More informationError analysis of a finite element method for the Willmore flow of graphs
Interfaces and Free Boundaries 8 6, 46 Error analysis of a finite element metod for te Willmore flow of graps KLAUS DECKELNICK Institut für Analysis und Numerik, Otto-von-Guericke-Universität Magdeburg,
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationDoubly Warped Products in Locally Conformal Almost Cosymplectic Manifolds
INTERNATIONAL ELECTRONIC JOURNAL OF GEOMETRY VOLUME 0 NO. 2 PAGE 73 8 207) Doubly Warped Products in Locally Conformal Almost Cosymplectic Manifolds Andreea Olteanu Communicated by Ion Miai) ABSTRACT Recently,
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationChapter 1. Density Estimation
Capter 1 Density Estimation Let X 1, X,..., X n be observations from a density f X x. Te aim is to use only tis data to obtain an estimate ˆf X x of f X x. Properties of f f X x x, Parametric metods f
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationLecture 21. Numerical differentiation. f ( x+h) f ( x) h h
Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationA New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions
Numerisce Matematik manuscript No. will be inserted by te editor A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions Wang Ming 1, Zong-ci Si 2, Jincao Xu1,3 1 LMAM, Scool of Matematical
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More informationDifferentiation. Area of study Unit 2 Calculus
Differentiation 8VCE VCEco Area of stud Unit Calculus coverage In tis ca 8A 8B 8C 8D 8E 8F capter Introduction to limits Limits of discontinuous, rational and brid functions Differentiation using first
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationarxiv: v1 [math.na] 11 May 2018
Nitsce s metod for unilateral contact problems arxiv:1805.04283v1 [mat.na] 11 May 2018 Tom Gustafsson, Rolf Stenberg and Jua Videman May 14, 2018 Abstract We derive optimal a priori and a posteriori error
More informationEffect of Numerical Integration on Meshless Methods
Effect of Numerical Integration on Mesless Metods Ivo Babuška Uday Banerjee Jon E. Osborn Qingui Zang Abstract In tis paper, we present te effect of numerical integration on mesless metods wit sape functions
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationA CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Number 3, September 2003 A CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION GH. MICULA, E. SANTI, AND M. G. CIMORONI Dedicated to
More informationAn approximation method using approximate approximations
Applicable Analysis: An International Journal Vol. 00, No. 00, September 2005, 1 13 An approximation metod using approximate approximations FRANK MÜLLER and WERNER VARNHORN, University of Kassel, Germany,
More informationarxiv: v1 [math.na] 20 Jul 2009
STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv:0907.3438v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationRIGIDITY OF TIME-FLAT SURFACES IN THE MINKOWSKI SPACETIME
RIGIDITY OF TIME-FLAT SURFACES IN THE MINKOWSKI SPACETIME PO-NING CHEN, MU-TAO WANG, AND YE-KAI WANG Abstract. A time-flat condition on spacelike 2-surfaces in spacetime is considered ere. Tis condition
More informationarxiv: v3 [math.na] 25 Jun 2017
ANALYSIS OF A HIGH ORDER UNFITTED FINITE ELEMENT METHOD FOR ELLIPTIC INTERFACE PROBLEMS CHRISTOPH LEHRENFELD AND ARNOLD REUSKEN arxiv:1602.02970v3 [mat.na] 25 Jun 2017 Abstract. In te context of unfitted
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
More informationWeierstraß-Institut. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stocastik im Forscungsverbund Berlin e.v. Preprint ISSN 0946 8633 Stability of infinite dimensional control problems wit pointwise state constraints Micael
More informationOn convexity of polynomial paths and generalized majorizations
On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationSome Error Estimates for the Finite Volume Element Method for a Parabolic Problem
Computational Metods in Applied Matematics Vol. 13 (213), No. 3, pp. 251 279 c 213 Institute of Matematics, NAS of Belarus Doi: 1.1515/cmam-212-6 Some Error Estimates for te Finite Volume Element Metod
More informationOverlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions
Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationMinimal surfaces of revolution
5 April 013 Minimal surfaces of revolution Maggie Miller 1 Introduction In tis paper, we will prove tat all non-planar minimal surfaces of revolution can be generated by functions of te form f = 1 C cos(cx),
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol 4, No, pp 86 84 c 004 Society for Industrial and Applied Matematics LEAST-SQUARES METHODS FOR LINEAR ELASTICITY ZHIQIANG CAI AND GERHARD STARKE Abstract Tis paper develops least-squares
More informationAdaptive Finite Element Method
39 Capter 3 Adaptive inite Element Metod 31 Introduction As already pointed out in capter 2, singularities occur in interface problems Wen discretizing te problem (221) wit inite Elements, te singularities
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationarxiv: v1 [math.dg] 4 Feb 2015
CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola
More informationContinuity. Example 1
Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)
More informationPOLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY
APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky
More informationON THE CONSISTENCY OF THE COMBINATORIAL CODIFFERENTIAL
TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 366, Number 10, October 2014, Pages 5487 5502 S 0002-9947(2014)06134-5 Article electronically publised on February 26, 2014 ON THE CONSISTENCY OF
More information1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist
Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x
More informationCONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION
CONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION SÖREN BARTELS AND ANDREAS PROHL Abstract. Te Landau-Lifsitz-Gilbert equation describes dynamics of ferromagnetism,
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationA SHORT INTRODUCTION TO BANACH LATTICES AND
CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationMath 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0
3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More informationOn the best approximation of function classes from values on a uniform grid in the real line
Proceedings of te 5t WSEAS Int Conf on System Science and Simulation in Engineering, Tenerife, Canary Islands, Spain, December 16-18, 006 459 On te best approximation of function classes from values on
More information1 Introduction to Optimization
Unconstrained Convex Optimization 2 1 Introduction to Optimization Given a general optimization problem of te form min x f(x) (1.1) were f : R n R. Sometimes te problem as constraints (we are only interested
More informationA finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations
ANZIAM J. 54 (CTAC2012) pp.c681 C698, 2013 C681 A finite element approximation for te quasi-static Maxwell Landau Lifsitz Gilbert equations Kim-Ngan Le 1 T. Tran 2 (Received 31 October 2012; revised 29
More informationBlanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS
Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating
More informationSection 15.6 Directional Derivatives and the Gradient Vector
Section 15.6 Directional Derivatives and te Gradient Vector Finding rates of cange in different directions Recall tat wen we first started considering derivatives of functions of more tan one variable,
More informationADAPTIVE MULTILEVEL INEXACT SQP METHODS FOR PDE CONSTRAINED OPTIMIZATION
ADAPTIVE MULTILEVEL INEXACT SQP METHODS FOR PDE CONSTRAINED OPTIMIZATION J CARSTEN ZIEMS AND STEFAN ULBRICH Abstract We present a class of inexact adaptive multilevel trust-region SQP-metods for te efficient
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationRecent Progress in the Integration of Poisson Systems via the Mid Point Rule and Runge Kutta Algorithm
Recent Progress in te Integration of Poisson Systems via te Mid Point Rule and Runge Kutta Algoritm Klaus Bucner, Mircea Craioveanu and Mircea Puta Abstract Some recent progress in te integration of Poisson
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
INTODUCTION ND MTHEMTICL CONCEPTS PEVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips of sine,
More informationCOMPUTATIONAL COMPARISON BETWEEN THE TAYLOR HOOD AND THE CONFORMING CROUZEIX RAVIART ELEMENT
Proceedings of ALGORITMY 5 pp. 369 379 COMPUTATIONAL COMPARISON BETWEEN E TAYLOR HOOD AND E CONFORMING OUZEIX RAVIART ELEMENT ROLF KRAHL AND EBERHARD BÄNSCH Abstract. Tis paper is concerned wit te computational
More informationSmoothed projections in finite element exterior calculus
Smooted projections in finite element exterior calculus Ragnar Winter CMA, University of Oslo Norway based on joint work wit: Douglas N. Arnold, Minnesota, Ricard S. Falk, Rutgers, and Snorre H. Cristiansen,
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationMultigrid Methods for Discretized PDE Problems
Towards Metods for Discretized PDE Problems Institute for Applied Matematics University of Heidelberg Feb 1-5, 2010 Towards Outline A model problem Solution of very large linear systems Iterative Metods
More informationIntegral Calculus, dealing with areas and volumes, and approximate areas under and between curves.
Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral
More informationarxiv: v1 [math.na] 27 Jan 2014
L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates
More informationarxiv: v1 [math.na] 20 Nov 2018
An HDG Metod for Tangential Boundary Control of Stokes Equations I: Hig Regularity Wei Gong Weiwei Hu Mariano Mateos Jon R. Singler Yangwen Zang arxiv:1811.08522v1 [mat.na] 20 Nov 2018 November 22, 2018
More information