Time-Fractional Diffusion with Mass Absorption in a Half-Line Domain due to Boundary Value of Concentration Varying Harmonically in Time

The time-fractional diffusion equation with mass absorption is studied in a half-line domain under the Dirichlet boundary condition varying harmonically in time. The Caputo derivative is employed. The solution is obtained using the Laplace transform with respect to time and the sin-Fourier transform with respect to the spatial coordinate. The results of numerical calculations are illustrated graphically.


Introduction
From a mathematical point of view, diffusion and heat conduction are described by the same equation of the parabolic type ∂u ∂t = a∆u (1) where a is the diffusivity coefficient, t denotes time, ∆ is the Laplace operator, and u stands for concentration in the case of diffusion and for temperature in the case of heat conduction. Ångström [1] was the first to consider Equation (1) under harmonic impact. In that case, sometimes the terms "oscillatory diffusion" or "diffusion waves" are used [2][3][4]. Introducing of oscillations into the diffusion equation can be done by two ways. The first possibility consists in considering the harmonic source term ∂u(x, t) ∂t = a ∂ 2 u(x, t) Nowacki [5,6] studied the equation in the domain −∞ < x < ∞ with δ(x) being the Dirac delta function under assumption u(x, t) = U(x) e iωt (4) and obtained the solution The square root √ iω is defined as a particular case of the general formula [7] (±iω) α = |ω| α e ±απi sign(ω)/2 .
Another possibility to introduce oscillations in the diffusion equation consists in imposing the harmonic boundary condition. For example, the diffusion equation is considered in the domain 0 < x < ∞ under condition as well as under Assumption (4). The solution has the following form: For the boundary condition the solution becomes [8,9] u(x, t) and for the boundary condition one obtains [10] u( In a medium with a chemical reaction or with heat absorption/release, in Equation (1) , there appears an additional linear term [11,12] The values of the coefficient b > 0 and b < 0 correspond to mass/heat absorption and mass/heat release, respectively. Equation (14) also describes mass or heat transport in a thin plate which lateral surfaces exchange mass or heat with surroundings [13] as well as bio-heat transfer [14][15][16]. In the case of one spatial variable, Equation (14) governs propagation of neuronal signals and is known as the cable equation [17,18].
The time-fractional counterpart of Equation (1) has the form [34][35][36] with Γ(α) being the gamma function. The Caputo derivative has the Laplace transform rule Here, the asterisk denotes the transform, and s is the Laplace transform variable. The equation can be regarded as the time-fractional generalization of the diffusion, bio-heat, and cable equations as well as time-fractional generalization of the Klein-Gordon equation [37][38][39][40][41][42][43].
In this paper, we study Equation (19) in a half-line domain under the Dirichlet boundary condition varying harmonically in time; the paper develops the results of [44].

The Statement and Solution of the Problem
The equation is considered in the domain 0 < x < ∞ under the boundary condition For later use in numerical calculations, it is appropriate to introduce the following nondimensional quantities: First, we investigate two particular cases of Equation (20) corresponding to the integer values of the order α under Assumption (4).

Bio-Heat Equation: Quasi-Steady-State Oscillations
In this case α = 1 and for the function U(x) we have equations and Using the sin-Fourier integral transform, we obtain where the tilde denotes the transform, and ξ is the transform variable. The inverse sin-Fourier transform yields for b > 0 and If the boundary condition is then where Figure 1 shows the dependence of the amplitude of oscillations described by Equation (29)

Klein-Gordon Equation: Quasi-Steady State Oscillations
For the hyperbolic Klein-Gordon equation (α = 2) under the boundary condition and Assumption (4), we obtain and Applying the sin-Fourier transform to Equation (33), we obtain The inversion of the transform depends on the relation between b and ω.
whereas, for b < ω 2 , Figure 2 shows the dependence of the nondimensional amplitude of oscillations described by Equations (36) and (37)   It should be noted the difference between the result of the present paper and the corresponding result of [44]. If Equation (31) is considered with the harmonic source term δ(x) e iωt in the domain −∞ < x < ∞ under Assumption (4), then there appears the resonance at ω = √ b (see Figure 3). In the case of the half-line domain with the harmonic Dirichlet boundary condition, there is no resonance.

Equation with Time-Fractional Derivative
It should be emphasized that the equation with time-fractional derivative under the boundary condition cannot be considered under Assumption (4). This is due to the formula [45] d α e λt dt α = λ α e λt γ(n − α, λt) where γ(a, x) is the incomplete gamma function [46] γ(a, x) = x 0 e −w w a−1 dw.
For equations with time-derivative of the fractional order, the initial conditions should be imposed; for example, and The Laplace transform with respect to time t and the sin-Fourier transform with respect to the spatial coordinate x yield Inversion of the integral transform using the convolution theorem results in the solution where the following formula [7,23] has been used. Here E α,β (z) is the Mittag-Leffler function in two parameters α and β:

Bio-Heat Equation: The Solution with Zero Initial Condition
Equation (44) for α = 1 takes the form The inverse Laplace transform results in u(ξ, t) = u 0 aξ and finally where erfc (x) is the complementary error function.

Klein-Gordon Equation: The Solution with Zero Initial Conditions
For α = 2, Equation (44) becomes The inverse sin-Fourier transform yields The inverse Laplace transform of exp −x √ s 2 + b/ √ a depends on the sign of b and reads [47]: (54) Here, J 1 (z) is the Bessel function of the first kind, I 1 (z) is the modified Bessel function of the first kind.
Using the convolution theorem, we obtain, for b > 0, and, for b < 0,

Conclusions
We have considered the time-fractional diffusion-wave equation with the Caputo fractional derivative of the order 0 < α ≤ 2 with mass absorption in a half-line domain under the Dirichlet boundary condition varying harmonically in time. The investigated equation can also be regarded as the time-fractional generalization of the bio-heat and Klein-Gordon equations. The Caputo derivative of the exponential function has more complicated form than the derivative of the integer order. As a consequence, the assumption that the solution can be represented as a product of a function of the spatial coordinate and the time-harmonic term without taking into account the initial conditions cannot be used. In the case of the standard Klein-Gordon equation, the solution of Equations (36) and (37) describes the quasi-steady-state oscillations, whereas the solution of Equations (55) and (56) also describes the transient process and has the wave front at x = √ at. When α approaches 2, the solution approximates this wave front.