1D Fourier equation with temperature dependent thermal properties

2 次查看（过去 30 天）
Maria Lorena Richiusa
评论： Maria Lorena Richiusa 2020-6-7
Hello! I am trying to solve the 1D Fourier's equation with temperature dependent thermal properties but I do not know how to get the temperature results from pdepe at each time step and use them to evaluate the material properties for the next time step. Can you help me sort this out? Find below the example.
In this case, for example, the thermal diffisivity is temperature dependent accordint to the law:
diffr=3.34E-10*T^4 - 9.37E-07*T^3 + 9.38E-04*T^2 - 4.80E-01*T + 2.09E+02
function [c,f,s] = pdex1pde(x,t,u,dudx,diffr,Tfront)
c = 1/diffr;
f = dudx;
s = 0;
end
Many thanks for the help!

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

回答（1 个）

Bill Greene
Bill Greene 2020-6-6
function [c,f,s] = pdex1pde(x,t,T,dudx,diffr,Tfront)
diffr=3.34E-10*T^4 - 9.37E-07*T^3 + 9.38E-04*T^2 - 4.80E-01*T + 2.09E+02;
c = 1/diffr;
f = dudx;
s = 0;
end
  3 个评论
显示 1更早的评论
Bill Greene
Bill Greene 2020-6-7
Yes, you are free to name that third argument whatever you like-- u, T, or any other legal matlab variable name.
Maria Lorena Richiusa
Maria Lorena Richiusa 2020-6-7
Thanks again for your answer. Could I define the temperature-dependent material property law on an external script that recalls the pdepe function, and passing this function as an argument of the pdepe? I have tried but in a wrong way, as the temperature is a result of the pdepe function.
Thanks a lot for your answer.

请先登录，再进行评论。

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by