Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd
Upload
46 views

Function Global: %solution For Dimensionless Concentration in Function of Eta

The document contains MATLAB code that models solute transport using a 1D convection-diffusion partial differential equation (PDE) under various conditions of the Biot and Fourier numbers (Bi, phi). It generates concentration profiles by running the PDE solver for different parameter values and plots the dimensionless concentration versus spatial and temporal variables. The code contains functions to define the PDE, initial and boundary conditions, and solves the PDE using the pdepe solver to generate concentration profiles as surfaces and curves.

Uploaded by

Pedro Henrique Magacho
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
46 views

Function Global: %solution For Dimensionless Concentration in Function of Eta

The document contains MATLAB code that models solute transport using a 1D convection-diffusion partial differential equation (PDE) under various conditions of the Biot and Fourier numbers (Bi, phi). It generates concentration profiles by running the PDE solver for different parameter values and plots the dimensionless concentration versus spatial and temporal variables. The code contains functions to define the PDE, initial and boundary conditions, and solves the PDE using the pdepe solver to generate concentration profiles as surfaces and curves.

Uploaded by

Pedro Henrique Magacho
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

function HW_2

global Bi phi Kv D n i j Concentration


Bi = [0.1 1 10 100];
phi = [0.1 1 3 10 100];
n = 0:0.001:1;
Kv = 0.03;
D = 0.00001;
k=1;
for i=1:4
for j=1:5
Conc
figure(1)
subplot(2,2,k)
plot(n, Concentration,'LineWidth',2)
grid on
set(gca,'fontsize',12)
hold on
end
hold off
title(sprintf('Dimensionless Concentration Profile for Biot =
%.2f',Bi(i)))
legend('\phi=0.1','\phi=1','\phi=3','\phi=10','\phi=100')
xlabel('\eta','fontsize',14)
ylabel('Concentration','fontsize',14)
k=k+1;
end
Conc_phi
pde
end
function Conc
global Bi phi Concentration i j n
Constant= Bi(i)/(-phi(j)*exp(-phi(j))+phi(j)*exp(phi(j))+Bi(i)*exp(phi(j))+Bi(i)*exp(phi(j)));
Concentration = Constant*(exp(-phi(j)*n)+exp(phi(j)*n)); %Solution for
Dimensionless Concentration in function of eta
end
function Conc_phi
figure(2)
f=0:0.01:4;
Conc = 2./(exp(-f)+exp(f));
plot(f,Conc,'LineWidth',2)
title('Dimensionless Concentration Profile in function of \phi')
xlabel('\phi','fontsize',14)
ylabel('Concentration','fontsize',14)
legend('\eta = 0')
set(gca,'fontsize',12)
grid on
end

function pde
m=0;
x = 0:0.01:1;
t = 0:0.01:1;
sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,t,x);
Concentration = sol(:,:,1);
figure(3)
surf(x,t,Concentration)
set(gca,'fontsize',12)
xlabel('\eta','fontsize',14)
ylabel('\tau','fontsize',14)
zlabel('Concentration','fontsize',14)
axis([0,1,0,1,0,1])
title('Transient Dimensionless Concentration Profile')
figure(4)
plot(t,Concentration(:,1),t,Concentration(:,51),'LineWidth',2)
title('Dimensionless Concentration Profile in function of \tau')
xlabel('\tau')
ylabel('Concentration')
legend('z = 0','z = L/2')
set(gca,'fontsize',12)
axis([0,1,0,0.25])
end
function [c,f,s] = pdex1pde(x,t,u,DuDx)
c = 1; %Building pde
f = DuDx; %Building pde
s = -9*u; %Building pde
end
function u0 = pdex1ic(x)
u0 = 0;%Initial Condition
end
function [pl,ql,pr,qr] = pdex1bc(xl,ul,xr,ur,t)
pl = 0; %Boundary condition no flux on the z = 0
ql = 1; %Boundary condition no flux on the z = 0
pr = ur-1; %Boundary condition Dimensionless Concentration = 1 for eta
= 1
qr = 0; %Boundary condition Dimensionless Concentration = 1 for eta =
1
end

Published with MATLAB R2015a

You might also like