## Mathematical Methods for Heat Transfer and Thermodynamic Analysis of Conductive, Convective, and Radiative Media

View this Special IssueResearch Article | Open Access

# Determination of One Unknown Thermal Coefficient through a Mushy Zone Model with a Convective Overspecified Boundary Condition

**Academic Editor:**Mohsen Torabi

#### 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 [1–9].

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 [15–40].

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 by**if 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 [10–12]: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 by**if 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 by**and given by (14), if and only if the parameters , , , , , , and satisfy inequalities and and**where 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 by**and given by (14), if and only if the parameters , , , , , , and satisfy inequalities , , and and**where 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 by**and is the unique solution to the equation**where is the real function defined by**if 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 ). *and 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 ). *and is the unique solution to the equation**where is the real function defined by**if 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.

#### References

- V. Alexiades and A. D. Solomon,
*Mathematical Modeling of Melting and Freezing Processes*, Hemisphere Publishing Corporation, Washington, DC, USA, 1993. - J. R. Cannon,
*The One-Dimensional Heat Equation*, Addison-Wesley, Menlo Park, Calif, USA, 1984. - H. S. Carslaw and J. C. Jaeger,
*Conduction of Heat in Solids*, The Clarendon Press, Oxford, UK, 1959. View at: MathSciNet - J. Crank,
*Free and Moving Boundary Problems*, The Clarendon Press, Oxford, UK, 1984. View at: MathSciNet - A. Fasano,
*Mathematical Models of Some Diffusive Processes with Free Boundaries*, vol. 11 of*MAT Serie A: Conferencias, Seminarios Y Trabajos de Matematica*, 2005. - S. C. Gupta,
*The Classical Stefan Problem. Basic Concepts, Modelling and Analysis*, Elsevier, Amsterdam, The Netherlands, 2003. View at: MathSciNet - J. M. Hill,
*One-Dimensional Stefan Problems: An Introduction*, Longman Higher Education, Essex, UK, 1987. - V. J. Lunardini,
*Heat Transfer with Freezing and Thawing*, Elsevier Science Publisher BV, 1991. - A. B. Tayler,
*Mathematical Models in Applied Mechanics*, Clarendon Press, Oxford, UK, 1986. View at: MathSciNet - A. D. Solomon, D. G. Wilson, and V. Alexiades, “A mushy zone model with an exact solution,”
*Letters in Heat and Mass Transfer*, vol. 9, no. 4, pp. 319–324, 1982. View at: Publisher Site | Google Scholar - D. A. Tarzia, “Determination of unknown thermal coefficients of a semi-infinite material for the one-phase Lamé-Clapeyron (Stefan) problem through the Solomon-Wilson-Alexiades' mushy zone model,”
*International Communications in Heat and Mass Transfer*, vol. 14, no. 2, pp. 219–228, 1987. View at: Publisher Site | Google Scholar - D. A. Tarzia, “Explicit solutions for the Solomon-Wilson-Alexiades's mushy zone model with convective or heat flux boundary conditions,” http://arxiv.org/abs/1503.02736. View at: Google Scholar
- D. Tarzia, “Relationship between Neumann solutions for two-phase Lamé-Clapeyron-Stefan problems with convective and temperature boundary conditions,”
*Thermal Science*, 12 pages, 2015. View at: Publisher Site | Google Scholar - D. A. Tarzia, “Determination of the unknown coefficients in the Lamé-Clapeyron Stefan problem (or one-phase Stefan problem),”
*Advances in Applied Mathematics*, vol. 3, no. 1, pp. 74–82, 1982. View at: Publisher Site | Google Scholar | MathSciNet - V. L. Baranov, A. A. Zasyad'ko, and G. A. Frolov, “Integro-differential method of solving the inverse coefficient heat conduction problem,”
*Journal of Engineering Physics and Thermophysics*, vol. 83, no. 1, pp. 60–71, 2010. View at: Publisher Site | Google Scholar - I. Boulaoued and A. Mhimid, “Determination of the diffusion coefficient of new insulators composed of vegetable fibers,”
*Thermal Science*, vol. 16, no. 4, pp. 987–995, 2012. View at: Publisher Site | Google Scholar - R. Das, S. C. Mishra, and R. Uppaluri, “Retrieval of thermal properties in a transient conduction-radiation problem with variable thermal conductivity,”
*International Journal of Heat and Mass Transfer*, vol. 52, no. 11-12, pp. 2749–2758, 2009. View at: Publisher Site | Google Scholar - Z.-C. Deng, L. Yang, J.-N. Yu, and G.-W. Luo, “An inverse problem of identifying the coefficient in a nonlinear parabolic equation,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 71, no. 12, pp. 6212–6221, 2009. View at: Publisher Site | Google Scholar | MathSciNet - A. Erdem, D. Lesnic, and A. Hasanov, “Identification of a spacewise dependent heat source,”
*Applied Mathematical Modelling*, vol. 37, no. 24, pp. 10231–10244, 2013. View at: Publisher Site | Google Scholar | MathSciNet - N. L. Gol'dman, “One-phase inverse stefan problems with unknown nonlinear sources,”
*Differential Equations*, vol. 49, no. 6, pp. 680–687, 2013. View at: Publisher Site | Google Scholar | MathSciNet - A. Hazanee, M. I. Ismailov, D. Lesnic, and N. B. Kerimov, “An inverse time-dependent source problem for the heat equation,”
*Applied Numerical Mathematics*, vol. 69, pp. 13–33, 2013. View at: Publisher Site | Google Scholar | MathSciNet - A. Hasanov and B. Pektas, “Identification of an unknown time-dependent heat source term from overspecified Dirichlet boundary data by conjugate gradient method,”
*Computers & Mathematics with Applications*, vol. 65, no. 1, pp. 42–57, 2013. View at: Publisher Site | Google Scholar | MathSciNet - E. Hetmaniok, I. Nowak, D. Slota, R. Witula, and A. Zielonka, “Solution of the inverse heat conduction problem with Neumann boundary condition by using the homotopy perturbation method,”
*Thermal Science*, vol. 17, no. 3, pp. 643–650, 2013. View at: Publisher Site | Google Scholar - N. Hryntsiv, “The inverse problem with free boundary for a weakly degenerate parabolic equation,”
*Journal of Mathematical Sciences*, vol. 183, no. 6, pp. 779–795, 2012. View at: Publisher Site | Google Scholar | MathSciNet - M. Hussein and D. Lesnic, “Determination of a time-dependent thermal diffusivity and free boundary in heat conduction,”
*International Communications in Heat and Mass Transfer*, vol. 53, pp. 154–163, 2014. View at: Publisher Site | Google Scholar - M. S. Hussein, D. Lesnic, and M. I. Ivanchov, “Simultaneous determination of time-dependent coefficients in the heat equation,”
*Computers & Mathematics with Applications*, vol. 67, no. 5, pp. 1065–1091, 2014. View at: Publisher Site | Google Scholar | MathSciNet - Y. Inatomi, F. Onishi, K. Nagashio, and K. Kuribayashi, “Density and thermal conductivity measurements for silicon melt by electromagnetic levitation under a static magnetic field,”
*International Journal of Thermophysics*, vol. 28, no. 1, pp. 44–59, 2007. View at: Publisher Site | Google Scholar - F. Kanca and M. I. Ismailov, “The inverse problem of finding the time-dependent diffusion coefficient of the heat equation from integral overdetermination data,”
*Inverse Problems in Science and Engineering*, vol. 20, no. 4, pp. 463–476, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - N. B. Kerimov and M. I. Ismailov, “An inverse coefficient problem for the heat equation in the case of nonlocal boundary conditions,”
*Journal of Mathematical Analysis and Applications*, vol. 396, no. 2, pp. 546–554, 2012. View at: Publisher Site | Google Scholar | MathSciNet - M. Lamvik and J. M. Zhou, “A novel technique for measuring the thermal conductivity of metallic materials during melting and solidification,”
*Measurement Science and Technology*, vol. 6, no. 7, pp. 880–887, 1995. View at: Publisher Site | Google Scholar - G. Li, D. Yao, Y. Wang, and X. Jia, “A nonlinear solute transport model and data reconstruction with parameter determination in an undisturbed soil-column experiment,”
*Mathematical Problems in Engineering*, vol. 2011, Article ID 679531, 14 pages, 2011. View at: Publisher Site | Google Scholar - T. Lu, W. Han, P. Jiang, Y. Zhu, J. Wu, and C. Liu, “A two-dimensional inverse heat conduction problem for simultaneous estimation of heat convection coefficient, fluid temperature and wall temperature on the inner wall of a pipeline,”
*Progress in Nuclear Energy*, vol. 81, pp. 161–168, 2015. View at: Publisher Site | Google Scholar - W. Peng, J.-M. Zhou, Y. Li, Y. Yang, and M.-R. Guo, “A dynamic technique for the measurement of thermal conductivity of molten salt based on cylindrical melting model,”
*Journal of Thermal Analysis and Calorimetry*, vol. 115, no. 2, pp. 1767–1777, 2014. View at: Publisher Site | Google Scholar - N. N. Salva and D. A. Tarzia, “A sensitivity analysis for the determination of unknown thermal coefficients through a phase-change process with temperature-dependent thermal conductivity,”
*International Communications in Heat and Mass Transfer*, vol. 38, no. 4, pp. 418–424, 2011. View at: Publisher Site | Google Scholar - E. A. Santillan Marcus and D. A. Tarzia, “Determination of one unknown thermal coefficient of a semi-infinite porous material through a desublimation problem with coupled heat and moisture flows,”
*JP Journal of Heat and Mass Transfer*, vol. 1, no. 3, pp. 251–270, 2007. View at: Google Scholar - R. S. Vajjha and D. K. Das, “Experimental determination of thermal conductivity of three nanofluids and development of new correlations,”
*International Journal of Heat and Mass Transfer*, vol. 52, no. 21-22, pp. 4675–4682, 2009. View at: Publisher Site | Google Scholar - A.-H. Wang, X.-G. Liang, and J.-X. Ren, “Constructal enhancement of heat conduction with phase change,”
*International Journal of Thermophysics*, vol. 27, no. 1, pp. 126–138, 2006. View at: Publisher Site | Google Scholar - Y. Yang, J.-N. Yu, and Z.-C. Deng, “An inverse problem of identifying the coefficient of parabolic equation,”
*Applied Mathematical Modelling*, vol. 32, no. 10, pp. 1984–1995, 2008. View at: Publisher Site | Google Scholar | MathSciNet - Y. Yang and J. Zhou, “An experimental technique for liquid/solid thermal conductivity measurements at the melting point,”
*International Journal of Thermophysics*, vol. 27, no. 1, pp. 184–194, 2006. View at: Publisher Site | Google Scholar - G. Zhou, B. Wu, W. Ji, and S. Rho, “Time- or space-dependent coefficient recovery in parabolic partial differential equation for sensor array in the biological computing,”
*Mathematical Problems in Engineering*, vol. 2015, Article ID 573932, 9 pages, 2015. View at: Publisher Site | Google Scholar | MathSciNet

#### Copyright

Copyright © 2015 Andrea N. Ceretani and Domingo A. Tarzia. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.