Journal of Applied Mathematics

Volume 2015, Article ID 375930, 9 pages

http://dx.doi.org/10.1155/2015/375930

## Explicit Solutions for the Solomon-Wilson-Alexiades’s Mushy Zone Model with Convective or Heat Flux Boundary Conditions

^{1}Departamento de Matemática, FCE, Universidad Austral, Paraguay 1950, S2000FZF Rosario, Argentina^{2}CONICET, Argentina

Received 14 June 2015; Revised 22 October 2015; Accepted 29 October 2015

Academic Editor: Mustafa Inc

Copyright © 2015 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.

#### Abstract

We complete the Solomon-Wilson-Alexiades’s mushy zone model (Solomon, 1982) for the one-phase Lamé-Clapeyron-Stefan problem by obtaining explicit solutions when a convective or heat flux boundary condition is imposed on the fixed face for a semi-infinite material. We also obtain the necessary and sufficient condition on data in order to get the explicit solutions for both cases which is new with respect to the original model. Moreover, when these conditions are satisfied, the two phase-change problems are equivalent to the same problem with a temperature boundary condition on the fixed face and therefore an inequality for the coefficient which characterized one of the two free interfaces of the model is also obtained.

#### 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 [1–9]. A review of a long bibliography on moving and free boundary problems for phase-change materials (PCM) for the heat equation is shown in [10]. The importance of obtaining explicit solutions to some free boundary problems was given in the works [11–26].

We consider a semi-infinite material, with constant thermal coefficients, that is initially assumed to be liquid at its melting temperature which is assumed to be equal to 0°C. At time , a heat flux or a convective boundary condition is imposed at the fixed face , and a solidification process begins where three regions can be distinguished [27, 28]:(H1)liquid region at the temperature 0°C, in , ;(H2)solid region at the temperature , in , ;(H3)mushy region at the temperature , in , . The mushy region is considered isothermal and we make the following assumptions on its structure:(H3i)the material contains a fixed portion (with ) of the total latent heat (see condition (3) in below);(H3ii)the width of the mushy region is inversely proportional to the gradient of temperature (see condition (4) below).

Following the methodology given in [27–29] and the recent one in [30], we consider a convective boundary condition in Sections 2 to 4 and a heat flux condition in Sections 5 and 6 at the fixed face , respectively. In both cases, we obtain explicit solutions for the temperature and the two free boundaries which define the mushy region. We also obtain, for both cases, the necessary and sufficient condition on data in order to get these explicit solutions given in Sections 2 and 5, respectively, which is new with respect to the original model when a temperature boundary condition at the face was imposed. Moreover, these two problems are equivalent to the same phase-change process with a temperature boundary condition on the fixed face studied in [27] and therefore an inequality for the coefficient which characterized one of the two free interfaces is also obtained in Sections 4 and 6. Moreover, in Section 3, we obtain the convergence of the solution of the phase-change process with a convective boundary condition to the solution given in [27] for a temperature boundary condition at the fixed face when the heat transfer coefficient goes to infinity, and we also give the order of the corresponding convergence when the coefficient that characterized the transient heat transfer at goes to infinity.

This paper completes the model given in [27] by considering two new boundary conditions (convective and heat flux) at the fixed face of the PCMs and obtaining explicit solutions for both cases when a restriction on data is satisfied.

#### 2. Explicit Solution with a Convective Boundary Condition

The phase-change process consists in finding the free boundaries and and the temperature such that the following conditions must be verified (problem ()):

Condition (6) represents a convective boundary condition (Robin condition) at the fixed face [31–33] with a heat transfer coefficient which is inversely proportional to the square root of the time [29, 30, 34, 35]. Now, we will obtain the solution of problem (1)–(6) when data satisfy the restriction (7).

Theorem 1. *If the coefficient satisfies the inequalitythen the solution of problem (1)–(6) is given bywithand the coefficient is given as the unique solution of the equationwhere the real functions and are defined by*

*Proof. *Taking into account that is a solution of the heat equation (3) [3], we propose as a solution of problem (1)–(6) the following expression:where the two coefficients and must be determined.

From condition (4), we deduce the expression (9) for the free boundary , where the coefficient must be determined. From conditions (6) and (2), we deduce the system of equationswhose solution is given byand then we get expression (8) for the temperature.

From condition (4), we deduce expression (10) for the interface and expression (11) for . From condition (3), we deduce (12) for the coefficient . Functions and have the following properties:Therefore, we deduce that (12) has a unique solution when the coefficient satisfies the inequalitythat is, inequality (7) holds.

Corollary 2. *If the coefficient satisfies inequality (7), then the temperature, defined by (8), verifies the following inequalities:*

*Proof. *From (8), we obtainMoreover, from (8) and (20), we also getthat is, (19) holds.

#### 3. Asymptotic Behavior When the Coefficient

Now, we will obtain the asymptotic behaviour of solution (8)–(12) of problem (1)–(6) when the heat transfer coefficient is large, that is, when . From a physical point of view, it must be convergent to the solution of the same problem with a temperature boundary condition at the fixed face given by (23).

For any coefficient satisfying inequality (7), we will denote the temperature and the two free boundaries and (defined in (8), (9), and (10), resp.) by , , and , respectively, with coefficients and . We will also denote by and the functions defined in (13). We have the following result.

Theorem 3. *One obtains the following limits:where , , and are the solutions of the following phase-change process with mushy region: (1)–(5) and instead of the boundary condition (6).*

*Proof. *The solution of problem (1)–(5) and (23) is given by [27]withand the coefficient given as the unique solution of the equationwhere the real function is defined by withThen,And, therefore, the limits (22) hold.

Now, by studying the real functions and as functions of the variable , we can obtain the order of the convergence of solution (8)–(12) of problem (1)–(6) to solution (24)–(28) of problem (1)–(5) and (23) when .

Theorem 4. *When the variable , one obtains the following estimations:*

*Proof. *As the variable , we can consider that and then solution (8)–(12) of problem (1)–(6) is well defined.

Function is an increasing function in variable ; therefore, function is a decreasing function in variable . Then, function is a decreasing function in variable , , which is convergent to as because (30) and (31) hold.

By using (13) and (31), we haveBy using (13) and (30), we haveTherefore, we haveThen, the estimation (33) holds andBy using (11) and (27), we getFinally, by using (8) and (24), we getand the thesis holds. In the particular case, when , we have

#### 4. Equivalence between the Mushy Zone Models with Convective and Temperature Boundary Conditions

We consider problem () defined by conditions (1)–(5) and temperature boundary conditionat the fixed face , whose solution was given in [27]. We have the following property.

Theorem 5. *If the coefficient satisfies inequality (7), then problem (), defined by conditions (1)–(6), is equivalent to problem (), defined by conditions (1)–(5) and (45), when the parameter in problem () is related to parameters and in problem () by the following expression:where the coefficient is given as the unique solution of (12) for problem () or as the unique solution of equationfor problem (), where the real function is defined by *

*Proof. *If the coefficient satisfies inequality (7), then the solution of problem () is given by (8)–(12). Taking into account thatwe can define problem () by imposing the temperature boundary condition (45) with data given in (46). By using this data in problem () and the method developed in [30], we can prove that the solutions of both problems () and () are the same and then the two problems are equivalent.

Corollary 6. *If the coefficient satisfies inequality (7), then the coefficient of the solid-mushy zone interface of problem () verifies the following inequality:Then,*

*Remark 7. *The real functions , defined in (30), and , defined in (48), are similar; the difference between them is the parameters or used in each definition.

#### 5. Explicit Solution with a Heat Flux Boundary Condition

Now, we will consider a phase-change process which consists in finding the free boundaries and and the temperature such that the following conditions must be verified (problem ()): conditions (1)–(5) and

Condition (52) represents the heat flux at the fixed face characterized by a coefficient which is inversely proportional to the square root of the time [34].

Theorem 8. *If the coefficient satisfies the inequalitythen the solution of problem (1)–(5) and (52) is given bywithand the coefficient given as the unique solution of the equationwhere the real function is defined by*

*Proof. *Following the proof of Theorem 1, we propose as a solution of problem (1)–(5) and (52) the following expression:where the two coefficients and must be determined.

From condition (2), we deduce expression (55) for the free boundary , with the coefficient to be determined. From conditions (2) and (52), we deduceand then we get expression (54) for the temperature.

From condition (4), we deduce expression (56) for the interface and expression (57) for . From condition (3), we deduce (58) for the coefficient . Since function has the following properties:we can deduce that (58) has a unique solution when the coefficient satisfies the inequalitywhich is inequality (53).

*Remark 9. *We have a relationship between (the lower limit for coefficient in order to have a phase-change process with a mushy region with a heat flux boundary condition at ) and (the lower limit for the coefficient in order to have a phase-change process with a mushy region with a convective boundary condition at ) given by

#### 6. Equivalence between the Mushy Zone Models with Heat Flux and Temperature Boundary Conditions

Following Section 4, we will now study the relationship between problems () and (). We have the following property.

Theorem 10. *If the coefficient satisfies inequality (53), then problem (), defined by conditions (1)–(5) and (52), is equivalent to problem (), defined by conditions (1)–(5) and (45), when the parameter in problem () is related to the parameter in problem () by the following expression:where the coefficient is given as the unique solution of (58) for problem () or as the unique solution of (47) for problem ().*

*Proof. *If the coefficient satisfies inequality (53), then the solution of problem () is given by (54)–(58). Taking into account thatwe can define problem () by imposing the temperature boundary condition (45) with the data given in (65). By using this data in problem () and the method developed in [30], we can prove that the solutions of both problems () and () are the same and then the two problems are equivalent.

Corollary 11. *If the coefficient satisfies inequality (53), then the coefficient of the solid-mushy zone interface of problem () verifies inequality (51) which is the same as that we have obtained through the equivalence between problems () and ().*

*Remark 12. *At last, for a suggestion of an anonymous referee, we will transform problem (), given by the equations and conditions (1)–(6), and inequality (7) in a dimensionless form. We define the following dimensionless change of variables:where is a characteristic length. Therefore, the equations and conditions (1)–(6) are transformed aswhere is the Stefan number and is the Biot number defined by the following expressions:Moreover, inequality (7) for the physical coefficient , which characterized the heat transfer coefficient in the boundary condition (6), is transformed in the following way: Therefore, limit in problem (1)–(6) in physical variables is equivalent to limit in problem (68) in dimensionless variables.

By using the results of this work, we can now obtain new explicit expression for the determination of one or two unknown thermal coefficients through a phase-change process with a mushy zone by imposing an overspecified convective boundary condition at the fixed face . This will complete and improve the results obtained previously in [28].

#### 7. Conclusions

The goal of this paper is to complete the solution of Solomon-Wilson-Alexiades’s model for a mushy zone model for phase-change materials when a convective or a heat flux boundary condition at the fixed face is imposed. In both cases, explicit solutions for the temperature and the two free boundaries which define the mushy region were obtained and, for both cases, the necessary and sufficient conditions on data in order to get these explicit solutions are also obtained which is new with respect to the original model when a temperature boundary condition at the face was imposed. Moreover, the equivalence of this two phase-change process with the one with a temperature boundary condition on the fixed face was obtained and an inequality for the dimensionless coefficient that characterizes the first free boundary is also given. On the other hand, the convergence of the phase-change process with mushy zone when the heat transfer coefficient goes to infinity was also obtained and the order of the convergence is also shown.

#### Nomenclature

Dimensionless coefficient that characterizes the transient heat transfer at (Biot number), defined in (70) | |

Specific heat, J/(kg°C) | |

Temperature at the fixed face , °C | |

Bulk temperature at the fixed face , °C | |

Coefficient that characterizes the transient heat transfer at , kg/(C°s^{5/2}) | |

Thermal conductivity, W/(m°C) | |

Latent heat of fusion by unit of mass, J/kg | |

Characteristic length, m | |

Phase-change process defined by conditions (1)–(6) | |

Phase-change process defined by conditions (1)–(5) and (45) | |

Phase-change process defined by conditions (1)–(5) and (52) | |

Coefficient that characterizes the transient heat flux at , kg/s^{5/2} | |

Position of the liquid-mushy zone interface at time , m | |

Dimensionless position of the liquid-mushy zone interface at time | |

Position of the solid-mushy zone interface at time , m | |

Dimensionless position of the solid-mushy zone interface at time | |

Stefan number, defined in (69) | |

Time, s | |

Temperature of the solid phase, °C | |

Spatial coordinate, m. |

*Greek Symbols*

Diffusivity coefficient, m^{2}/s | |

One of the two coefficients that characterizes the mushy zone, °C | |

One of the two coefficients that characterizes the mushy zone, being dimensionless | |

Coefficient that characterizes the free boundary in (45), being dimensionless | |

Density of mass, kg/m^{3} | |

Coefficient that characterizes the free boundary in (10), being dimensionless | |

Coefficient that characterizes the free boundary in (26), being dimensionless | |

Coefficient that characterizes the free boundary in (43), being dimensionless | |

Coefficient that characterizes the free boundary in (9), being dimensionless | |

Coefficient that characterizes the free boundary in (25), being dimensionless | |

Dimensionless time | |

Dimensionless temperature of the solid phase | |

Dimensionless spatial coordinate. |

#### Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors would like to thank an anonymous referee for his constructive comments which improves the readability of the paper. The present work has been partially sponsored by the Projects PIP no. 0534 from CONICET, Universidad Austral, Rosario, Argentina and AFOSR-SOARD Grant FA9550-14-1-0122.

#### References

- V. Alexiades and A. D. Solomon,
*Mathematical Modeling of Melting and Freezing Processes*, Hemisphere-Taylor & Francis, Washington, DC, USA, 1996. - J. R. Cannon,
*The One-Dimensional Heat Equation*, Addison-Wesley, Menlo Park, Calif, USA, 1984. - H. S. Carslaw and C. J. Jaeger,
*Conduction of Heat in Solids*, Clarendon Press, Oxford, UK, 1959. - J. Crank,
*Free and Moving Boundary Problems*, Clarendon Press, Oxford, UK, 1984. View at MathSciNet - A. Fasano, “Mathematical models of some diffusive processes with free boundary,”
*Matematica—Serie A*, vol. 11, pp. 1–128, 2005. View at Google Scholar - S. C. Gupta,
*The Classical Stefan Problem. Basic Concepts, Modelling and Analysis*, Elsevier, Amsterdam, The Netherlands, 2003. - V. J. Lunardini,
*Heat Transfer with Freezing and Thawing*, Elsevier, London, UK, 1991. - L. I. Rubinstein,
*The Stefan Problem*, American Mathematical Society, Providence, RI, USA, 1971. - A. B. Tayler,
*Mathematical Models in Applied Mechanics*, Clarendon Press, Oxford, UK, 1986. View at MathSciNet - D. A. Tarzia, “A bibliography on moving-free boundary problems for heat diffusion equation. The Stefan problem,”
*MAT Serie A*, vol. 2, pp. 1–297, 2000. View at Google Scholar - N. M. Arifin, R. Nazar, and I. Pop, “Similarity solution of Marangoni convection boundary layer flow over a flat surface in a nanofluid,”
*Journal of Applied Mathematics*, vol. 2013, Article ID 634746, 8 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - A. C. Briozzo and D. A. Tarzia, “Explicit solution of a free boundary problem for a nonlinear absorption model of mixed saturated-unsaturated flow,”
*Advances in Water Resources*, vol. 21, no. 8, pp. 713–721, 1998. View at Publisher · View at Google Scholar · View at Scopus - P. Broadbridge, “Solution of a nonlinear absorption model of mixed saturated-unsaturated flow,”
*Water Resources Research*, vol. 26, no. 10, pp. 2435–2443, 1990. View at Publisher · View at Google Scholar · View at Scopus - F. Font, S. L. Mitchell, and T. G. Myers, “One-dimensional solidification of supercooled melts,”
*International Journal of Heat and Mass Transfer*, vol. 62, no. 1, pp. 411–421, 2013. View at Publisher · View at Google Scholar · View at Scopus - Y. Huang and Y. Shang, “The Bäcklund transformations and abundant exact explicit solutions for a general nonintegrable nonlinear convection-diffusion equation,”
*Abstract and Applied Analysis*, vol. 2012, Article ID 489043, 11 pages, 2012. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - A. Hussanan, I. Khan, and S. Shafie, “An exact analysis of heat and mass transfer past a vertical plate with Newtonian heating,”
*Journal of Applied Mathematics*, vol. 2013, Article ID 434571, 9 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - D. Kaya, “Explicit solutions of generalized nonlinear Boussinesq equations,”
*Journal of Applied Mathematics*, vol. 1, no. 1, pp. 29–37, 2001. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - T. G. Myers and F. Font, “On the one-phase reduction of the Stefan problem with a variable phase change temperature,”
*International Communications in Heat and Mass Transfer*, vol. 61, pp. 37–41, 2015. View at Publisher · View at Google Scholar · View at Scopus - T. G. Myers, S. L. Mitchell, and F. Font, “Energy conservation in the one-phase supercooled Stefan problem,”
*International Communications in Heat and Mass Transfer*, vol. 39, no. 10, pp. 1522–1525, 2012. View at Publisher · View at Google Scholar · View at Scopus - M. F. Natale and D. A. Tarzia, “Explicit solutions to the two-phase Stefan problem for storm-type materials,”
*Journal of Physics A: Mathematical and General*, vol. 33, no. 2, pp. 395–404, 2000. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - C. Rogers, “On a class of reciprocal Stefan moving boundary problems,”
*Zeitschrift für Angewandte Mathematik und Physik*, vol. 66, no. 4, pp. 2069–2079, 2015. View at Publisher · View at Google Scholar · View at MathSciNet - N. N. Salva and D. A. Tarzia, “Explicit solution for a Stefan problem with variable latent heat and constant heat flux boundary conditions,”
*Journal of Mathematical Analysis and Applications*, vol. 379, no. 1, pp. 240–244, 2011. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - A. Sojoudi, A. Mazloomi, S. C. Saha, and Y. T. Gu, “Similarity solutions for flow and heat transfer of non-Newtonian fluid over a stretching surface,”
*Journal of Applied Mathematics*, vol. 2014, Article ID 718319, 8 pages, 2014. View at Publisher · View at Google Scholar · View at Scopus - A. D. Solomon, D. G. Wilson, and V. Alexiades, “Explicit solutions to change problems,”
*Quarterly of Applied Mathematics*, vol. 41, pp. 237–243, 1983. View at Google Scholar - V. R. Voller and F. Falcini, “Two exact solutions of a Stefan problem with varying diffusivity,”
*International Journal of Heat and Mass Transfer*, vol. 58, no. 1-2, pp. 80–85, 2013. View at Publisher · View at Google Scholar · View at Scopus - V. R. Voller, J. B. Swenson, and C. Paola, “An analytical solution for a Stefan problem with variable latent heat,”
*International Journal of Heat and Mass Transfer*, vol. 47, no. 24, pp. 5387–5390, 2004. View at Publisher · View at Google Scholar · View at Scopus - 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 · View at Google Scholar · View at Scopus - D. A. Tarzia, “Determination of unknown thermal coefficients of a semi-infinite material for the one-phase Lame-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 · View at Google Scholar - D. A. Tarzia, “An explicit solution for a two-phase unidimensional Stefan problem with a convective boundary condition at the fixed face,”
*MAT Serie A*, vol. 8, pp. 21–27, 2004. View at Google Scholar - D. A. Tarzia, “Relationship between Neumann solutions for two-phase Lamé-Clapeyron-Stefan problems with convective and temperature boundary conditions,”
*Thermal Science*, p. 3, 2015. View at Publisher · View at Google Scholar - P. M. Beckett, “A note on surface heat transfer coefficients,”
*International Journal of Heat and Mass Transfer*, vol. 34, no. 8, pp. 2165–2166, 1991. View at Publisher · View at Google Scholar · View at Scopus - S. D. Foss, “An approximate solution to the moving boundary problem associated with the freezing and melting of lake ice,”
*AIChE Symposium Series*, vol. 74, no. 174, pp. 250–255, 1978. View at Google Scholar · View at Scopus - H. Ching-Lun and S. Yen-Ping, “Perturbation solution for planar solidification of a saturated liquid with convection at the wall,”
*International Journal of Heat and Mass Transfer*, vol. 18, no. 12, pp. 1481–1483, 1975. View at Publisher · View at Google Scholar · View at Scopus - D. A. Tarzia, “An inequality for the coefficient
*σ*of the free boundary $s\left(t\right)=2\sigma \sqrt{t}$ of the Neumann solution for the two-phase Stefan problem,”*Quarterly of Applied Mathematics*, vol. 39, no. 4, pp. 491–497, 1981. View at Google Scholar · View at MathSciNet - S. M. Zubair and M. A. Chaudhry, “Exact solutions of solid-liquid phase-change heat transfer when subjected to convective boundary conditions,”
*Heat and Mass Transfer*, vol. 30, no. 2, pp. 77–81, 1994. View at Publisher · View at Google Scholar