注塑成型三维流动的数值模拟

建立了非等温、粘性、不可压缩、非牛顿流体流动的控制方程。为了避免同时求解耦合的压力场、速度场,本文通过修改Galerkin方法的变分方程,导出了关于压力场的拟Poisson方程,用迭代法独立地求解连续性方程、动量方程,并进行速度一粘度迭代求出最终的压力场、速度场。由于直接使用Galerkin方法求解能量方程容易引起温度场的振荡,本文采用隐式格式及“上风”法离散能量方程,用超松驰迭代法求解温度场的代数方程组。比较了模拟结果与等温管道流动的解析解及法兰的实际注射结果,算例表明本文方法可以预测注射成型流动过程中的一些重要特征。与传统Galerkin方法相比,本文方法可以减少内存,提高数值方法的稳定性。     

快速求解叶栅流动的边界元方法

本文利用边界元方法，通过ｋｒｉｃｈｈｏｆｆ变换将描述叶栅流动的控制方程转换成线性方程。并将广义ｋ－Ｊ条件与边界积分方程联立求解，避免了非线性项和叶片出气角的迭代计算。完成了一种快速求解任意迥转面叶栅流场的计算程序，实用表明与其它数值方法及实验结果符合较好，具有快带、简明、实用的特点。     

剧变截面圆管内渗流的数值计算方法

对于剧变截面圆管的渗流问题写出不可压缩渗流的基本方程组,对直接求解原始变量(速度和压力)的数值计算方法作出改进。先由非主流方向的运动方程计算压力,后由主流方向的运动方程计算主流方向的速度分量,再由连续性方程计算非主流方向的速度分量。这样可以避免在一般的求解原始变量方法中由连续性方程计算压力时出现的困难和麻烦。根据本方法和剧变截面圆管的特点,采用半交错不等距非正交贴体混合网格系。本文详细写出差分方程和迭代计算公式,对剧变截面圆管内的渗流算例进行数值计算。本方法的优点是简单和实用,在工程上具有较大的应用价值。     

非正交曲线坐标的叶轮机械三元跨声速流方程和特征理论

本文使用张量运算,导出非正交曲线坐标下三元弱守恒型方程、时间相关流面迭代方程及特征理论.它们对于求解叶轮机械三元带激波跨声速流场,进行周向、径向畸变的数值分析及颤振不稳定气动力计算等都是有价值的. 1.非正交曲线坐标下的三元弱守恒型方程 叶轮机械非正交曲线坐标三元流方程     

FV/MC混合算法求解轴对称钝体后湍流流场

介绍一种有限容积／Monte Carlo结合求解湍流流场的相容的混合算法．有限容积法求解Reynolds平均的动量方程和能量方程，Monte Carlo方法求解模化的脉动速度—频率—标量联合的PDF方程．将该算法发展到无结构网格，探讨了在无结构网格中实现两种方法的耦合，包括颗粒定位，颗粒场和平均场之间数据交换等问题．并以二维轴对称钝体后湍流流场作为算例，比较了计算结果与实验结果．     

高阶紧致格式求解二维粘性不可压缩复杂流场

提出了一种求解二维不可压缩复杂流场的高精度算法．控制方程为原始变量、压力Ｐｏｉｓｓｏｎ方程提法．在任意曲线坐标下，采用四阶紧致格式求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，时间推进采用交替方向隐式（ＡＤＩ）格式，在非交错网格上用松弛法求解压力Ｐｏｉｓｓｏｎ方程．对于复杂的流场，采用了区域分解方法，并在每一时间步对各子域实施松弛迭代使之能精确地反映非定常流场．利用该算法计算了二维受驱空腔流动，弯管流动和垂直平板的突然起动问题．计算结果与实验结果和其他研究者的计算结果相比较吻合良好．对于平板起动流动，成功地模拟了流场中旋涡的生成以及Ｋａｒｍａｎ涡街的形成     

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

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

基于涡粘性模型的翼型湍流流场熵产率计算

利用有限体积法实现了基于非正交同位网格的SIMPLE算法。基于熵分析方法,采用涡粘性模型求解湍流熵产方程,系统研究了湍流模型对二维翼型绕流流场熵产率的影响。通过计算NACA0012翼型在来流雷诺数为2.88×106时,0°攻角~16.5°攻角范围内的翼型表面压力系数分布和升阻力特性,验证了算法及程序的正确性。结果表明,选择不同湍流模型时,翼型流场熵产的计算结果存在差异,湍流耗散是引起流场熵产的主要原因;翼型流场的熵产主要发生在翼型前缘区、壁面边界层和翼型尾流区域,流场熵产率与翼型阻力系数线性相关;当产生分离涡时,粘性耗散引起的熵产下降。     

用鲁棒Riemann求解器和运动重叠网格计算旋翼粘性绕流

发展了一种基于鲁棒Riemann求解器和运动重叠网格技术计算直升机悬停旋翼流场的方法。基于惯性坐标系,悬停旋翼流场是非定常流场,控制方程为可压缩Reynolds平均Navier-Stoke方程,其对流项采用Roe近似Reimann求解器离散,使用改进的五阶加权基本无振荡格式进行高阶重构,非定常时间推进采用含牛顿型LUSGS子迭代的全隐式双时间步方法。为实施旋转运动和便于捕捉尾迹,计算采用运动重叠网格技术。计算得到的桨叶表面压力分布及桨尖涡涡核位置都与实验结果吻合较好。数值结果表明:所发展方法对桨尖涡具有较高的分辨率,对激波具有较好的捕捉能力,该方法可进一步推广到前飞旋翼粘性绕流的计算。     

数值求解不可压粘性流体定常运动的格林函数方法

本文提出了一种数值求解不可压粘性流体定常运动的格林函数方法.在本文中利用Stokes方程的基本解作为格林函数将求解不可压粘性流体定常运动的边值问题化为求解速度场和边界应力的非线性积分方程组,在解出速度场和边界应力后可直接计算流场中各点的压力;用有限元近似将积分方程离散化而进行数值求解。对于小雷诺数流动,只归结为求解边界积分方程,使求解区域减少一个维度。对于非线性问题,可用迭代方法求解,在每次迭代中只须解出边界点上的速度或应力。通过几个简单的算例,表明本文所提出的方法具有精度高、处理边界条件简单、通用性强的优点,并具有求解各种复杂流动的潜力。     

Conservative calculations of non-isentropic transonic flows

The classical potential formulation of inviscid transonic flows is modified to account for non-isentropic effects. The density is determined in terms of the speed as well as the pressure, which in turn is calculated from a second-order mixed-type equation derived via differentiating the momentum equations. The present model differs in general from the exact inviscid Euler equations since the flow is assumed irrotational. On the other hand, since the shocks are not isentropic, they are weaker and are placed further upstream compared to the classical potential solution. Furthermore, the streamline leaving the aerofoil does not necessarily bisect the trailing edge. Results for the present conservative calculations are presented for non-lifting and lifting aerofoils at subsonic and transonic speeds and compared to potential and Euler solutions.     

Lifting aerofoil calculation using the boundary element method

The bbundary integral formulation and boundary element method are extended to include lifting flow problems. This involves inclusion of a branch cut in the flow field and imposition of a Kutta condition to determine the circulation, Γ Additional boundary integral contributions arise from the cut surface. Techniques for calculating Γ are developed and we treat, in particular, a superposition procedure which permits very efficient computation. Numerical results are presented for an NACA0012 aerofoil at several angles of attack.     

Aeroacoustic Modelling of Low-Speed Flows

A new numerical algorithm for acoustic noise generation is developed. The approach involves two steps comprising an incompressible flow part and an inviscid acoustic part. The acoustic part can be started at any time of the incompressible computation. The formulation can be applied both for isentropic flows and non-isentropic flows. The model is validated for the cases of an isentropic pulsating sphere and non-isentropic flows past a circular cylinder. In the latter case the computations show that the generated acoustic field in addition to the dominant Strouhal frequency, f ₀, contains a slightly higher frequency, f ₂, and a modulating lower frequency, f ₁=f ₂−f ₀. Numerical experiments with different interpolation schemes, boundary conditions, etc., show that the appearance of these modes is not an artifact from the numerical discretization. Received 20 July 1998 and accepted 26 May 1999     

A new approach to the solution of the Navier-Stokes equations

In the present paper a numerical algorithm is given for solving a standard problem in fluid dynamics, that of inviscid, irrotational, incompressible flow over an arbitrary symmetric profile. The purpose of the paper is to propose an alternative approach to solve certain fluid dynamic flows. This paper may be thought of as the first of a possible series of papers solving new and fundamental problems. In a sense, this new approach asks the question: what is the simplest and most efficient method of solving the problem considered by finite difference methods. It is believed that the following algorithm answers this question. Standard second-order finite difference techniques, such as SLOR and ADI, are used to solve numerically a mixed boundary value problem comprised of a pair of elliptic partial differential equations with constant coefficients.     

水下爆炸非均熵二维定常流的三族特征线解法

针对二维定常可压缩超声速非等熵柱状流,提出一种特征线差分解法,通过在沿马赫线的相容方程中添加沿流线的熵变项以描述非等熵效应,得到等熵流和非等熵流均适用的三族特征线方程组。根据水下爆炸近场特点,建立无限长柱状装药的定常模型,将三族特征线方程组用有限差分法离散求解,通过构造合适的网格保证计算格式可以数值上收敛,由此编制程序并计算几种柱状炸药的水下爆炸近场冲击波。对比有限元模拟结果和实验结果发现,特征线差分法可以比较准确地捕捉冲击波形状并计算冲击波后流场,从而验证了所提出的三族特征线差分法的准确性。     

Comparison between the traction and pressure-imposed boundary conditions in finite element flow analysis

A formulation is developed to impose pressure-prescribed boundary conditions in the penalty finite element method. Some numerical experiments for the Poiseuille flow problem are performed to compare it with the conventional traction-prescribed boundary condition. Also the incorrectness of the traction-free outlet boundary condition for contained-flows is studied with explanatory numerical examples. Discussion is focused on the inlet and outlet boundary conditions to simulate fully developed flows. Finally, the three-dimensional flow in a bifurcated pipe is analysed with the proposed formulation.     

A design method for two-dimensional cascades of turbomachinery blades

A design method for two-dimensional cascades of turbomachinery blades is presented. A finite element potential flow program is extended to allow fluid to transpire through the blade surface, the displaced surface streamline defining a new blade geometry. The potential changes are related linearly to the transpired flow rates. New surface velocities may then be specified as a function of surface distance, in accordance with boundary layer considerations. Closure and smoothness of the new blade are successfully achieved, while large changes in the blade geometry are possible.     

含铝炸药爆轰驱动的非线性特征线模型

针对含铝炸药爆轰的非理想特性,提出了含铝炸药爆轰产物膨胀的局部等熵假设,建立含铝炸药爆轰驱动的非线性特征线模型,为研究含铝炸药爆轰产物的非等熵流动和膨胀做功提供了一种新的理论分析方法。设计了5、50 μm含铝炸药和含LiF炸药驱动0.5、1 mm厚金属板实验,通过激光位移干涉仪测试金属板运动的速度历程,再通过实验结果计算得到铝粉在爆轰产物中的反应度变化规律,结合含铝炸药爆轰产物的非线性特征线模型,理论计算了含铝炸药驱动金属板的速度历程。对比理论与实验结果,理论方法能够很好地描述铝粉二次反应对炸药做功能力的贡献,同时验证了含铝炸药爆轰驱动的非线性特征线模型的正确性。     

Numerical solution of viscous flows using integral equation methods

A formulation of the boundary element method for the solution of non-zero Reynolds number incompressible flows in which the non-linear terms are lumped together to form a forcing function is presented. Solutions can be obtained at low to moderate Reynolds numbers. The method was tested using the flow of a fluid in a two-dimensional converging channel (Hamel flow) for which an exact solution is available. An axisymmetric formulation is demonstrated by examining the drag experienced by a sphere held stationary in uniform flow. Performance of the method was satisfactory. New results for an axisymmetric free jet at zero Reynolds number obtained using the boundary element method are also included. The method is ideal for this type of free-surface problem.     

Measurement of confined supersonic, 2-d jet lengths using the hydraulic analogy

 In this study, inviscid, confined supersonic and subsonic stream interactions are studied using a hydraulic analog experiment. Specifically, the length of the first minimum location of the slipline is measured. Reasonable agreement is shown between a theoretical model and analog measurements. A method for extending the isentropic limited hydraulic analogy to model these inherently non-isentropic flows is also developed. Received: 28 October 1996/Accepted: 10 June 1997