COMPACT FINITE-DIFFERENCE SCHEMES FOR WEAKLY NON-LINEAR PROBLEMS AND BOUNDARY CONDITIONS IMITATING CAUCHY PROBLEM

  • V. A. Gordin National Research University «Higher School of Economics» & Hydrometeorological Centre of Russia
DOI 10.29006/1564-2291.JOR-2019.47(1).9
Keywords compact finite-difference scheme, test functions, accuracy order, Richardson extrapolation, discontinuous diffusion coefficient, transparent boundary conditions imitating cauchy problem

Abstract

Compact finite-difference schemes are well known and provide high accuracy order for differential equation with constant coefficients.
Algorithms for constructing compact schemes of the 4-th order for boundary value problems with variable (smooth or jump) coefficient are developed. For the diffusion equations with a smooth variable coefficient and the Levin – Leontovich equation, compact finite-difference schemes are also constructed and their 4-th order is experimentally confirmed. The method of constructing compact schemes of the 4-th order can be generalized to partial differential equations and systems with weak nonlinearity, for example, for the Fisher – Kolmogorov – Petrovsky – Piskunov equation, for the nonlinear Schrödinger equation or for the Fitzhugh – Nagumo system. For such nonlinear problems, a combination of simple explicit schemes and relaxation is used. Richardson’s extrapolation increases the order of the circuits to the 6-th.
To approximate multidimensional problems with discontinuous coefficients, for example, the two-dimensional stationary diffusion equation in inhomogeneous media, it is necessary to estimate the possible asymptotics of solutions in the vicinity of the boundary line’s breaks. To do this, we use generalized eigen-functions in the angle, which can be used as a set of test functions and build compact difference schemes approximating the problem on triangular grids with high order of accuracy. The asymptotics along the radius of these generalized eigen-functions (in polar coordinates in the vicinity of the vertex of the angle) have irrational indices which can be found from a special dispersion equation and which determine the indices of the corresponding Bessel functions along the radius.
For a number of difference schemes approximating the most important evolutionary equations of mathematical physics, it is possible to construct special boundary conditions imitating the Cauchy problem (ICP) on the whole space. These conditions depend not only on the original equation, but also on the type of the difference scheme, and even on the coefficients of the corresponding differential equation. The ICP conditions are determined with accuracy to a gauge. But the choice of this gauge turns out to be essential with numerical implementation. The role of rational approximations of the Pade – Hermite type of the symbol of the corresponding pseudo-differential operator is important. Examples of movie solutions of problems with ICP conditions for various finite-difference schemes approximating the basic mathematical physics equations, see https://cs.hse.ru/mmsg/transbounds.
The study was realized within the framework of the Academic Fund Program at the National Research University – Higher School of Economics (HSE) in 2016–2017 (grant No. 16-05-0069) and by the Russian Academic Excellence Project «5–100».

References


  1. Gordin V.A. O smeshannoj kraevoj zadache, imitiruyushchej zadachu Koshi. Uspekhi matem. nauk.(On Mixed Boundary Problem, Simulated Cauchy Problem. Advances of Mathematical Sciences), 1978, Vol. 33, No. 5, pp. 181–182 (Russian), pp. 189–190 (English).

  2. Gordin V.A. Primenenie proektorov v prognosticheskih skhemah. Trudy Gidrometcentra SSSR (Projectors Using in Forecasting Schemes. Proceedings of the USSR Hydrometeorological Center), 1978, No. 212, pp. 79–96.

  3. Gordin V.A. Granichnoe uslovie polnogo pogloshcheniya voln, vyhodyashchih iz prognosticheskoj oblasti dlya differencial’nogo uravneniya v chastnyh proizvodnyh. Trudy Gidrometcentra SSSR (Boundary Condition of Waves Full Absorption, which Go Away from Prognostic Area for Difference Equation in Partial Derivatives. Proceedings of the USSR Hydrometeorological Center), 1982, No. 242, pp. 104–120.

  4. Gordin V.A. Primenenie vektornoj approksimacii Pade dlya chislennogo resheniya ehvolyucionnyh prognosticheskih uravnenij. (Application of the Pade Vectorial-Approximation to the Numerical Solution of Evolutionary Forecasting Equations). Meteorology and Hydrology, 1982, No. 11, pp. 24–37 (Russian), pp.18–27 (English).

  5. Gordin V.A. Matematicheskie zadachi gidrodinamicheskogo prognoza pogody. Vychislitel’nye aspekty (Mathematical Problems of the Hydrodynamical Weather Forecasting. Numerical Aspects), Leningrad: Gidrometeoizdat, 1987, 264 p.

  6. Gordin V.A. Matematika, komp’yuter, prognoz pogody i drugie scenarii matematicheskoj fiziki (Mathematics, Computer, Weather Forecasting and Other Scenarios of Mathematical Physics), Moskva: Physmatlit, 2010, 2013 (in Russian).

  7. Gordin V.A. and Tsymbalov E.A. Compact difference schemes for the diffusion and Schrodinger equations. Approximation, stability, convergence, effectiveness, monotony. Journal of Computational Mathematics, 2014, Vol. 32, No. 3, pp. 348–370, http://www.global-sci.org/jcm/galley/JCMCR-14.pdf

  8. Gordin V.A. Differencial’nye i raznostnye uravneniya. Kakie yavleniya oni opisyvayut i kak ih reshat (Differential and Finite-Difference Equations. What are Differential and FiniteDifference Equations. What are Phenomena That They Described and How They are Solved), Publisher Moscow: House of Higher School of Economics, 2016.

  9. Gordin V.A. and Tsymbalov E.A. Raznostnaya skhema 4-go poryadka tochnosti dlya differencial’nogo uravneniya s peremennymi koehfficientami (4-th order difference scheme for differential equation with variable coefficients). Mathematical Models and Computer Simulations, 2018, Vol. 10, No. 1, pp. 3–14, (English).

  10. Gordin V.A. and Tsymbalov E.A. Kompaktnaya raznostnaya skhema dlya differencial’nogo uravneniya s kusochno-postoyannym koehfficientom (Compact difference scheme for a differential equation with a piecewise constant coefficient). Matematicheskoe modelirovanie, 2017, Vol. 29, No. 12, pp. 16–28.

  11. Gordin V.A. and Tsymbalov E.A. Compact difference schemes for weakly-nonlinear parabolic and Schrodinger-type equations and systems. 2017, arXiv preprint arXiv, 1712.05185.

  12. Gordin V.A. and Tsymbalov E.A. Compact difference scheme for parabolic and Schrodinger-type equations with variable coefficients. J. Comp. Phys., 2018, Vol. 375, pp. 1451–1468.

  13. Gordin V.A. and Shemendyuk A.A. «Transparent» Boundary Conditions for the Equation of Rod Transverse Vibrations. Preprint https://arxiv.org/submit/2447606/view 2018, submitted to Journal of Sound and Vibration.
Published
2019-05-28
Section
The XXII workshop of the Council of nonlinear dynamics of the Russian Academy of