Computational Analysis of the Stability of 2D Heat Equation on Elliptical Domain Using Finite Difference Method

Abstract

Aims: The aim and objective of the study to derive and analyze the stability of the finite difference schemes in relation to the irregularity of domain. Study Design: First of all, an elliptical domain has been constructed with the governing two dimensional (2D) heat equation that is discretized using the Finite Difference Method (FDM). Then the stability condition has been defined and the numerical solution by writing MATLAB codes has been obtained with the stable values of time domain. Place and Duration of Study: The work has been jointly conducted at the MUET, Jamshoro and QUEST, Nawabshah Pakistan from January 2019 to December 2019.  Methodology: The stability condition over an elliptical domain with the non-uniform step size depending upon the boundary tracing function is derived by using Von Neumann method. Results: From the results it was revealed that stability region for the small number of mesh points remains larger and gets smaller as the number of mesh nodes is increased. Moreover, the ranges for the time steps are defined for varied spatial step sizes that help to find the stable solution. Conclusion: The corresponding stability range for number of nodes N=10, 20, 30, 40, 50, and 60 was found respectively. Within this range the solution remains smooth as time increases. The results of this study attempt to provide the stable solution of partial differential equations on irregular domains.