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.     

Numerical boundary conditions for unsteady transonic flow calculations

In calculations of transonic flows it is necessary to limit the domain of computation to a size that is manageable by computers. At the boundary of the computational domain, boundary conditions are required to ensure a unique solution. Since wave solutions exist in the unsteady transonic flow field, incorrect boundary conditions may result in spurious reflections from the computational boundary. This may introduce errors into the solution. To prevent the spurious reflections, absorbing boundary conditions are often used on the computational boundary. In this paper we describe a method to derive absorbing boudary conditions for transonic calculations. We demonstrate both theoretically and numerically that the use of the absorbing boundary conditions will reduce the spurious reflections in the calculation.     

An efficient approach for calculation of pitching moment in nonlinear reduced frequency method at low Mach number transonic flows

In this research, an efficient methodology for calculation of pitching moment coefficient at low Mach number transonic flows by using the perturbed nonlinear reduced frequency approach is presented. The proposed approach uses the perturbation technique in the nonlinear frequency domain (NLFD) method to estimate the solution at high harmonics. In this approach, the density and velocity fields at high harmonics are perturbed about those at low harmonics. Perturbing the density and velocity fields, the semi‐linear form of the governing equations is obtained. The resulting solution vector and spatial operator are then approximated by discrete form of Fourier transformation and governing equations are solved by using the pseudo‐spectral approach. Numerical results show that the proposed approach predicts good pitching moment coefficient at low Mach number transonic flows with up to 50% savings in computational time. Copyright © 2011 John Wiley & Sons, Ltd.     

Time marching integral equation method for unsteady transonic flows

In this paper, we have proposed a time marching intregral equation method which does not have the limitation of the time linearized integral equation method in that the latter method can not satisfactorily simulate the shock-wave motions. Firstly, a model problem—one dimensional initial and boundary value wave problem is treated to clarify the basic idea of the new method. Then the method is implemented for 2-D and 3-D unsteady transonic flow problems. The introduction of the concept of a quasi-velocity-potential simplifies the time marching integral equations and the treatment of trailing vortex sheet condition. The numerical calculations show that the method is reasonable and reliable.     

全速解法在湍流跨音速流动中的应用

本文对不可压流动的常用算法SIMPLE算法进行了推广,使其能计算从亚音速到超音速一定马赫数范围内的流动,这里,可压缩流动和不可压缩流动的数值自满实现了统一,称为全速解法本文对全速解法在二维流动计算中的应用性进行了初步的研究,采用了非交错网格的有限体积方法对控制方程进行离散,并用动量插值法来求得连续方程中单元边界上的变量值,本文对全速解法在二维层流的计算效果进行了考核,而后又将此算法在湍流跨音速流动中应用,计算表明本方法是成功的,能够很好地反映各种马赫数下的流场特性。     

Projection conditions on the vorticity in viscous incompressible flows

The problem of establishing appropriate conditions for the vorticity transport equation is considered. It is shown that, in viscous incompressible flows, the boundary conditions on the velocity imply conditions of an integral type on the vorticity. These conditions determine a projection of the vorticity field on the linear manifold of the harmonic vector fields. Some computational consequences of the above result in two-dimensional calculations by means of the nonprimitive variables, stream function and vorticity, are examined. As an example of the application of the discrete analogue of the projection conditions, numerical solutions of the driven cavity problem are reported.     

Diffusion velocity for a three-dimensional vorticity field

A discussion is presented on the existence of a diffusion velocity for the vorticity vector that satisfies extensions of the Helmholtz vortex laws in a three-dimensional, incompressible, viscous fluid flow. A general form for the diffusion velocity is derived for a complex-lamellar vorticity field that satisfies the property that circulation is invariant about a region that is advected with the sum of the fluid velocity and the diffusion velocity. A consequence of this property is that vortex lines will be material lines with respect to this combined velocity field. The question of existence of diffusion velocity for a general three-dimensional vorticity field is shown to be equivalent to the question of existence of solutions of a certain Fredholm equation of the first kind. An example is given for which it is shown that a diffusion velocity satisfying this property does not, in general, exist. Properties of the simple expression for diffusion velocity for a complex-lamellar vorticity field are examined when applied to the more general case of an arbitrary three-dimensional flow. It is found that this form of diffusion velocity, while not satisfying the condition of circulation invariance, nevertheless has certain desirable properties for computation of viscous flows using Lagrangian vortex methods. The significance and structure of the noncomplex-lamellar part of the viscous diffusion term is examined for the special case of decaying homogeneous turbulence.     

Distinctive features of the transonic flow behind the trailing corner edge of an oversized cone-cylinder body

The results of an experimental study of the transonic flow behind the trailing corner edges of oversized cone-cylinder bodies of different dimensions are presented. Emphasis is placed on the induced transonic flow restructuring and the local aerodynamic forces accompanying the process, both steady and unsteady. The restructuring that occurs with increase in the freestream Mach number is studied.     

Steady-state solutions for some simple flows of generalized Burgers fluids

The steady-state solutions for three types of unsteady oscillating flows of generalized Burgers fluids are determined by means of the Fourier sine transforms. These solutions are also presented in equivalent forms in terms of elementary functions exp, sine, cosine, hyperbolic sine and hyperbolic cosine. The similar solutions for Burgers, Oldroyd-B, Maxwell, Second grade and Navier-Stokes fluids can be also obtained as limiting cases of our solutions.     

Transient flows of a second grade fluid

Exact analytical solutions for a class of unsteady unidirectional flows of an incompressible second-order fluid are constructed. The flows are generated impulsively from rest by motion of a plate or two plates or by sudden application of a pressure gradient. Expressions for velocity, flux and skin friction are obtained for both large and small times. It is found that large and small times solutions are dependent on the coefficient of viscoelasticity. The solutions corresponding to Newtonian fluids can be easily obtained from those for fluids of second order by letting the viscoelastic parameter to be zero.     

Transient solutions for three-dimensional lid-driven cavity flows by a least-squares finite element method

A time-accurate least-squares finite element method is used to simulate three-dimensional flows in a cubic cavity with a uniform moving top. The time- accurate solutions are obtained by the Crank-Nicolson method for time integration and Newton linearization for the convective terms with extensive linearization steps. A matrix-free algorithm of the Jacobi conjugate gradient method is used to solve the symmetric, positive definite linear system of equations. To show that the least-squares finite element method with the Jacobi conjugate gradient technique has promising potential to provide implicit, fully coupled and time-accurate solutions to large-scale three-dimensional fluid flows, we present results for three-dimensional lid-driven flows in a cubic cavity for Reynolds numbers up to 3200.     

Improved relaxation schemes for transonic potential calculations

A block relaxation scheme, grouped in a red-black ordering, is applied to transonic aerofoil calculations using body-fitted co-ordinates. The scheme is simple and easily vectorizable. Detailed comparisons with the approximate factorization method (AF2) are presented and it is shown that the new scheme is competitive in all cases considered. Transonic results, of engineering accuracy, on an 0-type grid of 149 × 30 points, are usually obtained within 200 iterations (? 40s on a Cyber 175).     

A vorticity stretching diagnostic for turbulent and transitional flows

Vorticity stretching in wall-bounded turbulent and transitional flows has been investigated by means of a new diagnostic measure, denoted by Γ, designed to pick up regions with large amounts of vorticity stretching. It is based on the maximum vorticity stretching component in every spatial point, thus yielding a three-dimensional scalar field. The measure was applied in four different flows with increasing complexity: (a) the near-wall cycle in an asymptotic suction boundary layer (ASBL), (b) K-type transition in a plane channel flow, (c) fully turbulent channel flow at Re _τ = 180 and (d) a complex turbulent three-dimensional separated flow. Instantaneous data show that the coherent structures associated with intense vorticity stretching in all four cases have the shape of flat ‘pancake’ structures in the vicinity of high-speed streaks, here denoted ‘h-type’ events. The other event found is of ‘l-type’, present on top of an unstable low-speed streak. These events (l-type) are further thought to be associated with the exponential growth of streamwise vorticity in the turbulent near-wall cycle. It was found that the largest occurrence of vorticity stretching in the fully turbulent wall-bounded flows is present at a wall-normal distance of y ⁺?=?6.5, i.e. in the transition between the viscous sublayer and buffer layer. The associated structures have a streamwise length of ~200–300 wall units. In K-type transition, the Γ-measure accurately locates the regions of interest, in particular the formation of high-speed streaks near the wall (h-type) and the appearance of the hairpin vortex (l-type). In the turbulent separated flow, the structures containing large amounts of vorticity stretching increase in size and magnitude in the shear layer upstream of the separation bubble but vanish in the backflow region itself. Overall, the measure proved to be useful in showing growing instabilities before they develop into structures, highlighting the mechanisms creating high shear region on a wall and showing turbulence creation associated with instantaneous separations.     

Fourier time spectral method for subsonic and transonic flows

The time accuracy of the exponentially accu-rate Fourier time spectral method (TSM) is examined and compared with a conventional 2nd-order backward differ-ence formula (BDF) method for periodic unsteady flows. In particular, detailed error analysis based on numerical com-putations is performed on the accuracy of resolving the local pressure coefficient and global integrated force coefficients for smooth subsonic and non-smooth transonic flows with moving shock waves on a pitching airfoil. For smooth sub-sonic flows, the Fourier TSM method offers a significant accuracy advantage over the BDF method for the predic-tion of both the local pressure coefficient and integrated force coefficients. For transonic flows where the motion of the discontinuous shock wave contributes significant higher-order harmonic contents to the local pressure fluctuations, a sufficient number of modes must be included before the Fourier TSM provides an advantage over the BDF method. The Fourier TSM, however, still offers better accuracy than the BDF method for integrated force coefficients even for transonic flows. A problem of non-symmetric solutions for symmetric periodic flows due to the use of odd numbers of intervals is uncovered and analyzed. A frequency-searching method is proposed for problems where the frequency is not known a priori. The method is tested on the vortex shedding problem of the flow over a circular cylinder.     

An inviscid model for the base pressure of transonic turbine cascade

An inviscid base pressure model for transonic turbine blade has been presented. It has been shown that for a given back pressure the base pressure at the trailing edge, and the profile loss of a turbine blade are fixed according to the model and the base pressure can be calculated with the help of an inviscid numerical scheme. A parameteric study on the model shows that a blade profile with positive curvature downstream of the throat is advantageous for generating less loss, whilst the worst situation is when the exit flow reaches the sonic condition.     

Anisotropic multistage finite element method for two dimensional viscous transonic flow in turbomachinery

A new method based on the anisotropic tensor force finite element and Taylor-Galerkin finite element is presented in the present paper. Its application to two-dimensional viscous transonic flow in turbomachinery improves the convergence rate and stability of calculation, and the results obtained agree well with the experimental measurements.     

Magnetohydrodynamic transient flows of a non-Newtonian fluid

Some exact solutions of the time-dependent partial differential equations are discussed for flows of an Oldroyd-B fluid. The fluid is electrically conducting and incompressible. The flows are generated by the impulsive motion of a boundary or by application of a constant pressure gradient. The method of Laplace transform is applied to obtain exact solutions. It is observed from the analysis that the governing differential equation for steady flow in an Oldroyd-B fluid is identical to that of the viscous fluid. Several results of interest are obtained as special cases of the presented analysis.     

Aerodynamic drag of bodies in transonic flow. Theory and applications to computational aerodynamics

Formulas for all the components of the aerodynamic drag (total, friction, inductive, wave, pressure, and heat-transfer) are uniformly derived as applied to flows governed by the Navier-Stokes and Reynolds equations. For flows of this type the definition of the aerodynamic drag components is refined and the physical basis of the chosen method of breaking up the total drag into components is discussed. Ways of calculating the aerodynamic drag components using the methods of computational aerodynamics are considered. On the basis of the refined formulas the drag components are calculated for flows around airfoils and a high-aspect-ratio wing in transonic flow.