New Stable, Explicit, Shifted-Hopscotch Algorithms for the Heat Equation

Abstract

Our goal was to find more effective numerical algorithms to solve the heat or diffusion equation. We created new five-stage algorithms by shifting the time of the odd cells in the well-known odd-even hopscotch algorithm by a half time step and applied different formulas in different stages. First, we tested 105 = 100,000 different algorithm combinations in case of small systems with random parameters, and then examined the competitiveness of the best algorithms by testing them in case of large systems against popular solvers. These tests helped us find the top five combinations, and showed that these new methods are, indeed, effective since quite accurate and reliable results were obtained in a very short time. After this, we verified these five methods by reproducing a recently found non-conventional analytical solution of the heat equation, then we demonstrated that the methods worked for nonlinear problems by solving Fisher’s equation. We analytically proved that the methods had second-order accuracy, and also showed that one of the five methods was positivity preserving and the others also had good stability properties.

A NSFD method for the singularly perturbed Burgers-Huxley equation

“…For future work, one may proceed with a higher order scheme for the problem under consideration for example by using the Crank-Nicolson method or those presented in [20,21]. We are currently working in this direction.…”

Comparison of the Performance of New and Traditional Numerical Methods for Long-Term Simulations of Heat Transfer in Walls with Thermal Bridges

“…We started working on new explicit schemes a couple of years ago for determining heat conduction in any number of spatial dimensions. In our original articles [33][34][35][36][37][38][39][40], the novel algorithms were theoretically and experimentally investigated. They were proven to be stable for the linear diffusion-or diffusion-reaction equation, and their theoretical order of convergence is also stated.…”

A New Parameter-Uniform Discretization of Semilinear Singularly Perturbed Problems

“…These problems arise in different fields of study such as fluid dynamics, magnetohydrodynamics, aerodynamics, oceanography, quantum mechanics, plasma dynamics, chemical reactions, and liquid crystal modeling [1][2][3]. As an example, the heat and mass transport phenomena [4] are described by singularly perturbed differential equations in which the diffusion coefficient is regarded as a perturbation parameter.…”