二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.      

运用一维太阳风MHD CE/SE模型预报行星际激波到达时间

在一维太阳风磁流体(MHD)数值模拟中, 应用时空守恒元和解元(Conservation Element/Solution Element, CE/SE)方法, 建立了一个行星际激波扰动传播模型(1D-MHD (CE/SE)模型), 用来预报行星际激波到达时间. 收集了1997年2月至2002年8月间的137个激波事件, 对激波到达地球轨道附近的传播时间进行了预测, 并将预报结果与STOA, ISPM, HAFv.2以及SPM模型所得结果进行比较. 对于相同的样本事件, 1D-MHD (CE/SE)模型给出的渡越时间平均绝对值误差并不大于其他4个模型, 且该模型预报的相对误差小于10 %的事件占25.6 %, 小于30 %的事件占69.3 %, 小于50 %的事件占87.6 %, 其预报精度与其他模型相比基本相当. 这表明该模型在空间天气的激波到达时间预报方面有潜在的应用价值.      

r型网格自适应在间断Galerkin有限元激波捕捉中的应用

间断Galerkin（DG）有限元方法因计算精度高、适用于非结构网格等特点得到广泛研究和应用，其在数值模拟包含强间断流场时存在残差收敛性和计算鲁棒性差问题，均匀分布的网格加剧这一问题并影响激波分辨率。针对该问题，发展了r型网格自适应方法，实现间断Galerkin有限元数值模拟过程中网格自适应加密。基于网格点归一化的压力值作为r型网格自适应中网格点移动驱动力的重要权值，并将网格自适应后的网格点位移变化量与网格点之间的初始位移之比作为驱动力的另一重要权值，实现网格沿激波方向各向异性自适应加密，并且激波附近网格点的相邻网格点同步向激波方向移动。发展了适合间断Galerkin有限元方法的Venkatakrishnan限制器。并列NACA0012翼型超声速算例及三维并列圆柱相互干扰算例结果表明:基于r型网格自适应的间断Galerkin有限元方法能够清晰锐利捕捉激波，提高模拟精度，具有良好的收敛性和鲁棒性。      

3点逐步r型自适应网格算法

为了根据某个流场参数梯度较合理地分布网格同时又不增加网格点的数量,提出一种基以总变差为基础的r型3点逐步自适应网格算法,通过3点逐步调整网格点的位置以减小网格线的扭曲,生成质量较好的计算网格.以绕二维圆柱和三维双椭球的超声速流动为算例,用Euler方程模拟,以流场的压力总变差来进行自适应,通过分别用所提出的自适应方法生成的自适应网格与用原始网格计算所得的流场结果对比,计算结果表明,用自适应网格计算所得的激波比用初始网格计算所得的激波薄,采用自适应方法所计算出的流场数值解具有更高的分辨率.      

不同厚度翼型动态失速涡运动数值研究

在低马赫数下, 对三种不同相对厚度的NACA系列基本厚度翼型在俯仰振荡运动中的动态失速现象进行了数值研究.数值模拟时采用双时间法和LU-SGS隐式解法相结合的模式求解了非惯性坐标系下的非定常纳维-司托克斯(NS)方程组, 空间离散采用Roe格式,并结合刚性动网格生成技术.采用全湍流计算,通过引入BL(Baldwin-Lomax)模型计入湍流影响.计算出的气动力迟滞曲线与实验结果变化趋势符合较好,表明了数值方法的有效性,同时通过比较和分析不同厚度翼型在轻失速和深失速下的流场结构,发现绕翼型的失速旋涡产生和发展规律是明显不同的.      

基于简单WENO-间断Galerkin的Euler方程自适应计算

为了得到Euler方程的高精度、高分辨率数值解,介绍了间断Galerkin方法、三角形单元上简单WENO限制器的基本原理以及基于自适应网格加密的激波捕捉方法。将简单WENO限制器-间断Galerkin方法应用到曲边四边形单元上,通过单元边界上高斯积分点的坐标来搜索相邻单元从而得到相邻单元的单元编号,实现了基于“问题单元”的局部网格加密自适应计算。对若干典型问题进行编程计算,结果表明,简单WENO限制器可以应用到曲边四边形单元上,且可适用于局部网格加密时具有“悬挂节点”的非结构网格上的激波捕捉。      

翼型波浪水面巡航地面效应数值模拟

数值求解非定常不可压缩流动的雷诺平均N-S(Navier-Stokes)方程和标准k-ε湍流模型,用VOR(Volume Of Fluid)方法模拟波浪水面,模拟了NACA2410翼型在波浪水面上方飞行的流场.研究了余弦波浪水面的生成方法,选取合适的计算网格和时间步长,避免了在生成规则余弦波浪水面时由于数值耗散使波面形状衰减.比较了固壁波浪与水面波浪计算结果的差异,发现固壁波浪的结果更接近余弦曲线分布.研究了波浪等级对翼型气动性能的影响,分析计算结果发现:在规则的余弦水面波浪上方飞行,翼型气动力呈现周期性,给出了一个周期内气动力的变化过程,以及波长和波幅对气动力平均值和波动幅度的影响规律.      

圆弧翼型跨声速流动的动态模态分析

跨声速翼型的激波周期性自激振荡会给机翼结构带来附加的脉动载荷,从而加剧飞行器表面结构的疲劳损伤。使用动态模态分解（DMD）方法研究了跨声速下绕厚度18%的对称双圆弧翼型的压力脉动场,分析了DMD提取的各阶主模态的频率特征、压力脉动的空间分布以及压力脉动随激波振荡的时间演化过程,并使用DMD模态进行流场重构。结果表明,DMD方法能准确捕捉流场各特征频率的模态,第1阶模态是激波抖振的主频,在激波的自激振荡过程中占主导作用,前7阶模态重构的流场损失函数降低至4%以内,误差主要分布于激波间断处。      

一种适用于浸入有限元方法的网格自适应方法

针对动边界流固耦合的数值模拟问题，基于浸入有限元方法提出了一种耦合流场特征和几何特征的笛卡儿网格局部加密自适应方法，克服了单个自适应指示因子无法精确捕捉固体运动的特征的不足。在耦合自适应策略中，分别以流场涡量和固体位置作为流场和几何信息指示因子来驱动网格自适应。通过方腔顶盖驱动圆盘流动算例，以圆盘体积守恒和特征点的运动轨迹验证耦合自适应方法的优势。计算结果表明:仅基于流动特征的自适应不能很好地保证圆盘的体积守恒；仅基于几何特征的自适应无法有效追踪圆盘的轨迹；而耦合自适应策略能同时较好地保证两项指标的计算精度，在保证总体计算自由度不变的情况下，圆盘区域速度散度2-范数降低了一个数量级，圆盘的轨迹误差2-范数降低了2个数量级。      

结冰模拟与AEDC方法有效性的验证

介绍了AEDC(Arnold Engineering Development Center)结冰缩比方法,并对其有效性问题,用数值仿真方法进行了验证.利用结冰仿真计算软件对NACA0012翼型及其1/2和1/3缩比模型、圆柱及其1/2缩比模型进行了结冰数值模拟,内容涉及流场、水滴撞击特性和冰生长计算等.缩比部件的结冰条件由AEDC方法确定.结果表明,翼型、圆柱与其缩比模型的结冰过程具有较好的相似性.说明AEDC方法科学合理,其确定的结冰缩比试验参数有效,同时也说明结冰模拟技术可用于结冰缩比相似性研究.     

迎风格式在二维非结构网格中的应用

将原来在结构网格中提出的AUSM+ (AdvectionUpstreamSplittingMethod)格式 ,推广应用于二维非结构网格中 .利用格林 高斯公式对控制体内的变量进行线性重构 ,获得空间二阶精度 ,并采用Barth型限制器以抑制数值解的振荡 .为提高计算效率 ,采用了自适应多重网格法 .最后给出了NACA0 0 1 2翼型和平板激波反射的几个Euler方程的算例 ,并进行了分析比较 .     

迎风格式在二维非结构网格中的应用

将原来在结构网格中提出的AUSM⁺(Advection Upstream Splitting Method)格式,推广应用于二维非结构网格中.利用格林-高斯公式对控制体内的变量进行线性重构,获得空间二阶精度,并采用Barth型限制器以抑制数值解的振荡.为提高计算效率,采用了自适应多重网格法.最后给出了NACA0012翼型和平板激波反射的几个Euler方程的算例,并进行了分析比较.     

用GAO-YONG湍流模型计算激波/边界层干扰

采用GAO-YONG可压缩湍流方程组,模拟了平板激波/湍流边界层干扰现象.运用SIMPLE算法求解方程组,并分别采用三阶精度的QUICK格式和中心格式离散对流项和扩散项.计算结果较好预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征: 分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波.对流场的时均参数与实验值进行了比较,计算得到的壁面压力分布、摩阻系数分布和速度型与实验值比较吻合很好.结果表明GAO-YONG可压缩湍流方程组能够高精度模拟平板激波/湍流边界层干扰流动.      

低速流预处理Roe格式中的数值粘性

使用Roe格式计算低马赫数流动,需采用预处理技术.在预处理技术中,需要选择适合的预处理参数,但有时很难找到,因而得不到合理的计算结果.通过数值实验,发现使用预处理方法时,格式的数值粘性对低速流动的模拟有显著的影响.提出了在预处理Roe格式中加入可调参数,用以控制数值粘性对计算结果的影响.带有可调节数值粘性参数的预处理Roe格式的数值粘性更接近低速流动的真实耗散,反映了低速流动的真实物理现象.对方腔流动和圆柱绕流的低雷诺数定常流动所做的数值实验,得到了比较满意的结果,表明这一改进是有效的.对同一低雷诺数不同低马赫数的定常圆柱绕流所做的数值实验,均得到了满意的结果,表明了方法的适用性.     

Application of ADER Scheme in MHD Simulation

The Arbitrary accuracy Derivatives Riemann problem method(ADER) scheme is a new high order numerical scheme based on the concept of finite volume integration,and it is very easy to be extended up to any order of space and time accuracy by using a Taylor time expansion at the cell interface position.So far the approach has been applied successfully to flow mechanics problems.Our objective here is to carry out the extension of multidimensional ADER schemes to multidimensional MHD systems of conservation laws by calculating several MHD problems in one and two dimensions: (ⅰ) Brio-Wu shock tube problem,(ⅱ) Dai-Woodward shock tube problem,(ⅲ) Orszag-Tang MHD vortex problem.The numerical results prove that the ADER scheme possesses the ability to solve MHD problem,remains high order accuracy both in space and time,keeps precise in capturing the shock.Meanwhile,the compared tests show that the ADER scheme can restrain the oscillation and obtain the high order non-oscillatory result.     

守恒律方程的导数人工压缩法

将高分辨率差分格式用于守恒律方程的导数方程,可以克服传统高分辨率格式在极值点精度退化的缺点,类似于Harten的人工压缩法,新方法称为导数人工压缩法.本方法既能提高间断的分辨率,又能提高极值点处的分辨率,是一种高分辨率低耗散低扩散格式.用单个守恒律方程带间断和多极值的初值问题和一维激波管问题进行了验证,比较了Harten-TVD,人工压缩,导数人工压缩方法的在间断和极值点的分辨率问题.     

大气波动传播问题的一种无反射数值边界格式

根据无反射边界条件的基本原则,考虑到高阶的数值边界格式可以减小边界误差和虚假反射,利用曲线拟合中最小二乘法的思想,提出了一种高阶的光滑拟合外推边界格式(SFEBS).因为大气运动的控制方程简化后可以归结为一个对流方程,所以作为边界格式之间的比较和检验,在一维情况下用对流方程和数值模拟中具有代表性的波包和激波作为算例,将其与传统的基于Taylor展开思想构造的边界条件(TEBS)进行了比较.计算结果表明,在高阶情况下,用SFEBS计算波动传播问题的虚假反射约为用同阶TEBS计算的1/6,说明高阶的SFEBS可以作为一种很好的无返射数值边界格式.为今后模拟大气波动提供了一种很好的数值边界格式.