FINITE ELEMENT APPROXIMATIONS AND THE DIRICHLET PROBLEM FOR SURFACES OF PRESCRIBED MEAN CURVATURE
|
|
- William Lambert
- 5 years ago
- Views:
Transcription
1 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 corresponding to surfaces of prescribed mean curvature and prove an optimal convergence estimate in te H 1 -norm. 1. H-Harmonic Maps Te numerical solution of te classical H-Plateau Problem consists of approximating disc-like surfaces wit prescribed boundary curve and prescribed mean curvature H. For a detailed discussion of te algoritms and teory see [6] for te case of zero mean curvature, and [7] for te constant mean curvature case. In tis paper we consider te associated H-Diriclet problem. Estimates for finite element approximations to solutions of general nonlinear elliptic systems are obtained in [4], using a continuity metod involving L estimates for te discrete problem. Here we give a muc sorter proof of te H 1 estimate, avoiding te need for L estimates and only assuming te discrete and smoot data are close in te H 1 sense. Our tecniques apply to a wide class of nonlinear systems. We treat te case of a non-polygonal and non-convex boundary and give te explicit dependence on te non-degeneracy constant of te smoot solution being approximated. Te arguments are prototypes of tose used in [7] for treating te more difficult case of te (free boundary) H-Plateau Problem. Te main tool for avoiding L norms in te present borderline case is te isoperimetric inequality due to Rado, see Remark 3.1. Tis work was partially supported by te Australian Researc Council. Trougout, ( R 2 ) is a bounded domain wit C 2 boundary. Function spaces will consist of functions defined over wit values in R 3 unless oterwise clear from context. Constants will depend on and oter quantities as indicated Matematics Subject Classification. Primary: 65N30; Secondary: 49Q05, 53A10. Key words and prases. finite elements, mean curvature, H-surface, H-armonic. 1
2 2 GERHARD DZIUK AND JOHN E. HUTCHINSON By H 1 is meant te H 1 seminorm, and by H 1 te full norm. Note tat by Poincaré s inequality, H 1 () is a norm on H 1 0(). For vectors a, b, c R 3, te triple product is defined by [a, b, c] =a b c. Tis is invariant under cyclic permutations of a, b and c, and antisymmetric wit respect to intercanging any two. It is te volume of te parallelopiped spanned by a, b and c. Definition 1.1. Suppose H is a real number. A function u H 2 (; R 3 ) is H-armonic wit boundary data u 0 H 2 (; R 3 )if (1.1) (1.2) u=2hu x u y u = u 0 a.e. in on Example 1.1. Let D be te closed unit disc in R 2. Let u 0 (x, y) =(x, y, 0): D R 3 wit 0 <H<1. Tere are two solutions of (1.1) and (1.2) obtained by mapping te unit disc D conformally, i.e. stereograpically projecting from a suitable point, onto te lower sperical caps obtained from eac of te two speres of radius 1/H (mean curvature H) wic contain te image of u 0. Tese solutions are called small or large depending on weter teir images do not, or do, contain a emispere. (We use tis example for test computations, see Tables 1 and 2.) If 1 <H< 0 ten one similarly obtains two solutions from te upper sperical caps. If H = 0 ten tere is exactly one solution, te map u(x, y) = (x, y, 0): D R 3. If H = 1 ten one obtains a solution by mapping onto te lower emispere of a spere of radius 1, and onto te upper emispere if H = 1. Equation (1.1) is te Euler-Lagrange system associated to te H- Diriclet integral D H (u) =D H (u;)= 1 (1.3) u 2 +2HV (u), 2 were V (u) =V(u;):= 1 (1.4) [u, u x,u y ] 3 can be tougt of as te signed volume of te cone over te origin obtained from te image of u. In fact, direct computation and integration by parts easily gives D H (u),ϕ = D H (u;),ϕ := d D dt H (u + tϕ) (1.5) t=0 = u ϕ +2H [ϕ, u x,u y ]
3 PRESCRIBED MEAN CURVATURE DIRICHLET PROBLEM 3 for u C 2 ( ; R 3 ) and ϕ C0( ; 2 R 3 ), and ence for u H 1 L (; R 3 ) and ϕ H0 1 L (; R 3 ) by a limit argument, for example see [10, Remark III.1.1]. If u H 1 L is stationary for D H, i.e. (1.6) u ϕ +2H [ϕ, u x,u y ]=0 for all ϕ H0 1 L, ten u is said to be a weak solution of (1.1). Example 1.1 is fairly typical. Arguing euristically, te energy functional D H (u) is cubic in u and tus one expects (generically) eiter two or no stationary points. In te former case one expects te smaller solution to be a local minimum and te larger solution to be unstable. Indeed, one as te following result due to te combined work of Heinz, Werner, Hildebrandt, Jäger, Wente, Brezis Coron, Struwe and Steffen. For detailed references see Struwe [10, 11]. Teorem 1.1. Assume u 0 H 1 L (; R 3 ) and H R satisfy u 0 L H 1. Ten tere exists u u 0 + H0 1 L suc tat D H (u) = min { D H (v) :v u 0 +H 1 0, v L H 1 }. Moreover, u L u 0 L ( ) and u is a weak solution to (1.1) and (1.2). If furtermore u 0 L H < 1 ( ) ten u is te unique local minimum of D H in u 0 + H0 1 L. Moreover, u is te unique weak solution of (1.1) and (1.2) wic satisfies ( ). Te function u is called te small solution of (1.1) and (1.2). Under te same assumption ( ) if H 0and u 0 is not constant, tere is also a second weak solution ū to (1.1) and (1.2) wic satisfies ū L > u 0 L. Any suc solution is called a large solution to (1.1) and (1.2). If u 0 H 2 (, R 3 ) ten any weak solution to (1.1) and (1.2) belongs to H 2 (, R 3 ). Remark Te large solution need not be unique, altoug one would expect tat tis is te generic situation. An example of Wente, [11, Example IV.3.7], gives a continuum of solutions for te unit disc and boundary data u 0 (x, y) =(x, 0, 0). See Fig. 1 for te image of te trivial small solution and of one of te large solutions on a relatively coarse grid. Rotation of u() around te u 1 -axis gives a continuum of solutions.
4 4 GERHARD DZIUK AND JOHN E. HUTCHINSON Figure 1. Wente s example (discrete approximations); small and one of a continuum of large solutions 2. Te existence of a large solution is obtained by a mountain pass type argument, see [2] and [11, Teorem III.4.8] Remark 1.2 (Nondegeneracy). We will be interested in approximating functions u H 2 (; R 3 ) wic are H-armonic and nondegenerate in te sense tat te second variation D H (u) as no zero eigenvalues. Tis is always true for small solutions, see [11, Lemma III.4.7]. More precisely, for u H 2 (; R 3 ) and ϕ, ψ H0(; 2 R 3 ) one first easily cecks by direct computation and integration by parts tat D H (u)(ϕ, ψ) =D H (u; )(ϕ, ψ) := 2 D s t H (u + tϕ + sψ) s=t=0 (1.7) = ϕ ψ +2H [u, ϕ x,ψ y ]+[u, ψ x,ϕ y ] (1.8) = ϕ ψ +2H [ψ, u x,ϕ y ]+[ψ, ϕ x,u y ], see [10, Remark III.1.1] and te paragrap following (3.2). From (1.8) D H (u) extends to a bounded symmetric bilinear functional on H 1 0, since u ϕ ϕ u L 4 ϕ L 4 ϕ L 2 c u H 2 ϕ 2 L 2 It follows tat te inner product H 1 induces a bounded self-adjoint linear operator 2 D H (u): H0 1 H0. 1 Te eigenvalues of 2 D H (u) are real, bounded below by λ 0 (say) and ave no accumulation point. Moreover, λ 0 = λ 0 (,H, u H 2), as follows from using (1.8) to estimate te Raleig-Ritz quotient. Te nondegeneracy constant λ of D H (u) is defined by λ = min { γ : γ is an eigenvalue of 2 D H (u) }
5 PRESCRIBED MEAN CURVATURE DIRICHLET PROBLEM 5 Ten 2 D H (u) is one-one and onto iff λ>0. Let (1.9) ϕ = ϕ + + ϕ denote te H 1 ortogonal decomposition of ϕ H 1 0 into members of te positive and negative spaces H + and H corresponding to te eigenvalues and eigenfunctions of 2 D H (u). Ten (1.10) D H (u)(ϕ, ϕ + ϕ ) λ ϕ 2 H 1 for all suc ϕ and λ is te largest suc real. From te eigenfunction equation, c.f. (1.8), togeter wit te estimate for λ 0, one obtains ϕ H 2 and (1.11) ϕ H 2 ν ϕ H 1 were ν = ν(,h, u H 2,d) wit d te dimension of H. 2. Discrete H-Harmonic Maps For >0 let T be a triangulation of by triangles T wose side lengts are bounded above by c for some c independent of and wose interior angles are bounded away from zero uniformly and independently of. Te intersection of any two different triangles is eiter empty, a common vertex, or a common edge. Let = T. T T Let X = { } u C 0 ( ; R 3 ):u T P 1 (T) T T, X 0 ={ϕ X :ϕ =0}, were P 1 (T ) is te set of polynomials over T of degree at most one. For some δ>0 and all sufficiently small, := { x R 2 : d(x, ) <δ}. If u H 2 () ten by te C 2 regularity of tere exists an extension of u to, also denoted by u, suc tat (2.1) u H 2 ( ) c u H 2 (). Definition 2.1. Te discrete H-Diriclet integral is defined by D H (u ; )= 1 u 2 +2HV (u ; ) 2 for u X.
6 6 GERHARD DZIUK AND JOHN E. HUTCHINSON It follows from (1.5) and (1.7) wit replaced by and a limit argument, or by direct computation and noting tat boundary integrals on internal edges cancel, tat (2.2) D H(u ; ),ϕ = u ϕ +2H [ϕ,u x,u y ], D H (u ; )(ϕ,ψ )= (2.3) ϕ ψ +2H [u,ϕ x,ψ y ]+[u,ψ x,ϕ y ], for u X and ϕ,ψ X 0. Motivated by (1.6), one as Definition 2.2. A function u X is discrete H-armonic if (2.4) u ϕ +2H [u x,u y,ϕ ]=0, for all ϕ X 0. We will prove te following. Teorem 2.1. Let u H 2 (; R 3 ) be H-armonic and u = u 0 on were u 0 H 2 (; R 3 ). Assume u is nondegenerate wit nondegeneracy constant λ. Let u 0 X and assume u 0 u 0 H 1 ( ) α. Ten tere exist constants 0 = 0 ( u H 2, u 0 H 2,α,,d,H,λ),ε 0 = ε 0 (H, λ), and c 0 = c 0 ( u H 2, u 0 H 2,α,H) suc tat if 0 < 0 ten: 1. Tere exists a unique discrete H-armonic function u suc tat u = u 0 on and 2. Moreover, u u H 1 ( ) ε 0 ; u u H 1 ( ) c 0 λ Proof of Main Teorem Wit u, u 0 and u 0 as in te main teorem define (3.1) J u = u 0 + I (u u 0 ) u 0 + X 0, were I is te standard nodal interpolation operator. Te proof of te main teorem will use te following quantitative version of te Inverse Function Teorem wit X = u 0 +X 0, X = X 0, Y = X 0 (te dual space of X 0 ), x 0 = J u, f = D H( ; ). Te proof of te lemma follows from tat in [1] pp Lemma 3.1. Let X be an affine Banac space wit Banac space X as tangent space, and let Y be a Banac space. Suppose x 0 X and
7 PRESCRIBED MEAN CURVATURE DIRICHLET PROBLEM 7 f C 1 (X,Y). Assume tere are positive constants α, β, δ and ε suc tat were f(x 0 ) Y δ, f (x 0 ) 1 L(Y,X) α 1, f (x) f (x 0 ) L(X,Y ) β for all x B ε (x 0 ), β<α, δ (α β)ε. Ten tere exists a unique x B ε (x 0 ) suc tat f(x )=0. Remark 3.1 (Te Volume Functional). A fundamental result wic is due to Wente [12] states tat for any u 0 H 1 L (; R 3 ) te functional V (and ence D H ) extends to an analytic functional on te affine space u 0 + H0(; 1 R 3 ). Tis is peraps surprising, since from (1.5) and (1.7) one migt expect bounds for te relevant integrals to also involve ϕ L and u L respectively. More generally, one as te following. For u, v, w H 1 L (; R 3 ) define te trilinear functional V (u, v, w) =V(u, v, w;)= 1 (3.2) [u, v x,w y ]+[u, w x,v y ]. 6 Note tat V (u) =V(u, u, u). Assume now tat at least one of u, v, w also belongs to H0(; 1 R 3 ). Ten V is invariant under cyclic permutations of its arguments, as follows from integration by parts in te C 2 case and in general by a limit argument, see [10, Remark III.1.1.iii]. Since V is invariant under permutation of its first two arguments, it ten follows it is invariant under any permutation of its arguments. Moreover under te same assumptions, from an argument similar to tat in [10, proof of Teorem III.2.3] wic uses an isoperimetric inequality due to Radó [8], one also as (3.3) V (u, v, w) c u H 1 () v H 1 () w H 1 (). Similar remarks and estimates apply if is everywere replaced by. Assume now u, v H 1 L (; R 3 ) and ϕ H 1 0 L (; R 3 ). It follows tat (3.4) (3.5) (3.6) V (u),ϕ =3V(u, u, ϕ) c u 2 H () ϕ 1 H 1 (), V (u)(v, ϕ) =6V(u, v, ϕ) c u H 1 () v H 1 () ϕ H 1 (), V ( )(u, v, ϕ) =6V(u, v, ϕ) c u H 1 () v H 1 () ϕ H 1 (). (In particular, if u H 1 L (; R 3 ) tese estimates allow one to define te integrals in (1.5) and (1.7) for arbitrary ϕ, ψ H 1 0(; R 3 ).) Similar results also old if is replaced by.
8 8 GERHARD DZIUK AND JOHN E. HUTCHINSON For te remainder of tis section, u is as in te Main Teorem. Extend u to as in (2.1) and restrict to as necessary. Bot te extension and restriction will also be denoted by u. Lemma 3.2. u J u H 1 ( ) c 1 were c 1 = c 1 ( u H 2 (), u 0 H 2 (),α). Proof. u J u H 1 ( ) = (u 0 u 0 ) (I (u u 0 ) (u u 0 )) H 1 ( ) u 0 u 0 H 1 ( ) + c u u 0 H 2 () c, were c = c( u H 2, u 0 H 2,α). Since u H 1 ( ) c u H 2 ( ) c u H 2 () by elementary estimates and (2.1), te result follows. Lemma 3.3. If ϕ X 0 ten D H(J u; ),ϕ c 2 ϕ H 1 ( ), were c 2 = c 2 ( u H 2 (), u 0 H 2 (),α,h). Proof. D H(J u; ),ϕ = ( D H(J u; ),ϕ D H(u; ),ϕ ) + D H(u; ),ϕ =: A + B From te Taylor series expansion for V ( ; ) and Remark 3.1 A = (J u u) ϕ +2H V (J u; ),ϕ V (u; ),ϕ J u u H 1 ( ) ϕ H 1 ( ) +2 H V (u; )(J u u, ϕ ) + H V (u; )(J u u, J u u, ϕ ) J u u H 1 ( ) ϕ H 1 ( ) + c H u H 1 ( ) J u u H 1 ( ) ϕ H 1 ( ) + c H J u u 2 H 1 ( ) ϕ H 1 ( ) c ϕ H 1 ( ), from Lemma 3.2 and (2.1), were c = c( u H 2 (), u 0 H 2 (),α,h). Also, B = u ϕ +2H [ϕ,u x,u y ] = ( u +2Hu x u y ) ϕ = ( u +2Hu x u y ) ϕ c ϕ L 2 ( ) c ϕ H 1 ( ) were c = c( u H 2,H), as follows from (2.1), a Sobolev imbedding teorem, and elementary calculus. Te required result follows.
9 PRESCRIBED MEAN CURVATURE DIRICHLET PROBLEM 9 Remark 3.2 (A Discrete Eigenspace Decomposition). If ϕ X 0 let ϕ also denote te zero extension to. Note tat ϕ / H0() 1 unless is convex. For tis reason define P : X 0 H0() 1 to be te H 1 projection, i.e. (Pϕ ) ϕ = ϕ ϕ for all ϕ H0(). 1 One as (3.7) ϕ H 1 ( ) c 1/2 ϕ H 1 ( ) (3.8) Pϕ H 1 () ϕ H 1 ( ) (3.9) ϕ Pϕ H 1 () c 1/2 ϕ H 1 ( ) To see (3.7) note tat ϕ is constant on any triangle T T and tat T ( ) c T. Inequality (3.8) is immediate, since P is just H 1 ortogonal projection onto H0(). 1 For (3.9) first note tat ϕ Pϕ H 1 () ϕ ϕ H 1 () for any ϕ H0(), 1 by ortogonality. Now coose ϕ by suitably deforming ϕ in a boundary strip. Let (Pϕ ) +,(Pϕ ) H0() 1 be te components of Pϕ as in (1.9). Note tat (Pϕ ) is smoot, and in particular (3.10) (Pϕ ) H 2 () ν ϕ H 1 ( ) since (Pϕ ) H 2 () ν (Pϕ ) H 1 () ν Pϕ H 1 () ν ϕ H 1 ( ) from (1.11), (1.9) and te H 1-ortogonality of (Pϕ ) + and (Pϕ ), and (3.8). Define a discrete analogue of (1.9) by ϕ ( ) = I (Pϕ ) X 0, ϕ (+) = ϕ, (3.11) ϕ = ϕ (+) + ϕ ( ). Taking te zero extension of (Pϕ ) and (Pϕ ) + to, and of ϕ ( ) and ϕ (+) to, we claim (Pϕ ) ϕ ( ) H 1 ( ) c ϕ H 1 ( ), (3.12) (Pϕ ) + ϕ (+) H 1 ( ) c 1/2 ϕ H 1 ( ), were c = c(ν) Proof of claim. (Pϕ ) ϕ ( ) H 1 ( ) c (Pϕ ) H 2 () cν ϕ H 1 ( ) from (3.10). Also (Pϕ ) H 1 ( ) c (Pϕ ) H 2 () cν ϕ H 1 ( ). Tis gives te first result.
10 10 GERHARD DZIUK AND JOHN E. HUTCHINSON For te second, (Pϕ ) + ϕ (+) H 1 ( ) Pϕ ϕ H 1 ( ) + (Pϕ ) ϕ ( ) H 1 ( ) c( 1/2 + ν) ϕ H 1 ( ) from te first result and (3.9). On, ϕ = ϕ (+) required estimate now follows from (3.9). = 0 and so te We also ave (3.13) ϕ ( ) H 1 ( ) (1 + c) ϕ H 1 ( ) ϕ (+) H 1 ( ) (1 + c 1/2 ) ϕ H 1 ( ) from (3.12), te ortogonal decomposition Pϕ =(Pϕ ) +(Pϕ ) + and (3.8). Tus (3.11) is an almost ortogonal decomposition for small. Lemma 3.4. If ϕ X 0 ten D H(J u; )(ϕ,ϕ (+) ) 3λ 4 ϕ 2 H 1 ( ) provided 1 were 1 = 1 ( u H 2 (), u 0 H 2 (),α,,d,h,λ). Proof. Since V ( ; ) is cubic, from (2.3) D H(J u; )(ϕ,ϕ (+) = D H (u; )(ϕ,ϕ (+) But ) 2HV (u; )(J u u, ϕ,ϕ (+) )+2HV (u; )(J u u, ϕ,ϕ (+) ) c J u u H 1 ( ) ϕ H 1 ( ) ϕ (+) H 1 ( ) c ϕ 2 H 1 ( ) ). were c = c( u H 2 (), u 0 H 2 (),α,ν,h), from Remark 3.1, and also from Lemma 3.2 and (3.13). Now were D H (u; )(ϕ,ϕ (+) E 1 c(1 + u L ) ϕ H 1 ( ) ϕ (+) )=D H (u; )(ϕ,ϕ (+) )+E 1 H 1 ( ) c ϕ 2 H 1 ( ) wit c = c( u H 2 (),ν,h), from (3.7) and (3.12), since (Pϕ ) + = (Pϕ ) = 0 in.
11 PRESCRIBED MEAN CURVATURE DIRICHLET PROBLEM 11 Also D H (u; )(ϕ,ϕ (+) ) = D H (u; )(Pϕ,(Pϕ ) + (Pϕ ) ) +D H (u; )(ϕ Pϕ,(Pϕ ) + (Pϕ ) ) from (1.10). But +D H (u; )(ϕ, (ϕ (+) λ Pϕ 2 H 1 () + E 2 + E 3 (Pϕ ) + ) (ϕ ( ) (Pϕ ) )) E 2, E 3 c(1 + u L ) 1/2 ϕ 2 H 1 ( ) from (3.9) and (3.8), and (3.12) respectively. It follows tat D H(J u; )(ϕ,ϕ (+) ) λ Pϕ 2 H 1 () c 1/2 ϕ 2 H 1 ( ) 3λ 4 ϕ 2 H 1 ( ) from (3.9), for 1 = 1 ( u H 2 (), u 0 H 2 (),α,,d,h,λ). Lemma 3.5. If v u 0 + X 0 and ϕ,ψ X 0 ten D H(v ; )(ϕ,ψ ) D H(J u; )(ϕ,ψ ) λ 4 ϕ H 1 ( ) ψ H 1 ( ) provided v J u H 1 ( ) ε 1 were ε 1 = ε 1 (H, λ). Proof. Tis follows from D H(v ; )(ϕ,ψ ) D H(J u; )(ϕ,ψ ) =2HV (J u; )(v J u, ϕ,ψ ), (3.6) and Lemma 3.2. Completion of proof of Main Teorem. We use Lemma 3.1 wit X = u 0 + X 0, X = X 0, Y = X 0, x 0 = J u, f = D H( ; ). Te norm on X 0 is H 1 ( ) and on X 0 is te corresponding dual norm. Note tat D H( ; ): u 0 +X 0 X0 wit derivative D H( ; ): u 0 +X 0 L(X 0,X0) using standard identifications. From Lemma 3.3 (3.14) D H(J u; ) c 2. From Lemma 3.4, D H(J u; ) is invertible and ( [D H(J u; )] 1 / 3λ 4 ϕ H 1 ( ) ϕ (+) H 1 ( ) ) 1
12 12 GERHARD DZIUK AND JOHN E. HUTCHINSON provided 1. But ϕ (+) H 1 ( ) (Pϕ ) + (Pϕ ) H 1 ( ) + ϕ (+) (Pϕ ) + H 1 ( ) + ϕ ( ) (Pϕ ) H 1 ( ) (1 + c 1/2 ) ϕ H 1 ( ) were c = c(ν), from (3.8) and (3.12). Hence (3.15) ( ) 1 [D H(J u; )] 1 λ 2 if 3 were 3 = 3 ( u H 2 (), u 0 H 2 (),α,,d,h,λ). Finally, from Lemma 3.5 (3.16) D H(v ; ) D H(J u; ) λ 4 if v J u H 1 ( ) ε 1 were ε 1 = ε 1 (H, λ). Take δ = c 2, α = λ/2, β = λ/4 and ε = ε 1. Ten from (3.14) (3.16) te ypoteses of Lemma 3.1 are satisfied provided 3, c 2 λ ε 4 1. Tis establises te first (uniqueness) part of te main teorem wit ε 0 = ε 1 and 0 = 0 ( 3,ε 1,λ,c 2 )= 0 ( u H 2, u 0 H 2,α,ν,H,λ). Taking δ = c 2, α = λ/2, β = λ/4 and ε = λ 1 c 0 te ypoteses of Lemma 3.1 are again satisfied from (3.14) (3.16) provided 3, c 2 1c 4 0. Tis establises te second (O() convergence) part of te main teorem wit 0 = 3 and c 0 =4c Numerical Results In Tables 1 and 2 we present te results of test computations for te explicitly known sperical solutions described in Example 1.1 wit H =0.5 and = B 1 (0). Denote by e te error between te continuous solution and te discrete solution in te cosen norm. For two successive grids wit grid sizes 1 and 2 te experimental order of convergence is eoc =ln e / 1 ln 1. e 2 2 Te test computations confirm te order 1 for te H 1 ()-norm and additionally sow te order 2 for te L 2 ()-norm. Figures 2 and 3 sow computational results wit = B 1 (0), H = 0.5 and boundary values u(e iφ ) = (cos(φ), sin(φ), (2 + 3) cos(2φ) 0.5 cos(6φ)) on a grid wit 8192 triangles. For better visibility te resulting surfaces are scaled, but te boundaries of te solution surfaces are te same. Figure 4 sows a solution for te annular domain = { x 1 < x < 2 } and boundary data wic give knotted boundary curves.
13 PRESCRIBED MEAN CURVATURE DIRICHLET PROBLEM 13 nodes level L 2 -error L 2 -eoc H 1 -error H 1 -eoc e e e e e e e e e e e e Table 1. Small solution, H =0.5 nodes level L 2 -error L 2 -eoc H 1 -error H 1 -eoc Table 2. Large solution, H =0.5 Figure 2. Small solution, H =0.5 References [1] M. Berger, Nonlinearity and Functional Analysis, Academic Press, [2] H. Brezis, J.-M. Coron, Multiple solutions of H systems and Rellic s conjecture, Comm. Pure Appl. Mat. 37, 1984,
14 14 GERHARD DZIUK AND JOHN E. HUTCHINSON Figure 3. Large solution, H =0.5 Figure 4. Annulus, H =0.5 [3] P. G. Ciarlet, Te Finite Element Metods for Elliptic Problems, Nort Holland, [4] M. Dobrowolski, R. Rannacer Finite element metods for nonlinear elliptic systems of second order, Mat. Nacr. 94, 1974, [5] G. Dziuk, J.E. Hutcinson, On te approximation of unstable parametric minimal surfaces, Calc. Var. 4, 1996, [6] G. Dziuk, J.E. Hutcinson, A finite element metod for approximating minimal surfaces, Preprint 4 Matematisce Fakultät Freiburg, 1996 CMA Mat. Res. Rep. 5, Australian National University, [7] G. Dziuk, J.E. Hutcinson, A finite element metod for approximating surfaces of prescribed mean curvature, in preparation.
15 PRESCRIBED MEAN CURVATURE DIRICHLET PROBLEM 15 [8] T. Radó, Te isoperimetric inequality and te Lebesgue definition of surface area, Trans. Amer. Mat. Soc. 61, 1947, [9] M. Rumpf, A. Scmidt et al, GRAPE, Grapics Programming Environment, Report 8, SFB 256, Bonn (1990). [10] M. Struwe, Plateau s Problem and te Calculus of Variations, Princeton University Press, [11] M. Struwe, Variational Metods, Springer Verlag, [12] H. C. Wente, An existence teorem for surfaces of constant mean curvature, J. Mat. Anal. Appl. 26, 1969, Institut für Angewandte Matematik, Universität Freiburg, Hermann Herder Str. 10, D Freiburg i. Br., GERMANY address: gerd@matematik.uni-freiburg.de Scool of Matematical Sciences, Australian National University, GPO Box 4, Canberra, ACT 0200, AUSTRALIA address: Jon.Hutcinson@anu.edu.au
THE 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 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 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 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 informationIsoparametric finite element approximation of curvature on hypersurfaces
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
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 informationA Finite Element Primer
A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.
More informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
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 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 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 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 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 informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
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 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 informationMath 161 (33) - Final exam
Name: Id #: Mat 161 (33) - Final exam Fall Quarter 2015 Wednesday December 9, 2015-10:30am to 12:30am Instructions: Prob. Points Score possible 1 25 2 25 3 25 4 25 TOTAL 75 (BEST 3) Read eac problem carefully.
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 informationA posteriori error estimates for non-linear parabolic equations
A posteriori error estimates for non-linear parabolic equations R Verfürt Fakultät für Matematik, Rur-Universität Bocum, D-4478 Bocum, Germany E-mail address: rv@num1rur-uni-bocumde Date: December 24 Summary:
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 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 informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
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 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 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 informationCHAPTER 4. Elliptic PDEs
CHAPTER 4 Elliptic PDEs One of te main advantages of extending te class of solutions of a PDE from classical solutions wit continuous derivatives to weak solutions wit weak derivatives is tat it is easier
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 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 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 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 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 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 informationTest 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 =
Test Review Find te determinant of te matrix below using (a cofactor expansion and (b row reduction Answer: (a det + = (b Observe R R R R R R R R R Ten det B = (((det Hence det Use Cramer s rule to solve:
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 informationConvexity and Smoothness
Capter 4 Convexity and Smootness 4.1 Strict Convexity, Smootness, and Gateaux Differentiablity Definition 4.1.1. Let X be a Banac space wit a norm denoted by. A map f : X \{0} X \{0}, f f x is called a
More informationCharacterization of reducible hexagons and fast decomposition of elementary benzenoid graphs
Discrete Applied Matematics 156 (2008) 1711 1724 www.elsevier.com/locate/dam Caracterization of reducible exagons and fast decomposition of elementary benzenoid graps Andrej Taranenko, Aleksander Vesel
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 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 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 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 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 information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationOn convergence of the immersed boundary method for elliptic interface problems
On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be 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 informationGeneric maximum nullity of a graph
Generic maximum nullity of a grap Leslie Hogben Bryan Sader Marc 5, 2008 Abstract For a grap G of order n, te maximum nullity of G is defined to be te largest possible nullity over all real symmetric n
More informationAPPROXIMATION 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 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 informationTHE IMPLICIT FUNCTION THEOREM
THE IMPLICIT FUNCTION THEOREM ALEXANDRU ALEMAN 1. Motivation and statement We want to understand a general situation wic occurs in almost any area wic uses matematics. Suppose we are given number of equations
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 informationSemigroups of Operators
Lecture 11 Semigroups of Operators In tis Lecture we gater a few notions on one-parameter semigroups of linear operators, confining to te essential tools tat are needed in te sequel. As usual, X is a real
More informationarxiv:math/ v1 [math.ca] 1 Oct 2003
arxiv:mat/0310017v1 [mat.ca] 1 Oct 2003 Cange of Variable for Multi-dimensional Integral 4 Marc 2003 Isidore Fleiscer Abstract Te cange of variable teorem is proved under te sole ypotesis of differentiability
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
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 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 informationSmoothness of solutions with respect to multi-strip integral boundary conditions for nth order ordinary differential equations
396 Nonlinear Analysis: Modelling and Control, 2014, Vol. 19, No. 3, 396 412 ttp://dx.doi.org/10.15388/na.2014.3.6 Smootness of solutions wit respect to multi-strip integral boundary conditions for nt
More informationFigure 1. Figure. Anticipating te use of two-level domain decomposition preconditioners, we construct a triangulation of in te following way. Let be d
Lower Bounds for Two-Level Additive Scwarz Preconditioners wit Small Overlap * Susanne C. Brenner Department of Matematics University of Sout Carolina Columbia, SC 908 Summary. Lower bounds for te condition
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 information1. Introduction. Consider a semilinear parabolic equation in the form
A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
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 informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
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 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 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 informationControllability of a one-dimensional fractional heat equation: theoretical and numerical aspects
Controllability of a one-dimensional fractional eat equation: teoretical and numerical aspects Umberto Biccari, Víctor Hernández-Santamaría To cite tis version: Umberto Biccari, Víctor Hernández-Santamaría.
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 information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationITERATIVE SOLUTION OF A CONSTRAINED TOTAL VARIATION REGULARIZED MODEL PROBLEM
ITERATIVE SOLUTION OF A CONSTRAINED TOTAL VARIATION REGULARIZED MODEL PROBLEM Abstract. Te discretization of a bilaterally constrained total variation minimization problem wit conforming low order finite
More informationPROJECTION-FREE APPROXIMATION OF GEOMETRICALLY CONSTRAINED PARTIAL DIFFERENTIAL EQUATIONS
PROJECTION-FREE APPROXIMATION OF GEOMETRICALLY CONSTRAINED PARTIAL DIFFERENTIAL EQUATIONS Abstract. We devise algoritms for te numerical approximation of partial differential equations involving a nonlinear,
More informationTopics in Generalized Differentiation
Topics in Generalized Differentiation J. Marsall As Abstract Te course will be built around tree topics: ) Prove te almost everywere equivalence of te L p n-t symmetric quantum derivative and te L p Peano
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 informationarxiv: v1 [math.na] 29 Jan 2019
Iterated total variation regularization wit finite element metods for reconstruction te source term in elliptic systems Micael Hinze a Tran Nan Tam Quyen b a Department of Matematics, University of Hamburg,
More informationFINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS
J. KSIAM Vol.22, No.2, 101 113, 2018 ttp://dx.doi.org/10.12941/jksiam.2018.22.101 FINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS DEOK-KYU JANG AND JAE-HONG PYO DEPARTMENT
More informationREGULARITY OF SOLUTIONS OF PARTIAL NEUTRAL FUNCTIONAL DIFFERENTIAL EQUATIONS WITH UNBOUNDED DELAY
Proyecciones Vol. 21, N o 1, pp. 65-95, May 22. Universidad Católica del Norte Antofagasta - Cile REGULARITY OF SOLUTIONS OF PARTIAL NEUTRAL FUNCTIONAL DIFFERENTIAL EQUATIONS WITH UNBOUNDED DELAY EDUARDO
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
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 informationarxiv: v3 [math.na] 15 Dec 2009
THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING M.A. EBRAHIMI, MICHAEL HOLST, AND EVELYN LUNASIN arxiv:91.4548v3 [mat.na] 15 Dec 9 ABSTRACT. In tis paper we investigate te use of te D Navier-Stokes-Voigt
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 informationBILAYER PLATES: MODEL REDUCTION, Γ-CONVERGENT FINITE ELEMENT APPROXIMATION AND DISCRETE GRADIENT FLOW
BILAYER PLATES: MODEL REDUCTION, Γ-CONVERGENT FINITE ELEMENT APPROXIMATION AND DISCRETE GRADIENT FLOW SÖREN BARTELS, ANDREA BONITO, AND RICARDO H. NOCHETTO Abstract. Te bending of bilayer plates is a mecanism
More informationNUMERICAL APPROXIMATION OF THE LANDAU-LIFSHITZ-GILBERT EQUATION AND FINITE TIME BLOW-UP OF WEAK SOLUTIONS
NUMERICAL APPROXIMATION OF THE LANDAU-LIFSHITZ-GILBERT EQUATION AND FINITE TIME BLOW-UP OF WEAK SOLUTIONS SÖREN BARTELS, JOY KO, AND ANDREAS PROHL Abstract. Te Landau-Lifsitz-Gilbert equation describes
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 informationf a h f a h h lim lim
Te Derivative Te derivative of a function f at a (denoted f a) is f a if tis it exists. An alternative way of defining f a is f a x a fa fa fx fa x a Note tat te tangent line to te grap of f at te point
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationFINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES
FINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES MICHAEL HOLST AND CHRIS TIEE ABSTRACT. Over te last ten years, te Finite Element Exterior Calculus (FEEC) as
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 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 informationMath Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim
Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z
More informationOSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix
Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra
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 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 informationInfluence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations
International Journal of Difference Equations ISSN 0973-6069, Volume 12, Number 2, pp. 281 302 (2017) ttp://campus.mst.edu/ijde Influence of te Stepsize on Hyers Ulam Stability of First-Order Homogeneous
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 informationSubdifferentials of convex functions
Subdifferentials of convex functions Jordan Bell jordan.bell@gmail.com Department of Matematics, University of Toronto April 21, 2014 Wenever we speak about a vector space in tis note we mean a vector
More informationMultigrid Methods for Obstacle Problems
Multigrid Metods for Obstacle Problems by Cunxiao Wu A Researc Paper presented to te University of Waterloo in partial fulfillment of te requirement for te degree of Master of Matematics in Computational
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 informationCELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION SCHEME WITHOUT FICTITIOUS DOMAINS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7, Number 1, Pages 1 9 c 010 Institute for Scientific Computing and Information CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More information