空化、超空化流动的数值模拟方法研究

基于结构化网格,运用可压缩流N—S方程及k-ε湍流模型对流场进行求解,在低压区域引入一种基于混合密度函数的空化模型对轴对称体的空化、超空化流动进行了数值模拟．通过将半球圆柱的计算结果与实验数据和前人的计算结果进行对比,验证了所发展的数值方法的可靠性．最后,采用非定常的数值方法,研究了钝头体射弹的空化、超空化流动特性,并模拟了其超空泡的发展过程．     

海堤越浪的数值模拟

基于RANS方程和两方程湍流模型,采用有限体积法,将人射波波场作为人工的分布源项加人动量方程,提出了适用于VOF方法的源造波一消波技术。通过对行波及驻波的计算,分别考察了数值波浪水槽前端及末端消波段的有效性。在本文建立的数值波浪水槽内对规则波在海堤上爬高和越浪过程进行了数值模拟,并将计算结果与现有实验结果进行了比较。验证计算结果表明,数值模拟结果较好地复演了海堤越浪过程。为了研究模型尺度对越浪量的影响,文中设计了一组满足重力相似但具有不同几何比尺的数值实验模型。系列数值实验结果表明,若按重力相似换算越浪量,计算结果与实验预报值间的偏差随模型比尺的增大和堤前波浪破碎强度的增强而增大,建议在进行越浪物理模型实验时需进一步考虑模型比尺对原型预报值的影响。     

基于Boltzmann模型方程的气体运动论统一算法研究

模型方程出发,研究确立含流态控制参数可描述不同流域气体流动特征的气体分子速度分布函数方程; 研究发展气体运动论离散速度坐标法, 借助非定常时间分裂数值计算方法和NND差分格式, 结合DSMC方法关于分子运动与碰撞去耦技术, 发展直接求解速度分布函数的气体运动论耦合迭代数值格式; 研制可用于物理空间各点宏观流动取矩的离散速度数值积分方法, 由此提出一套能有效模拟稀薄流到连续流不同流域气体流动问题统一算法. 通过对不同Knudsen数下一维激波内流动、二维圆柱、三维球体绕流数值计算表明, 计算结果与有关实验数据及其它途径研究结果(如DSMC模拟值、N-S数值解)吻合较好, 证实气体运动论统一算法求解各流域气体流动问题的可行性. 尝试将统一算法进行HPF并行化程序设计, 基于对球体绕流及类``神舟'返回舱外形绕流问题进行HPF初步并行试算, 显示出统一算法具有很好的并行可扩展性, 可望建立起新型的能有效模拟各流域飞行器绕流HPF并行算法研究方向. 通过将气体运动论统一算法推广应用于微槽道流动计算研究, 已初步发展起可靠模拟二维短微槽道流动数值算法; 通过对Couette流、Poiseuille流、压力驱动的二维短槽道流数值模拟, 证实该算法对微槽道气体流动问题具有较强的模拟能力, 可望发展起基于Boltzmann模型方程能可靠模拟MEMS微流动问题气体运动论数值计算方法研究途径.      

大气绕山体的运动及其输移特性的研究

对分层大气山体绕流的流动模式及扩散输移特性进行了数值模拟。采用隐式时间离散方法在贴体网格系统下求解雷诺平均的N-S方程,计算结果描述出大气流动的特征,证实了分层(以Froude数为特征参数,定义为F=U/NH,U为来流风速,N为Brunt-Vaisala频率,H为山体高度)变化对山体绕流流态的影响。数值结果表明:当Froude数大于4.0时,山体绕流的流动不再依赖于大气分层的变化。当Froude数介于4.0和1.0之间时,流场中出现了Lee波,并随着Froude数的进一步减小,流动分离发生及Lee波破碎现象。同时模型也预测了在各种流动模式下大气中夹杂着的污染物绕山体的传输特性,表明大气的分层现象对污染物的分布有着非常重要的影响。     

基于输运方程类空化模型的通气空泡流数值模拟

超空泡流动涉及多相、相变、湍流、可压缩性和非定常特性等复杂情况,亟待发展相应的精细数值模拟方法和计算软件.本文基于均质平衡流理论和输运方程类空化模型,提出了复杂气、汽、水多相的通气空泡流数学模型,采用有限体积方法和SIMPLE算法,建立了求解混合介质的雷诺平均N-S方程、湍流方程、各组分质量分数输运方程的数值模拟方法,自主开发了相应的三维计算程序.模拟了绕圆盘的轴对称通气超空泡流,计算结果与经典的近似公式、半经验公式及试验数据吻合良好,验证了数学模型和数值计算方法能够正确预报通气空泡流场的基本空泡形态和流场动力特性.     

高黏性流动有限元模拟的迭代稳定分步算法

分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中,已有的显式和半隐式分步算法由于 其显式本质,必须采用小时间步长计算,不但降低了计算效率,同时也常与为使分步算法稳 分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中,已有的显式和半隐式分步算法由于 其显式本质,必须采用小时间步长计算,不但降低了计算效率,同时也常与为使分步算法稳 定必须满足的最小时间步长要求冲突. 本文目的是构造一种含迭代格式的分步算法,它能在 保证精度的前提下大幅度地增大时间步长. 方腔流和平面Poisseuille流数值计算结果证实 了此特点,该方法被有效应用于充填流动过程的数值模拟.     

孤立波斜坡堤堤顶及后坡越浪流数值分析

孤立波往往被用来模拟海啸等浅水大波.运用基于RANS方程、VOF方法以及Goring造波方法建立了二维数值波浪水槽,进行了孤立波在简单斜坡堤上越浪过程的数值模拟.通过数值模拟,研究了不同相对堤顶宽及不同相对波高条件下对孤立波越堤流基本特征,包括堤顶前缘水层厚度及速度、堤顶末端水层厚度及速度、堤顶流厚度及速度沿程分布、后坡流厚度及速度沿程分布.依据数值实验结果,建立了堤顶流最大厚度分布及最大速度分布的经验公式.计算结果表明,堤顶流厚度沿程减小,而堤顶流速度则沿程增大;后坡流厚度先增大,随后再沿程缓慢减小;后坡流速度先减小,随后再沿程缓慢增大.     

数值模拟界面流方法进展

用数值模拟方法研究界面流,对于工程设计和理论分析复杂流动都是非常重要的.单纯用现有的高分辨格式无法清晰定出流体界面,故此对界面必须进行特殊的处理,例如界面跟踪(front tracking)方法．本文对界面流的各种数值模拟方法进行了综述．这些方法包括锋面跟踪法、移动网格法、基於粒子法、边界积方法、连续对流格式、体积跟踪方法和水平集的方法．我们着重介绍了现阶段处理流体界面问题使用较多的方法:体积跟踪方法和水平集方法．同时我们也对各种方法的优缺点作了评述.      

NUMERICAL INVESTIGATION ON FLOW STABILITY OF RAYLEIGH-B?NARD CONVECTION OF COLD WATER IN A RECTANGULAR CAVITY COOLED AND HEATED SYMMETRICALLY RELATIVE TO THE TEMPERATURE OF DENSITY MAXIMUM

为了解具有密度极值流体瑞利-贝纳德对流特有现象和规律,利用有限容积法对长方体腔内关于密度极值温度对称加热-冷却时冷水瑞利-贝纳德对流的分岔特性进行了三维数值模拟,得到了不同条件下的对流结构型态及其分岔序列,分析了密度极值特性、瑞利数、热边界条件以及宽深比对瑞利-贝纳德对流的影响. 结果表明:具有密度极值冷水瑞利-贝纳德对流系统较常规流体更加稳定,且流动型态及其分岔序列更加复杂;相同瑞利数下多种流型可以稳定共存,各流型在相互转变中存在滞后现象;随着宽深比的增加,流动更易失稳,对流传热能力增强;系统在导热侧壁时比绝热侧壁更加稳定,对流传热能力有所减弱;基于计算结果,采用线性回归方法,得到了热壁传热关联式.     

气体运动论数值算法在微槽道流中的应用研究

简要介绍基于Boltzmann模型方程的气体运动论数值算法基本思想及其对二维微槽道流动问题数值计算的推广，并阐述适用于微尺度流动问题的气体运动论边界条件数值处理方法。通过对压力驱动的二维微槽道流动问题进行数值模拟，将不同Knudsen数下的微槽道流计算结果分别与有关DSMC模拟值和经滑移流理论修正的N—S方程解进行比较分析，表明基于Boltzmann模型方程的气体运动论数值算法对微槽道气体流动问题具有很好的模拟能力。     

Preface

This special issue of PARTICUOLOGY is devoted to the first UK-China Particle Technology Forum taking place in Leeds, UK, on 1-3 April 2007. The forum was initiated by a number of UK and Chinese leading academics and organised by the University of Leeds in collaboration with Chinese Society of Particuology, Particle Technology Subject Group （PTSG） of the Institution of Chemical Engineers （IChemE）, Particle Characterisation Interest Group （PCIG） of the Royal Society of Chemistry （RSC） and International Fine Particle Research Institute （IFPRI）. The forum was supported financially by the Engineering and Physics Sciences Research Council （EPSRC） of United Kingdom,     

The 7th International Conference on Fluid Mechanics May 24-27,2015,Qingdao,China

正http://www.icfm7.org First Announcement and Call for PapersThe objective of International Conference on Fluid Mechanics(ICFM)is to provide a forum for researchers to exchange new ideas and recent advances in the fields of theoretical,experimental,computational Fluid Mechanics as well as interdisciplinary subjects.It was successfully convened by the Chinese Society of Theoretical and Applied Mechanics(CSTAM)in Beijing(1987,     

INFORMATION FOR CONTRIBUTORS

Contributions： The Journal, Acta Mechanica Solida Sinica, is pleased to receive papers from engineers and scientists working in various aspects of solid mechanics. All contributions are subject to critical review prior to acceptance and publication.     

Guide for authors

正Each of the sections below provides essential information for authors.We recommend that you take the time to read them before submitting a contribution to Acta Mechanica Sinica.We hope our guide to authors may help you navigate to the appropriate section.How to prepare a submission This document provides an outline of the editorial process involved in publishing a scientific paper in Acta Mechanica Sinica.     

基于鲁棒滤波的捷联导引头视线角速度估计方法

针对捷联导引头无法直接获取视线角速度等信息的问题,研究了鲁棒滤波在大气层外飞行器捷联导引头视线角速度估计中的应用。为了建立非线性滤波估计模型,考虑目标视线角速度的慢变特性,采用一阶马尔科夫模型建立了状态方程;推导了视线角速度的解耦模型,并建立了量测方程;考虑到实际应用中存在系统噪声统计特性失准的问题,基于Huber-Based鲁棒滤波方法,设计了视线角速度滤波器,并完成了基于Huber-Based滤波方法和扩展卡尔曼滤波方法的数学仿真。仿真结果表明Huber-Based滤波方法的视线角、视线角速度及视线角加速度估计精度分别达到0.1140'、0.1423'/s、0.0203'/s2,而扩展卡尔曼滤波方法的视线角、视线角速度及视线角加速度估计精度仅分别为0.6577'、0.6415'/s、0.0979'/s~2。仿真结果证明了该方法可以有效地估计出相对视线角速度等信息,并且在非高斯噪声的条件下,依然可获得较高的估计精度,具有一定的鲁棒性。     

Guide for authors

正Each of the sections below provides essential information for authors.We recommend that you take the time to read them before submitting a contribution to Acta Mechanica Sinica.We hope our guide to authors may help you navigate to the appropriate section.How to prepare a submission This document provides an outline of the editorial process involved in publishing a scientific paper in Acta Mechanica