Thanks to visit codestin.com
Credit goes to www.scribd.com

Scribd
Upload
54 views4 pages

Course Project Report

The document describes a numerical simulation of heat transfer between two parallel plates with different temperatures, where one plate is stationary and the other is moving. The goal is to compare the analytical and numerical solutions for the non-dimensional temperature profile and analyze the effect of the Eckert number. The finite difference method is used to numerically solve the governing equations and the results are plotted against the analytical solution and the error is calculated.

Uploaded by

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

Course Project Report

The document describes a numerical simulation of heat transfer between two parallel plates with different temperatures, where one plate is stationary and the other is moving. The goal is to compare the analytical and numerical solutions for the non-dimensional temperature profile and analyze the effect of the Eckert number. The finite difference method is used to numerically solve the governing equations and the results are plotted against the analytical solution and the error is calculated.

Uploaded by

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

COURSE PROJECT-1

Consider fully developed flow between two large parallel plates separated by a distance L. The
upper plate is moving with a velocity, U and the lower plate is kept stationary. The lower and
upper plate is maintained at temperaturesT 0and T1respectively. Numerically simulate this problem
using finite difference method and validate the results with analytical solution. Also plot the
effect of Eckert number on the non-dimensional temperature profile.

Fig 1. Schematic of the physical domain

Aim: To compare analytical and numerical solution of coutte flow where two plates are kept at
constant wall temperature.

Solution:
The problem is first solved analytically which yields the following temperature and
velocity profile. It should be noted that the profiles here are nondimensional velocity and
temperatures infact all variables have been nondimensionalised which are denoted by uppercase
characters.
The differential equation that governs the flow is
2

d U
=0
2
dY
Which yields the velocity profile to be linear and on applying no slip boundary conditions we get.

V =Y
The governing equations for temperature is
2

d
=EcPr dV
2
dY
dY

( )

For nondimensional variables

dV
dY

is 1 and hence the equation simplifies to

d
=EcPr
dY 2
Solving this numerically using FDM using the code given below in MATLAB

function HT_Course_Project
clc;
format short;
l=1;
ecpr=input('Enter the value of Brinkman number');
h=input('Enter the value of step size');
Space=[0:h:l];
temp=size(Space);
n=temp(1,2);
T=zeros(1,n);
T(1,1)=0;
T(1,n)=1;
A=zeros(n-2);
B=zeros(1,n-2);
B=transpose(B);
A(1,1)=-2/(h^2);
A(1,2)=1/(h^2);
B(1)=-ecpr;
A(n-2,n-2)=-2/(h^2);
A(n-2,n-3)=1/(h^2);
B(n-2)=-ecpr-1/(h^2);
for i=2:n-3
A(i,i-1)=1/(h^2);
A(i,i)=-2/(h^2);
A(i,i+1)=1/(h^2);
B(i)=-ecpr;
end
T(1,2:n-1)=Solution_of_Tridiagonal_Matrix(A,B)
Actual=-0.5*ecpr*Space.^2+(1+0.5*ecpr)*Space;
for i=1:n
error(i)=Actual(i)-T(i);
end
grid on;
plot(Space,T,'b',Space,Actual,'r',Space,error,'g');
legend('Numerical','Actual','Error','Location','northwest')
end

Comparison between analytical and numerical solution and effect of Brinkman number is shown in
following figures.

Nondimensional Temperature

Nondimensional Plate Length


Graph for actual and numerical solution and error between is shown in figure with EcPr=2

Nondimensional Temperature

Nondimensional Plate Length

Effect of different values of EcPr is shown in figure with EcPr=1,2 and 3

You might also like