Determination of two unknown thermal coefficients through an inverse one-phase fractional Stefan problem
We consider a semi-infinite one-dimensional phase-change material with two unknown constant thermal coefficients among the latent heat per unit mass, the specific heat, the mass density and the thermal conductivity. Aiming at the determination of them, we consider an inverse one-phase Stefan problem with an over-specified condition at the fixed boundary and a known evolution for the moving boundary. We assume that the phase-change process presents latent-heat memory effects by considering a fractional time derivative of order () in the Caputo sense and a sharp front model for the interface. According to the choice of the unknown thermal coefficients, six inverse fractional Stefan problems arise. For each of them, we determine necessary and sufficient conditions on data to obtain the existence and uniqueness of a solution of similarity type. Moreover, we present explicit expressions for the temperature and the unknown thermal coefficients. Finally, we show that the results for the classical statement of this problem, associated with , are obtained through the fractional model when .
Keywords: Inverse Stefan problems, Time fractional Caputo derivative, Unknown thermal coefficients, Explicit solutions, Wright and Mainardi Functions, Over-specified boundary condition.
Mathematics Subject Classification 2010: 26A33, 35C05, 35R11, 35R35, 80A22.
Determination of thermal coefficients for phase-change materials through inverse Stefan problems has been widely studied during the last decades [5, 6, 14, 27, 28, 29]. Especially, phase-change processes involving solidification or melting have been extensively studied because of their scientific and technological applications [1, 2, 4, 7, 8, 10, 18, 26, 31]. A rewiew of a long bibliography on moving and free boundary value problems for the heat equation can be consulted in . Recently, a new sort of Stefan problems including time-fracional derivatives have begun to be studied [13, 15, 33, 34, 35, 23, 24, 25]. Some references in fractional derivatives can be found at [16, 19, 22, 9, 17, 20]. In [23, 24, 25] there are considered free boundary value problems which are obtained by replacing the time derivative in a one-phase Stefan problem by a fractional derivative of order () in the Caputo sense , and explicit solutions of similarity type are given for the resultant fractional Stefan problems. A physical interpretation of the problems considered in [23, 24, 25] is given in . In that article, authors derive fractional Stefan problems for phase-change processes by substituting the local expression of the heat flux given by the Fourier law for a new non-local definition. They consider a heat flux given as a weighted sum of local fluxes back in time, which they express in terms of the Riemann-Liouville integral of order () of the local flux given by the Fourier law. They also explain how this change implies that the new model takes into consideration latent-heat memory effects in the evolution of the phase-change process and give to the parameter the physical meaning of being the strength of memory retention. This fractional model reduces to the classical Stefan problem when . The same occurs with the solutions of similarity type given in [23, 24, 25] in the sense that they converge to the similarity solutions of the classical Stefan problems with which they are related to, when . To the authors knowledge, the first use of inverse fractional Stefan problems for the determination of thermal coefficients has been done recently in . In that article the author studies the determination of one unknown thermal coefficient for a semi-infinite material through a fractional one-phase Stefan problem with an over-specified condition at the fixed boundary. Necessary and sufficient conditions on data to obtain the existence and uniqueness of solutions of similarity type are established, and explicit expressions for the the temperature of the material, the free boundary and the unknown thermal coefficient are given. Moreover, it has shown that results through the fractional model reduce to the results previously obtained in  for the determination of one unknown thermal coefficient using a classical inverse Stefan problem. Encouraged by [28, 32], we consider here the problem of determining two unknown thermal coefficients through an inverse fractional one-phase Stefan problem for which it is known the evolution of the free boundary. In order to have dimensional coherence in the time fractional heat equation as well as in the fractional Stefan condition, we have included two extra parameters in the model, which are such that:
In particular, and can be considered equal to 1 with the corresponding physical dimension (see below).
More precisely, we consider the following inverse problem for a one-phase melting process:
where the unknowns are the temperature [C] of the liquid phase and two thermal coefficients among:
|:||latent heat per unit mass|||
We assume that the remaining coefficients:
|:||temperature at the boundary|||
|:||coefficient characterizing the heat flux at|||
|:||strength of the memory retention||(dimensionless)|
|:||parameter required to have dimensional coherence in equation (2a)|||
|:||parameter required to have dimensional coherence in condition (2c)|||
involved in Problem (2), are all known (e.g. through a phase-change experiment). Aiming at the simultaneous determination of two unknown thermal coefficients, we consider that the time evolution of the sharp interface is also known. More precisely, we follow [23, 24, 25] in assuming that it is given by:
for any , where is the Gamma function defined by:
2 Solutions of similarity type
The triplet given by the temperature and two unknown thermal coefficients among , , and is a solution to the inverse fractional one-phase Stefan problem (2) if:
is defined in .
, where and .
is a continuous function over , where , except maybe in the point , in which it must be satisfied the following condition:
For all , exists .
The unknown thermal coefficients are positive real numbers.
The temperature and the two unknown thermal coefficients verify (2).
According to [23, 24], we know that the function given by (5) fulfill conditons 1 to 4 in Definition 2.1 and that it satisfies the fractional diffusion equation (2a). Noting that , it follows from condition (2d) that:
From this and condition (2b), we have that:
We observe that because is a strictly decreasing function in  and, as we have already noted, .
and the fractional Caputo derivative of a power function with positive exponent is given by :
it follows from the fractional Stefan condition (2c) that:
where is the Mainardi function, which is defined by :
and satisfies for all .
We have thus proved the following result:
3 Existence and uniqueness of solutions of similarity type. Formulae for the two unknown thermal coefficients.
In this section we will look for necessary and sufficient conditions on data to have existence and uniqueness of solution to Problem (2) for each possible choice of the two unknown thermal coefficients, as well as explicit formulae for them. Thanks to Theorem 2.1, it can be done through analysing and solving the system of equations (13) for each pair of unknown thermal coefficients. With the aim of organizing the main results of this section, we will write:
|Case 1: Determination of and||Case 4: Determination of and|
|Case 2: Determination of and||Case 5: Determination of and|
|Case 3: Determination of and||Case 6: Determination of and .|
For each , we introduce the real functions , and defined in by:
The following result will be useful all throughout this section:
For any , the real functions , and defined in (15) satisfy the following conditions:
Theorem 3.1 (Case 1: Determination of and ).
Isolating from equation (13b), we have that is given by (18a). Now, by combining this with equation (13a), it can be obtained that is given by (18b). It must be noted that the parameter involved in both (18a) and (18b) depends on . Nevertheless, it can be determined without making any reference to as follows. By replacing (18a) in the definition of given in (14), we have that must be a solution to equation (19). It follows from (17a) that the equation (19) admits a solution if and only if its RHS is between and . To complete the proof only remains to observe that this is equivalent to say that inequality (20) must hold and that, when this happens, equation (19) has an only one positive solution. ∎
Theorem 3.2 (Case 2: Determination of and ).
By following the same steps as in the demonstration of the Theorem 3.1, it can be shown that and must be given by (21), where the parameter should be a solution to equation (22). Since the Mainardi function is a strictly decreasing function from to in , we have that the equation (22) admits a solution if and only if its RHS is between and . This is equivalent to say that inequality (23) must holds. Moreover, when data satisfy (23), equation (22) has an only one positive solution. ∎
Theorem 3.3 (Case 3: Determination of and ).
By proceeding analogously to the proofs of the previous Theorems 3.1 and 3.2, we have that and should be given by (24), being a solution to equation (25). Since the RHS of the equation (25) is a positive number, it follows from (17b) that the equation (25) admits an only one solution for any set of data. ∎
Theorem 3.4 (Case 4: Determination of and ).
It is similar to the demonstration of Theorem 3.1. ∎
Theorem 3.5 (Case 5: Determination of and ).
It is similar to the demonstration of Theorem 3.1. ∎
Theorem 3.6 (Case 6: Determination of and ).
By following the same ideas as in the proof of the Theorem 3.1, we have that and must be given by (28), where should be a solution to equation (29). By noting that the RHS of this equation is a positive number, it follows from (17c) that the equation (29) admits an only one positive solution for any set of data. ∎
Table 1 summarizes the formulae obtained for the two unknown thermal coefficients and the condition that data must verify to obtain them, for each one of the six possible choices of the two unknown thermal coefficients among , , and in Problem (2) when the moving boundary is defined by (3).
|Case||Formulae for the unknown||Equation for||Condition|
|thermal coefficients||for data|
4 Convergence to the classic case when
When , the time fractional derivative of order in the Caputo sense of a function coincides with its classical time derivative. Then, if we allow to be equal to 1 in Problem (2) and we consider that case, we obtain that Problem (2) reduces to the classical inverse one-phase Stefan problem studied in . This problem, which we will refer to as Problem (2), is given by the classical diffusion equation:
the classical Stefan condition:
and conditions (2b), (2d) and (2e). Of course, to obtain (2a) and (2c) we have also considered and in (2a) and (2c), respectively. According to the physical interpretation given in , as increases its value to 1, the ability of memory retention of the system diminishes to the limit case of no memory retention corresponding to . In this context, classical Stefan problems are representing phase-change processes whose temporal evolution can be described in terms of local in time properties (absence of memory).
The determination of two unknown thermal coefficients through a classical inverse one-phase Stefan problem was done in . In that article, necessary and sufficient conditions on data to obtain existence and uniqueness of solution to Problem (2) are given, together with formulae for the unknown thermal coefficients. In several articles [23, 24, 25, 32] it has been proved the convergence when of the solution to a fractional Stefan problem with to the solution to the associated classical problem obtained by considering . Encouraged by those works, we are interested in this section in proving the convergence when of the results obtained in Section 3 to the ones given in .
In order to emphasize the dependence on of the formulae given in Theorems 3.1 to 3.6, we will mention it here explicitly. For example, if we are analysing the convergence of the solution to Problem (2) given in Theorem 3.1, we will refer to it as , and . We will also write to represent the coefficient defined by (14).
The next result will be useful in the following:
For each , the Wright and Mainardi functions are such that:
where is the error function, which is defined by:
See . ∎
Theorem 4.1 (Convergence related to Case 1).
On one hand, we have that the LHS of the equation (34) defines a strictly decreasing function from to in . On the other hand, since inequality (20) holds, the RHS of equation (20) is between and . Therefore, it follows that equation (34) has an only one positive solution. By introducing the parameter defined by (32), we can rewrite equation (34) as follows:
and see that the solution to the equation (19) is such that:
Theorem 4.2 (Convergence related to Case 2).
Since inequality (23) holds for all , we have that the following inequality also holds:
Theorem 4.3 (Convergence related to Case 3).
Since the LHS of equation (48) defines a strictly increasing function from to in and the RHS of equation (48) is a positive number, it follows that equation (48) has an only one positive solution. By introducing the parameter defined by (32), the equation (48) can be rewritten as:
from which we can see that the solution to the equation (25) is such that:
being the only one solution to equation (47). The rest of the proof runs as before. ∎
The following results can be proved in the same manner than the previous theorems in this section. Then, we have prefered to omit their proofs.