摘要边界条件设置为两端dT / dx = 0。你在右端设置T = 293。
一维热传导由PDE估计
79(30天)
显示旧的评论
我试图复制图3 (b)。
https://pubs.rsc.org/en/content/articlelanding/2016/RA/c6ra21902k
这是支持信息,展示了如何万博1manbetx计算。
https://www.rsc.org/suppdata/c6/ra/c6ra21902k/c6ra21902k1.pdf
看来,温度太高,除非我s方程除以~ 20。
当我试着模拟其他材料,入射激光的能量似乎会增加20倍。
我不知道为什么的原因。边界条件也影响呢?
如果你有时间,我真的很感激如果你能给我你的建议。
clearvars
clearvars变量
clearvars刚keepVariables
clearvars刚keepVariables变量
x = linspace (0、4000、4000);
t = linspace (0500500);
m = 0;
索尔= pdepe (m, @heatpde, @heatic @heatbc x, t);
colormap默认的
显示亮度图像(x, t,索尔);
爬(1023 [273]);
colorbar
包含(“距离x”,“翻译”,“乳胶”)
ylabel (“t”,“翻译”,“乳胶”)
标题(“热方程”)
函数[c、f、s] = heatpde (x, t, u, dudx)
%的材料属性
n = 1000000000;%到纳米级
一个= 1776000 / n;%的吸收系数
r = 0.0377;%反射
k = 112.4 / n / n;%导热系数
d = 2260000 / n / n / n;%质量密度
h = 0.72;%的热容
v = 26;%脉冲持续时间
b = 6;%脉冲形状的常量
我= 2000 / n / n;%激光强度
c = d * h;
f = k * dudx;
s =(我)/ 20 *功率(t / v, b) * (exp (b * (1 - t / v))) *(第一轮)* (exp ((——) * x));
结束
函数情况= heatic (x)
情况= 293;
结束
函数[pl, ql,公关,qr] = heatbc (xl, ul, xr, ur, t)
pl = 0;
ql = 1;
公关= ur - 293;
qr = 0;
结束
