High-order compact schemes for semilinear parabolic moving boundary problems
Introduction
In recent years, there has been an increasing number of moving boundary problems that occur frequently in physical, chemical, meteorological and biological modeling [1], [8], [17], [22], [24], [26], [28]. For example, the process of tumor growth has been formulated into different mathematical models, which can advance our understanding of cancer research [1], [17], [22]. Besides, the moving boundary problems also arise in the descriptions of the oriented movement of cells [8] and movement of the shoreline in a sedimentary ocean [28] and so on.
In terms of analysis of moving-boundary problems, Caboussat and Rappaz study the well-posedness of a one-dimensional free surface problem and obtain the existence and uniqueness results locally in time [5]. Moreover, Muntean and Böhm deal with a one-dimensional coupled system of semi-linear parabolic equations on the moving boundary and prove a global existence and uniqueness theorem of positive weak solutions [20]. Escher and Matioc re-express the problem of the growth of nonnecrotic tumors as an abstract evolution equation and prove local well-posedness in the small Hölder spaces context [10]. Besides, Sun et al. consider the global existence and uniqueness of a transient solution of a non-autonomous free boundary model for tumor growth [25]. However, as we all know, it is extraordinarily challenging to obtain every analytical solution of a moving boundary problem on account of the complexity of boundaries. Therefore, establishing effective and high-accuracy numerical methods is a must.
To the authors' best knowledge, some numerical methods have been described for the solution of time-dependent partial differential equations with moving boundaries, which can be classified in two ways. On the one hand, based on the moving mesh methods, there are many different numerical algorithms proposed by researchers, which can adapt to different problems with specific features [2], [3], [4], [7], [12], [14]. For example, Baines et al. [2], [3], [4] study moving finite element methods and propose moving mesh finite element algorithm for the adaptive solution of nonlinear diffusion equations with moving boundaries in one and two dimensions. Lee et al. [14] construct a velocity-based method for modeling avascular tumor growth, whose estimated order of convergence is mostly less than or equal to two. Huang et al. establish a moving-mesh method based on the geometric conservation law in [7] and study the stability of three moving-mesh finite-difference schemes in the norm for one-dimensional linear convection-diffusion equations in [12].
On the other hand, an alternative approach is to rescale the original problem with moving boundary to an equivalent problem on a certain fixed domain [6], [11], [13], [17], [18], [21], [25], [33], [38]. Kovács and Lubich study spatial and temporal discretizations of parabolic differential equations with dynamic boundary conditions in a weak formulation [13]. Besides, numerous finite difference schemes have been applied to solve moving boundary problems numerically [11], [17], [18], [21], [25]. Unfortunately, there is no numerical analysis of the convergence and stability of difference schemes in the above literature. Cao and Sun [6] establish a Crank-Nicolson scheme with the unconditional stability for the linear convection-diffusion equation with moving boundaries and prove that the scheme proposed is convergent in the maximum norm and the order of convergence is . Yuan et al. study a three-dimensional moving boundary problem on the compressible miscible (oil and water) displacement by a second-order upwind difference fractional steps scheme applicable to parallel computing [33]. Zheng et al. construct an efficient numerical method to solve a fractional reaction-diffusion model with a moving boundary and obtain the convergence and stability [38].
In terms of the compact difference methods widely used to solve partial differential equations numerically, there are a great number of works in the literature on the numerical solution of partial differential equations. The readers interested in this aspect are referred [9], [15], [16], [19], [29], [30], [31], [32], [34], [35], [36], [37] for details. However, to the authors' best knowledge, there are few works in the literature on the compact schemes with the higher-order convergence for partial differential equations with moving boundary conditions.
In reality, moving-boundary problems include both problems with unknown boundary, which are often called Stefan problems, and prescribed-boundary problems. In this work, we focus on the semilinear convection-diffusion equation with prescribed-boundary condition. where κ and γ are both positive constants, , and the boundary value condition can be rewritten as
(FP) can be regarded as an extension of the linear prescribed-boundary problem in [6]. Here, we assume that the initial value in (FP) is regular enough and the exterior force g has the second-order continuous derivative with respect to the first components in the -neighborhood of the solution v, where is a positive constant.
Besides, we assume that the boundary value Φ is piecewise smooth and the following consistency conditions hold: , i.e., . We suppose that and , for any . And the domain can be defined as , which is an interval in R for each .
In this paper, we develop a kind of compact scheme with high accuracy and high convergence rate to solve (FP) numerically and carry out the corresponding analysis. Firstly, using a linear transformation and an exponential transformation, we can establish a compact scheme for the equation (FP). Secondly, the unique solvability and the convergence are established. Thirdly, we extend the methodology to establish both a Crank-Nicolson scheme and a compact scheme for the two-dimensional moving boundary problem. Finally, we implement some numerical examples to support our theoretical results.
The content of this paper will be organized as follows. In Section 2, we obtain an equivalent system defined on a rectangular region. In Section 3, we establish a compact scheme with fourth-order accuracy in the spatial dimension and second-order accuracy in the temporal dimension, see Theorem 4.5. Besides, unique solvability and convergence of the established scheme are studied in Section 4. Furthermore, a Crank-Nicolson scheme and a compact scheme for a two-dimensional semilinear moving boundary problem are given as well in Section 5. Finally, some results of numerical examples are given in Section 6.
Section snippets
Problem reformation by a linear variable transformation and an exponential transformation
First and foremost, we introduce a linear transformation Denote and We can obtain where .
Furthermore, the exponential transformation [15] is introduced as
Derivation of a compact scheme
Before the presentation of the compact scheme, the uniform-grid discretization is used. Let be a uniform mesh on the interval with , where m is a positive integer and . Let be a uniform mesh on the interval with , where n is a positive integer and . Therefore, the rectangular region is subdivided by . Denote . Let
Convergence of the compact scheme
In this section, we analyze the solvability and convergence of the compact scheme (CS). Denote and introduce the following notations for :
The following lemmas are necessary for our main results. Lemma 4.1 [23], [27] If is a grid function defined on and satisfy , then Proof For the last inequality, we
Extension to two-dimensional semilinear moving boundary problem
Considering the two-dimensional semilinear moving boundary problem as along with the initial condition and moving boundary condition where are positive constants. , and can be defined as two domains, where and ,
Numerical examples
To verify our results in the previous sections, we implement four numerical examples, including one-dimensional linear problem (Case 1 and Case 2), one-dimensional semilinear problem (Case 3) and two-dimensional problem (Case 4). Although the compact scheme (CS) is established for semilinear parabolic moving boundary problems, it should be used to solve linear problems as well. That is to say that the compact scheme (CS) should apply equally when the equation (FP) reduces to a linear problem.
Conclusion
In summary, we establish the compact scheme for one-dimensional semilinear parabolic moving boundary problems and study the convergence. The proposed method has the convergence order of two in the temporal direction and the convergence order of four in the spatial direction in -norm, respectively. Furthermore, we extend the methodology to two-dimensional semilinear problems as well in this paper. With the help of the exponential transformation, we can avoid discussing a compact scheme on
Acknowledgement
This work was supported by the National Natural Science Foundation of China (Grant No. 11871435, 11471287, 11501514), the Natural Sciences Foundation of Zhejiang Province (Grant No. LY19A010026), Zhejiang Province “Yucai” Project (2019), Fundamental Research Funds of Zhejiang Sci-Tech University (Grant No. 2019Q072) and Zhejiang Provincial Department of Education Research Project (Grant No. Y201942553).
References (38)
- et al.
A history of the study of solid tumour growth: the contribution of mathematical modelling
Bull. Math. Biol.
(2004) - et al.
An stability analysis for finite difference solution of one-dimensional convection-diffusion problems on moving meshes
J. Comput. Appl. Math.
(2012) - et al.
A moving mesh approach for modelling avascular tumour growth
Appl. Numer. Math.
(2013) - et al.
Numerical simulation of blowup in nonlocal reaction-diffusion equations using a moving mesh method
J. Comput. Appl. Math.
(2009) - et al.
A compact-difference scheme for the Navier-Stockes equations in vorticity-velocity formulation
J. Comput. Phys.
(2000) - et al.
A moving-boundary problem for concrete carbonation: global existence and uniqueness of weak solutions
J. Math. Anal. Appl.
(2009) - et al.
Numerical solution for the model of vibrating elastic membrane with moving boundary
Commun. Nonlinear Sci.
(2007) - et al.
Analysis of two-dimensional freezing on the outside of a coolant carrying tube
Int. J. Heat Mass Transf.
(1981) - et al.
High-order compact exponential finite difference methods for convection-diffusion type problems
J. Comput. Phys.
(2007) - et al.
A high-order exponential ADI scheme for two dimensional time fractional convection-diffusion equations
Comput. Math. Appl.
(2014)
The compact and Crank-Nicolson ADI schemes for two-dimensional semilinear multidelay parabolic equations
J. Comput. Appl. Math.
Numerical solution of the time fractional reaction diffusion equation with a moving boundary
J. Comput. Phys.
Moving Finite Elements
A moving mesh finite element algorithm for the adaptive solution of time-dependent partial differential equations with moving boundaries
Appl. Numer. Math.
Velocity-based moving mesh methods for nonlinear partial differential equations
Commun. Comput. Phys.
Analysis of a one-dimensional free boundary flow problem
Numer. Math.
Maximum norm error estimates of the Crank-Nicolson scheme for solving a linear moving boundary problem
J. Comput. Appl. Math.
A moving-mesh method based on the geometric conservation law
SIAM J. Sci. Comput.
A free boundary problem for a class of parabolic type chemotaxis model
Kinet. Relat. Models
Cited by (1)
Stability estimate for a nonlinear coupled heat transfer model
2023, Applied Mathematics Letters