thermal conductivity, pde toolbox
조회 수: 15 (최근 30일)
이전 댓글 표시
Dear all,
I have a question concerning the pde tool box. Indeed, I would like to model the thermal behaviour of a cyclindrical battery. For said purpose, I would like to specify the thermal conductivities in the x and y directions. How can I specify said parameters using the following code?
specifyCoefficients(Thermal, 'm',0,...% 'd',(para.rho)/para.cp,...%rho*cp 'c',???,...%k 'a',0,...%0 'f',Heat)%source
Thanks in advance.
댓글 수: 0
채택된 답변
Ravi Kumar
2018년 10월 15일
Hi Juan,
In the ThermalModel based workflow you can use the function thermalProperties to define thermal conductivity, mass density, and specific heat, all of which can also be function spatial coordinates (x,y) and the temperature itself (nonlinear).
For example if your thermal conductivity is linear function of x & y, then you can define the thermal conductivity as a function:
function k = kFunc(region,state)
k = 1 + 0.1*region.x + 0.01*region.y;
end
You can then specify kFunc in thermalProperties function as in the example below:
model = createpde('thermal','transient');
importGeometry(model,'Block.stl');
thermalProperties(model,'ThermalConductivity',@kFunc,'MassDensity',1,'SpecificHeat',1)
generateMesh(model);
thermalBC(model,'Face',1,'Temperature',1)
thermalIC(model,0);
R = solve(model,[0,0.01])
pdeplot3D(model,'ColorMapData',R.Temperature(:,2))
Regards, Ravi
댓글 수: 3
Ravi Kumar
2018년 10월 15일
Hi Juan,
You should be able to use ThermalModel workflow in R2017a version of MATLAB.
Regards,
Ravi
추가 답변 (1개)
Ravi Kumar
2018년 10월 16일
Hi Juan,
Looks like you have orthotropic thermal conductivity, not spatially varying conductivity. in that case, you can just use the matrix form of thermal conductivity, k = [0.15, 0 ; 0, 30].
thermalProperties(model,'ThermalConductivity',k,'MassDensity',1,'SpecificHeat',1)
Note if the problem is in 3-D space, then you need to specify a 3x3 matrix of conductivity. You don't need a function for the case you describe now.
You can specify heat source as a function of time using the function fFunc. Note that you can specify such heat source as internalHeatSource, if it is a domain heat source. If it is a boundary heat input, then you need to use 'HeatFlux' option of thermalBC function.
Regards,
Ravi
댓글 수: 2
Ravi Kumar
2018년 10월 22일
Hi Juan,
I think you might have figured out answer for your question by now. Solver works by discretization of the domain. So, solver needs value of heat source at several points on the boundary. The coordinates points at which heat source should be provided by the user are passed to the function as location.x, location.y, etc. So what you have done is correct.
Regards,
Ravi
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!