冶金过程数值模拟复习资料参考——LoganZhou详解

1.数学模型、数学模拟、数值模拟、物理模拟数学模型：用数学语言描述实际现象。

数学模拟：通过数学模型对实际现象的描述和求解，实现对实际过程的数学再现。

数值模拟：需要计算机进行数值求解的数学模拟。

物理模拟：是通过实验室物理实验模拟真实物理过程的方法。

将实际地形物理的缩小模型置于实验体（如风洞、水槽等）内，在满足基本相似条件（包括几何、运动、热力、动力和边界条件相似）的基础上，模拟真实过程的主要特征，如空气动力规律和扩散规律。

由于所有相似条件不可能完全满足，针对具体要求恰当选取相似参数是实现物理模拟的关键。

物理模拟主要用于数值计算模式难于处理的复杂地形以及建筑物影响时的扩散研究。

物理模拟实验与现场实验相比条件易控制、重复，且省人力、物力，可进行较全面和规律性实验，是大气扩散研究的重要手段。

2、建立数学模型的方法步骤及意义数学模型的建立可以分成三个步骤：(a). 分解，将复杂的冶金反应过程分解为流体流动、传热、传质和化学反应等基本单元过程；(b). 简化，通过一系列假设对问题进行合理简化，简化原则包括：抓住主要矛盾而不失真实性，满足应用的精度要求，适应当前实验条件，以便进行模型识别和参数估算，适应现有计算机能力；(c). 推导，基于上述分解和假设，开始对所有的数学表达进行具体推导，构建模型的核心方程组，可以是线性、非线性，常微分、偏微分，甚至隐函数；1）、摸排阶段：根据实际生产提出的问题，了解问题背景、主次、本质、目标、建模方法；2）、梳理阶段：确定目标，限定过程现象，据此进行参量分析、确定问题的参数、搜集各种必要信息，明确已知量、未知量、自变量、因变量、主要量、次要量，以及这些量之间的关系和所属的基础理论范畴，同时做出误差允许范围内的简化处理，舍弃一些次要量，便于求解；3）、建模阶段：通过基本过程描述，进行理论依据的选取，根据文献收集，进行合适的数学模型的比较分析和选取，如果没有合适模型，需要自定义模型。

作用和意义：现有工艺：加深对过程的基本现象、反应机理的认识，为改善工艺过程和操作提供依据；探索设备、工艺过程和操作参数的变化对冶金效果的影响和变化规律及它们之间的定量关系，为优化工艺、改进设备、改善操作提供必要的数据和依据；实现对工艺过程的诊断和过程的自动控制；指导中试厂和现场实物实验的设计和规划，以节省开支等。

新工艺：对新设计工艺的可行性和灵活性做出准确的估计；对规划和设计实验室、中试厂或实物规模的实验提供指导；帮助评估中试厂或实物试验结果和进行比例放大；在一定条件下，可替代中试厂或现场实物进行开发性试验，以节省费用等。

3、通量传输速度的确定通量：在空间任意位置上，单位时间内通过垂直于运动方向上单位面积的物理量。

其本身是矢量性质，单位是“物理量单位/(m2·s)”。

层流laminar flow：流速小，流层不混湍流turbulent flow：流速大，流层混合剧烈，微团运动极不规则质量传输：有效扩散系数＝分子扩散系数+涡流扩散系数动量传输：有效粘度＝分子粘度+涡流粘度能量传输：有效导热系数＝分子导热系数+涡流导热系数4.1、控制方程-边界条件及初始条件只有具备足够数量的赋值（或关系式），微分方程才能有特解。

一般非稳态方程要求有一个初始条件。

方程数目根据变量个数确定，而边界条件数目则由方程中变量的导数阶次和个数共同决定，每个n阶导数需要n个边界条件。

一般边界条件取决于局部条件，边界条件的典型类别有：一类边界条件：直接给定边界上因变量的数值。

如研究流体流动时常设流体与固体边界无相互滑移，即固－液界面处u=0。

二类边界条件：边界上存在通量连续条件。

如分析钢锭模向外散热时有：式中，ε——锭模表面发射率；σ——斯芯藩－玻尔兹曼常数；T0——环境温度；Ta——模表面温度。

三类边界条件：直接给定边界传输通量。

如钢包中心线两侧钢液的动能耗散通量为零。

4.2、流体的3种湍流模型求解求μt （μt 为涡流粘度或表观湍流粘度）一般采用的有三种方法：普朗特混合长理论（零方程模型）、κ方程模型（单方程模型）和κ－ε双方程模型，以及后一种方法的变体。

①.混合长模型：对于粘性流体，如果湍流微团从某一层中由于脉动的作用而到达速度不同的另一层，微团运动过程中经历了l m 距离，且这一运动导致目标层的扰动，则随机湍流速度u t 可表达为： 。

其中，为x 方向时均速度在y 方向上的速度梯度的绝对值，l m 为特征混合长度。

进而，普朗特提出涡流粘度可以表达为 该模型被称为零方程模型，因为它是以代数方程表示特征量的。

优点：是不必求解与μt 有关的微分方程，只需要确定混合长（但相当困难） 缺点：该模型仅限于简单流场的描述，复杂流场（如环流）无能为力。

②.κ方程模型：已知μt =ρl m u t ，而湍流流动速度与湍流动能的平方根成正比，即湍流脉动速度u t 有式中κ为湍流脉动动能，于是，涡流粘度可以表达为： 式中C μ为经验常数，κ的数值由湍流能量衡算得到。

例如，在x 方向流动的湍流边界层中，如假设流动属稳态，其湍流能量守恒方程为式中C D ——流量系数，Pr κ——湍流动能的普朗特数。

优缺点：与零方程一样，κ方程模型同样需要首先确定混合长度，这就对研究诸如气体搅拌钢铁等环流过程构成了障碍。

③.κ－ε双方程模型 ：这一模型的基本出发点是以特征能量和特征耗散速率来表示对特征长度的函数关系，即： 式中κ——湍流脉动动能；ε——湍流脉动动能的耗散率。

只要确定湍流脉动动能κ和耗散速率ε，则涡流粘度有解。

由脉动动量方程推导可以得到描述κ和ε的偏微分方程为式中G κ为湍流脉动动能κ的产生速率，C μ、C 1及C 2都是经验常数，Pr κ、Pr ε分别为湍流动能和动能耗散速率的普朗特数。

一般取在κ－ε方程中包括速度项，可见湍流条件下求解速度场需要将连续性方程、运动方程及湍流κ－ε方程联立求解。

步骤：求解运动方程获得速度场→求解k-ε方程获得涡流粘度分布→以涡流粘度为前提，求解速度场温度场→收敛→结束5、数值模拟考虑湍流时，如何理解？ x t m u u l y ∂=∂2x t m t mu l u l yμρρ∂==∂t u ∝22212x y z u u u κ⎡⎤'''=++⎢⎥⎣⎦322Pr t x x yt D mu u u C x y y y y l κμκκκρκρμ⎛⎫⎛⎫⎛⎫∂∂∂∂∂+=-⨯+- ⎪ ⎪ ⎪∂∂∂∂∂⎝⎭⎝⎭⎝⎭322t and m l C μκεμρκε==()eff eff 12 and Pr Pr D D G CG C Dt Dt κκκεμμκεερκρερερεκ⎛⎫⎛⎫=∇⋅∇+-=∇⋅∇+- ⎪ ⎪⎝⎭⎝⎭121.45; 1.92; 0.09; Pr 1.0; Pr 1.3C CC μκε=====x u y ∂∂12t m C l μμρκ=对钢铁冶金而言，湍流特征的描述应当摆在突出位置，因为钢铁冶金过程所涉及的流动问题大多数是湍流问题，一是因为钢液、熔渣等高温流体都是高粘度流体，二是这些流体（包括很多气体）的流动都是高速流动，从它们的粘度和速度计算得到的雷诺数都非常高。

描述这些湍流的目的就是为了求解速度场，确切地说，就是通过给定适当的系数来描述湍流条件下的混合效果，以便进一步利用纳维尔－斯托克斯方程求解湍流速度场。

从模型设计者的角度，就是找到湍流条件下有效传输系数的途径。

6.1、当量直径、比表面积6.2、气—固相反应的数学模型一界面模型假设所有化学反应都发生在内核界面，而这显然是不合适的，实际的Fe2O3的还原过程是一个逐步还原过程，这必然涉及到还原程度不同的分层。

三界面未反应核模型就是认为存在三个反应界面层：Fe-FexO 、FexO-Fe3O4、Fe3O4-Fe2O3。

还原气体CO （或H2）通过边界层扩散到De3层（ Fe-FexO 层）界面，与FexO 反应消耗一部分还原剂，产物CO2（或H2O ）向外扩散，其余还原剂继续往中心扩散，通过De3到De2、De1…整个过程共有11个环节，其中扩散环节8个（DCO,界面层，DCO,Fe 层，DCO,FeO 层，DCO,Fe3O4层，DCO2,Fe3O4层， DCO2,FeO 层， DCO2,Fe 层， DCO2,界面层），界面反应环节3个。

7、气-液双模理论是当气液两相作相对运动时其接触界面两侧分别存在气体边界层气膜和液体边界层液膜。

气膜和液膜均属层流。

双膜理论的基本论点如下：1、相互接触的气、液两相流体间存在着稳定的相界面，界面两侧各有一个很薄的停滞膜，相界面两侧的传质阻力全部集中于这两个停滞膜内，吸收质以分子扩散方式通过此二膜层由气相主体进入液相主体；2、在相界面处，气、液两相瞬间即可达到平衡，界面上没有传质阻力，溶质在界面上两相的的组成存在平衡关系，即所需的传质推动力为零或气、液两相达到平衡。

3、在两个停滞膜以外的气、液两相主体中，由于流体充分湍动，不存在浓度梯度，物质组成均匀。

溶质在每一相中的传质阻力都集中在虚拟的停滞膜内。

8、冶金过程的模拟方法数值方法的计算流程图：偏微分方程的数学分类:9、离散化方法：区域离散化、控制方程离散化区域离散化：1）、内节点法：①先节点后界面；②界面位于相邻两节点的正中间；③代表控制容积的节点位置不一定在控制容积的正中心。

2）、外节点法：①先界面后节点；②节点位置在控制容积的正中心；③界面不一定位于相邻两节点的正中间位置。

注：当采用均匀离散的时候，外节点法＝内节点法【外节点法 & 内节点法之相关说明：1. 边界节点所代表的控制体不同。

外节点法位于非角顶上的边界节点代表了半个控制体；而内节点法则可看成厚度为零的控制体的代表。

2.在非均匀离散中，外节点法的界面是处于节点之间的中间位置，因此所计算的通量值精度较高，但节点位置不处于控制容积的正中心，所以用节点位置的值代替控制容积的值，会有一定误差；内节点法则正好相反。

3.求解区域内如果有物性突变，则内节点法处理这种突变界面比较容易，因为内节点法容易保证控制单元内的物性均一，比如模拟相变过程。

4.离散过程以尽量均匀为宗旨，应注意相邻控制体的厚度变化不宜过大，同一控制体各个方向的尺寸一般不宜相差太远，这些都容易降低计算准确度。

5. 为了书写和编程方便，必须对节点加上编号，原则上编号顺序是任意的，但习惯采用与坐标轴正方向相一致的顺序编号。

如用i表示x方向的节点的编号，j表示y轴方向，k表示z轴方向，因此，节点坐标为(xi,yj)的点可以直接写成(i,j)。

由于i,j和k的值都随坐标轴的前进而增加，所以有下列关系 Xi+△Xi=Xi+1; Xi-△Xi-1=Xi-1;Yj+△Yj=Yj+1;Yj-△Yj-1=Yj-1在节点(xi,yj)上的某传输量的值，如温度，也可以写成T(i,j)或Tij。