Next Article in Journal
A Roadside Unit Deployment Optimization Algorithm for Vehicles Serving as Obstacles
Next Article in Special Issue
Application of a Variant of Mountain Pass Theorem in Modeling Real Phenomena
Previous Article in Journal
Fault Diagnosis Analysis of Angle Grinder Based on ACD-DE and SVM Hybrid Algorithm
Previous Article in Special Issue
Non-Relativistic Treatment of the 2D Electron System Interacting via Varshni–Shukla Potential Using the Asymptotic Iteration Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advanced Analytic Self-Similar Solutions of Regular and Irregular Diffusion Equations

by
Imre Ferenc Barna
1,* and
László Mátyás
2
1
Wigner Research Centre for Physics, Konkoly–Thege Miklós út 29-33, H-1121 Budapest, Hungary
2
Department of Bioengineering, Faculty of Economics, Socio-Human Sciences and Engineering, Sapientia Hungarian University of Transylvania, Libertătii sq. 1, 530104 Miercurea Ciuc, Romania
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(18), 3281; https://doi.org/10.3390/math10183281
Submission received: 28 July 2022 / Revised: 5 September 2022 / Accepted: 6 September 2022 / Published: 9 September 2022

Abstract

:
We study the diffusion equation with an appropriate change of variables. This equation is, in general, a partial differential equation (PDE). With the self-similar and related Ansatz, we transform the PDE of diffusion to an ordinary differential equation. The solutions of the PDE belong to a family of functions which are presented for the case of infinite horizon. In the presentation, we accentuate the physically reasonable solutions. We also study time-dependent diffusion phenomena, where the spreading may vary in time. To describe the process, we consider time-dependent diffusion coefficients. The obtained analytic solutions all can be expressed with Kummer’s functions.

1. Introduction

The study of classical transport processes are crucial both from scientific and from engineering points. One of such process is diffusion, which is a quite general phenomena. It can be formulated for particles, which are the mass diffusion, or to energy transport, which is related to thermal conduction. This manuscript is an extension of certain studies, in which we present additional analytic solutions for the regular diffusion equation from symmetry considerations. In our last paper [1], after a historical overview, we presented a new class of analytic solutions derived with the help of the reduction mechanism provides a more detailed analysis of various self-similar and generalized self-similar solutions. Here, we show how the original self-similar trial function can be generalized in a kind of power law expansion. All the presented solutions are new and cannot be found elsewhere in the literature.
The whole diffusion phenomena is, in general, well introduced—with cases and studies—by the monograph of Ghez [2]. Embedding into the larger field of transport processes diffusion is discussed by John Newman and Vincent Battaglia in their series of lectures [3]. Gillespie and Seitariodu provided an introduction to the standard theoretical models for simple Brownian type diffusion [4] in 2013. The diffusion equation can be also related to the Wiener process [5]. The anomalous diffusion was analyzed with statistical methods by Weihua and co-workers [6]. Uchaikin investigated self-similar anomalous diffusion [7]. Ari Arapostathis et al. studied the ergodic control of diffusion processes in a monograph [8]. One may find applications in solid state physics, binary alloys, thin films, etc., in [9,10,11,12,13,14,15]. Defects and diffusion in nanotubes was summarized by Fisher [16]. Diffusion processes a peculiar material in ceramics was analyzed in the monograph by Pelleg [17]. Atomic diffusion processes is a scientific field which was presented in a monograph as well [18]. Diffusion processes are the starting points for reaction diffusion processes [19] or diffusion in porous media [20]. A vitally important application of such mathematical equations is the mathematical modeling of aircraft cabin fires [21,22]. The two dimensional diffusion equation, completed with certain reaction terms may lead to pattern formation, for instance the Turing patterns derived from Schnakenberg equations [23] or Brusselator model [24]. If the diffusion equation is completed with the simplest nonlinear term, the gradient of the variable on the second power we arrive to the Kardar–Parizi–Zwang (KPZ) equation which is the simplest successful model for surface growth phenomena. In two of our former studies we investigated the KPZ equation (with additional noise terms) with the self-similar [25] and the traveling waves Ansatz [26]. Interesting results have been obtained in the study of irregular diffusion [27,28]. Now, in the following we deepen this analysis and present additional analytic solutions with detailed parameter study. We present a series of solutions which are defined on the whole real axis having a decaying and spreading property with additional oscillations. Our main investigation tool, the self-similar Ansatz, helps us to build a link from the regular to irregular diffusion processes, which makes up the second part of the study.
Although the diffusion process can be studied in different dimensions, here, we consider only one Cartesian coordinate therefore the equation reads
C ( x , t ) t = D 2 C ( x , t ) x 2 ,
where C ( x , t ) is the distributions of the particle concentration in space and time and D is the diffusion coefficient. C ( x , t ) in the equation above is considered up to a constant, which means, that C ( x , t ) + C 0 may be also a solution, depending on the initial conditions. The function C ( x , t ) fulfills the necessary smoothness conditions with existing continuous first and second derivatives in respect to time and and space. From causality, the diffusion coefficient should be a positive real number ( D > 0 ) . Numerous physics textbooks gives us the derivation how the fundamental (the Gaussian) solutions can be obtained, e.g., [29,30].
The derivation and analysis of various analytic solutions of physical processes described by various mathematical equations, (e.g., algebraic, differential or partial differential) have crucial importance. As it was shown in our former paper [1] and as it will be shown here, there are far more solutions known for diffusion than the Gaussian and the error functions. We presented some relatively simple solutions (e.g., D t + x 2 / 2 ), other solutions which can derived with the general symmetry analysis method by Clarkson and Kruskal [31], the traveling profile method of Benhamidouche [32] or the self-similar Ansatz of Sedov [33]. Beyond the Gaussian and error functionsm most of our results are expressible with the Kummer special functions. In the following, we try additional two trial functions and enlarge the number of solutions known from the self-similar Ansatz. In the last part of the study, we investigate less regular diffusion processes where the diffusion coefficients gain temporal dependencies. The diffusion equation stands at the basis of more complex equations: in case on the r.h.s. beyond the second derivative, there is a function F ( C ) with certain properties, we can talk about the Kolmogorov–Petrovskii–Piskunov equation [34]. Explicitly, on the r.h.s., the term C ( 1 C ) yields the Fisher equation [35,36]. Environmental aspects can be found in [37]. The term p C + r C q in general means the Newel–Whitehead–Segel (NWS) equation [38,39] and the term C ( 1 C ) ( C α ) where 0 < α < 1 defines the general Zeldovich equation (or Huxley equation) which arises in combustion theory [40]. For certain non-trivial values of p , r and q, one may find exact solutions of the NWS equation [41,42]. Frank–Kamenetzkii used the exp ( a · e b C ) term [43] to explain thermal explosion.
Burgers used the C · C x (where the subscript stands for partial derivation) term to study turbulence [44]. Nariboli and Lin introduced the quadratic Burgers equation with the term of C 2 · C x [45]. Sachdev [46] modified the original Burgers equation and used a third order term of C 3 · C x . Later numerous generalization saw the light of sun by various authors, and the originality of the models are hard to identify. The generalized Huxley equation [47] has the source term of β C ( 1 C δ ) ( C δ γ ) where β ,   δ   , γ are free real parameters. Lastly, we mention the the Burgers–Huxley and the Burgers–Fisher equations [48]. The first has the source term of γ C x + β C ( 1 C δ ) ( C δ γ ) and the second of γ C x + β C ( 1 C ) . Our presented list is of course far from being complete.
Applications in different fields—for instance, plasma physics or condensed matter—one may find in [49,50,51].
We hope that this work may bring deeper understanding in the study of vapor diffusion [52,53], of the one dimensional convection- diffusion-reaction problem [54,55,56], and of diffusive aspects in different flows [57].
Diffusion processes can be coupled to fluid dynamics phenomenon to describe the double (or multiple) diffusive convection phenomena [58] where a convection is driven by two (or more) different density gradients described with different rates of diffusions.
Another way of generalization of diffusion is the application of fractional derivatives. First, consider when the time derivative is fractional. Such a study was conducted in the work of Wyss [59]. The solutions were exactly given and could be expressed with the Fox functions. For the mathematical details of Fox functions, see [60]. Regarding the space fractional diffusion equation, one may find studies in [61]. Comparison of our results to such functions could be the subject of future study.
The analysis and control of coupled neural networks can be conducted with the reaction diffusion term, as was given in the monograph of Wang et al. [62]. It is obvious, but we mention that the diffusion equation has the same form as the heat conduction equation. This field has a wealth of literature as well, from which we mention only two monographs [63,64]. Lastly, not to forget the field of continuously developing numerical methods of PDEs, it is worthwhile to mention the new results obtained by [65,66].

2. Theory and Results

2.1. Self-Similar Ansatz

We start the analysis with the self-similar Ansatz
C ( x , t ) = t α f x t β = t α f ( η ) ,
where α and β are the self-similar exponents being real numbers describing the decay and the spreading of the solution in time and space. These properties makes this Ansatz physically extraordinarily relevant; it was first introduced by Sedov [33]. Together with dimension-analysis method He applied it to various fluid mechanical problems such as incompressible heavy fluids or gas dynamics. Later Zeldovich [67] exhaustively investigated high temperature gas dynamical problems. The most recent and modern literature where the reader can learn the method and the possible connections to dimension-analysis was written by Barenblatt [68] in 2003. He addressed problems such as very intense concentrated flooding or flow in porous media.
For the present diffusion equation, after some trivial algebra (which can be found in our previous paper [1]) we get:
α = arbitrary real number , β = 1 / 2 ,
and there is a clear-cut time-independent ordinary differential equation (ODE) of
α f 1 2 η f = D f .
with the choice of α = 1 / 2 and setting the first integration constant to zero c 1 = 0 we obtain back the well-known Gaussian solution.
This is the so-called fundamental solution and sometimes referred to as source type solution—by mathematicians—because for t 0 the C ( x , 0 ) δ ( x ) .
The main goal of the present study is to evaluate and analyze the solutions for general ( α 1 / 2 ) real parameter values. All the subsequent results are completely new and cannot be found elsewhere in the scientific literature. Luckily, the ODE Equation (4) has closed-form analytic values, which are
f ( η ) = η · e η 2 4 D c 1 M 1 α , 3 2 , η 2 4 D + c 2 U 1 α , 3 2 , η 2 4 D ,
where M ( · , · , · ) and U ( · , · , · ) are the Kummer functions. As definition consider the series expansion of M
M ( a , b , z ) = 1 + a z b + ( a ) 2 z 2 ( b ) 2 2 ! + + ( a ) n z n ( b ) n n ! ,
with the ( a ) n = a ( a + 1 ) ( a + 2 ) ( a + n 1 ) , ( a ) 0 = 1 which is the so-called rising factorial or Pochhammer’s Symbol [60]. In our present case, b has a fix non-negative integer value, so none of the solutions have poles at b = n . For the Kummer function M, when the parameter a has negative integer numerical values ( a = m ), the solution is reduced to a polynomial of degree m for the variable z. In other cases, a m , we obtain a convergent infinite series for all values of a , b and z. There is a connecting formula between the two Kummer’s functions, U is defined from M via
U ( a , b , z ) = π s i n ( π b ) M [ a , b , z ] Γ [ 1 + a b ] Γ [ b ] z 1 b M [ 1 + a b , 2 b , z ] Γ [ a ] Γ [ 2 b ] ,
where Γ ( a ) is the Gamma function [60]. For exhaustive details (e.g., integral representation, recursion formula), see the NIST Handbook [60]. For non-negative integer α s, an alternative formulations of the result is possible in the form of
f ( η ) = e η 2 4 D c ˜ 1 H 2 α 1 η 2 D + c ˜ 2 · 1 F 1 1 2 α 2 , 1 2 ; η 2 4 D ,
where H a ( η ) is the Hermite polynomial and 1 F 1 ( · , · ; · ) is the hypergeometric function. The first part of the solution, the Hermite polynomials, form a complete orthonormal basis set on the + range with the Gaussian weight function and defines a Hilbert space (note that Hermite polynomials play an extraordinary role in quantum mechanics as the solution of the harmonic oscillator problem [69] which pioneered the way to second quantization or field theory). Unfortunately, the second term in Equation (8) if c 2 0 destroys the orthogonality causing over-completeness, which is an interesting feature, taking out the scope that the functional analysis might have far-reaching consequences, needing further in-depth investigation. We have an intuition that such a property could also help us say something new about turbulence. However, that is out of the scope of the present study, but is in line with our long-range scientific interest.
It is important to emphasize that there are four parameter ranges where the derived solutions behave qualitatively differently:
  • α < 0 , where the derived solutions are divergent for large η s, such solutions are nonphysical and out of the scope of our present analysis
  • α = 0 , the solution is zero in the origin and has an asymptotic positive finite value at asymptotic large η
  • 0 < α < 1 , the solution is zero in the origin has a local maxima and a decay to zero as η goes to infinity
  • 1 < α , the solutions are again zero in the origin then have a local maxima and a quick oscillatory decay to zero, larger α values mean more oscillations with more and more zero transitions.
We present here the form of f for the cases of α = 0 , 1 , 2 , 3 , 4 ; the last two cases are evaluated for the first time:
f ( η ) = e r f η 2 D , f ( η ) = κ 0 · η · e η 2 4 D , f ( η ) = κ 0 · η · e η 2 4 D · 1 1 6 D η 2 , f ( η ) = κ 0 · η · e η 2 4 D · 1 1 3 D η 2 + 1 60 1 D 2 η 4 , f ( η ) = κ 0 · η · e η 2 4 D · 1 1 2 D η 2 + 1 20 1 D 2 η 4 1 840 1 D 3 η 6 ,
the κ 0 is an arbitrary normalization constant and will be given later on. For completeness, the final concentration distributions are also provided; by inserting η = x / t 1 / 2 and the actual value of α , we obtain:
C ( x , t ) = e r f x 2 D t , C ( x , t ) = κ 1 x t 3 2 e x 2 4 D t , C ( x , t ) = κ 1 x t 5 2 e x 2 4 D t 1 x 2 6 D t , C ( x , t ) = κ 1 x t 7 2 e x 2 4 D t 1 x 2 3 D t + x 4 60 ( D t ) 2 , C ( x , t ) = κ 1 x t 9 2 e x 2 4 D t 1 x 2 2 D t + x 4 20 ( D t ) 2 x 6 840 ( D t ) 3 .
Figure 1 shows the given five shape functions. Functions with α > 0 clearly show a decaying and oscillatory behavior in terms of η . Figure 2 shows six C ( x , t ) s for different α s; for generality, we also show two solutions for half-integer α s. The quick decay and the slight oscillations are clear to see in all cases. Due to the linearity of the diffusion equation, any linear combination of Equation (5) is automatically a solution, also enriching the possible mathematical structure of the diffusion process.
To complete our investigation, we have to analyse which kind of initial and boundary value problems can be satisfied by these solutions. Our answer is twofold: Firstly, if one of the above solutions (10) is explicitly given, with a fixed numerical value of c 1 ,   c 2 ,   D and κ i , then we may fix an arbitrary time point as t 0 = t , and the corresponding C ( x ) curve can be defined as the initial problem. This curve is presented in Figure 2 for α = 1 . We focus on physically relevant boundary values, which mean, asymptotically for x ± , the solutions C ( ) 0 , to which one may arrive if α > 0 .
Secondly, even the general initial value problem can be handled with the help of the Green’s functions formalism. According to the standard theory of Green’s functions the solution of the diffusion Equation (1) can be obtained via the following convolution integral
C ( x , t ) = 1 π t + w ( x 0 ) G ( x x 0 ) d x 0
where ( x 0 ) contains the initial condition of the problem, C | t = 0 = w ( x 0 ) .
The Green function for diffusion is well-known and can be found in numerous mathematical textbooks, such as [70,71,72,73],
G ( x x 0 ) = e x p ( x x 0 ) 2 4 t D ,
which is the fundamental solution of diffusion.
It is now straightforward to use our new functions (10) in the convolution integral as a Green’s function to derive solutions for arbitrary initial conditions C | t = 0 = w ( x 0 ) . So, in this sense, we can give the most general Green’s functions
C ( x , t ) = κ n t α + w ( x 0 ) · ( x x 0 ) t 1 / 2 · e ( x x 0 ) 2 4 D t c 1 M 1 α , 3 2 , ( x x 0 ) 2 4 D t + c 2 U 1 α , 3 2 , ( x x 0 ) 2 4 D t d x 0
where κ n stands for the proper normalization. We see in the following that for some special forms of the initial conditions, such as polynomials, Gaussian, Sinus, or Cosines, the convolution integral can be conducted analytically. It is clear that, for α = 1 / 2 , we automatically obtain back the Gaussian Green’s function of (12).
Note that for positive integer values of α both Kummer’s U and Kummer’s M functions contain the same polynomials the only difference is just an overall normalization constant. It can be shown, as it was visualized in our last study [1], that for negative α s the shape functions f ( η ) are divergent for large η s. In the regime of 0 α 1 the solutions are similar to the “usual Gaussian” solutions, positive on the whole axis and goes to zero at infinite η . These solutions can be understood as different probability distribution functions and the corresponding expectation values, variance, higher moments, skewness, kurtosis, and other probabilities can be evaluated.
For α > 1 , the shape functions show some oscillations have finite negative values at some arguments and have a quick decay at large arguments η . The norm of the functions are still finite but, due to the oscillations such functions, cannot be interpreted as probability distributions.
To present a solution for a given initial value w ( x 0 ) , the numerical values of the normalization constants κ α have to be given. We apply the L 1 normalization with the definition of
1 = κ α 0 η · e η 2 4 D M 1 α , 3 2 , η 2 4 D d η .
arriving to the numerical values of: κ 1 , 2 , 3 , 4 = ( 1 2 , 3 2 , 5 2 , 7 2 ) . In some cases, the L 2 normalizations have to be used with the coefficients of κ 1 , 2 , 3 , 4 = 2 2 π 4 , 1 5 120 2 π 4 , 1 21 21 · 160 2 π 4 , 1 429 429 · 4480 2 π 4
As a rather physical example, we study the α = 1 ,   c 1 = 1 ,   c 2 = 0 case for the initial condition of
w ( x 0 ) = Heaviside ( x 0 1 ) · Heaviside ( 11 x 0 ) · ( [ x 0 6 ] 2 + 25 ) 25
which is an up-shifted upside down parabola cut at the x 0 = + 1 , + 11 visualized in Figure 3. Note that the w ( x 0 ) = 0 and the w ( x 0 = + ) = 0 boundary conditions are fulfilled. The solution obtained with the convolution integral of (13)
C ( x , t ) = 1 t · 2 t 25 π ( x 6 ) e r f x 1 2 t e r f x 11 2 t + 4 25 t 3 2 e 1 4 t x 1 2 e 1 4 t x 11 2
For a better understanding, this result is presented in Figure 3. Note that the solution function C ( x , t ) fulfills the C ( x 0 ) = 0 condition and the C ( x 0 + ) = 0 asymptotic condition as well. The convolution integral of (13) results a single-valued function of (16) therefore the derived solution is unique.

2.2. An Interesting Ansatz

To have a mathematically complete analysis of the possible Ansatz we have to investigate the “inverse self-similar Ansatz” in the form of:
C ( x , t ) = x α g t x β = x α g ( ω ) .
Now, the role of the temporal and spatial variables is interchanged. The physical interpretation of this new trial function is hard to see; until now, we cannot find any kind of reasonable physical explanation of this new Ansatz Equation (17). (To avoid confusion with the original Ansatz, we use the g ( ω ) notation for this case.) After having conducted the usual derivation and algebraic steps, we arrive at the relations of
α = arbitrary real number , β = 2 .
The obtained ODE looks similar but obviously contains more terms then the previous ones
4 D ω 2 g + ω g [ 2 α D 4 D + 2 ( α 1 ) D ] g + D α ( α 1 ) g = 0 .
The solutions for the shape functions can be evaluated with the help of the usual Kummer’s and exponential functions in the form of
g = c 1 e 1 4 D ω ω 5 4 25 4 α 2 + 4 α 4 M 9 4 + 25 4 α 2 + 4 α 4 , 25 4 α 2 + 4 α 2 , 1 4 D ω + c 2 e 1 4 D ω ω 5 4 25 4 α 2 + 4 α 4 U 9 4 + 25 4 α 2 + 4 α 4 , 25 4 α 2 + 4 α 2 , 1 4 D ω .
The two parameters of the Kummer functions should be real therefore α must lie in the interval of [ 1 2 26 2 , 1 2 26 2 ] which is approximately 2.1 < α < 3.1 . Our experience showed, that basically for any numerical α values the shape functions have a power-law dependence such as f ( η ) η n where 0 < n < 1 and the C ( x , t ) s are divergent at large x arguments. Therefore, we found no physically reasonable solutions, therefore present no figures for Equation (20).

2.3. A Generalization

At this point, it is straightforward to try the generalized form of the self-similar Ansatz
C ( x , t ) = a ( t ) · h x b ( t ) = a ( t ) · h ( ω ) ,
where all a , b and h are continuous real functions with existing continuous first temporal and second spatial derivatives and ω is the new reduced independent variable. Note that now the functions which are responsible for the time decay and spreading have a general form. Instead of the power law dependencies t α and t β , we apply a ( t ) and b ( t ) . By calculating the needed temporal and spatial derivatives and plugging back to the original diffusion equation, we arrive to the ODE of
a t h a b t b ω h = D a b 2 h ,
where prime means derivation in respect to ω and subscript t in respect to time. This equation should be an ODE for h ( ω ) ; therefore, the coefficients of h and h should be independent of time therefore should be equal to constants, so Equation (22) become:
h ω h = D h ,
It has the solution of
h ( ω ) = C 1 ω + C 2 2 D e ω 2 2 D + 2 π D ω · e r f 2 ω 2 D .
The equations of constraints have to be fulfilled as well:
b 2 a t = a · b · b t = a , b 2 0 .
The corresponding solutions can be easily obtained by direct integration and read:
a ( t ) = ± 2 t + c 1 , b ( t ) = c 2 · 2 t + c 1 .
Using the original definitions of the Ansatz (21) we can obtain the final solution in the form of:
C ( x , t ) = ( ± 2 t + c 1 ) · C 1 x c 2 · 2 t + c 1 + C 2 2 D e x 2 2 D c 2 2 ( 2 D t + c 1 ) + 2 π D · x c 2 · 2 t + c 1 · e r f 2 x 2 c 2 D ( 2 t + c 1 ) .
where c 1 ,   c 2 ,   C 1 ,   C 2 are integration constants. Choosing C 1 = c 1 = 0 , we obtain back the usual solution which is a sum of a Gaussian and and error function. Note that this is equivalent to the self-similar solution where α = β = 1 / 2 . It is instructive to see that a more general form of the Ansatz does not necessarily lead to a larger class of solutions. The functions a ( t ) and b ( t ) do not have additional freedom, the power laws however have two free parameters— α and β two real numbers—which expand the class of possible solutions. To the best of our best knowledge, this relatively simple derivation is not yet published or widely known in the scientific community. At this point, we have to note that, in our former study [1], we investigated the traveling-profile Ansatz from [32], which interpolates between the traveling wave and the self-similar trial functions in the form of C ( x , t ) = a ( t ) · h ( [ x b ( t ) / c ( t ) ] = a ( f ) h ( ω ) , where a ( t ) , b ( t ) and c ( t ) are arbitrary continuous functions, with existing first derivatives. The derived results were very similar to the Equation (5).

2.4. A Redefinition of Variables

In the following, we use the conjecture that the f function can be written as:
f ( η ) = η e η 2 4 D g ( η ) .
It is worth checking this form to derive possible new results. The derivative of the function f ( η ) is:
f ( η ) = e η 2 4 D g ( η ) η η 2 D e η 2 4 D + η e η 2 4 D g ( η ) .
The second derivative of function f reads as follows:
f ( η ) = e η 2 4 D η 2 D g ( η ) + g ( η ) η D g ( η ) + η η 2 4 D 2 g ( η ) η 2 2 D g ( η ) + g ( η ) η 2 2 D g ( η ) + η g ( η ) .
Inserting these functions into the equation, while keeping in mind that β = 1 2
α f 1 2 η f = D f ,
we obtain for g = g ( η )
α η g = 2 D g ( η ) η g η 2 2 g + η g D .
Reordering the terms leads to
η g + 2 g η 2 2 D g + ( α 1 ) η D g = 0 .
The solutions read
g ( η ) = c 1 M 1 α , 3 2 , η 2 4 D + c 2 U 1 α , 3 2 , η 2 4 D ,
note that these are the same solutions as for the self-similar Ansatz, only in a separated form. These two examples clearly show that there is a relatively large freedom to define an Ansatz, but only few of them lead to reasonable new solutions.

2.5. Using Various Series Expansions of f ( η )

As a possible generalization of the self-similar Ansatz, we may define the following infinite power series of
C ( x , t ) = i = 1 a i · t α i · ( f [ η ] ) i ,
where a i s and α i s are arbitrary real numbers. As first (and most logical case) just take the following two terms of:
C ( x , t ) = a t α f x t β + b t α f x t β 2 = a t α f ( η ) + b t α f ( η ) 2 ,
the role of α and β is still the same, and the role of a and b are to fix the ratio of the two components or—more importantly—to turn-on or turn-off one of them. It is clear that if we want to derive an ODE for the shape function the argument η should remain on the same first power.
After some relatively simple algebraic steps we got the usual constraints for the two exponents
α = arbitrary real number , β = 1 / 2 .
The derived ODE is
a α f a 2 η f b ( α f 2 η f f ) = D ( a f + 2 b [ f 2 + f f ] ) ,
for general α , the solution is cumbersome containing large number of Kummer’s M and Kummer’s U functions and given in the Appendix A at the end of the paper. However, for some given small values α = 1 ,   ± 1 / 2 ,   ± 1 ,   3 / 2 ,   5 / 2 , the results can be expressed with the help of Gaussians and with the error function. The overall and complete function test of the solutions of Equation (38) is a hard question due to the five parameters of { D , a , b , c 1 , c 2 } ( c 1 and c 2 stand for the integration constants).
For α = 1 / 2 and for arbitrary other real parameters the result reads the follows:
f = a e η 2 4 D 1 D ± e η 2 4 D a 2 e η 2 4 D 4 b c 1 D π D + 4 D b c 2 π D e r f 1 2 1 D η / D 2 b e η 2 4 D 1 D ,
where erf is the usual error function. Note that this solution let ( a = 0 & b ϵ R ) , but not the opposite case. Therefore, it is impossible to obtain back the fundamental or Gaussian solution.
Figure 4 shows the shape functions of Equation (39) for a = 2 ,   b = 1 and for a = 0 ,   b = 1 and the corresponding final C ( x , t ) s as well. Note the remarkable new feature when the f 2 ( η ) term is considered alone; the solutions have a compact support. In this sense, a linear PDE is reduced with a non-linear Ansatz to a non-linear ODE having non-linear properties. This solution is definitely unknown in the scientific community. We might begin to speculate about that even the regular diffusion equation could describe non-regular diffusion phenomena such as the porous media equation [20].
If we consider the quadratic term in (36) alone, we obtain a solution which is just a little more simplified than the one provided in the Appendix A. Therefore, we skip to presenting it.
After this idea, we may go a bit further, considering additional generalized forms such as
C ( x , t ) = t α f x t β + i = 1 a i · t α i · η i · ( f [ η ] ) i ,
Keeping only the first two terms, we arrive to
C ( x , t ) = t α f x t β + a t α η f x t β = t α f ( η ) + a t α η f ( η ) ,
which has the solutions of
f ( η ) = e η 2 4 D c 1 M 1 α , 3 2 , η 2 4 D η + c 2 U 1 α , 3 2 , η 2 4 D η a + b η .
Note that, for a , b > 0 , the only change is only a little different in terms of the scaling of the results. However, if the a · b < 0 , the solution has an obvious singularity, which makes it interesting, but nonphysical; therefore, for now, we ignore this. Last along this line of thought, we may try an Ansatz with the shape of:
C ( x , t ) = a t α f x t β + b t α f x t β 1 / 2 ,
which, unfortunately, provides no analytic solutions. Generally, higher order terms in power expansions provide higher degree non-linear second-oder ODEs, which rarely have analytic solutions.

2.6. Arbitrary Self-Similar Exponents

Lastly, we arrived at a question which leads us out of the problem of regular diffusion. We might ask (if nothing else, only for the sake of completeness): what does it mean when both α , β are arbitrary real numbers dictating the ODE of
α f β η f = D f .
where D is still the usual diffusion coefficient. The solutions remain similar
f ( η ) = η e β η 2 2 D c 1 M 2 β α 2 β , 3 2 , β η 2 2 D + c 2 U 2 β α 2 β , 3 2 , β η 2 2 D .
With some easy reasoning, we can find out the the original form of PDE. An additional t 2 β 1 time dependence has to be included to cancel the usual β = 1 / 2 constraint. Therefore, the starting PDE reads
C ( x , t ) t = D · t 2 β 1 · 2 C ( x , t ) x 2 .
The role of the α parameter—which is responsible for the temporal decay—could be easily investigated, as we saw above; however, now the role of the β —which is responsible for the spreading—was hidden until now. The solution functions are even functions; therefore, we only concentrate on positive arguments and investigate the Kummer M function only. From the formula of (45) two conditions are easy to notice. The first one is that β 0 this is due to the denominator of the first parameter in the Kummer function, and the second one is that for β < 0 the exponential multiplier function goes to infinity at large arguments. We ignore such nonphysical solutions.
The next figure, Figure 5a, presents various f ( η ) shape functions for α = 1 . Three cases can be distinguished:
  • When 2 β < α , the functions have zero transitions and show oscillatory behavior.
  • When 2 β = α , Kummer’s functions are equal to unity; hence, the solution is purely Gaussian, with the quickest possible decay to zero.
  • When 2 β > α , the larger the beta the lower the global maximum and the slower the decay at large argument. The numerical value of α is irrelevant if it is positive. It is interesting that for negative α s and for positive β s the total solution is again divergent at large arguments. For completeness, we show the C ( x , t ) for α = 1 and β = 1 on Figure 5b.
One can see, that for sufficiently large values of β , the shape function f has a maximum, which is followed by a relatively slow decay. For β = 1 / 4 , there is one root of the shape function f ( η ) , and correspondingly of the function C ( x , t ) . This shows that for values β smaller than 1 / 2 , the existence of nontrivial fluctuations in the value of C ( x , t ) is possible, which means that it may become smaller than the average of the background. If β is smaller than 1 / 4 , more roots of the shape function f ( η ) are possible, on Figure 6. These values of β shows a special behavior of the system, where the nontrivial diffusive effects may lead to temporal mass concentrations in case of mass diffusion.

3. Summary and Outlook

We investigated the regular diffusion equation with the usual self-similar Ansatz and discussed all the solutions which are achievable beyond the Gaussian one. Most of the solutions can be described with the help of Kummer’s M and Kummer’s U functions. In certain cases, the linear combination of these functions, which lead to solutions, are explicitly evaluated. For some special parameters, the solutions go over the Hermite polynomials. In the second part of the study, we presented additional solutions which all can be derived from different modifications of the original self-similar Ansatz. Such solutions are again far from being well-known among the scientific community and therefore have to be published and discussed in detail. At the end of our manuscript, we investigated a special diffusion process which has time-dependent diffusion coefficient. This time dependence may be associated to the parameter of spreading, and certain analytic forms can be obtained. Work is in progress to analyze spatial and temporal dependent diffusion equations, which will be the topic of out next study. In the long run, we would like to also investigate reaction–diffusion equations (even with non-constant diffusion coefficients). Our experience so far suggests that numerous existing models can have new solutions with interesting features.

Author Contributions

The corresponding author I.F.B. provided the original idea of the study, performed the analytic calculations, created the figures, and wrote some parts of the manuscript. L.M. evaluated certain formulas, checked all the analytic calculations, checked the spelling, improved the language of the final manuscript, and also collected some part of the cited references. The authors discussed the manuscript on a regular weekly basis. All authors have read and agreed to the published version of the manuscript.

Funding

I.F. Barna was supported by the NKFIH, the Hungarian National Research Development and Innovation Office.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that supports the findings of this study are all available within the article. At the request of readers, the appropriate Maple files are available.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

To be complete we give the exact form of the solution of Equation (38) for arbitrary α . For a better transparency we introduce the following four abbreviation:
M 1 α : = M 1 α , 3 2 , η 2 4 D , U 1 α : = U 1 α , 3 2 , η 2 4 D ,
and similary
M α : = M α , 3 2 , η 2 4 D , U α : = U α , 3 2 , η 2 4 D .
The solution formula is very elaborate but contains these four Kummer’s functions only, therefore the notation is applicable. Note that due to second degree of the ODE, two solutions exist:
f ( η ) = 1 2 2 a M 1 α U α + a M α U 1 α + 2 a α M α U 1 α ± 4 a 2 M 1 α U α 2 + M α M 1 α U α U 1 α { 4 a 2 + 8 a 2 α } + M α U 1 α { a 2 + 4 a 2 α } + 4 a 2 α 2 M α U 1 α + 8 c 2 b α M α U 1 α 2 8 c 1 b α M α M 1 α U 1 α + 8 c 2 b M 1 α U α U 1 α 8 c 1 b M 1 α U α + 4 c 2 b M α U 1 α 4 b c 1 M α M 1 α U 1 α 1 2 / ( b [ 2 α M α U 1 α + 2 M 1 α U α + M α U 1 α ] ) .

References

  1. Mátyás, L.; Barna, I.F. General Self-Similar Solutions of Diffusion Equation and Related Constructions. Rom. J. Phys. 2022, 67, 101–120. [Google Scholar]
  2. Ghez, R. Diffusion Phenomena; Dover Publication Inc.: Long Island, NY, USA, 2001. [Google Scholar]
  3. Newman, J.; Battaglia, V. The Newman Lectures on Transport Phenomena; Jenny Stanford Publishing: Singapore, 2021. [Google Scholar]
  4. Gillespie, D.T.; Seitaridou, E. Simple Brownian Diffusion; Oxford University Press: Oxford, UK, 2013. [Google Scholar]
  5. Bálint, P.; Gilbert, T.; Szasz, D.; Tóth, I. What mathematical billiards teach us about statistical physics? arXiv 2020, arXiv:2009.06284. [Google Scholar]
  6. Deng, W.; Hou, R.; Wang, W.; Xu, P. Modeling Anomalous Diffusion, From Statistics to Mathematics; World Scientific: Singapore, 2020. [Google Scholar]
  7. Uchaikin, V.V. Self-similar anomalous diffusion and Levy-stable laws. Phys. Uspekhi 2003, 46, 821–849. [Google Scholar] [CrossRef]
  8. Arapostathis, A.; Borkar, V.S.; Ghosh, M.K. Ergodic Control of Diffusion Processes; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  9. Divinski, S.V.; Aloke, P. Handbook of Solid State Diffusion; Elsevier: Amsterdam, The Netherlands, 2017. [Google Scholar]
  10. Dybkov, V.I. Reaction Diffusion and Solid State Chemical Kinetics; The IPMS Publications: Kyev, Ukraine, 2002. [Google Scholar]
  11. Shewmon, P. Diffusion in Solids; Springer: New York, NY, USA, 2016. [Google Scholar]
  12. Öchsner, A.; Murch, G.E.; Belova, I.E. Advanced Diffusion Processes and Phenomena; Trans Tech Publication: Stafa-Zurich, Switzerland, 2014. [Google Scholar]
  13. Martin, G. Atomic mobility in Cahn’s diffusion model. Phys. Rev. B 1990, 41, 2279–2283. [Google Scholar] [CrossRef]
  14. Beke, D.L.; Erdélyi, Z. Resolution of the diffusional paradox predicting infinitely fast kinetics on the nanoscale. Phys. Rev. B 2006, 73, 035426. [Google Scholar] [CrossRef]
  15. Balogh, Z.; Erdélyi, Z.; Beke, D.L.; Langer, G.A.; Csik, A.; Boyen, H.G.; Wiedwald, U.; Ziemann, P.; Portavoce, A.; Girardeaux, C. Transition from anomalous kinetics toward Fickian diffusion for Si dissolution into amorphous Ge. Appl. Phys. Lett. 2008, 92, 143104. [Google Scholar] [CrossRef]
  16. Fisher, D.J. Defects and Diffusion in Carbon Nanotubes; Defect and Diffusion Forum, Trans Tech Publication: Stafa-Zurich, Switzerland, 2014. [Google Scholar]
  17. Pelleg, J. Diffusion in Ceramics; Solid Mechanics and Its Applications; Springer: New York, NY, USA, 2016. [Google Scholar]
  18. Michaud, G.; Alecian, G.; Richer, J. Atomic Diffusion in Stars; Astronomy and Astrophysics Library; Springer: New York, NY, USA, 2015. [Google Scholar]
  19. Liehr, A.W. Dissipative Solitons in Reaction Diffusion Systems Volume. In Springer Series in Synergetics; Springer: New York, NY, USA, 2013; Volume 70. [Google Scholar]
  20. Vazquez, J.L. The Porous Medium Equation, Mathematical Theory; Oxford University Press: Oxford, UK, 2007. [Google Scholar]
  21. Galea, E.; Markatos, N. A review of mathematical modelling of aircraft cabin fires. Appl. Math. Model. 1987, 11, 162–176. [Google Scholar] [CrossRef]
  22. Díaz, J.L. Modeling of an aircraft fire extinguishing process with a porous medium equation. SN Appl. Sci. 2020, 2, 2108. [Google Scholar] [CrossRef]
  23. Iron, D.; Wei, J.; Winter, M. Stability analysis of Turing patterns generated by the Schnakenberg model. J. Math. Biol. 2004, 49, 358–390. [Google Scholar] [CrossRef]
  24. Golovin, A.A.; Matkowsky, B.J.; Volpert, V.A. Turing Pattern Formation in the Brusselator Model with Superdiffusion. SIAM J. Appl. Math. 2008, 69, 251–272. [Google Scholar] [CrossRef]
  25. Barna, I.F.; Bognár, G.; Guedda, M.; Hriczó, K.; Mátyás, L. Analytic Self-Similar Solutions of the Kardar–Parisi–Zhang Interface Growing Equation with Various Noise Terms. Math. Model. Anal. 2020, 25, 241–256. [Google Scholar] [CrossRef]
  26. Barna, I.F.; Bognár, G.; Mátyás, L.; Guedda, M.; Hriczó, K. Travelling-wave solutions of the Kardar–Parisi–Zhang interface growing equation with different kind of noise terms. AIP Conf. Proc. 2020, 2293, 280005. [Google Scholar] [CrossRef]
  27. Klages, R.; Barna, I.F.; Mátyás, L. Spiral modes in the diffusion of a granular particle on a vibrating surface. Phys. Lett. A 2004, 333. [Google Scholar] [CrossRef]
  28. Mátyás, L.; Barna, I.F. Geometrical origin of chaoticity in the bouncing ball billiard. Chaos Solitons Fractals 2010, 44, 1111–1116. [Google Scholar] [CrossRef]
  29. Crank, J. The Mathematics of Diffusion; Oxford Clarendon Press: Oxford, UK, 1975. [Google Scholar]
  30. Benett, T. Transport by Advection and Diffusion: Momentum, Heat and Mass Transfer; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar]
  31. Clarkson, P.A.; Kruskal, M.D. New similarity reductions of the Boussinesq equation. J. Math. Phys. 1989, 30, 2201–2213. [Google Scholar] [CrossRef]
  32. Benhamidouche, N. Exact solutions to some nonlinear PDEs, travelling profiles method. Electron. J. Qual. Theory Differ. Equ. 2008, 2008, 1–7. [Google Scholar] [CrossRef]
  33. Sedov, L.I. Similarity and Dimensional Methods in Mechanics; Academic Press: Boca Raton, FL, USA, 1959. [Google Scholar]
  34. Kolmogorov, A.N.; Petrovskii, I.G.; Piskunov, N.S. Selected Works of A. N. Kolmogorov; Kluver Academic Publisher: Norwell, Australia, 1991. [Google Scholar]
  35. Fisher, R.A. The wave of advance of advantageous genes. Ann. Eugen. 1937, 7, 353–369. [Google Scholar] [CrossRef]
  36. Al-Khaled, K. Numerical study of Fisher’s reaction—Diffusion equation by the Sinc collocation method. J. Comput. Appl. Math. 2001, 137, 245–255. [Google Scholar] [CrossRef]
  37. Tálos, K.; Páger, C.; Tonk, S.; Majdik, C.; Kocsis, B.; Kilár, F.; Pernyeszi, T. Cadmium biosorption on native Saccharomyces cerevisiae cells in aqueous suspension. Acta Univ. Sapientiae Agric. Environ. 2009, 1, 20–30. [Google Scholar]
  38. Newell, A.C.; Whitehead, J.A. Finite bandwidth, finite amplitude convection. J. Fluid Mech. 1969, 38, 279–303. [Google Scholar] [CrossRef]
  39. Segel, L.A. Distant side-walls cause slow amplitude modulation of cellular convection. J. Fluid Mech. 1969, 38, 203–224. [Google Scholar] [CrossRef]
  40. Zeldovich, Y.B.; Frank–Kamenetskii, D.A. Distant side-walls cause slow amplitude modulation of cellular convection. Acta Physicochim. 1938, 9, 341. [Google Scholar]
  41. Gliding, B.H.; Kersner, R. Traveling Waves in Nonlinear Diffusion-Convection-Reaction; Progress in Nonlinear Differential Equations and Their Applications; Birkhäuser: Basel, Switzerland, 2001. [Google Scholar]
  42. Nourazar, S.; Soori, M.; Nazari-Golshan, A. On The Exact Solution of Newell-Whitehead-Segel Equation Using the Homotopy Perturbation Method. Aust. J. Basic Appl. Sci. 2011, 5, 1400. [Google Scholar] [CrossRef]
  43. Frank–Kamenetskii, D.A. Calculation of thermal explosion limits. Acta. Phys. Chim USSR 1939, 10, 365. [Google Scholar]
  44. Burgers, J. A Mathematical Model Illustrating the Theory of Turbulence. Adv. Appl. Mech. 1948, 1, 171–199. [Google Scholar] [CrossRef]
  45. Nariboli, G.A.; Lin, W.C. A New Type of Burgers’ Equation. ZAMM—J. Appl. Math. Mech. Z. Angew. Math. Mech. 1973, 53, 505–510. [Google Scholar] [CrossRef]
  46. Sachdev, P.L.; Rao, C.S.; Enflo, B.O. Large-Time Asymptotics for Periodic Solutions of the Modified Burgers Equation. Stud. Appl. Math. 2005, 114, 307–323. [Google Scholar] [CrossRef]
  47. Inan, B. A New Numerical Scheme for the Generalized Huxley Equation. Bull. Math. Sci. Appl. 2016, 16, 105–111. [Google Scholar] [CrossRef]
  48. Ismail, H.N.; Raslan, K.; Rabboh, A.A.A. Adomian decomposition method for Burger’s—Huxley and Burger’s—Fisher equations. Appl. Math. Comput. 2004, 159, 291–301. [Google Scholar] [CrossRef]
  49. Negrea, M.; Petrisor, I.; Constantinescu, D. Aspects of the Diffusion of Electrons and Ions in Tokamak Plasma. Rom. J. Phys. 2010, 55, 1013. [Google Scholar]
  50. Heitjans, P.; Kärger, J. Diffusion in Condensed Matter; Methods, Materials, Models; Springer: New York, NY, USA, 2005. [Google Scholar]
  51. Pressyanov, D.; Dimitrov, D. The problem with temperature dependence of radon diffusion chambers with anti-thoron barrier. Rom. J. Phys. 2020, 65, 801. [Google Scholar]
  52. Jannot, Y.; Bal, H.M.; Degiovanni, A.; Moyne, C. Influence of heat transfer on the estimation of water vapor diffusion coefficient in transient regime. Int. J. Heat Mass Transf. 2021, 177, 121558. [Google Scholar] [CrossRef]
  53. Jiao, D.; Jiao, K.; Du, Q. Numerical investigations of vapor condensation and water transport in gas diffusion layers of PEMFC. Int. J. Heat Mass Transf. 2021, 177, 121543. [Google Scholar] [CrossRef]
  54. Kennedy, C.A.; Carpenter, M.H. Additive Runge–Kutta schemes for convection–diffusion–reaction equations. Appl. Numer. Math. 2003, 44, 139–181. [Google Scholar] [CrossRef]
  55. Sengupta, S.; Sengupta, T.K.; Puttam, J.K.; Vajjala, K.S. Global spectral analysis for convection-diffusion-reaction equation in one and two-dimensions: Effects of numerical anti-diffusion and dispersion. J. Comput. Phys. 2020, 408, 109310. [Google Scholar] [CrossRef]
  56. Jain, A.; Zhou, L.; Parhizi, M. Multilayer one-dimensional Convection-Diffusion-Reaction (CDR) problem: Analytical solution and imaginary eigenvalue analysis. Int. J. Heat Mass Transf. 2021, 177, 121465. [Google Scholar] [CrossRef]
  57. Guo, Z.; Zhao, T.S. Lattice Boltzmann model for incompressible flows through porous media. Phys. Rev. E 2002, 66, 036304. [Google Scholar] [CrossRef]
  58. Mojtabi, A.; Charrier-Mojtabi, M.C. Handbook of Porous Media; Dekker: Abingdon, VA, USA, 2000. [Google Scholar]
  59. Wyss, W. The fractional diffusion equation. J. Math. Phys. 1986, 27, 2782–2785. [Google Scholar] [CrossRef]
  60. Olver, F.W.J.; Lozier, D.W.; Boisvert, R.F.; Clark, C.W. (Eds.) NIST Handbook of Mathematical Functions; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  61. Yang, X.; Balenau, D.; Khan, Y.; Mohyud-Din, S.T. Local fractional variational iteration method for diffusion and wave equations on cantor sets. Rom. J. Phys. 2014, 59, 36. [Google Scholar]
  62. Wang, J.L.; Wu, H.N.; Huang, T.; Ren, S.Y. Analysis and Control of Coupled Neural Networks with Reaction-Diffusion Terms; Methods, Materials, Models; Springer: New York, NY, USA, 2018. [Google Scholar]
  63. Brush, G. The Kind of Motion We Call Heat; Elsevier: Amsterdam, The Netherlands, 1986. [Google Scholar]
  64. Lindhard, J.H.; Lindhard, J.H. A Heat Transfer Textbook; Phlogiston Press: Cambridge, UK, 2003. [Google Scholar]
  65. Jalghaf, H.K.; Kovács, E.; Majár, J.; Nagy, A.; Askar, A.H. Explicit Stable Finite Difference Methods for Diffusion-Reaction Type Equations. Mathematics 2021, 9, 3308. [Google Scholar] [CrossRef]
  66. Nagy, A.; Saleh, M.; Omle, I.; Kareem, H.; Kovács, E. New Stable, Explicit, Shifted-Hopscotch Algorithms for the Heat Equation. Math. Comput. Appl. 2021, 26, 61. [Google Scholar] [CrossRef]
  67. Zeldovich, Y.B.; Raizer, Y.P. Physics of Shock Waves and High-Temperature Hydrdodynamic Phemomena; Academic Press: Cambridge, MA, USA, 1966. [Google Scholar]
  68. Barenblatt, G.I. Scaling; Cambridge Texts in Applied Mathematics; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  69. Messiah, A. Quantum Mechanics; North-Holland Publishing Company: New York, NY, USA, 1961. [Google Scholar]
  70. Rother, T. Green’s Functions in Classical Physics; Lecture Notes in Physics; Springer International Publishing: New York, NY, USA, 2017; Volume 938. [Google Scholar]
  71. Kythe, P.K. Green’s Functions and Linear Differential Equations; Chapman & Hall/CRC Applied Mathrmatics and Nonliner Science; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  72. Greiner, W.; Reinhardt, J. Quantum Electrodynamics; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  73. Bronshtein, I.N.; Semendyayev, K.A.; Musiol, G.; Muehlig, H. Handbook of Mathematics; Springer: Wiesbaden, Germany, 2007. [Google Scholar]
Figure 1. Five evaluated shape functions f ( η ) in Equation (9). The gray, black, blue, red and green curves are for α = 0 , 1 , 2 , 3 and 4, respectively. Additional parameters κ 1 and D are set to unity.
Figure 1. Five evaluated shape functions f ( η ) in Equation (9). The gray, black, blue, red and green curves are for α = 0 , 1 , 2 , 3 and 4, respectively. Additional parameters κ 1 and D are set to unity.
Mathematics 10 03281 g001
Figure 2. The total solutions C ( x , t ) with the shape function of Equation (5) for six various α values. For α = 1 a possible initial condition is given for t = t 0 with the black curve. Additional parameters D = 2 ,   c 2 = 1 ,   c 2 = 0 are the same in all cases. Note that for a better comparison the same ranges are taken for the spatial and temporal variables in all six graphs.
Figure 2. The total solutions C ( x , t ) with the shape function of Equation (5) for six various α values. For α = 1 a possible initial condition is given for t = t 0 with the black curve. Additional parameters D = 2 ,   c 2 = 1 ,   c 2 = 0 are the same in all cases. Note that for a better comparison the same ranges are taken for the spatial and temporal variables in all six graphs.
Mathematics 10 03281 g002
Figure 3. (a) The initial condition of Equation (15), (b) the full solution C ( x , t ) of Equation (16).
Figure 3. (a) The initial condition of Equation (15), (b) the full solution C ( x , t ) of Equation (16).
Mathematics 10 03281 g003
Figure 4. (a) Equation ( 39 ) the black and red lines are for a = 2 ,   b = 1 and for a = 0 ,   b = 1 case for D = c 1 = c 2 = 1 , (b) the full solution C ( x , t ) for a = 2 ,   b = 1 , (c) the full solution C ( x , t ) for a = 0 ,   b = 1 .
Figure 4. (a) Equation ( 39 ) the black and red lines are for a = 2 ,   b = 1 and for a = 0 ,   b = 1 case for D = c 1 = c 2 = 1 , (b) the full solution C ( x , t ) for a = 2 ,   b = 1 , (c) the full solution C ( x , t ) for a = 0 ,   b = 1 .
Mathematics 10 03281 g004
Figure 5. (a) Functions of Equation (45) for α = 1 where black, red, blue, green and brown lines are for β = 1 / 4 ,   1 / 2 ,   1 ,   2 and 3, (b) The C ( x , t ) solution of Equation (46) for α = β = 1 .
Figure 5. (a) Functions of Equation (45) for α = 1 where black, red, blue, green and brown lines are for β = 1 / 4 ,   1 / 2 ,   1 ,   2 and 3, (b) The C ( x , t ) solution of Equation (46) for α = β = 1 .
Mathematics 10 03281 g005
Figure 6. (a) Functions of Equation (45) for α = 1 when β = 1 / 8 , (b) magnification of f ( η ) for values close to zero.
Figure 6. (a) Functions of Equation (45) for α = 1 when β = 1 / 8 , (b) magnification of f ( η ) for values close to zero.
Mathematics 10 03281 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Barna, I.F.; Mátyás, L. Advanced Analytic Self-Similar Solutions of Regular and Irregular Diffusion Equations. Mathematics 2022, 10, 3281. https://doi.org/10.3390/math10183281

AMA Style

Barna IF, Mátyás L. Advanced Analytic Self-Similar Solutions of Regular and Irregular Diffusion Equations. Mathematics. 2022; 10(18):3281. https://doi.org/10.3390/math10183281

Chicago/Turabian Style

Barna, Imre Ferenc, and László Mátyás. 2022. "Advanced Analytic Self-Similar Solutions of Regular and Irregular Diffusion Equations" Mathematics 10, no. 18: 3281. https://doi.org/10.3390/math10183281

APA Style

Barna, I. F., & Mátyás, L. (2022). Advanced Analytic Self-Similar Solutions of Regular and Irregular Diffusion Equations. Mathematics, 10(18), 3281. https://doi.org/10.3390/math10183281

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