IN 738LC材料Chaboche热粘塑性本构模型的隐式Euler格式

为了提高Chaboche模型的数值积分效率,推导了该模型的Eu ler隐式积分格式,针对IN 738LC材料,将该积分格式编写成FORTRAN程序,利用通用有限元分析软件M arc的用户子程序接口HYPELA 2,将该隐式积分格式与有限单元法结合起来。对实验进行了数值模拟,并将隐式算法和显式算法与实验结果作了对比,结果发现隐式算法是合理的,比显式算法收敛性要好。      

喷管化学反应流的对角化点隐式法求解

采用弱耦合对角化点隐式方法的ＭａｃＣｏｒｍａｃｋ格式,求解喷管的粘性化学反应流动,解决了控制方程组的刚性问题,又大大减少了常规点隐式方法中求逆矩阵的大量工作。对液氢液氧火箭发动机,采用６种组分、８个反应有限速率的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流场参数在喷管中的分布。计算结果表明,采用对角化点隐式方法求解化学反应喷管粘性流动能提高计算效率。     

显隐式分析方法在复合材料失效分析中的应用

针对复合材料失效等准静态问题,提出了显隐式联合的求解方法:在加载前期,非线性程度较低时,采用计算效率较高的隐式求解器求解;在材料出现损伤后,非线性程度较高时,改用条件收敛的显式求解器推进计算,中间通过数据传递将材料的应力应变水平及材料损伤状态由隐式求解器传递到显式求解器。该方法在保证计算精度的前提下,兼顾了隐式与显式分析的优点,有效提高了计算效率。     

双方程k-ωSST湍流模型的显式耦合求解及其在叶轮机械中的应用

双方程k-ω剪切应力输运(SST)湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax(BL)模型、Spalart-Allmaras(SA)模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。     

双方程 k-ω SST 湍流模型的显式耦合求解及其在叶轮机械中的应用简

 双方程k-ω剪切应力输运（SST）湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax（BL）模型、Spalart-Allmaras（SA）模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。     

用守恒型差分方法计算底部分离流场

本文用守恒型显式格式及交替方向隐式迭代计算二维底部流场。守恒显式格式是从全守恒的Navier-Stokes方程组出发,在Allen—Cheng格式的基础上引进了“调整因子”。这对放大稳定性要求的步长,加速收敛均有益处。本文还用交替方向隐式迭代计算了同样的流场。结果表明只要边界条件逼近一致性较好,交替方向隐式迭代比显式格式收敛快,节约40％机时。     

涡扇发动机加减速特性显式与隐式计算方法

为了探索一种便于进行涡扇发动机过渡态控制规律设计的性能计算模型,提出了基于部件法的涡扇发动机加减速的显式格式和隐式格式计算方法,该方法通过在发动机计算模型中直接给定喘振裕度限制值、燃烧室油气比限制值和涡轮进口总温的限制值,计算出最优的加减速特性,进而获得发动机的最优加减速控制规律。计算模型针对不同的给定值,选择了不同的燃烧室容积效应模型。证明了对一般的涡扇发动机,隐式格式计算模型中,给定压气机喘振裕度算法的解是唯一的。以某涡扇发动机在地面的加减速过程为例,按最优加减速控制规律计算,显式格式算法和隐式格式算法的结果误差小于1.3%.对给定高压转子转速加速率的加速特性也进行了验算,计算结果与最优加速过程的结果误差小于1.7%.本文提出的加减速特性计算方法可为涡扇发动机的过渡态开环和闭环控制规律设计提供便捷的手段。     

求解N-S方程的一个简单隐式算法

1.引言 近十几年来,求解分离流动的N-S方程的数值解法,取得了很大进展。显式有限差分解已发展得较为成熟。但由于稳定条件的限制,获得稳定解的计算时间较长。为了克服这一困难,发展了隐式方法,如隐式因式分解法和交替方向的隐式方法。但是这些方法都要求解三对角矩阵方程,同显式方法比较,程序复杂,每步计算时间较长。为了避免三对角矩阵迭代的困难,MacCormack提出了新的隐式算法,其计算迭代过程为两对角矩阵追杆。本文在[1]的基础上,提出了一种更为简单的隐式算法,其迭代过程     

矢通量分裂隐式有限体积法解超音速喷流强干扰流场

发展了用于轴对称流场,任意非正交曲线网格下的矢通量分裂隐式有限体积法;统一了二维和轴对称流场计算的表达式。该方法对粘性通量项进行了分解和归类,考虑了包括交叉导数项在内的所有项对隐式增量的贡献。对带中心喷流的流动计算表明该方法收敛速度和计算精度均优于ＭａｃＣｏｒｍａｃｋ显式格式。横向喷流强干扰流场的计算初步揭示了姿控发动机喷流产生间接推力的机理,由此,设计中可降低对主发动机额定推力的要求     

动能BGK格式模拟高雷诺数定常流的收敛性研究

为克服动能BGK格式计算耗时多、收敛慢的缺点,在格式中引入当地时间步长、隐式LU-SGS方法和多重网格技术,基于RAE2822翼型粘性跨声速绕流的数值模拟,对其收敛性的改善进行了研究。结果表明:上述3种加速手段对求解定常问题时的收敛性有明显的改善,达到定常解所需的迭代步数随着CFL数的提高而减少;因而求解某类问题时,隐式格式比显式格式更具优势;多重网格既适用于显式格式,也适用于隐式格式,加速性更加突出。     

全状态两相湍流燃烧模型及其检验

在欧拉-拉氏体系中,在传统点源两相湍流燃烧模型(TM)中考虑液滴带火燃烧状态,提出全状态两相湍流燃烧模型(FSM),并用甲醇-空气两相湍流射流火焰实验数据对两种模型进行检验,对比结果表明:FSM模型能够给出平均温度的双峰分布,得到更合理的预报结果.      

辅助动力装置环形回流燃烧室数值研究

利用Fluent软件计算辅助动力装置(auxiliary power unit,简称APU)环形回流燃烧室三维两相喷雾燃烧流场,研究不同进口温度和油气比对两相喷雾燃烧流场的影响,采用标准 k-ε 模型模拟湍流黏性,离散相模型(DPM)追踪油珠运动轨迹,燃烧模型采用非预混平衡化学反应模型.计算结果表明:随着进口温度和油气比的增加,燃烧室出口温度相应增加,但温度分布规律基本保持不变;计算结果与试验测量结果比较吻合,说明采用的数学模型和计算方法可用于预估实际APU回流燃烧室三维两相喷雾燃烧流场.      

氢/氧气-气燃烧与气-液燃烧特性对比研究

对氢/氧气-气和气-液同轴剪切式单喷注器进行了燃烧流动的仿真研究.采用带化学反应的湍流Navier-Stokes方程和颗粒轨道模型描述发动机内部喷雾两相燃烧流动过程,气相化学反应速率都由Ar-rhenius公式计算.对典型气-气燃烧和气-液燃烧仿真结果进行了比较,结果表明气-气燃烧完成长度相对气-液燃烧更长;并进行了同轴喷注器关键参数对两种燃烧流场的影响的仿真和分析比较,得到喷注流量和动量比均为影响两类型喷注器燃烧流场的关键因素,且这两因素对燃烧完成长度的影响趋势是完全相同的,而喷注速度对两类型喷注器燃烧流场影响程度都较小.      

基于入射空气温度的无焰燃烧气动特性数值分析

提出了基于入射空气温度的无焰燃烧（AITTC）新型燃烧室,CFD数值研究了AITFC模型燃烧室内的流体流动结构。对头部引射混合区的引射混合、空气流量分配、回流区内气一雾滴颗粒两相流行为、以及尾部稀释区混合效应等气动特性的数值研究结果表明,所设计的新型模型燃烧室能满足组织AITFC无焰燃烧所需要的气动力学条件。研究结果符合相关的理论和规律。研究结果能为无焰燃烧技术应用于燃气轮机燃烧室设计提供参考价值。     

两种亚网格燃烧模型的旋流扩散燃烧大涡模拟

用二阶矩亚网格（SOM—SGS）燃烧模型和文献中的涡旋破碎亚网格（EBU—SGS）燃烧模型，对甲烷一空气旋流扩散燃烧进行了大涡模拟，将二者得到的LES统计平均温度分布和实验结果以及用二阶矩燃烧模型的统观模拟（RANS—SOM）结果比较，表明LES—SOM和RANS—SOM的模拟结果都和实验符合较好，而LES—EBU的模拟结果和实验不符合，在不同区域内高估或者低估了燃烧温度。其原因是由于EBU模型不能有效地考虑有限反应动力学的作用。LES—SOM模拟的瞬态结果显示了旋流扩散火焰的湍流大涡结构不同于射流火焰的特点。     

双燃式冲压发动机中超燃燃烧室冷态流场数值模拟

采用Ｍａｃ Ｃｏｒｍ ａｃｋ 矢通分裂格式, 求解雷诺平均的Ｎ－Ｓ方程组, 采用Ｂａｌｄｗ ｉｎ－Ｌｏ－ｍ ａｘ 代数湍流模型和混合长度模型, 模拟超燃燃烧室内两股平行超声速来流的冷态流场, 取得了与试验相一致的结果, 计算表明回流区和射流流动参数对超燃燃烧室流场内波系的变化有着较大的影响。     

湍流燃烧模型对氢燃料超燃室流场模拟的影响

采用化学平衡的假定概率密度函数(PDF)模型和火焰面模型计算了德国宇航研究中心的超燃室反应流,计算结果与有限速率反应模型的和实验的结果进行了对比.使用有限体积法离散Favre平均的N-S方程,湍流模型采用k-ε模型.研究表明:(1)有限速率反应模型在喷氢孔近场,化学平衡的假定PDF模型在喷氢孔远场不能准确捕捉流场的细致结构,而火焰面模型对全流场预测较好,后两种模型的计算时间较有限速率反应模型节省约38%;(2)超燃室内湍流和燃烧相互作用不可忽略,从预测精度和计算效率来看,火焰面模型有较好的工程应用前景.      

用不同化学反应模型对煤油超声速燃烧的数值分析

采用单步不可逆有限速率化学反应模型,同时结合离散液滴模型、概率密度函数紊流燃烧模型的混合分数平衡化学反应模型,对煤油在双模态燃烧室内的燃烧进行了数值分析。通过比较数值计算结果与实验结果,说明计算所采用的模型是有效的。在给定的计算条件下,与单步不可逆有限速率化学反应模型相比,结合了离散液滴模型的混合分数化学反应模型更能准确地预测煤油在双模态燃烧室内的喷雾燃烧,并且集燃料喷射、混合及火焰稳定为一体的凹槽火焰稳定器增强了燃料的混合和燃烧,使燃烧效率在燃烧室出口达0.880 6。      

三维加力燃烧室两相湍流燃烧的数值模拟

本文对涡扇发动机三维加力燃烧室内的两相湍流燃烧过程进行了数值模拟。两相流模型采用双流体模型,湍流模型采用标准 k-ε模型,湍流燃烧采用涡旋破碎 ( EBU)模型,数值方法采用 LEAGAP算法。计算结果定性合理。      

双环预混旋流低污染燃烧室数值研究

利用Fluent软件计算双环预混旋流(TAPS)低污染燃烧室三维两相喷雾燃烧流场,研究两种燃烧室结构和两种喷油方式对流场与燃烧性能的影响,采用标准k-ε模型模拟湍流黏性,离散相模型追踪油珠运动轨迹,燃烧模型采用非预混平衡化学反应模型.计算结果表明:在进口条件不变情况下,改进燃烧室结构和喷油方式,能提高出口温度,同时可大幅降低出口污染物排放;在相同试验条件下,TAPS低污染燃烧室燃烧性能优于目前某在研发动机模型燃烧室.