THEORETICAL STUDY OF THREE-DIMENSIONAL NUMERICAL MANIFOLD METHOD

IntroductionNumerical manifold method is a new numerical method established on the basis of finitecover of manifold[1,2].By using continuous and non-continuous finite cover system,thenumerical manifold method includes the continuous and non-continuous as …     

Modes of vortex formation and transition to three-dimensionality in the wake of a freely vibrating cylinder

The present work is aimed to give some insight into the relation between vortex shedding modes and transition to three-dimensionality in the wake of a freely vibrating cylinder by establishing a numerical model and analyzing the relevant results of two- and three-dimensional simulations. The compressible flow past an elastically-mounted cylinder is solved by using the immersed boundary method (IB method). The cylinder is free to vibrate in the transverse direction with zero structure damping. The response of displacement amplitude is studied with the variation of reduced velocity. Whether P+S mode exists in three-dimensional flow and the occurrence of 2P mode is caused by flow transition from two-dimensional to three-dimensional are problems of concern. Both 2P and P+S wake modes are observed in two- and three-dimensional simulations. The numerical results indicate that the flow transition from two-dimensional to three-dimensional is coupled with the cylinder vibration in the synchronization/lock-in regime. The wake formation given by three-dimensional simulations suggests that the P+S mode might exist in reality when the flow is reverted to two-dimensional by vortex induced vibration (VIV) at Re=300–350. When Reynolds number increases to 425, the wake formation undergoes transition to three-dimensionality and 2P mode is observed. The effect of mass ratio on the flow transition to three-dimensionality is studied. The relationship between wake modes and aerodynamic forces is discussed.     

PDA measurements of droplet size and mass flux in the three-dimensional atomisation region of water jet in air cross-flow

The main objective of this paper is to develop a technique to measure the global droplet properties in the atomisation region of a water jet issuing in a high-speed air cross-flow. Knowledge of these global properties allows comparison of the break-up outcome of geometrically different water nozzles. This is achieved by extending a PDA system to enable measurements in three-dimensional droplet flows. First, the droplet size and the spatial droplet distribution are measured by the PDA method. The global droplet properties are then obtained by using the measured local mass flux as a weighting factor in integrating the local droplet size. To facilitate the measurement of mass flux in three-dimensional flows, the PDA method is extended so that the reference area for the mass flux is derived as a function of both the geometry of the measurement volume and the flow direction. In the present application of three-dimensional droplet flow (a water jet in air cross-flow), a simple method is developed to measure the three velocity components of droplets by means of a two-component PDA system. The paper outlines the measurement technique and the procedure of estimating the global droplet size and the global droplet size spectra from local droplet properties and local mass flux. Received: 26 July 1998/Accepted: 23 February 1999     

三向应力Mohr圆的真实构成及剪应力作用方向的确定

三向应力Mohr圆的构成在传统上是借助公式推证而得，并以平面图形来表示，缺乏三维的真实感和直观性。在应力应变分析中，对于平面应力问题，可以通过平面应力摩尔圆确定过一点不同斜面上的应力分量及其作用方向。对于三维问题，利用摩尔圆图解法可以确定某一斜面上的正应力和剪应力的数值，但不能表示剪应力的作用方向。剪应力的作用方向需要通过另外的图解方法来确定。本文分别从坐标系旋转和数值计算的角度解释了三向应力Mohr圆的构成过程，形象地说明了Mohr圆的物理本质。针对三向应力Mohr圆不能表示剪应力作用方向的问题，通过矢量运算，给出了剪应力作用方向的确定公式。     

Formation and evolution of the streamwise vortices in mixing layers

The evolution of the three-dimensional time-developing mixing layer is simulated numerically using the pseudo-spectral method. The initial perturbations used in this study consisted of the two-dimensional fundamental wave and the streamwise-invariant three-dimensional disturbance. A comparison of the formations of the streamwise vortices with different amplitude functions for three-dimensional disturbances is made. In one case the results are similar to that of Rogers & Moser^[1], whereas a different way in which the quadrupole forms and sudden expansion of the rib are observed in another case. The simulation also confirms that the stretching by the forming roller rather than Rayleigh centrifugal instability is responsible for the formation of the rib. Finally, numerical flow visualization results are presented. The project supported by the Zhejiang Province Natural Science Special Fund for Youth Scientistis' Cultivation.     

INVARIANT MANIFOLDS AND THEIR STABILITY IN A THREE－DIMENSIONAL MEASURE－PRESERVING MAPPING SYSTEMS

ＩＮＶＡＲＩＡＮＴＭＡＮＩＦＯＬＤＳＡＮＤＴＨＥＩＲＳＴＡＢＩＬＩＴＹＩＮＡＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＭＥＡＳＵＲＥ－ＰＲＥＳＥＲＶＩＮＧＭＡＰＰＩＮＧＳＹＳＴＥＭＳＬｉｕＪｉｅ（刘杰）（ＲｅｃｅｉｖｅｄＤｅｃ．１．１９９４;Ｃｏｍｍｕｎｉｃ...     

Analytic investigation of friction and heat transfer in the neighborhood of a three-dimensional stagnation point at low and moderate Reynolds numbers

A study is made of hypersonic three-dimensional flow of a viscous gas past blunt bodies at low and moderate Reynolds numbers with allowance for the effects of slip and a jump of the temperature across the surface. The equations of the three-dimensional viscous shock layer are solved by an integral method of successive approximation and a finite-difference method in the neighborhood of the stagnation point. In the first approximation of the method an analytic solution to the problem is found. Analysis of the obtained solution leads to the proposal of a simple formula by means of which the calculation of the heat flux to a three-dimensional stagnation point is reduced to the calculation of the heat flux to an axisymmetric stagnation point. A formula for the relative heat flux obtained by generalizing Cheng's well-known formula [1] is given. The accuracy and range of applicability of the obtained expressions are estimated by comparing the analytic and numerical solutions. Three-dimensional problems of the theory of a supersonic viscous shock layer at small Reynolds numbers were considered earlier in [2–5] in a similar formulation but without allowance for the effects of slip.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 143–150, March–April, 1988.     

Topological visualisation of focal structures in free shear flows

This paper describes a method for identifying and visualising the three-dimensional geometry of focal (vortex) structures in complex flows. The method is based primarily on the classification of the local topology as it is identified from the values of the velocity gradient tensor invariants. The identification of the local topology is reference frame invariant. Therefore, focal (vortex) structures can be unambiguously identified in these flows. A novel flow visualisation method is introduced whereby focal structures are rendered using a solid model view of the local topology. This new approach is applied to the identification of focal structures in three-dimensional plane mixing layer and plane wake flows.     

Numerical solution of 2D and 3D turbulent internal flow problems

The paper describes a method for solving numerically two-dimensional or axisymmetric, and three-dimensional turbulent internal flow problems. The method is based on an implicit upwinding relaxation scheme with an arbitrarily shaped conservative control volume. The compressible Reynolds-averaged Navier-Stokes equations are solved with a two-equation turbulence model. All these equations are expressed by using a non-orthogonal curvilinear co-ordinate system. The method is applied to study the compressible internal flow in modern power installations. It has been observed that predictions for two-dimensional and three-dimensional channels show very good agreement with experimental results.     

三维气泡与刚性壁面的相互作用研究

基于势流理论建立水下爆炸气泡运动三维模型,采用边界积分法求解拉普拉斯方程,得到气泡的变形及位置,并在计算过程中引入弹性网格技术,避免了因网格扭曲而导致的数值发散,进而模拟了刚性壁面附近三维气泡的动态特性。在数值模拟过程中,将本文计算值与实验数据进行对比分析,结果表明,计算值与实验数据吻合良好。在此基础上,分别模拟了弱浮力、强Bjerknes力,强浮力、弱Bjerknes力以及浮力与Bjerknes力相当时壁面附近气泡的运动特征,并将各种工况的计算结果与基于开尔文冲量理论(Kelvin Impulse)的Blake准则进行对比分析与讨论,得到了不同参数下气泡的运动特征。     

Development of finite-amplitude two-dimensional and three-dimensional disturbances in jet flows

The laminar-turbulent transition zone is investigated for a broad class of jet flows. The problem is considered in terms of the inviscid model. The solution of the initial-boundary value problem for three-dimensional unsteady Euler equations is found by the Bubnov-Galerkin method using the generalized Rayleigh approach [1–4]. The occurrence, subsequent nonlinear evolution and interaction of two-dimensional wave disturbances are studied, together with their secondary instability with respect to three-dimensional disturbances.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 8–19, September–October, 1985.     

Dissipative chaos, Shilnikov chaos and bursting oscillations in a three-dimensional autonomous system: theory and electronic implementation

A three-dimensional autonomous chaotic system is presented and physically implemented. Some basic dynamical properties and behaviors of this system are described in terms of symmetry, dissipative system, equilibria, eigenvalue structures, bifurcations, and phase portraits. By tuning the parameters, the system displays chaotic attractors of different shapes. For specific parameters, the system exhibits periodic and chaotic bursting oscillations which resemble the conventional heart sound signals. The existence of Shilnikov type of heteroclinic orbit in the three-dimensional system is proven using the undetermined coefficients method. As a result, Shilnikov criterion guarantees that the three-dimensional system has the horseshoe chaos. The corresponding electronic circuit is designed and implemented, exhibiting experimental chaotic attractors in accord with numerical simulations.     

Accurate measurement of three-dimensional deformations in deformable and rigid bodies using computer vision

Recently, digital-image-correlation techniques have been used to accurately determine two-dimensional in-plane displacements and strains. An extension of the two-dimensional method to the acquisition of accurate, three-dimensional surfacedisplacement data from a stereo pair of CCD cameras is presented in this paper. A pin-hole camera model is used to express the transformation relating three-dimensional world coordinates to two-dimensional computer-image coordinates by the use of camera extrinsic and intrinsic parameters. Accurate camera model parameters are obtained for each camera independently by (a) using several points which have three-dimensional world coordinates that are accurate within 0.001 mm and (b) using two-dimensional image-correlation methods that are accurate to within 0.05 pixels to obtain the computer-image coordinates of various object positions. A nonlinear, least-squares method is used to select the optimal camera parameters such that the deviations between the measured and estimated image positions are minimized. Using multiple orientations of the cameras, the accuracy of the methodology is tested by performing translation tests. Using theoretical error estimates, error analyses are presented. To verify the methodology for actual tests both the displacement field for a cantilever beam and also the surface, three-dimensional displacement and strain fields for a 304L stainless-steel compact-tension specimen were experimentally obtained using stereo vision. Results indicate that the three-dimensional measurement methodology, when combined with two-dimensional digital correlation for subpixel accuracy, is a viable tool for the accurate measurement of surface displacements and strains. Paper was presented at the 1989 SEM Spring Conference on Experimental Mechanics held in Cambridge, MA on May 29–June 1.     

Approximate solution of the problem of the three-dimensional boundary layer in an incompressible fluid

A method of successive approximations is proposed for the solution of the equations of the three-dimensional incompressible boundary layer on bodies of arbitrary shape. A coordinate system connected with the streamlines of the external nonviscous flow is used. It is assumed that the velocity across the external streamlines is small. When the intensity of secondary flow is low the equations describing the boundary layer in an incompressible fluid are reduced to a form analogous to the equations for the boundary layer on axially symmetrical bodies. An approximate analytical solution is obtained for the velocity and for the friction in the form of equations which can be used for any problems of a three-dimensional incompressible boundary layer. The method developed was applied to the problem of the three-dimensional boundary layer at a plate with a cylindrical obstacle in the presence of a slip angle.     

复合材料层合板的三维研究

本文提出了固支复合材料各向异性层合圆板受均布横向载荷作用下的满足三维弹性力学基本微分方程和边界条件的解析解答。文中采用一种发展的摄动方法进行求解，板中的每个应力和位移都展开为无量纲厚度参数ε的摄动级数，并采用二维板理论解答作为其相应三维摄动解答的一个基本解的形式，通过摄动方法逐级求解而获得完整的三维解答。文中以解析形式和数值形式给出了高精确度的三维应力和位移结果，结果表明，本文求解三维问题的解析方法是合理有效的。     

Spatial interaction of jets propagating in an accompanying supersonic flow

A numerical method of calculating a three-dimensional laminar supersonic underexpanded jet escaping into an accompanying supersonic flow is developed. The simplified Navier-Stokes equations for a steady-state three-dimensional flow are employed. Numerical calculations are carried out for several cases relating the outflow of jets from a four-nozzle assembly into an accompanying supersonic flow, and a number of the characteristics of three-dimensional flows of this kind are presented.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 88–93, November–December, 1972.     

Extension of e~N method to general three-dimensional boundary layers

In order to extend the e~N method to general three-dimensional boundary layers, the conservation law of the imaginary parts for the wave parameters with a fixed wave vector is deduced. The compatibility relationship(CR) and the general theory of ray tracing(RT), which have been extensively used in conservative systems, are applied to a general three-dimensional boundary layer belonging to non-conservative systems. Two kinds of e~N methods, i.e., the e~N -CR method and the e~N -RT method, are established.Both the two kinds of methods can be used to predict the evolutions of the spanwise wavenumber and the amplitude of the disturbances in general three-dimensional boundary layers. The reliability of the proposed methods is verified and validated by performing a direct numerical simulation(DNS) in a hypersonic general three-dimensional boundary layer over an aircraft model. The results are also compared with those obtained by other e~N methods, indicating that the proposed methods have great potential applications in improving the transition prediction accuracy in general three-dimensional boundary layers.     

Numerical solution of the transport equation for passive contaminants in three-dimensional complex terrains

This paper presents a numerical solution of the transport equation for heat and species in complex three-dimensional spaces. The solution domain of the equation is transformed into a cube, as also is the governing equation; the resultant equation is solved in the transformed space via a finite difference technique. The validity of the developed computer code is tested by predicting test cases for which either analytical or reliable experimental results exist. Results are also presented for the rate of convergence of the method and the computer storage requirements, from which the validity, the flexibility and the economy of the developed method are proved for flows in real three-dimensional complex terrains.     

Solution of pressure diffusion in radially composite reservoirs

This paper presents a new general method for solving the pressure-diffusion equation in cylindrically radial composite reservoirs, where the rock and fluid properties may change radially as a function ofr. Composite systems, such as formations with wellbore filtrate invasion and reservoirs with peripheral water encroachment, can be encountered as a result of drilling, secondary oil recovery, and water influx.The new solution method utilizes the reflection and transmission concept of electromagnetics to solve fluid flow problems in three-dimensional cylindrically radial reservoirs, where heterogeneity is in only one direction. The Green's function for a point source in a three-dimensional radially composite system is developed by using the reflection and transmission method. The method as well as the point source solution are sufficiently general that they may be applied to similar fluid flow and well testing problems involving single-phase flow.     

Numerical computation of a three-dimensional vortex sheet in a swirl flow

We consider an incompressible and inviscid fluid flow, called “swirl flow” that rotates around a certain axis in three-dimensional space. We investigate numerically the dynamics of a three-dimensional vortex sheet which is defined as a surface across which the velocity field of the swirl flow changes discontinuously. The vortex method and a fast summation method are implemented on a parallel computer. These numerical methods make it possible to compute the evolution of the vortex sheet for a long time and to describe the complex dynamics of the sheet.