Extension of high-order compact schemes to time-dependent problems (Q2762702)

The authors consider the model transient convection diffusion equation NEWLINE\[NEWLINE\frac{\partial\varphi}{\partial t}+u\frac{\partial\varphi}{\partial x}=a\;\frac{\partial^2\varphi}{\partial x^2}+S,\tag{1}NEWLINE\]NEWLINE where \(a\) is a strictly positive, \(u\) (for simplicity of exposition) is constant and \(S(x,t)\) is a smooth source term. For the associated steady state equation NEWLINE\[NEWLINE-a\;\frac{\partial^2\varphi}{\partial x^2}+u\;\frac{\partial\varphi}{\partial x}=S\tag{2}NEWLINE\]NEWLINE a higher-order compact scheme is constructed as follows: applying central differences to (2) and then repeatedly differentiate (2) to replace higher-order derivatives in the truncation error by lower-order derivatives that can be differenced gives a fourth-order compact difference scheme.NEWLINENEWLINENEWLINEThis approach is applied to the problem (1) replacing \(S\) by \(S-\frac{\partial\varphi}{\partial t}\). To the obtained semi-discrete system of ordinary differential equations the authors apply a two-level integration scheme and in terms of Courant and Reynolds (Pedet) numbers consider the von Neumann stability of the difference scheme. The same approach is applied for diffusion in \(2D\).NEWLINENEWLINENEWLINENext the implementation of the von Neumann boundary conditions is investigated. (The case of Dirichlet boundary conditions is not discussed because of its simplicity). A first-order system is also considered and several numerical examples are given. The bibliography contains thirty entries.