CFD经验漫谈+湍流初步

CFD经验漫谈及湍流模化初步CFD的基本流程前处理建模Autocad,CATIA,UG,Solidworks网格生成Gridgen,ICEM-CFD,Gridpro求解器设置计算参数后处理TecplotEnsight网格类型结构网格多块对接网格多块重叠网格非结构网格笛卡尔网格网格划分原则1.第一层网格距壁面的距离层流计算:y(1)~1/sqrt(Re)湍流计算:y+<1壁函数30-2002.网格法向间距拉伸率工程应用:r<1.25~1.3科学研究:r<1.2，ideal:r<1.1拉伸不是无休止的，deltay到边界层厚度/10停止(大约位于边界层中心)3.流向网格间距RANS计算:找感觉。个人经验，Re=10e.6,单位长的壁面附近布置80+点足够GridindependentStudy:粗中细三套网格，比较结果(特别是摩阻)LES计算：理想LES三方向的网格间距应相当法向与RANS相当，流向deltax+,deltaz+~8000(Spalart)，但不得超过边界层厚度/20GridRefinementStudy4.网格品质要求正交性物面网格线垂直于壁面空间网格线正交软件带有网格正交化功能 光滑性疏密过渡要平缓考虑流场特性流动变化剧烈的区域，如激波、分离区域局部加密ExamplesGridindependentStudyy+(1)大小的影响求解器求解参数时间步长隐式求解LUSGS,LU-ADI:CFL=5~10显式求解Runge-Kutta:Dt很小情况变化(Mach提高，高级湍流模型)，收敛困难手工减小时间步长，计算稳定后逐渐提高构造函数，使时间步长自动缓慢增长设置合理初场使用当地时间步技术限制器无激波，无需限制器，否则会引入不必要的粘性出现激波，限制器顺序minmod(高Mach),Smooth,VanAlbada,doubleminmod(低Mach)推荐应用多重网格，极大加速收敛，并促进计算稳定无粘通量空间计算方法RoeAUSM非定常计算Dt<delta/UmaxNaN错误非物理解修正StackOverflow错误数组别都放栈里，建议用动态数组建议学习一个开源代码，希望编程模块化，输入输出接口标准化，参数设置简单化计算环境使用Linux编译器使用intelfortran/C++如有可能，尽量抛弃Windows使用Linux，办公效率提高几十倍！对新生寄予厚望什么是湍流？LeonardodaVinci(1452-1519)在水中放置障碍物后观察流动现象。湍流：极不规则的流动现象。葛饰北斋(1760—1849)：神奈川巨浪歌川広重(1797—1858)：鸣门的风浪艺术家笔下的湍流云层自然界的湍流河流香烟的烟雾海浪减阻控制转捩：减少阻力燃烧强化湍流：加强扩散噪声抑制湍流：降噪气象预测大气湍流湍流的研究意义目的：预测和控制湍流雷诺实验湍流的形成雷诺数：层流湍流O.Reynolds(1842-1912)层流向湍流的发展过程

阶段一：当Re很小时，流体规则地流过圆柱

阶段二：

随着Re增加，尾涡不断增大

阶段三：经过Re阀值后，变为极不规则的运动

雷诺数：湍流的特征

特征一：不规则性射流

特征二：湍流脉动的尺度具有宽广的范围

特征三：拟序结构。剪切湍流存在拟序结构湍流脉动并非完全不规则的随机过程有序大尺度结构包含于不规则脉动中F.N.M.Brown,A.Roshko,1974湍流的分类均匀各向同性湍流剪切湍流：壁湍流：槽道流，圆管流，边界层湍流自由剪切湍流：射流，混合层，远场尾流平均剪切固壁均匀和各向同性平移几何构形平移加转动几何构形均匀：任意n点的空间几何构形在空间平移时，脉动量的n阶相关函数不变。均匀各向同性：任意n点的空间几何构形在空间平移和刚体转动时，脉动量的n阶相关函数不变。均匀：脉动量的任意阶统计相关函数和坐标系的平移无关。均匀各向同性：脉动量的任意阶统计相关函数和坐标系的平移和刚体转动无关。四阶统计相关函数不变均匀各向同性湍流格栅湍流，实验结果实验中采用的格栅均匀各向同性湍流的数值模拟结果壁湍流：槽道流，圆管流，边界层湍流自由剪切湍流：射流，混合层，远场尾流湍流的基本方程湍流的多尺度有结构不规则运动属于宏观流体运动，即使最小Kolmogorov尺度也远远大于分子自由程Reynolds(1886):连续介质假设下仍然满足牛顿定律，即N-S控制方程，并由此导出统计平均值满足的方程，即雷诺平均N-S方程湍流的基本方程Navier-StokesEquationReynoldsAverageO.Reynolds(1895)对上式平均，因为注意到因此雷诺应力本质不是应力，而是湍流对流的平均显现，