Numerical Integration of Functions with Endpoint Singularities and/or Complex Poles in 3D Galerkin Boundary Element Methods
|
|
- Andrew Morton
- 6 years ago
- Views:
Transcription
1 Publ. RIMS, Kyoto Univ. 4 (25), Numerical Integration of Functions with Endpoint Singularities and/or Complex Poles in 3D Galerkin Boundary Element Methods To Professor Masatake Mori for his outstanding contribution to numerical integration By Giovanni Monegato and Letizia Scuderi Abstract In this paper we propose special strategies to compute D integrals of functions having weakly or strong singularities at the endpoints of the interval of integration or complex poles close to the domain of integration. As application of the proposed strategies, we compute a four dimensional integral arising from 3D Galerkin boundary element methods (BEM) applied to hypersingular boundary integral equations.. Introduction In the computation of integrals and in the numerical solution of integral equations, one often has to deal with the numerical integration of functions with endpoint weak singularities, or with analytic functions having complex conjugate poles near the domain of integration. Nonlinear changes of variables have shown to be a very effective tool to face these problems. Indeed, in the case of weak singularities a proper change of variable can make the integrand function arbitrarily smooth. This, combined with the use of a standard rule Communicated by H. Okamoto. Received September 28, 24. Revised February 22, Mathematics Subject Classification(s): 65D32, 65R2, 65N38. This work was supported by the Ministero dell Istruzione, dell Università e della Ricerca of Italy and by the INDAM 23 project Modellistica Numerica per il Calcolo Scientifico e Applicazioni Avanzate. Dipartimento di Matematica, Politecnico di Torino, Corso Duca degli Abruzzi 24, 29 Torino, Italy. c 25 Research Institute for Mathematical Sciences, Kyoto University. All rights reserved.
2 87 Giovanni Monegato and Letizia Scuderi like the Gauss-Legendre one or the composite trapezoidal formula, allows to achieve high accuracy using a low number of abscissas (see [2], [8], [9], [], [4], [6]). Here we recall some of the known changes of variable, all having as co-domain the interval (, ), which have been proposed in the literature. For some comparisons (see [9]). In [8], the following polynomial transformation, which is a generalization of the one proposed in [6] to construct lattice rules for multiple integration on the hypercube, has been used: (.) ϕ (t) = (p + q ) (p )(q ) t u p ( u) q du, t, p,q. The integral in (.) could be computed, for any t (, ], using recurrence relations (see [9]). In any case it can be evaluated more efficiently, for the given values of p and q, bymeansofthen-point Gauss-Legendre rule, with n = p+q 2. Then, we recall a trigonometric transformation of the form t (.2) ϕ 2 (t) = (sin π 2 u)p (cos π 2 u)q du (sin, t, p,q. π 2 u)p (cos π 2 u)q du This latter generalizes a trigonometric transformation proposed in [4]. Also in this case the integral in (.2) could be computed for any t (, ] by using proper recurrence relations or by a direct computation using a Gauss-Legendre rule. A third transformation, of rational type, used in [], [3] in connection with the solution of certain integral equations, and in [4] to evaluate Hadamard finite-part integrals, is the following one: (.3) ϕ 3 (t) = t p t p, t, p,q. +( t) q The last transformation we recall is of exponential type: (.4) ϕ 4 (t) = 2 + ( π ) 2 tanh 2 sinh(t), <t<. It is the well-known double exponential (DE) transformation (see [6], [], [5]). In this paper we examine the use of nonlinear changes of variable for the numerical evaluation of some 4-dimensional integrals arising in the numerical solution of hypersingular boundary integral equations. In particular we consider
3 Numerical Integration in 3D Galerkin BEM 87 3D Galerkin boundary element methods, where a key issue is the efficient and fast computation of integrals of the form k(x, y) (.5) dy dx, j x y 3 i where i, j are triangular elements, possibly curved, and k(x, y) is bounded. When i and j are disjoint and not too close to each other, the computation of the above integral can be performed by standard rules. The difficult cases are those where (i) i and j are coincident, or share a common edge, or (ii) are disjoint but very close to each other. When i and j have only a common point, but are not close to each other at some other points, the computation can be classified standard. Many papers have been devoted to this topic (see, for example, [], [5], [2], [3]), although none of them seems to treat the case (ii), which is considered by some users the most difficult case. However in our opinion none of them are fully satisfactory because generally they either require some analytical calculation which is by non means trivial, or do not take into account the presence of complex poles which may effect adversely the numerical computation. Indeed the strategy adopted by these methods is that of reducing the computation of the multidimensional integral to the evaluation of nested D integrals of analytic functions. However this is still not satisfactory because it is by now well-known that when integrating an analytic function, the position of its poles is crucial for the performance of a quadrature rule. When an analytic function has poles which are outside the domain of integration but lie very close to it, and this is what happens in the computation of some of our integrals, its numerical integration by means of a quadrature rule can give very poor accuracy if one does not adopt a special strategy. Inthispaperweexaminethecaseof(.5)with i, j coinciding. The other cases of i, j having a common edge or disjoint but very close to each other will be treated in a following paper. We believe that our strategy can be applied successfully also to these cases. 2. One-dimensional Integration Rules There are two typical approaches to compute (.5). One uses the global cartesian coordinates (see [5]), the other uses local polar coordinates (see [2]). Both of them reduce (.5) to four nested one-dimensional integrals, whose integrand functions are analytical in the first case, and analytical or with logsingularities at the endpoints of the interval of integration in the second case.
4 872 Giovanni Monegato and Letizia Scuderi Unfortunately, in both cases the generated analytical functions have complex poles, depending upon outer variables or the form of the triangles, which may be arbitrarily close to the interval of integration. Therefore, for both approaches the Gaussian quadrature rules perform very poorly unless one considers a high number of quadrature nodes. Recalling that we have to compute four nested one-dimensional integrals, we can use only a few (let us say, at the most 6 ) quadrature nodes for each integral. Notice that quadrature nodes for each integral means a total of 4 quadrature nodes to compute (.5), which represents only one element of the system matrix arising from the Galerkin boundary element method. For this reason in this section we propose special strategies, which allow us to compute the integrals with the above-mentioned characteristics using very few Gaussian quadrature nodes. In particular we propose to use variable transformations of the type (.), either to make the integrand functions smoother at the endpoints of the interval of integration, or to move poles away from the interval. Then, the Gauss-Legendre rule is applied. In these situations the DE-rule, which is obtained by combining (.4) with the (truncated) trapezoidal rule (see [6]), has certainly a much higher (exponential) rate of convergence. However, in the particular application we will consider in Section 3, it is of key importance to obtain a good accuracy with a very small number of function evaluations. For such small values of n, let us say n, our approach appears more efficient. A typical case is an integral of a function having only a log-singularities at the endpoints of the interval of integration. For example, (2.) I = e x log x x By introducing in (2.) the change of variable (.), that here we write in the form (2.2) x = γ q,q, (t) := (q + q ) (q )(q ) t dx. x q ( x) q dx, q,q, and then applying the n-point Gauss-Legendre rule we obtain the improvements reported in Table. Notice that in (2.2) for q =wehave (2.3) γ q,, (t) =tq =: γ q (t) and for q = (2.4) γ,q, (t) = ( t)q =: γ q (t);
5 Numerical Integration in 3D Galerkin BEM 873 the choices (2.3) and (2.4) allow to smooth the singularities at the endpoints and, respectively. If q = q = q, weset (2.5) γ q, := γq,q,. n q = q =2 q =3 q = Table. Relative errors in I = eγ(t) log γ(t) γ(t) γ (t) dt, γ(t) =γ q, (t). In all tables the sign means that full relative accuracy (i.e. 4 significant digits in our case) has been achieved. Another situation that we will meet is that of integrals of the form (2.6) I 2 = f(x) x 2 + ε 2 dx with ε small and f(x) smooth. To compute efficiently this integral it is convenient to introduce preliminarily the simple change of variable x = γ q (t), q>, and then to apply the n-point Gauss-Legendre rule. For example, if we take f(x) =e x we obtain the relative errors of Table 2. ε = ε = 3 ε = 5 n q = q =2 q = q =4 q = q = Table 2. Relative errors in I 2 = qt q e tq t 2q +ε 2 dt. Since our D rules will be applied to a quadruple integral, it is important to obtain the required accuracy using a small number of nodes. In the integral
6 874 Giovanni Monegato and Letizia Scuderi (2.6), when ε is very small, for example in the above case ε = 5,thevalue of q one has to take makes the new integrand function so flat around the origin, that to achieve a good accuracy one has to take a value of n which is not any longer small. In this situation it is more convenient to proceed as follows. Write ( ε ) f(x) (2.7) I 2 = + ε x 2 + ε 2 dx; then apply the n-point Gauss-Legendre rule to the first integral on the righthand side, and to the second one after having introduced in it the change of variable x = γ q (t). Subdivision (2.7) is of key importance for the following reasons. In (2.6) the transformation x = γ q (t) leaves unchanged the interval of integration, but as q gets larger, the modulus of the poles tends to, and after a certain value of q some of the poles move again towards the interval of integration. In (2.7) this unwanted phenomenon does not arise. Indeed, in the first integral on the right-hand side the distance of the two poles from the interval of integration is equal to the length of this latter, while the length of the interval of integration of the second integral, after having introduced the transformation x = γ q (t), tends to zero as q tends to infinity. The effectiveness of this new approach is shown by the results in Table 3. ε = ε = 3 ε = 5 2n q =3 q = q =4 q = q =7 q = Table 3. Relative errors in I 2 = ε e x ε q x 2 +ε 2 dx + qt q e tq t 2q +ε 2 dt. Similarly, if the integral has the form b f(x) (2.8) I 3 = a (x 2 + ε 2 ) α dx with a<, b>andα>, then we proceed as in (2.7) by splitting (2.8) as follows ( ε ε ) b f(x) (2.9) I 3 = + + a ε ε (x 2 + ε 2 ) α dx. We apply the n-point Gauss-Legendre quadrature rule to each integral, after having introduced the change of variable x = t q in the first and in the third
7 Numerical Integration in 3D Galerkin BEM 875 integral to the right-hand side of (2.9). In the following Table 4, we have reported the relative errors obtained by this last numerical procedure when in (2.8) f, a =.5, b =.5 andα =3/2. In Table 4 and below GL m denotes the m-point Gauss-Legendre rule. ε = ε = 3 ε = 5 3n GL 3n q =2 GL 3n q =4 q = GL 3n q =7 q = Table 4. Relative errors in I 3 =( a q ε q q ε q + b + ε ε (x2 + ε 2 ) 3/2 dx, a =.5, b =.5. ) qtq (t 2q + ε 2 ) 3/2 dt The integral we have to compute could also be of the type (2.) I 4 = f(x) [(x r) 2 + ε 2 ] α dx with <r<andα>. Notice that also I 3 can be rewritten in this form and vice versa. In this case we can proceed as follows. Write first and then (2.) I 4 = r r ( r I 4 = + r r ) f(x) [(x r) 2 + ε 2 ] α dx, f( r r ( t)) [(t r) 2 +( r r ε)2 ] dt + f(x) α r [(x r) 2 + ε 2 ] α dx. This splitting is suggested by the remark that, when one uses Gaussian rules, poles with real parts close to the endpoints of the interval of integration are less adverse than poles with a real part in the central part of the same interval. Finally, we introduce, as for I 2, the change of variable t, x = u q, q, in both integrals and apply the same Gauss-Legendre rule to each of them. Since this strategy is also alternative to a splitting of the type (2.9), we have applied it to the same integral considered in Table 4. Therefore, in Table 5 we have reported the relative errors obtained by taking in (2.) f, r =.5 and α =3/2.
8 876 Giovanni Monegato and Letizia Scuderi ε = ε =5 2 2n GL 2n q = q =5 GL 2n q = q = Table 5. Relative errors in I 4 = r r +( r r r ε)2] 3/2 du + r q quq f( r r ( uq )) [ (u q r) 2 q quq f(u q )[(u q r) 2 + ε 2 ] 3/2 du, r =.5. From a comparison between the Tables 4 and 5, it appears that (2.) is more effective when ε is not too small; let us say of order greater than 2. Otherwise (2.9) seems to produce better results. The above strategy has been applied also to similar integrals with the algebraic polynomial (of degree 2) replaced by a corresponding trigonometric polynomial. In particular, we have also considered integrals of the type (2.2) I 5 = b a f(x) [(a cos(x)+b sin(x)) 2 +(a 2 cos(x)+b 2 sin(x)) 2 ] α dx. By z we denote the pair of complex conjugate poles of the trigonometric polynomial closest to the interval of integration. In this situation, we have applied the above strategy by choosing r = Re(z ). The results we have obtained by taking in (2.2) f, a =.5, b =,α =3/2 anda =,a 2 = b 2 =, b = c with varying c, are reported in Table 6. c = c =/2 c =/4 c =/8 Re(z ).55 Re(z ).72 Re(z ).77 Re(z ).78 Im(z ) ±.4 Im(z ) ±.24 Im(z ) ±.2 Im(z ) ±.6 2n GL 2n q =5 GL 2n q =5 GL 2n q =5 GL 2n q = Table 6. Relative errors in I 5 =.5 [(c sin(x))2 +(cos(x)+sin(x)) 2 ] 3/2 dx.
9 Numerical Integration in 3D Galerkin BEM 877 Also integrals of the form (2.3) I 6 = ε f(x) x with f(x) smoothandε> small, can be efficiently evaluated using a Gauss- Legendre rule, after having introduced the change of variable x = t q. Some relative errors concerning an integral of this type, with f(x) =e x, are reported in Table 7. dx, ε = ε = 3 ε = 5 n q = q =3 q = q =6 q = q = Table 7. Relative errors in I 6 = ε q qe tq t dt. In the introduction we have remarked that for small values of n, theapproach we have proposed gives slightly better results than the DE-rule. This is due to the fact that the stepsize h of the DE-rule is chosen to optimize the rate of convergence, that is, the behavior of the remainder term as n (see [5]). Indeed, as n the accuracy given by the DE-rule is superior to that given by our numerical integration approach. However, in the applications we are interested in, we ought to achieve a good accuracy by taking a small value of n. Therefore, in the case of the DErule, for a given (low) value of n we would like to choose the parameter h which minimizes the error. For the integrals I through I 5 of this section we have plotted the behavior of the error term as a function of h, forseveralvaluesofn. Inallcaseswehave examined, a proper choice of h has reduced substantially the error. Here we report a few examples, which are however representative of all cases we have examined. In all cases there exists an optimal choice of h which allows to reach very high accuracies using a small value of n. However we do not know a simple formula to compute such value, or a good approximation of it. We note that all figures have been obtained by subdividing the domain of h into equal parts. Therefore the actual minima could be even smaller than those appearing in the figures.
10 878 Giovanni Monegato and Letizia Scuderi 2 relative error stepsize h Figure. Integral I, n =3,h (.45,.65). 4 relative error stepsize h Figure 2. Integral I, n =5,h (.35,.55).
11 Numerical Integration in 3D Galerkin BEM relative error stepsize k (h=k+.48) x 3 Figure 3. Integral I, n =5,h (.4847,.4849). relative error stepsize h Figure 4. Integral I 4, r =.5, ε =, n =3,h (., 2).
12 88 Giovanni Monegato and Letizia Scuderi 8 relative error stepsize k (h=k+.355) x 5 Figure 5. Integral I 4, r =.5, ε =, n =3,h ( ,.35529). 6 7 relative error stepsize k (h=k+.48e 3) x 6 Figure 6. Integral I 4, r =.5, ε = 3, n =3,h (.48238e 3,.4824e 3).
13 Numerical Integration in 3D Galerkin BEM relative error stepsize k (h=k ) x 5 Figure 7. Integral I 5, ε =/2, n =3,h ( , ). 6 7 relative error stepsize k (h=k+2.4) x 3 Figure 8. Integral I 5, ε =/8, n =3,h (2.4832, ).
14 882 Giovanni Monegato and Letizia Scuderi The last integral we need to consider is (2.4) I 7 = = R(y) f(x) x defined as Hadamard finite-part, where y is an outer variable, which can make R(y) very close to zero. Since a change of variable does not leave unchanged the value of the integral, and we do not want to construct an ad hoc rule for each value of the variable y, we suggest to proceed as follows. If f() exists and can be easily computed, then write dx, (2.5) R(y) = f(x) x dx = R(y) f(x) f() x dx + f() log R(y) and apply the Gauss-Legendre to the integral on the right-hand side of (2.5). If however this is not the case, then write (2.6) R(y) = f(x) x dx = ( r = + R(y) r ) f(x) x dx. where r>has a suitable fixed value. To compute the first integral on the right-hand side we use a standard product rule for Hadamard finite part integrals, whose coefficients have been computed once for all ([7]). The second integral can be computed by the technique suggested for (2.3) with a large q (for example, q = ) when f is smooth. The choice of r may depend on f. Indeed, if f has a pair of complex poles with the real part fixed, independent of y and belonging to the interval of integration (,R(y)), it is convenient to take r equal to the real part of the poles. If the real part varies with y in a little range of (,R(y)) we can choose r as an intermediate value of the range; otherwise, it is convenient to use a product formula without splitting. Finally, if f is smooth, we suggest to choose r small. In the Tables 8 and 9 we have reported the relative errors obtained first by applying directly the product rule with n quadrature points (PR n ) to I 5 in (2.4), and then by splitting I 5 as in (2.6) and computing the integrals to the left-hand side as suggested above (PR n + GL n ). In particular, in Table 8 we report the relative errors corresponding to the choice f(x) = f (x) = /[(x r) 2 + ε 2 ]withr =. andε varying; in Table 9 we give those corresponding to f(x) =f 2 (x) =e x and r =. in (2.6).
15 Numerical Integration in 3D Galerkin BEM 883 ε = ε = 2 ε = 3 2n PR 2n PR n + GL n PR 2n PR n + GL n PR 2n PR n + GL n Table 8. Relative errors in I 7 = R(y) =( r + R(y) ) [(x.)2 +ε 2 ] r [(x.) 2 +ε 2 ] x dx x dx with r =., R(y) =.5. R(y) = 2 R(y) =.5 2n PR 2n PR n + GL n PR 2n PR n + GL n Table 9. Relative errors in I 7 = R(y) + (R(y)) r q q e x x qe tq t dt with r =., q =. dx = r e x x dx 3. The Four-dimensional Integral Theintegralwehavetocomputehasthefollowingform (3.) I = = d x = N i( x)n j (ȳ) x ȳ 3 dȳ, where for simplicity we assume to be a triangle with corners at (, ), (a,a 2 ), (b,b 2 )and{n i ( x)}, {N j (ȳ)} are the shape and the test functions, respectively. Our approach however easily extends to the case of a curved triangle, given by a smooth parametric representation. By the linear transformation x = Ax, ȳ = Ay with ( ) a b A = a 2 b 2 we map the domain into the reference triangle T having corners at (, ), (, ) (, ). In this case the integral (3.) can be replaced by
16 884 Giovanni Monegato and Letizia Scuderi (3.2) I T = A 2 = dx = N i(x)n j (y) A(x y) 3 dy, T T where A =det(a), N i (x) :=N i (Ax) andn j (y) :=N j (Ay). The value of I T is generally different from that of I ; however the sum of the values of the integrals defined over all boundary elements i is equal to the sum of the corresponding integrals of form (3.2). Therefore from now on we refer to the form (3.2). To compute (3.2) we split it as follows: [ (3.3) I T = A 2 N i (x) dx N j(y) N j (x) T T A(x y) 3 dy ] dy + N i (x)n j (x) dx = T T A(x y) 3 =: A 2 [I T, + I T,2 ] since the inner integral in I T,2 can be evaluated analytically and the strong singularity in I T, is of Cauchy type. As regards I T,, first we introduce the following polar coordinates (3.4) { y = x + r cos(ϑ) y 2 = x 2 + r sin(ϑ) where r R(ϑ), ϑ 2π. Thenwesplittheinterval(, 2π) intothree subintervals (ϑ,ϑ ), (ϑ,ϑ 2 )and(ϑ 2,ϑ 3 ), with ( ) x2 ϑ = π +arctan x ( ) x2 ϑ = arctan x (3.5) ( ) x2 ϑ 2 = π arctan x ( ) x2 ϑ 3 = π +arctan. For ϑ [ϑ l,ϑ l ], l =, 2, 3, we have r R l with (3.6) R l := R l (x,x 2,ϑ)= d l(x,x 2 ) p l (ϑ) x
17 Numerical Integration in 3D Galerkin BEM 885 where (3.7) d (x,x 2 )=x 2 d 2 (x,x 2 )= ( x x 2 ), 2 d 3 (x,x 2 )=x p (ϑ) = sin(ϑ) p 2 (ϑ) = [sin(ϑ)+cos(ϑ)]. 2 p 3 (ϑ) = cos(ϑ) Our integral I T, canthenbeexpressedasfollows: (3.8) I T, = 3 ϑl x dx N i (x,x 2 ) dx 2 Rl l= ϑ l dϑ = N j(x + r cos(ϑ),x 2 + r sin(ϑ)) N j (x,x 2 ) r 2 dr, f(ϑ) where f(ϑ) =[(a cos(ϑ)+b sin(ϑ)) 2 +(a 2 cos(ϑ)+b 2 sin(ϑ)) 2 ] 3/2. Notice that f(ϑ) vanishes when with α = a b + a 2 b 2 b 2 + b2 2 tan(ϑ) =α ± iβ,, β = a 2b a b 2 b 2 +. b2 2 In particular the poles of f(ϑ) depend on the corners of the triangle. Therefore their distance from the interval of integration with respect to the variable ϑ depends upon the shape of and could be very small. For example, in the case of a triangle, whose corners are at (, ), (a, ), (, ), we have approximately the following complex poles: (. + kπ) ±.4 ı when a =, (.85 + kπ) ±.24 ı when a =/2 and(.8 + kπ) ±.2 ı when a =/4, k ZZ. To proceed we set N r j (x,x 2 ; r, ϑ) = N j(x + r cos(ϑ),x 2 + r sin(ϑ)) N j (x,x 2 ) r
18 886 Giovanni Monegato and Letizia Scuderi and write (3.9) I T, = x dx N i (x,x 2 ) dx 2 ϑl Rl 3 [f(ϑ)] dϑ = N j r(x,x 2 ; r, ϑ) dr l= ϑ l r x 3 ϑl = dx N i (x,x 2 ) dx 2 [f(ϑ)] l= ϑ l [ Rl Nj r(x,x 2 ; r, ϑ) Nj r(x ],x 2 ;,ϑ) dr + Nj r (x,x 2 ;,ϑ)log(r l ) dϑ r Remark. Notice that in the case of {N j } linear, that is N j = N l j, j =, 2, 3, and (3.) N(x l,x 2 )=x, N2(x l,x 2 )=x 2, N3(x l,x 2 )= x x 2, we have N r (x,x 2 ; r, ϑ) =cos(ϑ), N2 r (x,x 2 ; r, ϑ) =sin(ϑ), N3 r (x,x 2 ; r, ϑ) = [sin(ϑ)+cos(ϑ)]. Therefore Nj r depends only upon ϑ and the integral over (,R l) in the first expression of I T, in (3.9) can be performed analytically. By denoting we write D j (x,x 2 ; r, ϑ) := N j r(x,x 2 ; r, ϑ) Nj r(x,x 2 ;,ϑ), r I T, = 3 IT,, l l=
19 Numerical Integration in 3D Galerkin BEM 887 where (3.) x ϑl IT, l = dx N i (x,x 2 ) dx 2 [f(ϑ)] ϑ l [ ] Rl D j (x,x 2 ; r, ϑ) dr + Nj r (x,x 2 ;,ϑ)log(r l ) dϑ, l =, 2, 3. and (3.2) we have Remark. In the case of {N j } quadratic, that is N j = N q j, j =,..., 6, N q (x,x 2 )=x (2x ), N q 2 (x,x 2 )=x 2 (2x 2 ), N q 3 (x,x 2 )=( x x 2 )(2( x x 2 ) ), N q 4 (x,x 2 )=4x x 2, N q 5 (x,x 2 )=4x 2 ( x x 2 ), N q 6 (x,x 2 )=4x ( x x 2 ), D (x,x 2 ; r, ϑ) =2cos 2 (ϑ), D 2 (x,x 2 ; r, ϑ) =2sin 2 (ϑ), D 3 (x,x 2 ; r, ϑ) =2[sin(ϑ)+cos(ϑ)] 2, D 4 (x,x 2 ; r, ϑ) =4sin(ϑ)cos(ϑ), D 5 (x,x 2 ; r, ϑ) = 4sin(ϑ)[sin(ϑ)+cos(ϑ)], D 6 (x,x 2 ; r, ϑ) = 4cos(ϑ)[sin(ϑ)+cos(ϑ)]. Therefore, since D j (x,x 2 ; r, ϑ) does not depend upon r, the integral over (,R l ) in (3.) can be performed analytically. Since R l takes different values for l =, 2, 3, the computation of I l T, have to be performed using different strategies according to these values. We first consider I T,, where we have set x 2 =( x ) x 2, x 2 (, ): (3.3) I T, = ϑ ( x ) dx N i (x, ( x ) x 2 ) d x 2 [f(ϑ)] [ ] R D j (x, ( x ) x 2 ; r, ϑ) dr + Nj r (x, ( x ) x 2 ;,ϑ)log( R ) dϑ ϑ
20 888 Giovanni Monegato and Letizia Scuderi with ) ( x ) x 2 ϑ = π +arctan(, x ϑ = arctan ( x 2 ), R = ( x ) x 2 sin(ϑ). Taking into account of the behaviors of the integrand functions of each integral in (3.3) and of the behaviors of the endpoints ϑ, ϑ as x, x 2,, the following smoothing transformations have been introduced x = γ q, (x), x (, ), x 2 = γ q 2 (y), y (, ), ϑ = γ q 3 ((b a)t + a), t (, ), whose smoothing exponents q l, l =, 2, 3, will be specified later. In the last change of variable (a, b) denotes the interval of integration with respect to the variable ϑ. If no real parts ϑ R of the poles of f(ϑ) belongto( ϑ, ϑ ), the interval of integration is ( ϑ, ϑ ). Otherwise we split this interval into two subintervals ( ϑ,ϑ R )and(ϑ R, ϑ ) and we apply the above change of variable to each of them, taking (a, b) ( ϑ,ϑ R )and(a, b) (ϑ R, ϑ ). Finally, for the variables x and x 2 the numerical integration is performed using the n-point Gauss-Legendre rule, while for the variable ϑ we have used the 2n-point Gauss-Legendre rule. InthecaseofIT, 2 we have (3.4) I 2 T, = ϑ2 ( x ) dx N i (x, ( x ) x 2 ) d x 2 [f(ϑ)] [ ] R 2 D j (x, ( x ) x 2 ; r, ϑ) dr + Nj r (x, ( x ) x 2 ;,ϑ)log( R 2 ) dϑ ϑ where ϑ = arctan ( x 2 ), ϑ 2 = π arctan ( ( x ) x 2 x R 2 = ( x )( x 2 ) sin(ϑ)+cos(ϑ). ),
21 Numerical Integration in 3D Galerkin BEM 889 In (3.4) we first introduce the smoothing transformations x = γ q 2, (x), x (, ), x 2 = γ q 22 (y), y (, ), ϑ = γ q 23, ((b a)t + a), t (, ), where, as above, (a, b) ( ϑ, ϑ 2 ) is split into ( ϑ,ϑ R ) (ϑ R, ϑ 2 ) whenever there exists a pole of f(ϑ) whose real part ϑ R ( ϑ, ϑ 2 ). Then we apply the same Gauss-Legendre product rule used in the previous case. The integral (3.5) I 3 T, = where ϑ3 ( x ) dx N i (x, ( x ) x 2 ) d x 2 [f(ϑ)] [ ] R 3 D j (x, ( x ) x 2 ; r, ϑ) dr + Nj r (x, ( x ) x 2 ;,ϑ)log( R 3 ) dϑ x ) ( x ) x 2 ϑ 3 = π +arctan(, ϑ 2 ϑ 2 = π arctan ( ( x ) x 2 x R 3 = cos(ϑ), is treated similarly to the previous cases. Here we set x = γ q 3 (x), x (, ), x 2 = γ q 32, (y), y (, ), ϑ = γ q 33, ((b a)t + a), t (, ), and the interval (a, b) ( ϑ 2, ϑ 3 ) is split into ( ϑ 2,ϑ R ) (ϑ R, ϑ 3 )ifϑ R ( ϑ 2, ϑ 3 ). In the case of linear and quadratic functions N j, a choice of the above smoothing exponents q l, q l2, q l3, l =, 2, 3, which leads to very satisfactory results, is reported in Table. However, for basis functions of higher order the choice x q l = q l2 = q l3 =2, l =, 2, 3, is quite satisfactory. As already stated earlier, the inner most integral in I T,2 defined in (3.3), can be evaluated analytically. Indeed, we have ),
22 89 Giovanni Monegato and Letizia Scuderi N i,n j linear I T, I 2 T, I 3 T, i j q q 2 q 3 q 2 q 22 q 23 q 3 q 32 q 33 : : : N i,n j quadratic IT, IT, 2 IT, 3 i j q q 2 q 3 q 2 q 22 q 23 q 3 q 32 q 33 : : : : : : Table. Thebestchoicefortheexponentsq lm, l, m =, 2, 3. (3.6) I T,2 = x dx N i (x,x 2 )N j (x,x 2 ) dx 2 y dy 2 dy {[a (x y )+b (x 2 y 2 )] 2 +[a 2 (x y )+b 2 (x 2 y 2 )] 2 } 3/2 = x [ ax 2 A 2 dx K i,j (x,x 2 ) + b(x 2 ) 2 +2cx (x 2 ) x (x + x 2 ) ax 2 + bx cx ] x 2 a(x ) bx cx 2(x ) dx 2 x x 2 x 2 (x + x 2 ) =: A 2 [I T,2 + I 2 T,2 + I 3 T,2], where K i,j (x,x 2 )=N i (x,x 2 )N j (x,x 2 ), a = a 2 + a 2 2, b = b 2 + b 2 2, c = a b + a 2 b 2. To compute IT,2 we proceed as follows: (3.7) I T,2 = dx x [ K i,j (x,x 2 ) ax 2 + b(x 2 ) 2 +2cx (x 2 ) b ( x x 2 ) K i,j (,x 2 ) ] a + b 2cx K i,j (x, x ) /[x (x + x 2 )] dx 2
23 Numerical Integration in 3D Galerkin BEM 89 dx x b K i,j (,x 2 ) dx 2 x + a + b 2c =: I, T,2 + I,2 T,2 + I,3 T,2 K i,j (x, x ) dx x dx 2 x + x 2 As regards I, T,2, first we introduce the changes of variable x = η and x 2 = ξ η +, with ξ and η ξ, and then we apply the Duffy transformation ξ = uv and η =( u)v, with u, v ; hence, we have (3.8) ξ [ I, T,2 = dξ K i,j (η, ξ η +) aη 2 + b(ξ + η) 2 2cη(ξ + η) bξk i,j (, ξ η +) ] a + b 2cηK i,j (η, η) /(ξη) dη = du [ K i,j (( u)v, v) a( u) 2 + b 2c( u) buk i,j (, v) ] a + b 2c( u)k i,j (( u)v, ( u)v) /[u( u)] dv. Moreover, [ x K (3.9) I,2 T,2 = i,j (,x 2 ) dx 2 ] K i,j(,x 2 ) dx 2 b dx x and (3.2) I,3 T,2 = a + b 2c K i,j (x, x )log( x ) dx. To compute IT,2 2, we first rewrite it as follows: x [ (3.2) IT,2 2 = dx K i,j (x,x 2 ) ax 2 + bx2 2 +2cx x 2 ax K i,j (x, ) ] bx 2 K i,j (,x 2 ) /(x x 2 ) dx 2 a b =: I 2, T,2 + I2,2 T,2 + I,2 T,2, K i,j (x, ) log( x ) dx [ x K i,j (,x 2 ) dx 2 ] K i,j(,x 2 ) dx 2 x dx
24 892 Giovanni Monegato and Letizia Scuderi being the latter integral equal to that defined in (3.9). We then introduce in I 2, T,2 the Duffy transformation x = uv and x 2 =( u)v, u, v ; we obtain (3.22) [ I 2, T,2 = du K i,j (uv, ( u)v) au 2 + b( u) 2 +2cu( u) auk i,j (uv, ) ] b ( u) K i,j (, ( u)v) /[u( u)] dv. Finally, we treat IT,2 3 in an analogous way to I T,2 ; therefore, first we write (3.23) I 3 T,2 = dx x [ K i,j (x,x 2 ) a(x ) 2 + bx cx 2(x ) a ( x x 2 ) K i,j (x, ) ] a + b 2cx 2 K i,j (x, x ) /[x 2 (x + x 2 )] dx 2 a a + b 2c K i,j (x, ) log( x ) dx =: I 3, T,2 + I2,2 T,2 + I,3 T,2, K i,j (x, x )log( x ) dx being the second and the third integral at the right-hand side of (3.23) equal to I 2,2 T,2 and I,3 T,2 defined in (3.2) and (3.2), respectively. Hence, we rewrite I 3, T,2 by introducing first the changes of variable x = ξ η + andx 2 = η, with ξ and η ξ, and then the Duffy transformation ξ = uv and η =( u)v, with u, v : (3.24) ξ [ I 3, T,2 = dξ K i,j ( ξ η +,η,) a(ξ + η) 2 + bη 2 2cη(ξ + η) bξk i,j ( ξ η +, ) ] a + b 2cηK i,j ( ξ η+,ξ+η) /(ξη) dη [ = du K i,j ( v, ( u)v) a + b( u) 2 2c( u) auk i,j ( v, ) ] a + b 2c ( u) K i,j ( v, v) /[u( u)] dv.
25 Numerical Integration in 3D Galerkin BEM 893 By taking into account equations (3.7)-(3.24), from (3.6) we derive that (3.25) I T,2 = A 2 [I, T,2 + I2, T,2 + I3, T,2 +2(I,2 T,2 + I2,2 T,2 + I,3 T,2 )]. The first three integrals are evaluated by using the n-point Gauss-Legendre formula, while the remaining ones are computed by using the same Gaussian rule, after having introduced our smoothing change of variable (2.4) with q =3. In Table we have reported the relative errors for I Ta obtained in the case of the triangle T a with corners at (, ), (a, ), (, ), with a varying in {, /2, /4} and by using linear functions N i,n j. i j n T T /2 T / Table. Relative errors for I Ta,whereT a has corners at (,), (a,), (,) and with N i,n j linear. In Table 2 we have reported the relative errors obtained for I Ta,onthe same triangle T a of Table and by using quadratic functions N i,n j. Remark. We explicitly remark that if we apply the above numerical procedure without introducing smoothing transformations, we generally obtain orders of relative accuracy halved with respect to those reported in the Tables and 2. Moreover, notice that for simplicity we have considered in Table only a choice for the smoothing exponents with varying {N j }, j =,..., 6.
26 894 Giovanni Monegato and Letizia Scuderi i j n T T /2 T /4 min max min max min max : : : : : : Table 2. Minimum and maximum relative errors for I Ta,whereT a has corners at (,), (a,), (,) and with N i,n j quadratic. Obviously this fact has involved in some particular cases a loss of precision; indeed, in these cases a suitable choice would have given a more accurate relative precision than that of the Tables and 2. References [] Andrä, H. and Schnack, E., Integration of singular Galerkin-type boundary element integrals for 3D elasticity problems, Numer. Math., 76 (997), [2] Diligenti, M. and Monegato, G., Integral evaluation in the BEM solution of (hyper)singular integral equations. 2D problems on polygonal domains, J. Comput. Appl. Math., 8 (997), [3] Elliott, D. and Prössdorf, S., An algorithm for the approximate solution of integral equations of Mellin type, Numer. Math., 7 (995), [4] Elliott, D. and Venturino, E., Sigmoidal transformations and the Euler-Maclaurin expansion for evaluating certain Hadamard finite-part integrals, Numer. Math., 77 (997), [5] Haas, M. and Kuhn, G., A symmetric Galerkin BEM implementation for 3D elastostatic problems with an extension to curved elements, Comp. Mech., 28 (22), [6] Korobov, N. M., Number-Theoretic Methods of Approximate Analysis, GIFL,Moscow, 963. [7] Monegato, G., Numerical evaluation of hypersingular integrals, J. Comput. Appl. Math., 5 (994), 9-3. [8] Monegato, G. and Scuderi, L., High order methods for weakly singular integral equations with non smooth input functions, Math. Comp., 67 (998), [9], Numerical integration of functions with boundary singularities, J. Comput. Appl. Math., Special Issue: Numerical Evaluation of Integrals, D. Laurie, R. Cools eds., 2 (999), 2-24.
27 Numerical Integration in 3D Galerkin BEM 895 [] Mori, M. and Sugihara, M., The double-exponential transformation in numerical analysis, J. Comput. Appl. Math., Numerical Analysis 2, Vol. V, Quadrature and orthogonal polynomials, 27 (2), [] Prössdorf, S. and Rathsfeld, A., Quadrature methods for strongly elliptic Cauchy singular integral equations on an interval, in: H. Dym (Ed.), The Gohberg Anniversary Collection, Vol. 2: Topics in Analysis and Operator Theory, Birkhäuser, Basel, (99), [2] Sauter, S. and Lage, C., Transformation of hypersingular integrals and black-box cubature, Math. Comp., 7 (2), [3] Sauter, S. and Schwab, C., Quadrature for hp-galerkin BEM in RI 3, Numer. Math., 78 (997), [4] Sidi, A., A new variable transformation for numerical integration, in: H. Braess, G. Hämmerlin, eds., Numerical Integration IV,Birkhäuser Verlag, Basel, ISNM 2 (993), [5] Sugihara, M., Optimality of the double exponential formula - functional analysis approach, Numer. Math., 75 (997), [6] Takahasi, H. and Mori, M., Double exponential formulas for numerical integration, Publ. RIMS, Kyoto Univ., 9 (974),
Dipartimento di Matematica
Dipartimento di Matematica L. SCUDERI ON THE COMPUTATION OF NEARLY SINGULAR INTEGRALS IN 3D BEM COLLOCATION Rapporto interno N. 14, giugno 2006 Politecnico di Torino Corso Duca degli Abruzzi, 24-10129
More informationb n x n + b n 1 x n b 1 x + b 0
Math Partial Fractions Stewart 7.4 Integrating basic rational functions. For a function f(x), we have examined several algebraic methods for finding its indefinite integral (antiderivative) F (x) = f(x)
More informationChapter 7: Techniques of Integration
Chapter 7: Techniques of Integration MATH 206-01: Calculus II Department of Mathematics University of Louisville last corrected September 14, 2013 1 / 43 Chapter 7: Techniques of Integration 7.1. Integration
More information(x 3)(x + 5) = (x 3)(x 1) = x + 5. sin 2 x e ax bx 1 = 1 2. lim
SMT Calculus Test Solutions February, x + x 5 Compute x x x + Answer: Solution: Note that x + x 5 x x + x )x + 5) = x )x ) = x + 5 x x + 5 Then x x = + 5 = Compute all real values of b such that, for fx)
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationPURE MATHEMATICS AM 27
AM Syllabus (014): Pure Mathematics AM SYLLABUS (014) PURE MATHEMATICS AM 7 SYLLABUS 1 AM Syllabus (014): Pure Mathematics Pure Mathematics AM 7 Syllabus (Available in September) Paper I(3hrs)+Paper II(3hrs)
More informationMath 181, Exam 2, Study Guide 2 Problem 1 Solution. 1 + dx. 1 + (cos x)2 dx. 1 + cos2 xdx. = π ( 1 + cos π 2
Math 8, Exam, Study Guide Problem Solution. Use the trapezoid rule with n to estimate the arc-length of the curve y sin x between x and x π. Solution: The arclength is: L b a π π + ( ) dy + (cos x) + cos
More informationFall 2013 Hour Exam 2 11/08/13 Time Limit: 50 Minutes
Math 8 Fall Hour Exam /8/ Time Limit: 5 Minutes Name (Print): This exam contains 9 pages (including this cover page) and 7 problems. Check to see if any pages are missing. Enter all requested information
More informationMathematics 136 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 19 and 21, 2016
Mathematics 36 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 9 and 2, 206 Every rational function (quotient of polynomials) can be written as a polynomial
More information(x x 0 )(x x 1 )... (x x n ) (x x 0 ) + y 0.
> 5. Numerical Integration Review of Interpolation Find p n (x) with p n (x j ) = y j, j = 0, 1,,..., n. Solution: p n (x) = y 0 l 0 (x) + y 1 l 1 (x) +... + y n l n (x), l k (x) = n j=1,j k Theorem Let
More informationDRAFT - Math 102 Lecture Note - Dr. Said Algarni
Math02 - Term72 - Guides and Exercises - DRAFT 7 Techniques of Integration A summery for the most important integrals that we have learned so far: 7. Integration by Parts The Product Rule states that if
More informationMethods of Integration
Methods of Integration Professor D. Olles January 8, 04 Substitution The derivative of a composition of functions can be found using the chain rule form d dx [f (g(x))] f (g(x)) g (x) Rewriting the derivative
More informationMathematics 104 Fall Term 2006 Solutions to Final Exam. sin(ln t) dt = e x sin(x) dx.
Mathematics 14 Fall Term 26 Solutions to Final Exam 1. Evaluate sin(ln t) dt. Solution. We first make the substitution t = e x, for which dt = e x. This gives sin(ln t) dt = e x sin(x). To evaluate the
More information8.3 Partial Fraction Decomposition
8.3 partial fraction decomposition 575 8.3 Partial Fraction Decomposition Rational functions (polynomials divided by polynomials) and their integrals play important roles in mathematics and applications,
More informationFast and accurate methods for the discretization of singular integral operators given on surfaces
Fast and accurate methods for the discretization of singular integral operators given on surfaces James Bremer University of California, Davis March 15, 2018 This is joint work with Zydrunas Gimbutas (NIST
More informationSOLUTIONS FOR PRACTICE FINAL EXAM
SOLUTIONS FOR PRACTICE FINAL EXAM ANDREW J. BLUMBERG. Solutions () Short answer questions: (a) State the mean value theorem. Proof. The mean value theorem says that if f is continuous on (a, b) and differentiable
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
More informationMIXED potential integral equation formulations have
1980 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 58, NO. 6, JUNE 2010 Evaluation of Weakly Singular Integrals Via Generalized Cartesian Product Rules Based on the Double Exponential Formula Athanasios
More informationAEA 2003 Extended Solutions
AEA 003 Extended Solutions These extended solutions for Advanced Extension Awards in Mathematics are intended to supplement the original mark schemes, which are available on the Edexcel website. 1. Since
More informationMath 230 Mock Final Exam Detailed Solution
Name: Math 30 Mock Final Exam Detailed Solution Disclaimer: This mock exam is for practice purposes only. No graphing calulators TI-89 is allowed on this test. Be sure that all of your work is shown and
More informationHere are brief notes about topics covered in class on complex numbers, focusing on what is not covered in the textbook.
Phys374, Spring 2008, Prof. Ted Jacobson Department of Physics, University of Maryland Complex numbers version 5/21/08 Here are brief notes about topics covered in class on complex numbers, focusing on
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 39, pp. 1-11, 1. Copyright 1,. ISSN 168-9613. ETNA TRIGONOMETRIC GAUSSIAN QUADRATURE ON SUBINTERVALS OF THE PERIOD GASPARE DA FIES AND MARCO VIANELLO
More informationt 2 + 2t dt = (t + 1) dt + 1 = arctan t x + 6 x(x 3)(x + 2) = A x +
MATH 06 0 Practice Exam #. (0 points) Evaluate the following integrals: (a) (0 points). t +t+7 This is an irreducible quadratic; its denominator can thus be rephrased via completion of the square as a
More informationPURE MATHEMATICS AM 27
AM SYLLABUS (2020) PURE MATHEMATICS AM 27 SYLLABUS 1 Pure Mathematics AM 27 (Available in September ) Syllabus Paper I(3hrs)+Paper II(3hrs) 1. AIMS To prepare students for further studies in Mathematics
More informationFormulas to remember
Complex numbers Let z = x + iy be a complex number The conjugate z = x iy Formulas to remember The real part Re(z) = x = z+z The imaginary part Im(z) = y = z z i The norm z = zz = x + y The reciprocal
More informationApproximation theory
Approximation theory Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 1 1.3 6 8.8 2 3.5 7 10.1 Least 3squares 4.2
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationLearning Objectives for Math 166
Learning Objectives for Math 166 Chapter 6 Applications of Definite Integrals Section 6.1: Volumes Using Cross-Sections Draw and label both 2-dimensional perspectives and 3-dimensional sketches of the
More informationFact: Every matrix transformation is a linear transformation, and vice versa.
Linear Transformations Definition: A transformation (or mapping) T is linear if: (i) T (u + v) = T (u) + T (v) for all u, v in the domain of T ; (ii) T (cu) = ct (u) for all scalars c and all u in the
More informationThe iterated sinh transformation
The iterated sinh transformation Author Elliott, David, Johnston, Peter Published 2008 Journal Title International Journal for Numerical Methods in Engineering DOI https://doi.org/10.1002/nme.2244 Copyright
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 informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationSET 1. (1) Solve for x: (a) e 2x = 5 3x
() Solve for x: (a) e x = 5 3x SET We take natural log on both sides: ln(e x ) = ln(5 3x ) x = 3 x ln(5) Now we take log base on both sides: log ( x ) = log (3 x ln 5) x = log (3 x ) + log (ln(5)) x x
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. M. A. and M.Sc. Scholarship Test. September 22, Time Allowed: 150 Minutes Maximum Marks: 30
NATIONAL BOARD FOR HIGHER MATHEMATICS M A and MSc Scholarship Test September 22, 2018 Time Allowed: 150 Minutes Maximum Marks: 30 Please read, carefully, the instructions that follow INSTRUCTIONS TO CANDIDATES
More informationHIGH ORDER METHODS FOR WEAKLY SINGULAR INTEGRAL EQUATIONS WITH NONSMOOTH INPUT FUNCTIONS
MATHEMATICS OF COMPUTATION Volume 67, Number 224, October 1998, Pages 1493 1515 S 0025-5718(98)01005-9 HIGH ORDER METHODS FOR WEAKLY SINGULAR INTEGRAL EQUATIONS WITH NONSMOOTH INPUT FUNCTIONS G. MONEGATO
More informationElliptic Problems for Pseudo Differential Equations in a Polyhedral Cone
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 9, Number 2, pp. 227 237 (2014) http://campus.mst.edu/adsa Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone
More informationRiemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E,
Tel Aviv University, 26 Analysis-III 9 9 Improper integral 9a Introduction....................... 9 9b Positive integrands................... 9c Special functions gamma and beta......... 4 9d Change of
More informationHybrid of Sigmoidal Transformations and Collocation Method for a Generalized Airfoil Equation
Journal of Numerical Mathematics and Stochastics, 3(1) : 80-95, 2011 http://www.jnmas.org/jnmas3-9.pdf JNM@S Euclidean Press, LLC Online: ISSN 2151-2302 Hybrid of Sigmoidal Transformations and Collocation
More informationConvergence and Stability of a New Quadrature Rule for Evaluating Hilbert Transform
Convergence and Stability of a New Quadrature Rule for Evaluating Hilbert Transform Maria Rosaria Capobianco CNR - National Research Council of Italy Institute for Computational Applications Mauro Picone,
More informationZERO DISTRIBUTION OF POLYNOMIALS ORTHOGONAL ON THE RADIAL RAYS IN THE COMPLEX PLANE* G. V. Milovanović, P. M. Rajković and Z. M.
FACTA UNIVERSITATIS (NIŠ) Ser. Math. Inform. 12 (1997), 127 142 ZERO DISTRIBUTION OF POLYNOMIALS ORTHOGONAL ON THE RADIAL RAYS IN THE COMPLEX PLANE* G. V. Milovanović, P. M. Rajković and Z. M. Marjanović
More informationFinal Exam 2011 Winter Term 2 Solutions
. (a Find the radius of convergence of the series: ( k k+ x k. Solution: Using the Ratio Test, we get: L = lim a k+ a k = lim ( k+ k+ x k+ ( k k+ x k = lim x = x. Note that the series converges for L
More informationPreliminary algebra. Polynomial equations. and three real roots altogether. Continue an investigation of its properties as follows.
978-0-51-67973- - Student Solutions Manual for Mathematical Methods for Physics and Engineering: 1 Preliminary algebra Polynomial equations 1.1 It can be shown that the polynomial g(x) =4x 3 +3x 6x 1 has
More informationy = x 3 and y = 2x 2 x. 2x 2 x = x 3 x 3 2x 2 + x = 0 x(x 2 2x + 1) = 0 x(x 1) 2 = 0 x = 0 and x = (x 3 (2x 2 x)) dx
Millersville University Name Answer Key Mathematics Department MATH 2, Calculus II, Final Examination May 4, 2, 8:AM-:AM Please answer the following questions. Your answers will be evaluated on their correctness,
More information(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2.
1. Complex numbers A complex number z is defined as an ordered pair z = (x, y), where x and y are a pair of real numbers. In usual notation, we write z = x + iy, where i is a symbol. The operations of
More informationDerivative and Integral Rules These are on the inside of the back cover of your text.
Derivative and Integral Rules These are on the inside of the back cover of your text. General Derivative Rule General Integral Rule d dx u(x) r = r u(x) r - 1 u(x) u(x)r u(x) dx = u(x) r1 r1 + C r U -1
More information10550 PRACTICE FINAL EXAM SOLUTIONS. x 2 4. x 2 x 2 5x +6 = lim x +2. x 2 x 3 = 4 1 = 4.
55 PRACTICE FINAL EXAM SOLUTIONS. First notice that x 2 4 x 2x + 2 x 2 5x +6 x 2x. This function is undefined at x 2. Since, in the it as x 2, we only care about what happens near x 2 an for x less than
More informationMA Spring 2013 Lecture Topics
LECTURE 1 Chapter 12.1 Coordinate Systems Chapter 12.2 Vectors MA 16200 Spring 2013 Lecture Topics Let a,b,c,d be constants. 1. Describe a right hand rectangular coordinate system. Plot point (a,b,c) inn
More informationPower Series. x n. Using the ratio test. n n + 1. x n+1 n 3. = lim x. lim n + 1. = 1 < x < 1. Then r = 1 and I = ( 1, 1) ( 1) n 1 x n.
.8 Power Series. n x n x n n Using the ratio test. lim x n+ n n + lim x n n + so r and I (, ). By the ratio test. n Then r and I (, ). n x < ( ) n x n < x < n lim x n+ n (n + ) x n lim xn n (n + ) x
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationNumerische Mathematik
Numer. Math. (28) 9:43 65 DOI.7/s2-7-25-7 Numerische Mathematik The superconvergence of Newton Cotes rules for the Hadamard finite-part integral on an interval Jiming Wu Weiwei Sun Received: 2 May 27 /
More informationarxiv:math.ca/ v2 17 Jul 2000
NECESSARY AND SUFFICIENT CONDITIONS FOR DIFFERENTIABILITY OF A FUNCTION OF SEVERAL VARIABLES. R.P. VENKATARAMAN, #1371, 13'TH MAIN ROAD, II STAGE, FIRST PHASE, B.T.M. LAYOUT, BANGALORE 560 076.
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 9, No. 3, pp. 07 6 c 007 Society for Industrial and Applied Mathematics NUMERICAL QUADRATURES FOR SINGULAR AND HYPERSINGULAR INTEGRALS IN BOUNDARY ELEMENT METHODS MICHAEL CARLEY
More informationHOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION)
HOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION) PROFESSOR STEVEN MILLER: BROWN UNIVERSITY: SPRING 2007 1. CHAPTER 1: MATRICES AND GAUSSIAN ELIMINATION Page 9, # 3: Describe
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationCalculus II Practice Test Problems for Chapter 7 Page 1 of 6
Calculus II Practice Test Problems for Chapter 7 Page of 6 This is a set of practice test problems for Chapter 7. This is in no way an inclusive set of problems there can be other types of problems on
More informationSolutions to Calculus problems. b k A = limsup n. a n limsup b k,
Solutions to Calculus problems. We will prove the statement. Let {a n } be the sequence and consider limsupa n = A [, ]. If A = ±, we can easily find a nondecreasing (nonincreasing) subsequence of {a n
More informationSection Taylor and Maclaurin Series
Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power
More informationMATHEMATICAL FORMULAS AND INTEGRALS
MATHEMATICAL FORMULAS AND INTEGRALS ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom Academic Press San Diego New York Boston London
More informationInfinite Series. 1 Introduction. 2 General discussion on convergence
Infinite Series 1 Introduction I will only cover a few topics in this lecture, choosing to discuss those which I have used over the years. The text covers substantially more material and is available for
More informationf(g(x)) g (x) dx = f(u) du.
1. Techniques of Integration Section 8-IT 1.1. Basic integration formulas. Integration is more difficult than derivation. The derivative of every rational function or trigonometric function is another
More informationExperimental mathematics and integration
Experimental mathematics and integration David H. Bailey http://www.davidhbailey.com Lawrence Berkeley National Laboratory (retired) Computer Science Department, University of California, Davis October
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More information1. COMPLEX NUMBERS. z 1 + z 2 := (a 1 + a 2 ) + i(b 1 + b 2 ); Multiplication by;
1. COMPLEX NUMBERS Notations: N the set of the natural numbers, Z the set of the integers, R the set of real numbers, Q := the set of the rational numbers. Given a quadratic equation ax 2 + bx + c = 0,
More informationOn the Chebyshev quadrature rule of finite part integrals
Journal of Computational and Applied Mathematics 18 25) 147 159 www.elsevier.com/locate/cam On the Chebyshev quadrature rule of finite part integrals Philsu Kim a,,1, Soyoung Ahn b, U. Jin Choi b a Major
More informationCandidates are expected to have available a calculator. Only division by (x + a) or (x a) will be required.
Revision Checklist Unit C2: Core Mathematics 2 Unit description Algebra and functions; coordinate geometry in the (x, y) plane; sequences and series; trigonometry; exponentials and logarithms; differentiation;
More informationGaussian interval quadrature rule for exponential weights
Gaussian interval quadrature rule for exponential weights Aleksandar S. Cvetković, a, Gradimir V. Milovanović b a Department of Mathematics, Faculty of Mechanical Engineering, University of Belgrade, Kraljice
More information1. Introduction. 2. Outlines
1. Introduction Graphs are beneficial because they summarize and display information in a manner that is easy for most people to comprehend. Graphs are used in many academic disciplines, including math,
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationReview of Integration Techniques
A P P E N D I X D Brief Review of Integration Techniques u-substitution The basic idea underlying u-substitution is to perform a simple substitution that converts the intergral into a recognizable form
More informationMath Real Analysis II
Math 432 - Real Analysis II Solutions to Homework due February 3 In class, we learned that the n-th remainder for a smooth function f(x) defined on some open interval containing is given by f (k) () R
More informationOR MSc Maths Revision Course
OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision
More informationWeek 2 Techniques of Integration
Week Techniques of Integration Richard Earl Mathematical Institute, Oxford, OX LB, October Abstract Integration by Parts. Substitution. Rational Functions. Partial Fractions. Trigonometric Substitutions.
More informationChapter 2: Differentiation
Chapter 2: Differentiation Winter 2016 Department of Mathematics Hong Kong Baptist University 1 / 75 2.1 Tangent Lines and Their Slopes This section deals with the problem of finding a straight line L
More informationThe Hardy-Littlewood Function: An Exercise in Slowly Convergent Series
The Hardy-Littlewood Function: An Exercise in Slowly Convergent Series Walter Gautschi Department of Computer Sciences Purdue University West Lafayette, IN 4797-66 U.S.A. Dedicated to Olav Njåstad on the
More information8.6 Partial Fraction Decomposition
628 Systems of Equations and Matrices 8.6 Partial Fraction Decomposition This section uses systems of linear equations to rewrite rational functions in a form more palatable to Calculus students. In College
More informationCourse Notes for Calculus , Spring 2015
Course Notes for Calculus 110.109, Spring 2015 Nishanth Gudapati In the previous course (Calculus 110.108) we introduced the notion of integration and a few basic techniques of integration like substitution
More informationPartial Fractions. June 27, In this section, we will learn to integrate another class of functions: the rational functions.
Partial Fractions June 7, 04 In this section, we will learn to integrate another class of functions: the rational functions. Definition. A rational function is a fraction of two polynomials. For example,
More informationn=0 ( 1)n /(n + 1) converges, but not
Math 07H Topics for the third exam (and beyond) (Technically, everything covered on the first two exams plus...) Absolute convergence and alternating series A series a n converges absolutely if a n converges.
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationTHE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK. Summer Examination 2009.
OLLSCOIL NA héireann, CORCAIGH THE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK Summer Examination 2009 First Engineering MA008 Calculus and Linear Algebra
More information(x 3)(x + 5) = (x 3)(x 1) = x + 5
RMT 3 Calculus Test olutions February, 3. Answer: olution: Note that + 5 + 3. Answer: 3 3) + 5) = 3) ) = + 5. + 5 3 = 3 + 5 3 =. olution: We have that f) = b and f ) = ) + b = b + 8. etting these equal
More informationOn the Lebesgue constant of subperiodic trigonometric interpolation
On the Lebesgue constant of subperiodic trigonometric interpolation Gaspare Da Fies and Marco Vianello November 4, 202 Abstract We solve a recent conjecture, proving that the Lebesgue constant of Chebyshev-like
More informationMath 23b Practice Final Summer 2011
Math 2b Practice Final Summer 211 1. (1 points) Sketch or describe the region of integration for 1 x y and interchange the order to dy dx dz. f(x, y, z) dz dy dx Solution. 1 1 x z z f(x, y, z) dy dx dz
More informationReal Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi
Real Analysis Math 3AH Rudin, Chapter # Dominique Abdi.. If r is rational (r 0) and x is irrational, prove that r + x and rx are irrational. Solution. Assume the contrary, that r+x and rx are rational.
More informationA PLANAR SOBOLEV EXTENSION THEOREM FOR PIECEWISE LINEAR HOMEOMORPHISMS
A PLANAR SOBOLEV EXTENSION THEOREM FOR PIECEWISE LINEAR HOMEOMORPHISMS EMANUELA RADICI Abstract. We prove that a planar piecewise linear homeomorphism ϕ defined on the boundary of the square can be extended
More informationPublic-key Cryptography: Theory and Practice
Public-key Cryptography Theory and Practice Department of Computer Science and Engineering Indian Institute of Technology Kharagpur Chapter 2: Mathematical Concepts Divisibility Congruence Quadratic Residues
More information2016 EF Exam Texas A&M High School Students Contest Solutions October 22, 2016
6 EF Exam Texas A&M High School Students Contest Solutions October, 6. Assume that p and q are real numbers such that the polynomial x + is divisible by x + px + q. Find q. p Answer Solution (without knowledge
More informationFinal Examination Solutions
Math. 5, Sections 5 53 (Fulling) 7 December Final Examination Solutions Test Forms A and B were the same except for the order of the multiple-choice responses. This key is based on Form A. Name: Section:
More informationSECTION A. f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes.
SECTION A 1. State the maximal domain and range of the function f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes. 2. By evaluating f(0),
More informationChapter 2: Differentiation
Chapter 2: Differentiation Spring 2018 Department of Mathematics Hong Kong Baptist University 1 / 82 2.1 Tangent Lines and Their Slopes This section deals with the problem of finding a straight line L
More informationChapter 8: Techniques of Integration
Chapter 8: Techniques of Integration Section 8.1 Integral Tables and Review a. Important Integrals b. Example c. Integral Tables Section 8.2 Integration by Parts a. Formulas for Integration by Parts b.
More informationarxiv: v1 [math.ca] 31 Dec 2018
arxiv:181.1173v1 [math.ca] 31 Dec 18 Some trigonometric integrals and the Fourier transform of a spherically symmetric exponential function Hideshi YAMANE Department of Mathematical Sciences, Kwansei Gakuin
More informationMaterial for review. By Lei. May, 2011
Material for review. By Lei. May, 20 You shouldn t only use this to do the review. Read your book and do the example problems. Do the problems in Midterms and homework once again to have a review. Some
More informationLINEAR SECOND-ORDER EQUATIONS
LINEAR SECON-ORER EQUATIONS Classification In two independent variables x and y, the general form is Au xx + 2Bu xy + Cu yy + u x + Eu y + Fu + G = 0. The coefficients are continuous functions of (x, y)
More informationarxiv: v1 [math.na] 14 Feb 2019
Closed-form evaluation of potential integrals in the Boundary Element Method Michael Carley (m.j.carley@bath.ac.uk) arxiv:192.551v1 [math.na] 14 Feb 219 February 15, 219 Abstract A method is presented
More informationStudy Guide/Practice Exam 3
Study Guide/Practice Exam 3 This study guide/practice exam covers only the material since exam. The final exam, however, is cumulative so you should be sure to thoroughly study earlier material. The distribution
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationGRE Math Subject Test #5 Solutions.
GRE Math Subject Test #5 Solutions. 1. E (Calculus) Apply L Hôpital s Rule two times: cos(3x) 1 3 sin(3x) 9 cos(3x) lim x 0 = lim x 2 x 0 = lim 2x x 0 = 9. 2 2 2. C (Geometry) Note that a line segment
More information