The Numerical Method of Lines for Partial （行部分的数值方法）.pdf

The Numerical Method of Lines for Partial Differential Equations by Michael B. Cutlip, University of Connecticut and Mordechai Shacham, Ben-Gurion University of the Negev The method of lines is a general technique for solving partial differential equations (PDEs) by typically using finite difference relationships for the spatial derivatives and ordinary differential equations for the time derivative. William E. Schiesser at Lehigh University has been a major proponent of the numerical method of lines, NMOL.1 This solution approach can be very useful with undergraduates when this technique is implemented in conjunction with a convenient ODE solver package such as POLYMATH.2 A Problem in Unsteady-State Heat Transfer3 This approach can be illustrated by considering a problem in unsteady-state heat conduction in a one-dimensional slab with one face insulated and constant thermal conductivity as discussed by Geankoplis.4 Unsteady-state heat transfer in a slab in the x direction is described by the partial differential equation 2 T T a 2 (1) t x 2 where T is the temperature in K, t is the time in s, and is the thermal diffusivity in m /s given by k/c . In this treatment, the thermal conductivity k in W/m·K, the density in p 3 kg/m , and the heat capacity c in J/kg·K are all considered to be constant. p Consider that a slab of material with a thickness 1.00 m is supported on a nonconducting insu