Scribd
Upload
77 views2 pages

Assignment 1: Parabolic PDE's: 1. 1D Heat Conduction

This document outlines an assignment to solve a 1D heat conduction problem using two numerical methods: an explicit FTCS method and an implicit Crank-Nicolson method. Students are asked to write code to implement the methods, compare results to analytical solutions, and analyze stability and error properties of the methods.

Uploaded by

Nurrohman Nurrohman
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
77 views2 pages

Assignment 1: Parabolic PDE's: 1. 1D Heat Conduction

This document outlines an assignment to solve a 1D heat conduction problem using two numerical methods: an explicit FTCS method and an implicit Crank-Nicolson method. Students are asked to write code to implement the methods, compare results to analytical solutions, and analyze stability and error properties of the methods.

Uploaded by

Nurrohman Nurrohman
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

King Abdulaziz University

Deprtment of Aerospace Engineering


AE-741 Computational Fluid Dynamics I.
Fall 2020

Assignment 1: Parabolic PDE’s

1. 1D Heat conduction
Write A computer code to determine the coefficients for the first order-derivative The heat
equation (1) governs the time-dependent temperature distribution in a homogeneous constant
property solid under conditions where the temperature varies only in one space dimension.
∂T ∂2 T
=ν 2 (1)
∂t ∂x
Consider a large wall of thickness L whose initial temperature is given by T (t, x) = c sin (πx/L)
If the faces of the wall continue to be held at 0◦ , then a solution for the temperature at t > 0,
0 ≤ x ≤ L is
−απ2 t
!  πx 
T x) = c exp
(t, sin (2)
L2 L
For this problem, let c = 100◦C, L = 1m, α = 0.02m2 /h. We will consider two methods of
solution:
(A) FTCS explicit method
(B) Crank-Nicolson implicit method

1.1. Part I:
- Classifiy the type of error for both methods by deriving and examining the modified equa-
tion.
1.2. Part II:
- Write a computer code to solve the problem described previously by methods A and B.
Make at least the following comparisons:
1. For ∆x = 0.1, ∆t = 0.1, compare the results from methods A and B and the exact solution
for t = 10h.
2. Repeat the aforementioned comparison after refining the space grid, that is, let ∆x = 0.01
(10 increments). Is the reduction in error as suggested by O[(∆x)2 ]? use L2 norm to show
this
3. For ∆x = 0.1, choose ∆t such that α∆t/(∆x)2 = 0.5 and compare the predictions of
methods A and B and the exact solution for t = 10h.
4. Demonstrate that method A does become unstable as α∆t/(∆x)2 exceeds 0.5. One sug-
gestion is to plot the centerline temperature versus time for α∆t/(∆x)2 ≈ 0.6 for 10 − 20h
of problem time.
5. For ∆x = 0.1, choose ∆t such that α∆t/(∆x)2 = 0.6 and compare the results of method B
and the exact solution for t ≈ 10h.
6. Increment α∆t/(∆x)2 to 2, then 3, etc., and repeat comparison (5) until the agreement
with the exact solution becomes noticeably poor.
7. Replace the second order discretization in FTCS with a fourth-order accurate discretiza-
tion.
8. Derive the stability condition
9. Repeat the above analysis (1 through 4) for The FTCS with the fourth-order accurate
discretization of the second order derivative.
Note:
• Comparison should be qualitative as well as quantitative
• Thorough analysis and discussion is expected.
• The format of the report should resemble the style of scientific papers.

You might also like