University of Cyprus. Abstract. singularities is analyzed. In this method, the solution is approximated by the
|
|
- Dale Caldwell
- 6 years ago
- Views:
Transcription
1 A Singular Function Boundary Integral Method for Elliptic Problems with Singularities hristos Xenophontos, Miltiades Elliotis and Georgios Georgiou Department of Mathematics and Statistics University of yprus P.O. Box Nicosia, yprus January 6, 005 Abstract A singular function boundary integral method for elliptic problems with boundary singularities is analyzed. In this method, the solution is approximated by the truncated asymptotic expansion for the solution near the singular point and the Dirichlet boundary conditions are weakly enforced by means of Lagrange multiplier functions. The resulting discrete problem is posed and solved on the boundary of the domain, away from the point of singularity. The main result of this paper is the proof of convergence of the method in particular, we show that the method approximates the generalized stress intensity factors, i.e. asymptotic expansion, at an exponential rate. convergence of the method are also presented. the coecients in the Numerical results illustrating the Key words: Boundary singularities boundary approximation methods Lagrange multipliers stress intensity factors. AMS Subject lassications: 65N5, 65N38, 65N, 65N30. Introduction In many engineering problems (e.g., fracture mechanics applications), governed by elliptic Partial Dierential Equations (PDEs), boundary singularities arise when there is a sudden orresponding author. xenophontos@ucy.ac.cy.
2 change in the boundary conditions and/or on the boundary itself. Singularities are known to aect adversely the accuracy and the convergence of standard numerical methods, such as nite element, boundary element, nite dierence and spectral methods. Grid renement and high-order discretizations are common strategies aimed at improving the convergence rate and accuracy of the above mentioned standard methods. If, however, the form of the singularity istaken into account and is properly incorporated into the numerical scheme, then a more eective method may be constructed (see, e.g., [7, 6]). (For a recent survey of treatment of singularities in elliptic boundary value problems see [8] and the references therein.) In the case of the two-dimensional Laplace equation, for example, the form of the singularity is visible through the asymptotic expansion for the solution u near the singular point. In polar co-ordinates (r ) centered at the singular point, u is given by [3] u(r )= X j= j r j j () () where j R are the unknown singular coecients, and j R j () are the eigenvalues and eigenfunctions of the problem, respectively, which are uniquely determined by the geometry and the boundary conditions along the boundaries sharing the singular point. The constants j are often called Generalized Stress Intensity Factors (GSIFs) or Flux IntensityFactors. In the case of elasticity problems with singularities, these constants are called Stress IntensityFactors (SIFs), and serve as a measure of the stress at which failure occurs. In most commonly used methods, such as the Finite Element Method (FEM), the SIFs are calculated as a post-solution operation (see, e.g., [4, 3, 4]). If, however, the goal of the computation is the calculation of the SIFs, then methods which calculate these quantities directly may be preferable. The Singular Function Boundary Integral Method (SFBIM), which we will analyze in this article, falls in the latter category. The SFBIM was originally developed by Georgiou et al. [7], and was subsequently rened and expanded by Elliotis et al. [8]{[0]. The method uses the leading terms in the local asymptotic expansion for the solution near the singular point as an approximation, while any Dirichlet boundary conditions are weakly enforced by means of Lagrange multiplier functions. The resulting problem is posed on the boundary of the domain, hence the dimension of the problem is reduced by one, leading to considerable computational cost reduction. We should also mention here the works of Li et al. [5]{[7], and Arad et al. [], who also developed similar methods. (See also [0] for a review of SIF evaluation and modeling of singularities in boundary integral methods.) The SFBIM has been successfully applied to a number of problems in solid and uid mechanics and excellent numerical results have been obtain thus far [8]{[]. In particular, it was observed that the method: (i) approximates the SIFs at an exponential rate of convergence, (ii) is very ecient, and (iii) compares extremely well with other accurate methods found in the literature. Our main goal in this article is to prove the observed convergence rates of the method.
3 The rest of the paper is organized as follows: In Section we present the formulation of the method for a two-dimensional Laplacian problem with a boundary singularity. In Section 3 we present the convergence analysis, and in Section 4 we comment onhow the method can be eciently implemented. Finally, Section 5 includes the results of some numerical computations illustrating the convergence of the method. Throughout the paper the usual notation H k () will be used for spaces containing functions on the domain R with having k generalized derivatives in L () the norm and seminorm on H k () will be denoted by kk k and jj k, respectively. Also, the letters and c will be used to denote generic positive constants independent ofany discretization parameters and possibly having dierent values in each occurrence. The Model Problem and its Formulation For simplicity we consider the Laplacian problem stated below and depicted graphically in Figure. Find u such that S = 4 S i. i= u = 0in = 0onS (3) u = 0onS (4) u = f(r )ons 3 = g(r )ons 4 (6) Figure : A two-dimensional Laplace equation problem with one boundary singularity. 3
4 In (){(6), denotes the Laplacian operator, the variables (r ) denote polar coordinates centered at O and the functions f and g are known. It is assumed that f g and the are such that there is only one boundary singularity ato. The local solution near O, isgiven by an asymptotic expansion of the form () [3]. Multiplying () by a test function v V (to be specied shortly) and integrating over, we get vu =0 and then, by means of Green's Theorem, we obtain ; rv ru =0: =0onS = g(r )ons 4,we further have Suppose v is chosen to satisfy rv ru ; S v = 0 in v=0ons ; S4 = vg: (7) one such choice being v r j j () (see eq. ()). Then (7) becomes Now, since u = f on S 3 wehave so adding this to (9), we get or, equivalently, Letting rv ru ; rv ru ; rv ru ; =0onS (8) S4 = vg: (u ; f) =0 ; S4 (u ; f) = S4 u = vg = (0) 4
5 the above equation becomes rv ru ; v ; u = S4 vg ; Hence, the variational problem to be solved reads: Find (u ) V V such that B(u v)+b(u v ) =F (v ) 8 (v ) V V () f: where B(u v)= b(u v ) =; F (v ) = S4 rv ru () v ; vg ; u (3) f: (4) The spaces V and V are chosen as follows. Let the trace space of functions in H () be denoted by H = (@) = u H () : L (@) : (5) With T : H ()! H = (@) denoting the trace operator, the norm of H = (@) is dened as n o k k = inf kuk uh : Tu= : (6) () Then, we dene H ;= (@) as the closure of H 0 (@) L (@) with respect to the norm k'k = ' : (7) k k sup H = (@) (See, e.g., [] for more details.) We then dene H () = u H () : uj S =0 (8) and we take V = H () V = H ;= (S 3 ) : (9) The discrete problem corresponding to (), then, reads: Find (u N h ) V N V h H H ;= (S 3 ) such that B(u N v)+b(u N v h )=F (v ) 8 (v ) V N V h (0) with B(u v) b(u v ) and F (v ) given by (){(4), and with V N V h nite dimensional spaces to be chosen shortly. 5
6 Remark The above formulation will be used in the analysis of the method. For the implementation, the following equivalent formulation will be used, in order to maximize computational eciency: From (9), using Green's formula once more, we get ; From (4) and (8) we obtain uv ; S4 = vg: S4 + ; S4 = vg () which is the equation that we discretize (see Section 4 ahead). Note that all integrals are one-dimensional, hence the computational cost is considerably reduced. Moreover, all integrations are carried out away from the point of singularity, where all integrands are smooth. 3 Error Analysis We begin by dening the nite dimensional spaces V N and V h which will be used in the approximate problem (0). First, with v i r i i () () denoting the singular functions, we dene the nite dimensional space V N = span fv i g N i= : (3) S n Next, let S 3 be divided into quasiuniform sections ; i i= ::: nsuch thats 3 = ; i= i. Let h i = j; i j and set h =max in h i. We assume that for each segment ; i, there exists an invertible mapping F i : I! ; i which maps the interval I =[; ] to ; i and we dene the space V h = h : h j ;i F ; i P p (I) i= ::: n (4) where P p (I) is the set of polynomials of degree p on I =[; ]: In practice, the representation of the boundary S 3 determines the mappings F i i.e.,ifs 3 is represented by a polynomial then an isoparametric mapping may be used, and if S 3 has some (general) parametric representation, then the blending map technique may be used (see h. 6 in []). We have the following theorem. 6
7 Theorem Let (u ) and (u N h ) be the solutions to () and (0), respectively. Suppose there exist positive constants c 0 c and, independent of N and h, suchthatthe following hold: Then B(v v) c 0 kvk and jb(u v)jc kuk kvk 8v V N (5) 9 0 6= v N V N s.t. ku ; u N k + k ; h k ;= with R + independent of N and h. h v N k h k ;= kv N k v kk ;= kvk inf vv N 8 h V h (6) 8v V N : (7) ku ; vk + inf V h k ; k ;= (8) This is similar to the classical result for saddle point problems (cf. [5]), and its proof for our problem appears in the appendix see also Theorem 6. in [5]. Before verifying that (5){(7) hold for the problem under consideration, we will make certain assumptions that will aid in the analysis. First, we note that for any function w which can be expressed in the form given by equation (), we can always write where w N = NX j= w = w N + r N j v j V N r N = We will assume that there exists a (0 ) such that jr N ja N NaN : X j=n + j v j : (9) (A) (A) If r< in (), assumptions (A), (A) hold trivially, since then, by () and the fact that j is harmonic, jr N j X j=n + j j j r j r N+ ; r an 7
8 with r<a<, and R + independent ofa and N: (The boundedness of the GSIFs was also used above.) X j=n + d dr j j j r j; = X j=n + j j j r j X j=n + d dr 8 < j j j : d dr r N+ ; r r 0 9 = j; d = d dr Na N : X j=n + j j j 8 < : r 0 j; d If r, one may partition the domain into subdomains in which separate approximations may be obtained, including one near O which isvalid for r <. The solution over the entire domain can then be composed by combining the solutions from each subdomain and properly dealing with their interactions across the interfaces separating them (see, e.g., [9]). Let us now verify that (5){(7) hold for the problem given by equation (0). First, note that B(v v)=jvj so that, by Poincare's inequality, By the auchy-schwartz inequality, B(v v) c 0 kvk 8v H () : (30) B(u v) c kuk kvk 8u v H () (3) so that (30) and (3) give (5). To verify (6), consider the following auxiliary problem (for which a unique solution exists): Find w such that where h Poincare's inequality) 9 = w =0in = h on S 3 (33) w =0onS =0onS [ S 4 (35) V h in (33). From (3) and (33) we obtain (using Green's formula and h w = with c 0 R + : Also, cf. [], = ww + jrwj = jwj c 0 kwk (36) k h k ;= kwk (37) ;= 8
9 so that by (36) and (37) with R + independent ofw and h. as given by equation (9). We have and also so that combining (38){(40) we get h w c 0 kwk kwk k hk ;= (38) h w N = Now, let w N V N be such thatw = w N + r N, h w ; h r N (39) h r N k h k ;= kr N k = k h k ;= kr N k (40) h w N kwk k h k ;= ; k h k ;= kr N k : (4) Now, using the reverse triangle inequality, we have kwk = kw N + r N k kw N k ;kr N k which along with (4), gives h w N kw N k ;kr N k k h k ;= ; k h k ;= kr N k kw N k k h k ;= ; ( + ) k h k ;= kr N k : (4) Since by assumption (A), r N converges to 0 exponentially, wehave and for N suciently large we may write 0 < kr Nk kw N k < kr N k kw N k ( + ) (43) where and are the constants from above. ombining (4) and (43) leads to h w N k hk ;= kw N k which in turn gives (6) once we replace w N by v N and = by. ondition (7) follows directly from the denition of the H ;= -norm (see also (40)). The preceding discussion leads to the following theorem. 9
10 Theorem 3 Let (u ) and (u N h ) be the solutions to () and (0), respectively. If H k (S 3 ) for some k, then there exists a positive constant, independent of N h and a (0 ), such that where m = minfk p +g: ku ; u N k + k ; h k ;= Proof. From Theorem we have ku ; u N k + k ; h k ;= inf vv N n pna N + h m p ;k o (44) ku ; vk + inf V h Now, inf ku ; vk ku ; w N k = kr N k vv N with w N r N given by (9). Using (A) and (A) we get inf vv N k ; k ;= : ku ; vk kr N k 0 + jr N j na N + p o Na N p Na N (45) with R + independent ofn. Next, let I be the p th ;order interpolant of on the partition f; i g n i= of S 3.Wehave (cf. [3]) k ; I k 0 h m p ;k kk k (46) where m =minfk p +g and >0 a constant independent ofh and p. Now, inf V h k ; k ;= k ; I k ;= k ; I k 0 (47) so that, since H k (S 3 ), (45){(47) give the desired result. Remark 4 The above theorem shows that if the number of singular functions N is increased then u N converges to u at an exponential rate. The theorem also shows that the convergence of h to can occur in one of three ways: (i) by keeping p xed andreducing h, (ii) by keeping h xed and increasing p, or (iii) by doing both. These three \options" loosely correspond to the three versions of the FEM, namely the h, p and hp versions (cf. [3]). Remark 5 Based on the above theorem, one may obtain the \optimal" matching between N and h, i.e. the relationship between the number of singular functions and the number of Lagrange multipliers used in the method, by choosing them in such a way so that the error in (44) is balanced. For example, in the case when p is kept xed and h! 0, we take h p+ p Na N. This leads to the following approximate expression for N: N (p +) ln h ln a : (48) 0
11 The approximation of the GSIFs is given by the following orollary. orollary 6 Let and u = u N = X j= NX j= j r j j () N j r j j () satisfy () and (0), respectively, with j N j denoting the true and approximate GSIFs. Then, there exists a positive constant, independent of N and a (0 ), such that j ; N j Proof. This is a direct consequence of (A) and the fact that (See also eq. (45).) 4 Implementation j ; N j a N : (49) ku ; un k 0 : We now give a brief description of the implementation in order to emphasize the properties of the method. As described in Remark, the equation which we discretize is We seek such that S4 + ; S4 = vg: u N = NX j= N j v j i u N S4 i u N N v i S4 = v i g (5) with v i given by (). The normal derivative ofu N above is replaced by the Lagrange multiplier function (cf. equations (0) and (4)) h = MX k= k k V h (S 3 ) (5)
12 where k R (to be determined), and V h (S 3 )=span f k g M. The boundary condition k= on S 3 is weakly enforced by requiring that i (f ; u N )=0 i= ::: M: (53) ombining (5){(53), we obtain the (N + M) (N + M) linear system K K ;!;! " ;!G;! # K T = (54) 0 F where ;! = N ::: N N T ;! =[ ::: M ] T and [K ] j k = [K ] j i = ; h ;!G i h ;!F i j i = = j v k S4 + v j j= ::: N k = ::: N v j i j= ::: N i= ::: M v j g j= ::: N f i i= ::: M: It is easily shown that the coecient matrix in (54) is symmetric. This matrix, however, is singular if N < M. Hence, we shouldchoose N larger than M, but not too large since for excessively large values of N the linear system (54) becomes ill-conditioned and the results obtained are unreliable. The relationship between N and h (hence M) described in Remark 4 should be our guide in deciding how large to choose these values { see also Section 5 ahead. As a nal remark in this section, we should point out that all integrals involved in the determination of the coecient matrix (and right hand side) in (54) are along the sides of the domain that do not contain the singularity. Moreover, they are one-dimensional and can be approximated by standard techniques, such as Gaussian quadrature. 5 Numerical Results We consider the model problem depicted in Figure, originally studied in [4]. The local asymptotic solution may be written as u = X j= j v j
13 where Figure : Geometry and boundary conditions for the model problem. v j = r (4j;3)=3 sin 3 (4j ; 3) : The numerical results presented here correspond to the following choices of all relevant parameters: The Lagrange multiplier function h used to impose the Dirichlet condition along S 3 is expanded in terms of quadratic basis functions i (see equation (5)), and boundary S 3 is divided into n quadratic elements of equal size. For the integration, boundary S 4 is also subdivided into n intervals of equal size. All integrals involved are calculated numerically by subdividing each interval above into 0 subintervals and using a 5-point Gauss-Legendre quadrature over each one. In computing the coecient matrix in (54), its symmetry is taken into account. To determine the relationship between the number of singular functions N, and the number of Lagrange multipliers M, we proceed as follows: Since for h we are using p = and h ==n, wehave M =n +: For the moment, we xn = 8 (say), which amounts to M = 7, and solve the linear system (54) for various values of N>M(e.g., N =9 3 :::). We concentrate only on the calculation of the rst GSIF N and record our results in Table. From the table we see that, for this choice for M, thevalue of a N is converged up to 4 signicant digits once N =35. Moreover, from (48) we have N (p +) ln h ln a ln a p + N =(p +)ln(=n) ln a ln p M ; ) a ln; p M; =(p +) ln a p M ; p+ N hence, using M =7andN = 35, we nd that a 0:89. With a known, we maynow compute the rest of the GSIFs and/or any other quantities of interest by choosing N and 3 )
14 N a N 9 : : : : : : : : : Table : Approximate GSIF a N with M = 7. M via equation (48). For example, for M =4(i.e. h ==0), we nd that N 60. Table below shows the converged approximate values for the rst ve GSIFs obtained using the SFBIM (with M =4andN =60),aswell as the hp-version of the FEM, which is considered to be the state-of-the-art for problems with singularities [] the values from reference [4] are also included for comparison. These results suggest that j a N j (SFBIM) a N j (hp-fem) a N j (Ref. [4]) : :79800 :80 0: : :699 3 ;0: ;0:03049 ;0: : : : : : :0009 Table : Approximate GSIFs a N j, j = ::: 5. the SFBIM can be an attractive (and often preferable) method for problems in which the GSIFs are the main goal of the computation. Figure 3 shows the convergence of the leading GSIFs witn N { in particular, the gure shows a semilog plot of the percentage relative error in a N j j = ::: 5. Since each error curve is essentially straight, this illustrates the exponential convergence of the method, as predicted by orollary 6. References [] M. Arad,. Yosibash, G. Ben-Dor and A. Yakhot, \omparing the ux intensity factors by a boundary element method for elliptic equations with singularities," ommun. Numer. Meth. Eng., Vol. 4, pp. 657{670 (998). 4
15 GSIFs using the SFBIM j = j = j = 3 j = 4 j = Percentage relative error N Figure 3: onvergence of a N j j = ::: 5 using the SFBIM. [] I. Babuska, \The nite element method with Lagrangian multipliers", Numer. Math., Vol. 0, pp. 79{9 (973). [3] I. Babuska and M. Suri, \The p and h-p versions of the nite element method, basic principles and properties", SIAM Review, Vol. 36, No. 4, pp. 578{63 (994). [4] S.. Brenner, \Multigrid methods for the computation of singular solutions and stress intensity factors I: orner singularities", Math. omp. Vol. 68, pp. 559{583 (999). [5] F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Methods, Springer-Verlag (99). [6] P. G. iarlet, The Finite Element Method for Elliptic Problems, lassics in Applied Mathematics 40, SIAM (00). [7] L. G. Olson, G.. Georgiou and W. W. Schultz, \An ecient nite element method for treating singularities in Laplace's equation", J. omp. Physics Vol. 96, pp. 39{ 40 (99). [8] M. Elliotis, G. Georgiou and. Xenophontos, \The solution of Laplacian problems over L-shaped domains with a singular function boundary integral method", 5
16 ommunications in Numerical Methods in Engineering, Vol. 8, No. 3, pp. 3{ (00) [9] M. Elliotis, G. Georgiou and. Xenophontos, \Solving Laplacian problems with boundary singularities: A comparison of a singular function boundary integral method with the p/hp version of the nite element method", accepted for publication in Appl. Math. omp., 004. [0] M. Elliotis, G. Georgiou and. Xenophontos, \Solution of the planar Newtonian stick-slip problem with a singular function boundary integral method", submitted, 004. [] M. Elliotis, G. Georgiou and. Xenophontos, \The singular function boundary integral method for a two-dimensional fracture problem", submitted, 004. [] G. Georgiou, A. Boudouvis, A. Poullikkas, \omparison of two methods for the computation of singular solutions in elliptic problems", J. omput. Appl. Math., Vol. 79, pp. 77{90 (997). [3] P. Grisvard, Singularities in Boundary Value Problems, Springer-Verlag, Berlin, 99. [4] T. Igarashi and I. Honma, \A boundary element method for potential elds with corner singularities", Appl. Math. Modelling, Vol. 0, pp. 847{85 (996). [5].. Li, ombined Methods for Elliptic Equations with Singularities, Interfaces and Innities, Kluwer Academic Plublishers, Boston, Amsterdam (998). [6].. Li, \Penalty combinations of Ritz-Galerkin and nite-dierence methods for singularity problems", J. omput. Appl. Math., Vol. 8, pp. {7 (997). [7].. Li, Y. L. han, G. Georgiou and. Xenophontos, \Special boundary approximation methods for Laplace's equation with singularities", accepted for publication in omp. Math. Appl., 004. [8].. Li and T. T. Lu, \\Singularities and treatments of elliptic boundary value problems," Math. omput. Model., Vol. 3, pp. 97{45 (000). [9].. Li, R. Mathon and P. Sermer, \Boundary methods for solving elliptic problems with singularities and interfaces", SIAM J. Num. Anal. Vol. 4, No. 3, pp. 487{498 (987). [0] N. K. Mukhopadhyay, S. K. Maiti and A. Kakodkar, \A review of SIF evaluation and modelling of singularities in BEM," omp. Mech., Vol. 5, pp. 358{375 (000). 6
17 [] J. Pitkaranta, \Boundary subspaces for the nite element method with Lagrange multipliers", Numer. Math. Vol. 33, pp. 73{89 (973). [] B. A. Szabo and I. Babuska, Finite Element Analysis, Wiley, 99. [3] B. A. Szabo and. Yosibash,\Numerical analysis of singularities in two dimensions. Part : omputation of the eigenpairs", Int. J. Numer. Methods Eng. Vol. 38, pp. 055{08 (995). [4] B. A. Szabo and. Yosibash,\Numerical analysis of singularities in two dimensions. Part : omputation of generalized ux/stress intensity factors", Int. J. Numer. Methods Eng. Vol. 39, pp. 409{434 (996). A Proof of Theorem Here, for completeness, we give a proof of Theorem. First, we note that since (u ) also satisfy (0), we have 8 (v ) V N V h B(u ; u N v)=;b(u ; u N v ; h )= v( ; h )+ Since u = f on S 3 and R (u N ; f) =08 V h,wehave u N = and thus the last integral in (55) is zero. Hence, Letting w = u N ; v v V N,we get B(u ; u N v)= (u ; u N ): (55) u 8 V h (56) v( ; h ): (57) B(u N ; v w) = B(u N ; u w)+b(u ; v w) =B(u ; v w) ; = B(u ; v w) ; By means of (56), the above equation becomes B(u N ; v w) =B(u ; v w) ; w( ; ) ; 7 w( ; ) ; w( ; h ) w( ; h ) ( h ; )(u ; v): (58)
18 Now, from (5), (7) and (58) we get c 0 ku N ; uk = c 0 kwk ku ; vk kwk + kwk k ; k ;= + ku ; vk k h ; k ;= nku ; vk + k ; k ;= kwk + ku ; vk k h ; k ;= o : This is an inequality of order : where x = kwk c 0 x bx + d b d>0 b= d = ku ; vk k h ; k ;= Therefore, we obtain the bound ku ; vk + k ; k ;= " k h ; k ;= + " ku ; vk ">0: x b + p b +4c 0 d c 0 or, equivalently, ku ; vk + k h ; k ;= + " ku ; vk kwk + "k h ; k ;= : (59) Next, using (6) with h = h ; we nd that there exists 0 6= v N V N such that k h ; k ;= R h ; )v N ( kv N k : (60) Also, it follows from (57) that ( h ; )v N = ( h ; )v N + ( ; )v N = B(u ; u N v N )+ ( ; )v N c ku ; u N k kv N k + kv N k k ; k ;= (6) so that by (60) and (6), o k h ; k ;= nku ; u N k + k ; k ;= o nku ; vk + kv ; u N k + k ; k ;= : 8
19 Since kv ; u N k = kwk,wehave from (59) k h ; k ;= n ku ; vk + k ; k ;= + " ; ku ; vk o + " k h ; k ;= with R + independent ofv and ": Letting " = weget o k h ; k ;= nku ; vk + k ; k ;= : (6) Using the triangle inequality, we immediately get k ; h k ;= k ; k ;= + k h ; k ;= o nku ; vk + k ; k ;= : (63) Similarly, using (59) and (63), we get ku ; u N k ku ; vk + kv ; u N k = ku ; vk + kwk o nku ; vk + k ; k ;= which gives the desired result. 9
Applied Mathematics and Computation 169 (2005)
Applied Mathematics and Computation 169 (2005) 485 499 www.elsevier.com/locate/amc Solving Laplacian problems with boundary singularities: a comparison of a singular function boundary integral method with
More informationNumerical Solutions of Laplacian Problems over L-Shaped Domains and Calculations of the Generalized Stress Intensity Factors
WCCM V Fifth World Congress on Computational Mechanics July 7-2, 2002, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Numerical Solutions of Laplacian Problems over L-Shaped Domains
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationNumerical Integration for Multivariable. October Abstract. We consider the numerical integration of functions with point singularities over
Numerical Integration for Multivariable Functions with Point Singularities Yaun Yang and Kendall E. Atkinson y October 199 Abstract We consider the numerical integration of functions with point singularities
More informationThe Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 29, 1999 Abstract This paper deals with the construction of wavelet approximatio
The Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 9, 1999 Abstract This paper deals with the construction of wavelet approximation spaces, in the framework of the Mortar method. We
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationt x 0.25
Journal of ELECTRICAL ENGINEERING, VOL. 52, NO. /s, 2, 48{52 COMPARISON OF BROYDEN AND NEWTON METHODS FOR SOLVING NONLINEAR PARABOLIC EQUATIONS Ivan Cimrak If the time discretization of a nonlinear parabolic
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationwith Applications to Elasticity and Compressible Flow Daoqi Yang March 20, 1997 Abstract
Stabilized Schemes for Mixed Finite Element Methods with Applications to Elasticity and Compressible Flow Problems Daoqi Yang March 20, 1997 Abstract Stabilized iterative schemes for mixed nite element
More informationXIAO-CHUAN CAI AND MAKSYMILIAN DRYJA. strongly elliptic equations discretized by the nite element methods.
Contemporary Mathematics Volume 00, 0000 Domain Decomposition Methods for Monotone Nonlinear Elliptic Problems XIAO-CHUAN CAI AND MAKSYMILIAN DRYJA Abstract. In this paper, we study several overlapping
More informationFinite element approximation on quadrilateral meshes
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2001; 17:805 812 (DOI: 10.1002/cnm.450) Finite element approximation on quadrilateral meshes Douglas N. Arnold 1;, Daniele
More informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
More informationElement diameter free stability parameters. for stabilized methods applied to uids
Element diameter free stability parameters for stabilized methods applied to uids by Leopoldo P. Franca Laboratorio Nacional de Computac~ao Cientica (LNCC/CNPq) Rua Lauro Muller 455 22290 Rio de Janeiro,
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationquantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner
Quantitative Justication of Linearization in Nonlinear Hencky Material Problems 1 Weimin Han and Hong-ci Huang 3 Abstract. The classical linear elasticity theory is based on the assumption that the size
More information/00 $ $.25 per page
Contemporary Mathematics Volume 00, 0000 Domain Decomposition For Linear And Nonlinear Elliptic Problems Via Function Or Space Decomposition UE-CHENG TAI Abstract. In this article, we use a function decomposition
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationon! 0, 1 and 2 In the Zienkiewicz-Zhu SPR p 1 and p 2 are obtained by solving the locally discrete least-squares p
Analysis of a Class of Superconvergence Patch Recovery Techniques for Linear and Bilinear Finite Elements Bo Li Zhimin Zhang y Abstract Mathematical proofs are presented for the derivative superconvergence
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationREPORTRAPPORT. Centrum voor Wiskunde en Informatica. Asymptotics of zeros of incomplete gamma functions. N.M. Temme
Centrum voor Wiskunde en Informatica REPORTRAPPORT Asymptotics of zeros of incomplete gamma functions N.M. Temme Department of Analysis, Algebra and Geometry AM-R9402 994 Asymptotics of Zeros of Incomplete
More informationAn Optimal Finite Element Mesh for Elastostatic Structural. Analysis Problems. P.K. Jimack. University of Leeds. Leeds LS2 9JT, UK.
An Optimal Finite Element Mesh for Elastostatic Structural Analysis Problems P.K. Jimack School of Computer Studies University of Leeds Leeds LS2 9JT, UK Abstract This paper investigates the adaptive solution
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationStability of implicit extrapolation methods. Abstract. Multilevel methods are generally based on a splitting of the solution
Contemporary Mathematics Volume 00, 0000 Stability of implicit extrapolation methods Abstract. Multilevel methods are generally based on a splitting of the solution space associated with a nested sequence
More informationConvergence of A Galerkin Method for 2-D Discontinuous Euler Flows Jian-Guo Liu 1 Institute for Physical Science andtechnology and Department of Mathe
Convergence of A Galerkin Method for 2-D Discontinuous Euler Flows Jian-Guo Liu 1 Institute for Physical Science andtechnology and Department of Mathematics University of Maryland College Park, MD 2742
More informationFind (u,p;λ), with u 0 and λ R, such that u + p = λu in Ω, (2.1) div u = 0 in Ω, u = 0 on Γ.
A POSTERIORI ESTIMATES FOR THE STOKES EIGENVALUE PROBLEM CARLO LOVADINA, MIKKO LYLY, AND ROLF STENBERG Abstract. We consider the Stokes eigenvalue problem. For the eigenvalues we derive both upper and
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationThe partition of unity quadrature in meshless methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2002; 54:987 1006 (DOI: 10.1002/nme.455) The partition of unity quadrature in meshless methods A. Carpinteri, G. Ferro
More informationMARY ANN HORN be decoupled into three wave equations. Thus, we would hope that results, analogous to those available for the wave equation, would hold
Journal of Mathematical Systems, Estimation, and Control Vol. 8, No. 2, 1998, pp. 1{11 c 1998 Birkhauser-Boston Sharp Trace Regularity for the Solutions of the Equations of Dynamic Elasticity Mary Ann
More information1e N
Spectral schemes on triangular elements by Wilhelm Heinrichs and Birgit I. Loch Abstract The Poisson problem with homogeneous Dirichlet boundary conditions is considered on a triangle. The mapping between
More information84 G. Schmidt (1.3) T : '! u := uj ; where u is the solution of (1.1{2). It is well known (cf. [7], [1]) that the linear operator T : H ( )! e H ( ) i
Numer. Math. 69: 83{11 (1994) Numerische Mathematik c Springer-Verlag 1994 Electronic Edition Boundary element discretization of Poincare{Steklov operators Gunther Schmidt Institut fur Angewandte Analysis
More informationBoundary-Layer Transition. and. NASA Langley Research Center, Hampton, VA Abstract
Eect of Far-Field Boundary Conditions on Boundary-Layer Transition Fabio P. Bertolotti y Institut fur Stromungsmechanik, DLR, 37073 Gottingen, Germany and Ronald D. Joslin Fluid Mechanics and Acoustics
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationTaylor series based nite dierence approximations of higher-degree derivatives
Journal of Computational and Applied Mathematics 54 (3) 5 4 www.elsevier.com/locate/cam Taylor series based nite dierence approximations of higher-degree derivatives Ishtiaq Rasool Khan a;b;, Ryoji Ohba
More informationMixed discrete least square meshless method for solution of quadratic partial di erential equations
Scientia Iranica A (014) 1(3), 49{504 Sharif University of Technology Scientia Iranica Transactions A: Civil Engineering www.scientiairanica.com Mixed discrete least square meshless method for solution
More informationMULTIGRID PRECONDITIONING FOR THE BIHARMONIC DIRICHLET PROBLEM M. R. HANISCH
MULTIGRID PRECONDITIONING FOR THE BIHARMONIC DIRICHLET PROBLEM M. R. HANISCH Abstract. A multigrid preconditioning scheme for solving the Ciarlet-Raviart mixed method equations for the biharmonic Dirichlet
More informationAn Accurate Semi-Analytic Finite Difference Scheme for Two-Dimensional Elliptic Problems with Singularities
An Accurate Semi-Analytic Finite Difference Scheme for Two-Dimensional Elliptic Problems with Singularities Z. Yosibash, M. Arad, A. Yakhot, G. Ben -Dor Pearlstone Center for Aeronautical Engineering Studies
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationA GLOBAL COLLOCATION METHOD FOR TWO-DIMENSIONAL RECTANGULAR DOMAINS
JOURNAL OF MECHANICS OF MATERIALS AND STRUCTURES Vol., No., A GLOBAL COLLOCATION METHOD FOR TWO-DIMENSIONAL RECTANGULAR DOMAINS CHRISTOPHER G. PROVATIDIS This paper proposes the use of a global collocation
More informationContents Introduction 3 Johnson's paradigm 4 3 A posteriori error analysis for steady hyperbolic PDEs 7 3. Symmetric hyperbolic systems : : : : : : :
Report no. 96/09 Finite element methods for hyperbolic problems: a posteriori error analysis and adaptivity a Endre Suli Paul Houston This paper presents an overview of recent developments in the area
More informationEvan Mitsoulis. School of Mining Engineering and Metallurgy. National Technical University of Athens. R.R. Huilgol
Cessation of Couette and Poiseuille ows of a Bingham plastic and nite stopping times Maria Chatzimina, Georgios C. Georgiou Department of Mathematics and Statistics, University of Cyprus P.O. Box 2537,
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationMath 497 R1 Winter 2018 Navier-Stokes Regularity
Math 497 R Winter 208 Navier-Stokes Regularity Lecture : Sobolev Spaces and Newtonian Potentials Xinwei Yu Jan. 0, 208 Based on..2 of []. Some ne properties of Sobolev spaces, and basics of Newtonian potentials.
More informationLeast-squares collocation meshless method
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2001; 51:1089 1100 Least-squares collocation meshless method Xiong Zhang ;, Xiao-HuLiu, Kang-ZuSong and Ming-Wan Lu
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationMODIFIED STREAMLINE DIFFUSION SCHEMES FOR CONVECTION-DIFFUSION PROBLEMS. YIN-TZER SHIH AND HOWARD C. ELMAN y
MODIFIED STREAMLINE DIFFUSION SCHEMES FOR CONVECTION-DIFFUSION PROBLEMS YIN-TZER SHIH AND HOWARD C. ELMAN y Abstract. We consider the design of robust and accurate nite element approximation methods for
More information114 EUROPHYSICS LETTERS i) We put the question of the expansion over the set in connection with the Schrodinger operator under consideration (an accur
EUROPHYSICS LETTERS 15 April 1998 Europhys. Lett., 42 (2), pp. 113-117 (1998) Schrodinger operator in an overfull set A. N. Gorban and I. V. Karlin( ) Computing Center RAS - Krasnoyars 660036, Russia (received
More informationUniversity of Hertfordshire Department of Mathematics. Study on the Dual Reciprocity Boundary Element Method
University of Hertfordshire Department of Mathematics Study on the Dual Reciprocity Boundary Element Method Wattana Toutip Technical Report 3 July 999 Preface The boundary Element method (BEM) is now recognised
More informationCONVERGENCE OF GAUGE METHOD FOR INCOMPRESSIBLE FLOW CHENG WANG AND JIAN-GUO LIU
MATHEMATICS OF COMPUTATION Volume 69, Number 232, Pages 135{1407 S 0025-571(00)0124-5 Article electronically published on March 24, 2000 CONVERGENCE OF GAUGE METHOD FOR INCOMPRESSIBLE FLOW CHENG WANG AND
More informationXIAOFENG REN. equation by minimizing the corresponding functional near some approximate
MULTI-LAYER LOCAL MINIMUM SOLUTIONS OF THE BISTABLE EQUATION IN AN INFINITE TUBE XIAOFENG REN Abstract. We construct local minimum solutions of the semilinear bistable equation by minimizing the corresponding
More informationRening the submesh strategy in the two-level nite element method: application to the advection diusion equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 39:161 187 (DOI: 10.1002/d.219) Rening the submesh strategy in the two-level nite element method: application to
More informationA NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES EQUATIONS
MATHEMATICS OF COMPUTATION Volume 74, Number 250, Pages 543 554 S 0025-5718(04)01711-9 Article electronically published on August 31, 2004 A NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationAn explicit nite element method for convection-dominated compressible viscous Stokes system with inow boundary
Journal of Computational and Applied Mathematics 156 (2003) 319 343 www.elsevier.com/locate/cam An explicit nite element method for convection-dominated compressible viscous Stokes system with inow boundary
More informationPositive Denite Matrix. Ya Yan Lu 1. Department of Mathematics. City University of Hong Kong. Kowloon, Hong Kong. Abstract
Computing the Logarithm of a Symmetric Positive Denite Matrix Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A numerical method for computing the logarithm
More informationSPECTRAL METHODS ON ARBITRARY GRIDS. Mark H. Carpenter. Research Scientist. Aerodynamic and Acoustic Methods Branch. NASA Langley Research Center
SPECTRAL METHODS ON ARBITRARY GRIDS Mark H. Carpenter Research Scientist Aerodynamic and Acoustic Methods Branch NASA Langley Research Center Hampton, VA 368- David Gottlieb Division of Applied Mathematics
More information0 Finite Element Method of Environmental Problems. Chapter where c is a constant dependent only on and where k k 0; and k k ; are the L () and H () So
Chapter SUBSTRUCTURE PRECONDITIONING FOR POROUS FLOW PROBLEMS R.E. Ewing, Yu. Kuznetsov, R.D. Lazarov, and S. Maliassov. Introduction Let be a convex polyhedral domain in IR, f(x) L () and A(x) be a suciently
More informationA priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces
A priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces A. Bespalov S. Nicaise Abstract The Galerkin boundary element discretisations of the
More informationK. BLACK To avoid these diculties, Boyd has proposed a method that proceeds by mapping a semi-innite interval to a nite interval [2]. The method is co
Journal of Mathematical Systems, Estimation, and Control Vol. 8, No. 2, 1998, pp. 1{2 c 1998 Birkhauser-Boston Spectral Element Approximations and Innite Domains Kelly Black Abstract A spectral-element
More informationSubstructuring Preconditioning for Nonconforming Finite Elements engineering problems. For example, for petroleum reservoir problems in geometrically
SUBSTRUCTURING PRECONDITIONING FOR FINITE ELEMENT APPROXIMATIONS OF SECOND ORDER ELLIPTIC PROBLEMS. I. NONCONFORMING LINEAR ELEMENTS FOR THE POISSON EQUATION IN A PARALLELEPIPED R.E. Ewing, Yu. Kuznetsov,
More informationVol. 29, No. 3, pp. 1{8, May WHERE BEST TO HOLD A DRUM FAST *
SIAM J. OPTIMIATION c 998 Society for Industrial and Applied Mathematics Vol. 29, No. 3, pp. {8, May 998 WHERE BEST TO HOLD A DRUM FAST * STEVEN J. COX y AND PAUL X. UHLIG z Abstract. Allowed to fasten
More information2 used to model the movement of water in groundwater reservoirs. However, it provides no information on the motion of air. If contaminant transport is
FULLY-DISCRETE FINITE ELEMENT ANALYSIS OF MULTIPHASE FLOW IN GROUNDWATER HYDROLOGY Zhangxin Chen* and Richard E. Ewingy Abstract. This paper deals with development and analysis of a fully discrete nite
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationNumerical Schemes from the Perspective of Consensus
Numerical Schemes from the Perspective of Consensus Exploring Connections between Agreement Problems and PDEs Department of Electrical Engineering and Computer Sciences University of California, Berkeley
More informationDynamic Green Function Solution of Beams Under a Moving Load with Dierent Boundary Conditions
Transaction B: Mechanical Engineering Vol. 16, No. 3, pp. 273{279 c Sharif University of Technology, June 2009 Research Note Dynamic Green Function Solution of Beams Under a Moving Load with Dierent Boundary
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationAbstract. A front tracking method is used to construct weak solutions to
A Front Tracking Method for Conservation Laws with Boundary Conditions K. Hvistendahl Karlsen, K.{A. Lie, and N. H. Risebro Abstract. A front tracking method is used to construct weak solutions to scalar
More informationby Ron Buckmire Mathematics Department Occidental College Los Angeles, CA , U.S.A.
A new nite-dierence scheme for singular dierential equations in cylinical or spherical coordinates by Ron Buckmire Mathematics Department Occidental College Los Angeles, CA 90041-3338, U.S.A. December
More informationIterative procedure for multidimesional Euler equations Abstracts A numerical iterative scheme is suggested to solve the Euler equations in two and th
Iterative procedure for multidimensional Euler equations W. Dreyer, M. Kunik, K. Sabelfeld, N. Simonov, and K. Wilmanski Weierstra Institute for Applied Analysis and Stochastics Mohrenstra e 39, 07 Berlin,
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationThe Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-Zero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In
The Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-ero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In [0, 4], circulant-type preconditioners have been proposed
More informationNumerische Mathematik
Numer. Math. (2003) 94: 195 202 Digital Object Identifier (DOI) 10.1007/s002110100308 Numerische Mathematik Some observations on Babuška and Brezzi theories Jinchao Xu, Ludmil Zikatanov Department of Mathematics,
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationThe meshless standard and hypersingular boundary node methods applications to error estimation and adaptivity in three-dimensional problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2001; 50:2233 2269 The meshless standard and hypersingular boundary node methods applications to error estimation and
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
More informationthe sum of two projections. Finally, in Section 5, we apply the theory of Section 4 to the case of nite element spaces. 2. Additive Algorithms for Par
ON THE SPECTRA OF SUMS OF ORTHOGONAL PROJECTIONS WITH APPLICATIONS TO PARALLEL COMPUTING PETTER E. BJRSTAD y AND JAN MANDEL z Abstract. Many parallel iterative algorithms for solving symmetric, positive
More information2 ZHANGXIN CHEN AND RICHARD EWING dierentiating the resulting pressure and multiplying it by the rough coecient. This approach generates a rough and i
FROM SINGLE-PHASE TO COMPOSITIONAL FLOW: APPLICABILITY OF MIXED FINITE ELEMENTS Zhangxin Chen and Richard E. Ewing Abstract. In this paper we discuss the formulation of the governing equations that describe
More informationMathematical Institute, University of Utrecht. The problem of estimating the mean of an observed Gaussian innite-dimensional vector
On Minimax Filtering over Ellipsoids Eduard N. Belitser and Boris Y. Levit Mathematical Institute, University of Utrecht Budapestlaan 6, 3584 CD Utrecht, The Netherlands The problem of estimating the mean
More informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationMathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr
HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL
More informationWALTER LITTMAN BO LIU kind of problem is related to the deformation of a membrane with strings attached to several sides of the polygon. In other word
THE REGULARITY AND SINGULARITY OF SOLUTIONS OF CERTAIN ELLIPTIC PROBLEMS ON POLYGONAL DOMAINS Walter Littman Bo Liu School of Math., University of Minnesota, Minneapolis, MN55455 Abstract. The regularity
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationLecture 2: Review of Prerequisites. Table of contents
Math 348 Fall 217 Lecture 2: Review of Prerequisites Disclaimer. As we have a textbook, this lecture note is for guidance and supplement only. It should not be relied on when preparing for exams. In this
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationPseudo-divergence-free element free Galerkin method for incompressible fluid flow
International Workshop on MeshFree Methods 23 Pseudo-divergence-free element free Galerkin method for incompressible fluid flow Y. Vidal () and A. Huerta (2) Abstract: Incompressible modelling in finite
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More information3rd Int. Conference on Inverse Problems in Engineering. Daniel Lesnic. University of Leeds. Leeds, West Yorkshire LS2 9JT
Inverse Problems in Engineering: Theory and Practice 3rd Int. Conference on Inverse Problems in Engineering June 13-18, 1999, Port Ludlow, WA, USA APP RETRIEVING THE FLEXURAL RIGIDITY OF A BEAM FROM DEFLECTION
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationLimit Analysis with the. Department of Mathematics and Computer Science. Odense University. Campusvej 55, DK{5230 Odense M, Denmark.
Limit Analysis with the Dual Ane Scaling Algorithm Knud D. Andersen Edmund Christiansen Department of Mathematics and Computer Science Odense University Campusvej 55, DK{5230 Odense M, Denmark e-mail:
More informationIN p-version AND SPECTRAL ELEMENT METHODS MARIO A. CASARIN
DIAGONAL EDGE PRECONDITIONERS IN p-version AND SPECTRAL ELEMENT METHODS MARIO A. CASARIN Abstract. Domain decomposition preconditioners for high-order Galerkin methods in two dimensions are often built
More informationANDREA TOSELLI. Abstract. Two-level overlapping Schwarz methods are considered for nite element problems
OVERLAPPING SCHWARZ METHODS FOR MAXWELL'S EQUATIONS IN THREE DIMENSIONS ANDREA TOSELLI Abstract. Two-level overlapping Schwarz methods are considered for nite element problems of 3D Maxwell's equations.
More information