Chemnitz Scientific Computing Preprints
|
|
- Shanon Price
- 5 years ago
- Views:
Transcription
1 G. Of G. J. Rodin O. Steinbac M. Taus Coupling Metods for Interior Penalty Discontinuous Galerkin Finite Element Metods and Boundary Element Metods CSC/11-02 Cemnitz Scientific Computing Preprints
2 Impressum: Cemnitz Scientific Computing Preprints ISSN ( : Preprintreie des Cemnitzer SFB393) Herausgeber: Professuren für Numerisce und Angewandte Matematik an der Fakultät für Matematik der Tecniscen Universität Cemnitz Postanscrift: TU Cemnitz, Fakultät für Matematik Cemnitz Sitz: Reicenainer Str. 41, Cemnitz ttp://
3 Cemnitz Scientific Computing Preprints G. Of G. J. Rodin O. Steinbac M. Taus Coupling Metods for Interior Penalty Discontinuous Galerkin Finite Element Metods and Boundary Element Metods CSC/11-02 Abstract Tis paper presents tree new coupling metods for interior penalty discontinuous Galerkin finite element metods and boundary element metods. Te new metods allow one to use discontinuous basis functions on te interface between te subdomains represented by te finite element and boundary element metods. Tis feature is particularly important wen discontinuous Galerkin finite element metods are used. Error and stability analysis is presented for some of te metods. Numerical examples suggest tat all tree metods exibit very similar convergence properties, consistent wit available teoretical results. CSC/11-02 ISSN September 2011
4 Contents 1. Introduction 1 2. Model Problem and Background Boundary Integral Equations Interior Penalty DG Finite Element Metods Coupling Metods New Coupling Metods First Metod: Non Symmetric Coupling Second Metod: Symmetric Tree Field Approac Tird Metod: Diriclet Based Coupling Stability and Error Analysis Numerical Examples Summary 18 A. Appendix 21 Autor s addresses: Günter Of Fakultät für Matematik, Tecnisce Universität Cemnitz, Reicenainer Str. 41, Cemnitz, Germany. of@tugraz.at Olaf Steinbac Institut für Numerisce Matematik, Tecnisce Universität Graz, Steyrergasse 30, 8010 Graz, Austria. o.steinbac@tugraz.at Gregory J. Rodin Institute for Computational Engineering and Sciences, Te University of Texas at Austin, Austin, TX, 78712, USA. gjr@ices.utexas.edu Mattias Taus Institute for Computational Engineering and Sciences, Te University of Texas at Austin, Austin, TX, 78712, USA. taus@ices.utexas.edu
5 1. Introduction Tis paper is concerned wit coupling metods for finite element and boundary element metods. Suc coupling metods are advantageous for problems wose domain involves an interior finite subdomain(s) embedded in an exterior unbounded subdomain, suc tat in te interior subdomain te governing partial differential equations are complex and require finite element metods, wereas in te exterior subdomain te governing partial differential equations are simple and can be solved using boundary element metods. Te coupling metods are well establised in te literature for classical (continuous) finite element and boundary element metods; we refer to [11] and te references given terein. Tis paper is motivated by applications tat require discontinuous Galerkin (DG) rater tan continuous finite element metods. Coupling metods involving DG finite element metods ave been analyzed by Bustinza, Gatica, Heuer, and Sayas [2, 3, 8, 9], wo establised tat essentially any boundary element metod can be combined wit any DG finite element metod, as long as one uses approximations continuous on te interface. For two dimensional problems, tis restriction can be removed if one combines DG finite element metods wit a particular Galerkin boundary element metod [8]. Te coupling metods considered by Bustinza, Gatica, Heuer, and Sayas are based on te symmetric formulation of boundary integral equations. In tis case, unique solvability of te coupling metod is a direct consequence of te unique solvability of te underlying finite element and boundary element systems. A disadvantage of symmetric boundary element metods is tat tey involve te ypersingular boundary integral operator tat not only precludes te use of basic collocation scemes but also requires functions continuous on te interface. Te latter restriction is particularly undesirable for coupling metods involving DG finite element metods in R 3 [8]. In tis paper, we present tree new metods tat allow for discontinuous functions on te interface, and terefore significantly simplify te coupling between DG finite element metods wit eiter Galerkin or collocation boundary element metods. Te first metod is based on te Jonson Nédélec coupling [13] extended to DG finite element metods. Tis metod admits bot collocation and Galerkin boundary element metods. However, te metod gives rise to non symmetric linear algebraic problems and its matematical foundations ave not been establised. Te second metod, wic combines a tree field approac [1] and a symmetric boundary integral formulation, addresses some of te drawbacks of te first metod, but it involves te ypersingular operator, and terefore it is limited to Galerkin boundary element metods. Tis metod gives rise to non symmetric but well structured linear algebraic problems tat can be solved almost as efficiently as symmetric ones. Following te coupling approac of DG and mixed finite element metods [10] te tird metod gives one two options. Te first option admits bot collocation and Galerkin scemes and results in non symmetric linear algebraic problems. Tis option as te advantage of aving a sound matematical foundation for te Galerkin sceme. Te second option is limited to Galerkin boundary element metods, but it results in symmetric 1
6 linear algebraic problems and as a sound matematical foundation. Te rest of te paper is organized as follows. In Section 2, we introduce a model problem and briefly outline relevant existing results necessary for presenting te new metods. In Section 3, we present te new coupling metods. In Section 4, we establis unique solvability and error estimates for one of te metods. In Section 5, we present numerical results indicating tat te proposed metods and teir variants ave very similar convergence properties, and tose properties are consistent wit available teoretical results. Te paper is concluded wit a brief summary. 2. Model Problem and Background For a bounded domain Ω R 3 wit a Lipscitz boundary Γ := Ω and a given volume density f L 2 (Ω), te model problem involves te partial differential equations te transmission conditions u i = f in Ω, u e = 0 in Ω c := R 3 \Ω, (2.1) u := u e u i = 0 and n u = n ( u e u i ) = 0 on Γ, (2.2) and te radiation condition ( ) 1 u e (x) = O x as x. (2.3) Here n denotes te outward unit normal vector on Γ. For our purposes, it is expedient to decouple te stated boundary value problem (2.1) (2.3) into an interior boundary value problem for te subdomain Ω and an exterior boundary value problem for te subdomain Ω c. We suppose tat te numerical treatment of te interior and exterior problems is based on DG finite element and boundary element metods, respectively. Our objective is to identify appropriate interior and exterior boundary value problems, boundary integral equations for te exterior problems, and discretization scemes on Γ Boundary Integral Equations Te Caucy data u e Γ and t e := (n u e ) Γ uniquely define a armonic function u e (x) for x Ω c via te representation formula, e.g. [24], u e (x) = U (x, y)t e (y)ds y + U (x, y)u e (y)ds y, (2.4) n y Γ Γ 2
7 wic satisfies te radiation condition (2.3), and were U (x, y) = 1 1 4π x y is te fundamental solution of te Laplace operator. Te Caucy data u e Γ and t e can be related to eac oter using boundary integral equations on Γ, all of wic follow from te exterior Calderon projection Here for x Γ (V t e )(x) = (K t e )(x) = Γ Γ ( ue t e ) ( 1 2 = I + K V ) ( ue 1 D I 2 K t e U (x, y)t e (y)ds y, (Ku e )(x) = n x U (x, y)t e (y)ds y, Γ ). (2.5) n y U (x, y)u e (y)ds y, (Du e )(x) = U (x, y)u e (y)ds y n x Γ n y denote te single layer, double layer, adjoint double layer, and te ypersingular boundary integral operators, respectively. Te mapping properties of tese operators are well establised, e.g. [7, 12, 14, 21, 24]. In particular, te single layer operator V : H 1/2 (Γ) H 1/2 (Γ) is bounded and H 1/2 (Γ) elliptic, and terefore invertible. Hence equations (2.5) imply te Diriclet to Neumann map t e = V 1 ( 1 2 I K)u e = [D + ( 12 I K )V 1 ( 12 ] I K) u e =: S ext u e. (2.6) Let us note tat bot representations of te Steklov Poincaré operator S ext : H 1/2 (Γ) H 1/2 (Γ) are self adjoint in te continuous setting. However, tese representations may ave different stability and symmetry properties in te discrete setting, e.g. [23]. Finally, let us mention tat te bilinear form induced by te ypersingular operator D allows for an alternative representation tat involves weakly singular surface integrals only [16]: Du, v Γ = 1 4π Γ Γ curl Γ u(y) curl Γ v(x) ds y ds x for all u, v H 1/2 (Γ) C(Γ), (2.7) x y were curl Γ is te surface curl operator. Tis representation is central to Galerkin boundary element metods involving te ypersingular operator D, as it allows one to represent te action of D in terms of te single layer operator V. Let us empasize tat (2.7) olds for continuous densities u and v only; oterwise (2.7) must include additional terms. 3
8 2.2. Interior Penalty DG Finite Element Metods In te proposed coupling metods, DG finite element metods are restricted to interior penalty metods [19, 20, 27] wic are well studied and widely applied. Let T = {T l } N Ω l=1 be a finite element mes of te interior domain Ω. For an element T T, we identify te boundary T, te diameter T, and te outward unit normal vector n T. We define te global mes size := max T T T. Te interior and exterior element faces of te finite element mes are defined as and E int := {e : T α, T β T : e = T α T β, α β}, E ext := {e : T T : e = T Γ}, respectively. For an interior face e = T α T β wit α < β, te jump and te average values of an element wise smoot function φ are defined as φ e := (φ Tα ) e (φ Tβ ) e and {φ} e := 1 ) ((φ Tα ) e + (φ Tβ ) e, 2 respectively, and te diameter of te face e is denoted by e. For s > 3 2 we introduce te broken Sobolev space V := { v L 2 (Ω) : v T H s (T ) T T }, and te semi discrete bilinear form for u, v V a DG (u, v) := u(x) v(x) dx T T {n u} e (x) v e (x) ds x (2.8) T e E int T e ξ {n v} e (x) u e (x) ds x + u e (x) v e (x) ds x, T e T e e e e E int e E int were ξ is a formulation parameter. In particular, te values ξ { 1, 0, 1} correspond to non symmetric, incomplete, and symmetric interior penalty DG finite element metods, respectively. Te parameters > 0 are required for stabilization. Te related energy norm is given by v 2 DG := v 2 L 2 (T ) + v 2 L 2 (e). (2.9) e T T Let e E int V := { v L 2 (Ω) : v T P p (T ) T T } = span{ϕi } M i=1 (2.10) denote te standard finite element space of local polynomials of degree p. For te coupling wit boundary element metods it is useful to consider a splitting of V = V Ω VΓ wit V Ω = span{ϕ i : ϕ i Γ = 0} M Ω i=1 and V Γ = span{ϕ i } M i=m Ω +1, (2.11) 4
9 were V Ω is te space spanned by te degrees of freedom in te interior of Ω, and V Γ te space spanned by te degrees of freedom on te interface Γ. In addition, let is Ṽ Γ = span { ϕ i } M M Ω +1 VΓ be te subspace of boundary basis functions continuous on Γ Coupling Metods In tis section, we briefly describe te coupling metod of Gatica, Heuer, and Sayas [8]. We regard tis metod as te current state of te art for coupling DG finite element and boundary element metods. To tis end, we consider te interior Neumann boundary value problem u i (x) = f(x) for x Ω, u i (x) = t i (x) for x Γ. n x In te context of te interior penalty metods, tis problem results in te variational problem of finding u i, V suc tat a DG (u i,, v ) = Ω f(x)v (x)dx + Γ t i (x)v (x)ds x for all v V. (2.12) By using te Neumann transmission condition, t i = t e, and te Diriclet to Neumann map (2.6), t e = S ext u e, we obtain a DG (u i,, v ) + S ext u e, v Γ = f, v Ω, u i = u e. As in te symmetric coupling of classical finite element and boundary element metods [6, 23], we use te symmetric representation (2.6) of te Steklov Poincaré operator S ext to obtain a DG (u i,, v ) + Du e, v Γ ( 1 2 I K )t e, v Γ = f, v Ω (2.13) were t e = V 1 ( 1 2 I K)u e H 1/2 (Γ) is te unique solution of te variational problem V t e, w Γ + ( 1 2 I K)u e, w Γ = 0 for all w H 1/2 (Γ). (2.14) For a Galerkin discretization of te variational problem (2.13) and (2.14), we introduce a finite dimensional ansatz space W = span{ψ k } N Γ k=1 H 1/2 (Γ), 5
10 and approximate u e by te Diriclet trace of u i, on Γ. Tis results in te variational problem of finding (u i,, t e, ) (V Ω ṼΓ ) W suc tat for all v V Ω ṼΓ, and a DG (u i,, v ) + Du i,, v Γ ( 1 2 I K )t e,, v Γ = f, v Ω (2.15) V t e,, w Γ + ( 1 2 I K)u i,, w Γ = 0 for all w W. (2.16) Tis variational problem was first proposed and analyzed in [8]. Since te ypersingular operator D requires te use of continuous basis functions, one must use te subspace V Ω ṼΓ instead of te general space V. Altoug tis restriction guarantees unique solvability and leads to optimal error estimates, it is incompatible wit te spirit of DG finite element metods. Furtermore, constructing te restricted subspace poses significant practical difficulties, especially for tree-dimensional problems, and terefore te entire approac may not be appealing to practitioners. Tis issue may be addressed by introducing additional Lagrange multipliers tat ensure continuity of u on Γ [8]. However, tis modification seems to be also cumbersome for tree-dimensional problems. Te variational problem (2.15) (2.16) is equivalent to te system of linear algebraic equations were K DG ΩΩ K DG ΩΓ K DG ΓΩ K ΓΓ DG ( 1 M 2 K ) 1 2 K V u Ω i u Γ i t e = f Ω f Γ 0, (2.17) KΩΩ DG[j, = a DG(ϕ i, ϕ j ) for i, j = 1,..., M Ω, K ΓΩ DG[j, = a DG( ϕ MΩ +i, ϕ j ) for i = 1,..., M M Ω, j = 1,..., M Ω, K ΩΓ DG[j, = a DG(ϕ i, ϕ MΩ +j) for i = 1,..., M Ω, j = 1,..., M M Ω, K ΓΓ DG[j, = a DG( ϕ MΩ +i, ϕ MΩ +j) for i, j = 1,..., M M Ω are te blocks of te DG finite element matrix. Furter, D [j, i] = D ϕ MΩ +i, ϕ MΩ +j Γ, M [l, i] = ϕ MΩ +i, ψ l Γ, V [l, k] = V ψ k, ψ l Γ, K [l, i] = K ϕ MΩ +i, ψ l Γ, for i, j = 1,..., M M Ω, k, l = 1,..., N Γ are te Galerkin boundary element matrices, and te rigt and side is given by fj Ω = f(x)ϕ j (x)dx for j = 1,..., M Ω, Ω fj Γ = f(x)ϕ MΩ +j(x)dx for all j = 1,..., M M Ω. Ω 6
11 Since te discrete single layer integral operator V is symmetric and positive definite, we can eliminate te discrete exterior Neumann datum t e to obtain te Scur complement system were ( K DG ΩΩ K DG ΩΓ K DG ΓΓ K DG ΓΩ + S sym ) ( u Ω i u Γ i ) ( ) f Ω =, (2.18) S sym = D + ( 1 2 M K )V 1 (1 2 M K ) (2.19) is a symmetric Galerkin boundary element approximation of te exterior Steklov Poincaré sym operator (2.6). Note tat S is positive definite for any coice of admissible basis functions ϕ i and ψ k, e.g. [23]. f Γ 3. New Coupling Metods In tis section, we present tree new coupling metods for te interior penalty DG finite element metods and boundary element metods. All tree metods allow one to use te standard space V of globally discontinuous finite element functions, wic significantly simplifies te implementation in comparison to te coupling metods tat require functions continuous on Γ. Some of te coupling metods admit bot collocation and Galerkin boundary element metods, wic is particularly useful for practitioners First Metod: Non Symmetric Coupling Tis metod simply extends Jonson Nédélec s coupling metod [13] involving classical finite element metods to DG finite element metods. Tis metod uses only te first integral equation in (2.5), and terefore te corresponding Steklov Poincaré operator is represented as t e = S ext u e = V 1 ( 1 2 I K)u e. By combining tis equation wit te Diriclet transmission condition u i = u e, we obtain te variational problem of finding (u i,, t e, ) V W suc tat a DG (u i,, v ) t e,, v Γ = f, v Ω for all v V, (3.1) V t e,, w Γ + ( 1 2 I K)u i,, w Γ = 0 for all w W. (3.2) Since te ypersingular operator D does not appear in tese equations, we can solve tem using te standard discontinuous finite element space V. 7
12 Te variational problem (3.1) and (3.2) is equivalent to te system of linear algebraic equations were K DG ΩΩ K DG ΩΓ K DG ΓΩ K DG ΓΓ M 1 2 M K V u Ω i u Γ i t e = K DG [j, i] = a DG (ϕ i, ϕ j ), f j = f, ϕ j Ω for i, j = 1,..., M. f Ω f Γ 0, (3.3) Te blocks of te stiffness matrix and te rigt and side vector are obtained according to te splitting (2.11). Te blocks M [l, i] = ϕ MΩ +i, ψ l Γ and K [l, i] = Kϕ MΩ +i, ψ l Γ for i = 1,..., M M Ω, l = 1,..., N Γ, and te block V as been already introduced in (2.17). Tis system of linear algebraic equations does not ave any apparent symmetry properties. By eliminating te discrete Neumann datum t e we obtain te Scur complement system ( ) K DG ( ΩΩ KΓΩ DG u Ω ) ( ) i f Ω KΩΓ DG KΓΓ DG + = u Γ Ŝns,G i f Γ, (3.4) were Ŝ ns,g = M V 1 (1 2 M K ) (3.5) is a non symmetric Galerkin boundary element approximation of te exterior Steklov Poincaré operator (2.6). In contrast to te symmetric approximation (2.19), Ŝns,G is in general not positive definite, and terefore a certain stability condition is necessary for positive definiteness. Tat condition can be satisfied wit a proper coice of te basis functions ϕ i Γ and ψ k, e.g. [23, 26]. Te stability condition is not necessary wen te coupling involves classical finite element metods [22, 25]. In tis paper, we simply conjecture tat te stability condition is also not necessary wen te coupling involves DG finite element metods. In Section 5, we present numerical examples supporting tis conjecture. Te singular boundary integral equation also admits collocation discretizations. Tis results in te approximate Steklov Poincaré operator were Ŝ ns,c = M V 1 ( 1 2 M K ), (3.6) V [l, k] = (V ψ k )(x l), M [l, i] = ϕ MΩ +i(x l), K [l, i] = (Kϕ MΩ +i)(x l) for i = 1,..., M M Ω, l = 1,..., N Γ are te entries of te collocation matrices, and x l are collocation nodes. 8
13 3.2. Second Metod: Symmetric Tree Field Approac Tis metod is based on coupling of te interior penalty DG finite element metods and symmetric boundary element metods. In contrast to te first metod, tis metod involves bot integral equations in (2.5). Te advantage of tis metod is tat its stability can be proved, and te resulting system of linear algebraic equations is block skew symmetric; tis algebraic structure can be advantageously exploited. Te drawback of te metod is tat it does not admit collocation scemes. Te metod as te structure similar to tat of te tree field domain decomposition metod of Brezzi and Marini [1]. Like in te first metod, we combine te variational problem (2.12) wit te Neumann transmission condition t i = t e : a DG (u i,, v ) t e, v Γ = f, v Ω for all v V. In contrast to te first metod, we insert te Diriclet transmission condition u i = u e into te first boundary integral equation in (2.5), but we do not use tis equation to eliminate u e : u i = u e = ( 1 2 I + K)u e V t e on Γ. To close te system of governing equations, we use te ypersingular boundary integral equation Du e + ( 1 2 I + K )t e = 0 on Γ. To proceed furter, we need to address te fact tat te ypersingular operator D and te double layer operator 1 I + K ave non trivial kernels. Terefore te exterior Diriclet 2 trace u e is not uniquely determined by te two boundary integral equations. To tis end we recall tat bot kernels coincide: u 0 (x) = 1, (Du 0 )(x) = 1 2 u 0(x) + (Ku 0 )(x) = 0 for x Γ a.e. Accordingly, we introduce te scaling condition u e, 1 Γ = 0, and te stabilized ypersingular operator D s defined as [17] D s u, v Γ = Du, v Γ + u, 1 Γ v, 1 Γ for all u, v H 1/2 (Γ). In addition to te trial space V, we introduce an ansatz space Q = span{φ i } M Γ i=1 H1/2 (Γ) C(Γ) 9
14 of continuous basis functions φ i. As a result we arrive at te variational problem of finding (u i,, t e,, u e, ) V W Q suc tat a DG (u i,, v ) t e,, v Γ = f, v Ω for all v V, (3.7) u i,, w Γ + V t e,, w Γ ( 1 2 I + K)u e,, w Γ = 0 for all w W, (3.8) ( 1 2 I + K )t e,, q Γ + D s u e,, q Γ = 0 for all q Q. (3.9) Tis variational problem is equivalent to te system of linear algebraic equations K DG ΩΩ K DG ΩΓ K DG ΓΩ K DG ΓΓ M M V 1M 2 K K D s, u Ω i u Γ i t e u e = were in addition to tose block matrices already used in (3.3) we ave f Ω f Γ 0 0, (3.10) D s, [j, i] = D s φ i, φ j Γ, M [l, i] = φ i, ψ l Γ, K [l, i] = Kφ i, ψ l Γ for i, j = 1,..., M Γ, l = 1,..., N Γ. By eliminating u e and t e from (3.10) we obtain te Scur complement system ( ) ( K DG ΩΩ KΓΩ DG u Ω ) ( ) i f Ω KΩΓ DG KΓΓ DG + = Stree u Γ i f Γ, (3.11) wit te symmetric tree field approximation of te exterior Steklov Poincaré operator S tree = M [ V + ( 1 2 M + K )D 1 s, (1 2 M + K )] 1 M. (3.12) 3.3. Tird Metod: Diriclet Based Coupling In te first two metods, te coupling involves passing te exterior Neumann data t e to te variational problem for te interior boundary value problem. In contrast, in te tird metod, te coupling is based on te Diriclet transmission condition u i = u e. Tis approac is similar to te coupling of DG and mixed finite element metods proposed by Girault, Sun, Weeler, and Yotov [10]. Te solution u i of te model problem (2.1) (2.3) coincides wit te solution of te Diriclet boundary value problem u i = f in Ω, u i = u e on Γ. 10
15 Following [10], te discrete variational problem corresponding to tis boundary value problem is to find u i, V suc tat for all v V â DG (u i,, v ) η u e, v n Γ u e v ds = f, v Ω, (3.13) e e were â DG (u, v) := a DG (u, v) n u, v Γ + η n v, u Γ + u v ds x. (3.14) e e Here η { 1, 0, +1} is a formulation parameter, similar to te formulation parameter ξ in (2.8). Let us combine te Diriclet to Neumann map (2.6) and te Neumann transmission condition t e = t i = n u i, so tat we obtain S ext u e = n u i on Γ. By adding a stabilization term to tis equation [10] and approximating u e by a function u e, Q, we obtain te variational problem of finding (u i,, u e, ) V Q suc tat â DG (u i,, v ) η n v, u e, Γ u e, v ds x = f, v Γ (3.15) e for all v V, and Su e,, q Γ + n u i,, q Γ + e [u e, u i, ]q ds x = 0 (3.16) e e for all q Q, were S is a boundary element approximation of te exterior Steklov Poincaré operator (2.6). Te variational problem (3.15) (3.16) is equivalent to te system of linear algebraic equations were K DG ΩΩ K DG ΩΓ K DG ΓΩ K ΓΓ DG B 1, B η, S + C K DG [j, i] = â DG (ϕ i, ϕ j ) u Ω i u Γ i u e = f Ω f Γ 0 for i, j = 1,..., M, (3.17) is te stiffness matrix of te modified discontinous Galerkin finite element metod wic is obtained according to te splitting (2.11). Furter, B r, [j, i] = r n ϕ MΩ +i, φ j Γ ϕ MΩ +i φ j ds x e e 11
16 for i = 1,..., M M Ω, j = 1,..., M Γ, and C [j, i] = φ i φ j ds x for i, j = 1,..., M Γ. e e Te matrix S can represent any of te following approximations of te exterior Steklov Poincaré operator (2.6), namely eiter te symmetric Galerkin approximation S sym,g = D + ( 1 2 M K )V 1 (1 2 M K ), (3.18) or te non symmetric Galerkin approximation S ns,g or te non symmetric collocation approximation S ns,c = M V 1 (1 2 M K ), (3.19) = M V 1 ( 1 2 M K ). (3.20) Te Scur complement form of (3.17) is ( K DG ΩΩ KΓΩ DG K DG ΩΓ K DG ΓΓ + S Diriclet ) ( u Ω i u Γ i ) ( f Ω = f Γ ), (3.21) wit te discrete representation of te exterior Steklov Poincaré operator S Diriclet = B η, [ S + C ] 1 B1,. (3.22) As before, te matrix S can be represented by eiter sym,g S, or ns,g S, or S ns,c. 4. Stability and Error Analysis In tis section, we establis unique solvability and error estimates for te governing equation (3.15) (3.16) of te tird metod, wit te provision tat te approximation S of te Steklov Poincaré operator (2.6) is stable. Let us associate te variational problem (3.15) (3.16) wit te bilinear form A(u i, u e ; v, q) = â DG (u i, v i ) η n v, u e Γ u e vds x (4.1) + Su e, q Γ + n u i, q Γ + e e e [u e u i ]qds x. e 12
17 Also we define te energy norm (v, q) 2 A := T T v 2 L 2 (T ) + e E int v 2 L 2 (e) + e e v q 2 L 2 (e) + q 2 H 1/2 (Γ). (4.2) Teorem 4.1 Let S be a stable boundary element approximation of te exterior Steklov Poincaré operator (2.6), Sq, q Γ c S 1 q 2 H 1/2 (Γ) for all q Q, (4.3) and η { 1, 0, +1}. Ten for sufficiently large stability parameters, te bilinear form (4.1) is elliptic, A(v, q ; v, q ) c A 1 (v, q ) 2 A for all (v, q ) V Q, (4.4) and te variational problem (3.15) (3.16) as a unique solution. Proof. For sufficiently large, te bilinear form a DG (, ) defined in (2.8) is elliptic [19]. Tat is, tere exists a constant c A 1 > 0 independent of suc tat a DG (v, v ) c A 1 v 2 DG for all v V, (4.5) were v 2 DG := T T v 2 L 2 (T ) + e E int e v 2 L 2 (e). By inserting (3.14) in (4.1) and using (4.3) and (4.5), we obtain te inequality A(v, q ; v, q ) = a DG (v, v ) + Sq, q Γ (1 η) n v, v q Γ + One can sow (see Remark 3.1 in [10]) tat c A 1 v 2 DG + c S 1 q 2 H 1/2 (Γ) (1 η) n v, v q Γ + v q, n v Γ 1 v 2 L 8 2 (T ) T T [v q ] 2 ds x e e [v q ] 2 ds x. e e e v q 2 L 2 (e), 13
18 wic implies te estimate A(v, q ; v, q ) ( c A 1 1 η ) v 2 L 8 2 (T ) + T T ( η ) 8 {( min 1 1 η 8 e E ext ), e E int v 2 L 2 (e) e v q 2 L 2 (e) + c S 1 q 2 H 1/2 (Γ) e ( c A 1 1 η ) }, c S 8 1 (v, q ) 2 A. Tus A(v, q ; v, q ) is elliptic for η { 1, 0, 1}, if we coose sufficiently large to ensure c A 1 1 η 8 > 0. Remark 4.1 If S corresponds to te symmetric Galerkin approximation (3.18), ten te stability estimate (4.3) olds for any pair of boundary element spaces Q and W, e.g. [23]. For te non symmetric approximations (3.19) and (3.20), an additional stability condition becomes necessary. Tat condition requires properly cosen boundary element spaces Q and W. In particular, tis can be acieved by constructing W on a finer mes in comparison to tat used for constructing Q, e.g. [26]. Lemma 4.2 Let (u i, u e ) be te weak solution of te model boundary value problem (2.1) (2.3) suc tat u i H s (Ω) wit s > 3 2 and (u i,, u e, ) be a solution of te variational problem (3.15) and (3.16). Ten te perturbed Galerkin ortogonality condition takes te form A(u i u i,, u e u e, ; v, q ) = ( S S ext )u e, q Γ for all V Q. (4.6) Proof. Teorem 3.2 in [10] implies â DG (u i, v) η n v, u e Γ Te solution (u i, u e ) satisfies S ext u e, q Γ + n u i, q Γ + Hence we conclude tat A(u i, u e ; v, q) + (S ext S)u e, q, Γ = f, v Ω Te solution (u i,, u e, ) satisfies u e vds x = f, v Γ for all v V. e e [u e u i ]qds x = 0 for all q H 1/2 (Γ). e e for all (v, q) V H 1/2 (Γ). A(u i,, u e, ; v, q ) = f, v Ω for all (v, q ) V Q. For a conforming approac, V Q V H 1/2 (Γ), te last two equations imply (4.6). 14
19 Teorem 4.3 Let (u i, u e ) be te weak solution of te model boundary value problem (2.1) (2.3) suc tat u i H s (Ω) wit s > 3 2. Let (u i,, u e, ) be te unique solution of te variational problem (3.15) and (3.16), and let S satisfy te approximation property ( S S ext )u e H 1/2 (Γ) c S inf t e w H w W 1/2 (Γ) wit t e = S ext u e. (4.7) Ten tere exists a constant C > 0 suc tat te quasi optimal error estimate (u i u i,, u e u e, ) 2 A C u e q 2 H 1/2 (Γ) + u i v 2 DG + e (u i v ) n 2 L 2 (e) + olds for all (v, q ) V Q. Proof. See Appendix. 1 + (u i v ) (u e q ) 2 L 2 (e) + c S inf t e w 2 e τ W H 1/2 (Γ) Remark 4.2 Te approximation property (4.7) olds unconditionally for te Galerkin approximations (3.18) and (3.19) [23, 24]. In contrast, for te collocation approximation (3.20), (4.7) olds if te approximation is stable. At present, stability of te approximation under general conditions is an open problem. Corollary 4.4 Let all assumptions of Teorem 4.3 old, and in addition 3/2 < s 2. If V is te space of discontinuous linear finite element functions, Q is te space of piecewise continuous linear boundary element basis functions, and W is te space of piecewise constant boundary element basis functions, ten tere exists a constant C suc tat (u i u i,, u e u e, ) A C s 1 [ u i H s (Ω) + t i H s 3/2 (Γ)]. (4.8) Remark 4.3 Of course if u i H 2 (Ω), te error estimate implies linear convergence rate. Tis result is straigtforward to generalize to approximations based on iger order polynomial basis functions. Suc approximations are meaningful for sufficiently large s only. Remark 4.4 Te Scur complement systems (2.18), (3.4), (3.11) and of (3.21) allow for a unified treatment of te coupling metods. In particular, tis unified structure allows one to exploit standard results pertaining to stability and error analysis. However, stability analysis of te non symmetric formulation (3.1) (3.2) requires additional considerations. Wile it is likely tat a successful treatment of te non symmetric formulation is possible along te lines proposed in [22, 25], it is not pursued in tis paper. Here, we limit our study of tis issue to presenting numerical results confirming expected teoretical results.. 15
20 5. Numerical Examples In tis section, we present numerical results confirming te teoretical results establised for te tird metod and suggesting tat te first two metods are stable and exibit expected rates of convergence in te energy and L 2 norms. In all examples, Ω is a unit spere, and f(x) = f(r, φ, θ) = 4(cos φ + sin φ) sin θ, were r, φ, θ are te sperical coordinates wose origin is at te spere center. Te exact solution for tis problem is u(x) = ũ(r, φ, θ) = 1 3 (cos φ + sin φ) sin θ { (4 3r)r for r < 1, r 2 for r > 1. Tis function is piecewise analytic, and terefore te numerical solutions are expected to ave te optimal rates of convergence. Numerical solutions were obtained using piecewise linear discontinuous finite element basis functions for V, piecewise linear continuous basis functions for Q, and piecewise constant basis functions for W. Te interface Γ was approximated using piecewise linear continuous basis functions; errors associated wit tis approximation can be estimated and controlled by using standard tecniques, e.g. [5, 15]. Te collocation and Galerkin boundary element metods were accelerated using te fast multipole metod; again te errors associated wit te use of multipole and local expansions were estimated and controlled following [18]. Te DG formulation parameters were cosen to be ξ = η = 1, and te DG stabilisation parameters were cosen to be = 5. In presenting results, we denote te number of finite elements T in Ω by N Ω, te number of elements on Γ by N Γ, and te number of nodes on Γ by M Γ. First, let us present numerical results for te tird metod (3.15) (3.16) in wic te exterior Steklov Poincaré operator is discretized using eiter te symmetric Galerkin sceme (3.18), or te non symmetric Galerkin sceme (3.19), or te non symmetric collocation sceme (3.20). In all cases, te error norm was computed as te modified energy norm [ê DG ()] 2 := T T (u i u i, ) 2 L 2 (T ) + + e E int (u i u i, ) (u e u e, ) 2 L 2 (e) + e e u i u i, 2 L 2 (e) 1 e u e u e, 2 L 2 (e) were te H 1/2 (Γ) norm in (4.2) was replaced by a weigted L 2 (Γ) norm, e.g. [4]. Te estimated order of convergence was computed as ( ) êdg ( l ) eoc := log 2, ê DG ( l+1 ) 16
21 were l refers to te refinement level, and l = 2 l+1. Table 1 contains numerical results for six meses (five refinement levels) for all tree scemes. Clearly te results support te notion tat te asymptotic order of convergence is linear. Table 2 mimics Table 1 except it is based on te L 2 norm e 2 () = u i u i, L2 (Ω) rater tan te energy norm. As expected, te asymptotic order of convergence is quadratic. Galerkin (3.18) Galerkin (3.19) Collocation (3.20) N Ω N Γ M Γ ê DG eoc ê DG eoc ê DG eoc Table 1: Errors and rates of convergence for te tird metod measured using te energy norm. Galerkin (3.18) Galerkin (3.19) Collocation (3.20) N Ω N Γ M Γ e 2 eoc e 2 eoc e 2 eoc Table 2: Errors and rates of convergence for te tird metod measured using te L 2 norm. Numerical results for te tree field formulation (3.7) (3.9) are given in Table 3, were we use te DG energy norm e DG () := u i u i, DG and te L 2 norm. Again we observe te expected linear and quadratic orders of convergence, respectively. Finally, Table 4 contains numerical results for te non symmetric approac (3.1) (3.2) for te Galerkin (3.5) and collocation (3.6) scemes. Numerical results indicate tat all versions of te tree metods performed similarly to eac oter. In particular, all of tem exibit linear order of convergence in te pertinent 17
22 Galerkin (3.10) N Ω N Γ M Γ e DG eoc e 2 eoc Table 3: Errors and rates of convergence for te second metod measured using te energy norm and L 2 norms. Galerkin (3.5) Collocation (3.6) N Ω N Γ M Γ e DG eoc e 2 eoc e DG eoc e 2 eoc Table 4: Errors and rates of convergence for te first metod measured using te energy norm and L 2 norms. energy norm and quadratic order of convergence in te L 2 (Ω) norm. Furtermore, te differences in te absolute values of te corresponding errors appear to be minimal. 6. Summary Tis paper introduces tree new coupling metods for te interior penalty DG finite element metods and boundary element metods. Te key advantage of te new metods is tat tey allow one to use discontinuous basis functions on te interface, wic is particularly important for coupling metods involving DG finite element metods. Te new coupling metods ave six variations associated wit different approximations of te Steklov Poincaré operator of te underlying boundary integral equations. We presented teoretical results pertaining to stability and error analysis for some of te versions of te metods, wereas establising suc results for oter versions can be done in a similar way. Numerical results suggest tat all versions perform very similar to eac oter, and exibit expected rates of convergence in bot energy and L 2 norms. 18
23 Acknowledgement Te autors would like to tank M. Neumüller for providing te DG solver wic was used for te DG part of all computations reported in tis paper. GJR acknowledges te financial support from te Department of Energy (grant DOE DE-FG02-05ER25701) and by te W. A. Tex Moncrief, Jr. Endowment In Simulation Based Engineering Sciences troug a Grand Callenge Faculty Fellowsip. MT acknowledges te financial support from te Fulbrigt Foundation. References [1] F. Brezzi, L. D. Marini: A tree field domain decomposition metod. Contemp. Mat. 157 (1994) [2] R. Bustinza, G. N. Gatica, F. J. Sayas: A LDG BEM coupling for a class of nonlinear exterior transmission problems. In: Numerical Matematics and Advanced Applications, Springer, Berlin, pp , [3] R. Bustinza, G. N. Gatica, F. J. Sayas: On te coupling of local discontinuous Galerkin and boundary element metods for non linear exterior transmission problems. IMA J. Numer. Anal. 28 (2008) [4] C. Carstensen, M. Maiscak, D. Praetorius, E. P. Stepan: Residual based a posteriori error estimate for ypersingular equation on surfaces. Numer. Mat. 97 (2004) [5] P. G. Ciarlet: Te Finite Element Metod for Elliptic Problems. Nort Holland, [6] M. Costabel: Symmetric metods for te coupling of finite elements and boundary elements. In: Boundary Elements IX (C. A. Brebbia, G. Kun, W. L. Wendland eds.), Springer, Berlin, pp , [7] M. Costabel: Boundary integral operators on Lipscitz domains: elementary results. SIAM J. Mat. Anal. 19 (1988) [8] G. N. Gatica, N. Heuer, F. J. Sayas: A direct coupling of local discontinuous Galerkin and boundary element metods. Mat. Comp. 79 (2010) [9] G. N. Gatica, F. J. Sayas: An a priori error analysis for te coupling of local discontinuous Galerkin and boundary element metods. Mat. Comp. 75 (2006) [10] V. Girault, S. Sun, M. F. Weeler, I. Yotov: Coupling discontinuous Galerkin and mixed finite element discretizations using mortar finite elements. SIAM J. Numer. Anal. 46 (2008)
24 [11] G. C. Hsiao, E. Scnack, W. L. Wendland: Hybrid coupled finite boundary element metods for elliptic systems of second order. Comput. Metods Appl. Mec. Engrg. 190 (2000) [12] G. C. Hsiao, W. L. Wendland: Boundary integral equations. Applied Matematical Sciences, vol. 164, Springer, Berlin, [13] C. Jonson, J. C. Nédélec: On te coupling of boundary integral and finite element metods. Mat. Comp. 35 (1980) [14] W. McLean: Strongly Elliptic Systems and Boundary Integral Equations. Cambridge University Press, [15] J. C. Nédélec: Curved finite element metods for te solution of singular integral equations on surfaces in R 3. Comp. Met. Appl. Mec. Engrg. 8 (1976) [16] J. C. Nédélec: Integral equations wit nonintegrable kernels. Int. Eq. Operator T. 5 (1982) [17] G. Of, O. Steinbac: A fast multipole boundary element metod for a modified ypersingular boundary integral equation. In: Analysis and Simulation of Multifield Problems (M. Efendiev, W. L. Wendland eds.), Lecture Notes in Applied and Computational Mecanics, vol. 12, Springer, Berlin, pp , [18] G. Of, O. Steinbac, W. L. Wendland: Te fast multipole metod for te symmetric boundary integral formulation. IMA J. Numer. Anal. 26 (2006) [19] B. Rivière: Discontinuous Galerkin metods for solving elliptic and parabolic equations. Frontiers in Applied Matematics, vol. 35, SIAM, Piladelpia, [20] B. Rivière, M. F. Weeler, V. Girault: Improved energy estimates for interior penalty, constrained and discontinuous Galerkin metods for elliptic problems. I. Comput. Geosci. 3 (2000) [21] S. A. Sauter, C. Scwab: Boundary Element Metods. Springer Series in Computational Matematics, vol. 39. Springer, Berlin, [22] F. J. Sayas: Te validity of Jonson Nedelec s BEM FEM coupling on polygonal interfaces. SIAM J. Numer. Anal. 47 (2009) [23] O. Steinbac: Stability estimates for ybrid coupled domain decomposition metods. Lecture Notes in Matematics, vol. 1809, Springer, Heidelberg, [24] O. Steinbac: Numerical Approximation Metods for Elliptic Boundary Value Problems. Finite and Boundary Elements. Springer, New York, [25] O. Steinbac: A note on te stable one equation coupling of finite and boundary elements. SIAM J. Numer. Mat. 49 (2011)
25 [26] W. L. Wendland: On asymptotic error estimates for te combined BEM an FEM. In: In Finite Element and Boundary Element Tecniques from Matematical and Engineering Point of View (E. Stein, W. L. Wendland eds.), CISM Lecture Notes, vol. 301, Springer, Wien, pp , [27] M. F. Weeler: An elliptic collocation finite element metod wit interior penalties. SIAM J. Numer. Anal. 15 (1978) A. Appendix In tis appendix we prove te error estimate given in Teorem 4.3. For arbitrary (v, q ) V Q, te triangle inequality implies (u i u i,, u e u e, ) A (u i v, u e q ) A + (u i, v, u e, q ) A. Te first term can be furter estimated by using standard tecniques, ence it remains to bound te second term. We use te ellipticity estimate (4.4) and te perturbed Galerkin ortogonality (4.6) to obtain c A 1 (u i, v, u e, q ) 2 A A(u i, v, u e, q ; u i, v, u e, q ) = A(u i, u i, u e, u e ; u i, v, u e, q ) + A(u i v, u e q ; u i, v, u e, q ) = (S ext S)u e, u e, q Γ + A(u i v, u e q ; u i, v, u e, q ). By combining tis equation wit (4.1) we obtain c A 1 (u i, v, u e, q ) 2 A (S ext S)u e, u e, q Γ + A(u i v, u e q ; u i, v, u e, q ) = (S ext S)u e, u e, q Γ + S(u e q ), u e, q Γ + a DG (u i v, u i, v ) n (u i v ), (u i, v ) (u e, q ) Γ + η n (u i, v ), (u i v ) (u e q ) Γ + [(u i v ) (u e q )][(u i, v ) (u e, q )]ds x e e 21
26 (S ext S)u e H 1/2 (Γ) u e, q H 1/2 (Γ) + c S 2 u e q H 1/2 (Γ) u e, q H 1/2 (Γ) +c A 2 u i v DG u i, v DG + n (u i v ) L 2 (e) (u i, v ) (u e, q ) L 2 (e) + + n (u i, v ) L 2 (e) (u i v ) (u e q ) L 2 (e) e (u i v ) (u e q ) L 2 (e) (u i, v ) (u e, q ) L 2 (e). By applying te weigted Hölder inequality we obtain c A 1 (u i, v, u e, q ) 2 A [ c (S ext S)u e 2 H 1/2 (Γ) + u e q 2 H 1/2 (Γ) + u i v 2 DG + e E ext e n (u i v ) 2 L 2 (e) e (u i v ) (u e q ) 2 L 2 (e) u e, q 2 H 1/2 (Γ) + u i, v 2 DG + e n (u i, v ) 2 L 2 (e) e (u i, v ) (u e, q ) 2 L 2 (e) Let T e T be a finite element containing a face e E ext. Ten for a finite element function v V we obtain n v L 2 (e) c e 1/2 v L 2 (T e), and Hence we obtain 1/2. 1/2 e n (u i, v ) 2 L 2 (e) c T T (u i, v ) 2 L 2 (T ). u e, q 2 H 1/2 (Γ) + u i, v 2 DG e (u i, v ) (u e, q ) 2 L 2 (e) e n (u i, v ) 2 L 2 (e) 1/2 c (u i, v, u e, q ) A, 22
27 and so c A 1 (u i, v, u e, q ) A [ c (S ext S)u e 2 H 1/2 (Γ) + u e q 2 H 1/2 (Γ) + u i v 2 DG + e n (u i v ) 2 L 2 (e) e Tis inequality implies te error estimate given in Teorem 4.3. (u i v ) (u e q ) 2 L 2 (e) 1/2. 23
28 Cemnitz Scientific Computing Preprints ISSN
Technische Universität Graz
Tecnisce Universität Graz Boundary element metods for Diriclet boundary control problems Günter Of, Tan Pan Xuan, Olaf Steinbac Bericte aus dem Institut für Numerisce Matematik Berict 2009/2 Tecnisce
More informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
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 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 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 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 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 informationBOUNDARY ELEMENT METHODS FOR POTENTIAL PROBLEMS WITH NONLINEAR BOUNDARY CONDITIONS
MATHEMATICS OF COMPUTATION Volume 70, Number 235, Pages 1031 1042 S 0025-5718(00)01266-7 Article electronically publised on June 12, 2000 BOUNDARY ELEMENT METHODS FOR POTENTIAL PROBLEMS WITH NONLINEAR
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
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 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 informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More 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 informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationarxiv: v1 [math.na] 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 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 informationPriority Program 1253
Deutsce Forscungsgemeinscaft Priority Program 53 Optimization wit Partial Differential Equations K. Deckelnick, A. Günter and M. Hinze Finite element approximation of elliptic control problems wit constraints
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationDomain decomposition methods via boundary integral equations
Domain decomposition methods via boundary integral equations G. C. Hsiao a O. Steinbach b W. L. Wendland b a Department of Mathematical Sciences, University of Delaware, Newark, Delaware 19716, USA. E
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 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 informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
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 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 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 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 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 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 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 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 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 informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More 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 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 informationKey words. Sixth order problem, higher order partial differential equations, biharmonic problem, mixed finite elements, error estimates.
A MIXED FINITE ELEMENT METHOD FOR A SIXTH ORDER ELLIPTIC PROBLEM JÉRÔME DRONIOU, MUHAMMAD ILYAS, BISHNU P. LAMICHHANE, AND GLEN E. WHEELER Abstract. We consider a saddle point formulation for a sixt order
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 informationA SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages 6 626 c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR
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 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 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 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 informationDownloaded 11/15/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 55, No. 6, pp. 2787 2810 c 2017 Society for Industrial and Applied Matematics EDGE ELEMENT METHOD FOR OPTIMAL CONTROL OF STATIONARY MAXWELL SYSTEM WITH GAUSS LAW IRWIN YOUSEPT
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More 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 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 informationSMAI-JCM SMAI Journal of Computational Mathematics
SMAI-JCM SMAI Journal of Computational Matematics Compatible Maxwell solvers wit particles II: conforming and non-conforming 2D scemes wit a strong Faraday law Martin Campos Pinto & Eric Sonnendrücker
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
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 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 informationHigh-Order Extended Finite Element Methods for Solving Interface Problems
Hig-Order Extended Finite Element Metods for Solving Interface Problems arxiv:1604.06171v1 [mat.na] 21 Apr 2016 Fei Wang Yuanming Xiao Jincao Xu ey words. Elliptic interface problems, unfitted mes, extended
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 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 informationEffect of Numerical Integration on Meshless Methods
Effect of Numerical Integration on Mesless Metods Ivo Babuška Uday Banerjee Jon E. Osborn Qingui Zang Abstract In tis paper, we present te effect of numerical integration on mesless metods wit sape functions
More informationFINITE ELEMENT APPROXIMATIONS AND THE DIRICHLET PROBLEM FOR SURFACES OF PRESCRIBED MEAN CURVATURE
FINITE ELEMENT APPROXIMATIONS AND THE DIRICHLET PROBLEM FOR SURFACES OF PRESCRIBED MEAN CURVATURE GERHARD DZIUK AND JOHN E. HUTCHINSON Abstract. We give a finite element procedure for te Diriclet Problem
More informationHYBRIDIZED GLOBALLY DIVERGENCE-FREE LDG METHODS. PART I: THE STOKES PROBLEM
MATHEMATICS OF COMPUTATION Volume 75, Number 254, Pages 533 563 S 0025-5718(05)01804-1 Article electronically publised on December 16, 2005 HYBRIDIZED GLOBALLY DIVERGENCE-FREE LDG METHODS. PART I: THE
More informationarxiv: v1 [math.na] 20 Nov 2018
An HDG Metod for Tangential Boundary Control of Stokes Equations I: Hig Regularity Wei Gong Weiwei Hu Mariano Mateos Jon R. Singler Yangwen Zang arxiv:1811.08522v1 [mat.na] 20 Nov 2018 November 22, 2018
More informationA symmetric mixed finite element method for nearly incompressible elasticity based on. biorthogonal system
A symmetric mixed finite element metod for nearly incompressible elasticity based on biortogonal systems Bisnu P. Lamicane and Ernst P. Stepan February 18, 2011 Abstract We present a symmetric version
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
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 informationArbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY
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 informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
More 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 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 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 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 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 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 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 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 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 informationERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS
ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS YONGYONG CAI, AND JIE SHEN Abstract. We carry out in tis paper a rigorous error analysis
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationA-posteriori error estimates for Discontinuous Galerkin approximations of second order elliptic problems
A-posteriori error estimates for Discontinuous Galerkin approximations of second order elliptic problems Carlo Lovadina, and L. Donatella Marini Abstract Using te weigted residual formulation we derive
More informationEXTENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS WITH APPLICATION TO AN AEROACOUSTIC PROBLEM
SIAM J. SCI. COMPUT. Vol. 26, No. 3, pp. 821 843 c 2005 Society for Industrial and Applied Matematics ETENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationarxiv: v2 [math.na] 5 Feb 2017
THE EDDY CURRENT LLG EQUATIONS: FEM-BEM COUPLING AND A PRIORI ERROR ESTIMATES MICHAEL FEISCHL AND THANH TRAN arxiv:1602.00745v2 [mat.na] 5 Feb 2017 Abstract. We analyze a numerical metod for te coupled
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationA SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS
A SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS CONSTANTIN BACUTA AND KLAJDI QIRKO Abstract. We investigate new PDE discretization approaces for solving variational formulations wit different types
More informationarxiv: v1 [math.na] 6 Dec 2010
MULTILEVEL PRECONDITIONERS FOR DISCONTINUOUS GALERKIN APPROXIMATIONS OF ELLIPTIC PROBLEMS WITH JUMP COEFFICIENTS BLANCA AYUSO DE DIOS, MICHAEL HOLST, YUNRONG ZHU, AND LUDMIL ZIKATANOV arxiv:1012.1287v1
More informationarxiv: v1 [math.na] 9 Mar 2018
A simple embedded discrete fracture-matrix model for a coupled flow and transport problem in porous media Lars H. Odsæter a,, Trond Kvamsdal a, Mats G. Larson b a Department of Matematical Sciences, NTNU
More informationComputing eigenvalues and eigenfunctions of Schrödinger equations using a model reduction approach
Computing eigenvalues and eigenfunctions of Scrödinger equations using a model reduction approac Suangping Li 1, Ziwen Zang 2 1 Program in Applied and Computational Matematics, Princeton University, New
More informationTHE DISCRETE PLATEAU PROBLEM: CONVERGENCE RESULTS
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 0025-5718XX0000-0 THE DISCRETE PLATEAU PROBLEM: CONVERGENCE RESULTS GERHARD DZIUK AND JOHN E. HUTCHINSON Abstract. We solve te problem of
More informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
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 informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationDepartment of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801
RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused
More informationCOMPUTATIONAL COMPARISON BETWEEN THE TAYLOR HOOD AND THE CONFORMING CROUZEIX RAVIART ELEMENT
Proceedings of ALGORITMY 5 pp. 369 379 COMPUTATIONAL COMPARISON BETWEEN E TAYLOR HOOD AND E CONFORMING OUZEIX RAVIART ELEMENT ROLF KRAHL AND EBERHARD BÄNSCH Abstract. Tis paper is concerned wit te computational
More 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 informationAPPROXIMATION BY QUADRILATERAL FINITE ELEMENTS
MATHEMATICS OF COMPUTATION Volume 71, Number 239, Pages 909 922 S 0025-5718(02)01439-4 Article electronically publised on Marc 22, 2002 APPROXIMATION BY QUADRILATERAL FINITE ELEMENTS DOUGLAS N. ARNOLD,
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 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 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 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 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 informationHigh-Order Energy and Linear Momentum Conserving Methods for the Klein-Gordon Equation
matematics Article Hig-Order Energy and Linear Momentum Conserving Metods for te Klein-Gordon Equation He Yang Department of Matematics, Augusta University, Augusta, GA 39, USA; yang@augusta.edu; Tel.:
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 informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More information