Numerical studies of slow viscous rotating flow past a sphere—1

The Navier-Stokes equations for a steady, viscous rotating fluid, rotating about the z-axis with angular velocity ω are linearized using the Stokes approximation. The linearized Navier-Stokes equations governing the axisymmetric flow can be written as three coupled partial differential equations for the stream function, vorticity and rotational velocity components. One parameter, R_eω = 2ωa²/v, enters the resulting equations. For R_eω « 1, the coupled equations are solved by the Peaceman-Rachford A.D.I. (Alternating Direction Implicit) method and the resulting algebraic equations are solved by the ‘method of sweeps’. Stream lines for ψ = 0·05, 0·2, 0·5 and magnitude of the vorticity vector z = 0·2 are plotted for R_eω = 0·1, 0·3, 0·5. Correction to the Stokes drag due to the rotation of fluid is calculated.     

A numerical study of pulsatile laminar flows in a pipe with a ring-type constriction

Numerical simulations have been carried out to study pulsatile laminar flows in a pipe with an axisymmetric ringtype constriction. Three types of pulsatile flows were investigated, namely a physiological flow, a pure sinusoidal flow and a non-zero mean velocity sinusoidal flow. The laminar flow governing equations were solved by the SIMPLE algorithm on a non-staggered grid and a modified Crank-Nicolson approximation was used to discretrize the momentum equations with respect to time. The maximum flow Reynolds numer (Re) is 100. The Womersley number (N_w) ranges from 0 to 50, with the corresponding Strouhal number (St) ranging from 0 to 3·98. The constriction opening ratio (d/D) and thickness ratio (h/D) are fixed at 0·5 and 0·1 respectively. Within the time period investigated, all these pulsatile flows include both forward and backward flows. The unsteady recirculation region and the recirculation points change in size and location with time. For N_w ≤ 1 and St≤ 1·56 x 10^?3 the three pulsatile flows have the same simple relation between the instantaneous flow rate and pressure loss (Δp) across the constriction and the pressure gradient in the axial direction (dp/dz) in the fully developed flow region. The phase angles between the flow rate and pressure loss and the pressure gradient are equal to zero. With increasing N_w and St, the phase angle between the flow rate and the dp/dz becomes larger and has its maximum value of 90° at N_w = 50 and St = 3·98. The three pulsatile flows also show different relations between the flow rate and the pressure gradient. The pure sinusoidal flow has the largest maximum pressure gradient and the non-zero mean velocity sinusoidal flow has the smallest. For larger N_w and St the fully developed velocity profiles in the fully developed flow region have a smaller velocity gradient along the radial direction in the central region. The maximum recirculation length increases for N_w ranging from 0 to 4·2, while this length becomes very small at N_w = 50 and St = 3·98. The deceleration tends to enlarge the recirculation region and this effect appears for N_w ≥ 3 and St ≥ 1·43×10^?2. Linear relations exist between the flow rate and the instantaneous maximum values of velocity, vorticity and shear stress.     

Numerical studies of slow viscous rotating flow past a sphere. III

The flow of steady incompressible viscous fluid rotating about the z-axis with angular velocity ω and moving with velocity u past a sphere of radius a which is kept fixed at the origin is investigated by means of a numerical method for small values of the Reynolds number Re_ω. The Navier–Stokes equations governing the axisymmetric flow can be written as three coupled non-linear partial differential equations for the streamfunction, vorticity and rotational velocity component. Central differences are applied to the partial differential equations for solution by the Peaceman–Rachford ADI method, and the resulting algebraic equations are solved by the ‘method of sweeps’. The results obtained by solving the non-linear partial differential equations are compared with the results obtained by linearizing the equations for very small values of Re_ω. Streamlines are plotted for Ψ = 0·05, 0·2, 0·5 for both linear and non-linear cases. The magnitude of the vorticity vector near the body, i.e. at z = 0·2, is plotted for Re_ω = 0·05, 0·24, 0·5. The correction to the Stokes drag as a result of rotation of the fluid is calculated.     

Force measurements on a circular cylinder in a Cryogenic Ludwieg-Tube using piezoelectric transducers

Force and pressure measurements were performed in a high Reynolds number facility, i.e. the Cryogenic Ludwieg-Tube (KRG). The balance based on multicomponent piezoelectric force transducers was applied totally in the cryogenic environment. The behaviour of the balance was tested ranging from ambient down to cryogenic conditions. As test cases, the flow around a circular cylinder was investigated at a sub- and supercritical Mach number. The highest possible Reynolds number in most nearly incompressible flow (Re = 5.8 × 10⁶) was achieved at the temperature of T = −150°C, the highest pressure possible, p ₀ = 10 bar, and the lowest attainable Mach number of Ma = 0.28. The results show that, in spite of the pulse operating mode of the tunnel, the steady and unsteady processes can be measured very well by means of a piezoelectric balance.     

Numerical modelling of confined flow past a cylinder of square cross-section at various orientations

Results are presented for the unsteady, two-dimensional flow and heat transfer due to a square obstruction of diameter d located asymmetrically between the parallel sliding walls of a channel with length-to-height ratio W/H = 6·44. Analysis is based on the numerical solution of spatially and temporally second-order accurate finite difference approximations of the transport equations expressed in curvilinear co-ordinates. Laminar, constant property flow is assumed for obstruction configurations in which the blockage ratio is d/H = 0·192, the nearest-wall distances are g/d = 0·2, 0·5 and 1, the orientation angles are α=0°, 10° and 20° and the Reynolds numbers are Re=100, 500, and 1000. Preparatory testing of the numerical procedure was performed for a variety of documented flows to verify its physiconumerical accuracy and obtain estimates of the residual grid-dependent uncertainties in the variables calculated. Heat transfer, drag and lift coefficients and Strouhal numbers for the present flow were finally calculated to within 4%–7% of their grid-dependent values using non-uniformly spaced grids consisting of (x=99, y=55) nodes. Above a critical value of the Reynolds number, which depends on the geometrical parameters, the flow is characterized by alternate vortex shedding from the obstruction top and bottom surfaces. Streamline, vorticity and particle streakline plots provide qualitative impressions of the unsteady vortical flow. Especially noteworthy are the extremes in the lift coefficient which ranges from large positive values for an obstruction with g/d=0·2 and α=10° to negative values for one with g/d=0·5 and α=0°. Both the drag and lift coefficients as well as the Strouhal number exhibit non-monotonic variations with respect to the parameters explored. Asymmetries in the obstruction location and orientation account for relatively large vortex-induced periodic variations in heat transfer, especially along the wall nearest the obstruction. Notable differences are also predicted for the heat transfer coefficients of the individual obstruction surfaces as a function of the orientation angle.     

Flow structures behind a vibrissa-shaped cylinder at different angles of attack: Complication on vortex-induced vibration

Complementary experimental studies have been conducted with a vibrissa-shaped cylinder at different angles of attack, through vortex-induced vibration (VIV) test in a wind tunnel, along with extensive measurements of wake dynamics in a water channel using time-resolved particle image velocimetry (TR-PIV). The VIV responses of an elastically mounted vibrissa-shaped cylinder are experimentally compared at various angles of attack in the range of θ = 0°–90°. At the reduced velocity of U₀/f₀D_h = 3–10 (f₀ being the system's natural frequency), the cross-flow displacement of the cylinders convincingly demonstrates that the vibrissa-shaped cylinder at a small angle of attack (θ ≤ 30°) is stable, and without appreciable displacement. Beyond θ = 30°, a prominent three-branched VIV response is readily identified, and increasing the angle of attack results in an upward shift of the synchronized region and a considerable intensification of the peak amplitude. Subsequently, TR-PIV measurements are made of the wake flow behind the vibrissa-shaped cylinder, to determine the spatio-temporally varying flow fields in two spanwise planes, i.e., the saddle and the nodal planes. Four systems with different angles of attack are chosen for comparison at Re_D = 1.8 × 10³, i.e., θ = 0°, 30°, 60° and 90°. In the two systems with θ = 0° and 30°, the wake regions feature weak velocity fluctuations in highly limited areas. However, increasing the angles of attack (to θ = 60° and 90°) gives rise to expanded recirculation zones, highly unstable flow reversals immediately behind the cylinder, and strengthened velocity fluctuations in the bulk wake regions. Cross-correlation of the fluctuating longitudinal velocities shows that at θ = 60° and 90° the energetic large-scale vortical structures form earlier, and they exert considerable influence on the near-wake fluid behind the cylinder. Finally, a sophisticated data-driven dynamic mode decomposition (DMD) process is used to extract the dominant unsteady structures in the four systems with different angles of attack. In the system with θ = 0°, two dominant DMD modes at frequencies St= 0.23 and 0.30 are identified in the saddle and the nodal planes, respectively, and those frequencies are St= 0.18 and 0.19 in the system with θ = 30°. The interaction between these dominant events at different frequencies tends to disrupt the formation of a strong vortex-shedding process. Therefore, the hydrodynamic force on the cylinder does not make a concerted contribution to suppressing the VIV behavior along the spanwise direction. In the systems with θ = 60° and 90°, the corresponding DMD modes exhibit much more synchronous, organized characteristics in the saddle and nodal planes, and unsteady events at the same frequencies are detected in both planes, reaching St = 0.14 (for θ = 60°) and 0.12 (for θ = 90°). These effects, along with the intensified vortex-shedding processes in the saddle and nodal planes, exert a concerted hydrodynamic force on the cylinder, causing it to start with an oscillatory state.     

A grid deformation technique for unsteady flow computations

A grid deformation technique is presented here based on a transfinite interpolation algorithm applied to the grid displacements. The method, tested using a two‐dimensional flow solver that uses an implicit dual‐time method for the solution of the unsteady Euler equations on deforming grids, is applicable to problems with time varying geometries arising from aeroelasticity and free surface marine problems. The present work is placed into a multi‐block framework and fits into the development of a generally applicable parallel multi‐block flow solver. The effect of grid deformation is examined and comparison with rigidly rotated grids is made for a series of pitching aerofoil test cases selected from the AGARD aeroelastic configurations for the NACA0012 aerofoil. The effect of using a geometric conservation law is also examined. Finally, a demonstration test case for the Williams aerofoil with an oscillating flap is presented, showing the capability of the grid deformation technique. Copyright © 2000 John Wiley & Sons, Ltd.     

Improvements to compressible Euler methods for low‐Mach number flows

In the present study improvements to numerical algorithms for the solution of the compressible Euler equations at low Mach numbers are investigated. To solve flow problems for a wide range of Mach numbers, from the incompressible limit to supersonic speeds, preconditioning techniques are frequently employed. On the other hand, one can achieve the same aim by using a suitably modified acoustic damping method. The solution algorithm presently under consideration is based on Roe's approximate Riemann solver [Roe PL. Approximate Riemann solvers, parameter vectors and difference schemes. Journal of Computational Physics 1981; 43 : 357–372] for non‐structured meshes. The numerical flux functions are modified by using Turkel's preconditioning technique proposed by Viozat [Implicit upwind schemes for low Mach number compressible flows. INRIA, Rapport de Recherche No. 3084, January 1997] for compressible Euler equations and by using a modified acoustic damping of the stabilization term proposed in the present study. These methods allow the compressible Euler equations at low‐Mach number flows to be solved, and they are consistent in time. The efficiency and accuracy of the proposed modifications have been assessed by comparison with experimental data and other numerical results in the literature. Copyright © 2000 John Wiley & Sons, Ltd.     

Periodic solution and chaotic behavior of a class of nonautonomic pendulum systems with large damping

In this paper the existence and uniqueness of the periodic solution is studied for a class of second order nonautonomic pendulum aystems x + ax + ψ(t)sinx=F(t) anil the parameter regions tor which the system in chaos is myestigated when ψ(t)=1-ελcosωt, F(t)=β +εμ(cosωt-ωsinωt) and the tamping coefficient a>0 is large. The result obtained generalize the corresponding conclusions of papers [1-8].     

The influence of acceleration and deceleration on shock wave movement on and around aerofoils in transonic flight

This paper examines the shock wave dynamics of a biconvex aerofoil in transonic flight during acceleration and retardation. The aerofoil has a cord length of 1 m and air at infinity is at 101.325 kPa and 300 K. Using Fluent as the CFD software, constant velocity (steady state) simulations were conducted at transonic Mach numbers. The aerofoil was then accelerated at 1041m/s² (106 g), starting at Mach 0.1, and decelerated at −1041m/s², starting at Mach 1.6, through the same range of Mach numbers using time-dependent (unsteady) simulations. Significant differences were found in the transonic region between the steady and the unsteady aerodynamic forces. Analysis of the flow field in this region showed that acceleration-dependent variations in the position of the shock wave on the surfaces of the aerofoil were the main reason for this. As very high accelerations were used in order to emphasize differences, which do not have many practical applications, simulations using accelerations lower than 9 g were also conducted in order to confirm the results. The acceleration-dependent behaviour of other shock waves around the aerofoil, such as the bow shock in front of the aerofoil and the trailing wave were also examined. The trailing wave followed behind the aerofoil changing position with different accelerations at the same Mach number.      

Off-design regimes of flow past waveriders based on isentropic compression flows

Supersonic off-design flow past waveriders on the M = 3 to 10 freestream Mach number range is numerically investigated. Configurations based on the flows behind plane shocks followed by isentropic flow compression are considered. The flow regimes are analyzed at the Mach numbers both smaller and greater than the design value M _d . The results are obtained by finite-volume solution of the Euler equations using higher-order Runge-Kutta TVD schemes.     

Computation of three-dimensional flow using the Euler equations and a multiple-grid scheme

The unsteady Euler equations are numerically solved using the finite volume one-step scheme recently developed by Ron-Ho Ni. The multiple-grid procedure of Ni is also implemented. The flows are assumed to be homo-enthalpic; the energy equation is eliminated and the static pressure is determined by the steady Bernoulli equation; a local time-step technique is used. Inflow and outflow boundaries are treated with the compatibility relations method of ONERA. The efficiency of the multiple-grid scheme is demonstrated by a two-dimensional calculation (transonic flow past the NACA 12 aerofoil) and also by a three-dimensional one (transonic lifting flow past the M6 wing). The third application presented shows the ability of the method to compute the vortical flow around a delta wing with leading-edge separation. No condition is applied at the leading-edge; the vortex sheets are captured in the same sense as shock waves. Results indicate that the Euler equations method is well suited for the prediction of flows with shock waves and contact discontinuities, the multiple-grid procedure allowing a substantial reduction of the computational time.     

Numerical analysis of wall pressure and heat flux fluctuations in shock-turbulent-boundary-layer interaction

The unsteady, compressible, Reynolds-averaged Navier-Stokes equations are solved numerically for an oblique shock-wave-induced turbulent boundary layer sepration. For the freestream Mach number 6 and the freestream Reynolds number 66·1 × 10⁶ m^?1, a time-dependent computation is performed, using MacCormack's explicit-implicit finite difference method with 82 × 42 grid points. A two-layer eddy viscosity turbulence model is employed in conjunction with a relaxation modification. Comparisons of the mean wall pressure and the mean heat transfer coefficient with the available experimental results are made and the evaluation of unsteady data for surface pressure and heat flux fluctuations is presented. It is found that the fluctuations in heat flux have qualitatively the same features as those of wall pressure but are different quantitatively.     

NUMERICAL STUDY OF EFFECTS OF PULSATILE AMPLITUDE ON UNSTEADY LAMINAR FLOWS IN RIGID PIPE WITH RING-TYPE CONSTRICTIONS

The effects of pulsatile amplitude on sinusoidal laminar flows through a rigid pipe with sharp-edged ring-type constrictions have been studied numerically. The parameters considered are: mean Reynolds number (Re) of the order of 100; Strouhal number (St) in the range 0·0–3·98; Womersley number (Nw) in the range 0·0–50·0. The pulsatile amplitude (A) varies in the range 0·0–2·0. The flow characteristics were studied through the pulsatile contours of streamline, vorticity, shear stress and isobars. Within a pulsatile cycle the relations between instantaneous flow rate (Q) and instantaneous pressure gradient (dp/dz) are observed to be elliptic. The relations between instantaneous flow rate (Q) and pressure loss (P_loss) are quadratic. Linear relations exist between instantaneous flow rate (Q) and maximum velocity, maximum vorticity and maximum shear stress. © by 1997 John Wiley & Sons, Ltd.     

A fast riemann solver with constant covolume applied to the random choice method

The Riemann problem for the unsteady one-dimensional Euler equations together with the constant-covolume equation of state is solved exactly. The solution is then applied to the random choice method to solve the general initial-boundary value problem for the Euler equations. The iterative procedure to find p*, the pressure between the acoustic waves, involves a single algebraic (non-linear) equation, all other quantities follow directly throughout the x–t plane, except within rarefaction fans where an extra iterative procedure is required. The solution is validated against existing exact results both directly and in conjunction with the random choice method.     

Prediction of oscillating thick cambered aerofoil aerodynamics by a locally analytic method

A complete first-order model and locally analytic solution method are developed to analyse the effects of mean flow incidence and aerofoil camber and thickness on the incompressible aerodynamics of an oscillating aerofoil. This method incorporates analytic solutions, with the discrete algebraic equations which represent the differential flow field equations obtained from analytic solutions in individual grid elements. The velocity potential is separated into steady and unsteady harmonic parts, with the unsteady potential further decomposed into circulatory and non-circulatory components. These velocity potentials are individually described by Laplace equations. The steady velocity potential is independent of the unsteady flow field. However, the unsteady flow is coupled to the steady flow field through the boundary conditions on the oscillating aerofoil. A body-fitted computational grid is then utilized. Solutions for both the steady and the coupled unsteady flow fields are obtained by a locally analytic numerical method in which locally analytic solutions in individual grid elements are determined. The complete flow field solution is obtained by assembling these locally analytic solutions. This model and solution method are shown to accurately predict the Theodorsen oscillating flat plate classical solution. Locally analytic solutions for a series of Joukowski aerofoils demonstrate the strong coupling between the aerofoil unsteady and steady flow fields, i.e. the strong dependence of the oscillating aerofoil aerodynamics on the steady flow effects of mean flow incidence angle and aerofoil camber and thickness.     

A high‐order discontinuous Galerkin solver for low Mach number flows

In this work, we present a high‐order discontinuous Galerkin method (DGM) for simulating variable density flows at low Mach numbers. The corresponding low Mach number equations are an approximation of the compressible Navier–Stokes equations in the limit of zero Mach number. To the best of the authors'y knowledge, it is the first time that the DGM is applied to the low Mach number equations. The mixed‐order formulation is applied for spatial discretization. For steady cases, we apply the semi‐implicit method for pressure‐linked equation (SIMPLE) algorithm to solve the non‐linear system in a segregated manner. For unsteady cases, the solver is implicit in time using backward differentiation formulae, and the SIMPLE algorithm is applied to solve the non‐linear system in each time step. Numerical results for the following three test cases are shown: Couette flow with a vertical temperature gradient, natural convection in a square cavity, and unsteady natural convection in a tall cavity. Considering a fixed number of degrees of freedom, the results demonstrate the benefits of using higher approximation orders. Copyright © 2015 John Wiley & Sons, Ltd.     

基于当地流活塞理论的气动弹性计算方法研究

发展了一种高效、高精度的超音速、高超音速非定常气动力计算 方法------基于定常CFD技术的当地流活塞理论. 运用当地流活塞理论计算非定常 气动力，耦合结构运动方程，实现超音速、高超音速气动弹性的时域模拟. 运用这 种方法计算了一系列非定常气动力算例和颤振算例，并和原始活塞理论、非定 常Euler方程结果作了比较. 由于局部地使用活塞理论假设，这种方法大大地克服 了原始活塞理论对飞行马赫数、翼型厚度和飞行迎角的 限制. 与非定常Euler方程方法相比，当地流活塞理论的效率很高.     

A high‐order discontinuous Galerkin method for all‐speed flows

In this article, we present a discontinuous Galerkin (DG) method designed to improve the accuracy and efficiency of steady solutions of the compressible fully coupled Reynolds‐averaged Navier–Stokes and k ? ω turbulence model equations for solving all‐speed flows. The system of equations is iterated to steady state by means of an implicit scheme. The DG solution is extended to the incompressible limit by implementing a low Mach number preconditioning technique. A full preconditioning approach is adopted, which modifies both the unsteady terms of the governing equations and the dissipative term of the numerical flux function by means of a new preconditioner, on the basis of a modified version of Turkel's preconditioning matrix. At sonic speed the preconditioner reduces to the identity matrix thus recovering the non‐preconditioned DG discretization. An artificial viscosity term is added to the DG discretized equations to stabilize the solution in the presence of shocks when piecewise approximations of order of accuracy higher than 1 are used. Moreover, several rescaling techniques are implemented in order to overcome ill‐conditioning problems that, in addition to the low Mach number stiffness, can limit the performance of the flow solver. These approaches, through a proper manipulation of the governing equations, reduce unbalances between residuals as a result of the dependence on the size of elements in the computational mesh and because of the inherent differences between turbulent and mean‐flow variables, influencing both the evolution of the Courant Friedrichs Lewy (CFL) number and the inexact solution of the linear systems. The performance of the method is demonstrated by solving three turbulent aerodynamic test cases: the flat plate, the L1T2 high‐lift configuration and the RAE2822 airfoil (Case 9). The computations are performed at different Mach numbers using various degrees of polynomial approximations to analyze the influence of the proposed numerical strategies on the accuracy, efficiency and robustness of a high‐order DG solver at different flow regimes. Copyright © 2014 John Wiley & Sons, Ltd.     

一种三阶有限体积法及其在欠膨胀射流激波结构数值模拟中的应用

以喷管出口欠膨胀射流为研究对象,在Lagrange坐标系下建立欠膨胀射流二维积分形式的流动方程。通过在单元交接面处进行三阶ENO(essentially nonoscillatory)格式插值,构造得到一种适用于求解该方程的三阶ENO有限体积法。采用该格式对一维Sod激波管算例和喷管出口欠膨胀射流进行数值计算。计算结果表明,该方法具有高精度、基本无振荡的特点,能很好地捕捉包含激波、滑移线以及三波交点等复杂流场波系结构。计算得到的波系结构中马赫盘的位置与实验结果吻合很好,相对误差小于1.1%。