Search

Nonlinear partial differential equations (PDEs) modelling unsteady boundary-layer flows are solved by the spectral relaxation method (SRM) and the spectral quasilinearization method (SQLM). The SRM and SQLM are Chebyshev pseudospectral based methods that have been successfully used to solve nonlinear boundary layer flow problems described by systems of ordinary differential equations. In this paper application of these methods is extended, for the first time, to systems of nonlinear PDEs that model unsteady boundary layer flow. The new extension is tested on two problems: boundary layer flow caused by an impulsively stretching plate and a coupled four-equation system that models the problem of unsteady MHD flow and mass transfer in a porous space. Numerous simulation experiments are conducted to determine the accuracy and compare the computational performance of the proposed methods against the popular Keller-box finite difference scheme which is widely accepted as being one of the ideal tools for solving nonlinear PDEs that model boundary layer flow problems. The results indicate that the methods are more efficient in terms of computational accuracy and speed compared with the Keller-box. [ABSTRACT FROM AUTHOR]

We introduce an adaptive numerical method for computing blow-up solutions for ODEs and well-known reaction-diffusion equations. The method is based on the implicit midpoint method and the implicit Euler method. We demonstrate that the method produces superior results to the adaptive PECE-implicit method and the MATLAB solver of comparable order. [ABSTRACT FROM AUTHOR]

Li et al. (SIAM J. Sci. Comput. 20:719–738, 1998) used the moving mesh partial differential equation (MMPDE) to solve a scaled Fisher's equation and the initial condition consisting of an exponential function. The results obtained are not accurate because MMPDE is based on a familiar arc-length or curvature monitor function. Qiu and Sloan (J. Comput. Phys. 146:726–746, 1998) constructed a suitable monitor function called modified monitor function and used it with the moving mesh differential algebraic equation (MMDAE) method to solve the same problem of scaled Fisher's equation and obtained better results. In this work, we use the forward in time central space (FTCS) scheme and the nonstandard finite difference (NSFD) scheme, and we find that the temporal step size must be very small to obtain accurate results. This causes the computational time to be long if the domain is large. We use two techniques to modify these two schemes either by introducing artificial viscosity or using the approach of Ruxun et al. (Int. J. Numer. Methods Fluids 31:523–533, 1999). These techniques are efficient and give accurate results with a larger temporal step size. We prove that these four methods are consistent for partial differential equations, and we also obtain the region of stability. [ABSTRACT FROM AUTHOR]