流动部件的建模.ppt

流动部件的建模 CDAJ-CHINA 主要内容 流动基本方程和解法 时间步长的计算 离散 普通管路的建模 分叉管的建模 连接部件的建模 边界条件的设定 流动的基本方程 连续方程： 能量方程： 动量方程： 流动的基本方程 求解的难点往往不是微分方程本身和编制程序本身而是难于获知解题所需的各种实用和合理的物理常数及参数值。而这些往往是各个公司的技术秘密。 针对具体的对象，通过大量的试验和产品的使用经验而获得的。有的时候，由于对于研究对象的物理规律还不是很清楚，或者情况过于复杂，就要利用经验公式和中间试验结果协助进行计算。 流动方程的解法 隐式： 计算结果依靠所有的网格和边界的条件建立非线性微分方程联立求解（迭代和插值）。一般采用比较多的计算步数。但是计算的稳定性要好于显式算法。所以这种算法用于比较长的时间的模拟。 显式：新的时间步长的结果只与控制体本身和相邻的节点有关，直接由前一个步长的结果求出不需要迭代。 时间步长的计算 显示算法： 这种方法的时间步长必须满足Courant条件： Courant数 =1，一般我们取0.8 m：系数（def＝1） 时间步长由GT-SUITE本身决定。 这种方法对于进排气系统和燃油喷射系统内的压力波动计算比较准确。而冷却系统中一般不考虑压力波动。 GT-FULE和大多数GT-POWER模型用这种算法。 时间步长的计算 隐式算法： GT-COOL和GT-POWER中暖机工况的排气系统计算用这种方法。这种方法计算速度比较快，但是也必须满足一些条件： 系统内参数波动比较小 （或者在计算中这种波动不重要) 系统内最大的Mach数小于0.3 （计算出的流量在节流出会高于实际流量） 时间步长由用户直接输入。如果时间步长过大，程序会自动警告（或者在*msg *out文件中显示) 在ImplicitSolnControl模块有比较详细的说明。 离散长度 推荐：进气系统： 离散长度＝0.40BORE 排气系统： 离散长度＝0.55BORE 在需要高频分析的时候要把离散长度变短。 例如: 声学分析和性能分析@6000rpm+ 声学分析的时候的离散长度应为性能分析的一半。 系统中不能有一到两个长度很短（远小于离散长度）的部件，这样会导致时间步长很小，计算时间过长。一般来讲可以把这些长度很短的部件合并到连接件中来考虑。 限制时间步长的部件可以在GT-POST的Time Step Restriction中看到，或者在* out文件中看到。 Pipe 摩擦损失： 在Run setup 中的Flow control可以通过定义 Global Friction Multiplier (constant term) 增加整个流动系统的摩擦。或者通过定义 Global Friction Multiplier (unsteady term) 研究流体流动脉冲对摩擦损失的影响。 Pipe 压力损失：主要是由于管路的弯曲、截面的变化引起的。程序中可以通过定义压力损失系数Cp来考虑。这里的压力损失系数不考虑壁面的摩擦。如果设置为Def,程序将对其进行自动的计算（以大量的试验数据为依据）。另外,这里的压力损失系数还不包括管路截面形状突变带来的损失，这些损失GT-SUITE也将自动计算。 传热 对于光滑表面的传热系数采用Colburn公式： 对于粗糙表面的传热系数需要进行修正： (因为表面粗糙度对传热影响很大） Pipe 传热系数 流体温度 壁面温度 总的散热量 壁面温度 Heat Conduction Object使用‘HeatCComps’ 用户输入：自然对流、强制对流以及辐射的相关参数。对于自然对流的空气其换热系数一般为5－20W/m^2*K。最大值通常在有风扇的情况下的值（但是不考虑风扇直接吹向表面的那部分）热辐射系数通常设置为0.5－0.8。发光物质的热辐射系数更小。 Pipe 温度可以进行稳态/瞬态的求解，但是在暖机过程中只能采用瞬态求解。在ThermalControl中控制。 壁面温度的求解可以进行开关控制，操作同上。 可以进行绝缘管和同心管的建模。在HeatCComps可以进行多层管的建模。第一行是最里面的一层，最后一行是最外面的一层。而且每层都可以采用不同的材料。（‘HeatCMaterials’） ‘HeatCFlange’可以进行从一个部件表面到另一个部件表面传热的建模（与流动平行）。特别是当两个部件间有较大的温度梯度的时候（比如：从进气歧管到缸头