Computation of the Eigenvalues of 3 D “ Charged ” Integral Equations

The Rayleigh-Ritz and the inverse iteration methods are used in order to compute the eigenvalues of 3D Fredholm-Stieltjes integral equations, i.e. 3D Fredholm equations with respect to suitable Stieltjes-type measures. Some applications are shown, relevant to the problem of computing the eigenvalues of a body charged by a finite number of masses concentrated on points, curves or surfaces lying in.


Introduction
The theory of Fredholm integral equations is strictly connected with the birth of functional analysis.A background of this theory can be found in classical books (see e.g.[1] [2]).For recent developments relevant to numerical computation of solutions, see [3].
In [4] [5], the Author shows that the Fredholm theory still holds considering the so called charged Fredholm integral equations, i.e.Fredholm equations with respect to a Stieltjes measure obtained by adding to the ordinary Lebesgue measure a finite combination of positive masses concentrated in arbitrary points of the considered interval.A topic which at present is included, as a particular case, in the theory of strictly positive compact operators.
In the one-dimensional case, the mechanical interpretation of these equations is connected with the problem of the free vibrations of a string charged by a finite number of cursors, and is related to an extension of the classical orthogonality property of eigensolutions, called the "Sobolev-type" orthogonality (see e.g.[6] [7]).
In preceding articles [8] [9] [10], the problem of numerical computation of the above mentioned eigenvalue problems was solved, by using the so called inverse iteration method, showing applications in one and two-dimensional cases.
In this article, after briefly recalling the main results about the theory of eigenvalues for charged Fredholm integral equations, we mention how to obtain (in the particular case of a symmetric and strictly positive operator), the lower and upper approximations of these eigenvalues by means of the Rayleigh-Ritz method [1], [11] and the Fichera orthogonal invariants method [11] [12] respectively.Then we conclude by showing that, even in the three-dimensional case, the lower approximations of the eigenvalues obtained by means of the Rayleigh-Ritz method can be improved by applying the inverse iteration method [13].Numerical computations relevant to the considered case are developed in the concluding section.

Lebesgue-Stieltjes Measures in a 3D Interval
Consider the interval  and the Hilbert space ( ) where [ ] ( ) 2 , , dM L dM U V U V ≡ and the subscripts " ( )d refer to the following definitions: (obviously ), so that the Stieltjes measure is the sum of the ordinary Lebesgue measure plus a finite sum of charges h m concentrated on It is worth noting that ( ) L  , a (complete) Hilbert space, which does not have singularities (i.e.

discontinuities) at points
and surfaces Σ  , ( ) , according to the usual condition for existence of the relevant Stieltjes integrals (see e.g.[14]).
Let us consider, for example, the computation of integrals with respect to the above Dirac-type measures: These formulas will be useful in the following.

The Eigenvalue Problem for 3D "Charged" Operators
Consider in ( ) , K P Q is a symmetric kernel, and ( δ denoting the usual Dirac-Delta function).According to the above positions, we have where we used the symmetry of the kernel, where we used the symmetry of the kernel, ( where we used the symmetry of the kernel, where we used the symmetry of the kernel, ( where we used the symmetry of the kernel, where we used the symmetry of the kernel, and lastly

Computation of the Eigenvalues for Charged Integral Equations
The computation of the eigenvalues of second kind Fredholm integral equations is usually performed by using the Rayleig-Ritz method [1] [11] for lower bounds, and the Fichera orthogonal invariants method [11] [12] for upper bounds.An alternative procedure, called the inverse iteration method, can be used to improve the lower approximations previously obtained by means of the Rayleigh-Ritz method.This approach was already considered in [13], and will be applied even in the present case.
We will not describe herewith the three above mentioned methods, because they are essentially independent of the dimension of the considered vibrating item (string, membrane or body).We refer for shortness to the above mentioned articles [8] [9] [11] [12].

Applications
Let , if .
The considered operator is compact and strictly positive, since it is connected with free vibrations of a body charged by a finite number of masses concentrated on points, curves, or surfaces contained in  .
, , e x y z x y z ρ − − − = (see Figure 5).Under this assumption, the approximate eigenvalues evaluated by using the Rayleigh-Ritz and inverse iteration methods are listed in Table 3, whereas the relevant eigenfunctions are shown in Figure 6.
( ) (see Figure 7).Under this assumption, the approximate eigenvalues evaluated by using the Rayleigh-Ritz and inverse iteration methods are listed in Table 4, whereas the relevant eigenfunctions are shown in Figure 8.   (see Figure 9).Under this assumption, the approximate eigenvalues evaluated by using the Rayleigh-Ritz and inverse iteration methods are listed in Table 5, whereas the relevant eigenfunctions are shown in Figure 10.

Numerical
(see Figure 11).Under this assumption, the approximate eigenvalues evaluated by using the Rayleigh-Ritz and inverse iteration methods are listed in Table 6, whereas the relevant eigenfunctions are shown in Figure 12.
(see Figure 13)..Under this assumption, the approximate eigenvalues evaluated by using the Rayleigh-Ritz and inverse iteration methods are listed in Table 8, whereas the relevant eigenfunctions are shown in Figure 16.
plus a finite sum of continuous charges belonging to the curves Σ ⊂  sum of continuous charges belonging to the surfaces Σ ⊂   , with densities ( )

2 .Figure 1 .Figure 2 .
Figure 1.Spatial distribution of the volume density function ( ) , , x y z ρ relevant to the

Figure 3 .Figure 4 .
Figure 3. Spatial distribution of the volume density function ( ) , , x y z ρ relevant to the et al.DOI: 10.4236/jamp.2017.5101712061 Journal of Applied Mathematics and Physics

Figure 5 .Figure 6 .
Figure 5. Spatial distribution of the volume density ( ) , , x y z ρ (a) and concentrated density functions (b) relevant to the numerical Example 3.

Figure 7 .Figure 8 .
Figure 7. Spatial distribution of the volume density ( ) , , x y z ρ (a) and concentrated density functions (b) relevant to the numerical Example 4.

Figure 9 .
Figure 9. Spatial distribution of the volume density ( ) , , x y z ρ (a) and concentrated density functions (b) relevant to the numerical Example 5.

Figure 11 .Figure 12 .
Figure 11.Spatial distribution of the volume density ( ) , , x y z ρ (a) and concentrated density functions (b) relevant to the numerical Example 6.

Figure 13 .Figure 14 .
Figure 13.Spatial distribution of the volume density ( ) , , x y z ρ (a) and concentrated density functions (b) relevant to the numerical Example 7.

Figure 15 .
Figure 15.Spatial distribution of the volume density ( ) , , x y z ρ (a) and concentrated density functions (b) relevant to the numerical Example 8.

Figure 17 .
Figure 17.Spatial distribution of the volume density ( ) , , x y z ρ (a) and concentrated density functions (b) relevant to the numerical Example 9.

Table 1 .
Approximate eigenvalues h µ  and ˆh µ ( 1, 2, , 6 h =  ) relevant to the numerical Example 1, as computed by means of the Rayleigh-Ritz and inverse iteration method, respectively.

Table 3 .
Approximate eigenvalues h µ  and ˆh µ ( 1, 2, , 6 h =  ) relevant to the numerical Example 3, as computed by means of the Rayleigh-Ritz and inverse iteration method, respectively.

Table 6 .
Approximate eigenvalues h µ  and ˆh µ ( 1, 2, , 6 h =  ) relevant to the numerical Example 6, as computed by means of the Rayleigh-Ritz and inverse iteration method, respectively.

Table 8 .
Approximate eigenvalues h µ  and ˆh µ ( 1, 2, , 6 h =  ) relevant to the numerical Example 8, as computed by means of the Rayleigh-Ritz and inverse iteration method, respectively.

Table 9 .
Approximate eigenvalues h µ  and ˆh µ ( 1, 2, , 6 h =  ) relevant to the numerical Example 9, as computed by means of the Rayleigh-Ritz and inverse iteration method, respectively.at The Antenna Company, and is funded by the Competitiveness Enhancement Program grant at Tomsk Polytechnic University.