Next Article in Journal
An Optimized Neural Network Acoustic Model for Porous Hemp Plastic Composite Sound-Absorbing Board
Previous Article in Journal
Lie Symmetries and Conservation Laws for the Viscous Cahn-Hilliard Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solutions to Three-Phase-Field Model for Solidification

1
School of Statistics and Applied Mathematic, Anhui University of Finance and Economics, Bengbu 233030, China
2
Institute of Quantitative Economics, Anhui University of Finance and Economics, Bengbu 233030, China
*
Author to whom correspondence should be addressed.
Submission received: 9 March 2022 / Revised: 10 April 2022 / Accepted: 13 April 2022 / Published: 22 April 2022

Abstract

:
Our research on applied mathematics is in line with the research scope of the journal. The phase-field model is applied to simulate the material and other areas. A phase-field model describing the non-isothermal solidification of an ideal multi-component alloy system is proposed in this paper. The time and space variation of a three-phase-field function and the governing equations of the temperature field are established. The global existence of weak solutions for three-dimensional parabolic differential equations is proved by the Faedo–Galerkin method. The existence of a maximum theorem is also extensively studied.

1. Introduction

Our research on applied mathematics is in line with the research scope of the journal. A multi-component alloy system is a kind of important material, especially in technology applications and technology. Therefore, it plays an important role in the formation of mechanical properties and the microstructure of materials. The phase-field model is studied in solid materials (cf. [1]). The multi-component in the alloy combines with the appearance of multi-phase, resulting in different phase transformations and different types of solidification. The solidification of binary alloys is the basis of studying the basic principles of the solidification process. The solidification of a multi-component system can be analyzed by the solidification characteristics of a binary system. In recent years, the multi-phase field (MPF) method has attracted extensive attention. This model is a generalization of Steinbach et al. (cf. [2]), which describes isothermal phase transitions of certain kinds of alloys. The multi-phase field method has been applied to the computer simulation of microstructure evolution (cf. [3,4]) and it was applied to simulate the dendrite solidification process in Fe-Cr-Ni-Mo-C (cf. [5]). It was also applied to the computer simulation of ice-formation in sea-water (cf. [6]). The classic solutions are proved in [7], where the authors set θ = e l ( u + v + w ) . In this paper, we study the non-isothermal solidification of ideal multi-component and multi-phase alloy systems. We set
θ = e l 1 4 ( 1 3 u 3 + u 2 v 1 3 v 3 u v 2 ) l 3 4 ( 1 3 u 3 + u 2 w 1 3 w 3 u w 2 ) l 2 4 ( 1 3 v 3 + v 2 w 1 3 w 3 v w 2 ) .
Thus, we allow a temperature, which is assumed to be a priori given the free energy functional F [ u , v , w , θ ] . This means that the model not only considers the phase transformation caused by the difference in solute, but also considers the phase transformation caused by temperature change.
We studied a system of partial differential equations simulating the evolution of three-phase boundary problems in seawater, ice, and snow, and proved that the system had a global solution in the case of a one-dimensional initial boundary value problem (cf. [8]), where we cannot especially have h ( u , v , w ) .
In this paper we prove the global solution by the Galerkin method and we conclude a maximum theorem.
The order parameter represents the solidification state of the alloy, which is equal to 0 in the solid phase and 1 in the liquid phase. The model is a generalization of the one introduced by Steinbach et al. in [2], for isothermal solidification/melting processes of certain kinds of alloys. Our model must satisfy the following partial differential equations:
u t k 1 ( w Δ u u Δ w ) k 3 ( v Δ u u Δ v ) = a 1 u w ( w u ) a 3 u v ( v u ) + θ ( l 1 u v + l 3 u w ) ,
v t k 2 ( w Δ v v Δ w ) k 3 ( u Δ v v Δ u ) = a 2 v w ( w v ) a 3 v u ( u v ) + θ ( l 1 u v + l 2 v w ) ,
w t k 1 ( u Δ w w Δ u ) k 2 ( v Δ w w Δ v ) = a 1 w u ( u w ) a 2 w v ( v w ) + θ ( l 3 u w l 2 v w ) ,
θ t + ( l 1 u v + l 3 u w ) u t + ( l 1 u v + l 2 v w ) v t + ( l 2 v w l 3 u w ) w t = D Δ θ .
for ( t , x ) ( 0 , T e ) × Ω . The boundary conditions and initial conditions are as follows:
u ( t , x ) = 0 , v ( t , x ) = 0 , w ( t , x ) = 0 , ( t , x ) [ 0 , T e ) × Ω ,
u ( 0 , x ) = u 0 , v ( 0 , x ) = v 0 ( x ) , w ( 0 , x ) = w 0 ( x ) , x Ω .
This system satisfies the second law of thermodynamics. Here, Ω R 3 is a bounded open domain. The function θ is the temperature. The phase-field functions u, v, and w are the fractions of two possible solid and liquid crystal states. For 3 phases we obtain the constraint u + v + w = 1 . For physical reasons, k 1 , k 2 , k 3 , a 1 , a 2 , a 3 , D are positive. In the free energy
F [ u , v , w , θ ] = Ω { k 1 2 | u w w u | 2 + k 2 2 | w v v w | 2 + k 3 2 | u v v u | 2 + ψ ^ ( u , v , w ) 1 4 θ ( l 1 ( 1 3 u 3 + u 2 v 1 3 v 3 u v 2 ) + l 3 ( 1 3 u 3 + u 2 w 1 3 w 3 u w 2 ) + l 2 ( 1 3 v 3 + v 2 w 1 3 w 3 v w 2 ) ) } d x ,
where
ψ ^ ( u , v , w ) = a 1 2 u 2 w 2 + a 2 2 v 2 w 2 + a 3 2 u 2 v 2 ,
we choose for ψ ^ C 2 ( R , [ 0 , ) ) , which represent the double well potential.
θ satisfies
θ = e l 1 4 ( 1 3 u 3 + u 2 v 1 3 v 3 u v 2 ) l 3 4 ( 1 3 u 3 + u 2 w 1 3 w 3 u w 2 ) l 2 4 ( 1 3 v 3 + v 2 w 1 3 w 3 v w 2 ) ,
where e is the local enthalpy and l 1 , l 2 , l 3 are the latent heat of fusion.
It is easy to see in the case of two-phase systems, u = 1 v , or u = 1 w , or v = 1 w , and u v = 1 , or u w = 1 , or v w = 1 .
We write Q T : = ( 0 , T e ) × Ω , where T e is a positive constant, and define
( υ , φ ) Z = Z υ ( y ) φ ( y ) d y ,
for Z = Ω . Since we must have u + v + w = 1 and u t + v t + w t = 0 , the model can be reduced to the following:
u t k 1 ( 1 v ) + k 3 v Δ u ( k 1 k 3 ) u Δ v = a 1 u ( 1 u v ) ( 1 2 u v ) a 3 u v ( v u ) + θ l 1 u v + l 3 u ( 1 u v ) ,
v t k 2 ( 1 u ) + k 3 u Δ v ( k 2 k 3 ) v Δ u = a 2 v ( 1 u v ) ( 1 u 2 v ) a 3 v u ( u v ) + θ l 1 u v + l 2 v ( 1 u v ) ,
θ t D Δ θ = l 1 u v + l 2 v ( 1 u v ) + 2 l 3 u ( 1 u v ) u t
l 1 u v + 2 l 2 v ( 1 u v ) + l 3 u ( 1 u v ) v t .
The boundary and initial conditions are therefore
u ( t , x ) = 0 , v ( t , x ) = 0 , θ ( t , x ) = 0 , ( t , x ) ( 0 , T e ) × Ω ,
u ( 0 , x ) = u 0 ( x ) , v ( 0 , x ) = v 0 ( x ) , θ ( 0 , x ) = θ 0 ( x ) , x Ω .
Definition 1.
Let ( u 0 , v 0 , θ 0 ) H 0 1 ( Ω ) × H 0 1 ( Ω ) × L 2 ( Ω ) . A triple ( u , v , θ ) with
u , v L ( 0 , T e ; H 0 1 ( Ω ) ) L 2 ( 0 , T e ; H 2 ( Ω ) ) ,
θ L ( 0 , T e ; L 2 ( Ω ) ) L 2 ( 0 , T e ; H 0 1 ( Ω ) ) ,
is a weak solution to problems(7)(10), if for all φ C 0 ( ( , T e ) × Ω ) ) , there hold
0 = ( u , φ t ) Q T e α ( u , φ ) Q T e β v , ( u φ ) a 1 u ( 1 u v ) ( 1 2 u v ) , φ Q T e a 3 ( u v ( v u ) , φ ) Q T e + ( u 0 , φ ( 0 ) ) Ω + θ l 1 u v + l 3 u ( 1 u v ) , φ Q T e ,
0 = ( v , φ t ) Q T e γ ( v , φ ) Q T e λ u , ( v φ ) a 2 v ( 1 u v ) ( 1 u 2 v ) , φ Q T e a 3 ( u v ( u v ) , φ ) Q T e + ( v 0 , φ ( 0 ) ) Ω + θ l 1 u v + l 2 v ( 1 u v ) , φ Q T e ,
0 = ( θ , φ t ) Q T e ( D θ , φ ) Q T e + ( θ 0 , φ ( 0 ) ) Ω ( l 1 ( 1 3 u 3 + u 2 v 1 3 v 3 u v 2 ) + l 3 1 3 u 3 + u 2 ( 1 u v ) 1 3 ( 1 u v ) 3 u ( 1 u v ) 2 + l 2 ( 1 3 v 3 + v 2 ( 1 u v ) 1 3 ( 1 u v ) 3 v ( 1 u v ) 2 ) , φ t ) Q T e ( l 1 ( 1 3 u 0 3 + u 0 2 v 0 1 3 v 0 3 u 0 v 0 2 ) + l 3 ( 1 3 u 0 3 + u 0 2 ( 1 u 0 v 0 ) 1 3 ( 1 u 0 v 0 ) 3 u 0 ( 1 u 0 v 0 ) 2 ) + l 2 ( 1 3 v 0 3 + v 0 2 ( 1 u 0 v 0 ) 1 3 ( 1 u 0 v 0 ) 3 v 0 ( 1 u 0 v 0 ) 2 ) , φ ( 0 ) ) Q T e ,
where α = k 1 ( 1 v ) + k 3 v > 0 , β = k 1 k 3 , γ = k 2 ( 1 u ) + k 3 u > 0 , λ = k 2 k 3 , w 0 = 1 u 0 v 0 .
Theorem 1.
For all ( u 0 , v 0 , θ 0 ) H 0 1 ( Ω ) × H 0 1 ( Ω ) × L 2 ( Ω ) , there exists a unique weak solution ( u , v , θ ) of problems(7)(12), which, in addition to(13)(14), satisfies
u t L 2 ( Q T e ) , u L 4 ( Q T e ) , v t L 2 ( Q T e ) , v L 4 ( Q T e ) , θ t L 2 ( 0 , T e ; H 1 ( Ω ) ) .
Notation. In the following sections, we use the letter C > 0 to indicate that it will change with the line. The L 2 ( Ω ) -norm is denoted by · .
The outline of this paper is as follows. In Section 2, we will prove the existence of solutions to the initial boundary-value problem of the nonlinear Equations (7)–(12) by the Faedo–Galerkin method.
In Section 3 we shall show that a maximum theorem holds.

2. Existence of Solutions

In this section, we construct the approximate solutions by the Galerkin method and derive the a priori estimates; then, we propose to send m and to show that a subsequence of our solutions u m , v m , θ m converges to a weak solution of problems (7)–(12).

2.1. Construction of Approximate Solutions

Let { ω k } k = 1 be a basis in H 0 1 ( Ω ) and ω k be a solution to the eigen-problem
Δ ω k = λ k ω k , i n Ω ω k = 0 , k = 1 , , o n Ω .
For a positive integer m, we will look for approximate solutions u m , v m , e m of the form
u m ( t ) = k = 1 m d m k ( t ) ω k , v m ( t ) = k = 1 m g m k ( t ) ω k , e m ( t ) = k = 1 m h m k ( t ) ω k , ( k = 1 , 2 , . . . , m ) ,
where we select the coefficients d m k ( t ) , g m k ( t ) , h m k ( t ) so that
d m k ( 0 ) = ( u 0 , ω k ) = δ m k , g m k ( 0 ) = ( v 0 , ω k ) = η m k , h m k ( 0 ) = ( θ 0 , ω k ) = ζ m k ,
and
( u m t , ω j ) α ( Δ u m , ω j ) β ( u m Δ v m , ω j ) = a 1 u m ( 1 u m v m ) ( 1 2 u m v m ) , ω j a 3 u m v m ( v m u m ) , ω j + ( ( e m l 1 4 ( 1 3 u m 3 + u m 2 v m 1 3 v m 3 u m v m 2 ) l 3 4 ( 1 3 u m 3 + u m 2 w m 1 3 w m 3 u m w m 2 ) l 2 4 ( 1 3 v m 3 + v m 2 w m 1 3 w m 3 v m w m 2 ) ) l 1 u m v m + l 3 u m w m , ω j ) ,
( v m t , ω j ) γ ( Δ v m , ω j ) λ ( v m Δ u m , ω j ) = a 2 ( v m ( 1 u m v m ) ( 1 u m 2 v m ) , ω j ) a 3 ( v m u m ( u m v m ) , ω j ) + ( ( e m l 1 4 ( 1 3 u m 3 + u m 2 v m 1 3 v m 3 u m v m 2 ) l 3 4 ( 1 3 u m 3 + u m 2 w m 1 3 w m 3 u m w m 2 ) l 2 4 ( 1 3 v m 3 + v m 2 w m 1 3 w m 3 v m w m 2 ) ) l 1 u m v m + l 2 v m w m , ω j ) ,
( e m t , ω j ) = D Δ ( θ m l 1 4 ( 1 3 u m 3 + u m 2 v m 1 3 v m 3 u m v m 2 ) l 3 4 ( 1 3 u m 3 + u m 2 w m 1 3 w m 3 u m w m 2 ) l 2 4 ( 1 3 v m 3 + v m 2 w m 1 3 w m 3 v m w m 2 ) , ω j ) , ( j = 1 , , m ) .
where w m = 1 u m v m . (22) and (24) comprise a system of nonlinear ordinary differential equations, and the nonlinear term is locally Lipschitz continuous.
Assuming u m , v m , e m have the structure (19), we note that
( u m ( t ) , ω j ) = ( k = 1 m d m k ( t ) ω k , ω j ) = Ω k = 1 m d m k ( t ) ω k ω j d x = k = 1 m d m k ( t ) Ω ω k ω j d x ,
( v m ( t ) , ω j ) = ( k = 1 m g m k ( t ) ω k , ω j ) = Ω k = 1 m g m k ( t ) ω k ω j d x = k = 1 m g m k ( t ) Ω ω k ω j d x ,
( e m ( t ) , ω j ) = ( k = 1 m h m k ( t ) ω k , ω j ) = Ω k = 1 m h m k ( t ) ω k ω j d x = k = 1 m h m k ( t ) Ω ω k ω j d x ,
We deal with other terms in the same way.
We introduce the vectors
D m = D m ( t ) = d m 1 ( t ) d m m ( t ) , G m = G m ( t ) = g m 1 ( t ) g m m ( t ) , H m = H m ( t ) = h m 1 ( t ) h m m ( t )
and F 1 ( D m T , G m T , H m T ) , F 2 ( D m T , G m T , H m T ) , F 3 ( D m T , G m T , H m T ) denote the nonlinear terms. Thus, we obtain a system of ordinal differential equations.
B D m + F 1 ( D m T , G m T , H m T ) = 0 ,
B G m + F 2 ( D m T , G m T , H m T ) = 0 ,
B H m + F 3 ( D m T , G m T , H m T ) = 0 ,
where
B = ( b i j ) = ( ω 1 , ω 1 ) ( ω 1 , ω m ) ( ω m , ω 1 ) ( ω m , ω m ) = ω 1 ω m ω 1 ω m d x .
Then, if we choose a vector X = ( x 1 , , x m ) that is not equal to zero, there holds for ( X ) T B X = Ω ( x 1 ω 1 + + x m ω m ) 2 d x > 0 ; otherwise, invoking that ω 1 , , ω m are linearly independent, we obtain X = 0 . B is a positive-definite matrix.
For the initial data, we make a smooth approximation
u 0 m = k = 1 m δ m k ω k u 0 , s t r o n g l y i n H 0 1 ( Ω ) ,
v 0 m = k = 1 m η m k ω k v 0 , s t r o n g l y i n H 0 1 ( Ω ) ,
e m = k = 1 m ζ m k ω k e 0 , s t r o n g l y i n H 0 1 ( Ω ) .
According to the existence theorem of local solutions to ordinary differential equations, there exist the solutions for a.e. 0 t t m . We extend t m = T e . Then, we obtain the global solutions. However, we need to show that the a priori estimates hold.

2.2. A Priori Estimates

Theorem 2.
There exists a constant C, depending on Ω , T e , such that
u m H 0 1 ( Ω ) 2 + v m H 0 1 ( Ω ) 2 + e m 2 + 0 T e u m L 4 ( Ω ) 4 d τ + 0 T e v m L 4 ( Ω ) 4 d τ + 0 T e ( u m H 2 ( Ω ) 2 + v m H 2 ( Ω ) 2 + e m H 0 1 ( Ω ) 2 ) d τ C .
for m = 1 , 2 , where u m L ( Ω ) , v m L ( Ω ) , w m L ( Ω ) are suitably small.
Proof. 
Multiplying (22)–(24) by d m k ( t ) , g m k ( t ) , and h m k ( t ) , respectively, summing up over j = 1 , , m , integrating over Ω , and adding and recalling (19), we find
1 2 d d t ( u m 2 + v m 2 + 1 D e m 2 ) + 1 2 γ ( γ α λ 2 ) u m 2 + 1 2 α ( γ α β 2 ) v m 2 + 1 2 e m 2 + 1 64 ( 128 a 1 33 | 3 a 1 a 3 | | 3 a 2 a 3 | 64 ϵ ) u m L 4 ( Ω ) 4 + 1 64 ( 128 a 2 33 | 3 a 2 a 3 | | 3 a 1 a 3 | 64 ϵ ) v m L 4 ( Ω ) 4 + 1 32 ( 32 a 1 + 64 a 3 + 32 a 2 15 | 3 a 1 a 3 | 15 | 3 a 2 a 3 | ) u m v m 2 C ( u m 2 + v 2 + 1 D e m 2 ) + C ( Ω ) ,
where α γ λ 2 > 0 , α γ β 2 , 128 a 1 33 | 3 a 1 a 3 | | 3 a 2 a 3 | 64 ϵ > 0 , 128 a 2 33 | 3 a 2 a 3 | | 3 a 1 a 3 | 64 ϵ > 0 , 32 a 1 + 64 a 3 + 32 a 2 15 | 3 a 1 a 3 | 15 | 3 a 2 a 3 | > 0 .
By using the Gronwall inequality, we obtain
( u m 2 + v m 2 + 1 D e m 2 ) + 1 2 γ ( γ α λ 2 ) 0 T e u m 2 d τ + 1 2 α ( γ α β 2 ) 0 T e v m 2 d τ + 1 2 0 T e e m 2 d τ + 1 64 ( 128 a 1 33 | 3 a 1 a 3 | | 3 a 2 a 3 | 64 ϵ ) 0 T e u m L 4 ( Ω ) 4 d τ + 1 64 ( 128 a 2 33 | 3 a 2 a 3 | | 3 a 1 a 3 | 64 ϵ ) 0 T e v m L 4 ( Ω ) 4 d τ + 1 32 ( 32 a 1 + 64 a 3 + 32 a 2 15 | 3 a 1 a 3 | 15 | 3 a 2 a 3 | ) 0 T e u m v m 2 d τ C ( T e ) .
Multiplying (22) and (23) by λ j d m k ( t ) and λ j g m k ( t ) , respectively, summing up over j = 1 , , m , then formally integrating over Ω and adding them, we obtain
1 2 d d t ( u m 2 + v m 2 ) + 1 2 γ ( α γ λ 2 γ ϵ ) Δ u m 2 + 1 2 α ( γ α β 2 α ϵ ) Δ v m 2 + 1 2 ( 6 a 1 | 3 a 1 a 3 | | 5 a 1 3 a 3 a 2 | 2 ) u m u m 2 + 1 2 ( 2 a 1 + 2 a 3 | 3 a 2 a 3 | | 5 a 1 3 a 3 a 2 | 2 ) v m u m 2 + 1 2 ( 6 a 2 | 3 a 2 a 3 | | 5 a 2 3 a 3 a 1 | 2 ) v m v m 2 + 1 2 ( 2 a 2 + 2 a 3 | 3 a 1 a 3 | | 5 a 2 3 a 3 a 1 | 2 ) u m v m 2 C ( u m 2 + v m 2 + e m 2 + 1 ) ,
where α γ λ 2 γ ϵ > 0 , γ α β 2 α ϵ > 0 , 6 a 1 | 3 a 1 a 3 | | 5 a 1 3 a 3 a 2 | 2 > 0 , 2 a 1 + 2 a 3 | 3 a 2 a 3 | | 5 a 1 3 a 3 a 2 | 2 > 0 , 2 a 2 | 3 a 2 a 3 | | 5 a 2 3 a 3 a 1 | 2 > 0 , 2 a 2 + 2 a 3 | 3 a 1 a 3 | | 5 a 2 3 a 3 a 1 | 2 > 0 .
By using the Gronwall inequality, we obtain
u m 2 + v m 2 + 1 2 γ ( α γ λ 2 γ ϵ ) 0 T e Δ u m 2 d τ + 1 2 α ( γ α β 2 α ϵ ) 0 T e Δ v m 2 d τ C ( T e ) .
It follows from (37) and (38) that
u m t L 2 ( Q T e ) + L 4 3 ( Q T e ) + v m t L 2 ( Q T e ) + L 4 3 ( Q T e ) + e m t L 2 ( 0 , T e ; H 1 ( Ω ) ) C ( T e ) .

2.3. Existence of Weak Solutions

Next we pass to limits as m , to build the weak solutions to the initial boundary-value problems (7)–(12).
Theorem 3
(Aubin–Lions). Let B 0 be a normed linear space embedded compactly into another normed linear space B, which is continuously embedded into a Hausdorff locally convex space B 1 , and 1 p < + . If v , v i L p ( 0 , t ; B 0 ) , i N , the sequence { v i } i N converges weakly to v in L p ( 0 , t ; B 0 ) , and { v i t } i N is bounded in L 1 ( 0 , t ; B 1 ) , then v i converges to v strongly in L p ( 0 , t ; B ) .
Lemma 1.
Let ( 0 , t ) × Ω be an open set in R + × R n . Suppose functions g n , g are in L q ( ( 0 , t ) × Ω for any given 1 < q < , which satisfy
g n L q ( ( 0 , t ) × Ω ) C , g n g a . e . i n ( 0 , t ) × Ω .
Then g n converges to g weakly in L q ( ( 0 , t ) × Ω ) .
Theorem 3 is a general version of the Aubin–Lions lemma valid under the weak assumption t v i L 1 ( 0 , t ; B 1 ) . This version, which we need here, is proved in [9]. A proof of Lemma 1 can be found in [9].
Lemma 2.
Problems (7)(12) have at least one weak solution ( u , v , e ) in the sense of Definition 1.1. Each of the weak solutions satisfies e t L 2 ( 0 , T e ; H 1 ( Ω ) ) and ( u t , v t ) L 2 ( 0 , T e ; L 2 ( Ω ) ) + L 4 3 ( Q T e ) .
Proof of Lemma 2.
According to the energy estimates (34), we see sequences { u m } m = 1 , { v m } m = 1 , { e m } m = 1 are bounded in L ( 0 , T e ; H 0 1 ( Ω ) ) L 2 ( 0 , T e ; H 2 ( Ω ) ) ,
L ( 0 , T e ; L 2 ( Ω ) ) L 2 ( 0 , T e ; H 0 1 ( Ω ) ) , respectively, and ( { u m t } m = 1 , { v m t } m = 1 ) , { e m t } m = 1 are bounded in L 2 ( 0 , T e ; L 2 ( Ω ) ) + L 4 3 ( Q T e ) , L 2 ( 0 , T e ; H 1 ( Ω ) ) , respectively.
Consequently there exist subsequences { u m l } l = 1 { u m } m = 1 , { v m l } l = 1 { v m } m = 1 , { e m l } l = 1 { e m } m = 1 and functions ( u , v ) L ( 0 , T e ; H 0 1 ( Ω ) ) L 2 ( 0 , T e ; H 2 ( Ω ) ) ,
e L 2 ( 0 , T e ; H 0 1 ( Ω ) ) L ( 0 , T e ; L 2 ( Ω ) ) , with ( u t , v t ) L 2 ( 0 , T e ; L 2 ( Ω ) ) + L 4 3 ( Q T e ) ,
e t L 2 ( 0 , T e ; H 1 ( Ω ) ) , such that
u m l u , w e a k l y i n L 2 ( 0 , T e ; H 2 ( Ω ) ) u m l * u , w e a k l y i n L ( 0 , T e ; H 0 1 ( Ω ) ) u m l t u t , w e a k l y i n L 2 ( 0 , T e ; L 2 ( Ω ) ) + L 4 3 ( Q T e ) v m l v , w e a k l y i n L 2 ( 0 , T e ; H 2 ( Ω ) ) v m l * v , w e a k l y i n L ( 0 , T e ; H 0 1 ( Ω ) ) v m l t v t , w e a k l y i n L 2 ( 0 , T e ; L 2 ( Ω ) ) + L 4 3 ( Q T e ) . e m l * e , w e a k l y i n L ( 0 , T e ; L 2 ( Ω ) ) e m l e , w e a k l y i n L 2 ( 0 , T e ; H 0 1 ( Ω ) ) e m l t e t , w e a k l y i n L 2 ( 0 , T e ; H 1 ( Ω ) )
There exist functions ξ , ϕ such that
u m 3 ξ , w e a k l y i n L 4 3 ( Q T e )
v m 3 ϕ , w e a k l y i n L 4 3 ( Q T e ) .
It remains to show that ξ = u 3 , ϕ = v 3 . To this end, we show first that u m u , v m v in L 2 ( 0 , T e ; H 1 ( Ω ) ) by applying Theorem 3. To apply this theorem, u t , v t have estimates u t L 2 ( 0 , T e ; L 2 ( Ω ) ) + L 4 3 ( Q T e ) L 4 3 ( Q T e ) L 1 ( 0 , T e ; L 4 3 ( Ω ) ) , v t L 2 ( 0 , T e ; L 2 ( Ω ) ) + L 4 3 ( Q T e ) L 4 3 ( Q T e ) L 1 ( 0 , T e ; L 4 3 ( Ω ) ) , e t L 2 ( 0 , T e ; H 1 ( Ω ) ) .
Applying Theorem 3 with B 0 = H 2 ( Ω ) , B = H 1 ( Ω ) , B 1 = L 4 3 ( Ω ) , and p = 2 we obtain
u m u , s t r o n g l y i n L 2 ( 0 , T e ; H 1 ( Ω ) ) v m v , s t r o n g l y i n L 2 ( 0 , T e ; H 1 ( Ω ) ) .
Therefore, we choose the other sequences from these sequences and denote them in the same way. They converge almost everywhere in Q T e . This implies the convergence u m 3 u 3 , v m 3 v 3 almost everywhere in Q T e . Using the embedding H 1 L 4 ( Ω ) and applying Lemma 1, we obtain ξ = u 3 , ϕ = v 3 . Then, we obtain the others in the same way. Equation (15) follows from these relations if we show that
(42) ( u 0 m , φ ( 0 ) ) Ω ( u 0 , φ ( 0 ) ) , (43) ( u m , φ t ) Q T e ( u , φ t ) Q T e , (44) ( u m , φ ) Q T e ( u , φ ) Q T e , (45) ( v m , ( u m φ ) ) Q T e ( v m , ( u φ ) ) Q T e , (46) u m v m ( v m u m ) , φ Q T e ( u v ( v u ) , φ ) Q T e , (47) ( θ m u m v m , φ ) Q T e ( θ u v , φ ) Q T e .
for m . Now, the relation (42) follows from (31), the relation (43) is a consequence of u m t L 4 3 ( Q T e ) , the relation (44) is consequence of u m L 2 ( 0 , T e ; H 2 ( Ω ) ) , the relation (45) is a consequence of u m u , u m u , v m v in L 2 ( Q T e ) , the relation (46) is a consequence of u m v m L 2 ( Q T e ) and u m u , v m v in L 2 ( Q T e ) , and the relation (47) is obtained from u m v m L 2 ( Q T e ) and θ m θ in L 2 ( Q T e ) . We obtain other terms in the same way as above. □

2.4. Uniqueness

In this subsection we show uniqueness of the solution of ( u , v , e ) that was obtained in Section 2.2 and Section 2.3.
Theorem 4.
Let ( u 0 i , v 0 i ) H 0 1 ( Ω ) , e 0 i L 2 ( Ω ) , i = 1 , 2 be given functions. Let ( u i , v i , e i ) be weak solutions of problems(7)(12)with ( u i , v i ) L ( 0 , T e ; H 0 1 ( Ω ) ) L 2 ( 0 , T e ; H 2 ( Ω ) ) , and e i L 2 ( 0 , T : H 0 1 ( Ω ) ) L ( 0 , T e ; L 2 ( Ω ) ) , i = 1 , 2 . Then
u ˜ 2 + v ˜ 2 + e ˜ 2 C ( 1 + T e e C 1 T e ) ( u ˜ 0 2 + v ˜ 0 2 + e ˜ 0 2 ) .
where u ˜ = u 1 u 2 , v ˜ = v 1 v 2 and e ˜ = e 1 e 2 . The constant C is independent of u i , v i , e i , u 0 i , v 0 i , e 0 i .
Proof. 
The result equation of the difference of (7) for u 1 and u 2 is multiplied by u ˜ . We obtain the others in the same way. Then, u ˜ , v ˜ , and e ˜ satisfy the inequality
1 2 ( u ˜ 2 + v ˜ 2 + 1 D e ˜ 2 ) + 1 2 γ ( γ α λ 2 ) u ˜ 2 + 1 2 α ( γ α β 2 ) v ˜ 2 + 1 4 e ˜ 2 C { ( u ˜ 0 2 + v ˜ 0 2 + e ˜ 0 2 ) + 0 t Ω ( u ˜ 2 + v ˜ 2 + e ˜ 2 ) } .
We use the Gronwall inequality, which yields
u ¯ 2 + v ¯ 2 + e ¯ 2 C ( 1 + T e e C 1 T e ) ( u ¯ 0 2 + v ¯ 0 2 + e ¯ 0 2 ) .
In particular, for u ¯ 0 = v ¯ 0 = e ¯ 0 = 0 , we obtain the uniqueness of the solution.

3. A Maximum Theorem

If we want the model to have physical meaning, we must prove that the order parameters u , v , w still change in the interval [ 0 , 1 ] during its evolution. We will study a maximum theorem.
Theorem 5.
Let ( u , v , w , θ ) be the solution of problems(7)(12)and u 0 [ 0 , 1 ] , v 0 [ 0 , 1 ] , w 0 [ 0 , 1 ] for each x Ω . Then such solution satisfies u [ 0 , 1 ] , v [ 0 , 1 ] , w [ 0 , 1 ] and u + v + w = 1 for each x Ω and t R + .
Proof. 
We only analyze Equation (7); Equation (8) is treated in a similar way. To prove that u ( t , x ) 0 , let us define
u = u , if u < 0 , 0 , if u 0 ,
so that u 0 and satisfies boundary and initial conditions
u ( t , x ) = 0 , t 0 , x Ω , u ( 0 , x ) = 0 , x Ω .
We multiply (7) by u and integrate in Ω to obtain
1 2 d d t Ω ( u ) 2 d x + α 2 Ω ( u ) 2 d x C Ω ( v ) 2 ( u ) 2 d x + C Ω θ ( u ) 2 d x a 1 Ω ( u ) 2 ( 1 u v ) ( 1 2 u v ) d x a 3 Ω ( u ) 2 v ( v u ) d x .
Using Gronwall’ s lemma, we obtain that ( u ) 2 = 0 a.e. in Ω for all t ( 0 , T e ) , and thus u 0 a.e. in Q T e . We obtain in a similar way as earlier that v 0 , w 0 a.e. in Q T e . Using a similar argument to prove that u ( t , x ) 1 , v ( t , x ) 1 , w ( t , x ) 1 , we define
u + = u 1 , if u > 1 , 0 , if u 1 ,
v + = v 1 , if v > 1 , 0 , if v 1 ,
w + = w 1 , if w > 1 , 0 , if w 1 ,
so that u + ( t , x ) 0 and satisfies boundary and initial conditions
u + ( t , x ) = 0 , v + ( t , x ) = 0 , w + ( t , x ) = 0 , t 0 , x Ω , u + ( 0 , x ) = 0 , v + ( t , x ) = 0 , w + ( t , x ) = 0 , x Ω .
Multiplying (1)–(3) by u + , v + , and w + , respectively, integrating in Ω , and adding the resulting equations, we have
1 2 d d t Ω u + ( 1 + u + ) d x + 1 2 d d t Ω v + ( 1 + v + ) d x + 1 2 d d t Ω w + ( 1 + w + ) d x + α Ω ( u + ) 2 d x + γ Ω ( v + ) 2 d x + k 1 u + k 2 ( 1 u ) Ω ( w + ) 2 d x | β | Ω Δ v u + ( 1 + u + ) d x + C Ω Δ u v + ( 1 + v + ) d x + C Ω Δ u w + ( 1 + w + ) d x a 1 Ω u + ( 1 + u + ) ( 1 u + v + ) ( 2 2 u + v + ) d x + C Ω θ u + ( 1 + u + ) d x a 3 ( Ω u + ( 1 + u + ) ( 1 + v + ) ( v + u + ) d x + C Ω θ v + ( 1 + v + ) d x + Ω v + ( 1 + v + ) ( 1 + u + ) ( u + v + ) d x ) + C Ω w + ( 1 + w + ) θ d x a 2 Ω v + ( 1 + v + ) ( 1 u + v + ) ( 2 2 v + u + ) d x a 1 Ω w + ( 1 + w + ) ( u + w + ) d x a 2 Ω w + ( 1 + w + ) ( 1 + v + ) ( v + w + ) d x .
By using Gronwall’s lemma, we obtain that u + = 0 , v + = 0 , w + = 0 a.e. in Ω for all t ( 0 , T e ) , and thus u 1 , v 1 , w 1 a.e. in Q T e . □

Author Contributions

Methodology, Y.T.; Writing—original draft, W.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Science Foundation of Anhui Province under grant No. 2108085MA04.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author gratefully acknowledges the help of Peicheng Zhu, Shanghai University. He has offered valuable suggestions in this study. The author would like to deeply thank all the reviewers for their insightful and constructive comments. This paper has been partly supported by the Grant 2108085MA04 of the key NSF of Anhui Province, China.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, L.; Mai, J.; Cheng, F. Weak solutions and simulations to a new phase-field model with periodic boundary. Appl. Anal. 2022, 101, 366–381. [Google Scholar] [CrossRef]
  2. Steinbach, I.; Pezzolla, F.; Nestler, B.; Seeßelberg, M.; Prieler, R.; Schmitz, G.J.; Rezende, J.L. A phase field concept for multiphase systems. Phys. D 1996, 94, 135–147. [Google Scholar] [CrossRef]
  3. Zhang, J.; Zheng, C.; Li, D. A multi-phase field model for static recrystallization of hot deformed austenite in a C-Mn steel. Acta Metall. Sin. (Engl. Lett.) 2018, 31, 208–215. [Google Scholar] [CrossRef]
  4. Zhu, C.S.; Wang, T.; Feng, L.; Lei, P.; Ma, F.L. Numerical study of growth competition between twin grains during directional solidification by using multi-phase field method. Chin. Phys. B 2022, 31, 028102. [Google Scholar] [CrossRef]
  5. Nomoto, S.; Mori, K.; Segawa, M. Solidification Simulation of Fe-Cr-Ni-Mo-C Duplex Stainless Steel Using CALPHAD-Coupled Multi-phase Field Model with Finite Interface Dissipation. In Proceedings of the 4th World Congress on Integrated Computational Materials Engineering (ICME 2017), Ypsilanti, MI, USA, 21–25 May 2017; pp. 283–292. [Google Scholar]
  6. Covello, V.; Abbà, A.; Bonaventyra, L.; Rocca, A.D.; Valdettaro, L. A multiphase model for the numerical simulation of ice-formation in sea-water. In Proceedings of the VII European Congress on Computational Methods in Applied Sciences and Engineering, Crete, Greece, 5–10 June 2016. [Google Scholar]
  7. Morelli, B.; Caretta, C.; Boldrini, J.L. Local existence of solutions of a three phase-field model for solidification. Math. Method. Appl. Sci. 2009, 32, I2. [Google Scholar]
  8. Tang, Y.; Zheng, L.; Luan, L. Solutions to a multi-phase model of sea ice growth. Open Math. 2021, 19, 1–19. [Google Scholar] [CrossRef]
  9. Lions, J. Quelques Methodes de Resolution des Problems aux Limites Non Lineaires; Dunod Gauthier-Villars: Pairs, France, 1969. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tang, Y.; Gao, W. Solutions to Three-Phase-Field Model for Solidification. Symmetry 2022, 14, 862. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14050862

AMA Style

Tang Y, Gao W. Solutions to Three-Phase-Field Model for Solidification. Symmetry. 2022; 14(5):862. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14050862

Chicago/Turabian Style

Tang, Yangxin, and Wei Gao. 2022. "Solutions to Three-Phase-Field Model for Solidification" Symmetry 14, no. 5: 862. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14050862

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