直方管内充分发展湍流的大涡模拟 第二部分

本文第一部分用大涡模拟的方法对直方管内Ｒｅ＝４．９×１０￣４的充分发展湍流运动进行了数值模拟。所得结果证实大涡模拟是高雷诺数复杂湍流运动数值模拟的一个有效手段。在本文第二部分用所得的数据库对湍流运动的结构和形态进行了分析，并且用流体质点运动图形显示了湍流运动的二次流。     

直方管内充分发展湍流的大涡模拟 第一部分

本文用大涡模拟的方法对直方管内Ｒｅ＝４．９×１０￣４的充分发展湍流运动进行数值模拟，所得结果清楚地显示了二次流的存在。统计平均量的数值结果也与实验结果吻合甚好，这表明大涡模拟是高雷诺数复杂湍流运动数值模拟的一个有效手段。     

燃烧室-涡轮耦合流动传热超大涡模拟研究

为了准确预测发动机燃烧室和涡轮耦合复杂气热环境下涡轮部件的流动和换热特性,应用基于BSL k-ω湍流模型的高精度超大涡模拟方法(VLES),以及SST k-ω雷诺平均湍流模型对雷诺数为Re=3.8×105的高压涡轮导叶在均匀进口条件以及燃烧室-涡轮耦合情况进口条件下的涡轮流动和传热进行了数值研究.耦合火焰面生成流型(F...     

基于人工神经网络的可压缩湍流大涡模拟模型

在国家数值风洞（NNW）工程项目的指导下，空间人工神经网络（SANN）模型被用于强可压缩湍流大涡模拟（LES）研究，其中流场的湍流马赫数分别为0.6、0.8、1.0。基于湍流的多尺度空间结构特性和人工神经网络方法发展的高精度空间神经网络（SANN）模型适用于不可压缩湍流和弱可压缩湍流。对于强可压缩湍流，流场中会出现激波结构，给大涡模拟带来了挑战。本文的研究结果表明：SANN模型适用于强可压缩湍流的大涡模拟。在先验分析中，SANN模型预测的亚格子应力和亚格子热流的相关系数超过0.995，远远高于梯度模型和近似反卷积模型等传统模型；传统模型的相对误差大于30%，而SANN模型在这方面有很大的改进，相对误差低于11%。在后验分析中，与隐式大涡模拟（ILES）、动态Smagorinsky模型（DSM）、动态混合模型（DMM）相比，SANN模型能更精确地预测能谱、各类湍流统计特性以及瞬态流动结构。因此，基于湍流多尺度空间结构特性的人工神经网络模型加深了人们对强可压缩湍流亚格子建模的认识，同时可以服务于NNW工程的流体力学模型构造。     

湍流大涡数值模拟进展

本文简要陈述湍流大涡数值模拟的原理、优点,着重讨论湍流大涡数值模拟方法的关键问题及其可能解决的途径,包括脉动的过滤、亚格子模型、近壁模型和标量湍流的大涡数值模拟中的特殊问题.文章强调大涡数值模拟中亚格子应力的本质是可解尺度湍流和不可解尺度湍流动量间的输运,并以作者最近提出的新型亚格子模型说明发展亚格子模型的正确途径.文章最后提出湍流大涡数值模拟近期需要迫切解决的问题和其他具有挑战性的方向.     

洁净室流场大涡模拟

应用湍流大涡模拟方法对矢流洁净室和全顶棚送风下部两侧回风式洁净室进行了数值模拟。计算中使用了隐含流线抑风耗散作用的具有三阶部分展开的Taylor-Galerkin离散格式，导出了经Gauss滤波后的基本方程组的有限元列式，应用涡粘性亚格子模式的湍流大涡模拟方法进行数值模拟。数值计算结果与实验结果符合较好，表明湍流大涡模拟能够较准确地捕捉到洁净室流场中旋涡的位置和尺度，在洁净室的数值仿真方面，与标准的K－ε湍流模式相比，湍流大涡模拟结果具有较高的可靠性。     

基于格子Boltzmann方法的方通道湍流的大涡模拟

基于格子Boltzmann方程的大涡模拟方法,对以摩擦速度、方通道水利直径为特征尺度,雷诺数为300的直方通道内湍流流动进行数值计算.利用多松弛时间格子Boltzmann方法来模拟流场的流动,切应力改善亚格子应力模型来模化滤波后的非封闭项.将模化后的亚格子应力与格子Boltzmann方法中的松弛时间相关联,使得松弛时间当地化,从而能够准确地模拟湍流.对湍流的平均流向速度、平均二次流速度以表征湍流强度的均方根速度以及不同截面流向瞬时涡做了计算和评估.计算结果与直接数值模拟、实验数据相吻合,证明了格子Boltzmann方法在计算通道湍流中的精度.      

三维槽道湍流的大涡模拟

利用大涡模拟方法模拟了雷诺数为13000的充分发展槽道湍流流动。为了比较网格和亚格子模型对计算结果的影响,采用了三种不同的网格数目,以及Smagorinsky模型和动力模型等两种亚格子模型,并与直接数值模拟结果进行了比较。结果表明,大涡模拟可以比较准确地给出速度分布、脉动速度均方根分布等湍流统计参数。在此基础上,对近壁区小尺度湍流结构进行了研究。     

燃气轮机模型燃烧室的大涡模拟

燃烧室内的燃油雾化、蒸发以及和空气进行混合过程对燃烧过程有重要影响。提出1种基于大涡模拟的数学模型来模拟燃烧室内燃料喷射、蒸发和混合过程。被空间滤波掉的亚网格尺度涡对大尺度涡的影响由求单方程SGS湍流模型进行模拟。采用拉格朗日法和蒙特卡洛技术对流场中的喷雾粒子进行采样跟踪,采样喷雾粒子在流场中作为点源项与气相进行质量、动量和能量的双向耦合。提出1个基于SGS湍流动能的双向耦合模型来模拟SGS脉动速度对喷雾粒子运动的影响以及喷雾相对SGS湍流动能的影响。通过对1个同轴模型燃烧室中的喷雾蒸发及混合过程的大涡模拟,将预测结果和试验值进行了比较,预测值和试验值吻合良好,验证了模型的可靠性。     

高阶气体动理学格式在湍流数值模拟中的应用

本文回顾了高阶气体动理学格式在湍流数值模拟中的应用。与传统的Riemann求解器相比,气体动理学格式可以提供时空耦合的演化过程,这对发展高精度格式十分重要。因此,基于两步四阶时间离散和高精度WENO重构,发展了具有四阶时间精度的气体动理学格式。该格式有更高的数值精度和稳定性,并且具有更好的处理复杂流动问题的能力。目前,两步四阶格式已经成功地应用到低雷诺数湍流直接数值模拟和高雷诺数工程湍流RANS模拟中,包括低速槽道湍流、超声速均匀各向同性衰减湍流、二维亚声速翼型湍流和三维跨声速翼身湍流等。数值结果表明该格式对湍流直接数值模拟和湍流RANS模拟具有高数值精度和高数值稳定性。下一步将利用高阶气体动理学格式研究更具有挑战性的可压缩湍流问题,例如超声速湍流边界层和激波边界相互作用等。     

风沙流中近床面沙粒三维运动的LES-DEM分析

采用自行开发的大涡模拟( Large-Eddy Simulation,LES)和离散单元法(Discrete Element Method,DEM)耦合的数值计算程序,对风沙气固两相三维运动进行了数值模拟.对近床面沙粒的起跳速度、起跳角、入射速度和入射角的概率分布,以及沙粒展向速度的分布做出统计分析.结果显示:本文所建立两相耦合数值模型,能够复现风沙输移中散体沙粒的运动特征；入射速度和起跳速度服从对数正态分布,并随着来流平均风速增大,概率分布呈现集中趋势；入射角、起跳角分别集中在10°和15°左右,表明沙粒动能向垂直方向转移的趋势；沙粒的展向速度的随来流风速增大而更加显著.     

航空发动机进气及叶栅通道内砂粒动力学特性分析

为提高发动机的吞砂、防砂能力,需要了解砂粒在发动机中运动规律及对发动机性能的影响。采用欧拉-拉格朗日法对E3发动机整流罩、风扇和增压级进行了砂粒与气流的耦合作用分析计算,加入颗粒碰撞模型和侵蚀模型。仿真分析结果表明:砂粒运动轨迹和侵蚀区域符合发动机使用规律。砂粒进入发动机与风扇叶片压力面发生碰撞,砂粒发生碰撞后大部分进入了外涵流道。砂粒碰撞位置主要集中在风扇叶片压力面处,风扇对颗粒运动轨迹的影响最显著,碰撞后的颗粒有明显的径向运动趋势,叶顶区域的颗粒富集程度较高。整流罩、风扇叶片压力面及外机匣壁面都发生较为严重的砂粒侵蚀现象,而增压级叶片砂粒侵蚀现象不明显。      

风沙两相流动光学测量及图像处理技术研究进展

介绍了自风沙动力学诞生以来所使用的各类光学测量方法,阐述了各种方法的优缺点和适用范围,针对最新光学测量技术在近代风沙两相流研究中的重要作用,着重介绍了PIV、高速数字图像测量及其图像处理技术在风沙两相流中的应用,指出光学测量在时间和空间的解析能力上的持续增长将为沙粒微观力学解析和风沙流宏观统计规律获取提供充足的支持。     

航空发动机砂尘吸入物静电监测仿真实验

提出了航空发动机砂尘吸入物静电监测的仿真实验方法,实验以软件ANSYS电磁场分析模块建立有限元模型为基础,模拟不同粒径、荷质比、运动速度及质量浓度情况下砂尘吸入物的静电感应信号,并分别从时域与频域对感应电荷与电压信号进行分析,研究砂尘吸入物的粒径及其他宏观参数与静电监测信号之间的关系,建立用以表征砂尘颗粒粒径大小的特征指标。基于IDMS(进气监测系统)感应电压信号功率谱密度分布建立了特征向量,并以其曼哈顿距离与欧氏距离作为特征指标表征砂尘颗粒粒径大小。经仿真实验验证:特征指标与砂尘粒径呈正相关,且不受砂尘吸入物运动速度及荷质比变化的影响,但受砂尘吸入物质量浓度影响。进一步研究将利用质量浓度对特征指标进行修正,并开展验证实验。      

减压器关闭过程内部流场的动态仿真和特性分析

针对自主研制的大流量减压器,采用动网格技术和流固耦合仿真技术,将气体的流动和活动组件的运动耦合起来,实现了关闭过程中内部流场的动态仿真,描述了关闭过程减压器内压力和速度变化过程,掌握了流场中波系、回流区的演化情况。同时,通过对流场中不同位置参数进行监测,分析了监测点各参数随时间的变化情况,以及流量和开度的实时变化情况,研究了基于动态流场的关闭过程动态特性,揭示了影响动态特性的内在因素。     

Aeolian Processes and their Effects on Understanding the Chronology of Mars

Aeolian (wind) processes can transport particles over large distances on Mars, leading to the modification or removal of surface features, formation of new landforms, and mantling or burial of surfaces. Erosion of mantling deposits by wind deflation can exhume older surfaces. These processes and their effects on the surface must be taken into account in using impact crater statistics to derive chronologies on Mars. In addition, mapping the locations, relative ages, and orientations of aeolian features can provide insight into Martian weather, climate, and climate history.     

基于SPH方法的圆柱形贮箱液体晃动研究

在飞行器燃油供给方式的选择和飞行姿态控制调整中,必须充分考虑贮箱内燃油的晃动情况及其对贮箱冲击力的影响。然而液体晃动问题的高度非线性使得其研究极具挑战性。采用光滑质点流体动力学方法对圆柱形贮箱内液体的晃动进行建模仿真,将晃动液体的控制微分方程组离散到一系列流体粒子上,从而将不同物理量函数的连续积分转换成支持域内粒子的离散化求和形式,通过求解线性方程组得到数值解。在固壁边界上布置一组虚粒子,以对内部粒子产生排斥力来模拟固壁边界条件,同时采用狄利克雷边界条件方法模拟自由液面,实现对自由液面的追踪。将仿真的液面变化与试验结果作了对比,吻合较好。得到了贮箱壁面上不同监测点的压强变化以及贮箱与液体之间的相互作用力,并说明了结果的准确性。数值计算方法对液体晃动问题的研究提供一定参考价值。     

旋流预混燃烧室流动混合的大涡模拟

采用大涡模拟的方法研究旋流预混燃烧室中的流动及燃料空气混合过程,选取亚格子湍动能模式为基准的涡黏模型,通过模拟涡核的运动过程定性的表征了拟序结构的发展;引入混合不均匀度的概率密度函数,研究入口边界条件尤其是入口速度对于涡团破碎及混合程度的影响规律;通过研究速度分量与燃料质量分数脉动的相关性,分析了不同方向的速度脉动对于混合过程的影响因子.这些均为进一步研究旋流预混燃烧室的燃烧不稳定性及污染物排放提供了理论上的指导.      

韧性材料的冲蚀损伤数值仿真综述

砂尘造成韧性材料的冲蚀损伤现象和损伤演化规律是冲蚀机理研究的重要内容，基于相应冲蚀试验和数值仿真的分析 方法是揭示冲蚀机理的重要手段。回顾了砂尘冲蚀损伤研究的缘起和发展经历，在介绍砂尘冲蚀损伤机理的基础上，总结了当前 数值仿真研究的3种主流方法：有限元方法、无网格方法和计算流体动力学方法。从数值仿真角度针对砂尘对韧性材料造成的破 坏模式、砂尘冲击姿态、砂尘冲击速度和形状、砂尘破碎对材料去除机制的影响等4个方面综述了影响冲蚀损伤机理的主要因素 和当前研究进展。砂尘对韧性材料造成的破坏模式通常分为材料凹陷、滑动、犁削、切削等4种不同机制；砂尘的冲击姿态决定了 材料的损伤模式；砂尘的冲击速度和砂尘的形状决定了材料的损伤严重程度；砂尘破碎现象对材料的去除机制有一定影响。