[2018年最新整理]PHOENICS软件培训.pptVIP

PHOENICS软件培训 PHOENICS软件的变量 因变量（dependent variables） 100个 自变量（independent variables）X Y Z t 辅助变量（auxiliary variables） RHO1, ENUL PHOENICS软件求解的方程： 其离散式为： aP?P=aN?N + aS?S + aE?E + aW?W + aH?H + aL?L + aT?T + sources 其中: aP = aN + aS + aE + aW + aH + aL + aT（by continuity) PHOENICS软件求解方法 1. Point by point： aPfPnew = aNfNold + aSfSold + aEfEold + aWfWold + aHfHold + aLfLold + aTfT + sources 收敛迟缓，适用于各变量强烈偶合 PHOENICS软件求解方法 2. Slabwise： aPfPnew - aNfNnew - aSfSnew - aEfEnew – aWfWnew = aHfHold + aLfLold + aTfT + sources Solve命令的默认解法，对抛物流有较大优越性 PHOENICS软件求解方法 3. Whole－field： PHOENICS中的边界条件（源项） PHOENICS默认(不设任何边界条件)： 所有边界对流体不产生流入（出）、无摩擦、绝热，表现出对称面或轴的特征；phoenics设置边界条件以源项的形式加入。 PHOENICS中的边界条件（源项） 传统源项，如压力梯度、粘性发热等已经在phoenics中内置，term命令控制； 有限容积法离散方程，源项线性表达为：Sbc=TC(V-?) ; T是几何因子,如面积或体积 对每个单元有： PHOENICS中的边界条件（源项） 每单元源项： PHOENICS中的边界条件（源项） PIL施加源项： where and when？ T、C、V？ Patch(name, type, ixf, ixl, iyf, iyl, izf, izl, itf, itl) Coval(name, variable, coefficient1, value1) Coval(name, variable, coefficient2, value2) variable是phoenics求解变量 PHOENICS中的边界条件（源项） 固定值边界：C赋为大数（fixval），V取边界值 PHOENICS中的边界条件（源项） 固定流率边界：C赋为小值，V取边界流率与C的商（C为fixflu，V取边界流率） PHOENICS中的边界条件（源项） 线性边界，如固定换热系数h(W/m2/K)与固定温度(K)的环境换热： Q = Sh(Tsur-Tp) 显然，T＝S，h＝C，V＝Tsur Patch(heatloss, area, ixf, ixl, iyf, iyl, izf, izl, itf, itl) Coval(heatloss, tem1, h, Tsur) PHOENICS中的边界条件（源项） 层流wall边界：静止壁面对流体的层流剪应力 可以令 用wall型边界： 则只需：C＝1.0 PHOENICS中的边界条件（源项） 流入边界：所有质量流边界都是线性引入连续方程，压力为变量 流入边界如果与内部压力无关，则Cm取为fixflu，Vm取所加质量流 固定压力边界：压力固定，与质量流无关，Cm取大值，Vm取边界压力。由于数值计算需要，C并不取fixval，而取1000左右 PHOENICS软件的构成（components） EARTH的组织： GROUND面向用户： Q1文件24GROUP TALK=T;RUN(1,1) GROUP 1. Run title and other preliminaries GROUP 2. Transience; time-step specification GROUP 3. X-direction grid specification GROUP 4. Y-direction grid specification GROUP 5. Z-direction grid specification GROUP 6. Body-fitted coordinates or grid distortion GROUP 7. Variables stored, solved named GROUP 8. Terms (in d