Next Article in Journal
Multitask Emergency Logistics Planning under Multimodal Transportation
Previous Article in Journal
A Study on the Dynamic Evolution and Regional Differences of Public Capital and Return to Capital in China
Previous Article in Special Issue
General Relativistic Space-Time with η1-Einstein Metrics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solutions for Multitime Reaction–Diffusion PDE

by
Cristian Ghiu
1,‡ and
Constantin Udriste
2,*,†,‡
1
Department of Mathematical Methods and Models, Faculty of Applied Sciences, University Politehnica of Bucharest, Splaiul Independentei 313, Sector 6, RO-060042 Bucharest, Romania
2
Department of Mathematics and Informatics, Faculty of Applied Sciences, University Politehnica of Bucharest, Splaiul Independentei 313, Sector 6, RO-060042 Bucharest, Romania
*
Author to whom correspondence should be addressed.
Second address: Academy of Romanian Scientists, Ilfov 3, Sector 5, RO-050044 Bucharest, Romania.
These authors contributed equally to this work.
Submission received: 26 August 2022 / Revised: 22 September 2022 / Accepted: 29 September 2022 / Published: 3 October 2022
(This article belongs to the Special Issue New Advances in Differential Geometry and Optimizations on Manifolds)

Abstract

:
A previous paper by our research group introduced the nonlinear multitime reaction–diffusion PDE (with oblique derivative) as a generalized version of the single-time model. This paper states and uses some hypotheses that allow the finding of some important explicit families of the exact solutions for multitime reaction–diffusion PDEs of any dimension that have a multitemporal directional derivative term. Some direct methods for determining the exact solutions of nonlinear PDEs from mathematical physics are presented. In the single-time case, our methods present many advantages in comparison with other known approaches. Particularly, we obtained classes of ODEs and classes of PDEs whose solutions generate solutions of the multitime reaction–diffusion PDE.

1. Introduction

Let μ R , μ 0 and let m , n N , m 2 , n 2 . We consider the functions
  • f : Ω 1 R , f of class C 1 , with Ω 1 R n , Ω 1 open;
  • for α { 1 , 2 , , m } , h α : Ω 2 R , h α of class C 1 , with Ω 2 R m , Ω 2 open;
    h ( t ) = h 1 ( t ) , h 2 ( t ) , , h m ( t ) , for all t = ( t 1 , t 2 , , t m ) Ω 2 .
These functions define the multitime reaction–diffusion PDE
h α ( t ) u t α = μ n u x n + f u , u x , , n 1 u x n 1 .
The present paper deals with new exact solutions admitted by multitime nonlinear reaction–diffusion PDEs (1) of a fairly general form that depends on one or more arbitrary functions. Particularly, solutions are also being built for simplified multitime reaction–diffusion PDEs of type (11). It is important to underline that exact solutions of mathematical physics PDEs, which contain arbitrary functions and, therefore, have a significant generality, is of great practical interest for evaluating the accuracy of various numerical and approximate analytical methods in order to solve corresponding initial boundary value problems.
The additional information we need is distributed as follows: (i) exact solutions of a model nonlinear single-time reaction–diffusion equation [1,2,3,4], single-time soliton solutions and fractional effects [5], single-time wave solutions [6], single-time Schrödinger solitons [7]; (ii) methods for constructing complex solutions of nonlinear PDEs using simpler solutions [8]; (iii) techniques to create multitime extensions of single-time ODEs and PDEs [9]; (iv) multitime Floquet theory [10]; (v) multitime wave functions [11,12]; (vi) multitime Boussinesq solitons [13], multitime Rayleigh solitons [14], multitime reaction–diffusion solitons [15]; (vii) discrete diagonal recurrences [16], discrete multiple recurrences [17], linear discrete multitime multiple recurrences [18];
The ingredients used to find the solutions of multitime reaction–diffusion PDEs (1) are: (i) the solutions of PDE (3) (i.e., first integrals of the ODE system (5)); (ii) the solutions of ODE (6); and a particular solution of PDE (4).
Here is the plan for the rest of the article. In Section 2, we describe how we obtain the solutions via the first integrals of an associated ODE system. In Section 3, we try to argue that, despite the obvious impossibility of developing a useful general theory that encompasses all, there nevertheless exists an ODE that generates solutions for PDE. Section 4 further develops this idea for simplified multitime reaction–diffusion PDEs. In Section 5, we discuss one of the most fundamental problems in PDE, namely, simplified multitime reaction–diffusion PDEs in a Riemannian setting. Finally, Section 6 tries to identify some of the main goals and open problems for exact solutions to PDEs.

2. Solutions via First Integrals

The first integrals of an associated symmetric ODE system are generators of solutions of the multitime reaction–diffusion PDE.
Definition 1.
Let Ω 0 R × Ω 2 be an open subset. A function u : Ω 0 R is called the solution of PDE (1) if
  • u is of class C 1 ;
  • for all ( x , t ) Ω 0 , there exist 2 u x 2 ( x , t ) , , n u x n ( x , t ) ;
  • the functions 2 u x 2 , , n 1 u x n 1 : Ω 0 R are continuous (if n = 2 , this condition does not appear);
  • for all ( x , t ) Ω 0 , we have u ( x , t ) , u x ( x , t ) , , n 1 u x n 1 ( x , t ) Ω 1 ;
  • function u verifies the relation (1) on Ω 0 .
From relation (1), it follows that function n u x n is continuous on the set Ω 0 .
Let k R be a fixed constant. We attached another PDE to multitime reaction–diffusion PDE (1)
h α ( t ) φ t α = μ n φ x n k φ x + f u , φ x , , n 1 φ x n 1 .
Let us consider that φ ( x , t ; c ) , with c = ( c 1 , , c r ) , is a solution of PDE (2), which depends on r parameters. Suppose that φ is a C 1 with respect to all arguments. We accept that 2 φ x 2 , , n 1 φ x n 1 are continuous with respect to all arguments (if n = 2 , this condition does not appear).
Let us consider the C 1 functions P 1 ( · ) , P 2 ( · ) , , P r ( · ) , P r + 1 ( · ) , as r + 1 solutions of the PDE
h α ( t ) w t α ( t ) = 0 .
We note P ( t ) = P 1 ( t ) , P 2 ( t ) , , P r ( t ) .
Suppose that v 0 ( · ) is a C 1 solution of the linear PDE
h α ( t ) v t α ( t ) = 1 .
Let us show that u ( x , t ) : = φ x + k v 0 ( t ) + P r + 1 ( t ) , t ; P ( t ) is a solution of PDE (1) (we assume that functions can be composed); in the case k = 0 , v 0 ( · ) is no longer considered. We denote σ ( x , t ) : = x + k v 0 ( t ) + P r + 1 ( t ) . Then
h α ( t ) u t α ( x , t ) = h α ( t ) φ x σ ( x , t ) , t ; P ( t ) k v 0 t α ( t ) + P r + 1 t α ( t ) + h α ( t ) φ t α σ ( x , t ) , t ; P ( t ) + h α ( t ) j = 1 r φ c j σ ( x , t ) , t ; P ( t ) P j t α ( t ) = k h α ( t ) v 0 t α ( t ) + h α ( t ) P r + 1 t α ( t ) φ x σ ( x , t ) , t ; P ( t ) + h α ( t ) φ t α σ ( x , t ) , t ; P ( t ) + j = 1 r φ c j σ ( x , t ) , t ; P ( t ) h α ( t ) P j t α ( t ) .
Since h α ( t ) v 0 t α ( t ) = 1 and h α ( t ) P j t α ( t ) = 0 (if 1 j r + 1 ), we obtain
h α ( t ) u t α ( x , t ) = h α ( t ) φ t α σ ( x , t ) , t ; P ( t ) + k φ x σ ( x , t ) , t ; P ( t ) .
μ n u x n ( x , t ) + f u ( x , t ) , u x ( x , t ) , , n 1 u x n 1 ( x , t ) = μ n φ x n σ ( x , t ) , t ; P ( t )
+ f u σ ( x , t ) , t ; P ( t ) , φ x σ ( x , t ) , t ; P ( t ) , , n 1 φ x n 1 σ ( x , t ) , t ; P ( t )
= h α ( t ) φ t α σ ( x , t ) , t ; P ( t ) + k φ x σ ( x , t ) , t ; P ( t ) ;
the last equality is true for the functions φ · , · ; c , which are solutions of PDE (2).
It follows that u is a solution of PDE (1).
Analogously, if P 1 ( t ) is a solution of PDE (3) and v 0 ( t ) is a solution of PDE (4), and if φ ( x , t ) is a solution of the PDE (2), then u ( x , t ) : = φ x + k v 0 ( t ) + P 1 ( t ) , t is a solution of PDE (1) (assuming that the composition of functions can be performed).
Therefore, we proved the following four propositions.
Proposition 1.
Let k R be a fixed constant. Let φ : Ω R be a solution of PDE (2), with Ω R × Ω 2 , Ω open; let P 1 : D 1 R be a C 1 solution of PDE (3), with D 1 Ω 2 , D 1 open; let v 0 : D 0 R be a C 1 solution of PDE (4), with D 0 Ω 2 , D 0 open.
Let Ω 0 : = ( x , t ) R × ( D 0 D 1 ) | x + k v 0 ( t ) + P 1 ( t ) , t Ω .
Suppose that set Ω 0 is nonvoid (it is obvious that Ω 0 is an open set).
Then, function u : Ω 0 R ,
u ( x , t ) = φ x + k v 0 ( t ) + P 1 ( t ) , t , f o r   a l l ( x , t ) Ω 0 ,
is a solution of PDE (1).
Proposition 2.
Let k R be a fixed constant. Let Ω R × Ω 2 × R r , Ω open, and let φ : Ω R be a C 1 function with respect to all arguments. Suppose there exist 2 φ x 2 , , n 1 φ x n 1 continuous functions with respect to all arguments (if n = 2 , this condition does not appear).
Denote p r 3 ( Ω ) : = c R r | t h e r e   e x i s t s ( x , t ) R × Ω 2 , s u c h   t h a t ( x , t , c ) Ω . Suppose that, for all c p r 3 ( Ω ) , function φ ( · , · ; c ) : ( x , t ) R × Ω 2 | ( x , t , c ) Ω R is a solution of PDE (2).
Let P 1 , , P r , P r + 1 : D 1 R be C 1 solutions of PDE (3), with D 1 Ω 2 , D 1 open; let v 0 : D 0 R be a C 1 solution of PDE (4), with D 0 Ω 2 , D 0 open.
Let Ω 0 : = ( x , t ) R × ( D 0 D 1 ) | x + k v 0 ( t ) + P r + 1 ( t ) , t ; P 1 ( t ) , , P r ( t ) Ω .
Suppose the Ω 0 is nonvoid (it is obvious that Ω 0 is an open set).
Then, function u : Ω 0 R ,
u ( x , t ) = φ x + k v 0 ( t ) + P r + 1 ( t ) , t ; P 1 ( t ) , , P r ( t ) , f o r   a l l ( x , t ) Ω 0 ,
is a solution of PDE (1).
Proposition 3.
Let φ : Ω R be a solution of the PDE (1), with Ω R × Ω 2 , Ω open; let P 1 : D 1 R be a C 1 solution of PDE (3), with D 1 Ω 2 , D 1 open.
Denote Ω 0 : = ( x , t ) R × D 1 | x + P 1 ( t ) , t Ω .
Suppose the set Ω 0 is nonvoid (it is obvious that Ω 0 is an open set).
Then, function
u : Ω 0 R , u ( x , t ) = φ x + P 1 ( t ) , t , f o r   a l l ( x , t ) Ω 0 ,
is a solution of PDE (1).
Proposition 4.
Let Ω R × Ω 2 × R r , Ω open, and let φ : Ω R be a C 1 function with respect to all arguments. Suppose there exist the functions 2 φ x 2 , , n 1 φ x n 1 , which are continuous with respect to all arguments (if n = 2 , this condition does not appear).
Denote p r 3 ( Ω ) : = c R r | t h e r e e x i s t s ( x , t ) R × Ω 2 , s u c h t h a t ( x , t , c ) Ω . Suppose that, for all c p r 3 ( Ω ) , the function φ ( · , · ; c ) : ( x , t ) R × Ω 2 | ( x , t , c ) Ω R is a solution of PDE (1).
Let P 1 , , P r , P r + 1 : D 1 R be C 1 solutions of PDE (3), with D 1 Ω 2 , D 1 open.
Denote Ω 0 : = ( x , t ) R × D 1 | x + P r + 1 ( t ) , t ; P 1 ( t ) , , P r ( t ) Ω .
Suppose that set Ω 0 is nonvoid (it is obvious that Ω 0 is open set).
Then, function u : Ω 0 R ,
u ( x , t ) = φ x + P r + 1 ( t ) , t ; P 1 ( t ) , , P r ( t ) , f o r   a l l ( x , t ) Ω 0 ,
is a solution of PDE (1).
The C 1 solutions of PDE (3) are the first integrals of the ODE system
d t 1 d τ ( τ ) = h 1 t ( τ ) d t 2 d τ ( τ ) = h 2 t ( τ ) d t m d τ ( τ ) = h m t ( τ ) ,
written in symmetric form d t 1 h 1 ( t ) = d t 2 h 2 ( t ) = = d t m h m ( t ) .
We use the following remark (similar to the one in paper [19]), in which it is shown how PDE (3) solutions are obtained, i.e., first integrals of ODEs system (5).
Remark 1.
Let t 0 Ω 2 , such that h ( t 0 ) ( 0 , 0 , , 0 ) ; hence there exists α 0 { 1 , 2 , , m } , that satisfies h α 0 ( t 0 ) 0 . Let F 1 , F 2 , …, F m 1 : V R , be C 1 functional independent first integrals of the ODE system (5) (in the sense specified in [19]), with V Ω 2 , t 0 V , V open and connected such that, for all t V we have h α 0 ( t ) 0 .
Denote F ( t ) : = F 1 ( t ) , F 2 ( t ) , , F m 1 ( t ) , t V .
Any solution of PDE (3), defined in a neighborhood of t 0 has the form w ( t ) = E F ( t ) , for all t W 0 , where E : U R is a C 1 function, U R m 1 , U open; t 0 W 0 V , W 0 open; F ( t ) U , for all t W 0 (see the results from Section 4.1, of the paper [19]).
Remark 2.
It follows that functions P j ( · ) , which appeared above in this section, have (in a neighborhood of t 0 ) the form specified in Remark 1; i.e., P j ( t ) = E j F ( t ) , for all t W j , where E j : U j R is a C 1 function, U j R m 1 , U j open; t 0 W j V , W j open; F ( t ) U j , for all t W j (for each j { 1 , 2 , , r + 1 } ).
We can consider that functions P 1 , P 2 , , P r + 1 have the same domain of definition D, choosing D = W 1 W 2 W r + 1 . We notice that D is an open set, t 0 D V , and, for all t D , we have F ( t ) U 1 U 2 U r + 1 .

3. Solutions via an Adapted ODE

We associate to the PDE (1) an n-th order ODE
μ y ( n ) k y + f y , y , y , , y ( n 1 ) = 0
(which depends on the parameter k R ).
For k = 0 , Equation (6) becomes
μ y ( n ) + f y , y , y , , y ( n 1 ) = 0 .
It can be seen that if y ( · ) is a solution of Equation (6) (respectively of Equation (7)), then the function φ ( x , t ) : = y ( x ) is a solution of PDE (2) (respectively of PDE (1)). From Propositions 1, 2, 3, and 4, four other statements are immediately obtained.
Proposition 5.
Let k R be a fixed constant. Let ψ : I R be a solution of ODE (6), with I R , I open interval; let P 1 : D 1 R be a C 1 solution of PDE (3), with D 1 Ω 2 , D 1 open; let v 0 : D 0 R be a C 1 solution of PDE (4), with D 0 Ω 2 , D 0 open.
Denote Ω 0 : = ( x , t ) R × ( D 0 D 1 ) | x + k v 0 ( t ) + P 1 ( t ) I .
Suppose that set Ω 0 is nonvoid (obviously, Ω 0 is an open set).
Then, function u : Ω 0 R ,
u ( x , t ) = ψ x + k v 0 ( t ) + P 1 ( t ) , f o r   a l l ( x , t ) Ω 0 ,
is a solution of PDE (1).
Proposition 6.
Let k R be a fixed constant, and let I R , I be an open interval. Suppose G R r , G is open.
Let ψ : I × G R be a C 1 function with respect to all arguments. Suppose there exist 2 ψ x 2 , , n 1 ψ x n 1 continuous functions (with respect to all arguments) (if n = 2 , this condition does not appear). Suppose that, for all c G , function ψ ( · ; c ) : I R is a solution of ODE (6).
Suppose P 1 , , P r , P r + 1 : D 1 R are C 1 solutions of PDE (3), with D 1 Ω 2 , D 1 open; let v 0 : D 0 R be a C 1 solution of PDE (4), with D 0 Ω 2 , D 0 open.
We define Ω 0 : = ( x , t ) R × ( D 0 D 1 ) | x + k v 0 ( t ) + P r + 1 ( t ) ; P 1 ( t ) , , P r ( t ) I × G .
Suppose that set Ω 0 is nonvoid (obviously, set Ω 0 is open).
Then, function u : Ω 0 R ,
u ( x , t ) = ψ x + k v 0 ( t ) + P r + 1 ( t ) ; P 1 ( t ) , , P r ( t ) , f o r   a l l ( x , t ) Ω 0 ,
is a solution of PDE (1).
Proposition 7.
Let ψ : I R be a solution of ODE (7), with I R , I open interval; let P 1 : D 1 R be a C 1 solution of PDE (3), with D 1 Ω 2 , D 1 open.
Denote Ω 0 : = ( x , t ) R × D 1 | x + P 1 ( t ) I .
Suppose that set Ω 0 is nonvoid (obviously, Ω 0 is an open set).
Then, function u : Ω 0 R ,
u ( x , t ) = ψ x + P 1 ( t ) , f o r   a l l ( x , t ) Ω 0 ,
is a solution of PDE (1).
Proposition 8.
Let I R , I be an open interval, and G R r , G open.
Let ψ : I × G R be a C 1 function with respect to all arguments. Suppose there exist functions 2 ψ x 2 , , n 1 ψ x n 1 , which are continuous with respect to all arguments (if n = 2 , this condition does not appear). Suppose that, for all c G , function ψ ( · ; c ) : I R is a solution of ODE (7).
We introduce P 1 , , P r + 1 : D 1 R as C 1 solutions of PDE (3), with D 1 Ω 2 , D 1 open.
Denote Ω 0 : = ( x , t ) R × D 1 | x + P r + 1 ( t ) ; P 1 ( t ) , , P r ( t ) I × G .
Suppose that set Ω 0 is nonvoid (obviously, Ω 0 is an open set).
Then, function u : Ω 0 R ,
u ( x , t ) = ψ x + P r + 1 ( t ) ; P 1 ( t ) , , P r ( t ) , f o r   a l l ( x , t ) Ω 0 ,
is a solution of PDE (1).
For λ R , c = ( c 1 , , c n ) Ω 1 , k R , let η ( · ; λ ; c ; k ) be maximal solution y ( · ) of Equation (6), which verifies y ( λ ) = c 1 , y ( λ ) = c 2 , y ( λ ) = c 3 , ,   y ( n 1 ) ( λ ) = c n . The domain of definition of this solution is an open interval I ( λ ; c ; k ) , with λ I ( λ ; c ; k ) .
Let A : = ( x ; λ ; c ; k ) R × R × Ω 1 × R | x I ( λ ; c ; k ) . Set A is open. The functions η , η x , 2 η x 2 , , n 1 η x n 1 : A R are of class C 1 . For any ( x ; λ ; c ; k ) A , we have η ( x ; λ ; c ; k ) , η x ( x ; λ ; c ; k ) , 2 η x 2 ( x ; λ ; c ; k ) , , n 1 η x n 1 ( x ; λ ; c ; k ) Ω 1 .
For any λ R , c = ( c 1 , , c n ) Ω 1 , k R , we have ( λ ; λ ; c ; k ) A , and
η ( λ ; λ ; c ; k ) = c 1 , η x ( λ ; λ ; c ; k ) = c 2 , 2 η x 2 ( λ ; λ ; c ; k ) = c 3 , , n 1 η x n 1 ( λ ; λ ; c ; k ) = c n .
For any λ R , c Ω 1 , k R , we have I ( λ ; c ; k ) = λ + I ( 0 ; c ; k ) , and
η ( x ; λ ; c ; k ) = η ( x λ ; 0 ; c ; k ) , for all x I ( λ ; c ; k ) ,
or, more generally, I ( λ + ξ ; c ; k ) = ξ + I ( λ ; c ; k ) , and
η ( x + ξ ; λ + ξ ; c ; k ) = η ( x ; λ ; c ; k ) , for all x I ( λ ; c ; k ) , for all ξ R .
(For the interval ( α , β ) , the set ξ + ( α , β ) is the interval ( ξ + α , ξ + β ) . )
Let x 0 R , c 0 Ω 1 , k 0 R . Since ( x 0 ; x 0 ; c 0 ; k 0 ) A , and A is an open set, it follows that there exists ρ > 0 , and there exist the open and bounded intervals Λ , L, and there exists the open set G Ω 1 , such that: x 0 Λ , c 0 G , k 0 L , and ( x 0 ρ , x 0 + ρ ) × Λ × G × L A .
Proposition 9.
Let x 0 R , t 0 Ω 2 , with h ( t 0 ) ( 0 , 0 , , 0 ) .
Denote ρ R , ρ > 0 . Let L R , L be an open and bounded interval, and let G Ω 1 , G open, such that: ( x 0 ρ , x 0 + ρ ) × { x 0 } × G × L A .
Let F 1 , F 2 , …, F m 1 : V R , be C 1 , functional independent first integrals for ODEs system (5), with V Ω 2 , t 0 V , V open. Let F ( t ) : = F 1 ( t ) , F 2 ( t ) , , F m 1 ( t ) , t V .
We consider the C 1 functions E 1 , E 2 , , E n + 1 : U R , E n + 2 : U L , with U R m 1 , U open. Suppose that, for all s U , we have E 1 ( s ) , E 2 ( s ) , , E n ( s ) G .
Suppose that there exists D V , D open, with t 0 D , such that F ( D ) U .
Let v 0 : D 0 R be a C 1 solution of PDE (4), with D 0 Ω 2 , D 0 open, t 0 D 0 .
For all j { 1 , 2 , , n 1 , n , n + 2 } , and for all t D , we denote P j ( t ) : = E j F ( t ) .
For all t D , we denote P n + 1 ( t ) : = E n + 1 F ( t ) E n + 1 F ( t 0 ) .
For all t D 0 , we take v ˜ 0 ( t ) : = v 0 ( t ) v 0 ( t 0 ) .
Then, there exists D 1 D D 0 , D 1 open, with t 0 D 1 , such that:
( i ) for all x x 0 ρ 2 , x 0 + ρ 2 , for all t D 1 , we have
x + P n + 2 ( t ) · v ˜ 0 ( t ) + P n + 1 ( t ) ( x 0 ρ , x 0 + ρ ) ;
( i i ) the function u : x 0 ρ 2 , x 0 + ρ 2 × D 1 R ,
u ( x , t ) = η x + P n + 2 ( t ) · v ˜ 0 ( t ) + P n + 1 ( t ) ; x 0 ; P 1 ( t ) , , P n ( t ) ; P n + 2 ( t ) , f o r   a l l ( x , t ) x 0 ρ 2 , x 0 + ρ 2 × D 1 ,
is a solution of PDE (1).
Proof. 
( i ) Let I : = ( x 0 ρ , x 0 + ρ ) , I 0 : = x 0 ρ 2 , x 0 + ρ 2 .
Because the interval L is bounded, there exists M > 0 , such that, for any k L , we have | k | M .
Since the functions v ˜ 0 and P n + 1 are continuous in t 0 , and v ˜ 0 ( t 0 ) = P n + 1 ( t 0 ) = 0 , it follows that there exists R > 0 , such that, for any t D D 0 B ( t 0 , R ) , the following relations hold
| v ˜ 0 ( t ) | < ρ 2 ( M + 1 ) a n d | P n + 1 ( t ) | < ρ 2 ( M + 1 )
(set B ( t 0 , R ) is an open ball of center t 0 and radius R, included in R m ).
Denote D 1 : = D D 0 B ( t 0 , R ) ; D 1 is open, t 0 D 1 , and D 1 D D 0 .
Let us show that
x + k v ˜ 0 ( t ) + P n + 1 ( t ) I , f o r   a l l x I 0 , f o r   a l l t D 1 , f o r   a l l k L .
We note x I 0 , t D 1 , k L . We use the obvious inequality
| x + k v ˜ 0 ( t ) + P n + 1 ( t ) x 0 | | x x 0 | + | k | · | v ˜ 0 ( t ) | + | P n + 1 ( t ) | .
From x I 0 , we deduce | x x 0 | < ρ 2 . Since k L , we have | k | M . Further, using inequalities (9), we obtain
| x x 0 | + | k | · | v ˜ 0 ( t ) | + | P n + 1 ( t ) | < ρ 2 + M · ρ 2 ( M + 1 ) + ρ 2 ( M + 1 ) = ρ .
It follows that | x + k v ˜ 0 ( t ) + P n + 1 ( t ) x 0 | < ρ , i.e., x + k v ˜ 0 ( t ) + P n + 1 ( t ) I .
We proved statement (10); from this statement ( i ) follows, since for all t D 1 , we have P n + 2 ( t ) = E n + 2 F ( t ) L , and if in (10) we choose k = P n + 2 ( t ) , we find (8).
( i i ) For each k L , we consider the function ψ k : I × G R ,
ψ k ( x ; c ) = η ( x ; x 0 ; c ; k ) , for all ( x ; c ) I × G .
For all c G , function ψ k ( · ; c ) : I R is a solution of Equation (6).
The functions P 1 ( · ) , P 2 ( · ) , , P n ( · ) , P n + 1 ( · ) , defined on D 1 , are C 1 solutions of PDE (3). Function v ˜ 0 : D 0 R is a C 1 solution of PDE (4).
Denote Ω 0 ( k ) = ( x , t ) R × D 1 | x + k v ˜ 0 ( t ) + P n + 1 ( t ) ; P 1 ( t ) , , P n ( t ) I × G .
Since for all t D 1 , we have P 1 ( t ) , P 2 ( t ) , , P n ( t ) G , we deduce that
Ω 0 ( k ) = ( x , t ) R × D 1 | x + k v ˜ 0 ( t ) + P n + 1 ( t ) I .
According to Proposition 6, function θ k ( · , · ) : Ω 0 ( k ) R ,
θ k ( x , t ) = ψ k x + k v ˜ 0 ( t ) + P n + 1 ( t ) ; P 1 ( t ) , , P n ( t ) , f o r   a l l ( x , t ) Ω 0 ( k ) ,
is a solution of PDE (1) (for all k L ).
From statement (10) we find that, for all k L , we have I 0 × D 1 Ω 0 ( k ) ; it follows that we can define the function
φ : I 0 × D 1 × L R , φ ( x , t ; k ) = θ k ( x , t ) , for all ( x , t ; k ) I 0 × D 1 × L , i.e.
φ ( x , t ; k ) = ψ k x + k v ˜ 0 ( t ) + P n + 1 ( t ) ; P 1 ( t ) , , P n ( t ) , or
φ ( x , t ; k ) = η x + k v ˜ 0 ( t ) + P n + 1 ( t ) ; x 0 ; P 1 ( t ) , , P n ( t ) ; k , for all ( x , t ; k ) I 0 × D 1 × L .
For any k L , function θ k ( · , · ) is a solution of PDE (1). We deduce that for any k L , function φ ( · , · ; k ) : I 0 × D 1 R is a solution of PDE (1).
Function P n + 2 ( · ) , defined on D 1 , is a C 1 solution of PDE (3).
We apply Proposition 4 for function φ ; we have r = 1 and Ω = I 0 × D 1 × L . Instead of P 1 ( · ) from Proposition 4, we set P n + 2 ( · ) , and P r + 1 ( · ) will be the null function.
Let Ω 0 = ( x , t ) R × D 1 | x , t ; P n + 2 ( t ) I 0 × D 1 × L .
Since, for all t D 1 , we have P n + 2 ( t ) L , and hence, Ω 0 = I 0 × D 1 .
According to Proposition 4, function u : I 0 × D 1 R ,
u ( x , t ) = φ x , t ; P n + 2 ( t ) , f o r   a l l ( x , t ) I 0 × D 1 ,
is a solution of PDE (1). We obtain statement ( i i ) . □

4. Simplified Multitime Reaction–Diffusion PDE

The simplified multitime reaction–diffusion PDE, in the direction h = ( h 1 , h 2 , , h m ) , is a PDE of the form (see [15])
h α ( t ) u t α = μ 2 u x 2 + g ( u ) ,
with g : J R a C 1 function, J R , J open interval.
PDE (11) is an equation of type (1), with n = 2 ; we have Ω 1 = J × R , and f ( q 1 , q 2 ) = g ( q 1 ) , for all ( q 1 , q 2 ) J × R .
Let us determine some solutions for a particular case of PDE (11), namely
t 1 u t 1 + + t m 1 u t m 1 + t m u t m = μ 2 u x 2 a u 3 + b u 2 , t m 0 ,
where μ , a and b are real constants, μ > 0 , a > 0 , b 0 (compare with soliton solutions, [15]).
In this case, the ODEs system (5) becomes
d t 1 d τ ( τ ) = t 1 ( τ ) d t 2 d τ ( τ ) = t 2 ( τ ) d t m 1 d τ ( τ ) = t m 1 ( τ ) d t m d τ ( τ ) = t m ( τ ) .
Functions F 1 , F 2 , , F m 1 : R m 1 × ( R { 0 } ) R ,
F α ( t ) = t α t m , f o r a l l t R m 1 × ( R { 0 } ) , f o r   a l l α { 1 , 2 , , m 1 } ,
are m 1 functional independent first integrals of the ODEs system (13).
We note F ( t ) = F 1 ( t ) , F 2 ( t ) , , F m 1 ( t ) = t 1 t m , t 2 t m , , t m 1 t m , t R m 1 × ( R { 0 } ) .
Let E : U R be any C 1 function ( U R m 1 , U open).
The PDE (4) is written
t 1 v t 1 + + t m 1 v t m 1 + t m v t m = 1 .
The function v 0 : R m 1 × ( R { 0 } ) R ,
v 0 ( t ) = ln | t m | , f o r   a l l t R m 1 × ( R { 0 } ) ,
is a particular solution for PDE (14).
Equations (6) and (7), associated to PDE (12), are
μ y k y a y 3 + b y 2 = 0
μ y a y 3 + b y 2 = 0 .
For statements ( i ) , ( i i ) and ( i i i ) , the below are true.
( i ) For k = b 2 μ a , the function y 1 ( x ) = 1 x 2 μ a is a solution for ODE (15).
For k = b 2 μ a , the function y 2 ( x ) = 1 x 2 μ a is a solution for ODE (15).
( i i ) For k = b μ 2 a , the function y 3 ( x ) = b a · 1 1 + c e b x 2 μ a is a solution for ODE (15).
For k = b μ 2 a , the function y 4 ( x ) = b a · 1 1 + c e b x 2 μ a is a solution for ODE (15).
Solutions y 3 and y 4 depend on the real constant c.
( i i i ) Functions y 5 ( x ) = b a · 2 + sinh b x μ a 2 2 + sinh b x μ a , y 6 ( x ) = b a · 2 + sinh b x μ a 2 2 + sinh b x μ a ,
y 7 ( x ) = 12 b μ 9 a μ 2 b 2 x 2 are three solutions for ODE (16).
Using statement ( i ) and Proposition 5, we deduce that the functions
u 1 ( x , t ) = 1 x b 2 μ a ln | t m | + E F ( t ) 2 μ a = 2 μ a x b 2 μ ln | t m | + a E F ( t ) ,
u 2 ( x , t ) = 1 x + b 2 μ a ln | t m | + E F ( t ) 2 μ a = 2 μ a x + b 2 μ ln | t m | + a E F ( t )
are solutions of PDE (12).
Using statement ( i i ) and Proposition 6 (with r = 1 and choosing P r + 1 ( · ) as the null function), we deduce that the functions
u 3 ( x , t ) = b a · 1 1 + E F ( t ) e b 2 μ a x b μ 2 a ln | t m | = b a · 1 1 + E F ( t ) e b x 2 μ a · | t m | b 2 2 a ,
u 4 ( x , t ) = b a · 1 1 + E F ( t ) e b 2 μ a x + b μ 2 a ln | t m | = b a · 1 1 + E F ( t ) e b x 2 μ a · | t m | b 2 2 a
are solutions for PDE (12).
Using statement ( i i i ) and Proposition 7, we deduce that functions
u 5 ( x , t ) = b a · 2 + sinh b x + b E F ( t ) μ a 2 2 + sinh b x + b E F ( t ) μ a , u 6 ( x , t ) = b a · 2 + sinh b x + b E F ( t ) μ a 2 2 + sinh b x + b E F ( t ) μ a ,
u 7 ( x , t ) = 12 b μ 9 a μ 2 b 2 x + E F ( t ) 2 = 12 b μ 9 a μ 2 b x + b E F ( t ) 2
are solutions for PDE (12).
Since function E is arbitrary, we can replace E, which appears in the expression of u 1 and u 2 , with E a ; we can replace E, which appears in the expression of u 5 and u 6 , with μ a b E ; we can replace E, which appears in the expression of u 7 , with E b .
Solutions u 1 , u 2 , u 3 , u 4 , u 5 , u 6 , u 7 of PDE (12) become
u 1 ( x , t ) = 2 μ a x b 2 μ ln | t m | + E F ( t )
u 2 ( x , t ) = 2 μ a x + b 2 μ ln | t m | + E F ( t )
u 3 ( x , t ) = b a · | t m | b 2 2 a | t m | b 2 2 a + E F ( t ) e b x 2 μ a
u 4 ( x , t ) = b a · | t m | b 2 2 a | t m | b 2 2 a + E F ( t ) e b x 2 μ a
u 5 ( x , t ) = b a · 2 + sinh b x μ a + E F ( t ) 2 2 + sinh b x μ a + E F ( t )
u 6 ( x , t ) = b a · 2 + sinh b x μ a + E F ( t ) 2 2 + sinh b x μ a + E F ( t )
u 7 ( x , t ) = 12 b μ 9 a μ 2 b x + E F ( t ) 2 .
For each u j , the domain of definition is chosen as an open set in R × R m 1 × ( R { 0 } ) , such that F ( t ) U and the respective denominator do not cancel. For example, for u 1 , the domain of definition is
( x , t ) R × R m 1 × ( R { 0 } ) | F ( t ) U , a x b 2 μ ln | t m | + E F ( t ) 0 .

5. The Simplified Multitime Reaction–Diffusion PDE in Riemannian Setting

From the physical point of view, it would be more important to further study a simplified multitime reaction–diffusion PDE in a Riemannian setting, which is still an open problem. For that, let ( M , g ) be a smooth, compact Riemannian manifold of dimension n (particularly, n = 2 , 3 ) without a boundary. In the Riemannian setting, the simplified multitime reaction–diffusion PDE is
i h α ( x , t ) u t α + μ Δ g u + f ( u ) = 0 , ( x , t ) M × R m .
Let d x g denote the volume element of the compact Riemannian manifold ( M , g ) and
| u | 2 = u u ¯ , | g u | g 2 = g u , g u ¯ g .
This simplified multitime reaction–diffusion PDE is mainly concerned with the interface between Riemannian geometry and heat theory, but it leads, in a natural way, to questions of functional analysis related to the theory of operators on Hilbert spaces. In some respects, these problems are similar to those studied in the standard Euclidean case, but depending on the Riemannian metric g these might go beyond and provide new aspects to the problem.
Open problems (i) What happens if we consider Riemannian metrics that vary periodically in multitime? (ii) Is a diffusion process determined by its intrinsic Riemannian metric?

6. Conclusions

The problem of exact solutions for PDEs was resuscitated by the publication of the book [20]. Although this book refers to geometrical PDEs (Exact Solutions of Einstein’s Field Equations), there are many other PDEs with a physical sense (Broer–Kaup PDE, Burgers PDE, Fokker–Planck PDE, Hamilton–Jacobi equation PDE, Hamilton–Jacobi–Bellman equation PDE, Diffusion PDE, Klein–Gordon PDE, Korteweg–de Vries PDE, Maxwell PDEs, Navier–Stokes PDEs, Schrödinger PDE, Wave equation, etc.) for which we need exact solutions. Particularly, solitons are exact solutions. In order to show the advantages of the work on exact solutions, one can point out that certain solutions have played very important roles in the solving of physical problems and hence finding exact solutions that have a practical interpretation is of fundamental importance.
Our techniques are similar to the one in paper [19]. There are practical as well as theoretical reasons for studying exact solutions of the multitime reaction–diffusion PDEs (see also single-time case [3,5,6,7,8,21,22,23]). In some cases, the reader can be surprised by the enormous number of exact solutions. In other cases, there are problems with a unique exact solution, for instance, the Kerr and Schwarzschild solutions for the final collapsed state of massive bodies. The single-time reaction–diffusion PDE or the multitime reaction–diffusion PDE can be used as mathematical models in a wide series of domains. One example is represented by many processes during embryonic development, which involves the transport and reaction of molecules, or transport and proliferation of cells, within growing tissues.
Due to the possibilities offered by symbolic computation software such as Maple or Mathematica [24], the computation of the exact solutions of PDEs has become easier and more stimulating for mathematicians and other scientists.
In Maple, the basic Calling Sequences are:
pdsolve(PDE, f, HINT = hint, INTEGRATE, build);
pdsolve(PDEsystem, funcs, HINT, otheroptions).

Author Contributions

Conceptualization, C.G. and C.U.; methodology, all authors; writing—original draft preparation, C.G. and C.U.; writing—review and editing, C.G. and C.U.; validation, C.U. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are grateful to the reviewers for interesting proposals to improve our point of view. Best wishes go to Oana-Maria Pastae, “Constantin Brancusi” University of Tg-Jiu, for the very conscientious revision of the language in this manuscript. Warm thanks go to Valeriu Prepelita (University Politehnica of Bucharest) for proofreading the manuscript regarding duplication problems.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cherniha, M. New exact solutions of nonlinear reaction-diffusion equations. Rep. Math. Phys. 1998, 41, 339–349. [Google Scholar] [CrossRef]
  2. Samsonov, A.M.; Gursky, V.V. Exact solutions to a nonlinear reaction-diffusion equation and hyperelliptic integrals inversion. J. Phys. A Math. Gen. 1999, 32, 6573. [Google Scholar] [CrossRef]
  3. Okoya, S.S. Some exact solutions of a model nonlinear reaction-diffusion equation. Int. Commun. Heat Mass Transfer 1996, 23, 1043–1052. [Google Scholar] [CrossRef]
  4. Cherniha, M. New ansätze and exact solutions for nonlinear reaction-diffusion equations arising in Mathematical Biology. Symmetry Nonlinear Math. Phys. 1997, 1, 138–146. [Google Scholar]
  5. Bashar, M.H.; Inc, M.; Islam, S.R.; Mahmoud, K.H.; Akbar, M.A. Soliton solutions and fractional effects to the time-fractional modified equal width equation. Alexandria Eng. J. 2022, 61, 12539–12547. [Google Scholar] [CrossRef]
  6. Islam, S.R.; Khan, S.; Arafat, S.Y.; Akbar, M.A. Diverse analytical wave solutions of plasma physics and water wave equations. Results Phys. 2022, 40, 105834. [Google Scholar] [CrossRef]
  7. Arafat, S.Y.; Khan, K.; Islam, S.R.; Rahman, M.M. Parametric effects on paraxial nonlinear Schrödinger equation in Kerr media. Chin. J. Phys. 2022, in press. [CrossRef]
  8. Aksenov, A.V.; Polyanin, A.D. Methods for constructing complex solutions of nonlinear PDEs using simpler solutions. Mathematics 2021, 9, 345. [Google Scholar] [CrossRef]
  9. Matei, L.; Pirvan, M.; Petrescu, L.L. Techniques to create multitime extensions of single-time ODEs and PDEs. Int. J. Appl. Math. Inform. 2013, 7, 135–142. [Google Scholar]
  10. Ghiu, C.; Tuliga, R.; Udriste, C.; Tevy, I. Floquet theory for multitime linear diagonal recurrence. UPB Sci. Bull. Ser. A 2016, 78, 3–8. [Google Scholar]
  11. Lienert, M.; Petrat, S.; Tumulka, R. Multi-Time Wave Functions, an Introduction; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  12. Petrat, S.; Tumulka, R. Multi-time equations, classical and quantum. Proc. R. Soc. A 2014, 470, 20130632. [Google Scholar] [CrossRef] [PubMed]
  13. Matei, L.; Udriste, C.; Ghiu, C. Multitime Boussinesq solitons. Int. J. Geom. Methods Mod. Phys. 2012, 9, 1250031. [Google Scholar] [CrossRef]
  14. Matei, L.; Udriste, C. Construction of multitime Rayleigh solitons. UPB Sci. Bull. Ser. A 2014, 76, 29–43. [Google Scholar]
  15. Udriste, C.; Petrescu, L.; Matei, L. Multitime reaction-diffusion solitons. Balkan J. Geom. Appl. 2012, 17, 115–128. [Google Scholar]
  16. Ghiu, C.; Tuliga, R.; Udriste, C.; Tevy, I. Discrete diagonal recurrences and discrete minimal submanifolds. Balkan J. Geom. Appl. 2015, 20, 49–64. [Google Scholar]
  17. Ghiu, C.; Tuliga, R.; Udriste, C. Discrete multiple recurrence. UPB Sci. Bull. Ser. A 2019, 81, 25–38. [Google Scholar]
  18. Ghiu, C.; Tuliga, R.; Udriste, C. Linear discrete multitime multiple recurrence. UPB Sci. Bull. Ser. A 2020, 82, 29–40. [Google Scholar]
  19. Ghiu, C.; Udriste, C.; Petrescu, L.L. Families of solutions of multitemporal nonlinear Schrödinger PDE. Mathematics 2021, 9, 1995. [Google Scholar] [CrossRef]
  20. Stephani, H.; Kramer, D.; Maccallum, M.; Hoenselaers, C.; Herlt, E. Exact Solutions of Einstein’s Field Equations, Cambridge Monographs on Mathematical Physics; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  21. Polyanin, A.D. Construction of exact solutions in implicit form for PDEs: New functional separable solutions of non-linear reaction–diffusion equations with variable coefficients. Int. J. Non-Linear Mech. 2019, 111, 95–105. [Google Scholar] [CrossRef]
  22. Simpson, M.J. Exact solutions of linear reaction-diffusion processes on a uniformly growing domain: Criteria for successful colonization. PLoS ONE 2015, 10, e0117949. [Google Scholar] [CrossRef] [Green Version]
  23. Feng, Z. On explicit exact solutions to the compound Burgers-KdV equation. Phys. Lett. A 2002, 293, 57–66. [Google Scholar] [CrossRef]
  24. Cheb-Terrab, E.S.; von Bulow, K. A computational approach for the analytical solving of Partial Differential Equations. Comput. Phys. Commun. 1995, 90, 102–116. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ghiu, C.; Udriste, C. Solutions for Multitime Reaction–Diffusion PDE. Mathematics 2022, 10, 3623. https://0-doi-org.brum.beds.ac.uk/10.3390/math10193623

AMA Style

Ghiu C, Udriste C. Solutions for Multitime Reaction–Diffusion PDE. Mathematics. 2022; 10(19):3623. https://0-doi-org.brum.beds.ac.uk/10.3390/math10193623

Chicago/Turabian Style

Ghiu, Cristian, and Constantin Udriste. 2022. "Solutions for Multitime Reaction–Diffusion PDE" Mathematics 10, no. 19: 3623. https://0-doi-org.brum.beds.ac.uk/10.3390/math10193623

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop