Numerical Solution of Parabolic Partial Differential Equation by Using Finite Difference Method
In the real world, many physical problems like heat equation, wave equation, Laplace equation and Poisson equation are modeled by partial differential equations (PDEs). A PDE of the form ut = α uxx, (α > 0) where x and t are independent variables and u is a dependent variable; is a one-dimensional heat equation. This is an example of a prototypical parabolic equation. The heat equation has analytic solution in regular shape domain. If the domain has irregular shape, computing analytic solution of such equations is difficult. In this case, we can use numerical methods to compute the solution of such PDEs. Finite difference method is one of the numerical methods that is used to compute the solutions of PDEs by discretizing the given domain into finite number of regions. Here, we derived the Forward Time Central Space Scheme (FTCSS) for this heat equation. We also computed its numerical solution by using FTCSS. We compared the analytic solution and numerical solution for different homogeneous materials (for different values of diffusivity α). There is instantaneous heat transfer and heat loss for the materials with higher diffusivity (α) as compared to the materials of lower diffusivity. Finally, we compared simulation results of different non-homogeneous materials.
All right reserved. No part of this Journal may be reproduced in any form or by any electronic or mechanical means, including information storage and retrieval system, without permission in writing from the publisher, except by a reviewer who may quote brief passage in a review. The views and interpretation in this journal are those of author(s) and they are not attributable to the NPS.