Numerical solution of an integral equation arising in the problem of cruciform crack using Daubechies scale function

This paper is concerned with obtaining approximate numerical solution of a classical integral equation of some special type arising in the problem of cruciform crack. This integral equation has been solved earlier by various methods in the literature. Here, approximation in terms of Daubechies scale function is employed. The numerical results for stress intensity factor obtained by this method for a specific forcing term are compared to those obtained by various methods available in the literature, and the present method appears to be quite accurate.


Introduction
Integral equations occur naturally in many areas of mathematical physics. Many engineering and applied science problems arising in water waves, potential theory and electrostatics are reduced to solving integral equations. The problem of finding out the crack energy and distribution of stress in the vicinity of a cruciform crack leads to the integral equation where This is an integral equation of some special type since the kernel L(x, t) has singularity at (0, 0) only. f(x) is a prescribed function relating to the internal pressure given by Since the cracks are in the shape of a cross, the problem is known as the cruciform crack problem. Of interest here, is the stress intensity factor (1) which is directly proportional to stress intensity at the crack trip. How the integral equation (1.1) occurs in the problem of cruciform crack is explained by Stallybrass [8] who solved the integral equation in a closed form using Wiener-Hopf technique and provided the numerical results for the stress intensity factor. Rooke and Sneddon [6] solved this integral equation approximately by using an expansion in terms of Legendre functions and obtained numerical results which are very close to those of Stallybrass [8] although the convergence is slow. The two methods appear to be somewhat elaborate. The integral equation (1.1) has also been solved numerically by various other methods from time to time. For example, Elliot [2] employed the method of Sigmoidal transformation to obtain approximate solution for the case f (x) = 1 . It is not obvious if this method is useful for other forms of f(x). Tang and Li [9] solved the integral equation approximately by employing Taylor series expansion for the unknown function and obtained very accurate numerical estimates for the stress x −1 .

3
intensity factor. They made use of Cramer's rule in the mathematical analysis so that if one increases the number of terms in the approximation the calculation becomes unwieldy so as to make the method unattractive. Bhattacharya and Mandal [1] solved the integral equation approximately by two different methods, one is based on expansion of the unknown function in terms of Bernstein polynomials and the other is based on expansion in terms of rationalized Haar functions. Singh and Mandal [7] also solved it by using Legendre multi-wavelets. All these methods provide numerical results for (1) which are very close to exact results given by Stallybrass [8]. Expansion in terms of Bernstein polynomials or Haar functions or Legendre multi-wavelets suggest expansion in terms of other functions such as Daubechies scale functions since these provide a somewhat new tool in the numerical solution of integral equations.
In this paper, Daubechies scale functions are employed to expand the unknown function (x) . K-Daubechies scale function is employed to find approximate solution of integral equation taking K = 3 . It may be noted that K = 1 corresponds to Haar wavelets. As the result can be improved taking larger value of K, so the results obtained by using K-Daubechies scale function are better than the results using the rationalized Haar functions. Though Legendre multi-wavelets give satisfactory results, K-Daubechies scale function has some interesting features like compact support, fractal nature and no explicit form at all resolutions. Only the knowledge of the low-pass filter coefficients in two-scale relation is required throughout the calculation. For these reasons, Daubechies scale function is used as an efficient and new mathematical tool to solve integral equations. At x = 1 , the expansion of (x) reduces to a finite expansion because most of Daubechies scale functions vanish. Actually, the integral equation (1.1) produces a system of linear equations in the unknown coefficients. After solving this linear system, the unknown function (x) is evaluated at x = 1 so as to obtain numerically the value of the stress intensity factor. For different values of in the expression of internal pressure f(x) given by (1.3), (1) is obtained and compared to known results available in the literature. It is found that the method is quite accurate as the approximate values of (1) obtained by the present method are seen to differ negligibly from exact values.

Basic properties of Daubechies scale function and wavelets
Daubechies discovered a whole new class of compactly supported orthogonal wavelets, which is generated from a single function (x) , known as Daubechies scale or refinable function. K-Daubechies scale function (K ≥ 1) has 2K scaling coefficients and has compact support [0, 2K − 1] . It may be noted that K = 1 corresponds to the Haar wavelets.
Using the explicit form of f(x) in (1.3) and the two-scale relation (2.1), the expression in (3.4) reduces to the form Now, using the Gauss-type quadrature rule with complex nodes and weights for integrals involving Daubechies scale function (cf. Panja and Mandal [5]), we obtain where The determination of the nodes x i and weights w i is described by Panja and Mandal [5].
The basic trick for the calculation of the integral (3.5) is described by Kessler et al. [3] and Panja and Mandal [4]. If (x) is the scale function with compact support [0, 2K − 1] (K ≥ 1) , then it produces a system of orthonormal basis sn given by (2.5). From    [4]. Now, the calculation of (3.6) is described. Using the twoscale relation (2.1), from (3.6), we obtain (3.13) (3.14) where I s m,n is given by the relation    If −(2K − 2) ≤ m, n ≤ 0 , then Θ m,n,l 1 ,l 2 (x, t) in (3.17 ) has singularity at (0, 0) , and for these values of m, n, the values of I s m,n cannot be determined using the relation (3.18). If −(2K − 2) ≤ m, n ≤ 0 , the recursion relation for I s m,n is obtained as Again using the Gauss quadrature rule involving the Daubechies scale function, (3.19) is reduced to the form Here, h l j for j = 1, 2 (l j = 0, 1, 2 … , 2K − 1) are the low-pass filters. Basic trick for calculating weights w with a program in MATHEMAT-ICA has been discussed by Panja and Mandal [5]. Also, I s m,n = 0 for m or n ≤ −(2K − 1) or m or n ≥ 2 s . We present here the numerical values of I s m,n for( K = 3)-Daubechies scale functions taking s = 3 for those values of m and n for which Θ m,n,l 1 ,l 2 (x, t) has singularity at (0, 0). Table 1 shows the values of I s m,n for N = 5, whereas Table 2 shows the values of I s m,n for N = 7. (3.20)

Numerical results
A comparison between the numerical values of (1) obtained here by using Daubechies scale functions and exact results of Stallybrass [8] is given in Table 3 for different values of = 1, 2, 3, … , 10.
The table shows the exact values of (1) ( = 1, 2, … , 10) according to Stallybrass [8] and the results obtained by the present method with their relative errors.
The numerical results are displayed in Fig. 1(a-e). For the sake of clarity, five figures are drawn, wherein the stress intensity factor (1) is depicted against the parameter for different integral values. In each figure, (1) obtained from Stallybrass's [8] exact result is denoted by the symbol " ◻ ," and (1) obtained from other approximate methods is shown. The figures are self-explanatory. However, as the result obtained by Sigmoidal transformation method is only for = 1 , this is not shown here. From these figures, it is obvious that all the methods including the present method provide very accurate results.

Conclusion
Here, a numerical scheme based on expansion in terms of K-Daubechies scale function is employed for obtaining approximate numerical estimates of an integral equation of some special type arising in the classical problem of cruciform crack in elasticity. Comparison between the numerical results obtained by the present method with the exact results obtained by Stallybrass [8] shows that the method is quite accurate. The method works nicely for moderate values of K (e.g., K = 3 ). The results can be further improved taking larger values of K (K > 3).