Abstract

A semi-infinite material under a solidification process with the Solomon-Wilson-Alexiades mushy zone model with a heat flux condition at the fixed boundary is considered. The associated free boundary problem is overspecified through a convective boundary condition with the aim of the simultaneous determination of the temperature, the two free boundaries of the mushy zone and one thermal coefficient among the latent heat by unit mass, the thermal conductivity, the mass density, the specific heat, and the two coefficients that characterize the mushy zone, when the unknown thermal coefficient is supposed to be constant. Bulk temperature and coefficients which characterize the heat flux and the heat transfer at the boundary are assumed to be determined experimentally. Explicit formulae for the unknowns are given for the resulting six phase-change problems, besides necessary and sufficient conditions on data in order to obtain them. In addition, relationship between the phase-change process solved in this paper and an analogous process overspecified by a temperature boundary condition is presented, and this second problem is solved by considering a large heat transfer coefficient at the boundary in the problem with the convective boundary condition. Formulae for the unknown thermal coefficients corresponding to both problems are summarized in two tables.

1. Introduction

Heat transfer problems with a phase-change such as melting and freezing have been studied in the last century due to their wide scientific and technological applications. Some books in the subject are [19].

In this paper we consider a phase-change process for a semi-infinite material, which is characterized by , that is initially assumed to be liquid at its melting temperature (which without loss of generality we assume to be equal to C). We consider this material under a solidification process with the presence of a zone where solid and liquid coexist, known as “mushy zone,” with a heat flux boundary condition imposed at the fixed face . We follow [10, 11] in considering three different regions in this type of solidification process:(1)Liquid region at temperature : .(2)Solid region at temperature : .(3)Mushy region at temperature : . and are the functions that characterize the free boundaries of the mushy zone. We also follow [10] in making the following assumptions on the structure of the mushy zone, which is considered as isothermal:(1)The material contains a fixed portion of the total latent heat per unit mass (see condition (3)).(2)Its width is inversely proportional to the gradient of temperature (see condition (4)).Thermal coefficients involved in the solidification process are assumed to be constant. They are: latent heat by unit mass,: thermal conductivity,: mass density,: specific heat,: one of the two coefficients which characterize the mushy zone,: one of the two coefficients which characterize the mushy zone,: coefficient that characterizes the heat flux at ,: coefficient that characterizes the heat transfer at ,: bulk temperature at .

We suppose that five of the six thermal coefficients , , , , , and of the solid phase are known and that, by means of a change of phase experiment (solidification of the material at its melting temperature), we are able to measure the quantities , , and .

Encouraged by the recent works [12, 13] and with the aim of the simultaneous determination of temperature , the two free boundaries and , and one unknown thermal coefficient among , , , , , and , we impose an overspecified boundary condition [2] which consists of the specification of a convective condition at the fixed face (see condition (7)) of the material undergoing the phase-change process. This leads us to the following free boundary problem:

This problem was first studied in [11] with a temperature boundary condition at instead of the convective condition (7) considered in this paper. Moreover, the determination of one unknown thermal coefficient for the one-phase Lamé-Clapeyron-Stefan problem with an overspecified heat flux condition at the fixed face without a mushy zone was done in [14]. Other papers related to determination of thermal coefficients are [1540].

The goal of this paper is to obtain the explicit solution to the phase-change process (1)–(7) with one unknown thermal coefficient independent of position and time and the necessary and sufficient conditions on data in order to obtain an explicit formula for the unknown thermal coefficient. In addition, we are interested in analysing the relationship between problem (1)–(7) and the phase-change process given by (1)–(6) besides the Dirichlet boundary condition overspecified at given by (31) (see below). In particular, we are interested in solving the problem with Dirichlet boundary condition through problem with convective boundary condition when large values of the coefficient that characterizes the heat transfer at are considered.

The organization of the paper is as follows. In Section 2 we prove a preliminary result where necessary and sufficient conditions on data for the phase-change process (1)–(7) are given in order to obtain the temperature and the two free boundaries and . Based on this preliminary result, in Section 3, we present and solve six different cases for the phase-change process (1)–(7) according to the choice of the unknown thermal coefficient among , , , , , and . In Section 4 we discuss the relationship between the phase-change process (1)–(6) with the Dirichlet boundary condition (31) and the same process with the convective boundary condition (7). We show that temperature , free boundaries and , and the explicit formula for the unknown thermal coefficient , , , , , or for the phase-change process (1)–(6) with the Dirichlet condition (31) can be obtained through the phase-change process with convective condition given by (1)–(7) when tends to . Explicit formulae for the unknown thermal coefficient for problems (1)–(7) and (1)–(6) and (31), besides restrictions on data that guarantees their validity, are summarized in Tables 1 and 2, respectively.

2. Explicit Solution to the Phase-Change Process

The following lemma represents the base on which the work in this section will be structured.

Lemma 1. The solution to problem (1)–(7) is given byif and only if the parameters involved in problem (1)–(7) satisfy the following two equations:where represents the thermal diffusivity.

Proof. The kind of phase-change processes considered in this paper has the following general solution [1012]:where coefficients , , , and depend on the particular specifications of the phase-change process.
In order to have the solution to problem (1)–(7), we impose conditions (2)–(4), (6) and (7) on (11) and obtain that coefficients , , and must be given bywhich corresponds to solution (8), and the parameters involved in the problem must satisfy (9) and (10).

As a consequence of Lemma 1, we know that we can solve the phase-change process (1)–(7) with one unknown thermal coefficient through the determination of parameter that characterizes one of the two free boundaries of the mushy zone and the unknown thermal coefficient among , , , , , and . In addition, we also know from Lemma 1 that we can do that by solving the system of (9)-(10).

3. Explicit Formula for the Unknown Thermal Coefficient

In this section we present and solve six different cases for the phase-change process (1)–(7) according to the choice of the unknown thermal coefficient among , , , , , and .

With the aim of organizing our work, we classify each case by making reference to the coefficients which is necessary to know in order to solve it (see Lemma 1):Case : determination of and .Case : determination of and .Case : determination of and .Case : determination of and .Case : determination of and .Case : determination of and .In addition, with the goal of making our presentation more readable, in the following statements and proofs we introduce several functions. We name these functions with a subscript according the case where they arise.

Theorem 2 (Case : determination of and ). If in problem (1)–(7) one considers the thermal parameter as an unknown, then its solution is given by (8) with and given byif and only if the parameters , , , , , and satisfy the following two inequalities:

Proof. Due to properties of the error function, it follows that a necessary and sufficient condition for the existence and uniqueness of a positive solution to (10) iswhich is equivalent to inequalities and . In that case, the positive solution to (10) is given by (14). Finally, it follows from (9) that is the positive thermal coefficient given by (13).

Theorem 3 (Case : determination of and ). If in problem (1)–(7) one considers the thermal parameter as an unknown, then its solution is given by (8) with given byand given by (14), if and only if the parameters , , , , , , and satisfy inequalities and andwhere the real function is defined by

Proof. As we see in the proof of Theorem 2, a necessary and sufficient condition that guarantees the existence and uniqueness of solution to (10) is that inequalities and hold, and, in that case, the coefficient is given by (14).
On the other hand, it follows from (9) that is given by (16). This coefficient is positive if and only ifwhere is the real function defined in (17). Taking into account the expression of given in (14), we have that inequality (18) is equivalent to inequality .

Theorem 4 (Case : determination of and ). If in problem (1)–(7) one considers the thermal parameter as an unknown, then its solution is given by (8) with given byand given by (14), if and only if the parameters , , , , , , and satisfy inequalities , , and andwhere is the real function defined in (17).

Proof. Conditions and and the expression of given in (14) arise in the same way as that in the precedent proofs. On the other hand, it follows from (9) that is given by (19), being the real function defined in (17). This coefficient is positive if and only ifTaking into account the expression of given in (14), we have that inequality (20) is equivalent to inequality . Finally, we have that given in (19) is less than 1 if and only if , which, as we see in the proof of Theorem 3, is equivalent to condition .

Theorem 5 (Case : determination of and ). If in problem (1)–(7) one considers the thermal parameter as an unknown, then its solution is given by (8) with given byand is the unique solution to the equationwhere is the real function defined byif and only if the parameters , , and satisfy inequality .

Proof. The system of (9)-(10) is equivalent towhere is the real function defined in (23). A necessary condition for existence of solution to this system is that inequality holds. Then, if we assume that holds, we immediately obtain that is given by (21). To complete the proof it only remains to demonstrate that (22) admits a unique positive solution. This follows from the fact that is an increasing function such that and .

Theorem 6 (Case : determination of and ). If in problem (1)–(7) one considers the thermal parameter as an unknown, then its solution is given by (8) with given byand is the unique solution to (22), if and only if the parameters , , and satisfy inequality .

Proof. It is similar to the proof of Theorem 5.

Theorem 7 (Case : determination of and ). If in problem (1)–(7) one considers the thermal parameter as an unknown, then its solution is given by (8) with given byand is the unique solution to the equationwhere is the real function defined byif and only if the parameters , , and satisfy inequality and

Proof. The system of (9)-(10) is equivalent towhere is the real function defined in (28). A necessary condition for existence of solution to this system is that inequality holds. Then, if we assume that holds, we immediately obtain that is given by (26). To complete the proof it only remains to demonstrate that (27) admits a unique positive solution. Since is an increasing function such that and , it follows that a necessary and sufficient condition for existence (and uniqueness) of solution to (27) is that which is equivalent to inequality .

Table 1 summarizes the results of this section, corresponding to 6 cases.

4. The Phase-Change Process with Large Heat Transfer Coefficient

A similar phase-change process to (1)–(7) with one unknown thermal coefficient has been studied in [11]. In that paper, the author considers a fusion process with a mushy zone given by (1)–(6) overspecified with a temperature boundary condition and obtains the explicit solution to some cases. Encouraged by [11], we consider the solidification process (1)–(6) with one unknown thermal coefficient overspecified with the Dirichlet boundary condition:We can see this condition as the limit case of the convective boundary condition (7) when the heat transfer coefficient tends to . From a physical point of view, overspecifying the phase-change process (1)–(6) by imposing the convective boundary condition (7) seems to be more appropriate than imposing the Dirichlet condition (31). This section is devoted to showing that the temperature , the free boundaries and , and the explicit formula for the unknown thermal coefficient , , , , , or for the phase-change process with Dirichlet boundary condition given by (1)–(6) and (31) can be obtained through the phase-change process with convective boundary condition given by (1)–(7) when tends to .

We begin with a result related to the solution to the phase-change process (1)–(6) and (31). This result may be shown in much the same manner as Lemma 1; thus we do not give its proof here.

Lemma 8. The solution to problem (1)–(6) and (31) is given by (8); that is,where represents the thermal diffusivity, if and only if the parameters involved in problem (1)–(6) and (31) satisfy the following two equations:

From Lemma 8 we know that, in order to have temperature , the free boundaries and , and the unknown thermal coefficient , , , , , or for problem (1)–(6) and (31), it is enough to find the unknown thermal coefficient and the parameter that characterizes the free boundary . Proceeding analogously to the work done in [11] or in Section 3, we can obtain the thermal coefficient , , , , , or and the parameter for problem (1)–(6) and (31). Formulae for those quantities, besides restrictions on data that guarantee their validity, are summarized in Table 2 (restrictions on data and definitions on functions mentioned in Table 2 are listed as follows).

List of restrictions on data for problem (1)–(6) and (31) mentioned in Table 2 iswhere is the unique positive solution to the equation being the real function defined in (17);where is the unique positive solution to the equation being the real function defined in (17);

List of definitions of functions related to problem (1)–(6) and (31) mentioned in Table 2 is being the real function defined in (17);

On the other hand, it is not difficult to verify that formulae and restrictions on data given in Table 2 correspond to formulae and restrictions on data given in Table 1 for tending to . This fact, besides Lemmas 1 and 8, allows us to conclude that we can solve the phase-change process (1)–(6) with one unknown thermal coefficient overspecified by the Dirichlet condition (31) through the phase-change process (1)–(7), which is overspecified by the more physically appropriate convective boundary condition (7), when the heat transfer coefficient tends to .

5. Conclusions

In this paper, we consider a semi-infinite material under a solidification process with a mushy zone caused by an initial heat flux boundary condition, when the thermophysical parameters involved in the phase-change process are assumed to be constant. We solve the associated free boundary problem overspecified with a convective boundary condition and obtain the temperature, the two free boundaries of the mushy zone and one thermal coefficient among the latent heat by unit mass, the thermal conductivity, the mass density, the specific heat, and the two coefficients that characterize the mushy zone, when the bulk temperature and the coefficients that characterize the heat flux and the heat transfer at the boundary are assumed to be known. As a consequence, we give formulae for the temperature, the two free boundaries, and the unknown thermal coefficient, besides necessary and sufficient conditions on data, in order to obtain them. In addition, we present the relationship between the phase-change process studied in this paper and another similar phase-change process which is overspecified by a Dirichelt boundary condition. From this relationship, we solve the problem with the Dirichlet condition by considering a large heat transfer coefficient in the problem with the convective condition. In this way, we solve the phase-change process overspecified with a temperature boundary condition through the more physically appropriate phase-change problem overspecified with a convective boundary condition. We summarize explicit formulae for the unknown thermal coefficient for both problems in Tables 1 and 2.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This paper has been partially sponsored by the Project PIP no. 0534 from CONICET-UA (Rosario, Argentina) and AFOSR-SOARD Grant FA 9550-14-1-0122. The authors like to thank two anonymous referees for the helpful comments.