基于虚网格的格心ADT搜索法

针对重叠网格寻点问题,提出了一种直接在格心网格下操作,基于虚网格的ADT(Alternating Digital Tree)搜索方法.通过对计算网格进行扩展,建立虚边界网格,解决了格心网格因覆盖区域不完整而无法直接建立ADT数据结构的问题.搜索结果即为包含解的格心单元集合,可直接对结果列表遍历以得到合理贡献单元,故完全摒弃了可靠性差的Stencil Walk方法.由虚网格的定义,使寻点在边界附近的处理更为灵活,可以准确给出边界附近贡献单元的有效信息,同时简化了虚网格体系的构建.扩展了重叠边界类型,构建的搜索空间能完整覆盖网格范围,解决了对称面重叠问题.算例研究表明:该方法可靠性好,边界处理能力强,有效提高了重叠网格方法对大型复杂网格的解算能力.     

加热器喷管热-流耦合传热分析

针对加热器喷管中复杂的气、固、液多相流动传热问题,建立了三维热流耦合换热计算模型,分别对燃气、冷却剂和喷管室壁建立不同的控制方程,将辐射热量作为源项加入到方程中,进行流动和传热的耦合计算.采用此方法对美国AEDC(Arnold Engineering Development Center)喷管的流动传热过程进行了计算,数值计算结果与试验结果吻合较好.在此基础上对某超燃冲压发动机试验台水冷式加热器喷管的换热问题进行了三维数值模拟,并定量分析了辐射换热对加热器喷管壁面温度分布的影响.结果表明:冷却水流量取2.0 kg/s时,加热器喷管气壁最高温度为660K,膜温度为430 K,加热器能可靠冷却,其热效率满足试验要求,对于含有H2O和CO2这样的高温燃气,辐射热量对喷管壁面温度分布有较大影响,必须引入到温度场的求解之中.     

结构/非结构网格中的数值优化计算

非结构网格因其在处理外形方面所具有的特殊优越性而倍受关注.在结构/非结构网格中做了粘性流动的数值优化计算,优化计算包括单目标和双目标气动优化,以及反设计计算,即从初始翼型不断修型而实现给定翼型的外形及压力分布的目标.算例及其结果分析证明了该方法用于翼型的修型和设计计算的有效性.      

应用于DSMC方法的二维Voronoi网格技术

研究了Voronoi网格技术并将其应用于直接模拟蒙特卡洛DSMC(Direct Simulation Monte Carlo)计算.基于Dirichlet镶嵌与Voronoi图理论,Voronoi网格利用特征点表征网格单元,具有建立粒子与网格单元之间映射关系的独特算法,适合于DSMC方法的统计特点.在剔除过于靠近边界的特征点以及必要情况下边界细化的基础上,通过区分由边界节点表征的非完整Voronoi网格单元以及由计算区域内镶嵌点表征的完整Voronoi网格单元,解决了Voronoi网格的二维边界匹配问题.Voronoi网格技术支持自适应DSMC计算.映射效率对比表明,Voronoi网格的DSMC计算效率高于三角形网格,低于多级直角网格.通过MEMS微喷管流动数值模拟,验证了Voronoi网格技术在DSMC方法中的有效性.      

矩形通道中亚尺度肋片的流动换热数值分析

对装有不同结构亚尺度肋片矩形通道的流动和换热进行了数值模拟,获得了通道流场、温度场分布以及平均努塞尔数的基本特征,对各种肋片通道的换热特性进行了对比分析.计算结果表明:亚尺度肋片距主肋片根部越近越有利于散热;在远离热面区域,亚尺度肋片的比表面积越大传热效果越好.当扩展表面积相同时,亚尺度肋片的长宽比越大换热效果越好;计算结果及分析揭示了控制肋片最优几何形状的统一原则——广义温度梯度均匀化原则,并依照此原则确定了可以强化换热的较优的肋片结构.计算分析表明广义温度梯度均匀化原则可以控制对流换热过程的传热强化.      

燃烧室三维贴体网格的生成方法研究

研究了燃烧室三维贴体网格人机交互生成的方法.使用该方法可以生成包括扩压器、内外环通道和火焰筒在内的燃烧室流场计算所需的整个网格.可以对网格的疏密和正交性进行调整.可以处理火焰筒上的主燃孔,掺混孔和气膜孔.并能够输出边界曲面上网格点的类型,用于流场计算.网格边界线的输入和边界线上的网格点分布采用人机交互方式生成.使用VISUAL BASIC语言编写程序界面.把整个网格生成及图形显示功能集成在一个软件中.     

2D约束Delaunay剖分生成表面模型的表面网格

提出一种生成表面模型的表面非结构化网格的方法.该方法仅假定表面模型的每个面片在取定适当的投影平面时可看成是一个单值函数.对表面模型的每个面片,首先在其相应的投影面上进行二维约束Delaunay剖分,然后对网格顶点进行插值.由于采用联动剖分的策略使得单独生成的各个表面网格在共同边界处匹配.与通常的生成表面网格的网格前沿法相比,本方法无需预先对域的边界进行离散,边界的离散体现在约束Delaunay剖分中恢复限定线段的边界细分过程中,减少了用户干预.通过合理地指定2D约束Delaunay三角化时网格单元尺寸分布函数,可以有效控制表面网格对表面模型的逼近精度和实现自适应的表面网格.实验结果证明了本算法的有效性.     

网格结构对Euler方程解的影响

 采用2种不同对接分区网格技术讨论了不同的网格结构对绕偏副翼的三维机翼流动的Euler方程解的影响.一种是部分不连续的网格(PDG);另一种是完全连续的网格(CCG),分别在PDG和CCG中采用Van Leer的矢通量分裂格式和Jameson的中心差分格式求解了Euler方程.算例表明,2种网格各有其优缺点,且都能得到有用的结果,但为了高效地取得更好的计算结果,必须根据流动和几何外形选取恰当的计算网格.     

复杂外形网格生成技术及应用

以某型双立尾战斗机为例,给出了实际复杂型号飞行器结构网格生成的一种分块方案及实现过程,用一种矢性三次多项式插值的方法生成相邻块公共交界面的网格,用Hilgenstock方法修正源项,求解椭圆型方程生成块内空间网格,使网格线与边界面夹角、与边界面间距控制为期望的值.在所生成的双立尾战斗机分块结构网格中求解Euler方程,获得了合理的全机气动力系数、翼面压力分布和绕流图谱结果.     

非线性问题的MPS无网格算法

给出了修正的MPS(Moving Particle Semi-implicit)无网格方法,指出恰当选取权函数可避免使用虚节点来处理边界,从而解决复杂边界处虚节点难设定的问题.以非线性Burgers方程为算例分析了其数值特性,结果表明:以样条函数为权函数时不采用虚节点的结果与使用虚节点的结果一致;对一维问题,光滑权函数性能明显优于非光滑权函数;而对二维问题,近边界处加权域不对称性的影响则不很显著.与精确解或由传统方法所得结果的对比分析表明该方案可行.该方法物理概念清晰,易于编程实现.      

龙卷风维持特性的探索

利用计算流体力学（CFD）方法建立了模拟龙卷风的装置模型。基于龙卷风在平面上的速度型,拟合成函数关系式,在三维CFD程序中设置初始速度场以及相关的边界条件,匹配不同的初始温度场进行了关于龙卷风维持和发展的一系列数值模拟。重现并利用漩涡管的Ranque-Hilsch效应以及汇流换热原理分析了本文2种不同的温度型对龙卷风内部流场变量随时间推移所发生的变化,同时揭示了龙卷风产生并维持下去的一种可能性,即冷热气流汇流换热生成温度场可能是产生和维持龙卷风的直接原因,并在此基础上提出使龙卷风削弱甚至消亡的方法,即破坏温度场。      

面向气热耦合的涡轮叶片计算域模型建模方法

针对冷却结构复杂的涡轮叶片在气热耦合数值模拟中计算域模型建模效率低、模型质量不稳定以及对数值模拟适应性差等问题,分析面向制造的涡轮叶片模型及建模方法,结合气热耦合数值模拟对涡轮叶片数值模拟的实际需求,提出一种面向气热耦合的涡轮叶片计算域建模方法。通过外型冷却特征自动定位生成算法,创建涡轮叶片冷气域部分;通过自适应管道求交算法、边界自动匹配算法,创建能够适应不同叶型的涡轮叶片燃气域部分;此外,在创建流体域的过程中提取气热耦合数值模拟所需的关键几何特征及非几何信息,并与冷气域、燃气域、叶片实体集成,完成气热耦合计算域模型的生成。基于以上研究开发涡轮叶片气热耦合计算域模型快速建模系统,验证所提方法的有效性。      

矩形微槽乙醇水溶液传热特性数值模拟

由于影响微槽传热特性的因素较多,实验研究结果相互矛盾,使得实验研究和数值仿真相结合成为进一步深入认识微槽传热特性的有效手段.在实验研究的基础上,采用有限体积法对矩形微槽中体积浓度30％的乙醇水溶液传热特性进行了三维数值模拟. 矩形微槽采用两种不同对称结构和热边界条件,进口采用流动完全发展边界条件.对不同对称结构和物性条件(常物性和变物性)的计算结果进行了比较,并与实验结果进行了对比.研究表明,对所研究微槽结构和工质,采用流动完全发展进口条件,以及随温度变化的热物性参数,数值计算结果与实验结果基本吻合.      

矩形微槽内FC-72的单相流动和换热实验研究

微槽换热器是解决航空电子设备元器件热控制问题的一条有效途径.以新型电子设备冷却液FC-72为工质,利用4种不同结构尺寸的矩形微槽,进行了微槽内单相流动和强迫对流换热性能实验研究;分析了工质流速、过冷度以及微槽结构等对换热性能的影响,给出了适用于层流和紊流的单相流动阻力特性实验关联式和单相强迫对流换热实验关联式,并对实验结果进行了不确定度分析.结果表明,在矩形微槽换热器中,FC-72的单相强迫对流换热能力能够满足中高强度的电子设备散热要求.      

两相控温型储液器进出流量的瞬态数值模拟

两相控温型储液器是泵驱两相流体回路（MPTL）系统中的一个重要部件,承担着工质存储、供给、气液分离及精密控温的作用。采用Navier-Stokes方程建立了MPTL系统瞬态模拟的仿真模型,可用于研究热源功率变化时储液器与主回路的动态传热和传质特性。通过仿真与试验对比发现,数值模型的流量误差在±10%以内,验证了模型的有效性和准确度。数值模拟结果表明:热源开关机时,储液器与主回路发生工质交换,气液两相的温度和压力受到影响,系统的流阻也受到影响;随着热源功率的增加,工质交换速率和交换总量随之增加,储液器内气液两相的温度和压力变化趋势随之增大。该模型可用于研究不同工作条件下的流量、温度和干度的变化特性,指导MPTL设计,并在系统搭建前预测系统特性。      

旋转状态下涡轮叶片前缘的流动与换热

用数值模拟的方法对旋转状态下涡轮叶片前缘冷却结构进行了数值研究,该结构由进气腔、叶片尾缘块和前缘块构成,对此结构不同的旋转速度情况进行了计算,根据计算结果分析了旋转对涡轮叶片前缘流动与换热的影响.计算结果表明,旋转状态下带气膜出流的冲击流动中,前尾缘冲击面的换热随着转速的增加而减小,且尾缘冲击面的换热比前缘冲击面的换热要好;同时前尾缘冲击面换热的差别随着转速的增加将越来越小.     

环肋翅片管相变传热特性及分形优化

基于管壳式相变换热器,通过数值模拟方法研究了35号石蜡在矩形环肋翅片管外的融化传热特性,建立肋片单元三维模型研究热流体温度、肋片高度参数对传热过程的影响,并探究了分形结构翅片的优化传热性能。结果表明：矩形环肋翅片管外相变融化过程分为传热速率差异明显的3段,适用于不同功率需求;提高热流体入口温度和相变材料温差近似等比增强总传热功率;肋片高度由10 cm分别提高至12.5 cm、15 cm时,总融化时间分别缩短42.89%、71.96%,增强传热,但功率重量比降低;分形结构优化翅片总融化时间缩短41.95%,功率提高,为相变翅片管的优化设计提供参考。     

毛细管束中层流粘性振荡流强化传热的若干问题

本文是作者从事振荡管流换热研究的总结。研究结果表明:流体在管内作轴向振荡运动时,可以极大地强化轴向的传热。本文给出了计及管壁效应的振荡管流换热的精确解析解,公布了变Pr数的实验结果。在此基础上,讨论了影响振荡管流轴向换热的6个无因次参数的物理含义,并指出其中的某些无因次参数在物理本质上就是傅里叶数的倒数的平方根,这将有利于理解振荡管流换热这种非定常换热过程。由本文的精确解析解,经过简化,可以得到国内外已公布的近似解。实验数据和解析解预言的基本规律符合很好。