Scribd
Upload
27 views1 page

α = " cos + β cos 1 + β #: T = T  2 X 1 − (−1) nπ sinh sinh sin nπx L 

This document provides instructions for a computational fluid dynamics assignment to: 1) Obtain the steady-state temperature distribution on a 2D rectangular plate with specified boundary conditions. 2) Discretize the governing equation and solve using Jacobi, Gauss-Seidel, and SOR methods. 3) Compute the error between analytical and numerical solutions and the number of iterations for each method.

Uploaded by

Aditya AP
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
27 views1 page

α = " cos + β cos 1 + β #: T = T  2 X 1 − (−1) nπ sinh sinh sin nπx L 

This document provides instructions for a computational fluid dynamics assignment to: 1) Obtain the steady-state temperature distribution on a 2D rectangular plate with specified boundary conditions. 2) Discretize the governing equation and solve using Jacobi, Gauss-Seidel, and SOR methods. 3) Compute the error between analytical and numerical solutions and the number of iterations for each method.

Uploaded by

Aditya AP
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/ 1

NATIONAL INSTITUTE OF TECHNOLOGY CALICUT

Department of Mechanical Engineering

ME4027 Computational Fluid Dynamics, Assignment-1

It is required to obtain steady-state temperature distribution on a 2D rectangular


plate. Governing equation (G.E.) and boundary conditions are:

∂ 2T ∂ 2T
∇2 T = + =0
∂x2 ∂y 2

y = 0, T = T1 ; x = 0, T = T2 ; y = H = 2m, T = T3 ; x = L = 1m, T = T4 where

T1 = 150◦ C, T2 = T3 = T4 = 0. Analytical solution is given by

   
∞ nπ(H−y)
X 1 − (−1)n sinh L
 nπx 
T = T1 2 nπH
 sin 
n=1
nπ sinh L
L

consider upper limit of summation ∞ = n = 110. Discretize the G.E. and solve by

Jacobi, Gauss-Seidel and Successive over-relaxation (SOR) at a convergence level

of max = 0.01, where the error () in successive iterations (l + 1 & l) is defined as

i=N −1
X−1
j=M
l+1 l
= |Ti,j − Ti,j |
j=2
i=2
p
2 − 2 (1 − α) ∆x
Use N = 21, M = 41 and ωoptimal = for SOR, where β = and
α ∆y

"
π
 2 π
 #2
cos N −1 + β cos M −1
α=
1 + β2
Compute error in analytical and numerical solutions, and number of iterations for all.

You might also like