Method of Infinite System of Equations for Problems in Unbounded Domains

Many problems of mechanics and physics are posed in unbounded or infinite domains. For solving these problems one typically limits them to bounded domains and find ways to set appropriate conditions on artificial boundaries or use quasi-uniform grid that maps unbounded domains to bounded ones. Differently from the above methods we approach to problems in unbounded domains by infinite system of equations. In this paper we present starting results in this approach for some one-dimensional problems. The problems are reduced to infinite system of linear equations. A method for obtaining approximate solution with a given accuracy is proposed. Numerical experiments for several examples show the effectiveness of the offered method.


Introduction
Many problems of mechanics and physics are posed in unbounded or infinite domains, for example, heat transport problems in infinite or semi-infinite bar, aerosol propagation in atmosphere, problem of ocean pollution, wave propagation in unbounded media, and problem of computing the potential distribution due to a source of current in or on the surface of the Earth.For solving these problems one usually restricts oneself to treat the problem in a bounded domain and try to use available efficient methods for finding exact or approximate solutions in the restricted domain.But there arise some questions: which size of restricted domain is enough and how to set conditions on artificial boundary for obtaining approximate solution with good accuracy?The simplest way to do this is to transfer boundary condition on infinity to the artificial boundary.This raw way may lead to large deviation of approximate solution from the solution of the original problem.Therefore, instead of transferring boundary condition on infinity without changes one tries to set appropriate conditions on artificial boundary.This is a direction of researches that attracts the attention of many specialists in the fields of mathematics, mechanics, and physics see 1-4 .[16] Infinite system of linear equations with infinite number of unknowns is the system of the form

Infinite System of Equations
A sequence of numbers x 1 , x 2 , . . . is called a solution of the system 2.1 if after the substitution of these numbers into the right hand sides we obtain convergent series and all equalities are satisfied.Suppose that the system has solutions.Then the solution found by the method of successive approximation or simple iteration with the zero starting approximation x 0 i 0 i 1, 2, . . . is called the main solution of the system.
We say that a solution x 1 , x 2 , . . . of the system tends to zero if x i → 0 as i → ∞.We are interested in the system for which

2.3
The system 2.1 is called regular system if it is called completely regular if there exists a constant θ > 0 such that Assume that there exists a number K such that the free members b i satisfy the condition

Quasi-Uniform Grid [8, 10]
Let x ξ be strictly monotone smooth function of the argument ξ ∈ 0, 1 .The grid ω N x i x i N , i 0, 1, . . ., N , 2.9 with x 0 a, x 1 ∞ is called quasi-uniform grid on a, ∞ .In this case the last node x N of the grid is on the infinity.

Journal of Applied Mathematics
Example of quasi-uniform grids are the grids , i 0, 1, . . ., N tangential grid . 2.10

Method of Infinite System for Stationary Problems
We will present in details the method of infinite system on the model problem of heat conductivity in a semi-infinite bar under the usual assumptions that the functions k x , d x , and f x are continuous and In the case if k, d are constants and f x has a compact support 0, L one can easily find the exact artificial boundary condition at x L with the help of Dirichlet-to-Neumann map 4 .If f x has no compact support but has special form such that it is possible to find a particular solution of the equation then an exact artificial boundary condition can be established.In general case when k, d, and f only satisfy the condition 3.2 one cannot find the exact condition at x L if the problem is restricted to finite interval 0, L .Below we consider this case.
In order to solve the problem 3.1 -3.2 we introduce the uniform grid {x i ih, i 0, 1, . ..} and consider the difference scheme where Now we rewrite the difference scheme in the form of customary three-point difference equations

3.6
Here Putting

3.8
we can rewrite the system 3.6 in the canonical form of infinite system as follows: 3.9 We have ρ 0 1 − p 0 − q 0 1 and from 3.7 , 3.8 it follows 0 < p i , q i < 1;

3.10
Thus, the system 3.9 is regular.More precisely, it is completely regular because it is easy to verify that where 3.12 Now consider r i /ρ i .From 3.8 and 3.10 we have r i /ρ i f i /d i , and from the assumptions 3.2 it follows that f i /d i → 0. Therefore, there exists a constant K * such that |f i | ≤ K * d i for any i.Consequently, the conditions of Theorem 2.3 are satisfied and the solution of the infinite system 3.8 can be found by the truncation method.
A question that arises here is to which size do we need to truncate the infinite system for obtaining approximate solution with a given accuracy?Below we give answer to this question.
Following the progonka method or Thomas algorithm which is a special form of the Gauss elimination 17 for tridiagonal system of equations we shall seek the solution of 3.9 in the form where coefficients are calculated by the formulas: For the coefficients α i and β i there hold the following assertions. ii Proof.i Analogously as for a finite system of three-point difference equations see 17 it is possible to prove by induction that 0 < α i < 1 i 2, 3, . . . .Hence, from the condition y i → 0 and from 3.13 it follows that β i → 0 as i → ∞. ii Since 0 < α i < 1 is proved, from 3.10 and 3.11 we have Further, taking into account 3.8 and 3.9 we obtain Then there holds the estimate 0 < α i < μ < 1, i 2, 3, . . . .

3.19
The proof of the lemma is complete.
Theorem 3.2.Given an accuracy ε > 0. If starting from a natural number N there holds then for the deviation of the solution of the truncated system

3.21
compared with the solution of the infinite system 3.9 there holds the following estimate Proof.Denote z i y i − y i .Then it is easy to see that z i satisfies the infinite system 23 where b i 0, i 0, . . ., N,
Remark 3.3.The above theorem permits us in the process of computation of the sweep coefficients 3.14 to determine when to truncate the infinite system 3.9 for guarantee that the solution of the infinite system less than given ε.
Below in order to illustrate the effectiveness of using Theorem 3.2 we consider the following.

3.26
This problem has the exact solution u x 1/ 1 x 2 .Construct the infinite system 3.9 and truncate it when Theorem 3.2 is satisfied.The solution of the truncated system is compared with the exact solution.The results of computation on the uniform grid with step h 0.1 and h 0.05 are given in the Tables 1 and 2, where N is the size of the system, that is, automatically truncated, error max 0 i N y i − u i , u i u x i .

3.27
The graphs of the exact and approximate solutions are given in Figure 1, meanwhile the graph of the right-hand side is given in Figure 2.  In the process of computation we observe that the coefficients β i tend to zero very fast and the coefficients α i have the trend of approaching to 1 but the ratio β i / 1 − α i tend to zero rather fast, too.The graphs of the coefficients α i , β i , and the ratio β i / 1 − α i are given in Figures 3, 4, and 5, respectively.Remark 3.5.The problem in semi-infinite bar for the more general than 3.1 equation − ku ru du f x , x > 0, is treated in similar way if for it we construct according to Samarskii 17 a monotone difference scheme.Remark 3.6.In the case if instead of Dirichlet boundary condition at left end-point there is given Neumann boundary condition then the problem can be treated in similar way.Remark 3.7.If the functions k x , d x , and f x in 3.1 are even functions then the problem on the whole real axis with vanishing conditions, that is, u ±∞ 0, is reduced to the problem on 0, ∞ with the boundary condition u 0 0.

Parabolic Equation on Semi-Infinite Bar
In this section we apply the infinite system technique proposed in the previous section to an initial-boundary problem for parabolic equation.a First we consider the heat conductivity problem with constant coefficient k > 0

4.1
This problem has the exact solution Using purely implicit difference scheme on uniform grid with spatial stepsize h and time stepsize τ we reduce the problem 4.1 to infinite system on each time layer j 1: −ry where r kτ/h 2 , i, j are indexes of grid nodes in space and in time.
The system 4.3 is treated in a similar way as for the system 3.6 .It is easy to show that at each time layer j 1 for the sweep coefficients α i and β i we have 0 < α i < r/ 1 r and β i → 0. So, we can use Theorem 3.2 to truncate the infinite system.
In order to demonstrate the advantage of the infinite system method over the quasiuniform grid method used in 10, 12 we perform computations by these two methods: infinite system on uniform grid and finite system on quasi-uniform grid x i i/N − i, i 0, 1, . . ., N with N 50.Since the density of quasi-uniform nodes is very sparse for i ≥ 25 the profiles of the approximate solution are broken lines.Figures 6 and 7 give the profiles of solution computed by the two methods with k 10, ε 0.001.From these figures we see that the infinite system method gives better results.
b The stationary problem of air pollution caused by a point source with constant power Q located at point 0, H is reduced to the problem where ϕ is the concentration of aerosol, u is the wind velocity in the x-direction, w g is the falling velocity of pollutants by gravity, σ ≥ 0 is the coefficient of transformation, νthe vertical diffusion coefficient, and α ≥ 0-the coefficient characterizing the reflection and absorption of the bedding surface 18, 19 .
The numerical solution of the above problem on uniform grid using the infinite system method was studied in 18 , where a theorem similar to Theorem 3.2 was proved.
An interesting fact was established in 20 .It is that if limiting the problem of air pollution to the domain 0 ≤ z ≤ Z with finite Z then artificial boundary condition ϕ x, Z 0 causes an undershoot and the Neumann condition ∂ϕ/∂z x, Z 0 causes an overshoot in numerical solution when compared with the solution of the problem with the boundary condition ϕ x, ∞ 0.

The Equation of Complex Type
The advantage of the infinite system method in comparison with the quasi-uniform grid method is revealed more clearly in the problems simulating wave phenomena, when high  density of grid nodes is near left endpoint and nodes are gradually sparse in the right cannot well describe waves running to the right.To show this fact we consider a problem describing ion wave in stratified incompressible fluid x, 0 0.  For solving this problem, as in 9 , we set φ ∂ 2 u/∂x 2 − u.Then the problem is decomposed into two consecutive problems of second order

5.2
On the uniform grid with step sizes h, τ replace the above differential problem by the difference schemes  From the explicit scheme 5.3 it is easy to calculate φ j 1 .After that we use the infinite system technique for the scheme 5.4 for computing u j 1 .This algorithm is applied to some examples in comparison with the quasi-uniform grid method.
Example 5.1.In this example we take initial conditions be homogeneous and the left boundary condition u 0, t arctan 2 10t • sin 0.3t .
The profiles u x, t with t 1, 2, . . ., 10 computed by the infinite system method for h 0.1 are given in Figure 8 and by the quasi-uniform grid method are given in Figure 9.  Example 5.2.The left boundary condition is zero, the initial condition is u x, 0 x 3 e −x , ∂u ∂t x, 0 0.

5.5
The propagation of the initial profile after the time t 4, 6, 8, 10, and 12 computed by the infinite system method and by the quasi-uniform grid method are given in Figures 10 and 11, respectively.
From  we see that the approximate solution computed by the method of infinite system describes the wave propagation better than the solution obtained by the method of quasi-uniform grid.

Concluding Remarks
In this paper we propose and investigate the infinite system method for solving several one-dimensional stationary and nonstationary problems, where the keystone is to determine when truncate the infinite system for assuring to obtain approximate solution with a given accuracy.This method reveals advantage over the quasi-uniform grid method, proposed in 2001 by Russian mathematicians in time-dependent problems, especially in problems of wave propagation.
In combination with the alternating directions method the proposed method can be applied to two-dimensional problems in semi-infinite and infinite strips.
The development of the method for solving other two-dimensional and threedimensional problems is the direction of our research in the future.

1 Figure 2 :
Figure 2: Function of right hand side.

Figure 6 :
Figure 6: u x, t j by infinite system.