翼型空化绕流数值研究

空化是发生在流体机械上的复杂过程,理论研究遇到很大困难.本文引入合适的空化数值模型,将空腔界面近似为自由面,用界面构造精度较高的流体体积方法求解空腔位置,通过直接求解原始变量的NavierStokes方程,数值模拟了无界域中空化在翼型上发生、发展和脱落的周期过程;并分析了空化产生对翼型表面的压力分布、翼型收到的阻力和升力的影响.结果表明,空化出现在翼型上表面;由于空化的产生,翼型表面压力分布不稳定,导致升力、阻力和流场压力出现波动,这是实际中产生噪声和损失的主要原因.     

非结构化网格下大范围波浪的折绕射计算

从含流缓坡方程出发,推导出光程函数方程。采用基于非结构化网格下的有限体积法对光程函数方程和波作用守恒方程进行数值离散和联合求解,从而构建了一个考虑绕射的计算近岸大范围波浪传播过程的数值模型。模型的空间步长不再受制于波长的限制,同时非结构化网格可以很好地拟合复杂岸线变形。模型分别通过了施奈尔定律、直立防波堤后的波浪绕射和圆形浅滩上的波浪变形验证,结果表明,该数值模式能有效模拟复杂边界条件和大范围水域下近岸波浪的传播过程中的折射和绕射等变形。     

基于非结构化同位网格的SIMPLE算法

通过基于非结构化网格的有限体积法对二维稳态Navier—Stokes方程进行了数值求解。其中对流项采用延迟修正的二阶格式进行离散;扩散项的离散采用二阶中心差分格式;对于压力-速度耦合利用SIMPLE算法进行处理;计算节点的布置采用同位网格技术,界面流速通过动量插值确定。本文对方腔驱动流、倾斜腔驱动流和圆柱外部绕流问题进行了计算,讨论了非结构化同位网格有限体积法在实现SIMPLE算法时,迭代次数与欠松弛系数的关系、不同网格情况的收敛性、同结构化网格的对比以及流场尾迹结构。通过和以往结果比较可知,本文的方法是准确和可信的。     

电磁力控制翼型绕流分离的增升减阻效率研究

电磁力可有效对流体流动进行控制,增升减阻,抑制流动分离,制约其推广应用的瓶颈为控制效率问题.为提高其控制效率,基于翼型绕流的电磁力控制,对电磁力增升减阻的控制效率问题进行数值研究. 根据能量守恒定律,推导电磁力控制能耗的比,基于升力和阻力计算节省能量. 定义电磁力的控制效率为能量节省与电磁力控制所需能耗的比值,研究不同工况下电磁力增升减阻的控制效率. 发现在控制开始阶段,电磁力能量主要消耗在增加边界层流体的动能上,电磁力控制效率非常低,但电磁力控制效率会随着电磁力工作时间的增长而增加;电磁力控制效率随着来流速度的增加呈指数下降;通过增加电磁力激活板的输入能量可增强电磁力的控制效果,但无法明显增加其控制效率.      

翼型绕流的电磁力控制

将表面包覆电磁激活板的翼型,按一定的攻角,置于流动的弱电介质溶液中,电磁激活板可产生作用于流体的切向电磁力(Lorentz力),从而改变流体边界层的结构. 在转动水槽中,对翼型绕流及电磁力控制下的绕流形态进行了实验研究. 结果表明,未加电磁力时,前缘涡的脱落点是不确定的,与流场具体条件有关,而后缘涡仅在尖角处脱落. 前缘涡与后缘涡相互影响,并周期性的脱体,在尾部形成涡街. 施加电磁力后,当力的方向与流动方向相同时,可以在一定程度上抑制分离,消除涡街,其效果与减小攻角类似. 加反向电磁力时,则相当于加大攻角,在翼型体的背风面形成涡街.      

均匀来流中旋转圆柱黏性绕流的数值研究

从不可压非定常Ｎ－Ｓ方程出发,首次数值求解了均匀来流中圆柱作周向旋转振荡的黏性绕流问题。探讨了旋转角速度振幅、振荡频率及Ｒｅ数等因素对流场结构及其非定常演化过程的影响,并根据计算结果,给出了在旋转振动频率。速度振幅平面内流场涡结构的分区图。     

均匀来流中旋转圆柱黏性绕流的数值研究

从不可压非定常Ｎ－Ｓ方程出发，首次数值求解了均匀来流中圆柱作周向旋转振荡的黏性绕流问题。探讨了旋转角速度振幅、振荡频率及Ｒｅ数等因素对流场结构及其非定常演化过程的影响，并根据计算结果，给出了在旋转振动频率。速度振幅平面内流场涡结构的分区图。     

非结构混合网格高超声速绕流与磁场干扰数值模拟

对均匀磁场干扰下的二维钝头体无粘高超声速流场进行了基于非结构混合网格的数值模拟.受磁流体力学方程组高度非线性的影响及考虑到数值模拟格式的精度,目前在此类流场的数值模拟中大多使用结构网格及有限差分方法,因而在三维复杂外形及复杂流场方面的研究受到限制.本文主要探索使用非结构网格(含混合网格)技术时的数值模拟方法.控制方程为耦合了Maxwell方程及无粘流体力学方程的磁流体力学方程组,数值离散格式采用Jameson有限体积格心格式,5步Runge-Kutta显式时间推进.计算模型为二维钝头体,初始磁场均匀分布.对不同磁感应强度影响下的高超声速流场进行了数值模拟,并与有限的资料进行了对比,得到了较符合的结果.     

一种基于三角形非结构化网格SIMPLE算法的程序设计

为了深入研究解决流动问题的基本算法,提出了一种三角形非结构化网格上SIMPLE算法的计算机程序设计。利用FORTRAN语言的特点,采用自定义变量存储控制单元以及界面信息。单元变量中定义了指向相邻单元和界面的指针,通过对网格文件所提供拓扑信息的读取对这些指针进行赋值,从而建立计算所需的基本数据结构。这样,就避免了采用数组直接存储网格拓扑信息,大大增加了程序的可读性、拓展性和二次开发潜力。通过顶盖驱动流、后台阶流以及圆柱绕流几个经典计算流体力学算例的验证,该程序计算收敛良好,精度较高,所得结果符合物理实际。因而,可在此基础上开发更具实际应用价值的三维程序。     

平面叶栅跨音流场气动特性的数值研究

通过对模型方程的分析，给出了一种新的隐格式构造思想。将它运用到关通量分裂格式中，可得到无近似因子分解、无矩阵运算的高效二阶精度隐式矢通量分裂差分格式，并用来直接求解时间平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组。数值计算标明：该方法具有精度高、稳定性好、计算量少、收敛快等优点，在平面叶栅跨音流场的计算中，较好地捕获了激波，与实验比较，结果令人满意。     

平流层双轴椭球体飞艇绕流场的数值分析

本文在CFD软件FLUENT6.3的计算平台上采用时均Navier—Stokes方程与Realizable k-ε湍流模型对平流层双轴椭球体飞艇进行了数值模拟。着重分析了与飞艇轴线垂直的截面的周向摩擦阻力系数和压力分布情况,飞艇背部的分离流动特点,物面剪切层分离流动规律,空间流态,以及它们随攻角的变化规律。结果发现,飞艇背部的分离流动是基本对称的。在小攻角下只存在闭式分离,而随着攻角的增大也出现了开式分离。分离区随着攻角的增大而增大,分离线的起始点位置随着攻角的增大向上游移动。15°攻角时没有发现二次分离,30°,40°,45°攻角下都发现了二次分离。在相同的攻角下伴随流动沿飞艇轴线向下游发展时,分离流动的区域在横向截面上增大了。     

Controlling Transonic Flow around Airfoils by Means of Local Pulsed Addition of Energy

The influence of local pulsed-periodic addition of energy into a supersonic region on the flow structure and wave drag of an airfoil in transonic flow regimes is considered by methods of mathematical modeling. The study reveals significant prospects of the considered method of controlling airfoil performance in transonic flow regimes, including wave-drag reduction.     

Numerical simulation of fluid flows using an unstructured finite volume method with adaptive tri‐tree grids

A tri‐tree grid generation procedure is developed together with a finite volume method on the unstructured grid for solving the Navier–Stokes equations. A hierarchic numbering system for the data structure is used. The grid is adapted by adding and removing cell elements dependent on the vorticity magnitude. A special treatment is developed to ensure good quality triangular elements around the cylinder boundary. The adopted finite volume method is based on the cell‐centred scheme. The pressure–velocity coupling is treated using the SIMPLE algorithm. A modified QUICK scheme for unstructured grids is derived. The developed method is used to simulate the flow past a single and multiple cylinders at low Reynolds number. The obtained results are in good agreement with the published data. Copyright © 2002 John Wiley & Sons, Ltd.     

Parallel DNS with Local Grid Refinement

A parallel finite volume method for unstructured grids is used for a direct numerical simulation of the flow around a sphere at Re = 5000 (based on the sphere diameter and undisturbed velocity). The observed flow structures are confirmed by visualization experiments. A quantitative analysis of the Reynolds averaged flow provides a data base for future model evaluations.     

陷落腔的流场数值计算研究

本文对有限分析解法的近似计算方法及其对计算精度的影响进行了数值计算和分析讨论；然后采用有限分析解法，对前、后壁壁缘高度不同的三种腔内流场运动发展过程及腔口剪切层运动情况作了数值计算，得到的腔口剪切层上腔口随边相互作用情况与对应腔的流场显示实验结果进行了比较，两者一致，计算结果表明，采用本文发展的近似方法计算可靠，可以较大地提高计算效率。     

NUMERICAL STUDY ON THE FLOW AROUND A CIRCULAR CYLINDER WITH SURFACE SUCTION OR BLOWING USING VORTICITY-VELOCITY METHOD

A vorticity-velocity method was used to study the incompressible viscous fluid flow around a circular cylinder with surface suction or blowing. The resulted high order implicit difference equations were effeciently solved by the modified incomplete LU decomposition conjugate gradient scheme (MILU-CG). The effects of surface suction or blowing ' s position and strength on the vortex structures in the cylinder wake, as well as on the drag and lift forces at Reynoldes number Re = 100 were investigated numerically. The results show that the suction on the shoulder of the cylinder or the blowing on the rear of the cylinder can effeciently suppress the asymmetry of the vortex wake in the transverse direction and greatly reduce the lift force; the suction on the shoulder of the cylinder, when its strength is properly chosen , can reduce the drag force significantly , too.     

Flows around moving bodies using a dynamic unstructured overset-grid method

In this article, a computational fluid dynamics algorithm is presented for simulations of complex unsteady flows around rigid moving bodies using an unstructured overset-grid method. For this purpose, a highly automated, three-dimensional, tetrahedral, unstructured overset-grid method is developed with one-cell-width overlapping zone in order to model the arbitrary geometries for steady and unsteady flow simulations. A method has been described to obtain the inter-grid boundaries of the one-cell-wide overlapping zone shared by a background grid and a minor grid. In the overset-grid methodology, vector intersection algorithm and bounding box techniques have been utilised. The mesh refinement and overset-scheme conservation studies proved the accuracy and efficiency of the method developed here. The applications of the developed algorithms were also performed through simulations that included complex internal flows around a flow-control butterfly valve as well as flows in an internal combustion engine with a moving piston. Lastly, validations with experimental data were conducted for both steady and unsteady flows around rigid bodies with relative motions.     

平板网架结构分析的超级条元法

本文针对复杂杆系结构提出一种超级有限条的分析方法。这种方法基于半连续半离散思想将各杆件结点自由度归结为整体结构自由度；进一步又将整体自由度简化为少数结线自由度上机计算。这将为复杂结构体系开拓一条简便、实用的工程分析途径。     

Numerical diffusion for flow‐aligned unstructured grids with application to estuarine modeling

The benefits of unstructured grids in hydrodynamic models are well understood but in many cases lead to greater numerical diffusion compared with methods available on structured grids. The flexible nature of unstructured grids, however, allows for the orientation of the grid to align locally with the dominant flow direction and thus decrease numerical diffusion. We investigate the relationship between grid alignment and diffusive errors in the context of scalar transport in a triangular, unstructured, 3‐D hydrodynamic code. Analytical results are presented for the 2‐D anisotropic numerical diffusion tensor and verified against idealized simulations. Results from two physically realistic estuarine simulations, differing only in grid alignment, show significant changes in gradients of salinity. Changes in scalar gradients are reflective of reduced numerical diffusion interacting with the complex 3‐D structure of the transporting flow. We also describe a method for utilizing flow fields from an unaligned grid to generate a flow‐aligned grid with minimal supervision. Copyright © 2013 John Wiley & Sons, Ltd.