Posteriori Analysis of a Finite Element Discretization for a Penalized Naghdi Shell
|
|
- Lesley Rice
- 6 years ago
- Views:
Transcription
1 International Journal of Difference Equations ISSN , Volume 8, Number 1, pp (2013) ttp://campus.mst.edu/ijde Posteriori Analysis of a Finite Element Discretization for a Penalized Nagdi Sell Nora Tabouce Laboratoire de Matématiques Appliquées et Modélisation Université 8 Mai 1945 de Guelma, Guelma, Algeria tabmoufida@yaoo.fr Abstract We consider a penalized Nagdi model in Cartesian coordinates for linearly elastic sells wit little regularity. A posteriori analysis of te discrete problem leads to te construction of error indicators, wic satisfy optimal estimates. We describe a mes adaptivity strategy relying on tese indicators and we present a numerical experiment tat confirms its efficiency. AMS Subject Classifications: 74K25, 74S05. Keywords: Nagdi s sell model, finite element approximation, a posteriori error estimates, residuals. 1 Introduction Nagdi s model is a linear elastic sell. Te formulation of te model used ere was introduced in [4, 6]. A posteriori analysis is now an important tool for improving te efficiency of te discretization. We refer to [7, 8] for te first works concerning a plate model and to [2] for sell models. Te first aim of a posteriori analysis is mes adaptivity. Indeed, a muc smaller number of degrees of freedom are needed to obtain a given accuracy wen te final mes is adapted to te solution and te construction of suc mes relies on error indicators, wic only depend on te discrete solution, and ence can be computed in an explicit way and often in a non expensive way. A posteriori estimate proves tat tese indicators provide a good representation of te local error see [11] for a detailed presentation of all tis. Here we perform a posteriori analysis of te discretization, relying on a penalized version studied in [5] and prove upper and lower bounds for te error as a function of residual type indicators. Finally, we describe Received June 29, 2012; Accepted January 27, 2013 Communicated by Sandra Pinelas
2 112 N. Tabouce te strategy tat is used for te adaptivity mes. Numerical experiments are in good agreement wit te analysis. Te article is organized as follows. We first briefly recall te geometry of te midsurface and Nagdi sell formulation given in [4, 6]. Tis formulation involves te infinitesimal rotation vector, a vector unknown tat is tangent to te midsurface. In Section 3, we recall te penalized version of Nagdi s model intended to approximate te above mentioned tangency. Section 4 is devoted to te a posteriori analysis of te finite element discretization. In Section 5 we present te adaptivity strategy and numerical experiments. 2 Presentation of te Model Greek indices and exponents take teir values in te set {1, 2}, wile Latin indices and exponents belong to te set {1, 2, 3}. Let (e 1, e 2, e 3 ) be te canonical ortonormal basis of R 3. We denote by u v te inner product of R 3, u = u u te associated Euclidean norm and u v te vector product of u and v. Let be a bounded connected domain of R 2. We consider a sell of midsurface S = ϕ(), were ϕ W 2, (, R 3 ) is a one-to-one mapping suc tat te two vectors a α (x) = α ϕ(x) are linearly independent at eac point x of. We let a 3 = a 1 a 2 be te unit normal vector on te midsurface a 1 a 2 at point ϕ(x). Te vectors a i (x) define te local covariant basis at point ϕ(x). Te contravariant basis a i (x) is defined by a i a j = δ j i, were δj i is de Kronecker symbol. We let a(x) = a 1 (x) a 2 (x) 2, so tat a(x) is te area element of te midsurface in te cart ϕ. Te first and te second fundamental forms of te surface are given in covariant components by a αβ (x) = a α (x) a β (x) and b αβ (x) = a 3 (x) β a α (x). Te contravariant components of te first fundamental form a αβ (x) = a α (x) a β (x). Te lengt element l on te boundary is given by a αβ τ α τ β, (τ 1, τ 2 ) being te covariant coordinates of a unit vector tangent to. Let a αβρσ L () be te elasticity tensor. We consider ere te case of a omogeneous, isotropic material wit Young modulus E > 0 and Poisson ratio ν, 0 ν 1, were tese components are given by 2 a αβρσ = E 2(1 + ν) (aαρ a βσ + a ασ a βρ ) + Eν 1 ν 2 aαβ a ρσ. (2.1) In tis context, te covariant components of te cange of metric tensor read γ αβ (u) = 1 2 ( αu a β + β u a α ), (2.2) te covariant components of te cange of curvature tensor read χ αβ (u, r) = 1 2 ( αu β a 3 + β u α a 3 + α r a β + β r a α ), (2.3)
3 Posteriori Analysis of a Finite Element Discretization 113 and te components of te cange of transverse sear tensor read δ α3 (u, r) = 1 2 ( αu a 3 + r a α ). (2.4) We assume tat te boundary of te cart domain is divided into two parts: γ 0 on wic te sell is clamped and γ 1 = \γ 0 on wic te sell is subjected to applied traction and moment. Let us now consider te function space V() introduced in [4,6]: V() = {V = (v, s) [ H 1 (, R 3 ) ] } 2, s a3 = 0 in, v = s = 0 on γ 0. (2.5) Tis space is endowed wit te natural Hilbert norm V V = ( v 2 H 1 (;R 3 ) + s 2 H 1 (;R 3 ) ) 1 2. (2.6) We now recall te variational formulation of te problem corresponding to te linear Nagdi model wit data (f, N, M) L 2 (; R 3 ) L 2 (γ 1 ; R 3 ) L 2 (γ 1 ; R 3 ): find U = (u, r) V() suc tat V V, a(u, V ) = L(V ), (2.7) were te bilinear form a(, ) is defined by ] a(u, V ) = {ea [γ αβρσ αβ (u)γ ρσ (v) + e2 12 χ αβ(u)χ ρσ (V ) and te linear form L( ) is given by L(V ) = f v a dx + +2e E } a 1 + ν aαβ δ α3 (U)δ β3 (V ) dx, (2.8) γ 1 (N v + M s) l dτ. (2.9) Te data f, N, M represent a given resultant force density, an applied traction density and an applied moment density, respectively. In te above formulas, te tickness of te sell, denoted by e, is assumed to be constant and positive. We refer to [4,6] for te proof of te following results: Te form L is continuous on V() and its norm satisfies ( ) L c f L 2 (;R 3 ) + N L 2 (γ 1 ;R 3 ) + M L 2 (γ 1 ;R 3 ). (2.10) Tere exists a constant c > 0 suc tat V V(), a(v, V ) c V 2 V(). (2.11) Problem (2.7) admits a unique solution U in V().
4 114 N. Tabouce 3 Penalized Version We consider te penalized Nagdi problem introduced in [5], in wic te unknowns are te displacement u and te rotation r, elements of te space H 1 (; R 3 ) witout any ortogonality constraint on r. Te relaxed function space is defined by X = { V = (v, s) H 1 (, R 3 ) 2, v = s = 0 on γ 0 }, (3.1) equipped wit te norm defined in (2.6), wic is now denoted by X(). Teorem 3.1. Let p R be suc tat 0 < p 1. Let f L 2 (, R 3 ), N L 2 (γ 1, R 3 ) and M L 2 (γ 1, R 3 ). Ten tere exists a unique solution to te following problem: find U p = (u p, r p ) X suc tat were Proof. See [5]. V X, a(u p, V ) + 1 p b(r p a 3 ; s a 3 ) = L(V ), (3.2) b(λ; µ) = α λ α µ dx. (3.3) Remark 3.2. We note tat te quantity a(u, V ) can be written in anoter form wic seems more appropriate for implementation, since it uncouples te two components v and s of te test function V = (v, s). Indeed, we introduce te contravariant components of te stress resultant, n ρσ (u) = ea αβρσ γ αβ (u), (3.4) of te stress couple, m ρσ (U) = e3 12 aαβρσ χ αβ (U), (3.5) and of te transverse sear force, We also ave: t β (U) = e E 1 + ν aαβ δ α3 (U). (3.6) χ ρσ (V ) = θ ρσ (v) + γ ρσ (s) wit θ ρσ (v) = 1 2 ( ρv σ a 3 + σ v ρ a 3 ). (3.7) Tus, te bilinear form a (U p, V ) can be rewritten as ( ) a (U p, V ) = n ρσ (u p )γ ρσ (v) + m ρσ (U p )θ ρσ (v) + t β (U p ) β v a 3 a dx ( ) + m ρσ (U p )γ ρσ (s) + t β (U p )s a β a dx. (3.8)
5 Posteriori Analysis of a Finite Element Discretization 115 To go furter, using tis new form togeter wit te symmetry properties n ρσ (u p ) = n σρ (u p ) and m ρσ (U p ) = m σρ (U p ), and by integration by parts in problem (3.2), we obtain te strong formulation of te penalized problem: ρ ((n ρσ (u p )a σ + m ρσ (U p ) σ a 3 + t ρ (U p )a 3 ) a) = f a in, ρ (m ρσ (U p )a σ a) + t β 1 (U p )a β a p ρρ(r p a 3 )a 3 = 0 in, u p = r p = 0 on γ 0, ν ρ (n ρσ (u p )a σ + m ρσ (U p ) σ a 3 + t ρ (U p )a 3 ) (3.9) a = Nl on γ 1, ν ρ (m ρσ 1 (U p )a σ a + p ρ(r p a 3 )a 3 ) = Ml on γ 1. Te discretization tat we intend to study is constructed by te Galerkin metod from problem (3.2). We refer to [5, 5.1] for more details. 4 A Posteriori Analysis of te Discrete Problem Let (T ) be a regular affine family of triangulations wic covers te domain and P k (K) denote te space of restrictions to K, element of T, of polynomials wit total degree k. Te discrete space is given by X = { V = (v, s ) C 0 (; R 3 ) 2, V K P 1 (K), v = s = 0 on γ 0 }. (4.1) Te discrete problem consists to find U p, = (u p,, r p, ) X suc tat V X, a(u p,, V ) + 1 p b(r p, a 3, s a 3 ) = L(V ). (4.2) Tis problem as a unique solution [5, T. 3.1]. Te a posteriori analysis of problem (4.2) relies on te residual equation a(u p U p,, V ) + 1 p b((r p a 3 r p, a 3 ); s a 3 ) = L(V V ) a(u p,, V V ) 1 p b (r p, a 3 ; (s a 3 s a 3 )), (4.3) valid for all V X() and for all V X. Te construction of error indicators from tese equations requires approximations of te data and of te coefficients [2, 11]. 4.1 Approximation of te Data Let E 1 denote te set of edges of elements of T wic are contained in γ 1. We consider an approximation f of f in Z and approximations N and M of N and M in Z 1, were te spaces Z and Z 1 are defined by Z = { g L 2 () 3 ; K T, g K P 0 (K) 3}, Z 1 = { E L 2 (γ 1 ) 3 ; e E 1, E e P 0 (e) 3}. (4.4)
6 116 N. Tabouce 4.2 Approximation of te Coefficients We denote by a αβ, aαβρσ, ( a) and l, te approximations of a αβ, a αβρσ, a and l, respectively, in te space M wic is given by M = { χ H 1 (); K T, χ K P 1 (K) }. (4.5) Similarly, we consider approximations a k of te vectors a k and d α of te α a 3 in te space M. We also agree to denote by γ αβ( ), χ αβ( ) and δ α3( ) te approximations of te tensors introduced in (2.2) to (2.4). For instance, γ αβ( ) is given by γ αβ(u) = 1 2 ( ) α u a β + β u a α. (4.6) Tis leads to te definition of te approximate linear form L (V ) = f v( a) dx + (N v + M s)l dτ, (4.7) γ 1 and also of approximate bilinear forms { ] a (U, V ) = ea αβρσ [γ αβ (u) γ ρσ (v) + e2 12 χ αβ (U) χ ρσ(v ) +2e E } 1 + ν aαβ δ α3 (U) δβ3 (V ) ( (4.8) a) dx, ( ) b (r a 3, s a 3 ) = α r a 3 + r d α α s a 3 dx. (4.9) It is easy to ceck tat V X() and V X one as L(V V ) a (U p,, V V ) 1 p b (r p, a 3 ; (s a 3 s a 3 )) = (L L )(V V ) + L (V V ) (a a )(U p,, V V ) a (U p,, V V ) 1 p (b b )(r p, a 3 ; (s a 3 s a 3 )) (4.10) 1 p b (r p, a 3 ; (s a 3 s a 3 )). To go furter, we recall some standard notations: (i) E K denotes te set of edges of K wic are not contained in γ 0 and E 1 K te set of elements of E K wic are contained in γ 1 ; (ii) for eac e E K, ν = (ν 1, ν 2 ) is a unit vector normal to e, wit te furter assumption tat, wen e belongs to E 1 K, ν is outward to ;
7 Posteriori Analysis of a Finite Element Discretization 117 (iii) for eac e E K, e stands for te lengt of e; (iv) for eac e E K \E 1 K, [ ] e denotes te jump troug e; (v) K is te union of triangles of T tat sare an edge wit K; (vi) K is te union of triangles of T tat intersect K. We recall tat from [3, Teorem IX.3.11 and Corollary IX.3.12] tere exists a Clément type operator R wic maps Hγ 1 0 () into M γ 0 = M Hγ 1 0 and satisfies, for all functions χ Hγ 1 0 (), eac K T and eac e of K wic is not contained in γ 0, χ R χ L 2 (K) + K χ R χ H 1 (K) c K χ H 1 ( K ), χ R χ L 2 (e) c 1 2 χ H 1 ( K ). (4.11) Te idea is to take V equal to (R v, R s) and χ equal to R χ in (4.3). From [2, Lemma 3.3, Lemma 3.4 and Lemma 3.5], we define te quantities linked to te local approximation error on te data: for eac K T, = K f f L 2 (K) 3 + ( ) 1 2 e N N L 2 (e) 3 + M M L 2 (e) 3, (4.12) ε (d) K e E 1 K and also from te global approximation error on te coefficients ε (c) ( = a ( a) L () l l L (γ 1 ) + sup 1 α,β,ρ,σ 2 + sup 1 k 3 a αβρσ a αβρσ ak a k L () L () 3 + sup 1 α 2 + sup 1 α,β 2 a αβ a αβ L () ) α a 3 d α L L. () 3 (4.13) We are now in a position to prove te a posteriori error estimate. In order to state it, we define te error indicators. We use Remark 3.2 to write a(u, V ) and observe tat a similar form olds for a (U, V ), wit te obvious notation for te quantities n ρσ ( ), mρσ ( ), t β ( ) and θ ρσ( ) (in comparison wit (3.4) to (3.7), all coefficients are replaced by teir approximations). For eac K T, te error indicator η K is defined by wit η K = η K1 + η K2 (4.14) η K1 = f K ( ( a) + ρ (n ρσ (u p,)a σ + m ρσ (U p,)d σ + t ρ (U p,) a 3 )( a) ) L 2 (K) [ ( 2 e νρ n ρσ (u p,)a σ + m ρσ (U p,)d σ + t ρ (U ) ] p,)a 3 ( a) e E K \E 1 K + e L 2 (e) e N l ν ρ ( n ρσ (u p,)a σ + m ρσ (U p,)d σ + t ρ (U p,)a 3) ( a) L 2 (e) 3 e E 1 K (4.15)
8 118 N. Tabouce and η K2 = ρ ( K m ρσ (U p,)a σ( ) a) t β (U p,)a β( a) + 1 ( ) ρ ( ρ (r p, a 3 )a 3 ρ (r p, a 3 )d p ρ) L 2 (K) 3 + [ 1 2 e ν ρ m ρσ (U p,)a σ( a) + 1 ] p ν ρ ρ (r p, a 3 )a 3 e E K \EK 1 e L 2 (e) e M l ν ρ m ρσ (U p,)a σ( a) 1 p ν ρ ρ (r p, a 3 )a 3. L 2 (e) 3 e E 1 K (4.16) Note tat tese indicators are of residual type and easy to compute since tey only involve polynomial functions. 4.3 Te Main Results Teorem 4.1. For any data (f, N, M) in L 2 (; R 3 ) L 2 (γ 1 ; R 3 ) L 2 (γ 1 ; R 3 ), te following a posteriori error estimate between te solution U p of problem (3.2) and te solution U p, of problem (4.2) olds: ( ( ) ) 1 2 U p U p, X() c ηk 2 + ε (d)2 K + ε (c). (4.17) K T Proof. We give an abridged proof. From te ellipticity property (2.11), using te residual equation (4.3) by replacing its second member by (4.10), we ten use te triangle inequality tat, combined wit [2, Lemma 3.3, Lemma 3.4 and Lemma 3.5], leads to te bound of te following quantity: L (V V ) a (U p,, V V ) 1 p b (r p, a 3 ; (s a 3 s a 3 )) wit V = (R v, R s). We can write L (V V ) = L ((v R v), 0) + L (0, (s R s)), and a (U p,, V V ) = a (U p,, (v R v, 0)) + a (U p,, (0, s R s), 1 p b (U p,, V V ) = 1 p b (r p, a 3 ; (v R v, 0)) + 1 p b (r p, a 3 ; (0, s a 3 R s a 3 )).
9 Posteriori Analysis of a Finite Element Discretization 119 Note tat b (r p, a 3 ; (v R v, 0)) = 0. It remains to bound te following two quantities: and A 2 = A 1 = sup s H 1 γ 0 (;R 3 ) L (v R v, 0) a (U p,, (v R v, 0)) sup v Hγ 1 0 (;R 3 ) v H 1 (;R 3 ) 1 s H 1 (;R 3 ) { L (0, s R s) a (U p,, (0, s R s)) 1 p b (r p, a 3 ; (0, (s R s) a 3 )) Setting w = v R v and using te symmetry properties of te n ρσ ( ) and mρσ ( ), we ave L (v R v, 0) a (U p,, (v R v, 0)) = f w( a) dx + N wl dτ n ρσ (u p,) ρ w a σ + m ρσ (U p,) ρ w d σ + t β (U p,) β w a 3)( a) dx. By cutting te integrals on into te sum of integrals on te K in T and integrating by parts on eac K, we derive L (v R v, 0) a (U p,, (v R v, 0)) = f w( a) dx + N w l dτ γ 1 }. γ 1 + K T ( K K ρ ((n ρσ (u p,)a σ + m ρσ (U p,)d σ + t β (U p,)a 3)( a) ) w dx ν ρ (n ρσ (u p,)a σ + m ρσ (U p,)d σ + t β (U p,)a 3)( a) w dτ Using te Caucy Scwarz inequality combined wit (4.11) leads to ). (4.18) A 1 c ( K T η 2 K1 Setting t = s R s, and using te same arguments as previously, we arrive to ) 1 2. Tis concludes te proof. A 2 c ( K T η 2 K2 ) 1 2.
10 120 N. Tabouce Teorem 4.2. For any data (f, N, M) L 2 (; R 3 ) L 2 (γ 1 ; R 3 ) L 2 (γ 1 ; R 3 ), te following bounds old for all indicators defined in (4.15) (4.16): ( ) 1 ( ) 2 2 η Ki c U p U p, X(K ) + ε (d) K + ε (c), i = 1, 2. (4.19) K K Proof. Tis is an abridged proof of te estimate for η K1. We write η K1 in compact form as η K1 = K F L 2 (K) e [G ] e L 2 (e) e N l G L. 2 (e) 3 e E K \E 1 K e E 1 K We take R v = 0 in (4.18) and w = v equal to { F ψ w = K on K, 0 on \K, were ψ K denotes te bubble function on K. Tus, all te terms on te rigt-and side of (4.18) vanis but te integral on K. We observe tat function w (tus F ) on K is a polynomial of degree 3. From te appropriate inverse inequality [3, Proposition VII.4.1] togeter wit (4.10) and (4.3) combined wit [2, Lemma 3.3, Lemma 3.4 and Lemma 3.5], leads to K F L 2 (K) 3 c ( U p U p, X(K) + ε (d) K + ) ε(c) K. (4.20) Similarly, to bound η K1 for any edge e sared by two elements K and K, we take w in (4.18) equal to { P e,k ([G ] e ψ e ) on k K, K }, w = { 0 on \ K K }, were ψ e is te bubble function on e and P e,k is te lifting operator introduced in [3, Lemma XI 2.7] from polynomials on e vanising on e into polynomials on K vanising on K\e, constructed by an affine transformation from a fixed lifting operator on te reference triangle. Finally, for eac e in EK, 1 we take te function w in (4.18) equal to { Pe,k ((N w = l G )ψ e ) on K, 0 on \K, and tis gives te bound for η K1.
11 Posteriori Analysis of a Finite Element Discretization Te Adaptivity Strategy and Numerical Experiment Now we present te adaptivity strategy and some numerical experiments. 5.1 Adaptivity Strategy (i) Construct a first mes T 0. Set i = 0. (ii) Solve te numerical problem on T i. Let u i denote te solution. (iii) Compute η Ki (u i ) on K i T i. (iv) If te global estimator is small enoug, ten stop. (v) Oterwise, compute te new step of mes Ki+1 : Ki+1 = 1 2 K i if η Ki (u i ) T OL, oterwise Ki+1 = l Ki, l 1. Here T OL = 1 n ti K i T i were n ti is te number of triangles of T i. (vi) Generate te new mes T i Numerical Experiment and return to (ii). η Ki (u i ), Te numerical experiment tat we present as been performed on te finite element code FreeFem++, see [10]. Tree-dimensional visualization of te deformed sells uses Medit, a free mes visualization software available at ttp:// jussieu.fr/ frey/logiciels/medit.tml. Hyperbolic paraboloid sell Te reference domain is te square { = (x, y); x + y } 2b, (5.1) as illustrated in [1, and 2.4.2] and te cart ϕ is defined by ( c ) T ϕ(x, y) = x, y, 2b 2 (x2 y 2 ). (5.2) We coose b = 50 cm, c = 10 cm and e = 0.8 cm. We assume tat te sell is clamped on γ 0 = and tat it is subjected to uniform pressure q = 0.01 kp/cm 2. Te mecanical data are E = kp/cm 2, ν = 0.4. Te reference value for tis test is te
12 122 N. Tabouce normal displacement at te center C(0, 0) of te sell. Its value computed by various metods is of cm; see [1]. We take p = 10 3 E, see [5]. 2(1 + ν) Results wit mes adaptation Iteration Degrees of freedom U 3 (C) Results witout mes adaptation Degrees of freedom U 3 (C) Note tat te resolution of te discrete problem wit mes adaptation gives te convergence to te solution after te tird iteration wit a number of degrees of freedom equal to but te resolution of te problem witout mes adaptation gives an approximate solution wit an error of order and a number of degrees of freedom was lower tan calculated by mes adaptation. Figure 5.1 presents te initial and te final adapted meses according to te strategy described above. Figure 5.2 and 5.3 present te over-deformed sell. Figure 5.1: Te initial and adapted meses
13 Posteriori Analysis of a Finite Element Discretization 123 Figure 5.2: Te over-deformed top side Figure 5.3: Te over-deformed bottom side References [1] M. Bernadou, Métodes d éléments finis pour les problèmes de coques minces, Recerces en Matématiques Appliquées 33, Masson, [2] C. Bernardi, A. Blouza, F. Hect and H. Le Dret, A posteriori analysis of finite element discretization of a Nagdi sell model, Laboratoire J-L-Lions de Paris VI, preprint R08048, [3] C. Bernardi, Y. Maday and F. Rapetti, Discrétisations variationnelles de problèmes aux limites elliptiques, Matématiques & Applications (Berlin), 45, Springer, Berlin, [4] A. Blouza, Existence et unicité pour le modèle de Nagdi pour une coque peu régulière, C. R. Acad. Sci. Paris Sér. I Mat. 324 (1997), no. 7,
14 124 N. Tabouce [5] A. Blouza, F. Hect and H. Le Dret, Two finite element approximations of Nagdi s sell model in Cartesian coordinates, SIAM J. Numer. Anal. 44 (2006), no. 2, [6] A. Blouza and H. Le Dret, Nagdi s sell model: existence, uniqueness and continuous dependence on te midsurface, J. Elasticity 64 (2001), no. 2-3, [7] C. Carstensen, Residual-based a posteriori error estimate for a nonconforming Reissner-Mindlin plate finite element, SIAM J. Numer. Anal. 39 (2002), no. 6, [8] C. Carstensen and J. Scöberl, Residual-based a posteriori error estimate for a mixed Reißner-Mindlin plate finite element metod, Numer. Mat. 103 (2006), no. 2, [9] P. J. Frey and P.-L. George, Maillages, applications aux éléments finis, Hermès, [10] F. Hect and O. Pironneau, FreeFem++, [11] R. Verfürt, A Review of a posteriori error estimation and adaptive mesrefinement tecniques, Wiley and Teubner, 1996.
A PRIORI AND A POSTERIORI OF A LINEAR ELLIPTIC PROBLEM WITH DYNAMICAL BOUNDARY CONDITION.
A PRIORI AND A POSTERIORI OF A LINEAR ELLIPTIC PROBLEM WITH DYNAMICAL BOUNDARY CONDITION. Toufic El Arwadi, Séréna Dib, Toni Saya To cite tis version: Toufic El Arwadi, Séréna Dib, Toni Saya. A PRIORI
More informationEXISTENCE OF SOLUTIONS FOR A CLASS OF VARIATIONAL INEQUALITIES
Journal of Matematics and Statistics 9 (4: 35-39, 3 ISSN: 549-3644 3 doi:.3844/jmssp.3.35.39 Publised Online 9 (4 3 (ttp://www.tescipub.com/jmss.toc EXISTENCE OF SOLUTIONS FOR A CLASS OF ARIATIONAL INEQUALITIES
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 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 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 information3. Using your answers to the two previous questions, evaluate the Mratio
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0219 2.002 MECHANICS AND MATERIALS II HOMEWORK NO. 4 Distributed: Friday, April 2, 2004 Due: Friday,
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
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 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 informationPart VIII, Chapter 39. Fluctuation-based stabilization Model problem
Part VIII, Capter 39 Fluctuation-based stabilization Tis capter presents a unified analysis of recent stabilization tecniques for te standard Galerkin approximation of first-order PDEs using H 1 - conforming
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 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 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 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 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 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 information, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim
Mat 132 Differentiation Formulas Stewart 2.3 So far, we ave seen ow various real-world problems rate of cange and geometric problems tangent lines lead to derivatives. In tis section, we will see ow to
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 informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationGuaranteed and robust a posteriori bounds for Laplace eigenvalues and eigenvectors: a unified framework
Guaranteed and robust a posteriori bounds for Laplace eigenvalues and eigenvectors: a unified framework Eric Cancès, Geneviève Dusson, Yvon Maday, Benjamin Stamm, Martin Voralík To cite tis version: Eric
More informationGradient Descent etc.
1 Gradient Descent etc EE 13: Networked estimation and control Prof Kan) I DERIVATIVE Consider f : R R x fx) Te derivative is defined as d fx) = lim dx fx + ) fx) Te cain rule states tat if d d f gx) )
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 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 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 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 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 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 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 informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
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 informationFE approximation for an hybrid Naghdi equations for shells with G 1 -midsurface
FE approximation for an hybrid Naghdi equations for shells with G 1 -midsurface Barbouche Refka To cite this version: Barbouche Refka. FE approximation for an hybrid Naghdi equations for shells with G
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 informationMore on generalized inverses of partitioned matrices with Banachiewicz-Schur forms
More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,
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 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 informationON A PROBLEM OF SHALLOW WATER TYPE
24-Fez conference on Differential Equations and Mecanics Electronic Journal of Differential Equations, Conference 11, 24, pp. 19 116. ISSN: 172-6691. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu
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 informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
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 informationLecture 10: Carnot theorem
ecture 0: Carnot teorem Feb 7, 005 Equivalence of Kelvin and Clausius formulations ast time we learned tat te Second aw can be formulated in two ways. e Kelvin formulation: No process is possible wose
More information6. Non-uniform bending
. Non-uniform bending Introduction Definition A non-uniform bending is te case were te cross-section is not only bent but also seared. It is known from te statics tat in suc a case, te bending moment in
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
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 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 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 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 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 informationA SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION
A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable
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 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 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 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 informationEXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD
MATHEMATICS OF COMPUTATION Volume 68, Number 228, Pages 1379 1396 S 0025-5718(99)01093-5 Article electronically publised on February 24, 1999 EXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationNEW UNKNOWNS ON THE MIDSURFACE OF A NAGHDI S SHELL MODEL
An. Şt. Univ. Ovidius Constanţa Vol. 12(2), 2004, 115 126 NEW UNKNOWNS ON THE MIDSURFACE OF A NAGHDI S SHELL MODEL Liliana Gratie To Professor Dan Pascali, at his 70 s anniversary Abstract In this work,
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
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 informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
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 informationTHE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225
THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:
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 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 informationHigher Derivatives. Differentiable Functions
Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.
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 informationA MULTILEVEL PRECONDITIONER FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD
A MULTILEVEL PRECONDITIONER FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD KOLJA BRIX, MARTIN CAMPOS PINTO, AND WOLFGANG DAHMEN Abstract. In tis article we present a multilevel preconditioner for
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 informationAndrea Braides, Anneliese Defranceschi and Enrico Vitali. Introduction
HOMOGENIZATION OF FREE DISCONTINUITY PROBLEMS Andrea Braides, Anneliese Defrancesci and Enrico Vitali Introduction Following Griffit s teory, yperelastic brittle media subject to fracture can be modeled
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 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 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 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 informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationA Weak Galerkin Method with an Over-Relaxed Stabilization for Low Regularity Elliptic Problems
J Sci Comput (07 7:95 8 DOI 0.007/s095-06-096-4 A Weak Galerkin Metod wit an Over-Relaxed Stabilization for Low Regularity Elliptic Problems Lunji Song, Kaifang Liu San Zao Received: April 06 / Revised:
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 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 informationEfficient algorithms for for clone items detection
Efficient algoritms for for clone items detection Raoul Medina, Caroline Noyer, and Olivier Raynaud Raoul Medina, Caroline Noyer and Olivier Raynaud LIMOS - Université Blaise Pascal, Campus universitaire
More informationTHE INF-SUP STABILITY OF THE LOWEST ORDER TAYLOR-HOOD PAIR ON ANISOTROPIC MESHES arxiv: v1 [math.na] 21 Oct 2017
THE INF-SUP STABILITY OF THE LOWEST ORDER TAYLOR-HOOD PAIR ON ANISOTROPIC MESHES arxiv:1710.07857v1 [mat.na] 21 Oct 2017 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL Abstract. Uniform LBB conditions are
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 informationSolving Continuous Linear Least-Squares Problems by Iterated Projection
Solving Continuous Linear Least-Squares Problems by Iterated Projection by Ral Juengling Department o Computer Science, Portland State University PO Box 75 Portland, OR 977 USA Email: juenglin@cs.pdx.edu
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 informationarxiv: v1 [math.na] 19 Mar 2018
A primal discontinuous Galerkin metod wit static condensation on very general meses arxiv:1803.06846v1 [mat.na] 19 Mar 018 Alexei Lozinski Laboratoire de Matématiques de Besançon, CNRS UMR 663, Univ. Bourgogne
More informationA domain decomposition convergence for elasticity equations
A domain decomposition convergence for elasticity equations A. Cakib a, A. Ellabib b, A. Nacaoui c, a Département de Matématiques Appliquées et Informatique FST de Beni-Mellal, Université Sultan My Slimane
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 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 informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips
More informationAn Implicit Finite Element Method for the Landau-Lifshitz-Gilbert Equation with Exchange and Magnetostriction
SMA An Implicit Finite Element Metod for te Landau-Lifsitz-Gilbert Equation wit Excange and Magnetostriction Jonatan ROCHAT under te supervision of L. Banas and A. Abdulle 2 ACKNOWLEDGEMENTS I would like
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
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 informationDerivatives of Exponentials
mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function
More informationChapter 10. Function approximation Function approximation. The Lebesgue space L 2 (I)
Capter 1 Function approximation We ave studied metods for computing solutions to algebraic equations in te form of real numbers or finite dimensional vectors of real numbers. In contrast, solutions to
More informationDecay of solutions of wave equations with memory
Proceedings of te 14t International Conference on Computational and Matematical Metods in Science and Engineering, CMMSE 14 3 7July, 14. Decay of solutions of wave equations wit memory J. A. Ferreira 1,
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More informationUniform estimate of the constant in the strengthened CBS inequality for anisotropic non-conforming FEM systems
Uniform estimate of te constant in te strengtened CBS inequality for anisotropic non-conforming FEM systems R. Blaeta S. Margenov M. Neytceva Version of November 0, 00 Abstract Preconditioners based on
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
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 informationAnalytic Functions. Differentiable Functions of a Complex Variable
Analytic Functions Differentiable Functions of a Complex Variable In tis capter, we sall generalize te ideas for polynomials power series of a complex variable we developed in te previous capter to general
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 information