Defining a geometry to apply Heat flux - PDEPE Toolbox

Question

Hung Nguyen 2020-3-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/509005-defining-a-geometry-to-apply-heat-flux-pdepe-toolbox

评论： Hung Nguyen 2020-3-4

Hi everyone,

I am making a heat transfer model that shows the impact of laser beam on a small, spherical particle using PDEPE in Matlab. Diameter of the particle is around 60 um, while the laser spot radius is 250 um.

I have created a model after reading the Matlab's example (https://www.mathworks.com/help/pde/ug/heat-conduction-in-multidomain-geometry-with-nonuniform-heat-flux.html).

thermalBC(thermalmodel,'Face',1,'HeatFlux',11713803.81)

If I define the heat flux as above, PDEPE model will assume that the heat flux is apply to all spots around the sphere. However, in my case, the heat flux only apply to half of the sphere - as the laser beam only shine through one side of the particle. How can I define the laser impact area, in order to define the area that the heat flux apply to?

Thank you everyone.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Ravi Kumar 2020-3-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/509005-defining-a-geometry-to-apply-heat-flux-pdepe-toolbox#answer_418515

Exactly as in the example you pointed to, in which you will find:

thermalBC(thermalmodel,'Face',4,'HeatFlux',@externalHeatFlux,'Vectorized','on');

if you look at the externalHeatFlux, it applies heat flux only to certain region based on the value of phi. Modify this function to suite your needs.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Ravi Kumar 2020-3-4

Qflux(ids) = theta(ids).^2.*(pi - theta(ids)).^2.*phi(ids).^2.*(pi - phi(ids)).^2;

The above line defines value of heat flux as RHS expression. In you case, once you know the ids and defined a zero vector of Qflux, as in the example you can assign the values to the ids as:

Qflux(ids) = 11713803.81

Regards,

Ravi

Hung Nguyen 2020-3-4