Vortex motion in the early stages of unsteady flow around a circular cylinder

In this paper, processes in the early stages of vortex motion and the development of flow structure behind an impulsively-started circular cylinder at high Reynolds number are investigated by combining the discrete vortex model with boundary layer theory, considering the separation of incoming flow boundary layer and rear shear layer in the recirculating flow region. The development of flow structure and vortex motion, particularly the formation and development of secondary vortex and a pair of secondary vortices and their effect on the flow field are calculated. The results clearly show that the flow structure and vortices motion went through a series of complicated processes before the symmetric main vortices change into asymmetric: development of main vortices induces secondary vortices; growth of the secondary vortices causes the main vortex sheets to break off and causes the symmetric main vortices to become “free” vortices, while a pair of secondary vortices is formed; then the vortex sheets, after breaking off, gradually extend downstream and the structure of a pair of secondary vortices becomes relaxed. These features of vortex motion look very much like the observed features in some available flow field visualizations. The action of the secondary vortices causes the main vortex sheets to break off and converts the main vortices into free vortices. This should be the immediate cause leading to the instability of the motion of the symmetric main vortices. The flow field structure such as the separation position of boundary layer and rear shear layer, the unsteady pressure distributions and the drag coefficient are calculated. Comparison with other results or experiments is also made.

Receptivity To Free-Stream Disturbance Waves For Hypersonic Flow Over A Blunt Cone

A high-order shock-fitting finite difference scheme is studied and used to do direction numerical simulation (DNS) of hypersonic unsteady flow over a blunt cone with fast acoustic waves in the free stream, and the receptivity problem in the blunt cone hypersonic boundary layers is studied. The results show that the acoustic waves are the strongest disturbance in the blunt cone hypersonic boundary layers. The wave modes of disturbance in the blunt cone boundary layers are first, second, and third modes which are generated and propagated downstream along the wall. The results also show that as the frequency decreases, the amplitudes of wave modes of disturbance increase, but there is a critical value. When frequency is over the critial value, the amplitudes decrease. Because of the discontinuity of curvature along the blunt cone body, the maximum amplitudes as a function of frequencies are not monotone.

Abel Inversion of Axially-Symmetric Shock Wave Flows

Finite-fringe interferograms produced for axisymmetric shock wave flows are analyzed by Fourier transform fringe analysis and an Abel inversion method to produce density field data for the validation of numerical models. For the Abel inversion process, we use basis functions to model phase data from axially-symmetric shock wave structure. Steady and unsteady flow problems are studied, and compared with numerical simulations. Good agreement between theoretical and experimental results is obtained when one set of basis functions is used during the inversion process, but the shock front is smeared when another is used. This is because each function in the second set of basis functions is infinitely differentiable, making them poorly-suited to the modelling of a step function as is required in the representation of a shock wave.

Numerical investigation of dynamic stall of an oscillating aerofoil

The flow structure around an NACA 0012 aerofoil oscillating in pitch around the quarter-chord is numerically investigated by solving the two-dimensional compressible N-S equations using a special matrix-splitting scheme. This scheme is of second-order accuracy in time and space and is computationally more efficient than the conventional flux-splitting scheme. A 'rigid' C-grid with 149 x 51 points is used for the computation of unsteady flow. The freestream Mach number varies from 0.2 to 0.6 and the Reynolds number from 5000 to 20,000. The reduced frequency equals 0.25-0.5. The basic flow structure of dynamic stall is described and the Reynolds number effect on dynamic stall is briefly discussed. The influence of the compressibility on dynamic stall is analysed in detail. Numerical results show that there is a significant influence of the compressibility on the formation and convection of the dynamic stall vortex. There is a certain influence of the Reynolds number on the flow structure. The average convection velocity of the dynamic stall vortex is approximately 0.348 times the freestream velocity.

三维直叶栅非定常流动的并行计算研究

流动分离直接关系到压气机运行的安全性与效率，对分离流动的研究是叶轮机械真实流动研究中的一个重大课题。本文针对三维压气机单转子叶片中截面所构成的三维直叶栅跨音速分离流开发了通用数值计算程序。该程序基于B-L湍流模型及高精度差分方法。多种工况的数值计算显示本程序结果与实验值吻合比较理想，验证了程序的正确性。10°攻角下分离区脉动压力的频谱与实验结果的数量级吻合，说明本程序能够较好地模拟大攻角分离流这种非定常复杂流动，为了提高计算规模及计算速度，作者对程序进行了并行化并针对微机机群系统进行了并行优化，实际计算表明本程序具有较高的并行效率。

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Receptivity to free-stream disturbance waves for blunt cone axial symmetry hypersonic boundary layer

Based on high-order compact upwind scheme, a high-order shock-fitting finite difference scheme is studied to simulate the generation of boundary layer disturbance waves due to free-stream waves. Both steady and unsteady flow solutions of the receptivity problem are obtained by resolving the full Navier-Stokes equations. The interactions of bow-shock and free-stream disturbance are researched. Direct numerical simulation (DNS) of receptivity to free-stream disturbances for blunt cone hypersonic boundary layers is performed.

西峰油田特低渗透油藏渗流机理及开发对策研究

Micro-pore-throat, micro-fracture and low permeability are the most obvious characters of Xifeng ultra-low permeability reservoir, and threshold pressure gradient and medium deformation during the period of oilfield developing results non-linear seepage feature of the formation liquid flowing in the porous medium underground. It is impossible to solve some problems in the ultra-low permeability reservoir development by current Darcy filtration theory and development techniques. In the view of the characters of ultra-low permeability and powerful-diagenesis and fracture up-growth, the paper quantitatively characterizes of through-going scope for reservoir parameters together with some materials such as similarity field outcrop, rock core, drilling, well logging and production dynamic, which provides geological base for further development adjustment. Based on the displacement experiment of different kinds of seepage fluid and oil-water two phases, this paper proves the relationship between threshold pressure gradient and formation permeability in experiment and theory, which is power function and its index is about -1. The variation rule and the mechanism of oil-water two phases threshold pressure gradient are studied. At the same time, based on the experiment of medium deformation, the variation rule of formation physical property parameters and the deformation mechanism are researched, and the influential factors on the medium deformation are analyzed systematically. With elastic unsteady filtration theory, nonlinear mathematical models of the steady and unsteady flow of single phase as well as horizontal well flow and oil-water two phases flow are deduced with the influence of nonlinear factors including threshold pressure gradient and media deformation. The influences of nonlinear factors upon well deliverability and reservoir pressure distribution as well as the saturation variation pattern of oil-water front are analyzed. By means of the researches such as reasonable well pattern, reasonable well array ration, artificial fracture length optimization advisable water flood timing and feasibility of advanced water flooding, it is necessary to find out effective techniques in order to improve development result of this kind of reservoir. This research result develops and improves on low-velocity nonlinear seepage theory, and offers ways to study similar kind of reservoir; it is meaningful to the development of the ultra-low permeability oil and gas reservoir.

Numerical simulation of axisymmetric unsteady incompressible flow by a vorticity-velocity method

A new numerical method for solving the axisymmetric unsteady incompressible Navier-Stokes equations using vorticity-velocity variables and a staggered grid is presented. The solution is advanced in time with an explicit two-stage Runge-Kutta method. At each stage a vector Poisson equation for velocity is solved. Some important aspects of staggering of the variable location, divergence-free correction to the velocity held by means of a suitably chosen scalar potential and numerical treatment of the vorticity boundary condition are examined. The axisymmetric spherical Couette flow between two concentric differentially rotating spheres is computed as an initial value problem. Comparison of the computational results using a staggered grid with those using a non-staggered grid shows that the staggered grid is superior to the non-staggered grid. The computed scenario of the transition from zero-vortex to two-vortex flow at moderate Reynolds number agrees with that simulated using a pseudospectral method, thus validating the temporal accuracy of our method.

Energy conversions and storage caused by an unsteady poloidal flow in active solar regions

In this paper we discuss coupling processes between a magnetic field and an unsteady plasma motion, and analyze the features of energy storage and conversions in active region. It is pointed out that the static force-free field is insufficient for a discussion of storage processes, and also the pure unsteady plasma rotation is not a perfect approach. In order to analyze the energy storage, we must consider the addition of poloidal plasma motion. The paper shows that because the unsteady poloidal flow is added and coupling occurs between the magnetic field and both the toroidal and the poloidal plasma flows, an unsteady process is maintained which changes the force-free factor with time. Hence, the energy in the lower levels can be transferred to the upper levels, and a considerable energy can be stored in the active region. Finally, another storage process is given which is due to the pure poloidal flow. The article shows that even if there is no twisted magnetic line of force, the energy in the lower levels may still be transferred to the upper levels and stored there.

Computation of Viscous Incompressible Flow using Pressure Correction Method on Unstructured Chimera Grid

In this paper, a pressure correction algorithm for computing incompressible flows is modified and implemented on unstructured Chimera grid. Schwarz method is used to couple the solutions of different sub-domains. A new interpolation to ensure consistency between primary variables and auxiliary variables is proposed. Other important issues such as global mass conservation and order of accuracy in the interpolations are also discussed. Two numerical simulations are successfully performed. They include one steady case, the lid-driven cavity and one unsteady case, the flow around a circular cylinder. The results demonstrate a very good performance of the proposed scheme on unstructured Chimera grids. It prevents the decoupling of pressure field in the overlapping region and requires only little modification to the existing unstructured Navier–Stokes (NS) solver. The numerical experiments show the reliability and potential of this method in applying to practical problems.

Rheological flow from a die and painting on a moving solid wall

Die swell is an important, phenomenon. in polymer processing, and is explained usually by rheological properties of the fluid. Because of the nonuniform of temperature distribution on the free surface of the liquid jet, the thermo capillary convection driven by surface tension gradient exists. The rheological fluid flowing out of a die and painting on a moving solid wall is studied by the numerical finite element method of a two-dimensional and unsteady model in the present paper, and both the rheological effect of a non-Newtonian fluid and the thermocapillary effect are considered. The results show that both,effects; will enlarge the cross-section of the fluid jet, and the rheological effect of non-Newtonian fluid dominates the process in general.

Drag on an ellipsoid particle in free-molecular plasma flow

Based on the analysis of molecular gas dynamics, the drag and moment acting on an ellipsoid particle of revolution X-2/a(2) + Y-2/a(2) + Z(2)/c(2) = 1, as an example of nonspherical particles, are studied under the condition of free-molecular plasma flow with thin plasma sheaths. A nonzero moment which causes nonspherical particle self-oscillation and self-rotation around its own axis in the plasma flow-similar to the pitching moment in aerodynamics-is discovered for the first time. When the ratio of axis length c/a is unity, the moment is zero and the drag formula are reduced to the well-known results of spherical particles. The effects of the particle-plasma relative velocity, the plasma temperature, and the particle materials on the drag and moment are also investigated.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.