TransientThermalResults
瞬态热解决方案和派生的数量
描述
一个TransientThermalResults对象包含温度和梯度值的形式方便绘图和后处理。
温度及其梯度计算节点的三角形或四面体网格生成的generateMesh。节点出现在温度值温度财产。解决方案*中出现SolutionTimes财产。这三个组件的温度梯度节点出现在XGradients,YGradients,ZGradients属性。您可以提取解决方案和梯度值指定的时间指标温度,XGradients,YGradients,ZGradients。
插入温度及其梯度定制网格(例如,指定的meshgrid),使用interpolateTemperature或evaluateTemperatureGradient。
热解决方案的评估热通量在节点或任意空间位置、使用evaluateHeatFlux。评估集成指定边界热流率正常,使用evaluateHeatRate。
创建
解决一个瞬态热使用问题解决函数。这个函数返回一个瞬态热作为一个解决方案TransientThermalResults对象。
属性
对象的功能
evaluateHeatFlux |
评价热通量在节点或任意空间位置的热解决方案 |
evaluateHeatRate |
评估集成指定边界热流率正常 |
evaluateTemperatureGradient |
评估温度梯度在任意空间位置的热解决方案 |
interpolateTemperature |
插入温度在热的结果在任意空间位置 |
例子
解决瞬态热模型
解决一个二维瞬态热的问题。
关于这个问题创建一个瞬态热模型。
thermalmodel = createpde(热=“瞬态”);
创建几何和包括在模型中。
于SQ1 = [3;4;0;3;3;0;0;0;3;3); D1 = [2; 4; 0.5; 1.5; 2.5; 1.5; 1.5; 0.5; 1.5; 2.5]; gd = [SQ1 D1]; sf =“于SQ1 + D1”;ns = char (“于SQ1”,“D1”);ns = ns ';dl = decsg (gd、科幻、ns);geometryFromEdges (thermalmodel dl);pdegplot (thermalmodel EdgeLabels =“上”FaceLabels =“上”)xlim ([-1.5 - 4.5]) ylim([-0.5 - 3.5])轴平等的
正方形,分配这些热性能:
导热系数是
质量密度是
比热是
thermalProperties (thermalmodel ThermalConductivity = 10,…MassDensity = 2,…SpecificHeat = 0.1,…脸= 1);
钻石的地区分配这些热性能:
导热系数是
质量密度是
比热是
thermalProperties (thermalmodel ThermalConductivity = 2,…MassDensity = 1,…SpecificHeat = 0.1,…面对= 2);
假设菱形区域密度的热源 。
internalHeatSource (thermalmodel 4面= 2);
施加一个恒定的温度为0°C的平方板。
thermalBC (thermalmodel温度= 0,边缘= [1 2 7 8]);
设置初始温度为0°C。
thermalIC (thermalmodel 0);
生成网格。
generateMesh (thermalmodel);
这个问题的动力学是非常快。温度大概在0.1秒到达稳定状态。捕捉最活跃的动力的一部分,设置时间的解决方案logspace (2, 1, 10)。这个命令返回10对数间隔的解决方案*在0.01和0.1之间。
tlist = logspace (2, 1, 10);
解方程。
thermalresults =解决(thermalmodel tlist);
情节的解决方案与等温线用等高线图。
T = thermalresults.Temperature;msh = thermalresults.Mesh;要看更多有关憩苑pdeplot (, XYData = T(:, 10),轮廓=“上”ColorMap =“热”)
解决暂时的轴对称模型
分析传热与圆形截面杆和内部热代通过简化二维模型的三维轴对称模型。
创建一个瞬态热模型求解轴对称问题。
thermalmodel = createpde (“热”,“transient-axisymmetric”);
二维矩形带的模型x维度延伸轴的外表面,其对称性y维度的实际长度延伸杆(从- - - - - -1.5米到1.5米)。创建几何通过指定坐标的四个角。对于轴对称模型,工具箱假设转动轴是垂直轴穿过r= 0。
g = decsg ([3 4 0 0。2。2 -1.5 1.5 1.5 -1.5]”);
包括几何模型中。
geometryFromEdges (thermalmodel g);
情节的几何边缘标签。
图pdegplot (thermalmodel,“EdgeLabels”,“上”)轴平等的
杆是由这些材料的热性能。
k = 40;%导热,W / (m * C)ρ= 7800;^ %密度,公斤/米3cp = 500;%比热,W * s / C(公斤*)q = 20000;%热源,W / m ^ 3
指定导热系数、质量密度和比热的材料。
thermalProperties (thermalmodel“ThermalConductivity”、钾、…“MassDensity”ρ,…“SpecificHeat”,cp);
指定内部热源和边界条件。
internalHeatSource (thermalmodel q);
定义边界条件。没有热量转移的方向正常对称轴(边缘1)。您不需要更改默认这条边的边界条件。边2是保持在一个恒定的温度T= 100°C。
thermalBC (thermalmodel“边缘”2,“温度”,100);
指定外边界上的对流边界条件(3)边缘。外边界周围的温度是100°C和传热系数 。
thermalBC (thermalmodel“边缘”3,…“ConvectionCoefficient”,50岁,…“AmbientTemperature”,100);
杆的底部的热通量(4)边缘 。
thermalBC (thermalmodel“边缘”4“HeatFlux”,5000);
指定杆的初始温度是零。
thermalIC (thermalmodel 0);
生成网格。
generateMesh (thermalmodel);
次计算瞬态解的解决方案t = 0来t = 50000秒。
tfinal = 50000;tlist = 0:100: tfinal;结果=解决(thermalmodel tlist);
画出温度分布t = 50000秒。
T = result.Temperature;图pdeplot (thermalmodel,“XYData”T(:,结束),…“轮廓”,“上”)轴平等的标题(sprintf ([瞬态温度的…在最后一次(% g秒)”),tfinal))
