arxiv: v1 [math.na] 8 Feb 2013

Size: px

Start display at page:

Download "arxiv: v1 [math.na] 8 Feb 2013"

Clifton Smith
5 years ago
Views:

1 A mixed method fo Diichlet poblems with adial basis functions axiv: v1 [math.na] 8 Feb 2013 Nobet Heue Thanh Tan Abstact We pesent a simple discetization by adial basis functions fo the Poisson equation with Diichlet bounday condition. A Lagangian multiplie using piecewise polynomials is used to accommodate the bounday condition. This simplifies pevious attempts to use adial basis functions in the inteio domain to appoximate the solution and on the bounday to appoximate the multiplie, which technically equies that the mesh nom in the inteio domain is significantly smalle than that on the bounday. Numeical expeiments confim theoetical esults. Key wods: scaled adial basis functions, finite elements, Diichlet condition AMS Subject Classification: 65N30, 65N12, 65N15 1 Intoduction Radial basis functions RBFs) have been used successfully [11] to solve patial diffeential equations with Neumann o Robin bounday conditions. When Diichlet conditions ae consideed they must be appoximated in an appopiate way. In the case of RBFs this causes analytical difficulties since taces of such functions have epesentations that change continuously with the position of thei centes. Theefoe, coesponding discete spaces on the bounday) do not have a clea stuctue that could be used fo analysis. Such difficulties do not appea when teating natual bounday conditions. Lagangian multiplies consisting of RBFs have been studied in [4]. Due to a weake invese popety of RBFs as compaed to that of finite elements, the condition imposed on the mesh noms in the inteio domain and on the bounday is too estictive. In this pape we suggest to use RBFs in the domain and finite elements on the bounday. We also analyze the influence of scaling of RBFs on the eo estimate. Scaling of RBFs can be used to avoid too much ovelap which is essential fo the conditioning of stiffness matices. The idea to couple RBFs with finite elements fo the Lagangian multiplie is poposed in [7] whee we solve a hypesingula integal equation. In that pape, in ode to deal with the popety that the solution of the integal equation can be extended by zeo, we have Patially suppoted by CONICYT-Chile though FONDECYT-poject and Anillo ACT1118 ANANUM), and UNSW FRG Gant numbe PS Facultad de Matemáticas, Pontificia Univesidad Católica de Chile, Santiago, Chile. mailto: nheue@mat.puc.cl School of Mathematics and Statistics, The Univesity of New South Wales, Sydney 2052, Austalia. mailto: Thanh.Tan@unsw.edu.au 1

2 to conside the Lagangian multiplie on an extended domain. In the cuent situation it suffices to use multiplies on the bounday of the domain whee the poblem is set. Moeove, not as in [7] whee the Lagangian multiplie does not have any physical meaning, in the poblem to be consideed in this pape it epesents the nomal deivative of the solution. We pove an inf-sup condition which allows an estimate fo the appoximation of this multiplie; we did not succeed in poving this condition fo the poblem consideed in [7]. The pape is oganized as follows. In Section 2 we intoduce the model poblem and the mixed vaiational fomulation. In Section 3 the finite-dimensional spaces and discetization ae intoduced. Section 4 is a evisit of appoximation popeties of RBFs whee we extend pevious esults so that they can be used in the cuent study. Section 5 pesents the main esult of the pape, namely an a pioi eo estimate fo the appoximation. This analysis is caied out afte we pove discete ellipticity of the Diichlet bilinea fom in the case that no mass tem is pesent) and a discete inf-sup condition of the bilinea fom involving the Lagangian multiplie. Section 6 pesents numeical expeiments that confim the eo estimate. Thoughout the pape, the notation a b indicates that thee exists a constant C > 0 independent of discetization o scaling paametes h X, k, and involved functions except whee othewise noted) such that a Cb. Similaly we use a b, and a b means that a b a. 2 Model poblem and mixed fomulation Conside the model poblem u + κu = f in Ω, u = g on Γ, 2.1) whee Ω R d d = 2, 3) is a polygonal d = 2) o polyhedal d = 3) domain with bounday Γ, κ 0 is a constant, and whee f L 2 Ω) and g H 1/2 Γ) ae given functions. When κ = 0 then f and g must satisfy the usual compatibility condition. The solution u of 2.1) will be found in the weak sense by using a standad mixed vaiational fomulation. Defining av, w) = v w + κvw) dx v, w H 1 Ω), Ω 2.2) bv, µ) = vµ ds v H 1 Ω), µ H 1/2 Γ), Γ we can easily see that thee exists a constant C > 0 such that av, w) C v H 1 Ω) w H 1 Ω) v, w H 1 Ω), bv, µ) C v H 1 Ω) µ H 1/2 Γ) v H 1 Ω), µ H 1/2 Γ). 2.3) A vaiational fomulation of 2.1) is fomulated as: Find u, λ) H 1 Ω) H 1/2 Γ) satisfying au, v) + bv, λ) = F v) v H 1 Ω), bu, µ) = Gµ) µ H 1/2 2.4) Γ), 2

3 whee F v) = The following esult is well known. Ω fv dx and Gµ) = Γ gµ ds. Poposition 2.1. Thee exists a unique solution u, λ) H 1 Ω) H 1/2 Γ) to the poblem 2.4). Moeove, λ = u/ n whee n is the outwad nomal vecto on Γ. 3 Discetization with RBFs and finite elements We fist define the finite-dimensional space that appoximates u in 2.4). Let Φ : R d R be a adial basis function whose Fouie tansfom Φ satisfies Φω) 1 + ω 2 ) τ, ω R d, 3.1) whee τ > d/2. We conside the scaled adial basis functions Φ x) := d Φx/), > 0, x R d, so that Φ ω) ω 2 ) τ, ω R d. 3.2) The native space associated with Φ is defined by N Φ = { v D R d ) : Rd vω) 2 Φ ω) dω < } whee D R d ) is the space of distibutions defined in R d. This space is equipped with an inne poduct and a nom defined by vω)ŵω) v, w Φ = dω and v Φ := R d Φ ω) Rd vω) 2 Φ ω) dω ) 1/2. Unde the assumption 3.1), the native space N Φ is isomophic to the Sobolev space H τ R d ) with equivalent nom v ) ) τ/2 L2 R d ). Given a set of quasi-unifom centes X = {x 1,..., x N } Ω with mesh nom h X = sup x Ω min 1 j N x x j 2, we define H X, = span{φ 1,..., Φ N }, 3.3) whee Φ i x) = Φ x x i ) fo i = 1,..., N. Note that since the nodes can be nea to o even on the bounday Γ of Ω, the suppots of the scaled adial basis functions ae not necessaily subsets of Ω.) The solution u to 2.4) will be appoximated by u X H X,. Fo the appoximation of the Lagangian multiplie λ we use functions not necessaily continuous) which ae piecewise polynomials of degee p 0 defined on a quasi-unifom patition T k of the bounday Γ of Ω: Λ k := {µ : Γ R µ T P p T T k }. 3.4) 3

4 Hee, k is the mesh size of T k, and P p is the space of polynomials of degee at most p. Using these discete spaces, the Galekin scheme with adial basis functions and Lagangian multiplies fo the appoximate solution of 2.4) is: find u X, λ k ) H X, Λ k satisfying au X, v) + bv, λ k ) = F v) v H X,, 3.5) bu X, µ) = Gµ) µ Λ k. 4 Appoximation popety of scaled RBFs Fo any intege m 0 and eal > 0 we denote the nom of the scaled Sobolev space H m Ω) by v H m := 1/2 2 α D α v 2 L 2 Ω) α m with multi-index α = α 1,..., α d ) and α = α α d. We note that H 0 = H 0 Ω). Fo s 0, m) with intege m > 0 we define the scaled intepolation space H s := [H 0 Ω), H m ] θ, whee θ = s/m. Hee we employ the so-called eal K-method; see [2]. The intepolation nom in H s can be epesented as follows. Let S : DomS ) H m H 0 Ω) be an unbounded linea opeato defined by S v, w H 0 Ω) = v, w H m v DomS ), w H m. It is clea that S is self-adjoint and positive. Thus thee exists Λ := S 1/2 H 0 Ω) satisfying v, w H m = Λ v, Λ w H 0 Ω). : H m The inne poduct and nom in H s Ω) can now be epesented as v, w H s Λ = θ v, Λ θ w and v H s = Λ θ v H 0 Ω). 4.1) H 0 Ω) The Sobolev spaces H s fom a Hilbet scale with the following popety. Lemma 4.1. Let s 1 and s 2 be non-negative eal numbes, and let s 0 = s 1 + s 2 )/2. Then fo any v H s 1 H s 0 thee holds v H s 1 = sup w DΩ)\{0} v, w s H 0 w s H 2 whee DΩ) is the space of all functions which ae estictions on Ω of infinitely diffeentiable functions in R d. 4

5 Poof. Let m be an intege not less than max{s 1, s 2 }. We may assume that max{s 1, s 2 } > 0. Then H s i = [H 0 Ω), H m ] θi, whee θ i = s i /m, i = 0, 1, 2. Fo any v H s 1 H s 0, thee holds v H s 1 = sup z DΩ)\{0} v, z H s 1 Ω) z s H 1 Fo each z DΩ) we define w = Λ θ 1 θ 2 z. Then by noting 4.1), the elation θ 0 = θ 1 + θ 2 )/2, and the self-adjointness of Λ θ, we obtain z s H 1 = Λθ 1 z H 0 Ω) = Λ θ 2 w H 0 Ω) = w s H 2. Ω) and v, z s H 1 = Λ θ 1 v, Λ θ 1 z H0Ω) = = Λ θ 0 v, Λ θ 0 w Λ θ 1 v, Λ θ 2 w H 0 Ω) = v, w H s 0. H 0 Ω) Theefoe, v H s 1 = sup w DΩ)\{0} v, w s H 0 w s H 2 and the lemma is poved. Fo any v H τ, let I X v denote its intepolant in the space H X,, i.e., I X v H X, satisfies I X vx j ) = vx j ), j = 1,..., N. Fo a domain in two dimensions d = 2), the following appoximation popety is poved in [7, Lemma 5.3]; see also the poof of this lemma. The same aguments apply also to the case d = 3. Lemma 4.2. Let assumption 3.1) be satisfied. Then fo any v H τ thee holds v I X v H s Cs, τ) hx ) τ s v I X v H τ Cs, τ) fo 0 s τ and 0 < 0 with 0 > 0 abitay but fixed. hx ) τ s v H τ When v is smoothe, the eo bound can be extended by using the technique developed in [8], [12], and modified in [9] fo a sphee. Lemma 4.3. Assume that 3.1) holds. Let T be an opeato defined by T ψx) := Φ x y)ψy) dy, x R d. R d Then fo any s R 5

6 i) T is an isomophism fom H s τ R d ) onto H s+τ R d ) and satisfies T ψ H s+τ R d ) ψ H s τ R d ) ; ii) Fo any ψ H τ R d ) and ξ N Φ thee holds T ψ, ξ Φ = ψ, ξ H 0 R d ), i.e., T is the adjoint of the embedding opeato of the native space N Φ into H 0 R d ). Poof. The lemma follows fom T ψω) = Φ ω) ψω) ω 2 ) τ ψω), ω R d. With the help of the above lemma, we now extend the eo bound in Lemma 4.2 when v is smoothe. Lemma 4.4. Let the assumption 3.1) be satisfied. Then fo any s, t R satisfying 0 s τ t 2τ, if v H t Ω) then the following estimate holds v I X v H s Cs, t, τ) fo 0 < 0 with 0 > 0 abitay but fixed. hx Poof. Conside fist the case when t = 2τ. Let E : H σ H σ R d ) ) t s v H t be an -unifomly bounded extension opeato fo any σ > 0; cf. [7, Lemma 5.1]. Fo any v H 2τ Ω), since I X Ev = I X v = EI X v on Ω, one finds that Ev I X Ev is an extension of v I X v. Theefoe, the popety that I X is an othogonal pojection in N Φ yields v I X v 2 H τ Ev I XEv 2 H τ R d ) Ev I XEv 2 Φ Ev I X Ev, Ev Φ. It follows fom Lemma 4.3 that thee exists ψ H 0 R d ) such that and Theefoe T ψ = Ev, ψ H 0 R d ) Ev H 2τ R d ) Ev I X Ev, Ev Φ = Ev I X Ev, ψ H 0 R d ). v I X v 2 H τ Ev I XEv H 0 R d ) Ev H 2τ R d ) v I X v H 0 v H 2τ. 4.2) 6

7 By applying Lemma 4.2 with s = 0 and using the above inequality we obtain v I X v 2 H 0 Cτ) hx ) 2τ v I X v 2 H τ Cτ) hx ) 2τ v I X v H 0 v H 2τ. Thus the equied estimate is poved fo s = 0 and t = 2τ. It also follows fom 4.2) that the equied estimate holds fo s = τ and t = 2τ. By using intepolation we deduce the estimate fo 0 s τ and t = 2τ. Intepolation between ) 2τ s hx v I X v H s v H 2τ and v I X v H s hx ) τ s v H τ yields the equied estimate fo 0 s τ t 2τ, poving the lemma. To extend Lemma 4.4 to include less smooth functions, i.e. v H t Ω) fo 0 t < τ, we use Lemma 4.1. Lemma 4.5. Assume that 3.1) holds. Fo any v HΩ) t with 0 t 2τ thee exists z X H X, satisfying ) t s hx v z X H s Cs, t, τ) v H t 4.3) Ω) fo 0 s min{t, τ} and 0 < 0 with 0 > 0 abitay but fixed. Poof. We only need to pove the lemma fo 0 t < τ. Conside fist the case when τ/2 t < τ. Then 0 2t τ < t. Let P t : H t Ω) H X, be the pojection defined by It can be seen that P t v, z H t = v, z H t z H X,. 4.4) P t v v H t v H t. 4.5) If s [0, 2t τ] then τ 2t s < 2τ. Hence, fo any w H 2t s it follows fom Lemma 4.4 that ) t s hx w I X w H t w H 2t s. 4.6) By using Lemma 4.1 and 4.4) 4.6) we obtain P t v v H s = sup w DΩ)\{0} = sup w DΩ)\{0} hx P t v v, w H t w H 2t s P t v v, w I X w H t ) t s v H t. w H 2t s 7

8 In paticula, thee holds P t v v H 2t τ hx ) τ t v H t Ω). Hence, fo s [2t τ, t] by noting 4.5) and using intepolation we obtain the equied estimate, and thus pove 4.3) fo 0 s t and τ/2 t < τ. By successively consideing the case τ/4 t < τ/2, then τ/8 t < τ/4, etc., and using the same agument, we finish the poof of the lemma. We ae now able to deive the appoximation popety of H X, in non-scaled Sobolev noms. Lemma 4.6. Assume that 3.1) holds. Fo any v H t Ω) with 0 t 2τ thee exists z X H X, satisfying v z X H s Ω) Cs, t, τ) ht s X t v H t Ω) fo 0 s min{t, τ} and 0 < 0 with 0 > 0 abitay but fixed. Poof. Fist we note that since 0, 0 ], fo any function v and any positive intege m thee hold v H 0 Ω) = v H 0 and v H m v H m Ω) m v H m. By intepolation it follows that fo 0 s m v H s v H s Ω) s v H s. The equied esult is then a consequence of the above inequalities and Lemma Eo estimate In this section we pove ou main esult, Theoem 5.3, which establishes quasi-optimal convegence of ou mixed method and convegence odes. Of couse, appoximation popeties depend on the egulaity of solutions. To keep things simple, we assume that, fo smooth data, we have standad elliptic egulaity limited by the smoothness of Γ. Moe pecisely, let δ be such that the solution u of 2.1) fo any κ 0 and any sufficiently smooth data f, g, satisfies δ 1, 2] : u H δ Ω). 5.1) We estict ou consideations to egulaity no moe than H 2 Ω) since we ae inteested in non-smooth poblems and to simplify esults when appoximation spaces use piecewise polynomials of highe degee. Of couse, in two dimensions Ω being a polygon, thee holds { 2 if Ω is convex, δ = 2π/ω if Ω is non-convex, with ω π, 2π) being the angle of the lagest e-entant cone of the bounday Γ in case of non-convex Ω. The chaacteization of δ fo a polyhedal domain is a bit moe involved and not given hee. 8

9 We use standad Babuška-Bezzi theoy to pove the main esult. In ode to do so we now pove ellipticity of the bilinea fom a, ) on V X, := {v H X, ; bv, µ) = 0 µ Λ k } when κ = 0 and an inf-sup condition fo the bilinea fom b, ). Lemma 5.1. Fo k sufficiently small thee holds v H 1 Ω) v H 1 Ω) v V X,. Poof. The poof is standad cf. [5, Lemma 4.5]) and is given fo convenience of the eade. Let v V X, be given and decomposed as v = v 0 + D with D = Ω 1 Ω v dx so that Ω v 0 dx = 0. Hee, Ω denotes the measue of Ω. Thee holds v 2 H 1 Ω) = v 0 2 H 1 Ω) + Ω D2. 5.2) Now let µ L 2 Γ) \ {0} and its L 2 Γ)-pojection Πµ onto Λ k be given. By duality, a standad appoximation esult and the tace theoem, we find that thee holds that is, bd, µ) = bv v 0, µ) = bv, µ Πµ) bv 0, µ) D v H 1/2 Γ) k1/2 µ L2 Γ) + v 0 L2 Γ) µ L2 Γ) ) k 1/2 v H 1 Ω) + v 0 H 1 Ω) µ L2 Γ), bd, µ) sup k 1/2 v µ L 2 Γ)\{0} µ H 1 Ω) + v 0 H 1 Ω). 5.3) L2 Γ) Now, using 5.2), 5.3) and Poincaé-Fiedichs inequality, we find that v 2 H 1 Ω) v 0 2 H 1 Ω) + k v 2 H 1 Ω) v 2 H 1 Ω) + k v 2 H 1 Ω). Selecting k small enough finishes the poof. Lemma 5.2. Assume that 3.1) and the elliptic egulaity 5.1) hold. Suppose that the paametes h X, k and ae chosen such that K t h X, k, ) := ht 1 X k t 1 t is sufficiently small 5.4) with t = δ when δ < 2, o with some t [1, 2) when δ = 2. Then thee exists a positive constant α, independent of h X, k and, except fo thei elation via K t, such that thee holds α µ k H 1/2 Γ) sup v X H X, \{0} bv X, µ k ) v X H 1 Ω) µ k Λ k. 9

10 Poof. Fo any µ k Λ k, conside the poblem w + w = 0 in Ω, w n = µ k on Γ. 5.5) Since µ k H s Γ) fo any s < 1/2, thee exists a unique vaiational solution w H 1 Ω) of 5.5) with egulaity estimate limited by assumption 5.1) depending on the geomety) w H min{δ,3/2+s} Ω) µ k H min{δ 3/2,s} Γ) s [ 1/2, 1/2); 5.6) see e.g. [6]. Moeove, it is shown in [1, Theoem 2.7] that bw, µ k ) µ k 2 H 1/2 Γ). 5.7) On the othe hand, since 1 < δ 2 < 2τ we can invoke Lemma 4.6 to obtain, fo some z X H X,, { w z X H 1 Ω) ct) ht 1 X t = δ if δ < 2, t w H t Ω) with t [1, 2) if δ = 2. By using 5.6), the invese popety and the assumption 5.4), we deduce that µ k H t 3/2 Γ) k1 t µ k H 1/2 Γ) t [1, 2), w z X H 1 Ω) ht 1 X t k t 1 µ k H 1/2 Γ) = K th X, k, ) µ k H 1/2 Γ) 5.8) with t as given if δ = 2) o chosen if δ < 2) fo 5.4). This inequality, assumption 5.4) and 5.6) with s = 1/2 give z X H 1 Ω) µ k H 1/2 Γ). 5.9) On the othe hand 2.3), 5.7) and 5.8) yield fo K t h X, k, ) small enough i.e., This togethe with 5.9) yields bz X, µ k ) = bw, µ k ) + bz X w, µ k ) µ k 2 H 1/2 Γ) µ k H 1/2 Γ) w z X H 1 Ω) 1 K t h X, k, )) µ k 2 H 1/2 Γ) µ k 2 H 1/2 Γ), µ k H 1/2 Γ) bz X, µ k ). µ k H 1/2 Γ) µ k H 1/2 Γ) bz X, µ k ) z X H 1 Ω) sup v X H X, \{0} bv X, µ k ), v X H 1 Ω) poving the lemma. The following theoem is ou main esult. It poves the quasi-optimal convegence of ou mixed RBF appoximation of the Diichlet poblem with finite element Lagangian multiplie. 10

11 Theoem 5.3. Let us assume that 3.1) and 5.4) hold and conside adius paametes > 0 which ae bounded. In the case that κ = 0 in 2.1), we additionally assume that k is small enough. Then thee exists a unique solution u X, λ k ) H X, Λ k to the poblem 3.5). Moeove, let f and g be sufficiently smooth so that u, λ) is the solution to 2.4) with u H δ Ω) and λ H δ 3/2 Γ) with δ 1, 2], cf. 5.1). Then u u X H 1 Ω) + λ λ k H 1/2 Γ) inf u v H v H 1 Ω) + inf λ µ k H X, µ k Λ 1/2 Γ) k hδ 1 X δ u H δ Ω) + k δ 1 λ H δ 3/2 Γ), whee the implicitly appeaing constants depend on the constant C in 2.3), α in Lemma 5.2, and the ellipticity constant of bilinea fom a, ). Poof. Fo κ > 0, the fist bound of quasi-optimal convegence follows fom standad Babuška-Bezzi theoy cf. [3, Coollay ]) by making use of the continuity and ellipticity of the bilinea fom a, ), and continuity and continuous as well as discete inf-sup condition of the bilinea fom b, ), cf. Lemma 5.2. In the case that κ = 0 we use Lemma 5.1 to obtain V X, -ellipticity of av, v) = Ω v 2 dx. Then the esult follows the same way. To show the second estimate fo κ 0 we use Lemma 4.6 and the appoximation popety of Λ k. Depending on the egulaity of the solution, the eo estimate above can be simplified. Coollay 5.4. Let the assumptions of Theoem 5.3 be satisfied. If δ > 3/2 then h δ 1 u u X H 1 Ω) + λ λ k H 1/2 Γ) X δ + k δ 1) u H δ Ω). Poof. In the case u H δ Ω) with δ > 3/2, the nomal deivative λ = u/ n can be defined in the standad way nomal component tace of weak gadient) so that λ H δ 3/2 Γ) u H δ Ω). The assetion then is diect consequence of Theoem Numeical esults We conside the model poblem 2.1) with Ω = 0, 1) 0, 1), κ = 0 and f, g such that ux 1, x 2 ) = x x2 2, i.e. u Hδ Ω) with δ = 2. The nodes of X ae distibuted unifomly on Ω including nodes on the bounday. We use scaled adial basis functions with the adial basis functions defined in [10]. We conside the two cases τ = 1.5 and τ = 2.5 which coespond to C 0 and C 2 -functions, espectively. They ae otations of univaiate polynomials of degees 2 and 5, espectively. We have implemented the method by numeical integation with an ovekill of numbe of integation nodes. With δ = 2 assumption 5.4) equies that the atio be small enough fo k 1 ɛ 2 ɛ some ɛ > 0. We simply choose k h X, / moe pecisely an intege appoximation to 11 h 1 ɛ X

12 1/k smalle than o equal to 1 since the length of the sides of Ω is 1). In this way, fo fixed, K t is fixed and 5.4) is not guaanteed. Howeve, ou numeical esults do not show stability poblems that might be caused by a violation of the inf-sup condition) in the ange of unknowns unde consideation. Fo fixed, the eo estimate deived in Coollay 5.4 gives an uppe bound u u X H 1 Ω) + λ λ k H 1/2 Γ) h X + k. In the gaphs below we plot the individual eos u u X H 1 Ω) and λ λ k L2 Γ) on a double logaithmic scale. Fo the latte eo, which is measued in the L 2 athe than H 1/2 -nom, we expect a educed convegence like k 1/2. Both expected eo tems, h X labeled as h) and k 1/2, ae also given in the plots multiplied by 10 to shift them close to the coesponding eo cuves). Figues 1 and 2 show the esults fo = 0.2 with τ = 1.5 and τ = 2.5, espectively. Figues 3 and 4 show the coesponding esults fo educed adius = 0.1. In all the cases thee is some pe-asymptotic ange and the eos behave as expected fo lage numbe of unknowns h 10 k^{1/2} H1-eo L2-eLag) dimh_{x,}) Figue 1: Eos fo = 0.2 and τ = 1.5. Refeences [1] I. Babuška. The finite element method with Lagangian multiplies. Nume. Math., /73), [2] J. Begh and J. Löfstöm. Intepolation Spaces: An Intoduction. Spinge-Velag, Belin, [3] S. C. Benne and L. R. Scott. The Mathematical Theoy of Finite Element Methods. Spinge, Belin,

13 h 10 k^{1/2} H1-eo L2-eLag) dimh_{x,}) Figue 2: Eos fo = 0.2 and τ = 2.5. [4] Y. Duan and Y.-J. Tan. A meshless Galekin method fo Diichlet poblems using adial basis functions. J. Comput. Appl. Math., ), [5] G. N. Gatica, M. Healey, and N. Heue. The bounday element method with Lagangian multiplies. Nume. Methods Patial Diffeential Equations, ), [6] P. Gisvad. Elliptic Poblems in Nonsmooth Domains. Pitman, Boston, [7] N. Heue and T. Tan. Radial basis functions fo the solution of hypesingula opeatos on open sufaces. Computes and Mathematics with Applications, ), [8] R. Schaback. Impoved eo bounds fo scatteed data intepolation by adial basis functions. Math. Comp., ), [9] T. Tan, Q. T. Le Gia, I. H. Sloan, and E. P. Stephan. Bounday integal equations on the sphee with adial basis functions: Eo analysis. Appl. Nume. Math., ), [10] H. Wendland. Eo estimates fo intepolation by compactly suppoted adial basis functions of minimal degee. J. Appox. Theoy, ), [11] H. Wendland. Meshless Galekin methods using adial basis functions. Math. Comp., ), [12] H. Wendland. Scatteed Data Appoximation. Cambidge Univesity Pess, Cambidge,

14 h 10 k^{1/2} H1-eo L2-eLag) dimh_{x,}) Figue 3: Eos fo = 0.1 and τ = h 10 k^{1/2} H1-eo L2-eLag) dimh_{x,}) Figue 4: Eos fo = 0.1 and τ =

Measure Estimates of Nodal Sets of Polyharmonic Functions

Measure Estimates of Nodal Sets of Polyharmonic Functions Chin. Ann. Math. Se. B 39(5), 08, 97 93 DOI: 0.007/s40-08-004-6 Chinese Annals of Mathematics, Seies B c The Editoial Office of CAM and Spinge-Velag Belin Heidelbeg 08 Measue Estimates of Nodal Sets of