scholarly journals Verification of a mathematical model for the solution of the Stefan problem using the mushy layer method

2021 ◽  
Vol 2021 (3) ◽  
pp. 119-125
Author(s):  
R.S. Yurkov ◽  
◽  
L.I. Knysh ◽  

The use of solar energy has limitations due to its periodic availability: solar plants do not operate at night and are ineffective in dull weather. The solution of this problem involves the introduction of energy storage and duplication systems into the conversion loop. Among the energy storage systems, solid–liquid phase transition modules have significant energy, ecologic, and cost advantages. Physical processes in modules of this type are described by a system of non-stationary nonlinear partial differential equations with specific boundary conditions at the phase interface. The verification of a method for solving the Stefan problem for a heat-storage material is presented in this paper. The use of the mushy layer method made it possible to simplify the classical mathematical model of the Stefan problem by reducing it to a nonstationary heat conduction problem with an implicit heat source that takes into account the latent heat of transition. The phase transition is considered to occur in an intermediate zone determined by the solidus and liquidus temperatures rather than in in infinite region. To develop a Python code, use was made of an implicit computational scheme in which the solidus and liquidus temperatures remain constant and are determined in the course of numerical experiments. The physical model chosen for computer simulation and algorithm verification is the process of ice layer formation on a water surface at a constant ambient temperature. The numerical results obtained allow one to determine the temperature fields in the solid and the liquid phase and the position of the phase interface and calculate its advance speed. The algorithm developed was verified by analyzing the classical analytical solution of the Stefan problem for the one-dimensional case at a constant advance speed of the phase interface. The value of the verification coefficient was determined from a numerical solution of a nonlinear equation with the use of special built-in Python functions. Substituting the data for the physical model under consideration into the analytical solution and comparing them with the numerical simulation data obtained with the use of the mushy layer method shows that the results are in close agreement, thus demonstrating the correctness of the computer algorithm developed. These studies will allow one to adapt the Python code developed on the basis of the mushy layer method to the calculation of heat storage systems with a solid-liquid phase transition with account for the features of their geometry, the temperature level, and actual boundary conditions.

Author(s):  
Tomáš Roubíček

The classical Stefan problem, concerning mere heat-transfer during solid-liquid phase transition, is here enhanced towards mechanical effects. The Eulerian description at large displacements is used with convective and Zaremba-Jaumann corotational time derivatives, linearized by exploiting the additive Green-Naghdi’s decomposition in (objective) rates. In particular, the liquid phase is a viscoelastic fluid while creep and rupture of the solid phase is considered in the Jeffreys viscoelastic rheology exploiting the phase-field model, exploiting a concept of slightly (so-called “semi”) compressible materials. The $L^1$-theory for the heat equation is adopted for the Stefan problem relaxed by allowing for kinetic superheating/supercooling effects during the solid-liquid phase transition. A rigorous proof of existence of week solutions is provided for an incomplete melting, exploiting a time-discretisation approximation.


Sign in / Sign up

Export Citation Format

Share Document