Next Article in Journal
How to Find the Equivalence Classes in a Set of Linear Codes in Practice?
Previous Article in Journal
On the Analytical Solution of the SIRV-Model for the Temporal Evolution of Epidemics for General Time-Dependent Recovery, Infection and Vaccination Rates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Solid–Liquid Phase Interface Dynamics in an Undercooled Melt with a Solid Wall

by
Ekaterina A. Titova
and
Dmitri V. Alexandrov
*
Laboratory of Multi-Scale Mathematical Modeling, Department of Theoretical and Mathematical Physics, Ural Federal University, Lenin Ave., 51, Ekaterinburg 620000, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(2), 327; https://doi.org/10.3390/math12020327
Submission received: 19 December 2023 / Revised: 17 January 2024 / Accepted: 18 January 2024 / Published: 19 January 2024

Abstract

:
A new boundary integral equation for the interface function of a curved solid/liquid phase interface propagating into an undercooled one-component melt is derived in the presence of a solid wall in liquid. Green’s function technique is used to transform a purely thermal boundary value problem to a single integro-differential equation for the interface function in two- and three-dimensional cases. It is shown that a solid wall represents an additional source of heat and melt undercooling can be negative in the vicinity of the wall. The new boundary integral equation has a limiting transition to previously developed theory in the absence of a solid wall.

1. Introduction

Directional and volumetric crystallization are the basis of many technological solidification processes and are often found in nature (e.g., freezing of water and various solutions, solidification of magma) [1,2,3,4,5,6,7]. Mathematical models of such processes are based on the description of heat and mass transfer in solid and liquid phases separated by a moving crystallization front. This interphase boundary is generally curvilinear and moves according to a certain time-dependent law, which is determined by solving the problem. The mathematical model of a purely thermal problem was first formulated by Josef Stephan [8,9,10,11], and all problems with a moving boundary of phase transformation now bear his name. Note that, in general, the unsteady Stefan problem has no exact analytical solution, and the solution of each individual model strongly depends on the geometry of the crystallization domain and boundary conditions. The solution to such problems is usually constructed using approximate analytical (e.g., the method of differential series [12,13]) and numerical (e.g., the method of fixed boundaries [14]) methods.
One of the fruitful approaches to solving the Stefan problem is the method of boundary integral equation (BIE), first proposed by Nash and Glicksman [15,16]. The idea of this method is to derive an integro-differential equation for the interface function determining the position and velocity of the crystallization front. Note that such an equation was derived in [15,16] for a purely thermal problem. The BIE that describes the thermal concentration problem was derived by Alexandrov and Galenko [17] in the case of parabolic mass transfer in the melt/solution. Further, they considered the scenario of local non-equilibrium (fast) crystallization described by a hyperbolic impurity diffusion equation and derived the corresponding BIE [18]. Note that the high-speed BIE has a limiting transition to the low-speed BIE [18]. An important step in the BIE theory was the case of convective fluid flows considered in [19,20]. Note that the BIE allows us to study the morphological stability of interfacial boundaries of a certain shape [21,22,23], derive the selection criterion of stable growth mode for dendritic crystals [24,25,26], as well as study the shape and evolution of various patterns [27].
In this paper, we derive a new BIE for the case when a single-component melt or solution crystallizes under the significant influence of a solid wall in liquid. Examples of such processes are (i) dendritic growth in a mold, (ii) freezing of water in lakes of shallow depth, (iii) and solidification of lava in magma chambers. The new BIE has a limiting transition to the previously developed theory in the absence of a solid wall and is found below for 2D and 3D cases.

2. 2D Boundary Integral Equation

Let us consider the evolution of a curved solid/liquid phase transition boundary in an undercooled one-component liquid in the presence of a solid wall (Figure 1). Below, we assume that the one-component melt is motionless, and the BIE is written out in a reference frame that is moving with a constant (V) toward the liquid phase. Also, we assume that the temperature conductivity coefficient D T is the same constant in both phases. In addition, the solid wall is smooth and planar; no particle nucleation occurs in the bulk melt or at the solid wall. The temperature field in the solid and liquid phases satisfies the temperature conductivity equation in a reference frame that is moving with a constant velocity (V) [28]
D T 2 T T t + V T z = 0 ,
where T is the temperature, D T is the thermal conductivity, t is the time and V is the constant velocity of moving Cartesian coordinate systems x, y, and z fixed on the growing crystal (steady-state growth rate). The boundary conditions at the solid/liquid interface take the form [17,20,22,25]
T = T i = T 0 d c K Q c p β ˜ V + ζ t , D T T s o l i d T l i q u i d · d s = Q c p V + ζ t d 2 x ,
where d c is the anisotropy capillary length, K is the average interface curvature, Q is the latent heat of crystallization, c p is the heat capacity, β ˜ is the kinetic coefficient, T i is the interfacial temperature, T 0 is the phase transition temperature for a flat crystallization front, d s is the surface area vector element directed towards the liquid phase, and ζ ( x , t ) is the interface function. Let us especially emphasize that ζ ( x , t ) represents a curved line in 2D whereas ζ ( x , y , t ) is a crystal surface in 3D. Note that the boundary conditions (2) represent the Gibbs–Thomson and heat balance conditions, respectively. The average interface curvature K in a two-dimensional case is determined as
K ( x , t ) = 2 ζ / x 2 1 + ζ / x 2 3 / 2 .
Here, we consider how a solid wall z = z 0 kept at a fixed temperature
T ( z 0 ) = T 1 = const .
influences the boundary integral equation. To do this, we use Green’s function G ( p | p 1 ) to derive an integro-differential equation for the interface function ζ x , t . Green’s function satisfies the following equation and boundary condition [28,29]
G ( p | p 1 ) t 1 + D T 1 2 G ( p | p 1 ) V G ( p | p 1 ) z 1 = δ ( p p 1 ) , G ( x , z 0 , t | x 1 , z 1 , t 1 ) = 0 ,
where p = ( x , z , t ) and p 1 = ( x 1 , z 1 , t 1 ) .
Applying the method of images to Green’s function for an infinite domain [29], we construct Green’s function for the problem under question with a solid wall as follows
G ( p | p 1 ) = 1 4 π D T ( t t 1 ) exp ( x x 1 ) 2 + ( z z 1 + V ( t t 1 ) ) 2 4 D T ( t t 1 ) exp ( x x 1 ) 2 + ( z 2 z 0 + z 1 V ( t t 1 ) ) 2 4 D T ( t t 1 ) .
It is easily seen that Green’s function (6) satisfies the boundary value problem (5). Note that the first term in expression (6) coincides with Green’s function for the unbounded problem [30] (when a solid wall is absent).
Next, multiplying Equation (1) taken at point p 1 by G ( p | p 1 ) and then subtracting Equation (5) multiplied by T ( p 1 ) , we obtain
D T G ( p | p 1 ) 1 2 T T ( p 1 ) 1 2 G ( p | p 1 ) G ( p | p 1 ) T t 1 T ( p 1 ) G ( p | p 1 ) t 1 + V G ( p | p 1 ) T z 1 + T ( p 1 ) G z 1 = δ ( p p 1 ) T ( p 1 ) .
Now, we integrate Equation (7) over time and volume
D T t + ε Λ 1 G ( p | p 1 ) 1 2 T T ( p 1 ) 1 2 G ( p | p 1 ) d t 1 d x 1 d z 1 Λ 1 G ( p | p 1 ) T ( p 1 ) t 1 t 1 t + ε d x 1 d z 1 + V t + ε G ( p | p 1 ) T ( p 1 ) z 1 z 1 z 0 d x 1 d t 1 = T ( p ) ,
where ε > 0 is an infinitely small parameter and Λ 1 is an integration domain without the interface ζ ( x , t ) ( < x 1 < , < z 1 < z 0 ).
The second integral in Equation (8) vanishes since Equation (6) gives zero in the limit t 1 , and G ( p | p 1 ) = 0 when t 1 > t (causality condition). Considering the third integral in (8) with allowance for T ( z 1 z 0 ) = T 1 = const . and T ( z 1 ) = const . , we are able to integrate Green’s function as follows
V t + ε G ( p | p 1 ) d x 1 d t 1 = exp ( z 1 z ) V D T , z > z 1 , 1 , z < z 1 . exp ( z 1 + z 2 z 0 ) V D T , z 2 z 0 + z 1 < 0 , 1 , z 2 z 0 + z 1 > 0 .
Here, the coordinate z 1 ( z 0 , ) . At z 1 , both terms on the right-hand side of expression (9) are vanishing. Considering the boundary z 1 = z 0 and taking into account that the total domain of the problem under question is bounded by the solid wall, we arrive at
V t + ε G ( p | p 1 ) T ( p 1 ) z 1 z 1 z 0 d x 1 d t 1 = T 1 1 exp ( z z 0 ) V D T .
Now, we apply Green’s second identity to the first term in (8)
t + ε Λ 1 G ( p | p 1 ) 1 2 T T ( p 1 ) 1 2 G ( p | p 1 ) d t 1 d x 1 d z 1 = t + ε S 1 G ( p | p 1 ) 1 T T ( p 1 ) 1 G ( p | p 1 ) · d S 1 d t 1 ,
where S 1 is a boundary encompassing the volume Λ 1 in 3D. By substituting (10) and (11) into (8), we have
D T t + ε S 1 G ( p | p 1 ) 1 T T ( p 1 ) 1 G ( p | p 1 ) · d S 1 d t 1 + T 1 1 exp ( z z 0 ) V D T = T ( p ) .
The integral containing 1 G ( p | p 1 ) vanishes due to the continuity of temperature T [28,30]. Approaching the point p to the interface and using the heat balance condition (2), we come to
Q c p t + ε G ( p | p 1 ) V + ζ 1 t 1 d t 1 d x 1 + T 1 1 exp ( ζ z 0 ) V D T = T i ,
where ζ 1 = ζ ( x 1 , t 1 ) .
The final step is to use the Gibbs–Thomson condition (2), which gives
Δ d c K + T 2 exp ( ζ z 0 ) V D T β V + ζ t = t + ε G ( p | p 1 ) V + ζ 1 t 1 d t 1 d x 1 .
Here, Green’s function G ( p | p 1 ) is defined by expression (6), β = β ˜ c p / Q and we introduce the dimensionless undercooling Δ = ( T 0 T 1 ) c p / Q and temperature T 2 = T 1 c p / Q . Note that Equation (14) transforms to the corresponding BIE in the case of an unbounded domain at z 0 , i.e., (see, among others, [15,30])
Δ d c K β V + ζ t = 1 4 π D T t + ε 1 t t 1 exp ( x x 1 ) 2 + ( ζ ζ 1 + V ( t t 1 ) ) 2 4 D T ( t t 1 ) V + ζ 1 t 1 d t 1 d x 1 .
As a special note, Equation (15) has a particular solution for a parabolic/paraboloidal dendrite previously discussed by Brener and Mel’nikov [25]. These dendritic shapes are the exact solutions of the BIE (15) in 2D and 3D cases when the interface curvature K is negligible or constant.

3. 2D Stationary Crystallization

Let us write the BIE (14) in dimensionless coordinates, where the dendrite tip diameter ρ plays the role of the length scale, and z 0 , ζ , and x ( ζ 1 and x 1 ) are now dimensionless (for the sake of simplicity, we use here the same designations). In addition, we use the dimensionless time and spatial coordinates in the BIE with time scale ρ / V . As a result, we have
Δ d c ρ K + T 2 exp ( ζ z 0 ) 2 P T β V 1 + ζ t = P T 2 π 0 exp P T 2 τ ( x x 1 ) 2 exp P T 2 τ ( ζ ζ 1 + τ ) 2 exp P T 2 τ ( ζ 2 z 0 + ζ 1 τ ) 2 1 + ζ 1 ( x 1 , t τ ) t d τ τ d x 1 .
Here, τ = t t 1 and P T = ρ V / ( 2 D T ) is the Péclet number. Let us consider the frequently occurring case of steady-state crystal growth with a constant velocity. Under such conditions, the interface function ζ = ζ ( x ) does not depend on time and the BIE (16) transforms to
Δ d c ρ K + T 2 exp ( ζ z 0 ) 2 P T β V = P T 2 π × d x 1 exp P T ( ζ 1 ζ ) 0 d τ τ exp P T 2 ( x x 1 ) 2 + ( ζ ζ 1 ) 2 τ + τ d x 1 exp P T ( ζ 1 + ζ 0 ) 0 d τ τ exp P T 2 ( x x 1 ) 2 + ( ζ 0 + ζ 1 ) 2 τ + τ
with ζ 0 = ζ 2 z 0 . Replacing the variable τ in (17) by y 1 and y 2
τ ( x x 1 ) 2 + ( ζ ζ 1 ) 2 = y 1 , τ ( x x 1 ) 2 + ( ζ 0 + ζ 1 ) 2 = y 2 ,
and taking into account that the integrals over y 1 and y 2 are even functions, we obtain
Δ d c ρ K + T 2 exp ( ζ z 0 ) 2 P T β V = P T π × d x 1 exp P T ( ζ 1 ζ ) 0 d y 1 exp P T 2 ( x x 1 ) 2 + ( ζ ζ 1 ) 2 e y 1 + e y 1 d x 1 exp P T ( ζ 1 + ζ 0 ) 0 d y 2 exp P T 2 ( x x 1 ) 2 + ( ζ 0 + ζ 1 ) 2 e y 2 + e y 2 .
By using the definition of the modified Bessel function, we now come to the BIE describing the steady-state crystal growth
Δ d c ρ K + T 2 exp ( ζ z 0 ) 2 P T β V = P T π × d x 1 exp P T ( ζ 1 ζ ) K 0 P T ( x x 1 ) 2 + ( ζ ζ 1 ) 2 d x 1 exp P T ( ζ 1 + ζ 0 ) K 0 P T ( x x 1 ) 2 + ( ζ 0 + ζ 1 ) 2 .

4. Parabolic Dendrite

As a special case, we consider the growth of a parabolic dendrite with the dimensionless interface function
ζ = x 2 2 , ζ 1 = x 1 2 2 ,
Substituting the interface function (21) into the BIE (14) and integrating Green’s function over the unbounded region, we have
Δ d c K ρ + T 2 exp P T ( x 2 + 2 z 0 ) β V = π P T exp ( P T ) erfc ( P T ) 0 P T 2 π τ exp P T 2 τ ( x x 1 ) 2 + 2 z 0 + x 2 + x 1 2 2 + τ 2 d τ d x 1 .
Figure 2 and Figure 3 show the melt undercooling determined by Equation (22). As can be seen, the wall represents an additional source of heat and the undercooling can be negative at small distances from the wall. In this case, we have melting instead of crystallization. It should be emphasized that melt undercooling depends only on the Péclet number if there is no solid wall and if interface curvature is constant or negligible. Contrary to this, melt undercooling is a function of both the Péclet number and the distance between the crystal and the solid wall in a semi-bounded problem. In addition, this distance depends on the x-coordinate. It means that the interface temperature and melt undercooling are also dependent on x.

5. 3D Boundary Integral Equation

Here, we consider a more general and realistic case of the three-dimensional problem defined by Equations (1), (2), and (4). Green’s function for the 3D problem reads as
G 3 D ( p | p 1 ) = D T π V 2 ( t t 1 ) 3 / 2 exp | x x 1 | 2 + ( z z 1 + V ( t t 1 ) ) 2 4 D T ( t t 1 ) exp | x x 1 | 2 + ( z 2 z 0 + z 1 V ( t t 1 ) ) 2 4 D T ( t t 1 ) .
Here, p = ( x , y , z , t ) , p 1 = ( x 1 , y 1 , z 1 , t 1 ) , x = ( x , y ) , and x 1 = ( x 1 , y 1 ) . The 3D boundary integral equation can be obtained by analogy with the 2D case. In dimensionless coordinates, it takes the form of
Δ d c ρ K + T 2 exp ( ζ 2 z 0 ) P T β V 1 + ζ t = P T 2 π 3 / 2 0 d τ τ 3 / 2 exp P T 2 τ | x x 1 | 2 + ( ζ ζ 1 + τ ) 2 exp P T 2 τ | x x 1 | 2 + ( ζ 2 z 0 + ζ 1 τ ) 2 1 + ζ 1 t 1 d x 1 d y 1 .

6. Paraboloid of Revolution

The steady-state growth of a needle 3D dendrite can be modeled by the BIE with a paraboloid of revolution as an interface function
ζ = x 2 + y 2 2 , ζ 1 = x 1 2 + y 1 2 2 .
In the dimensionless time and spatial coordinates defined in Section 4, the BIE (14) for the phase transition interface (25) transforms to
Δ d c K ρ + T 2 exp P T ( x 2 + y 2 + 2 z 0 ) β V = P T exp ( P T ) 1 exp ( P T η ) η d η P T 2 π 3 / 2 0 d τ τ 3 / 2 × d x 1 d y 1 exp P T 2 τ | x x 1 | 2 + 2 z 0 + x 2 + y 2 + x 1 2 + y 1 2 2 + τ 2 .

7. Conclusions

In this paper, we derive a new boundary integral equation (BIE) for the interface function describing the motion of an interfacial boundary in a semi-bounded domain of supercooled melt adjacent to a solid wall of constant temperature. Such a temperature regime can be achieved in practice by artificially maintaining a constant temperature at the boundary (e.g., on the walls of an ingot mold, an ice surface, or the bottom of a lava chamber). This new BIE is derived using Green’s function technique for both 2D and 3D cases (formulas (14) and (24), respectively). The derived BIE transitions to the previously known theory if we take the distance between the solid/liquid interface and the solid wall to infinity. The new BIE has been analyzed for the growth of a parabolic dendrite. The solution shows that melt undercooling near a solid wall can become negative because such a wall reflects heat. This in turn leads to the melting of the crystal instead of solidification. An important fact is that the total melt undercooling depends on both the Péclet number and the distance between a curved crystal surface and the solid wall. This means that this undercooling is different at different points on the crystal surface. This means that the dendrite tip slows down as it approaches the wall due to a smaller driving force (melt supercooling). Therefore, the dendrite shape (the surface that envelops its tip and secondary branches) becomes wider near a solid wall. If we consider the growth of several dendrites together (the growth of a dendritic forest), a crystal evolving close to a given dendrite can act as a solid wall (heat source). Therefore, the growth of dendritic branches in the corresponding spatial direction will be retarded. This conclusion is confirmed by experimental data and computer simulations of dendritic growth [31], where the near-complete disappearance of secondary branches of neighboring dendrites was recorded.
In summary, the following new results have been obtained:
  • A new BIE for a curved solid/liquid interface propagating into a single-component supercooled melt with a solid wall of constant temperature has been derived.
  • An analytical solution of this BIE for a parabolic dendrite is obtained. This solution demonstrates that melt undercooling near a solid surface can be negative when the wall reflects heat and melting occurs.
  • It has been shown that the solid wall leads to a smaller driving force, which decelerates the dendrite tip motion and enlarges the crystal shape.
It is of interest to generalize this theory to the effect of a curved solid wall or several solid walls in the melt, modeling the real geometry of an ingot mold [32,33]. An important step is also the generalization of the theory for binary melts, where the impurity redistribution upstream of the solid/liquid interface as well as impurity absorption by the solid phase play a decisive role in the crystallization process and the properties of the solidified phase. Since improved properties of materials are often observed at high crystal growth rates, generalization is also required for local nonequilibrium solidification processes with hyperbolic mass transfer in liquid [18,34,35,36,37,38,39]. The BIE should also be generalized for convective solidification in the presence of a solid wall [19,20,40,41,42,43]. The BIE describing the morphology of solid phase protrusions in the two-phase region can also be used to determine the solid phase fraction in this region and, consequently, how the heat and mass transfer coefficient depends on the solid phase fraction. Therefore, it is reasonable to combine the present approach with the two-phase region theory [44,45,46,47,48,49,50] for more accurate modeling of the directional solidification process. The influence of these effects is important for studying the dynamics of a curved crystallization front in the presence of a solid wall. For this purpose, it is necessary to use the boundary integral method under consideration (developed in the presence of a solid wall) and the theories of the aforementioned papers (elaborated on in the absence of a solid wall).

Author Contributions

Conceptualization, E.A.T. and D.V.A.; methodology, E.A.T. and D.V.A.; software, E.A.T.; validation, E.A.T. and D.V.A.; formal analysis, E.A.T. and D.V.A.; investigation, E.A.T. and D.V.A.; resources, D.V.A.; writing—original draft preparation, E.A.T. and D.V.A.; writing—review and editing, E.A.T. and D.V.A.; visualization, E.A.T.; supervision, E.A.T. and D.V.A.; project administration, E.A.T. and D.V.A.; funding acquisition, E.A.T. and D.V.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Ministry of Science and Higher Education of the Russian Federation (project no. FEUZ-2023-0022).

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kurz, W.; Fisher, D.J.; Trivedi, R. Progress in modelling solidification microstructures in metals and alloys: Dendrites and cells from 1700 to 2000. Int. Mater. Rev. 2019, 64, 311–354. [Google Scholar] [CrossRef]
  2. Kurz, W.; Rappaz, M.; Trivedi, R. Progress in modelling solidification microstructures in metals and alloys. Part II: Dendrites from 2001 to 2018. Int. Mater. Rev. 2021, 66, 30–76. [Google Scholar] [CrossRef]
  3. Worster, M.G. Solidification of an alloy from a cooled boundary. J. Fluid Mech. 1986, 167, 481–501. [Google Scholar] [CrossRef]
  4. Wettlaufer, J.S.; Worster, M.G.; Huppert, E. The phase evolution of young sea ice. Geophys. Res. 1997, 24, 1251–1254. [Google Scholar] [CrossRef]
  5. Worster, M.G.; Huppert, H.E.; Sparks, R.S.J. The crystallization of lava lakes. J. Geophys. Res. 1993, 98, 15891–15901. [Google Scholar] [CrossRef]
  6. Makoveeva, E.V. Steady-state crystallization with a mushy layer: A test of theory with experiments. Eur. Phys. J. Spec. Top. 2023, 232, 1165–1169. [Google Scholar] [CrossRef]
  7. Herlach, D.; Galenko, P.; Holland-Moritz, D. Metastable Solids from Undercooled Melts; Elsevier: Amsterdam, The Netherlands, 2007. [Google Scholar]
  8. Wettlaufer, J.S. The Stefan problem: Polar exploration and the mathematics of moving boundaries. In Die Zentralanstalt für Meteorologie und Geodynamik, 1851–2001, 150 Jahre Meteorologie und Geophysik in Osterreich; Styria Verlag: Graz, Austria, 2001; pp. 420–435. [Google Scholar]
  9. Rubinstein, L.I. The Stefan Problem; American Mathematical Society: Providence, RI, USA, 1971. [Google Scholar]
  10. Meirmanov, A.M. The Stefan Problem; De Gruyter Expositions in Mathematics; De Gruyter: Berlin, Germany, 1992. [Google Scholar]
  11. Alexandrov, D.V.; Ivanov, A.A. The Stefan problem of solidification of ternary systems in the presence of moving phase transition regions. J. Exper. Theor. Phys. 2009, 108, 821–829. [Google Scholar] [CrossRef]
  12. Lubov, B.Y. The Theory of Crystallization in Large Volumes; Nauka: Moscow, Russia, 1975. [Google Scholar]
  13. Alexandrov, D.V. Nucleation and evolution of spherical crystals with allowance for their unsteady-state growth rates. J. Phys. A Math. Theor. 2018, 51, 075102. [Google Scholar] [CrossRef]
  14. Samarskii, A.A.; Vabishchevich, P.N. Computational Heat Transfer, Vol.1, Mathematical Modelling; Wiley: Chichester, UK, 1995. [Google Scholar]
  15. Nash, G.E. Capillary-limited, steady state dendritic growth. Part I. Theoretical development. NRL Rep. 1974, 7679. Available online: https://apps.dtic.mil/sti/citations/AD0780781 (accessed on 1 December 2023).
  16. Nash, G.E.; Glicksman, M.E. Capillary-limited steady-state dendritic growth—I. Theoretical development. Acta Metall. 1974, 22, 1283–1290. [Google Scholar] [CrossRef]
  17. Alexandrov, D.V.; Galenko, P.K. Boundary integral approach for propagating interfaces in a binary non-isothermal mixture. Physica A 2017, 469, 420–428. [Google Scholar] [CrossRef]
  18. Alexandrov, D.V.; Galenko, P.K. Selected mode for rapidly growing needle-like dendrite controlled by heat and mass transport. Acta Mater. 2017, 137, 64–70. [Google Scholar] [CrossRef]
  19. Saville, D.A.; Beaghton, P.J. Growth of needle-shaped crystals in the presence of convection. Phys. Rev. A 1988, 37, 3423–3430. [Google Scholar] [CrossRef] [PubMed]
  20. Titova, E.A.; Alexandrov, D.V. The boundary integral equation for curved solid/liquid interfaces propagating into a binary liquid with convection. J. Phys. A Math. Theor. 2022, 55, 055701. [Google Scholar] [CrossRef]
  21. Ben Amar, M.; Pelcé, P. Impurity effect on dendritic growth. Phys. Rev. A 1989, 39, 4263–4269. [Google Scholar] [CrossRef] [PubMed]
  22. Bouissou, P.; Pelcé, P. Effect of a forced flow on dendritic growth. Phys. Rev. A 1989, 40, 6673–6680. [Google Scholar] [CrossRef] [PubMed]
  23. Alexandrov, D.V.; Galenko, P.K.; Herlach, D.M. Selection criterion for the growing dendritic tip in a non-isothermal binary system under forced convective flow. J. Cryst. Growth 2010, 312, 2122–2127. [Google Scholar] [CrossRef]
  24. Barbieri, A.; Langer, J.S. Predictions of dendritic growth rates in the linearized solvability theory. Phys. Rev. A 1989, 39, 5314–5325. [Google Scholar] [CrossRef]
  25. Brener, E.A.; Mel’nikov, V.I. Pattern selection in two-dimensional dendritic growth. Adv. Phys. 1991, 40, 53–97. [Google Scholar] [CrossRef]
  26. Brener, E.A.; Mel’nikov, V.A. Two-dimensional dendritic growth at arbitrary Peclet number. J. Phys. Fr. 1990, 51, 157–166. [Google Scholar] [CrossRef]
  27. Pelce, P.; Pomeau, Y. Dendrites in the small undercooling limit. Stud. Appl. Math. 1986, 74, 245–258. [Google Scholar] [CrossRef]
  28. Langer, J.S. Studies in the theory of interfacial stability—II. Moving symmetric model. Acta Metall. 1977, 25, 1121–1137. [Google Scholar] [CrossRef]
  29. Morse, P.M.; Feshbach, H. Methods of Theoretical Physics; McGraw-Hill: New York, NY, USA, 1953. [Google Scholar]
  30. Langer, J.S.; Turski, L.A. Studies in the theory of interfacial stability—I. Stationary symmetric model. Acta Metall. 1977, 25, 1113–1119. [Google Scholar] [CrossRef]
  31. Hu, M.; Sun, C.; Fang, H.; Zhu, M. Competitive dendrite growth during directional solidification of a transparent alloy: Modeling and experiment. Eur. Phys. J. E 2020, 43, 16. [Google Scholar] [CrossRef] [PubMed]
  32. Kermanpur, A.; Eskandari, M.; Purmohamad, H.; Soltani, M.A.; Shateri, R. Influence of mould design on the solidification of heavy forging ingots of low alloy steels by numerical simulation. Mater. Des. 2010, 31, 1096–1104. [Google Scholar] [CrossRef]
  33. Abootorabi, A.; Korojy, B.; Jabbareh, M.A. Effect of mould design on the Niyama criteria during solidification of CH3C 80t ingot. Ironmak. Steelmak. 2019, 47, 722–730. [Google Scholar] [CrossRef]
  34. Galenko, P.K.; Danilov, D.A. Local nonequilibrium effect on rapid dendritic growth in a binary alloy melt. Phys. Lett. A 1997, 235, 271–280. [Google Scholar] [CrossRef]
  35. Galenko, P.K.; Danilov, D.A. Model for free dendritic alloy growth under interfacial and bulk phase nonequilibrium conditions. J. Cryst. Growth 1999, 197, 992–1002. [Google Scholar] [CrossRef]
  36. Galenko, P.K.; Danilov, D.A. Selection of the thermodynamically stable regime of rapid solidification front motion in an isothermal binary alloy. J. Cryst. Growth 2000, 216, 512–526. [Google Scholar] [CrossRef]
  37. Zhang, L.; Danilova, E.V.; Steinbach, I.; Medvedev, D.; Galenko, P.K. Diffuse-interface modeling of solute trapping in rapid solidification: Predictions of the hyperbolic phase-field model and parabolic model with finite interface dissipation. Acta Mater. 2013, 61, 4155–4168. [Google Scholar] [CrossRef]
  38. Lebedev, V.G.; Abramova, E.V.; Danilov, D.A.; Galenko, P.K. Phase-field modeling of solute trapping: Comparative analysis of parabolic and hyperbolic models. Int. J. Mater. Res. 2010, 101, 473–479. [Google Scholar] [CrossRef]
  39. Kharchenko, D.; Lysenko, I.; Galenko, P.K. Fluctuation effects on pattern selection in the hyperbolic model of phase decomposition. Stoch. Differ. Equ. 2011, 97, 97–127. [Google Scholar]
  40. Wolff, F.; Viskanta, R. Solidification of a pure metal at a vertical wall in the presence of liquid superheat. Int. J. Heat Mass Trans. 1988, 31, 1735–1744. [Google Scholar] [CrossRef]
  41. Gau, C.; Viskanta, R.C. Melting and solidification of a pure metal on a vertical wall. J. Heat Transfer 1986, 108, 174–181. [Google Scholar] [CrossRef]
  42. Galenko, P.K.; Funke, O.; Wang, J.; Herlach, D.M. Kinetics of dendritic growth under the influence of convective flow in solidification of undercooled droplets. Mater. Sci. Eng. A 2004, 375–377, 488–492. [Google Scholar] [CrossRef]
  43. Roshchupkina, O.; Shevchenko, N.; Eckert, S. Observation of dendritic growth under the influence of forced convection. IOP Conf. Ser. Mater. Sci. Eng. 2015, 84, 012080. [Google Scholar] [CrossRef]
  44. Huppert, H.E. The fluid mechanics of solidification. J. Fluid Mech. 1990, 212, 209–240. [Google Scholar] [CrossRef]
  45. Chiareli, A.O.P.; Huppert, H.E.; Worster, M.G. Segregation and flow during the solidification of alloys. J. Cryst. Growth 1994, 139, 134–146. [Google Scholar] [CrossRef]
  46. Peppin, S.S.L.; Aussillous, P.; Huppert, H.E.; Worster, M.G. Steady-state mushy layers: Experiments and theory. J. Fluid Mech. 2007, 570, 69–77. [Google Scholar] [CrossRef]
  47. Makoveeva, E.V.; Alexandrov, D.V.; Ivanov, A.A.; Alexandrova, I.V. Desupersaturation dynamics in solutions with applications to bovine and porcine insulin crystallization. J. Phys. A Math. Theor. 2023, 56, 455702. [Google Scholar] [CrossRef]
  48. Makoveeva, E.V.; Ivanov, A.A. Analysis of an integro-differential model for bulk continuous crystallization with account of impurity feeding, dissolution/growth of nuclei and removal of product crystals. Math. Meth. Appl. Sci. 2023. [Google Scholar] [CrossRef]
  49. Makoveeva, E.V. Mathematical modeling of the crystal growth process in a binary system. AIP Conf. Proc. 2020, 2313, 030058. [Google Scholar]
  50. Nizovtseva, I.G.; Starodumov, I.O.; Pavlyuk, E.V.; Ivanov, A.A. Mathematical modeling of binary compounds with the presence of a phase transition layer. Math. Meth. Appl. Sci. 2021, 44, 12260–12270. [Google Scholar] [CrossRef]
Figure 1. A scheme of solid/liquid interface propagation into an undercooled liquid in the presence of a solid wall z = z 0 .
Figure 1. A scheme of solid/liquid interface propagation into an undercooled liquid in the presence of a solid wall z = z 0 .
Mathematics 12 00327 g001
Figure 2. Undercooling at different distances from the flat wall, x = 0 , T 1 = 1680 K, β = d c = 0 .
Figure 2. Undercooling at different distances from the flat wall, x = 0 , T 1 = 1680 K, β = d c = 0 .
Mathematics 12 00327 g002
Figure 3. Undercooling at different distances from the flat wall at fixed P T = 2 × 10 4 , T 1 = 1680 K, β = d c = 0 .
Figure 3. Undercooling at different distances from the flat wall at fixed P T = 2 × 10 4 , T 1 = 1680 K, β = d c = 0 .
Mathematics 12 00327 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Titova, E.A.; Alexandrov, D.V. The Solid–Liquid Phase Interface Dynamics in an Undercooled Melt with a Solid Wall. Mathematics 2024, 12, 327. https://doi.org/10.3390/math12020327

AMA Style

Titova EA, Alexandrov DV. The Solid–Liquid Phase Interface Dynamics in an Undercooled Melt with a Solid Wall. Mathematics. 2024; 12(2):327. https://doi.org/10.3390/math12020327

Chicago/Turabian Style

Titova, Ekaterina A., and Dmitri V. Alexandrov. 2024. "The Solid–Liquid Phase Interface Dynamics in an Undercooled Melt with a Solid Wall" Mathematics 12, no. 2: 327. https://doi.org/10.3390/math12020327

APA Style

Titova, E. A., & Alexandrov, D. V. (2024). The Solid–Liquid Phase Interface Dynamics in an Undercooled Melt with a Solid Wall. Mathematics, 12(2), 327. https://doi.org/10.3390/math12020327

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop