Research article Topical Sections

Modeling of the continuous casting process of steel via phase-field transition system. Fractional steps method

  • Here we consider the phase field transition system (a nonlinear system of parabolic type), introduced by G. Caginalp to distinguish between the phases of the material that is involved in the solidification process. On the basis of the convergence of an iterative scheme of fractional steps type, a conceptual numerical algorithm is elaborated in order to approximate the solution of the nonlinear parabolic problem. The advantage of such approach is that the new method simplifies the numerical computations due to its decoupling feature. The finite element method (fem) in 2D is used to deduce the discrete equations and numerical results regarding the physical aspects of solidification process are reported. In order to refer the continuous casting process, the adequate boundary conditions was considered.

    Citation: Costică Moroşanu. Modeling of the continuous casting process of steel via phase-field transition system. Fractional steps method[J]. AIMS Mathematics, 2019, 4(3): 648-662. doi: 10.3934/math.2019.3.648

    Related Papers:

    [1] Mingliang Liao, Danxia Wang, Chenhui Zhang, Hongen Jia . The error analysis for the Cahn-Hilliard phase field model of two-phase incompressible flows with variable density. AIMS Mathematics, 2023, 8(12): 31158-31185. doi: 10.3934/math.20231595
    [2] Lakhlifa Sadek, Tania A Lazǎr, Ishak Hashim . Conformable finite element method for conformable fractional partial differential equations. AIMS Mathematics, 2023, 8(12): 28858-28877. doi: 10.3934/math.20231479
    [3] Jaeyong Choi, Seokjun Ham, Soobin Kwak, Youngjin Hwang, Junseok Kim . Stability analysis of an explicit numerical scheme for the Allen-Cahn equation with high-order polynomial potentials. AIMS Mathematics, 2024, 9(7): 19332-19344. doi: 10.3934/math.2024941
    [4] Martin Stoll, Hamdullah Yücel . Symmetric interior penalty Galerkin method for fractional-in-space phase-field equations. AIMS Mathematics, 2018, 3(1): 66-95. doi: 10.3934/Math.2018.1.66
    [5] Yasir Nawaz, Muhammad Shoaib Arif, Wasfi Shatanawi, Muhammad Usman Ashraf . A new unconditionally stable implicit numerical scheme for fractional diffusive epidemic model. AIMS Mathematics, 2022, 7(8): 14299-14322. doi: 10.3934/math.2022788
    [6] Saulo Orizaga, Ogochukwu Ifeacho, Sampson Owusu . On an efficient numerical procedure for the Functionalized Cahn-Hilliard equation. AIMS Mathematics, 2024, 9(8): 20773-20792. doi: 10.3934/math.20241010
    [7] Jing Li, Linlin Dai, Kamran, Waqas Nazeer . Numerical solution of multi-term time fractional wave diffusion equation using transform based local meshless method and quadrature. AIMS Mathematics, 2020, 5(6): 5813-5838. doi: 10.3934/math.2020373
    [8] Zhichao Fang, Ruixia Du, Hong Li, Yang Liu . A two-grid mixed finite volume element method for nonlinear time fractional reaction-diffusion equations. AIMS Mathematics, 2022, 7(2): 1941-1970. doi: 10.3934/math.2022112
    [9] Alessandra Jannelli, Maria Paola Speciale . On the numerical solutions of coupled nonlinear time-fractional reaction-diffusion equations. AIMS Mathematics, 2021, 6(8): 9109-9125. doi: 10.3934/math.2021529
    [10] Mubashara Wali, Sadia Arshad, Sayed M Eldin, Imran Siddique . Numerical approximation of Atangana-Baleanu Caputo derivative for space-time fractional diffusion equations. AIMS Mathematics, 2023, 8(7): 15129-15147. doi: 10.3934/math.2023772
  • Here we consider the phase field transition system (a nonlinear system of parabolic type), introduced by G. Caginalp to distinguish between the phases of the material that is involved in the solidification process. On the basis of the convergence of an iterative scheme of fractional steps type, a conceptual numerical algorithm is elaborated in order to approximate the solution of the nonlinear parabolic problem. The advantage of such approach is that the new method simplifies the numerical computations due to its decoupling feature. The finite element method (fem) in 2D is used to deduce the discrete equations and numerical results regarding the physical aspects of solidification process are reported. In order to refer the continuous casting process, the adequate boundary conditions was considered.


    We consider a material in a region ΩRn, n3, which may be in either of two phases, e.g., solid and liquid (see Figure 1). Let us denote by

    u(t,x)=θ(t,x)θM,    (t,x)Q=(0,T)×Ω, T>0,
    Figure 1.  A material Ω exists in two phases. The dotted lines indicate possible thickness of the continuous region.

    the reduced temperature distribution, where θ(t,x) represent the temperature of the material and θM is the melting temperature (the temperature at which solid and liquid may coexist in equilibrium, separated by an interface).

    In the following we will describe the framework of our problem. So, let us consider the interface as a continuous region, more vast (in which the liquid can coexist with the solid) and of finite thickness, in which the change of phase occurring continuously.

    The following nonlinear parabolic system

    {Cρtu+2tφ=kΔuinQ, αξtφ=ξΔφ+12ξ(φφ3)+sξu (1.1)

    with the non-homogeneous Cauchy-Neumann boundary conditions

    {kνu+hu=w1(t,x)onΣ=(0,T]×Ω,ξνφ=w2(t,x) (1.2)

    and with the initial conditions

    u(0,x)=u0(x),φ(0,x)=φ0(x)onΩ, (1.3)

    represents the mathematical model called the phase field transition system, introduced by G. Caginalp (see [3] and references therein) to model the transition between the solid and liquid phase in melting/solidification process to a matter occupying a region Ω, while:

    u(t,x) - represents the reduced temperature distribution in Q;

    φ(t,x) - is the phase function (the order parameter) used to distinguish between the states (phases) of material which occupies the region Ω at every time t[0,T];

    Cρ=ρV; ρ - the density, V - the casting speed;

    ,k,α,ξ,h are physical parameters representing, respectively: the latent heat, the thermal conductivity, the relaxation time, the measure of the interface thickness, the heat transfer coefficient;

    sξ=m[S]E2σTE, a bounded and positive quantity, expressed by positive and bounded physical parameters: m=11(2F(s))12ds, F(s)=14(s21)2, [S]E - the entropy difference between phases per volume, σ - the interfacial tension, TE - the equilibrium melting temperature (see Caginalp & Chen [4]);

    w1(t,x)W112p,21pp(Σ), p2, - is a given function: the temperature of the surrounding at Ω for each time t[0,T] (can also be interpreted as boundary control);

    w2(t,x)W112p,21pp(Σ) - is a given function;

    u0,φ0W22pp(Ω), with kνu0+hu0=w1(0,x) and ξνφ0=w2(0,x);

    The model (1.1)–(1.3) represents a refinement of the classical Stefan problem (see [21,23,24]) in two phases by adding a new nonlinear parabolic equation, derived from the Euler-Lagrange equations for the free energy in Landau-Ginzburg field theory. This new mathematical model reflects more accurately the physical phenomenon of solidification, like: superheating, supercooling, etc.

    Different other nonlinearities capable to come out the complexity of the physical phenomena (the effect of surface tension, separating zone of solid and liquid states, etc) have been proposed by several authors (see Cârjă, Miranville & Moroşanu [5], Kenmochi & Niezgódka [8], Miranville & Moroşanu [9], Moroşanu [16], Moroşanu & Motreanu [20], Penrose & Fife [22] and Temam [24]). The general nonlinear term in Moroşanu & Motreanu [20], is (possibly) non-convex and non-monotone and cover a large class of nonlinearities, including the known cases as well as other new relevant situations. Moreover, different types of boundary conditions on Σ can be associated to (1.1) (see Moroşanu [16] for more details).

    In the present section we will investigate the solvability of the first boundary value problems of the form (1.1)–(1.3) in the class W1,2p(Q), p2.

    The main result of this section establishes the dependence of the solution u(t,x), φ(t,x) in the nonlinear parabolic system (1.1)–(1.3) on the terms w1(t,x), w2(t,x) in the right-hand side of (1.2).

    Theorem 2.1 Problem (1.1)(1.3) has a unique solution (u,φ) with uW1,2p(Q) and φW1,2ν(Q), where ν=min{q,μ}, qp2. In addition (u,φ) satisfies

    uW1,2p(Q)+φW1,2ν(Q)C[1+u0W22pp(Ω)+φ032pW22qq(Ω)+w1W112p,21pp(Σ)+w2W112p,21pp(Σ)], (2.1)

    where the constant C depends on |Ω| (the measure of Ω), T, n, p, q and physical parameters, but is independent of u, φ, w1 and w2.

    Moreover, given any number Md>0, if (u1,φ1) and (u2,φ2) are solutions to (1.1)(1.3) for the same initial conditions, corresponding to the data w11,w21, w12,w22 W112p,21pp(Σ), respectively, such that φ1Lν(Q), φ2Lν(Q)Md, then the estimate below holds

    u1u2W1,2p(Q)+φ1φ2W1,2ν(Q)C[w11w21W112p,21pp(Σ)+w12w22W112p,21pp(Σ)], (2.2)

    where the constant C depends on |Ω|, T, Md, n, p, q and physical parameters, but is independent of u1,u2, φ1,φ2, w11, w21, w12 and w22.

    Proof. The basic tools in the analysis of the problem (1.1) (see [16] and references there in) are the Leray-Schauder degree theory, the Lp-theory of linear and quasi-linear parabolic equations, as well as the Lions and Peetre embedding Theorem, which ensures the existence of a continuous embedding W1,2p(Q)Lμ(Q), where the number μ is defined as follows

    μ={if   p>32,anypositivenumber3pif   p=32.

    The proof of Theorem 2.1 was given in Moroşanu [16] noting that there formulation differs from this by certain physical parameters, which implies different values for the constant C in (2.1) and (2.2). Moreover, corresponding to different boundary conditions (including nonlinear and nonhomogeneous boundary conditions), similar results were proved in Cârjă, Miranville & Moroşanu [5] and Miranville & Moroşanu [9].

    Corollary 2.2 Under hypotheses H0 and H2 in [20] the problem (1.1) possesses a unique solution (u,φ)W1,2p(Q)×W1,2p(Q).

    Proof. Let w11=w21=w12=w22=w in the Theorem 2.1. Then (2.2) shows that the conclusion of the corollary is true.

    The aim of this section is to use the fractional steps method in order to approximate the solution of the system (1.1)–(1.3), whose uniqueness is guaranteed by Corollary 2.2. To do that, let us associate to the time-interval [0,T] the equidistant grid of length ε=TM, for any integer M1. Then, the following approximating scheme can be written in order to approximate the solution of the nonlinear boundary value problem (1.1)–(1.3):

    {ρVtuε+2tφεkΔuε=0 in Qεi,kνuε+huε=w1(t,x) on Σεi, (3.1)
    {αξtφεξΔφε=12ξφε+sξuε in Qεi,ξνφε=w2(t,x) on Σεi,φε+(iε)=z(ε,φε(iε)), (3.2)

    where Qεi=[iε,(i+1)ε]×Ω, Σεi=[iε,(i+1)ε]×Ω, i=0,,M1, and z(t,φε(iε)) is the solution of the Cauchy problem

    {z(s)+12ξz(s)3=0    s[0,T],z(0)=φε(iε),    φε(0,x)=φ0(x), (3.3)

    computed at s=ε, for i=0,,M1. Here φε+(iε)=limtiεφε(t) and φε(iε)=limtiεφε(t).

    For later use, we set: W=L2(0,T;H1(Ω))W1,2([0,T];(H1(Ω))).

    Definition 3.1. By weak solution to the problem (1.1)–(1.2) we mean a pair of functions (u,φ)W×W which satisfy (1.1)–(1.2) in the following sense

    Qt(ρVu+2φ)ψdxdt+kQuψdxdt+hΣuψdγdt=Σw1ψdγdt, (3.4)
    Q(αξtφζ+ξφζ12ξ(φφ3)ζsξuζ)dxdt=Σw2ζdγdt, (3.5)

    for all (ψ,ζ)L2(0,T;H1(Ω))×L2(0,T;H1(Ω)), together with the initial conditions (1.3). In (3.4) and (3.5) we have denoted by the symbol Q the duality between L2(0,T;H1(Ω)) and L2(0,T;(H1(Ω))).

    The main result of this section is the following

    Theorem 3.2. Assume that u0,φ0W22pp(Ω) with kνu0+hu0=w1(0,x) and ξνφ0=w2(0,x). Let (uε,φε) be the solution of the approximating scheme (3.1)(3.3). Then, for ε0, we have

    (uε(t),φε(t))(u(t),φ(t))stronglyinL2(Ω) (3.6)

    for any t[0,T], where u,φW1,2([0,T];L2(Ω))L(0,T;H2(Ω)) is the weak solution to the problem (1.1)(1.3).

    Proof.(see [2]) The proof is based on compactness methods. As a matter of fact it turns out from Theorem 3.2 that if u0,φ0L2(Ω), then the weak solution (u(t),φ(t)) of the system (1.1)–(1.3) is a strong solution, i.e., it is absolutely continuous in t on [0,T] and satisfies a.e. the system (1.1)–(1.3). So Theorem 3.2 can be also viewed as a constructive way to prove the existence in (1.1)–(1.3).

    The result in Theorem 3.2 remains true by replacing the boundary condition (1.2) with kνu+hu=w(t)g(x) and ξνφ=0 (see [6,7,10,11,12,13,14,15,16,17,18,19]).

    The Cauchy problem (3.3) has the solution

    z(ε,φε(iε,x))=|φε(iε,x)|ξξ+ε(φε(iε,x))2,    i=0,,M1, (3.7)

    and then the general algorithm to compute the approximate solution by means of fractional steps method consist in the following sequence (i denotes the time level)

    Begin algfrac

      i:=0 uε,0=u0, φε,0=φ0 from the initial conditions (1.3);

      For i:=0 to M1 do

        Compute z(ε,φε(iε,x)) from (3.7);

        φε+:=z(ε,φε(iε,x));

         Compute φε,i+1,uε,i+1 solving the linear system (3.1)-(3.2);

      Endfor;

    End.

    A comparison between the fractional steps method and the standard iterative Newton method can be found in Moroşanu [16].

    The finite element method (fem in short) is a general method for approximating the solution of boundary value problems for partial differential equations. This method derives from the Ritz (or Gelerkin) method (see Axelson & Barker [1]), characteristic for the finite element method being the chose of the finite dimensional space, namely, the span of a set of finite element basis functions.

    The steps in solving a boundary value problem using fem are:

    P0. (D) The direct formulation of the problem;

    P1. (V) A variational (weak) formulation for problem (D);

    P2.     The construction of a finite element mesh (triangulation);

    P3.     The construction of the finite dimensional space of test

        function (called finite element basis functions);

    P4. (Vnn) A discrete analogous of (V);

    P5.     Assembly of the system of linear equations;

    P6.     Solve the system in P5.

    The finite element method is used in the sequel in order to deduce the discrete state equations. A conceptual numerical algorithm of fractional step type is then formulated to approximate the weak solution corresponding to (3.1)–(3.2), that is:

    (ρVuεt+2φεt,ψ)+k(uε,ψ)+hΩuεψdxdy=Ωw1ψdxdy, (4.1)
    ψH1(Ω),a.e. in (iε,(i+1)ε),
    αξ(φεt,ζ)+ξ(φε,ζ)12ξ(φε,ζ)=sξ(uε,ζ)+Ωw2ζdxdy, (4.2)
    ζH1(Ω),a.e. in (iε,(i+1)ε),

    together with the initial conditions

    u(0,x)=u0(x),φ(0,x)=φ0(x),xΩ.

    By (,) we have denoted the scalar product in L2(Ω).

    Let ε=T/M be the time step size. We assume that ΩR2 is a polygonal domain. Let Tr be the triangulation (mesh) over Ω and ˉΩ=KTrK, and let Nj=(xk,yl),j=¯1,nn, be the nodes associated to Tr. Denoting by Vnn the corresponding finite element space to Tr, then the basic functions {bj}nnj=1 of Vnn are defined by

    bj(Ni)=δji,i,j=¯1,nn,

    and so

    Vnn=span{b1,b2,...,bnn}.

    For i=¯1,M, we denote by ui and φi the Vnn interpolant of uε and φε, respectively. Then ui,φiVnn and

    ui(x,y)=nnl=1uilbl(x,y)i=¯1,M, (4.3)
    φi(x,y)=nnl=1φilbl(x,y)i=¯1,M, (4.4)

    where uil=uε(ti,Nl), φil=φε(ti,Nl), i=¯1,M, l=¯1,nn are the unknowns to be computed.

    Using in addition an implicit (backward) finite difference scheme in time, we introduce now the discrete equations corresponding to (4.1)–(4.2) as follows (see [7,12,13,14,15,16]) for more explanations)

    {Ruil+2Bφil+εhFRuil=B(ρVui1l+2φi1l+εwi1,l1)SφilsξεBuil=B(αξφi1l+εwi1,l2),    i=¯1,M (4.5)

    where uil and φil, l=¯1,nn, are the vectors of unknowns for time level i.

    From the initial conditions (1.3) we have

    u0(x,y)not=u0(x,y)=nnl=1u0(Nl)bl(x,y),φ0(x,y)not=φ0(x,y)=nnl=1φ0(Nl)bl(x,y), (4.6)

    and then from (4.6) we get (see (4.3)–(4.4))

    u0l=u0(Nl)l=¯1,nn,φ0l=φ0(Nl)l=¯1,nn. (4.7)

    The numerical algorithm to compute the approximate solution by fractional steps method can be obtained from the following sequence (again, i denotes the time level)

    Begin algfrac_fem

      i:=0  Compute u0l,φ0l, l=¯1,nn from (4.7);

      Choose wi,l1=w1(ti,Nl), wi,l2=w2(ti,Nl), NlΩ, i=¯0,M1,l=¯1,nn;

      For i:=1 to M do

        Compute zl=z(,Nl), l=¯1,nn from (3.7);

        φi1:=zl, l=¯1,nn;

       Compute uil, φil, l=¯1,nn, solving the linear system (4.5);

      Endfor;

    End.

    The convergence result established by Theorem 2.1 guaranty that the approximate solution computed by the conceptual algorithm algfrac_fem is in fact the approximate solution of the nonlinear parabolic system (1.1)–(1.3).

    The aim of this section is to present an industrial implementation of conceptual algorithm algfrac_fem established in the preview section (in fact an implementation of the numerical model stated by the linear system (4.5) (see P5)).

    From the thermo-kinetics point of view, the solidification and cooling, as well as the simultaneous heating in a continuous casting process of steel, represents a very complex problem of non-stationary heat and mass transfer (see [25]). To solve at present such a problem it is impossible without numerical models of the temperature field and computer technology. In the sequel we will present in short a continuous casting machine, including the key phenomena of interest and a new numerical model that we have used in the settlement of the problem mentioned above.

    The continuous casting process in the metallurgy. In a modern steel casting machine (its essential features are illustrated in Figure 2), the molten metal is tapped from a ladle into a copper mold (crystallizer). Here, the water-cooled walls of the mold (the primary cooling zone) extract heat what leads to solidify a shell that contains the liquid pool. Below the mold (the secondary cooling zone), the product is supported by rollers and is cooled down by water sprays that extract heat from the surface, and, eventually, the core becomes fully solid when the metallurgical length increase at 12÷13m. After the end of secondary cooling zone the product is cooled only by radiation (traiber). Finally, the continuous-casting product must be cut into the optimum lengths (cutting) to achive a maximum yield of metal.

    Figure 2.  Schematic representation of a continuous casting machine.

    The application of the numerical model (4.5) to the continuous casting process, requires experimental research and measurements of operational parameters at MTC2 from Mittal Steel S.A. Galaţi, as well as laboratory research. So, the most important input data in order to do this are (in round bracket we have written the value used by our Matlab program to do the numerical simulations):

    ● the casting speed (V=12.5 mm/s);

    ● physical parameters:

        the density (ρ=7850 kg/m3),

        the latent heat (=65.28kcal/kg),

        the relaxation time (α=1.0e+2),

        the length of separating zone (ξ=.5),

        the coefficients of heat transfer (h=32.012),

        T=44s;

    ● the boundary conditions: (w1(ti,Nl), NlΩ,  i=¯0,M1,l=¯1,nn, in the primary cooling zone (see Figure 4)), (w1(ti,Nl), NlΩ,  i=¯0,M1,l=¯1,nn, in the secondary cooling zone (see Figure 8)) and w2(ti,Nl)=0, NlΩ,  i=¯0,M1,l=¯1,nn;

    Figure 3.  The triangulation Tr over Ω = [0,650]×[0,220].
    Figure 4.  a) the values wi,l1 on the mobile part b) the values wi,l1 on the immobile part – the primary cooling zone.
    Figure 5.  a) the approximate temperature u1, b) the approximate function φ1.
    Figure 6.  The approximate temperature u5.
    Figure 7.  The approximate temperature u20.
    Figure 8.  The values wi,l1 on the fix part - the secondary cooling zone.

    ● dimensions of cristallizer (650 x 1900 x 220), in mm;

    ● the casting temperature (u0=15300C);

    ● the termal conductivity k(u):

    k(u) = [20 100 200 300 400 500 600 700 800 850 900 1000 1100 1200 1600;

        1.43e-5 1.42e-5 1.42e-5 1.42e-5 1.42e-5 9.5e-6 9.5e-6 9.5e-6 8.3e-6...

        8.3e-6 8.3e-6 7.8e-6 7.8e-6 7.4e-6 7.4e-6].

    Numerical experiments. In Figure 3 it can be seen the number of nodes associated to the mesh Tr in the x1 and x2 – axis directions of one half of a rectangular profile. Considering the symmetrical heat removal from the continuous casting (CC) according to the vertical symmetry axis of the rectangular profile, only a half of the cross-section is used in the computation program.

    The numerical model (4.5) uses the temperatures w(t,x,y),t[0,T], (x,y)Ω measured by the thermocouples; the values are illustrated in the Figure 4.

    Figures 57 represents the approximate solution ui, φi (see (4.3), (4.4)), corresponding to different moments of time (i=1, i=5 and i=M).

    A close examination of the Figures 57 tell us the dimension of the solid and liquid zone resulting by runing the Matlab computation program developed on the basis of the conceptual algorithm algfrac_fem.

    Moreover (see [17]), the shape of the graphs shows the stability and accuracy of the numerical results obtained by implementing the fractional steps method.

    The most interesting aspect that we can observe analysing the Figures 67 are the supercooling and superheating phenomenon.

    ● The solidification model that we have considered in this work consist in a nonlinear system of two parabolic differential equations [4]. This new mathematical description of the real phenomenon reflects more accurately the physical aspects, like: superheating, supercooling (see Figures 7 and 11, for example), the effects of surface tension, separating zone of solid and liquid states, etc.

    Figure 9.  The approximate temperature u20, nkx1=10.
    Figure 10.  The approximate temperature u20, nkx1=20.
    Figure 11.  The approximate temperature u40, nkx1=20.

    ● From numerical point of view, the main difficulty in treating the phase field transition system (1.1) is due to the presence of the nonlinear equation corresponding to phase function φ. Thus it is intensely motivated the work in finding more efficient algorithms in order to compute numerically the solution of such system. A scheme of fractional steps type is considered in this sense. This numerical method avoids the iterative process required by the classical methods (e.g., Newton's type approaches) in passing from a time level to another. Numerical tests show that the fractional steps method is faster (CPU-time spent is very small) and the stability and accuracy are higher (see [16,17]) than the Newton's methods.

    ● The distribution of the temperature and the thickness of the solidifying shell, calculated with the numerical model (4.5) obtained following this technique, show that it really is (see Figures 57). New fundamental material properties can also be extracted by analysing the implementation of the numerical model 4.5 (see Figure 7).

    This model is able to simulate the temperature field of a CCM (Continuous Casting Machine) as a whole or any of its parts. In addition, the program elaborated may be used for different slab profiles. The industrial implementation of the numerical model enable us the analysis of the temperature field of the slab when it passes through primary, secondary and traiber zone. The Figures 911 display the calculated isotherms in the secondary cooling zone, while the Figure 12 display the temperature curves on the mobile part, corresponding to the final time level tM.

    Figure 12.  The temperature on the mobile part in the secondary zone: M=40, nkx1=20.

    ● In order to refer the continuous casting process, the following boundary conditions was considered: νu+hu=w1(t,x,y), where h is the heat transfer coefficient and the given function w1(t,x,y) represents the temperature of the surrounding at t[0,T] and (x,y)Ω, ΩR2.

    Generally, the numerical method considered here can be used to approximate the solution of a nonlinear parabolic equation (system) containing a general nonlinear part.

    ● The numerical solution of the phase field transition system of solidification, approximated by this numerical scheme, can be considered as an admissible one for the corresponding boundary optimal control problem (see [16]), formulated in order to study the optimization of the continuous casting process. The numerical results presented in Figures 9-11 illustrate the accuracy of the numerical model (the influence of the density of the net: M = 20, nkx1 = 10, nky1 = 10; M = 20, nkx1 = 20, nky1 = 10; M = 40, nkx1 = 20, nky1 = 10). So, it is strongly motivated to investigate in further the numerical stability of this new numerical model taking into account all parameters (steel casting parameters, physical parameters, net parameters, etc.)

    A detailed discussions on the errors produced by the fractional steps method, illustrate the influence of time and space parameters as well as of all physical parameters (see [17]).

    The author declares that there is no conflicts of interest in this paper.



    [1] O. Axelson and V. Barker, Finite element solution of boundary value problems, Academic Press, 1984.
    [2] T. Benincasa and C. Moroşanu, Fractional steps scheme to approximate the phase-field transition system with non-homogeneous Cauchy-Neumann boundary conditions, Numer. Func. Anal. Opt., 30 (2009), 199-213. doi: 10.1080/01630560902841120
    [3] G. Caginalp, An analysis of a phase field model of a free boundary, Arch. Ration. Mech. An., 92 (1986), 205-245. doi: 10.1007/BF00254827
    [4] G. Caginalp and X. Chen,Convergence of the phase field model to its sharp interface limits, Eur. J. Appl. Math., 9 (1998), 417-445. doi: 10.1017/S0956792598003520
    [5] O. Cârjă, A. Miranville and C. Moroşanu, On the existence, uniqueness and regularity of solutions to the phase-field system with a general regular potential and a general class of nonlinear and non-homogeneous boundary conditions, Nonlinear Anal-Theor, 113 (2015), 190-208.
    [6] G. Iorga, C. Moroşanu and S. C. Cocindău, Numerical simulation of the solid region via phase field transition system, Metal. Int., 13 (2008), 91-95.
    [7] G. Iorga, C. Moroşanu and I. Tofan, Numerical simulation of the thickness accretions in the secondary cooling zone of a continuous casting machine, Metal. Int., 14 (2009), 72-75.
    [8] N. Kenmochi and M. Niezgόdka, Evolution systems of nonlinear variational inequalities arising from phase change problems, Nonlinear Anal-Theor, 22 (1994), 1163-1180. doi: 10.1016/0362-546X(94)90235-6
    [9] A. Miranville and C. Moroşanu, On the existence, uniqueness and regularity of solutions to the phase-field transition system with non-homogeneous Cauchy-Neumann and nonlinear dynamic boundary conditions, Appl. Math. Model., 40 (2016), 192-207.
    [10] A. Miranville and C. Moroşanu, Analysis of an iterative scheme of fractional steps type associated to the nonlinear phase-field equation with non-homogeneous dynamic boundary conditions, Discrete Cont. Dyn-S, 9 (2016), 537-556. doi: 10.3934/dcdss.2016011
    [11] C. Moroşanu, Approximation of the phase-field transition system via fractional steps method, Numer. Func. Anal. Opt., 18 (1997), 623-648. doi: 10.1080/01630569708816782
    [12] C. Moroşanu, et al., Report Stage II/2006, CEEX program no. 84/2005.
    [13] C. Moroşanu, Fractional steps method for approximation the solid region via phase field transition system, 6-th International Conference APLIMAT2007, Bratislava, 6-9 Feb. 2007.
    [14] C. Moroşanu, et al., Report Stage III/2007, CEEX program no. 84/2005.
    [15] C. Moroşanu, Approximation of the solid region in the continuous casting process of steel via phase-field transition system, 6th European Conference on Continuous Casting, Riccione, Italy, 3-6 Jun., 1-6, 2008.
    [16] C. Moroşanu, Analysis and optimal control of phase-field transition system: Fractional steps methods, Bentham Science Publishers, 2012.
    [17] C. Moroşanu, Qualitative and quantitative analysis for a nonlinear reaction-diffusion equation, ROMAI J., 12 (2016), 85-113.
    [18] C. Moroşanu and A. Croitoru, Analysis of an iterative scheme of fractional steps type associated to the phase-field equation endowed with a general nonlinearity and Cauchy-Neumann boundary conditions, J. Math. Anal. Appl., 425 (2015), 1225-1239. doi: 10.1016/j.jmaa.2015.01.033
    [19] C. Moroşanu, I. Crudu, G. Iorga, et al. Research Concerning the Evolution of Solidification Front via Phase-Field Transition System, CEx05-D11-Prog.no., 84 (2008), IFA Bucharest.
    [20] C. Moroşanu and D. Motreanu, A generalized phase field system, J. Math. Anal. Appl., 237 (1999), 515-540. doi: 10.1006/jmaa.1999.6467
    [21] O. A. Oleinik, A method of solution of the general Stefan problem,.Dokl. Akad. Nauk SSSR, 135 (1960), 1354-1357.
    [22] O. Penrose and P. C. Fife, Thermodynamically consistent models of phase-field type for kinetics of phase transitions, International Symposium on Physical Design, 43 (1990), 44-62.
    [23] L. I. Rubinstein, The Stefan problem, Transl. Math. Monographs, 27, American Mathematical Society, Providence, Rhode Island, 1971.
    [24] R. Temam, Infinite-dimensional dynamical systems in mechanics and physics, Vol. 68.of Applied Mathematical Sciences, Springer-Verlag, New York, second edition, 1997.
    [25] B. G. Thomas, Continuous Casting: Modeling, The Encyclopedia of Advanced Materials, (J. Dantzig, A. Greenwell, J. Mickalczyk, eds.), Pergamon Elsevier Science Ltd., Oxford, UK, 2 (2001), 8p.
  • This article has been cited by:

    1. Anca Croitoru, Gabriela Tănase, On a nonlocal and nonlinear second-order anisotropic reaction-diffusion model with in-homogeneous Neumann boundary conditions, 2023, 16, 1937-1632, 75, 10.3934/dcdss.2022155
    2. Ana-Maria Moșneagu, On some local and nonlocal reaction-diffusion models with Robin boundary conditions, 2023, 16, 1937-1632, 104, 10.3934/dcdss.2022161
    3. Cǎtǎlin Stoicescu, On a nonlocal and nonlinear second-order anisotropic reaction-diffusion system with in-homogeneous Cauchy-Neumann boundary conditions. Applications on epidemic infection spread, 2023, 16, 1937-1632, 126, 10.3934/dcdss.2022172
    4. Silviu Dumitru Pavăl, Alex Vasilică, Alin Adochiei, Qualitative and quantitative analysis of a nonlinear second-order anisotropic reaction-diffusion model of an epidemic infection spread, 2023, 16, 1937-1632, 33, 10.3934/dcdss.2022094
    5. Mouhamadou Samsidy Goudiaby, Ben Mansour Dia, Longtime behavior of a semi-implicit scheme for Caginalp phase-field model, 2021, 12, 25900374, 100213, 10.1016/j.rinam.2021.100213
    6. Vasile Berinde, Alain Miranville, Costicǎ Moroşanu, A qualitative analysis of a second-order anisotropic phase-field transition system endowed with a general class of nonlinear dynamic boundary conditions, 2023, 16, 1937-1632, 148, 10.3934/dcdss.2022203
    7. Tudor Barbu, Alain Miranville, Costică Moroşanu, On a Local and Nonlocal Second-Order Boundary Value Problem with In-Homogeneous Cauchy–Neumann Boundary Conditions—Applications in Engineering and Industry, 2024, 12, 2227-7390, 2050, 10.3390/math12132050
    8. Constantin Fetecau , Costică Moroşanu, Dorin-Cătălin Stoicescu, Fractional Steps Scheme to Approximate the Phase Field Transition System Endowed with Inhomogeneous/Homogeneous Cauchy-Neumann/Neumann Boundary Conditions, 2023, 12, 2075-1680, 1098, 10.3390/axioms12121098
    9. Constantin Fetecău, Costică Moroşanu, Fractional Step Scheme to Approximate a Non-Linear Second-Order Reaction–Diffusion Problem with Inhomogeneous Dynamic Boundary Conditions, 2023, 12, 2075-1680, 406, 10.3390/axioms12040406
    10. Gabriela Tănase, A first-order fractional–steps–type method to approximate a nonlinear reaction–diffusion equation with homogeneous Cauchy–Neumann boundary conditions, 2024, 0, 1937-1632, 0, 10.3934/dcdss.2024002
    11. Ana-Maria Moşneagu, Iulian Stoleriu, An adaptive time-step method for a nonlocal version of the Allen-Cahn equation, 2024, 0, 1937-1632, 0, 10.3934/dcdss.2024219
    12. Ana-Maria Moşneagu, Gabriela Tănase, A second-order fractional-steps method to approximate a nonlinear reaction-diffusion equation with non-Homogeneous Robin boundary conditions, 2025, 0, 1937-1632, 0, 10.3934/dcdss.2025028
    13. Mihail-Cristian Heghea, Silviu-Dumitru Pavăl, Vasile-Ion Manta, Reaction-diffusion applications to data dynamic clustering and partitioning, 2020, 0, 1937-1632, 0, 10.3934/dcdss.2025042
  • Reader Comments
  • © 2019 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6082) PDF downloads(1023) Cited by(13)

Figures and Tables

Figures(12)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog