[A-10]乘员舱热舒适分析技术

CAE－成本节约和技术革新的关键技术乘员舱热舒适性分析技术～STAR-CCM+的MRT、湿度、PMV的预测～CDAJ-China技术部唐连伟相羽俊輝(CDAJ)-2-概述对于住宅、办公室、学校和大厅等建筑物内的舒适性分析的已经很常见了，近来对于乘员舱内热舒适性的需求越来越多。乘员舱内空间一般相对建筑物内空间狭小，人员更容易受到来自窗外光照的影响，同时空调口吹风的分布对乘员舱内的环境的影响也很大，因此舱内热舒适性的评价是设计人员必须考虑的因素，也是提升设计附加值的方法。在此我们介绍如何使用STAR-CCM+进行乘员舱热舒适性分析评价。-3-目录1.热舒适性计算模型介绍2.舒适性评价使用的指标3.各指标在STAR-CCM+中的设定4.范例的展示5.参考文献6.致谢-4-•计算模型概要–对成员舱内的驾驶员、副驾驶员进行PMV/PPD的舒适性预测。•基础解析内容–稳态–湍流（Realizablek-ε）–多组分气体（空气，水蒸气）–密度：满足理想气体状态方程–重力–辐射（S2S，灰体，日照）–PassiveScalar•网格类型–TrimMesh1.舒适性计算模型介绍-5-•几何说明驾驶员副驾驶员前窗玻璃后窗玻璃侧窗玻璃1.舒适性计算模型介绍-6-•边界条件速度入口流量分配出口1.0[m/s]25.0[℃]质量分数Air[0.9902197]H2O[0.0097803]→相对湿度50[%]流出比1.01.舒适性计算模型介绍-7-•边界条件车体无滑移25.0[℃]反射比/发射率/透射比=0.2/0.8/0.0玻璃（前/后/侧）无滑移25.0[℃]反射比/发射率/透射比=0.1/0.4/0.5成员（驾驶员，副驾驶）无滑移绝热反射比/发射率/透射比=0.8/0.2/0.01.舒适性计算模型介绍-8-•日照设定1.舒适性计算模型介绍-9-•网格数：888,336•边界层：5.0mm(2层）•全局尺寸：100mm•最小尺寸：10mm1.舒适性计算模型介绍-10-2.舒适性评价使用的指标通常的舒适性评价对应以下这些评价指标：-预期平均通感[PredictedMeanVote(PMV)]-不适人员比例[Percentageofdissatisfiedpersons(PDD)]-平均辐射温度[MeanRadiantTemperature(MRT)]-相对湿度[RelativeHumidity(RH)]-平均空气龄[MeanAgeofAir(MAA)]-11--预期平均通感：PredictedMeanVote(PMV)-PMV是由丹麦工科大学的Fanger教授在1967年发表的热舒适方程式所提出的温冷感指标。其中涉及室温、平均辐射温度、相对湿度、平均风速、着衣量和人体作业量六个变量。使用从-3(冷感)～0(中间感受)～+3（热感）共七个阶段的舒适性评价。-ISO标准中，PMV在±0.5以内状态为推荐热舒适状态。2.舒适性评价使用的指标-12--预期平均通感：PredictedMeanVote(PMV)-定义式：-各变量说明-M：体表面的代谢率[W/m2]，其基本单位1Met=58.15[W/m2]。睡眠时约0.8[Met]，办公室工作时约1.2[Met]。平均成人体表面机约1.7[m2]，此时办公室工作的场合，代谢量约120[W]。-W：人体所做的机械功[W/m2]。通常取0。-H：体表面的对流、辐射的热损失[W/m2]-Ec：皮肤上的散热和出汗散热量[W/m2]-Cres：呼吸中的显热损失量[W/m2]-Eres：呼吸中的潜热损失量[W/m2]()()[]resrescMECEHWMePMV−−−−−×+=−028.0303.0036.02.舒适性评价使用的指标-13-()()[]resrescMECEHWMePMV−−−−−×+=−028.0303.0036.0fcl：裸露面积与着衣面积的比例。Icl：服装的热阻[m2K/W]tcl：衣服的表面温度[℃]tr：平均辐射温度[℃]Var：风速[m/s]hc：对流换热系数[W/m2K]ta：室温[℃]-预期平均通感：PredictedMeanVote(PMV)-H：体表面的对流、辐射的热损失2.舒适性评价使用的指标-14-()()[]resrescMECEHWMePMV−−−−−×+=−028.0303.0036.0pa：水蒸气分压力[Pa]-预期平均通感：PredictedMeanVote(PMV)-Ec：皮肤上的散热量-Cres：呼吸中的显热损失量-Eres：呼吸中的潜热损失量2.舒适性评价使用的指标-15--不适人员比例：Percentageofdissatisfiedpersons(PDD)-PPD指得是，对在一特定环境感到不适人员的比例。-定义式-ISO标准中，PDD小于10%的状态为推荐热舒适状态。()242179.003353.095100PMVPMVePDD+−−=2.舒适性评价使用的指标-16--平均辐射温度：MeanRadiantTemperature(MRT)-平均辐射温度指的是，人体实际周围全方向受到的辐射热等价至一个均一温度的假象空间的所具有的表面温度，近似等于周边壁面的平均表面温度。-定义式2.舒适性评价使用的指标-17--相对湿度：RelativeHumidity(RH)-相对湿度指得是，水蒸气分压力同当前温度下饱和水蒸气压力之比。-定义式-水的饱和蒸汽压计算公式使用JIS采用的SONNTANG公式-水蒸气的分压力即，(压力)×(水蒸气的摩尔分数)satOHPPRH2=()()()()]ln433502.25673952.12711193.22409642.219385.6096exp[21TTeTeTPsat+−+−−++−=−PH2O：水蒸气分压力[Pa]Psat：饱和水蒸气压[Pa]T：水蒸气温度[K]vOHXPP⋅=2P：圧力[Pa]Xv：水蒸气的摩尔分数2.舒适性评价使用的指标-18--平均空气龄：MeanAgeofAir（MAA）-空气龄指得是空气由进气口到达室内某一位置的移动时间。这种瞬态一样扩散假定的平均到达时间，即被称作平均空气龄。2.舒适性评价使用的指标-19-3.各指标在STAR-CCM+中的设定-预期平均通感：PredictedMeanVote(PMV)-不适人员比例：Percentageofdissatisfiedpersons(PDD)-对PMV，PDD的计算，可以在得到计算结果后通过在STAR-CCM+中运行Java宏文件来处理而获得。-20--平均辐射温度：MeanRadiantTemperature(MRT)1.Continua＞Physics1＞SelectModels…依次选择，“Radiation”Æ“Surface-to-surfaceRatiation”Æ“GrayThermalRadiation”Æ“SolarLoads”2.选择“Surface-to-surfaceRatiation”（S2S)的场合，“BoundaryMeanRadiantTemprature”就可以作为标准计算变量被选择了。①依次选择辐射，表面-表面辐射，灰体辐射，日照负荷②选择BoundaryMeanRadiantTemperature3.各指标在STAR-CCM+中的设定-21--相对湿度：RelativeHumidity(RH)1.Continua＞Physics1＞SelectModels…依次选择：“Multi-ComponentGas”Æ“Non-reacting”2.Continua＞Physics1＞Models1＞Multi-ComponentGas＞GasMixture＞右键点击GasComponents＞NewGasComponent＞选择‘Air(Air)’，重复前面步骤，选择‘H2O(Water)’。①选择多组分气体、无反应②气体成分选择Air，H2O3.各指标在STAR-CCM+中的设定-22--相对湿度：RelativeHumidity(RH)3.Tool＞FieldFunction＞New制作新的“场函数”用来监控“水蒸气分压力”4.右键点击UserFieldFunction1＞Rename…更名为“P_H2O”5.在Properties栏中，设定“FunctionName”、“Dimensions”和“Definition”“FunctionName”：PH2O“Dimensions”：Pressure=1“Definition”：$AbsolutePressure*$Species1MoleFraction※SpeciesX的X对应的数字指的是多组分气体组分定义的编号，这里需要根据H2O的Scalar编号一致。③新的FieldFunction作成④改名为“P_H2O”⑤定义函数3.各指标在STAR-CCM+中的设定-23--相对湿度：RelativeHumidity(RH)6.Tool＞FieldFunction＞New创建新的“场函数”--“饱和水蒸汽压”。7.右键点击UserFieldFunction1＞Rename…更名为“P_SAT”。8.在Properties栏中，设定“FunctionName”、“Dimensions”和“Definition”“FunctionName”：PSAT“Dimensions”：Pressure=1“Definition”：exp((-6096.9385*pow($Temperature,-1))+(21.2409642)-(2.711193e-2*$Temperature)+(1.673952e-5*pow($Temperature,2))+2.433502*log($Temperature))※饱和水蒸汽压力[Pa]=SONNTANG公式（以温度为变量的函数）⑦改名为“P_SAT”⑧定义函数3.各指标在STAR-CCM+中的设定⑥新的FieldFunction作成-24--相对湿度：RelativeHumidity(RH)9.Tool＞FieldFunction＞New创建新的“场函数”--“相对湿度”10.右键点击UserFieldFunction1＞Rename…更名为“RH(%)”11.在Properties栏中，设定“FunctionName”、“Dimensions”和“Definition”“FunctionName”：RH“Dimensions”：无维数“Definition”：100*$PH2O/$PSAT⑨新的FieldFunction作成⑩改名为“RH(%)”⑪定义函数3.各指标在STAR-CCM+中的设定-25--平均空气龄：MeanAgeofAir（MAA）1.Continua＞Physics1＞SelectModels…选择“PassiveScalar”2.Region＞PhysicsCondition＞PassiveScalarSourceOption在属性栏中钩选被动标量源3.Region＞PhysicsValues＞PassiveScalarSource下属性栏中Method选择FieldFunction4.Region＞PhysicsValues＞PassiveScalarSource＞FieldFunction下属性栏中ScalarFunction选择Density①选择PassiveScalar②钩选PassiveScalarSource③选择FieldFunction④ScalarFunction选择Density3.各指标在STAR-CCM+中的设定-26--平均空气龄：MeanAgeofAir（MAA）5.Tool＞FieldFunction＞New创建新的“场函数”--“平均空气龄”6.右键点击UserFieldFunction1＞Rename…更名为“MeanAgeofAir”7.在Properties栏中，设定“FunctionName”、“Dimensions”和“Definition”“FunctionN