Computational Analysis of the Stability of 2D Heat Equation on Elliptical Domain Using Finite Difference Method
Asian Research Journal of Mathematics,
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.
- Modeling and simulation
- stability analysis
- computational analysis
- finite difference method
- elliptical domain
- heat equation.
How to Cite
[Retrieved in May 2016]
Olsen-Kettle L. Numerical solution of partial differential equations; 2011.
Isaacson E, Keller HB. Analysis of Numerical Methods. 1994;541.
Han Chen, Cho Hong Min, Frederick Gibou. A Supra-convergent finite difference scheme for the poisson and heat equations on irregular domains and Non-Graded Adaptive Grids; 2006.
Peter McCorquodale, Phillip Colella, Hans Johansen, A cartesian grid embedded boundary method for the heat equation on irregular domains.
Verena Hora, Peter Gruber. Parallel Numerical Solution of 2-D Heat Equation, Parallel Numerics. 2005;47-56.
Gerald W. Recktenwald. Finite-difference approximations to the heat equation; 2011.
Kazufumi Ito, Zhilin Li, Yaw Kyei. Higher-order, Cartesian grid based finite difference schemes for elliptic equations on irregular domains. SIAM Journal on Scientific Computing. 2005;27(1): 346-367.
Izadian J, Karamooz N. New method for solving poisson equation on irregular domains. Applied Mathematical Sciences. 2012;6(8):369–380.
Gavete L, Benito JJ, Urena F. Generalized finite differences for solving 3D elliptic and parabolic equations. Applied Mathematical Modelling. 2016;40(2):955-965.
Gilberto E. Urroz. Convergence, stability, and consistency of finite difference schemes in the solution of partial differential equations; 2004.
Courant R, Friedrichs K, Lewy H. Über die partiellen Differenzengleichungen der mathematischen Physik", Mathematische Annalen (in German). 1928;100(1):32-74.
Courant R, Friedrichs K, Lewy H. On the partial difference equations of mathematical physics. AEC Research and Development Report, NYO-7689, New York: AEC Computing and Applied Mathematics; 1956.
Courant R, Friedrichs K, Lewy H. On the partial difference equations of mathematical physics. IBM Journal of Research and Development. 1928;11(2):215-34.
Charney JG, Fjörtoft R, Von Neumann J. Numerical integration of the barotropic vorticity equation. In The Atmosphere—A Challenge. American Meteorological Society, Boston, MA. 1990;267-284.
Lax PD, Richtmyer RD. Survey of the stability of linear finite difference equations. Comm. Pure Appl. Math. 1956;9:267-293 MR 79204.
Crank J, Nicolson P. A practical method for numerical evaluation of solutions of partial differential equations of heat conduction type. Proc. Camb. Phil. Soc. 1947;43.
Yuste SB. Weighted average finite difference methods for fractional diffusion equations. Journal of Computational Physics. 2006;216:264–274.
Yuste SB, Acedo L. An explicit finite difference method and a new von neumann-type stability analysis for fractional diffusion equations. SIAM Journal on Numerical Analysis. 2005;42(5):1862-1874.
Bilbao S. Modeling of complex geometries and boundary conditions in finite difference/finite volume time domain room acoustics simulation. In IEEE Transactions on Audio, Speech, and Language Processing. 2013;21(7):1524-153.
Jincheng Ren, Zhi-zhong Sun, Xuan Zhao. Compact difference scheme for the fractional sub-diffusion equation with Neumann boundary conditions. Journal of Computational Physics. 2013; 232(1):456-467.
Ercília Sousa, Can Li. A weighted finite difference method for the fractional diffusion equation based on the Riemann–Liouville derivative. Applied Numerical Mathematics. 2015;90:22-37.
Ercília Sousa. A second order explicit finite difference method for the fractional advection diffusion equation, Computers & Mathematics with Applications. 2012;64(10):3141-3152.
Sweilam NH, Khader MM, Mahdy AMS. Crank-nicolson finite difference method for solving time-fractional diffusion equation. Journal of Fractional Calculus and Applications. 2012;2(2):1-9.
Quintana-Murillo J, Yuste SB. A finite difference method with non-uniform time steps for fractional diffusion and diffusion-wave equations. The European Physical Journal Special Topics. 2013;222(8): 1987–1998.
Norihiro Watanab, Olaf Kolditz. Numerical stability analysis of two-dimensional solute transport along a discrete fracture in a porous rock matrix, Water Resources Research. 2015;51(7):5855–5868.
Jiequan Li, Zhicheng Yang. The von Neumann analysis and modified equation approach for finite difference schemes. Applied Mathematics and Computation. 2013;225:610-621.
Ehlers W, Zinatbakhsh S, Markert B. Stability analysis of finite difference schemes revisited: A study of decoupled solution strategies for coupled multifield problems. Int. J. Numer. Meth. Engng. 2013; 94:758–78.
Baudouin L, Seuret A, Gouaisbaut F, Dattas M. Lyapunov stability analysis of a linear system coupled to a heat equation. IFAC-PapersOnLine. 2017;50(1):11978-11983.
Konangi S, Palakurthi NK, Ghia U. Von Neumann stability analysis of first-order accurate discretization schemes for one-dimensional (1D) and two-dimensional (2D) fluid flow equations. Computers & Mathematics with Applications. 2018;75(2):643-665.
Jahangir K, Rehman SU, Ahmad F, Pervaiz A. Sixth-Order Stable Implicit Finite Difference Scheme for 2-D Heat Conduction Equation on Uniform Cartesian Grids with Dirichlet Boundaries. Punjab Univ. J. Math. 2019;51(5):27-42.
Singh AK, Bhadauria BS. Finite difference formulae for unequal sub-intervals using Lagrange’s interpolation formula. Int. J. Math. Anal. 2009;3(17):815.
Abstract View: 1892 times
PDF Download: 1229 times