基于可压缩SIMPLE算法的叶栅通道湍流流场的数值模拟
基于Fluent内燃机压缩过程数值模拟-孙宝庆
(a) smoothing 参数设置
(b) remeshing 参数设置
图 4 动网格参数设置
图 5 活塞运动参数设置
(a)5oCA 网格
(b)60oCA 网格
(c)120oCA 网格 图 6 不同曲轴转角下动网格模型
(d)180oCA 网格
4
模拟结果分析
4.1 速度场分析
如图 7 所示，本文给出了，在不同曲轴转角下，速度场云图。由图 7 可知，在不同曲 轴转角下速度场云图变化较大，在 0-180oCA 进气过程中，活塞由上止点到下止点运动 140mm,速度最大值变化较大，在 60oCA 时，最大速度为 237ms/s,在 120oCA 时最大速度为 148m/s，而到 180oCA 最大速度为 55.8m/s。由此可得在进气过程中，速度由零逐渐增大， 而后逐渐降低，在进气末尾速度逐渐降低。且由云图可知，在运动过程中存在着漩涡。这 更加有利有在发动机后期燃烧过程中，燃气混合。
[参考文献] [1] 杨忠敏.车用动力柴油化的技术研究与开发方向[J].柴油机设计与制造，2004,(01) [2] 王东,黄震,张有等.螺旋进气道三维流场数值模拟[J].计算机辅助工程， 2007,(09) [3] 曹永悌. 基于 Fluent 对康明斯 ISM370 柴油掺氢发动机燃烧过程的数值模拟[C].郑州： 华北水利水电大学， 2012 [4] 张国庆. 内燃机进排气过程数值模拟[C].太原：中北大学,2013
1
前言
气缸是发动机压缩过程中最主要的工作部件，对缸内气流流动数值模拟也显的格外重 要，同时由于目前随着柴油机转速逐渐提高，柴油机的动力、经济性指标和环保指标的提 高也越来越受到人们的重视。 而这些性能指标的提高在很大程度上受限于燃料的燃烧程度。 气流动状态则对柴油机燃烧室内的混合气生成的质量有着很大的影响。进气道的形状对柴 油机的进气性能有着很大影响，进而也影响着整机的性能 随着有限元技术的成熟，在包括汽车发动机在内的几乎所有工程领域得到愈来愈广泛 的运用 [1] 。有限元技术在汽车发动机零部件设计上的应用也越发成熟。计算流体动力学的 应用可以提高直喷柴油机进气道设计的整体性能，探求气道内每一点气体流动的变化，使 我们了解相应流场的微观结构。同时可以缩短研发周期，大大推动柴油机工业的发 [2]。
压气机跨音速叶栅叶型的数值模拟
求解器采用耦合求解器隐式方案 ， 根据湍流模型
收 稿 日期 ：０ ８— ６—０ ２０ ０ ３ 作 者 简介 ： 张广 （９ ３一 ， ， 龙 江 哈 尔滨 人 ， 读 硕 士 ， １８ ）男 黑 在 主要 研
对弯度． 最大相对厚度 ｅ不变的基础上 ， 厂 ＿ 、 将最大 弯度 位置 、 最大 厚度 位 置 向后 缘移 动 ， 过 这 通
摘
要： 以抑制叶型吸力面激波强度 、 附面层分离 ， 降低叶型损失为 目标 ， 对某型 压气机跨 音速 叶
栅 叶型进行气 动优化 。运用数值模拟 的方法对 比研究优化后 的叶栅 同原型叶栅在相 同工况下 的 气动性能 ， 概述了气动性能产生变化 的机 理。与优化 前对 比发现 ， 优化后 的叶 型具 有低损 失 ， Ｔ
采用 迎风 格式 的有 限体积 法对 二维 定 常 Ⅳ 一 ｓ方程 进行 数 值 求 解 ， 维 流 场 计 算 区域 采 用 Ｈ 二
型 网格 ， 网格 数 为 ９ （ 距 方 向 ） ×３ ０ 流 动 方 ０栅 ８（
２ 气 动 优 化 思 路
借 鉴高 亚音 速叶栅 以及 超音 速 叶栅 叶型特点 对 该跨 音速 叶 型进 行 优 化设 计 ， 索 验 证 叶 型最 探
马赫数 。叶栅 前后 延伸 部分 均采用 周期 性边 界条 件 。叶盆和 叶背均 采用 绝热 的无 滑移壁 面条 件 。
研究 方 向。本 文通 过 求 解 可 压缩 雷 诺 平 均 Ｎ—ｓ
方程 ， 数值 模拟 了压 气机 叶栅跨 音速 流 动情况 。
１ 计 算方 法 与边 界 条 件
作 范 围大 的特 点 。
关键 词 ： 跨音速压气机 ； 叶型 ； 气动优化 ； 数值计算
（完整版）渐缩渐扩管的无粘与可压缩流动模拟（fluent）
（完整版）渐缩渐扩管的⽆粘与可压缩流动模拟（fluent）渐缩渐扩管的流动是计算流体⼒学模拟的经典问题之⼀。
在这类流动中，激波的出现是流动中可压缩效应的体现。
精确的激波模拟是CFD 研究的热点之⼀。
为了更好捕捉压⼒梯度，需要采⽤较细的⽹格并结合合适的数值模拟和格式。
很多实际模拟中，局部⽹格的⾃适应会很有帮助。
⼀、Double Precision 的选择如果⼏何模型包含多个通过⼩直径管道相互连接的体，⽽某⼀个区域的压⼒特别⼤（因为⽤户只能设定⼀个总体的参考压⼒位置），此时，双精度求解器可能更能体现压差带来的流动（如渐缩渐扩管的⽆粘与可压缩流动模拟）。
⼆、操作条件对于可压缩流动的计算，推荐在Define/Operating Conditions 中将⼯作压⼒设为0，从⽽，最⼩化由于压⼒脉动⽽引起的误差。
三、压强进⼝边界条件压强进⼝边界条件⽤于定义进⼝流体的压强，可⽤于不可压缩和可压缩流动。
当进⼝压强已知，⽽流动速度或流量未知时，可使⽤压强进⼝边界条件。
压强进⼝边界条件也可⽤于定义外部或⾮受限流动的“⾃由边界”。
定义总压与静压⾸先道总压（p ）与静压（0p ）的关系如下（根据伯努⼒积分）：2012p p v ρ=+在Momentum （动量）选项卡内，Reference Frame 为参考值，有绝对值（Absolute ）与相对于临近区域值（Relative to Adjacent Zone ）两个选项供选择，⼀般保持默认的绝对值，Gauge Total Pressure （表总压）⽂本框中输⼊总压的值。
静压在FLUENT 中被称为Supersonic/Initial Gauge Pressure （超⾳速/初始表压），如果进⼝流动是超⾳速的或者是准备压强进⼝边界条件进⾏计算的初始化⼯作，则必须定义静压。
在流场为亚⾳速时，FLUENT 将忽略Supersonic/Initial Gauge Pressure （超⾳速/初始表压）的输⼊数据，⽽⽤驻点参数求出静压。
直通式迷宫密封内可压缩流场的CFD数值模拟_
3 二维迷宫密封数值模拟
密封模型如图 1所示 。
图 1 密封模型示意
数值模拟模型中 [1 ] , L = B = 3. 21mm, C = 0. 2 mm, R = 75 mm,齿厚 t = 0. 2mm; 出口静压为一个 大气压 1. 01325 ×105 Pa。入口参数 :M a = 0. 65, R
齿间距 L 对流量系数的影响如图 7所示 。齿 间距 L 对流量系数的影响较弱 ,随着 L 的增加 ,流 量系数线性递减 。
图 4 速度矢量云图
分析密封齿单个几何参数对泄漏量的影响 , 其他参数及边界条件保持不变 。齿隙 C 和齿厚 t 对流量系数的影响如图 5、6所示 。
图 5 齿隙 C对流量系数的影响
μ t
———湍流粘度系数
,μt
=ρCμ K2
/ε
其中 C1ε = 1. 44; C2ε = 1. 92; Cμ = 0. 09;
σ k
= 1.
0;σε
= 1.
3
近壁面处理采用标准壁面函数 。
2. 2 边界条件 以尽可能的接近物理真实为原则 , 进口边界
条件为滞止总压力 ;出口边界条件为静压 ,取值一 个标准大气压 。壁面为无滑移条件 。
综合以上分析 ,对齿型进行改进 ,把原齿齿顶 削去一角 ,形成 45°的齿尖 ,有 A、B 两种齿型可选 (见图 8) 。这样既能保证齿的整体强度 ,又可尽 量增强节流效果 。
前后两齿分别依次采用 A、B 两种齿型 ,对各 种组合进行计算 ,部分结果如表 1所示 。
28 FLU ID MACH INERY Vol136, No106, 2008
·
mi
=μ1
μ
i2
i
离心式空气压缩机叶轮内湍流的数值模拟
离心式空气压缩机叶轮内湍流的数值模拟
龚宝龙;李剑峰;贾秀杰;李方义;王光存
【期刊名称】《石油机械》
【年(卷),期】2013(41)9
【摘要】应用粘性标准两方程湍流模型、三维雷诺时均N-S方程,采用非定常数值模拟方法,研究了带分流叶片离心式空气压缩机叶轮内部的流场分布情况.研究结果表明,压力面的静压大于相应位置吸力面的压力,最大静压位于压力面后缘;最大相对速度位于叶片入口,长叶片吸力面有回流现象,短叶片对遏止回流有明显作用;长叶片压力面前缘压力非定常波动最严重,后缘其次,叶片中部的压力波动最弱.该项研究成果可应用于叶轮结构的优化设计及指导叶轮的再制造.
【总页数】4页(P37-40)
【作者】龚宝龙;李剑峰;贾秀杰;李方义;王光存
【作者单位】山东大学机械工程学院高效洁净机械制造教育部重点实验室;山东大学机械工程学院高效洁净机械制造教育部重点实验室;山东大学机械工程学院高效洁净机械制造教育部重点实验室;山东大学机械工程学院高效洁净机械制造教育部重点实验室;山东大学机械工程学院高效洁净机械制造教育部重点实验室
【正文语种】中文
【中图分类】TE924
【相关文献】
1.离心式叶轮内三维紊流场的数值模拟 [J], 胡志伟;姜培正
2.双流道泵叶轮内湍流的数值模拟 [J], 刘厚林;袁寿其;施卫东;马皓晨;袁建平
3.扭曲叶片离心泵叶轮内三维湍流数值模拟 [J], 付强;朱荣生;李维斌
4.双流道叶轮内湍流的三维数值模拟 [J], 刘厚林;袁寿其;施卫东;孔繁余
5.带分流叶片的离心泵叶轮内三维不可压湍流场的数值模拟 [J], 何有世;袁寿其;郭晓梅;袁建平;丛小青;黄良勇
因版权原因，仅展示原文概要，查看原文内容请购买。
力学流体simple算法完成
力学流体simple算法完成一、简介力学流体（Computational Fluid Dynamics，CFD）是一种通过数值模拟来研究流体运动的方法。
它可以帮助工程师和科学家了解流体如何在各种条件下移动和相互作用，以及如何在不同的环境中优化设计。
其中simple算法是一种常用的求解Navier-Stokes方程的方法，本文将对其进行详细介绍。
二、simple算法基础1. Navier-Stokes方程Navier-Stokes方程是描述流体运动的基本方程之一，它包含了三个主要部分：质量守恒方程、动量守恒方程和能量守恒方程。
其中，质量守恒方程描述了质量在空间和时间上的变化；动量守恒方程描述了力对流体速度的影响；能量守恒方程描述了温度在空间和时间上的变化。
2. 离散化方法为了求解Navier-Stokes方程，需要将其离散化为有限个离散点上的代数形式。
这可以通过有限差分、有限元或谱方法等技术实现。
其中，有限差分法是最常用的离散化方法之一。
3. simple算法原理simple算法是一种求解Navier-Stokes方程的迭代方法，它采用了分步求解的策略。
simple算法包括三个主要步骤：预处理、压力修正和速度修正。
其中，预处理步骤用于确定初始速度场和压力场；压力修正步骤用于计算压力场的变化；速度修正步骤用于更新速度场。
通过迭代这三个步骤，可以逐渐收敛到稳定的解。
三、simple算法实现1. 离散化首先需要将Navier-Stokes方程离散化为有限差分形式。
例如，在二维情况下，可以将动量守恒方程离散化为：$ \frac{u_{i+1,j}^{n+1}-u_{i,j}^{n}}{\Delta t}+\frac{u_{i,j+1}^{n}-u_{i,j-1}^{n}}{2\Delta y}-\frac{\nu}{\Delta x^2}(u_{i+1,j}^{n+1}-2u_{i,j}^{n+1}+u_{i-1,j}^{n+1})=-\frac{p_{i+1,j}-p_{i,j}}{\Delta x}$其中，$u$表示水平方向的速度，$v$表示垂直方向的速度，$p$表示压力，$\nu$表示粘性系数。
叶轮内部三维不可压湍流数值模拟
叶轮内部三维不可压湍流数值模拟
黄道见
【期刊名称】《农机化研究》
【年(卷),期】2005(000)006
【摘要】一般来说,离心泵叶轮内的流动是三维的湍流流动,叶轮的旋转和表面曲率效应以及随之而来的哥氏力和离心力,使叶轮内的流动极其复杂,致使内部流场测试困难.随着计算机技术的迅速发展,叶轮内流数值模拟研究相当活跃.为此,将计算流体力学(CFD)技术应用于叶轮设计,基于Navier-Stokes方程和标准κ-ε紊流模型,依据三维数值模拟的结果,优化与叶轮设计相关的几何参数,使叶轮内的流态接近于理想流态,从而保证叶轮具有良好的性能.
【总页数】3页(P203-204,207)
【作者】黄道见
【作者单位】江苏大学,流体工程技术中心,江苏,镇江,212013
【正文语种】中文
【中图分类】TH302
【相关文献】
1.低脉冲浆泵叶轮的三维造型及其内部的湍流数值模拟研究 [J], 胡庆喜;黄运贤;黄放辉;陈中豪
2.低比速离心泵叶轮内三维不可压湍流场计算 [J], 陈池;袁寿其;郑铭
3.双流道污水泵叶轮内部三维湍流流动的数值模拟 [J], 张静;齐学义;侯华;王建森
4.带分流叶片的离心泵叶轮内三维不可压湍流场的数值模拟 [J], 何有世;袁寿其;郭
晓梅;袁建平;丛小青;黄良勇
5.无过载离心泵叶轮内三维不可压湍流场计算 [J], 袁寿其;陈池;郑铭;李昳
因版权原因，仅展示原文概要，查看原文内容请购买。
【国家自然科学基金】_叶栅通道_基金支持热词逐年推荐_【万方软件创新助手】_20140730
推荐指数 3 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2Hale Waihona Puke 10年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
科研热词 气冷涡轮 黏性连续伴随方程 高负荷扩压叶栅 约束 燃气涡轮 熵增 激波 湍流模型 气膜冷却 气热耦合 气流角 气动优化 数值模拟 控制理论 换热系数准则公式 局部边界层吸除 实验研究 吸气量 吸气位置 参数化建模 冲压叶栅 冲击冷却
2011年 科研热词 数值模拟 附面层 激波 抽吸流 压气机 高负荷压气机叶栅 频谱分析 非轴对称端壁 非定常流动特性 遗传算法 转差 跨声速叶栅 航空航天推进系统 脱落频率 端壁运动 相干效应 燃气涡轮 湍流模型 涡轮设计 涡轮 流场结构 气热耦合 换热系数准则式 损失系数 总压损失 尾迹涡 对转技术 叶顶间隙 叶片造型 叶栅吹风实验 叶栅 可控涡设计 压力可控涡 优化设计 二次流 piv技术 推荐指数 5 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
推荐指数 3 3 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2013年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43
基于Fluent流场数值仿真的管路流量计算_张功晖
Hydraulics Pneumatics&Seals/No.12.2010基于Fluent流场数值仿真的管路流量计算张功晖1黎志航2周志鸿1(1.北京科技大学土木与环境工程学院，北京100083; 2.广东肇庆爱龙威机电有限公司，广东肇庆526238)摘要：利用Fluent三维单精度求解器，对管路内的三维稳态流场进行仿真，利用后处理工具得到管路体积流量，并将Fluent数值仿真计算的体积流量结果与实测结果进行对比，数值仿真计算结果得到实际测量实验的验证。
关键词：Fluent；管路；流量中图分类号：TH138.52文献标识码：A文章编号：1008-0813（2010）12-0041-03Air-passage Structure Improving of Pneumatic ElectromagneticValve Based on Flow Field Simulation withing FluentZHANG Gong-hui1LI Zhi-hang2ZHOU Zhi-hong1(1.Civil&Environment Engineering School of University of Science and Technology Beijing,Beijing100083,China; 2.Guangdong Zhaoqing L&V Co.,Ltd.,Zhaoqing526238,China)Abstract:This thesis applies Fluent single-precision solver calculate the volumetric flow rate by simulating3D steady flow field of the pipeline,and compares the calculated flow rate and the actual measured result.Key Words:fluent；pipeline；volumetric flow rate0提出问题广东肇庆爱龙威公司构建了如图1所示的管路，管路由一段长为L1=500mm、管内径为D1=4mm的塑料管AB，与一个长度为L2=40.14mm、孔径为D1=1.25mm 的不锈钢零件BC连接而成。
机翼绕流可压缩流动的模拟实验报告
机翼绕流可压缩流动的模拟实验报告一、实验目的学会使用Fluent进行机翼扰流可压缩流动的设置和求解方法二、实验步骤1.读入网格文件，检查网格质量，确认最小体积为正值，并进行网格缩放。
2.设置求解器参数，保持默认参数设置，激活Spalart-Allmaras模型，选择Spalart-Allmaras模型以模拟粘性单方程。
3.定义材料属性,选择ideal-gas作为流体，粘度选择sutherland。
4.设置操作条件，设置工作压力为零，以最小化压力脉动引起的误差。
5.设置边界条件，设置压力远场边界条件，包括表压101325Pa、马赫数0.4066、湍流粘度比10等。
6.求解与后处理，初始化流动，激活残差图，设置求解控制参数，开始求解，直到达到收敛。
7.设置参考值，为升力、阻力和力矩的参考值，用于无量纲化；检查升力系数，设置求解控制参数，指定收敛标准，设置升力系数与阻力系数的监视器，再次迭代求解。
8.检查翼型附近Y+的值，通过显示Y+的分布来检查网格，确保壁面附近的网格满足Y+值要求。
9.自适应网格，如果Y+值不在推荐范围内，使用自适应命令细化壁面附近的网格。
10.后处理：检查Y+值，显示速度矢量和压力分布，保存工况和数据文件。
三、结果分析1.实验过程中得到的云图/矢量图/流线图/监控变量的图…图1计算域网格显示图 2 残差值图 3 升力系数值图4 阻力系数图图 5 壁面Y+图 6 局部网格展示2.练习与讨论中获得的结果的各种图…(1)图1显示了一个流体动力学模拟的网格图，绿色线框代表了计算域的网格结构，在翼型前后方以及上下做了局部加密。
本次实验模拟的是亚声速流场，流场范围设定为弦长的50-100倍。
(2)图2是一个残差图，显示了在迭代过程中不同物理量（如连续性、x-速度、y-速度、能量和湍值）的残差变化。
这些残差随着迭代次数的增加而减小，模拟正在逐渐收敛。
(3)图3显示了升力系数（Cl）随迭代次数的变化，在初始迭代后迅速下降并趋于稳定。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
关键词 ： 流动控 制方 程 ； 差分格式 ； 网格 分布 ； 数值模拟 分类号 ：Ｈ ５ Ｔ ４ 文献标识码 ： Ａ 文章编号 ： ０ —８ ４ ２０ ）４０ ６ －３ １ １５８ （０ ７ ０ － ６０ ０ ２
Ｎｕ ｒ ａ ｎ ｅｔ ａｉｎｏ ｕｂ ｌｎ ｌｗ ｆｒＣ ｓａ ｅＰ ｓａｅＢ ｓｄ ｍｅ ｃｌＩｖ ｓ ｇｔ ｆＴ ｒｕｅ ｔＦｏ ｏ ａｃ ｄ ａｓｇ ａｅ ｉ ｉ ｏ ｏ ｏ ｒｓｉｌ Ｉ Ｌ ｒｔｍｅｉ ｎ Ｃ ｍｐｅｓｂｅＳＭＰ Ｅ Ａ ｈ ｔ ｉ ｃ
周红梅 ， 苏莫 明， 任 雄
（ 西北 工业 大学 ， 西安 ７ｏ ７ ） １０ ２
摘要 ： 开发 了二维任 意曲线坐标系下 、 采用 同位 网格 布局 、 能计 算跨音速 可压缩流动 的 ＳＭＰ Ｅ方法 计算程 序。为 Ｉ Ｌ 提高精度 ， 采用 Ｑ Ｉ Ｋ格式 、 Ｕ 格式在不同 的网格分布下计算 。为 了验证本 文方法对 低速 到跨音速 范围 内的粘 ＵＣ ＣＩ 性流动的数值分析效果 ， 对几个典型算例进行了数值试验 ， 计算结果 与文献 和实验数据吻合较好。
ａ ＋
∞
ａ ｖ
一 一
ａ ∞
＋ｐＪ ｏ（
’
一２ ） ｖ
ａ ＋
ｘ
＋ｐｏ ｏ ＋２ ） ￣（Ｊ ｃ ｕ
Ｐ—ｐ ８
１ 控制方程与 湍流模 型
将笛 卡尔坐标系（ Ｙ 转 化为任 意 曲线坐 标系 （ ｒ ， Ｘ， ） ，／ ）
采用标 准湍 流模 型 ， 则连续方程 ， 动量方 程 ， 湍动能方程 和
Ｚ Ｕ ｎ － ｉＳ － ｎ ， Ｎ ｘｏ ｇ ＨＯ Ｈｏ ｇｍｅ ， Ｕ Ｍｏｍｉｇ ＲＥ ｉｎ
（ ｏ ｈ ｅｔｒ ｏ ｔ ｈ ｉ ｌＵ ｉｅｓｙ Ｘ ’ｎ７ ７ ， ｈｎ ） Ｎ ｒ ｗ ｓ ｎＰ ｌ ｅ ｎ ａ ｎｖｒｉ ， ｉａ １ ０ ２ Ｃ ｉａ ｔ ｅ ｙｃ ｃ ｔ ０
和不可压缩 流来 建立 和发 展差分 格式 和计算 方法 的。在计
Ｊ
ｑ ＝ 。 ， ＋ “ ，
＝
ｙ ‘， ｙ
算 流体力学 三种 主要 的对 控制 方程 的离 散方 式 中 ， ＳＭ— 以 Ｉ Ｐ 方 法为代表的有限体积法近些 年被广泛应用 ， Ｅ Ｌ 它使在采 用 分离式求解方法 时各类 变量 能同步 的加 以改进 ， 从而提高 收敛速 度 ， 因此 已成为一种可 以计算具有 任何流速 流动 的数
ｌ ｎｅ ｂｔ ｒ ｕｖ ｏｒｉａ ｙｔ ｆｗ ｕ ｄｒ ｒｉａｙｃ ｒｅｃｏｄｎｔ ｓｓ ｍ．Ｆｒｉ ｒｖ ｇｃｍ ｕ ｔ ｎｌｐｅｉｏ ，ｔｅｈｇ ｅｏ ｔ ｎ ｓｈｍ Ｕ Ｃ ｏ ａ ｒ ｅ ｅ ｏ ｏｉ ｏ ｐ ｔ ｉ ａ ｒ ｓ ｎ ｈ ｉ ｒｓｌ ｉ ｃｅ ｅＱ ＩＫ ｍｐ ｎ ａｏ ｃ ｉ ｈ ｕｏ
维普资讯
第４ ９卷 第 ４期
２０ ０ ７年 ８月
汽
轮
机
技
术
Ｖｏ ． ９ Ｎｏ ４ １４ ．
Ｔ ＲＢ Ｎ Ｅ Ｕ Ｉ Ｅ ｒ ＣＨＮＯＬ ＯＧＹ
Ａｕ ． ０ ７ ｇ２ ０
基 于 可 压缩 ＳＭ Ｌ Ｉ Ｐ Ｅ算 法 的 叶栅 通 道 湍 流 流场 的数 值模 拟
Ａｂ ｔ ａ ｔ Ａ ａｃ ｌｔ ｎ ｒ ｇ ａ ｂ ｓ ｄ ｏ Ｉ Ｌ ｒ ｈ ｔ ｓ ｄ ｖ ｌｐ ｄ ｔ ａｃ ｌｔ ｏ ｒｓ ｉ ｌ ｎ ｎ ｏ ｒ ｓ ｉｌ ｓｒ ｃ ： ｃ ｕ ａｉ ａ ｐ ｏ ｒｍ ａ ｅ ｎ Ｓ ＭＰ ｌ ｏｌ Ｅ ａ ｔ ｍｅｉ ｉ ｅ ｅ ｏ ｅ ｏ ｃ ｌｕ ａ ｅ ｃ ｍｐ ｅ ｓ ｅ ａ ｄ ｉ ｃ ｍｐ ｅ ｓｂ ｅ ｉ ｃ ｂ
逆变分 量 （ ， ） Ｊｃｂａ 、 、ａｏｉ Ｊ 变换 系 数 ｑ 扩 散 系数 ｎ
０ 前
言
， 分别 由下式 给出 ：
＝
Ｊ 蟛 ＋ （
ｙ 一
．
）
随着计算机技术 突飞猛进 的发 展 ， 算流体力 学在流动 计 领域 的应用 日益 广泛 。一般 而言 ， 数值计 算是针对 可压缩流
值方法。本 文将 ＳＭＰ 算法 应用 于可 压缩 流场 ， 双圆弧 Ｉ Ｌ Ｅ 对 平面叶栅通道 内部流 场流动进行 数值计算 的
目的 。
Ｓ ＝
［ ｔ １
０
，＋＂ Ｏ ｋ
方 程 中 的源 项 可用 向 量 表 示 为