扭矩转化器（T／C）内部叶片流道流场的二维PIV实验研究

本文利用二维ＰＩＶ技术研究了汽车无级调连装置中扭矩转化器内部叶片流道流场的复杂流动,实验结果表明泵和透平在不同转速下,静子叶片流道中的流场体现出不同的流动特征,文中还对实验误差进行了分析,提出了减少误差的措施,通过叶栅流动理论分析说明了本文结果的可靠性。     

Accurate 3D viscous incompressible flow calculations with the FEM

A second-order-accurate (in both time and space) formulation is developed and implemented for solution of the three-dimensional incompressible Navier–Stokes equations involving high-Reynolds-number flows past complex configurations. For stabilization, only a fourth-order-accurate artificial dissipation term in the momentum equations is used. The finite element method (FEM) with an explicit time-marching scheme based on two-fractional-step integration is used for solution of the momentum equations. The element-by-element (EBE) technique is employed for solution of the auxiliary potential function equation in order to ease the memory requirements for matrix. The cubic cavity problem, the laminar flow past a sphere at various Reynolds numbers and the flow around the fuselage of a helicopter are successfully solved. Comparison of the results with the low-order solutions indicates that the flow details are depicted clearly even with coarse grids. © 1997 John Wiley & Sons, Ltd.     

二维翼型亚临界定常流动的有限元计算

本文使用了曲四边形Ｈｅｒｍｉｔｅ等参元，通过变分有限元法，计算了二维翼型亚临界定常位势流动中的压力分布，数值计算结果与试验结果符合较好。文中还使用了有限元法中的局部线性化理论，插值出元素中心点密度去处理迭代过程中元素系数阵的计算，其结果与通常的有限元法计算相一致，但计算时间却大大减少了。     

弯曲渠道中的振荡流和二次流

本文以弯道理想不可压缩流动的流线和等势线为曲线坐标系给出了弯道内不可压缩粘性振荡流动的一阶解以及边界层流型下的二阶解,以直角弯道为例对弯道內的粘性振荡流动以及二次流的影响进行了讨论。     

A PRESSURE-CORRECTION METHOD FOR THE SOLUTION OF INCOMPRESSIBLE VISCOUS FLOWS ON UNSTRUCTURED GRIDS

An unstructured grid, finite volume method is presented for the solution of two-dimensional viscous, incompressible flow. The method is based on the pressure-correction concept implemented on a semi-staggered grid. The computational procedure can handle cells of arbitrary shape, although solutions presented herein have been obtained only with meshes of triangular and quadrilateral cells. The discretization of the momentum equations is effected on dual cells surrounding the vertices of primary cells, while the pressure-correction equation applies to the primary-cell centroids and represents the conservation of mass across the primary cells. A special interpolation scheme s used to suppress pressure and velocity oscillations in cases where the semi-staggered arrangement does not ensure a sufficiently strong coupling between pressure and velocity to avoid such oscillations. Computational results presented for several viscous flows are shown to be in good agreement with analytical and experimental data reported in the open literature.     

粘性流体动力学有限元变分原理

本文把建立有限元变分原理的一种新方法“N>2直接方法”从固体力学推广到流体力学,并用该方法把粘性流体动力学的广义功率消耗原理和广义变分原理发展成为有限元变分原理。还在论证中发现,相邻有限元交界面上的应力协调条件会自然地满足而无需引进任何拉民乘子。本文还建立了混合杂交非协调元的变分原理和广义变分原理,它解除了全部协调性约束条件和其它的边界性约束条件,但是并不增加待定的拉氏乘子,因此使有限元计算得到简化。本文结果可以作为粘性流体动力学有限元计算的基础定理。     

UNIFIED COMPUTATION OF FLOW WITH COMPRESSIBLE AND INCOMPRESSIBLE FLUID BASED ON ROE''''S SCHEME

A unified numerical scheme for the solutions of the compressible and incompressible Navier-Stokes equations is investigated based on a time-derivative preconditioning algorithm. The primitive variables are pressure, velocities and temperature. The time integration scheme is used in conjunction with a finite volume discretization. The preconditioning is coupled with a high order implicit upwind scheme based on the definition of a Roe's type matrix. Computational capabilities are demonstrated through computations of high Mach number, middle Mach number, very low Mach number, and incompressible flow. It has also been demonstrated that the discontinuous surface in flow field can be captured for the implementation Roe's scheme.     

AN EULERIAN FINITE ELEMENT METHOD FOR TIME-DEPENDENT FREE SURFACE PROBLEMS IN HYDRODYNAMICS

A numerical method based on the finite element method is presented for simulating the two-dimensional transient motion of a viscous liquid with free surfaces. For ease of numerical treatment of the free surface expressed by a multiple-valued function, the marker particle method is employed. Numerous virtual particles are spread over all regions occupied by liquid. They move about on a fixed finite element mesh with the liquid velocity at their positions. These particles contribute nothing to the dynamics of the liquid and only serve as markers of liquid regions. The velocity field within liquid regions is calculated by solving the Navier– Stokes equations and the equation of continuity by the finite element method based on quadrilateral elements. A detailed discussion is given of the methodological problems arising in the implementation of the marker particle method on an unstructured finite element mesh and of the solutions to these problems. The proposed method is demonstrated on three sample problems: the broken dam problem, the impact of a falling liquid drop on a still liquid and the entry of a rigid block into water. Good agreement has been obtained in the comparison of the present numerical results with available experimental data.     

FEM ANALYSIS FOR INFLUENCES OF A FAULT ON COUPLED T-H-M-M PROCESS IN DUAL-POROSITY ROCK MASS

For the case in which a large geological structure like fault existing within the surrounding rock mass in the near field of a repository for high-level radioactive nuclear waste, one kind of coupled thermo-hydro-mechanical-migratory model of dual-porosity medium for saturatedunsaturated ubiquitous-joint rock mass was established. In the present model, the seepage field and the concentration field are double, but the stress field and the temperature field are single, and the influences of sets, spaces, angles, continuity ratios, stiffness of fractures on the constitutive relationship of the medium can be considered. At the same time, a two-dimensional program of finite element method was developed. Taking a hypothetical nuclear waste repository located at a rock mass being unsaturated dual-porosity medium as a calculation example, the FEM analysis for thermo-hydro-mechanical-migratory coupling were carried out under the condition of radioactive nuclide leaking for the cases with and without a fault, and the temperatures, pore pressures, flow velocities, nuclide concentrations and principal stresses in the rock mass were investigated. The results show that the fracture water in the fault flows is basically along the fault direction, and its flow velocity is almost three orders of magnitude higher than that of fracture water in rock mass; the nuclide concentration in the fault is also much higher than that without fault, and the nuclides move along the fault faster; moreover, the fault has obvious influences on the pore pressures and the principal stresses in the rock mass.     

带尾涡面的三维简化航天机亚音速有限元计算

本文采用Ｇａｌｅｒｋｉｎ有限元方法以全位势方程为控制方程计算了航天机三维简化模型的亚声速流动。为模拟实际流动,在后部加一尖劈形后体,并拖出一尾涡面。为局部超场 速区解的稳定,采用人工密度修正。对密度和环量进行双重迭代,得到了亚声速下三维航天机的压分分布和气动力系数。     

UNSTEADY DETACHED SEPARATION FROM A CIRCULAR CYLINDER PERFORMING ROTATIONAL OSCILLATIONS IN A UNIFORM VISCOUS INCOMPRESSIBLE FLOW

The origination of detached separation is studied on the basis of a numerical solution of the full Navier–Stokes equations. Fluxes of vorticity with different signs generated with twice the frequency of cylinder oscillation move from the cylinder to the outer surface of a detached liquid layer in the form of concentric rings. Near the critical layer between the attached layer and the main flow these rings are torn and crimped to the regions of separated vortices of the corresponding sign. The form of detached separated vortices is similar to that of vortices originating from a stationary circular cylinder in a uniform flow. Transition of the flow to a non-symmetric form with Karman vortex street generation at a Reynolds number (based on the radius) greater than 17 is revealed. This critical Reynolds number is smaller than that for a stationary circular cylinder in a viscous stream (where Re=20 has been determined to be a critical value) and corresponds to the Reynolds number extrapolated from the critical value for the stationary cylinder by increasing the cylinder radius by the attached layer thickness. The vorticity flux from the cylinder surface immediately into the separation region decreases as the frequency of cylinder oscillation increases. Violation of the flow potentiality in the detached separation region is the main cause of the vorticity generation on the outer surface of the attached liquid layer. © 1997 John Wiley & Sons, Ltd.