The Solution of Fractal Diffusion Retrospective Problem
1Physics and Mathematics faculty, Penza state university, Penza, Russia
4. Analogues of the Hermite Functions on Piecewise-Homogeneous Real Axis
6. Retrospective Problem for Iteration Heat Equation
7. Power Function with Discontinuous Coefficients and its Application
8. Retrospective Problem for the System of the Diffusion Equations
9. Retrospective Problem for Fractal System of the Diffusion Equations
Abstract
In this article we study the retrospective inverse problem. The retrospective inverse problem consists of in the reconstruction of a priori unknown initial condition of the dynamic system from its known final condition. Existence and uniqueness of the solution is proved.
Keywords: Hermite function, retrospective problem, integral equation, fractal diffusion
Applied Mathematics and Physics, 2013 1 (3),
pp 60-66.
DOI: 10.12691/amp-1-3-3
Received September 15, 2013; Revised September 23, 2013; Accepted September 29, 2013
Copyright © 2013 Science and Education Publishing. All Rights Reserved.Cite this article:
- Yaremko, O., and E. Mogileva. "The Solution of Fractal Diffusion Retrospective Problem." Applied Mathematics and Physics 1.3 (2013): 60-66.
- Yaremko, O. , & Mogileva, E. (2013). The Solution of Fractal Diffusion Retrospective Problem. Applied Mathematics and Physics, 1(3), 60-66.
- Yaremko, O., and E. Mogileva. "The Solution of Fractal Diffusion Retrospective Problem." Applied Mathematics and Physics 1, no. 3 (2013): 60-66.
| Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |
1. Introduction
In this article we study the retrospective inverse problem. The retrospective Inverse problem consists of in the reconstruction of a priori unknown initial condition of the dynamic system from its known final condition. The direct problem of heat conductivity is well-posed; the inverse problem is not well-posed. In mathematics the vast majority of inverse problems set not well-posed - small perturbations of the initial data (observations) can correspond to arbitrarily large perturbations of the solution. The French mathematician Jacques Hadamard in 1939 defined, the problem is called correct or well-posed problem if a solution exists, the solution is unique, the solution’s behavior hardly changes when there’s a slight change in the initial condition. If at least one of these three conditions is not fulfilled, problems are termed ill-posed or not well-posed. The most often in the case of ill-posed problems of the third condition are violated the condition of the stability of solutions. In this case, there is a paradoxical situation: the problem is mathematically generated, but the solution cannot be obtained by conventional methods. A classic example of ill-posed problem is retrospective problem for heat equation on the real axis. Mathematically retrospective problem leads to a Fredholm integral equation of the first kind:
![]() | (1) |
in which
- is the initial distribution of the temperature field,
- is the distribution of the fields in the moment of time t. As shown in [1], the solution of equation (1) expressed by the formula:
![]() | (2) |
2. Problem Statement
In the inverse problem of heat conductivity the initial distribution of sources is unknown. The initial distribution of sources generates the specified temperature distribution in an infinite piecewise-homogeneous rod. In Mathematical statement of the problem consists in finding a solution separatist system (n+1) equations of parabolic type
![]() | (3) |
on the initial conditions
![]() | (4) |
on the boundary conditions
![]() | (5) |
and the coupling conditions
![]() | (6) |
![]() |
here
- unknown function,
- set function,
![]() |
![]() |
- given real number, in which the condition of unlimited solvability of the problem considered fulfilled [2].
The solution to problem (3)-(6) is of the form:
![]() | (7) |
Where
- influence function [2] of the mixed boundary value problem.
Retrospective problem for the heat equation in the case of infinite piecewise-homogeneous rod consists in the determination of the unknown initial distribution of sources
, which generates the specified temperature distribution
in the moment of time
.
3. Transformation Operators
Method of transformation operators is used to solve the problem [2]. Necessary definitions from [2, 11, 12]. The direct
f and inverse
transformation operators are set equalities:
![]() |
![]() |
Here
- are the eigenfunctions [13[13],14][14] of the direct and coupling Sturm–Liouville problems for the Fourier operator in piecewise-homogeneous axis In. Eigenfunction
![]() |
is a solution of the system of separate differential equations
![]() |
on the coupling conditions
![]() |
on the boundary conditions
![]() |
Similarly eigenfunction
![]() |
is a solution of the system of separate differential equations
![]() |
With the coupling conditions
![]() |
where
![]() |
on the boundary conditions
![]() |
Let for some
of the considered boundary value problems have nontrivial solutions
, in this case the number
is called the eigenvalue [13[13],14][14], corresponding solutions
- is called the eigenfunctions of the direct and coupling Sturm–Liouville problems, respectively. In the further we shall adhere to the following normalization of eigenfunctions:
![]() |
4. Analogues of the Hermite Functions on Piecewise-Homogeneous Real Axis
Let define analogues of the system the Hermite functions on piecewise-homogeneous real axis:
![]() |
![]() |
where
- the system of the classical orthogonal Hermite functions [1].
Lemma 1. Functions
form biorthogonal system of functions by piecewise-homogeneous real axis.
Proof. We have the equality:
![]() |
We change the integrals of places, we get:
![]() |
On the decomposition theorem, we have:
![]() |
Consequently,
![]() |
5. Main Result
The problem of determining the initial distribution of the temperature field
mathematically leads to the separate system of integral equations:
![]() | (8) |
![]() |
Method of transformation operators applicable to solving separate system of integral equations (8).
Theorem 1. If the function
and for her the condition
![]() |
that the separate system of integral equations (8) has a unique solution
(definition
[6]), is according to the formula:
![]() |
where
![]() |
Proof . Let’s apply the transformation operator
to separate system of integral equations (8). As a result come to a model integral equation (1). Let’s apply the operator
in both parts of the obtained equality (9); as a result, taking into account the continuity of the operator
, we find the unknown distribution of temperature:
![]() | (9) |
Let’s calculate numbers 
![]() |
from the definition of the operator conversion of
the equality follows:
![]() |
Thus,
![]() |
6. Retrospective Problem for Iteration Heat Equation
6.1. Homogeneous CaseLet
- decision iteration heat equation
![]() |
limited on the set
![]() |
with initial conditions
![]() |
where the designation is accepted:
![]() |
and boundary conditions
![]() |
![]() |
We consider the problem of determining of initial values
of the system of functions
that is on known to the solution and its derivatives up to order
at time 
Let
- solutions of model heat equations:
![]() | (10) |
with initial conditions
![]() | (11) |
boundary conditions
![]() | (12) |
Direct check shows, that the system of functions dependent.
and
.
Lemma 2. The system of functions
and
is connected among themselves by equalities:
![]() | (13) |
- is operator matrix:
![]() |
- is matrix of size
where the diagonal connecting elements
consists of units, and all other elements are equal to zero.
Theorem 2. If the condition
is fulfilled for the system of functions
, then:
![]() | (14) |
and
values
must be taken from the formula (13).
Proof. Let's apply the formula (9) to each of
the initial conditions
in the model problems (10)-(12). As a result we come to (14).
Let the initial distribution is unknown, it generates the specified distribution at time
![]() |
it generates the specified distribution at time
:
![]() |
- is the decision of the separate system of differential equations
![]() | (15) |
![]() |
the specified distribution at time
of the function derivatives
up to order
including:
![]() |
![]() |
initial conditions
![]() | (16) |
boundary conditions
![]() |
![]() | (17) |
and coupling conditions
![]() | (18) |
are satisfied,
Here
-are set real numbers at which the condition of unlimited solvability of the problem (26.15)-(26.18) is fulfilled. Search of the decision leads to the separate system of the integrated equations:
![]() | (19) |
Here
![]() |
![]() |
![]() |
![]() |
- are the eigenfunctions of the corresponding to direct and dual Sturm–Liouville problems [5]. They are connected with this mixed boundary problem (15)-(18). Application of the transformation operator
leads the studied task to the uniform case, that is to the problem (10)-(12).
Truth of the result is established by the scheme proposed in the proof of theorem 2:
Theorem 3. If the function
for each of them the condition is fulfilled
![]() |
that system of the separate integrated equations (26.19) has the unique solution
the solution is found by the formula:
![]() | (20) |
where
![]() |
![]() |
- is the Fourier transform on the piecewise-homogeneous axis with the coupling points [5].
7. Power Function with Discontinuous Coefficients and its Application
We consider the Fourier transform of the Delta function [15]
![]() |
consequently
![]() |
We find as a consequence
![]() |
Let’s define analog of the power function as follows
![]() |
We find equality from the definition of the transformation operators
![]() |
The last equality means that the power function with discontinuous coefficients is obtained by the action of the transformation operator to the power function.
Theorem 4. The ratio connects the generalized power function and differentiation
![]() |
Proof. We have a chain of equalities
![]() |
8. Retrospective Problem for the System of the Diffusion Equations
Let’s return to the solution of the separate system of integral equations (8) in space of the generalized functions [16]
![]() | (21) |
![]() |
We will convert the found solution. An analogue of the series for the function
is of the form
![]() | (22) |
From decomposition theorem have
![]() | (23) |
Let’s find the decomposition of the generalized series for eigenfunction from definition of the generalized power function
![]() |
Let’s substitute this decomposition in a formula (22) and let’s integrate term by term. We come to the formula (23) in which
![]() | (24) |
Let’s substitute decomposition of eigenfunction
in the generalized power series in a formula (21) we will receive
![]() |
where designation is accepted
![]() |
Thus, the solution of the retrospective problem (21) is obtained.
Remark 1. The generating function for the
is the form
![]() |
Corollary 1. If to choose
then we obtain
![]() |
functions
-are the classical Hermite polynomials.
9. Retrospective Problem for Fractal System of the Diffusion Equations
Retrospective problem for fractal system of diffusion equations in the space of generalized functions
leads to the separatist system of integral equations:
![]() | (25) |
![]() |
where
- the Mittag-Leffler function [9].
We get the solution of the fractal retrospective problem repeating reasoning’s from paragraph 7.
![]() |
where designation is accepted
![]() |
You can find the generating functions for
:
![]() | (26) |
We find the explicit expression for the functions
. We find decomposition of the left part of the formula (26) in the series
![]() |
We get the expression, comparing the two views 
![]() |
If
, the formula takes the form
![]() |
Define a "fractal" generalization of the Hermite polynomials
![]() |
then the solution of the fractal retrospective problem has the form
![]() |
Corollary 2. In the hyperbolic case
the solution of the retrospective problem has the form
![]() |
Where
![]() | (27) |
Proof. We replace
for
in formula (16) for the solution of the direct problem. As a result, we get the formula
![]() |
We use the formula (26)
![]() |
In the end we find the solution of the Cauchy problem for the hyperbolic equation
![]() |
Thus,
![]() |
Remark. If, as an example, take
![]() |
the
and, so,
![]() |
10. The Inverse Dirichlet Problem for a Half Plane
Solution of the inverse Dirichlet problem for the right half-plane has the form:
![]() | (28) |
Let’s repeat the above reasoning. Let’s receive expressions for analogues of Hermite polynomials:
![]() |
As a result for the solution of the inverse Dirichlet problem we receive the representation in the form of the sum of the Taylor series:
![]() |
Corollary 3. If the function
admit continued with the real axis of the complex plane as a whole, the
![]() |
Example 2. Let
, then
therefore, we find
![]() |
11. Conclusion
In this article the formal solution of the retrospective problem is provided. The third aspect in determining the well-posed problem is not taken into account. Theorem of existence and uniqueness of solution are given. From the analysis of the formula (9): the solution of the retrospective heat problem with discontinuous coefficients is received by replacement in the final result of the Hermite functions [1] on the Hermite functions with discontinuous coefficients, defined in the article. The derivatives
are necessary to replace on
. The noticed analogy allows hoping on the possibility of obtaining the solution of problems of mathematical physics in which the Hermite functions with discontinuous coefficients.
References
| [1] | F.M. Mors, G. Fishbah, Methods of theoretical physics, 1958. | ||
In article | |||
| [2] | Yaremko, O.E. Matrix integral Fourier transforms for problems with discontinuous coefficients and transformation operators (2007) Doklady Mathematics, 76 (3), pp. 323-325. | ||
In article | CrossRef | ||
| [3] | O.M. Alifanov, Inverse problems of heat exchange, M, 1988, p. 279. | ||
In article | |||
| [4] | O.M. Alifanov, B.A. Artyukhin, S.V. Rumyancev, The extreme methods of solution of ill-posed problems, M, 1988, p. 288. | ||
In article | |||
| [5] | J.V. Beck, V. Blackwell, C.R. Clair, Inverse Heat Conduction. Ill-Posed Problems , M, 1989, p. 312. | ||
In article | |||
| [6] | V.K. Ivanov, V.V. Vasin, V.P. Tanana, Theory of linear ill-posed problems and its applications, M, 1978, p. 206. | ||
In article | |||
| [7] | M.M. Lavrentev, Some ill-posed problems of mathematical physics, Novosibirsk, AN SSSR,1962, p. 92. | ||
In article | |||
| [8] | A.N. Tikhonov, V. Ya. Arsenin, Methods of solution of ill-posed problems, M,1979, p. 288. | ||
In article | |||
| [9] | M.M. Dzhrbashyan, Integral Transforms and Representations of Functions in the Complex Domain, M, 1966. | ||
In article | |||
| [10] | Tikhonov, A. N.; Arsenin, V. Y. (1977). Solutions of Ill-Posed Problems. New York: Winston. | ||
In article | |||
| [11] | Yaremko O.E. Transformation operator and boundary value problems Differential Equation. Vol.40, No. 8, 2004, pp.1149-1160. | ||
In article | CrossRef | ||
| [12] | Bavrin, I.I., Yaremko, O.E. Transformation Operators and Boundary Value Problems in the Theory of Harmonic and Biharmonic Functions (2003) Doklady Mathematics, 68 (3), pp. 371-375. | ||
In article | |||
| [13] | Aldrich, John (2006), "Eigenvalue, eigenfunction, eigenvector, and related terms", in Jeff Miller (Editor), Earliest Known Uses of Some of the Words of Mathematics, retrieved 2006-08-22. | ||
In article | |||
| [14] | Bowen, Ray M.; Wang, Chao-Cheng (1980), Linear and multilinear algebra, Plenum Press, New York, NY. | ||
In article | |||
| [15] | Arfken, G. B.; Weber, H. J. (2000), Mathematical Methods for Physicists (5th ed.), Boston, MA: Academic Press. | ||
In article | PubMed | ||
| [16] | Yu. V. Egorov (1990). "A contribution to the theory of generalized functions". Russ. Math. Surveys (Uspekhi Mat. Nauk)45 (5): 1-49. | ||
In article | |||
OPEN ACCESS
PEER-REVIEWED













































































































In article
CiteULike
Delicious

