×

Numerical methods for nonlinear second-order hyperbolic partial differential equations. II: Rothe’s techniques for 1-D problems. (English) Zbl 1141.65380

Summary: Two families of Rothe’s methods that are based on the discretization of the time variable and keeping the spatial one continuous for the solution of second-order hyperbolic equations with damping and nonlinear source terms are presented. The first family is based on time integration and results in a Volterra integro-differential equation which, upon approximating certain time integrals, can be written as an ordinary differential equation in space. Upon discretizing the spatial derivatives by means of finite difference formula, this family results in nonlinear algebraic equations which contain exponential terms that depend on the time step, and its accuracy is shown to be lower than that of an explicit second-order accurate discretization in both space and time of second-order hyperbolic equations.
The second family of methods is based on the discretization of the time derivatives and time linearization, and results in linear ordinary differential equations in space. Upon freezing the coefficients of these equations, one can integrate analytically the resulting linear ordinary differential equations to obtain piecewise exponential solutions which are continuous in space, and three-point finite difference equations which depend exponentially on the time step, spatial grid size, and the diffusion, relaxation, damping and reaction times. The finite difference equations are shown to result in non-diagonally dominant matrices unless the time step is smaller than the characteristic relaxation, diffusion, damping and reaction times.
To avoid this problem, two Rothe’s techniques that do not account for the Jacobian of the reaction terms in the ordinary differential operator are developed, and it is shown that the Rothe’s techniques belonging to the second family are as accurate as the linearly implicit methods presented in Part I [ibid. 190, No. 1, 722–756 (2007; Zbl 1123.65087)] provided that the relaxation time is smaller than the critical time and the solution does not oscillate in space. When the relaxation time is greater than the critical one, Rothe’s methods have been found to be less accurate than the finite difference techniques of Part I [loc. cit.]. The Rothe’s techniques belonging to the second family are generalized to systems of nonlinear second-order hyperbolic equations and mixed systems of parabolic and second-order hyperbolic partial differential equations.

MSC:

65M20 Method of lines for initial value and initial-boundary value problems involving PDEs
35L70 Second-order nonlinear hyperbolic equations
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs

Citations:

Zbl 1123.65087
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] J.I. Ramos, Numerical methods for nonlinear second-order hyperbolic partial differential equations. I - Time-linearized finite difference methods for 1-D problems, Appl. Math. Comput., in press.; J.I. Ramos, Numerical methods for nonlinear second-order hyperbolic partial differential equations. I - Time-linearized finite difference methods for 1-D problems, Appl. Math. Comput., in press. · Zbl 1123.65087
[2] Rothe, E., Zweidimensionale parabolische Randwertaufgaben als eindimensionaler Randwertaufgaben, Math. Ann., 102, 652-670 (1930) · JFM 56.1076.02
[3] Ashyralyev, A.; Sobolevskii, P. E., Well-Posedness of Parabolic Difference Equations (1994), Birkhäuser Verlag: Birkhäuser Verlag Basel, Boston, Berlin · Zbl 1077.39015
[4] Martı´nez, V.; Marquina, A.; Donat, R., Shooting methods for one-dimensional diffusion-absorption problems, SIAM J. Numer. Anal., 31, 572-589 (1994) · Zbl 0821.65062
[5] Gilbarg, D.; Trudinger, N. S., Elliptic Partial Differential Equations of Second Order (1983), Springer-Verlag: Springer-Verlag Berlin · Zbl 0691.35001
[6] Ashyralyev, A.; Sobolevskii, P. E., New Difference Schemes for Partial Differential Equations (2004), Birkhäuser Verlag: Birkhäuser Verlag Basel, Boston, Berlin · Zbl 1354.65212
[7] Strikwerda, J. C., Finite Difference Schemes and Partial Differential Equations (1989), Wadsworth & Brooks: Wadsworth & Brooks Boston · Zbl 0681.65064
[8] Kikuchi, N.; Kac˘ur, J., Convergence of Rothe’s method in Hölder spaces, Appl. Math., 48, 353-365 (2003) · Zbl 1099.65079
[9] Cimrák, I., Comparison of Broyden and Newton methods for solving nonlinear parabolic equations, J. Electr. Eng., 52, 48-52 (2001) · Zbl 1057.35003
[10] Lang, J., High-resolution self-adaptive computations on chemical reaction-diffusion problems with internal boundaries, Chem. Eng. Sci., 51, 1055-1070 (1996)
[11] Kadalbajoo, M. K.; Sharma, K. K.; Awasthi, A., A parameter-uniform implicit difference scheme for solving time-dependent Burgers’s equations, Appl. Math. Comput., 170, 1365-1393 (2005) · Zbl 1084.65081
[12] Aksan, E. N.; Özdeş, A., A numerical solution of Burgers’s equation, Appl. Math. Comput., 156, 395-402 (2004) · Zbl 1061.65085
[13] Aksan, E. N., A numerical solution of Burgers’s equation by finite element method constructed on the method of discretization of time, Appl. Math. Comput., 170, 895-904 (2005) · Zbl 1103.65104
[14] Abbasbandy, S.; Darvishi, M. T., A numerical solution of Burgers’s equation by modified Adomian method, Appl. Math. Comput., 163, 1265-1272 (2005) · Zbl 1060.65649
[15] Schmidt, F.; Yevick, D., Discrete transparent boundary conditions for Schrödinger-type equations, J. Comput. Phys., 134, 96-107 (1997) · Zbl 0878.65111
[16] Tokar, M. Z.; Kalupin, D.; Pilipenko, D., Numerical solution of transport equations for plasmas with transport barriers, Comput. Phys. Commun., 175, 30-35 (2006) · Zbl 1196.76062
[17] Ramos, J. I., On the numerical treatment of an ordinary differential equation arising in one-dimensional non-Fickian diffusion problems, Comput. Phys. Commun., 170, 231-238 (2005) · Zbl 1196.65163
[18] Mueller, J. L.; Shores, T. S., A new sinc-Galerkin method for convection-diffusion equations with mixed boundary conditions, Comput. Math. Appl., 47, 803-822 (2004) · Zbl 1058.65083
[19] Grossmann, C.; Noack, A., Smoothing and Rothe’s method for Stefan problems in enthalpy form, J. Comput. Appl. Math., 138, 347-366 (2002) · Zbl 0994.65099
[20] Samec, N.; Škerget, Integral formulation of a diffusive-convective transport equation for reacting flows, Eng. Anal. Bound. Elem., 28, 1055-1060 (2004) · Zbl 1130.76382
[21] John, V.; Matthies, G.; Rang, J., A comparison of time-discretization/linearization approaches for the incompressible Navier-Stokes equations, Comput. Meth. Appl. Mech. Eng., 195, 5995-6010 (2006) · Zbl 1124.76041
[22] Sheu, T. W.H.; Wang, S. K.; Lin, R. K., An implicit scheme for solving the convection-diffusion-reaction equation in two dimensions, J. Comput. Phys., 164, 123-142 (2000) · Zbl 0995.76065
[23] Sheu, T. W.H.; Lin, R. K., Development of a convection-diffusion-reaction magnetohydrodynamic solver on non-staggered grids, Int. J. Numer. Meth. Eng., 45, 1209-1233 (2004) · Zbl 1060.76619
[24] Ramos, J. I., On diffusive methods and exponentially fitted techniques, Appl. Math. Comput., 103, 69-96 (1999) · Zbl 0929.65057
[25] Chen, G. Q.; Gao, Z.; Yang, Z. F., A perturbational \(h^4\) exponential finite difference scheme for the convective diffusion equation, J. Comput. Phys., 104, 129-139 (1993) · Zbl 0774.65073
[26] Tian, Z. F.; Ye, Y. B., A fourth-order compact ADI method for solving two-dimensional unsteady convection-diffusion problems, J. Comput. Appl. Math., 198, 268-286 (2007) · Zbl 1104.65086
[27] Campo, A.; Amon, C. H., Remarkable improvement of the Lévêque solution for isoflux heating with a combination of the transversal method of lines (TMOL) and a computer-extended Fröbenius power series, Int. J. Heat Mass Transfer, 48, 2110-2116 (2005) · Zbl 1189.76384
[28] Ingber, M. C.; Tanskin, J. A., A parallel domain decomposition BEM for diffusion problems with surface reactions, Eng. Anal. Bound. Elem., 29, 878-885 (2005) · Zbl 1182.76911
[29] Garbey, M., On some applications of the superposition principle with Fourier basis, ISIAM J. Sci. Comput., 22, 1087-1116 (2000) · Zbl 0980.65109
[30] Schiesser, W. E., The Numerical Method of Lines Integration of Partial Differential Equations (1991), Academic Press: Academic Press San Diego · Zbl 0763.65076
[31] Sinha, R. K., Finite element approximations with quadrature for second-order hyperbolic equations, Numer. Meth. PDEs, 18, 537-559 (2002) · Zbl 1011.65068
[32] Liskovets, O. A., The method of lines (review), Differen. Equat., 1, 1308-1323 (1965) · Zbl 0229.65065
[33] Baker, G. A.; Bramble, J. H., Semidiscrete and single step fully discrete approximations for second order hyperbolic equations, RAIRO Anal. Numér./Numer. Anal., 13, 75-100 (1979) · Zbl 0405.65057
[34] van der Houwen, P. J., Modified Nyström methods for semi-discrete hyperbolic differential equations, SIAM J. Numer. Anal., 18, 1081-1097 (1981) · Zbl 0498.65039
[35] Alonso-Mallo, I.; Cano, B.; Moreta, M. J., Order reduction and how to avoid it when explicit Runge-Kutta-Nyström methods are used to solve linear partial differential equations, J. Comput. Appl. Math., 176, 293-318 (2005) · Zbl 1069.65101
[36] Chen, L.-Q.; Zhao, W.-J.; Zu, J. W., Simulations of transverse vibrations of an axially moving string: a modified difference approach, Appl. Math. Comput., 166, 596-607 (2005) · Zbl 1329.74301
[37] L.W. Nagel, SPICE: A computer program to simulate semiconductor circuits, Memo No. CB/ERL M520, Electronic Research Laboratory, University of California, Berkeley, California, 1975.; L.W. Nagel, SPICE: A computer program to simulate semiconductor circuits, Memo No. CB/ERL M520, Electronic Research Laboratory, University of California, Berkeley, California, 1975.
[38] PSPICE, version 6.0, Microsim Corporation, Irvine, California, 1994.; PSPICE, version 6.0, Microsim Corporation, Irvine, California, 1994.
[39] Ramos, J. I., Equivalence of \(C^0\) and \(C^1\) methods for ODEs, Appl. Math. Comput., 167, 1331-1338 (2005) · Zbl 1082.65552
[40] Ramos, J. I., Exponential methods for one-dimensional reaction-diffusion equations, Appl. Math. Comput., 170, 380-398 (2005) · Zbl 1082.65572
[41] Ramos, J. I., A piecewise-analytical method for singularly perturbed parabolic problems, Appl. Math. Comput., 161, 501-512 (2005) · Zbl 1061.65088
[42] Ramos, J. I., Linearly-implicit, approximate factorization, exponential methods for multi-dimensional reaction-diffusion equations, Appl. Math. Comput., 174, 1609-1633 (2006) · Zbl 1090.65109
[43] Malachowski, M. J., Transition probability coefficients and the spatial stability of finite difference schemes for the diffusion and Telegrapher’s equations, Int. J. Numer. Model., 15, 243-249 (2002) · Zbl 0996.65097
[44] Tzou, D. Y., Macro- to Microscale Heat Transfer (1997), Taylor & Francis: Taylor & Francis Washington, DC
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.