Steel corrosion plays a central role in different technological fields. In this article, we consider a simple case of a corrosion phenomenon which describes a pure iron dissolution in sodium chloride. This article is devoted to prove rigorously that under rather general hypotheses on the initial data, the solution of this iron dissolution model converges to a self-similar profile as $t\rightarrow +\infty$. We will do so for an equivalent formulation as presented in the book of Avner Friedman about parabolic equations (Friedman (1964) Partial Differential Equations of Parabolic Type, Prentice-Hall, Inc., Englewood Cliffs, NJ.). In order to prove the convergence result, we apply a comparison principle together with suitable upper and lower solutions.

We present an efficient, accurate computational method for a coordinate-free model of flame front propagation of Frankel and Sivashinsky. This model allows for overturned flames fronts, in contrast to weakly nonlinear models such as the Kuramoto–Sivashinsky equation. The numerical procedure adapts the method of Hou, Lowengrub and Shelley, derived for vortex sheets, to this model. The result is a nonstiff, highly accurate solver which can handle fully nonlinear, overturned interfaces, with similar computational expense to methods for weakly nonlinear models. We apply this solver both to simulate overturned flame fronts and to compare the accuracy of Kuramoto–Sivashinsky and coordinate-free models in the appropriate limit.

The work in this paper concerns the study of different approximations for one-dimensional one-phase Stefan-like problems with a space-dependent latent heat. It is considered two different problems, which differ from each other in their boundary condition imposed at the fixed face: Dirichlet and Robin conditions. The approximate solutions are obtained by applying the heat balance integral method (HBIM), the modified HBIM and the refined integral method (RIM). Taking advantage of the exact analytical solutions, we compare and test the accuracy of the approximate solutions. The analysis is carried out using the dimensionless generalised Stefan number (Ste) and Biot number (Bi). It is also studied the case when Bi goes to infinity in the problem with a convective condition, recovering the approximate solutions when a temperature condition is imposed at the fixed face. Some numerical simulations are provided in order to assert which of the approximate integral methods turns out to be optimal. Moreover, we pose an approximate technique based on minimising the least-squares error, obtaining also approximate solutions for the classical Stefan problem.

We study the extended Stefan problem which includes constitutional supercooling for the solidification of a binary alloy in a finite spherical domain. We perform an asymptotic analysis in the limits of large Lewis number and small Stefan number which allows us to identify a number of spatio-temporal regimes signifying distinct behaviours in the solidification process, resulting in an intricate boundary layer structure. Our results generalise those present in the literature by considering all time regimes for the Stefan problem while also accounting for impurities and constitutional supercooling. These results also generalise recent work on the extended Stefan problem for finite planar domains to spherical domains, and we shall highlight key differences in the asymptotic solutions and the underlying boundary layer structure which result from this change in geometry. We compare our asymptotic solutions with both numerical simulations and real experimental data arising from the casting of molten metallurgical grade silicon through the water granulation process, with our analysis highlighting the role played by supercooling in the solidification of binary alloys appearing in such applications.

In this paper, we focus on the process of mass transfer in the rotating disk apparatus formulated as a Stefan problem with consideration given to both the hydrodynamics of the process and the specific chemical reactions occurring in the bulk. The wide range in the reaction rates of the underlying chemistry allows for a natural decoupling of the problem into a simplified set of weakly coupled convective–reaction–diffusion equations for the slowly reacting chemical species and a set of algebraic relations for the species that react rapidly. An analysis of the chemical equilibrium conditions identifies an expansion parameter and a reduced model that remains valid for arbitrarily large times. Numerical solutions of the model are compared to an asymptotic analysis revealing three distinct time scales and chemical diffusion boundary layer that lies completely inside the hydrodynamic layer. Formulated as a Stefan problem, the model generalizes the work of Levich (Levich and Spalding (1962) Physicochemical hydrodynamics, vol. 689, Prentice-Hall Englewood Cliffs, NJ) and will help better understand the natural limitations of the rotating disk reaction vessel when consideration is made for the reacting chemical species.

We critically compare the practicality and accuracy of numerical approximations of phase field models and sharp interface models of solidification. Here we focus on Stefan problems, and their quasi-static variants, with applications to crystal growth. New approaches with a high mesh quality for the parametric approximations of the resulting free boundary problems and new stable discretizations of the anisotropic phase field system are taken into account in a comparison involving benchmark problems based on exact solutions of the free boundary problem.

In this paper, a meshless regularization method of fundamental solutions is proposed for a two-dimensional, two-phase linear inverse Stefan problem. The numerical implementation and analysis are challenging since one needs to handle composite materials in higher dimensions. Furthermore, the inverse Stefan problem is ill-posed since small errors in the input data cause large errors in the desired output solution. Therefore, regularization is necessary in order to obtain a stable solution. Numerical results for several benchmark test examples are presented and discussed.

In this paper, we present a new model developed in order to analyze phenomena which arise in the solidification of binary mixtures using phase-field method, which incorporates the convection effects and the action of magnetic field. The model consists of flow, concentration, phase field and energy systems which are nonlinear evolutive and coupled systems. It represents the non-isothermal anisotropic solidification process of a binary mixture together with the motion in a melt with the applied magnetic field. To illustrate our model, numerical simulations of the influence of magnetic-field on the evolution of dendrites during the solidification of the binary mixture of Nickel-Copper (Ni-Cu) are developed. The results demonstrate that the dendritic growth under the action of magnetic-field can be simulated by using our model.

Phase change in ice-water systems in the geometry of horizontal cylindrical annulus with constant inner wall temperature and adiabatic outer wall is modeled with an enthalpy-based mixture model. Solidification and melting phenomena under different temperature conditions are analyzed through a sequence of numerical calculations. In the case of freezing of water, the importance of convection and conduction as well as the influence of cold pipe temperature on time for the complete solidification are examined. As for the case of melting of ice, the influence of the inner pipe wall temperature on the shape of the ice-water interface, the flow and temperature fields in the liquid, the heat transfer coefficients and the rate of melting are analyzed. The results of numerical calculations point to good qualitative agreement with the available experimental and other numerical results.

A thermal lattice Boltzmann method (LBM) for two-phase fluid flows in nucleate pool boiling process is proposed. In the present method, a new function for heat transfer is introduced to the isothermal LBM for two-phase immiscible fluids with large density differences. The calculated temperature is substituted into the pressure tensor, which is used for the calculation of an order parameter representing two phases so that bubbles can be formed by nucleate boiling. By using this method, two-dimensional simulations of nucleate pool boiling by a heat source on a solid wall are carried out with the boundary condition for a constant heat flux. The flow characteristics and temperature distribution in the nucleate pool boiling process are obtained. It is seen that a bubble nucleation is formed at first and then the bubble grows and leaves the wall, finally going up with deformation by the buoyant effect. In addition, the effects of the gravity and the surface wettability on the bubble diameter at departure are numerically investigated. The calculated results are in qualitative agreement with other theoretical predictions with available experimental data.