Previous Articles    

ANALYSIS ON A NUMERICAL SCHEME WITH SECOND-ORDER TIME ACCURACY FOR NONLINEAR DIFFUSION EQUATIONS

Xia Cui1, Guangwei Yuan1, Fei Zhao2,3   

  1. 1. Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, P. O. Box 8009-26, Beijing 100088, China;
    2. College of Science, North China University of Technology, Beijing 100144, China;
    3. Graduate School of China Academy of Engineering Physics, Beijing 100088, China
  • Received:2020-03-12 Revised:2020-05-27 Online:2021-09-15 Published:2021-10-15
  • Supported by:
    This work is supported by the National Natural Science Foundation of China (11871112, 11971069, 11971071, U1630249), Yu Min Foundation and the Foundation of LCP. The authors would like to thank the anonymous referees for their helpful suggestions to enhance the paper.

Xia Cui, Guangwei Yuan, Fei Zhao. ANALYSIS ON A NUMERICAL SCHEME WITH SECOND-ORDER TIME ACCURACY FOR NONLINEAR DIFFUSION EQUATIONS[J]. Journal of Computational Mathematics, 2021, 39(5): 777-800.

A nonlinear fully implicit finite difference scheme with second-order time evolution for nonlinear diffusion problem is studied. The scheme is constructed with two-layer coupled discretization (TLCD) at each time step. It does not stir numerical oscillation, while permits large time step length, and produces more accurate numerical solutions than the other two well-known second-order time evolution nonlinear schemes, the Crank-Nicolson (CN) scheme and the backward difference formula second-order (BDF2) scheme. By developing a new reasoning technique, we overcome the difficulties caused by the coupled nonlinear discrete diffusion operators at different time layers, and prove rigorously the TLCD scheme is uniquely solvable, unconditionally stable, and has second-order convergence in both space and time. Numerical tests verify the theoretical results, and illustrate its superiority over the CN and BDF2 schemes.

CLC Number: 

[1] C. Cancès and C. Guichard, Numerical analysis of a robust free energy diminishing finite volume scheme for parabolic equations with gradient structure, Found. Comput. Math., 17:6(2016), 1525-1584.
[2] Z.X. Chen and R.E. Ewing, Fully discrete finite element analysis of multiphase flow in groundwater hydrology, SIAM J. Numer. Anal., 34:6(1997), 2228-2253.
[3] X. Cui, G.W. Yuan and J.Y. Yue, Numerical analysis and iteration acceleration of a fully implicit scheme for nonlinear diffusion problem with second-order time evolution, Numer. Meth. Part. D. E., 32(2016), 121-140.
[4] X. Cui and J Y. Yue, Property analysis and quick solutions for nonlinear discrete schemes for conservative diffusion equation, Math. Numer. Sin., 37:3(2015), 227-246.
[5] W.L. Dai and P.R. Woodward, Numerical simulations for nonlinear heat transfer in a system of multimaterials, J. Comput. Phys., 139:1(1998), 58-78.
[6] W.W. Dai and A.J. Scannapieco, Interface- and discontinuity-aware numerical schemes for plasma 3-T radiation diffusion in two and three dimensions, J. Comput. Phys., 300(2015), 643-664.
[7] W.W. Dai and A.J. Scannapieco, Second-order accurate interface- and discontinuity-aware diffusion solvers in two and three dimensions, J. Comput. Phys., 281(2015), 982-1002.
[8] J. Douglas JR. and T. Dupont, Galerkin methods for parabolic equations, SlAM J. Numer. Anal, 7:4(1970), 575-626.
[9] J. Droniou, R. Eymard, T. Gallouët, C. Guichard and R. Herbin, The Gradient Discretisation Method, Math matiques et Applications 82, Springer International Publishing AG, part of Springer Nature, Cham, Switzerland, 2018.
[10] J. Droniou, R. Eymard and R. Herbin, Gradient schemes:generic tools for the numerical analysis of diffusion equations, ESAIM:Math. Model. Numer. Anal., 50:3(2016), 749-781.
[11] R. Eymard, T. Gallouët and R. Herbin, Finite Volume Methods. In Techniques of Scientific Computing, Part Ⅲ, Handbook of Numerical Analysis, VⅡ, P. G. Ciarlet and J. L. Lions, (ed.), North-Holland, Amsterdam, (2000), 713-1020.
[12] D.A. Knoll, W.J. Rider and G.L. Olson, Nonlinear convergence, accuracy, and time step control in non-equilibrium radiation diffusion, J. Quant. Spectrosc. Radiat. Transfer, 70(2001), 25-36.
[13] Y.X. Liu and C.W. Shu, Analysis of the local discontinuous Galerkin method for the drift-diffusion model of semiconductor devices, Sci. China Math., 59:1(2016), 115-140.
[14] B.Y. Li, H.D. Gao and W.W. Sun, Unconditionally optimal error estimates of a Crank-Nicolson Galerkin method for the nonlinear thermistor equations, SIAM J. Numer. Anal., 52:2(2014), 933-954.
[15] B.Y. Li and W.W. Sun, Error analysis of linearized semi-implicit Galerkin finite element methods for nonlinear parabolic equations, Int. J. Numer. Anal. Model., 10:3(2013), 622-633.
[16] B.Y. Li and W.W. Sun, Unconditional convergence and optimal error estimates of a Galerkinmixed FEM for incompressible miscible flow in porous media, SIAM J. Numer. Anal., 51:4(2013), 1959-1977.
[17] D.F. Li and J.L. Wang, Unconditionally optimal error analysis of Crank-Nicolson Galerkin FEMs for a strongly nonlinear parabolic system, J. Sci. Comput., 72:2(2017), 892-915.
[18] C.C. Ober and J.N. Shadid, Studies on the accuracy of time-integration methods for the radiationdiffusion equations, J. Comput. Phys., 195(2004), 743-772.
[19] G.L. Olson, Efficient solution of multi-dimensional flux-limited nonequilibrium radiation diffusion coupled to material conduction with second-order time discretization, J. Comput. Phys., 226(2007), 1181-1195.
[20] R.M. Rauenzahn, V.A. Mousseau and D.A. Knoll, Temporal accuracy of the nonequilibrium radiation diffusion equations employing a Saha ionization model, Comput. Phys. Commun., 172(2005), 109-118.
[21] D.Y. Shi, F.N. Yan and J.J. Wang, Unconditionally superclose analysis of a new mixed finite element method for nonlinear parabolic equations, J. Comput. Math., 37:1(2019), 1-17.
[22] W.W. Sun and Z.Z. Sun, Finite difference methods for a nonlinear and strongly coupled heat and moisture transport system in textile materials, Numer. Math., 120(2012), 153-187.
[23] M.F. Wheeler, A priori L2 error estimates for Galerkin approximations to parabolic partial differential equations, SIAM J. Numer. Anal., 10:4(1973), 723-759.
[24] G.W. Yuan, Uniqueness and stability of difference solution with nonuniform meshes for nonlinear parabolic systems, Math. Numer. Sin., 22:2(2000), 139-150.
[25] G.W. Yuan and X.D. Hang, Acceleration methods of nonlinear iteration for nonlinear parabolic equations, J. Comput. Math., 24:3(2006), 412-424.
[26] G.W. Yuan, X.D. Hang, Z.Q. Sheng and J.Y. Yue, Progress in numerical methods for radiation diffusion equations, Chinese J. Comput. Phys., 26:4(2009), 475-500.
[27] J.Y. Yue and G.W. Yuan, Picard-Newton iterative method with time step control for multimaterial non-equilibrium radiation diffusion problem, Commun. Comput. Phys., 10:4(2011), 844-866.
[28] G.D. Zhang, J.J. Yang and C.J. Bi, Second order unconditionally convergent and energy stable linearized scheme for MHD equations, Adv. Comput. Math., 44(2018), 505-540.
[29] F. Zhao, X. Cui and G.W. Yuan, Iterative acceleration methods with second-order time accuracy for nonlinear diffusion equations, Adv. Comput. Math., 46:7(2020), https://doi.org/10.1007/s10444-020-09756-4.
[30] B.Y. Zhou and D.F. Li, Newton linearized methods for semilinear parabolic equations, Numer. Math. Theor. Meth. Appl., 13(2020), 1-18.
[31] Y.L. Zhou, Applications of Discrete Functional Analysis to the Finite Difference Method, International Academic Publishers, Beijing, 1990.
[32] Y.L. Zhou, Difference schemes with nonuniform meshes for nonlinear parabolic system, J. Comput. Math., 14:4(1996), 319-335.
[1] Baiying Dong, Xiufeng Feng, Zhilin Li. AN L SECOND ORDER CARTESIAN METHOD FOR 3D ANISOTROPIC INTERFACE PROBLEMS [J]. Journal of Computational Mathematics, 2022, 40(6): 882-912.
[2] Darko Volkov. A STOCHASTIC ALGORITHM FOR FAULT INVERSE PROBLEMS IN ELASTIC HALF SPACE WITH PROOF OF CONVERGENCE [J]. Journal of Computational Mathematics, 2022, 40(6): 955-976.
[3] Yanping Chen, Qiling Gu, Qingfeng Li, Yunqing Huang. A TWO-GRID FINITE ELEMENT APPROXIMATION FOR NONLINEAR TIME FRACTIONAL TWO-TERM MIXED SUB-DIFFUSION AND DIFFUSION WAVE EQUATIONS [J]. Journal of Computational Mathematics, 2022, 40(6): 936-954.
[4] Mingming Zhao, Yongfeng Li, Zaiwen Wen. A STOCHASTIC TRUST-REGION FRAMEWORK FOR POLICY OPTIMIZATION [J]. Journal of Computational Mathematics, 2022, 40(6): 1004-1030.
[5] Rong Zhang, Hongqi Yang. A DISCRETIZING LEVENBERG-MARQUARDT SCHEME FOR SOLVING NONLIEAR ILL-POSED INTEGRAL EQUATIONS [J]. Journal of Computational Mathematics, 2022, 40(5): 686-710.
[6] Wei Zhang. STRONG CONVERGENCE OF THE EULER-MARUYAMA METHOD FOR A CLASS OF STOCHASTIC VOLTERRA INTEGRAL EQUATIONS [J]. Journal of Computational Mathematics, 2022, 40(4): 607-623.
[7] Yaolin Jiang, Zhen Miao, Yi Lu. WAVEFORM RELAXATION METHODS FOR LIE-GROUP EQUATIONS* [J]. Journal of Computational Mathematics, 2022, 40(4): 649-666.
[8] Xiaonian Long, Qianqian Ding. A SECOND ORDER UNCONDITIONALLY CONVERGENT FINITE ELEMENT METHOD FOR THE THERMAL EQUATION WITH JOULE HEATING PROBLEM [J]. Journal of Computational Mathematics, 2022, 40(3): 354-372.
[9] Siyuan Qi, Guangqiang Lan. STRONG CONVERGENCE OF THE EULER-MARUYAMA METHOD FOR NONLINEAR STOCHASTIC VOLTERRA INTEGRAL EQUATIONS WITH TIME-DEPENDENT DELAY [J]. Journal of Computational Mathematics, 2022, 40(3): 437-452.
[10] Xiaoyu Wang, Ya-xiang Yuan. STOCHASTIC TRUST-REGION METHODS WITH TRUST-REGION RADIUS DEPENDING ON PROBABILISTIC MODELS [J]. Journal of Computational Mathematics, 2022, 40(2): 294-334.
[11] Mohammed Harunor Rashid. METRICALLY REGULAR MAPPING AND ITS UTILIZATION TO CONVERGENCE ANALYSIS OF A RESTRICTED INEXACT NEWTON-TYPE METHOD [J]. Journal of Computational Mathematics, 2022, 40(1): 44-69.
[12] Yang Chen, Chunlin Wu. DATA-DRIVEN TIGHT FRAME CONSTRUCTION FOR IMPULSIVE NOISE REMOVAL [J]. Journal of Computational Mathematics, 2022, 40(1): 89-107.
[13] Qianqian Chu, Guanghui Jin, Jihong Shen, Yuanfeng Jin. NUMERICAL ANALYSIS OF CRANK-NICOLSON SCHEME FOR THE ALLEN-CAHN EQUATION [J]. Journal of Computational Mathematics, 2021, 39(5): 655-665.
[14] Lu Zhang, Qifeng Zhang, Hai-wei Sun. A FAST COMPACT DIFFERENCE METHOD FOR TWO-DIMENSIONAL NONLINEAR SPACE-FRACTIONAL COMPLEX GINZBURG-LANDAU EQUATIONS [J]. Journal of Computational Mathematics, 2021, 39(5): 708-732.
[15] Yong Liu, Chi-Wang Shu, Mengping Zhang. SUB-OPTIMAL CONVERGENCE OF DISCONTINUOUS GALERKIN METHODS WITH CENTRAL FLUXES FOR LINEAR HYPERBOLIC EQUATIONS WITH EVEN DEGREE POLYNOMIAL APPROXIMATIONS [J]. Journal of Computational Mathematics, 2021, 39(4): 518-537.
Viewed
Full text


Abstract