Scribd
Upload
4 views2 pages

Lampiran A Program Matlab: Ta %CM %CM %oc L cm-1 Oc-1 L cm-2 Oc-1 L cm-1 Oc-1

Uploaded by

Gusti Handayani
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
4 views2 pages

Lampiran A Program Matlab: Ta %CM %CM %oc L cm-1 Oc-1 L cm-2 Oc-1 L cm-1 Oc-1

Uploaded by

Gusti Handayani
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

LAMPIRAN A

PROGRAM MATLAB

A.1 Main Program 1 (Soal No.1)


%DATA
L = 10; %cm
D = 2.5; %cm
T1 = 250;%oC
k1 = 0.4; %cal cm-1 oC-1
h1 = 0.02; %cal cm-2 oC-1
k2 = 0.005; %cal cm-1 oC-1
Nu = h1*D/k1;
h2 = Nu*k2/D;
Tu = 30;
Nz = 100;
z = linspace(0,L,Nz+1);
dz = L/Nz;
n = 6/dz+1; %posisi pada z=6
A = zeros(Nz-1,Nz-1);
B = zeros(Nz-1,1);

%Calculation
A(1,1) = 1;
B(1) = T1;
A(Nz+1,Nz-1) = -k2;
A(Nz+1,Nz) = 4*k2;
A(Nz+1,Nz+1) = -3*k2-2*dz*h2;
B(Nz+1) = 2*dz*h2*Tu;
for i = 2:n-1
A(i,i-1) = 1;
A(i,i) = -(4*h1*dz^2/(k1*D)+2);
A(i,i+1) = 1;
B(i) = -(4*h1*dz^2/(k1*D)*Tu);
end
for i = n+1:Nz
A(i,i-1) = 1;
A(i,i) = -(4*h2*dz^2/(k2*D)+2);
A(i,i+1) = 1;
B(i) = -(4*h2*dz^2/(k2*D)*Tu);
end

A(n,n-2) = k1/k2;
A(n,n-1) = -4*k1/k2;
A(n,n) = 3*(k1/k2+1);
A(n,n+1) = -4;
A(n,n+2) = 1;
B(n) = 0;
T=A\B;

%PLOTTING
figure (1)
plot(z,T)
title('Profil Suhu Batang Pada Berbagai Posisi')
xlabel('Posisi ( cm )')
ylabel('Suhu ( ^0C )')

A.2 Subroutine 2 (Soal No.2)


function dvdt=Pensil07(t,v,Nx,n,miuw,miuo,dx)
dvdt=zeros(Nx+1,1);
v(1) = (-v(3)+4*v(2))/3;
v(end) = 8;
v(n) = (4*v(n+1)-v(n+2)-miuw/miuo*(-4*v(n-1)+v(n-2)))./(3*(miuw/miuo+1));
for i=2:n-1
dvdt(i)=miuw*(v(i+1)-2*v(i)+v(i-1))/dx^2;
end
for i=n+1:Nx
dvdt(i)=miuo*(v(i+1)-2*v(i)+v(i-1))/dx^2;
end
end

A.3 Main Program 2 (Soal No.2)


%DATA
miuw = 1; %cP
miuo = 3; %cP
L = 10; %cm
Nx = 100;
x = linspace(0,L,Nx+1);
dx = L/Nx;
n = 6/dx+1; %posisi pada x=6
tf = 30;
Nt = 60;
tspan = linspace(0,tf,Nt);
IC = zeros(1,Nx+1);

%Solver
[t,v] = ode15s(@Pensil07,tspan,IC,[],Nx,n,miuw,miuo,dx);
v(:,1) = (-v(:,3)+4*v(:,2))/3;
v(:,end) = 8;
v(:,n) = (4*v(:,n+1)-v(:,n+2)-miuw/miuo*(-4*v(:,n-1)+v(:,n-
2)))./(3*(miuw/miuo+1));

%PLOTTING
figure(1)
imagesc(x,t,v)
title('Profil Kecepatan Fluida')
xlabel('Posisi ( m )')
ylabel('Waktu ( detik )')
colormap jet
colorbar

You might also like