MATLAB PDE BC'S

1 次查看（过去 30 天）
Mr.DDWW
Mr.DDWW 2022-3-25
回答： Mr.DDWW 2022-3-27
clc;clear all;close all;
L = 1;
x = linspace(0,L,75);
t = linspace(0,1,75);
m = 1;
sol = pdepe(m,@heatpde,@heatic,@heatbc,x,t);
sol1=1-sol;
figure(1)
surf(x,t,sol1);
xlabel('y/b');
zlabel('(T-T_0)/(T_1-T_0)');
title('Fig 12.1-1'); grid on;
function [c,f,s] = heatpde(x,t,u,dudx)
c = 1;
f = dudx;
s = 0;
end
function u0 = heatic(x)
u0 = 1;
end
function [pl,ql,pr,qr] = heatbc(xl,ul,xr,ur,t)
% left Bc = ul
pl = ul;
ql = 0;
% right BC= ur
pr = ur;
qr = 0;
end
I am a pde code. I am having a problem using the BC's in image. can you please help me.

请先登录，再回答此问题。

回答（2 个）

Torsten
Torsten 2022-3-25
function [pl,ql,pr,qr] = heatbc(xl,ul,xr,ur,t)
Nu = 1.0;
% left Bc = ul
pl = 0;
ql = 1;
% right BC= ur
pr = Nu*ur;
qr = 1;
end
  4 个评论
显示 2更早的评论
Mr.DDWW
Mr.DDWW 2022-3-26
Well, I am supposed to obtain the numerical solution from the image by changing the size of the mesh
Torsten
Torsten 2022-3-26
If it's the equation from the image you are trying to solve, you'll have to set m=0 instead of m=1 in your code.
If you want to solve the problem for different spatial meshes, change the "75" in
x = linspace(0,L,75);
to a different number.

请先登录，再进行评论。

Mr.DDWW
Mr.DDWW 2022-3-27
It is a slab. So the symmetry (m) = 1

类别

在 Help Center 和 File Exchange 中查找有关 PDE Solvers 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by