基于格子Boltzmann方法表面形貌对微通道对流换热的影响

为推动微通道冷却技术在航空发动机中的应用,进一步提高航空发动机性能,本文通过扫描电子显微镜、粗糙度扫描仪等多方式扫描得到实际加工物体表面形貌,结合分形理论进行处理,得到真实、光滑、分形插值、康托尔集四类圆形截面管模型,管径100μm~400μm;采用格子Boltzmann方法进行数值模拟,首先通过与文献中实验结果的对比,验证了此方法的正确性,然后对空气在这四类圆形截面微通道中的对流换热特性进行了数值模拟,计算雷诺数Re范围80~640;结果发现：表面形貌是微通道中对流换热不可忽略的因素,相对粗糙度越大越有利于换热,真实扫描管的换热性能要比相同条件下光滑管高2.22%;分形理论可用于微通道表面形貌的构建,在相同迭代次数下,分形插值模型与真实扫描形貌换热性能相差0.03%,而康托尔集模型与真实扫描形貌换热性能相差1.88%,分形插值模型较康托尔集模型更能真实地反映实际形貌;微通道中,随着雷诺数的增加Nu将不再是常数,而是有增加的趋势。     

台阶后紊流场的大涡模拟

用大涡模拟法对台阶后复杂的紊流场进行了数值研究。用有限体积法求解滤波N－S方程，亚格子紊流模型采用结构模型。数值计算中，对流项，扩散项和时间分别采用三阶SHARP中式，二阶中心差分格式和三阶Runge-Kutta格式，给出了台阶后紊流场的时间统计特性，瞬时特性和紊流强度分布，并与试验结果进行了比较。     

圆柱尾流与游鱼间相互作用的数值计算研究

为了预测鱼类游泳问题的流动演变,开发了一种基于浸入边界格子玻尔兹曼方法的流动求解器.采用基于隐式边界力校正框架的灵活迭代算法降低计算成本和内存.在确定边界节点的速度修正后动量力由简单的直接力公式描述,无需复杂的积分计算.使用所提出的流动求解器,分析了鱼引起的动态失速涡与非定常流动中的传入涡之间的流体动力学相互作用.数值...     

不同亚格子尺度模型对气动噪声仿真精度的影响

大涡模拟（Large Eddy Simulation，LES）湍流模型在汽车外流场及气动噪声模拟计算分析中有着广泛的应用，高精度的非定常流场计算结果是有效进行气动噪声仿真分析的前提和基础。采用现代简易汽车模型（HSM模型），基于不同亚格子模型的LES对其进行非定常流场计算分析以及声场计算，并利用模型内部监测点的实验数据和仿真计算的结果进行对比验证。结果表明，使用WMLES模型可以得到更为准确的结果。     

不可压缩流动的高精度非标准格子玻尔兹曼方法

首先回顾了高精度非标准格子玻尔兹曼方法的发展历程,基于高精度通量重构格式,发展了一种通量重构格子玻尔兹曼方法（FRLBM）。采用两种求解方法：一种是将碰撞项隐式处理,直接求解离散速度玻尔兹曼方程（直接法）;另一种是先执行碰撞步,再求解纯对流方程（分步法）。通过收敛性研究,比较了这两种方法的精度和稳定性。研究结果表明,在小时间步长下两种方法误差近似,都能取得高阶精度;然而当时间步长增大,直接法误差几乎不变,分步法误差出现明显上升。由此表明,当取得近似误差时,直接法可以采用较大时间步长,计算效率更高,且直接法的稳定性略占优。接着通过模拟顶盖驱动方腔流验证了FRLBM捕捉流场细节的能力,并且比较了基于半隐格式的显式方法和一阶、二阶及三阶隐式-显式Runge-Kutta格式的时间离散,在不同雷诺数下的最大允许Courant-Friedrichs-Lewy数,数值结果表明二阶隐式-显式Runge-Kutta格式效果最优。最后数值模拟了圆柱绕流,验证了FRLBM计算复杂外形绕流的可靠性。     

亚格子模型对泵喷推进器宽带非定常力预测的影响研究

为了研究亚格子模型对泵喷推进器非定常流动与宽带非定常力预报结果的影响,采用分块结构化网格建立了模型尺度下艇后泵喷推进器的计算模型,并进行了大涡模拟数值计算。从艇尾非定常流场特征量和推进器转子脉动载荷两个方面对比了三种不同亚格子模型计算结果的差异,并分析了泵喷内部流动与转子非定常力间的内在联系。研究结果表明：三种亚格子模型均能得到含有叶频宽带峰的轴向推力谱,且整体趋势相近。泵喷转子上游的湍流强度和尺度的分布对亚格子模型较敏感,其中Smargrinsky-Lilly模型得到的湍流强度较强,尺度较大,该模型下的湍流谱在非平衡区量级较大,但由低频向高频的衰减较快,并且预测到的分离区范围大,导叶尾缘脱落涡分散,导致转子上游来流空间分布不均程度较强。对于转子叶片上的载荷脉动,Smargrinsky-Lilly模型预测的推力谱中线谱成分明显,并且叶频处宽带谱峰“陡峭”,而WALE模型和KET模型的结果宽带特性较强,对比非定常推力测试结果可知,WALE模型和KET模型更适于该问题宽带非定常力预测。     

导热-辐射耦合传热的多尺度分析和数值模型

复合材料高温传热特性的准确预测对飞行器热防护结构的设计有重要意义,相关模型也是国家数值风洞工程中多相多介质计算模型的重要组成部分。针对周期性结构复合材料的高温传热问题,利用多尺度渐进分析方法,对包含导热方程和辐射传输方程的导热-辐射耦合传热模型开展了研究。建立了表征单元尺度模型及宏观平均导热-辐射耦合传热模型,获得了材料宏观等效导热系数与表征单元模型间的关系,发现宏观等效辐射吸收和散射等系数可通过表征单元内的体积平均求取。结合有限容积方法与格子Boltzmann方法,建立了复合材料导热-辐射耦合传热多尺度数值模型。采用二维常物性材料传热过程的模拟验证了多尺度模型的有效性,结果表明所建立的多尺度模型能够对温度场给出准确高效的计算结果。该方法有助于为复合材料传热特性的预测提供数值手段。     

采用格子Boltzmann方法对低雷诺数平面射流的流场分析(英文)

采用格子Boltzmann方法对雷诺数(基于射流出口宽度和出口最大速度)为42和65的二维不可压缩平面射流进行了数值模拟。研究结果表明,在射流的自相似区,射流中心轴线速度按-1/3指数规律衰减并且射流按2/3指数率扩展,与理论分析及实验结果吻合。在射流的远场区域,中心轴线上瞬态速度的时间历程与频谱分析展示出速度的周期变化特性,并且周期性分量引起的相互作用不应被忽略;相比较而言,速度的周期性分量在近场区不明显。     

高超声速湍流高效模拟算法

采用了直接过滤的Navier-Stokes(N-S)方程组对高速可压缩湍流进行研究。针对高超声速湍流的非线性流动特性,对N-S方程直接过滤推导了大尺度湍流流场的控制方程,更精确地反映高速湍流的可压缩性,建立了可压缩湍流的大涡模拟TDM模型。使用传统的Smagorinsky模型的非线性推广,采用基于非Favre过滤的超声速可压缩湍流的大涡模拟模型,应用Caylay-Hamilton定理,建立可压缩湍流大涡模拟的非线性亚格子模型,并发展为动力学模式,模型中的两个常数通过当地流场动态的确定,消除了可调经验常数的影响。针对构造的高超声速湍流大涡模型开发相应的高效并行算法。     

一种剪切修正的涡黏性亚格子模型及其在槽道湍流中的应用

提出了一种利用壁面附近强剪切对涡黏系数进行修正的亚格子模型,称为剪切修正的涡黏模型.该模型能使涡黏系数在壁面附近自动满足壁面修正,在脉动为零的地方,涡黏系数自动为零,保证了壁面附近涡黏系数趋于零的特性.用该模型对槽道湍流壁面无吹吸、有吹吸,较低和较高雷诺数的情况进行了数值模拟,与直接数值模拟的结果进行了对比,结果吻合较好.