哈尔滨工业大学(威海)本科-终期概要.ppt

学生: 巩向涛 指导老师: 王富强 槽式太阳能热发电系统中内置扰流柱群管式集热器的数值研究 终期答辩 背景和意义 模型介绍 抛物面槽式太阳能集热系统 扰流柱 网格划分及无关性验证 圆管类模型的高质量O型网格 不规则模型的 非结构四面体网格 + 结构网格 与 非结构网格 融合之后再 周期操作 = 网格无关性验证的评判标准：流固耦合面的热通量（整体平均和扰流柱区域平均） 选取170万左右网格数目进行计算 边界条件和udf编译 DEFINE_SPECIFIC_HEAT(fluid_cp, T, Tref, h, yi)/*定义流体的比热容*/ { cell_t c; Thread *t; real Temp = C_T(c,t); return 2014.22+3.86884*(Temp-273)+2.05029*10e-3*(Temp-273)*(Temp-273)- 1.12621*10e-5*(Temp-273)*(Temp-273)*(Temp-273)+3.86282*10e-8*(Temp-273)*(Temp-273)*(Temp-273)*(Temp-273); } DEFINE_PROPERTY(densitydaoreyou,c,t)/*流体密度随温度变化*/ { real T = C_T(c,t); return -0.696982*(T-273)-0.000131384*(T-273)*(T-273)-0.00000209079*(T-273)*(T-273)*(T-273)+776.257; } 物性参数表 拟合公式 C语言代码 边界条件和udf编译 蒙特卡罗光线追踪法 得出热流密度 热流密度（W） 汇聚光照射 q=469687.5×z+45307 y?0, 0≤z≤0.016m q=-1789444.4×z+81453.1 y?0, 0.016m≤z≤0.0196m q=-4030789.5×z+125383.5 y?0, 0.0196m≤z≤0.031m q=429 y?0, 0.031m≤z≤0.035m q=-469687.5×z+45307 y?0, 0≥z≥-0.016m q=1789444.4×z+81453.1 y?0, -0.016m≥z≥-0.0196m q=4030789.5×z+125383.5 y?0, -0.0196m≥z≥-0.031m q=429 y?0, -0.031m≥z≥-0.035m 热流密度（WW） 非汇聚光照射 q=25714.29×z+31.5 y0, -0.035m≤z≤0 q=-25714.29×z+31.5 y0, 0≤z≤0.035m 整理成函数式 DEFINE_PROFILE(lianxi2_heatflux,thread,position) /*高斯热流*/ { real x[ND_ND],y,z; face_t f; begin_f_loop(f,thread) { F_CENTROID(x,f,thread); z=x[2]; if (-0.035=zz=0) F_PROFILE(f,thread,position) = 25714.29*z+31.5; else if (0=zz=0.035) F_PROFILE(f,thread,position) = -25714.29*z+31.5; else F_PROFILE(f,thread,position) = 1000; } end_f_loop(f,thread) } 写进C语言代码 模型对比验证 文章Thermal analysis of solar receiver pipes with superheated steam中的实验数据 出口截面最低温度 出口截面最高温度 实验 云图 模拟 云图 出口壁面温度云图（case6） 数值模拟的算例 光管的研究 光管管壁最大温差随Re变化 光管面平均温度随着坐标变化 结论：光管的温度沿着流动方向升高，在出口处达到最高。 结论：管壁壁面温度分布不均匀和热流密度分布相似，Re越大，温度分布相对越均匀。 光管的研究 Re=11151（左）和Re=19332（右）管内流体湍流动能分布 不同雷诺数下流体平均TKE 结论：通过拟合出的计算式可以看出，雷诺数越大，流体的湍动越大，TKE也越大，Nu上升，换热性能越好，范宁摩擦因子呈现下降。 不同雷诺数下的Nu和f 1979Re9107： 9108Re17286： 内置扰流柱群管的研究 Re=7063时，光管和内置扰流柱群管内TKE分布图 进口速度为0.1m/s时，管内壁面温度分布云图 结