摘要： 本文采用数值模拟方法研究二维前台阶流动问题。着重研究了层流和湍流两种情况下前台阶分离流的不同现象，给出了雷诺数为100时的层流流动结果，以及湍流情形下（re=10000）的dns(直接数值模拟)分析和标准k-e模型的模拟结果，并将计算结果与dns进行比较分析。
关键词： 前台阶 dns k-e模型数值模拟。
近年来，人们在研究数值求解n-s方程方面取得了很大进展。由于计算机技术的飞速发展，相对于已经发展得较成熟的雷诺平均模型（rans）,众多研究者倾向于直接求解数值模拟（即dns）。本文采用dns和标准k-e模型两种方法模拟研究二维前台阶湍流流动问题，以及在雷诺数相对较小的层流问题。
二维前台阶分离流的数值研究已经有了很大进展，文章采用不可压n-s方程对这一问题进行研究。
流体掠过一个直角表面的二维流场，在工程实践中有重要意义。图1是气流掠过前台阶二维流场的模型示意图。整个流场大致分为五个区域：
主流，初始边界层，边界层分离区，台阶前部回流区和台阶后部回流区。边界层从主边开始形成和发展，在台阶以前开始分离，在下底面与迎流面之间形成第一个回流区，气流掠过台阶后在上部平面又形成第二个回流区。与后台阶流场相比，更为复杂。
图1 气流掠过前台阶的二维流场。
1.主流；2.初始边界层； 3.边界层分离区；
4.台阶前部回流区； 5.台阶后部回流区)
本文计算的流动模型为经典的前台阶流动。由于其相对简单的外形允许采用较密的网格。为了简化物理问题以及为与dns结果比较，计算采取改变不同的运动粘性系数从而达到改变雷诺数大小的目的，台阶高度为0.
4m（h为特征长度），来流温度为常温t=293，来流速度u=1m/s，压强为常压p=0.1mpa。因此开展了以下三个算例：
1） 层流情形：re=100，
2） 湍流dns：re=10000，
3） 湍流k-e模型与（2）条件相同。
在计算中，为了使得前缘，物面及台阶角点附近有足够多的网格点数，且便于边界处理，本文在前缘、台阶附近有足够多的均匀密网格。
图2为层流情形下的速度分布图。
图3为层流情形下的压力分布图。
图4为湍流dns情形下的速度分布图。
图5为湍流dns情形下的压力分布图。
图6为湍流k-e情形下的速度分布图。
图7为湍流k-e情形下的压力分布图。
由上述三个算例的压力、速度分布图可知，在层流情形下，压力的变化在流场中不是很剧烈，仅在台阶处产生较小的压力梯度；从湍流dns和k-e模型比较看，dns结果明显的更能够反映湍流流场中的脉动特性，而k-e模型是对流场中的速度进行平均后的结果，与实际情况仍然存在一定的误差；但从压力场看，层流结果与k-e模型模拟结果更符合实际情况。
图8为三种情形**阶后壁面上的压力分布图。
a） 层流情形下。
b） 湍流dns
c） 标准k-e模型。
由壁面上的压力分布图可知，在台阶附近压力均出现剧烈变化，尤其是层流和湍流k-e模型模拟得到的均是压力急剧减小，然后随着流动压力又逐渐回升，说明在台阶前缘附近产生了强烈的涡旋结构。而dns结果与这两种情形相差较大，在台阶的壁面处整个压力分布均处在低压状态，说明dns计算得到了明显的边界层分离区，这是与实验观察相吻合的，更符合实际情况。
对于前台阶分离流动，在流动速度很高的情况下，流体可压缩性对流动的影响很大，台阶角点附近，这一影响尤为明显。用不可压 n-s方程来模拟此流动是远远不够的层流。就目前的研究情况来看，用可压缩n-s方程对前台阶分离流数值模拟的文献很少，这方面工作显得很薄弱，也是需要进一步研究的工作。
同时流体粘性和台阶的前缘效应等也是函待解决的问题。
fluent大作业
1 确定求解器。选择用于进行cfd计算的求解器，本例选择fluent5 62 创建圆环。在gambit中创建一个圆环直径为1，外圆半径为10的圆环 3 创建立方体 4 移动立方体。5 分割圆环。6 删除3 4圆环和立方体。7 建立弯管直段。8 移动弯管直段。9 整合弯管和直段。10 边界层的设定。1...
Fluent结课作业
结课作业。学生姓名 王一帆。学号 2201300293 专业 动力工程。题目 层流无热传导槽道问题的模拟。授课教师 郑建祥。日期 2014年5月14日。本题描述层流无热传导槽道问题的fluent模拟，即流动为层流，无热传导，不需考虑其他特殊因素。1.启动gambit软件。2.创建面。3.网格划分。4...
fluent笔记
discretization离散 node values节点值，coarsen粗糙 refine细化 curvature曲率，x wall shear stress 壁面切应力的x方向。strain rate应变率。1 求解器 solver 分为分离方式 segeragated 和耦合方式 coup...