Unsteady RANS method for surface ship boundary layer and wake and wave field

Results are reported of an unsteady Reynolds‐averaged Navier–Stokes (RANS) method for simulation of the boundary layer and wake and wave field for a surface ship advancing in regular head waves, but restrained from body motions. Second‐order finite differences are used for both spatial and temporal discretization and a Poisson equation projection method is used for velocity–pressure coupling. The exact kinematic free‐surface boundary condition is solved for the free‐surface elevation using a body‐fitted/free‐surface conforming grid updated in each time step. The simulations are for the model problem of a Wigley hull advancing in calm water and in regular head waves. Verification and validation procedures are followed, which include careful consideration of both simulation and experimental uncertainties. The steady flow results are comparable to other steady RANS methods in predicting resistance, boundary layer and wake, and free‐surface effects. The unsteady flow results cover a wide range of Froude number, wavelength, and amplitude for which first harmonic amplitude and phase force and moment experimental data are available for validation along with frequency domain, linear potential flow results for comparisons. The present results, which include the effects of turbulent flow and non‐linear interactions, are in good agreement with the data and overall show better capability than the potential flow results. The physics of the unsteady boundary layer and wake and wave field response are explained with regard to frequency of encounter and seakeeping theory. The results of the present study suggest applicability for additional complexities such as practical ship geometry, ship motion, and maneuvering in arbitrary ambient waves. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical Investigation of Unsteady Transitional Flows in Turbomachinery Components Based on a RANS Approach

This paper presents results of the numerical simulation of periodically unsteady flows with focus on turbomachinery applications. The unsteady CFD solver used for the simulations is based on the Reynolds averaged Navier–Stokes equations. The numerical scheme applies an extended version of the Spalart–Allmaras one-equation turbulence model coupled with a transition correlation. The first example of validation consists of boundary layer flow with separation bubble on a flat plate, both under steady and periodically unsteady main flow conditions. The investigation includes a variation of the major parameters Strouhal number, amplitude, and Reynolds number. The second, more complex test case consists of the flow through a cascade of turbine blades which is influenced by wakes periodically passing over the cascade. The computations were carried out for two different blade loadings. The results of the numerical simulations are discussed and compared with experimental data in detail. Special emphasis is given to the investigation of boundary layers with regard to transition, separation and reattachment under the influence of main flow unsteadiness. This revised version was published online in July 2006 with corrections to the Cover Date.     

Unsteady incompressible flows past two cylinders in tandem and staggered arrangements

A stabilized finite element formulation is employed to study incompressible flows past a pair of cylinders at Reynolds numbers 100 and 1000 in tandem and staggered arrangements. Computations are carried out for three sets of cylinder arrangements. In the first two cases the cylinders are arranged in tandem and the distance between their centres is 2·5 and 5·5 diameters. The third case involves the two cylinders in staggered arrangement. The distance between their centres along the flow direction is 5·5 diameters, while it is 0·7 diameter in the transverse direction. The results are compared with flows past a single cylinder at corresponding Reynolds numbers and with experimental observations by other researchers. It is observed that the qualitative nature of the flow depends strongly on the arrangement of cylinders and the Reynolds number. In all cases, when the flow becomes unsteady, the downstream cylinder, which lies in the wake of the upstream one, experiences very large unsteady forces that may lead to wake-induced flutter. The Strouhal number, based on the dominant frequency in the time history of the lift coefficient, for both cylinders attains the same value. In some cases, even though the near wake of the two cylinders shows temporal periodicity, the far wake does not. © 1997 John Wiley & Sons, Ltd.     

The complementary RANS equations for the simulation of viscous flows

A complementary set of Reynolds‐averaged Navier–Stokes (RANS) equations has been developed for steady incompressible, turbulent flows. The method is based on the Helmholtz decomposition of the velocity vector field into a viscous and a potential components. In the complementary RANS solver a potential solution coexists with a viscous solution with the purpose of contributing to a fastest decay of the viscous solution in the far field. The proposed complementary RANS equations have been validated for steady laminar and turbulent flows. The computational results show that the complementary RANS solver is able to produce less grid‐dependent solutions than a conventional RANS solver. Copyright © 2005 John Wiley & Sons, Ltd.     

流向振荡柱体尾流控制研究进展

详细介绍了近几年采用尾部喷射、隔离板和小窄条控制件等3 种方法对流向振荡柱体绕流旋涡脱落的抑制情况. 在研究范围内存在非锁频和3种锁频旋涡脱落模式. 风洞实验表明, 尾部喷射对这4 种模式都有抑制效果,窄条控制件对非锁频和2种锁频模式具有抑制效果, 而隔离板仅对非锁频和1 种锁频模式有效. 在不同流动和振荡条件下找出了每种方法的有效控制区, 研究了减阻和减少脉动力的效果, 并探讨了控制机理.      

Comparison of the RANS and PDF methods for air-particle flows

Two simulation methods, namely Reynolds-Averaged Navier–Stokes (RANS) equations, and Probability Distribution Function (PDF) are currently widely used for the modeling of multiphase flows. These two approaches are supplemented with appropriate closure equations that take into account all the pertinent forces and interaction effects on the solid particles, such as: particle–turbulence interactions; turbulence modulation; particle–particle interactions; particle–wall interactions; gravitation, drag and lift forces. The two methods have been used in order to simulate the turbulent particulate flow in upward pipes. The flow domain in all cases was a cylindrical pipe and the computations were carried for upward pipe flow. Monodisperse as well as polydisperse mixtures of particles have been considered. In general, the average velocity results obtained from the two methods are in close agreement, because the methods predict well the average velocity distribution of the carrier fluid as well as the solids. Thus, the differences in the average axial velocities predicted by the methods are not substantial. Differences in the turbulence intensity are more significant. A comparison of the numerical results obtained shows the relative importance of retaining the diffusion terms in both the axial and radial directions in the RANS method. Also the comparisons of the results show the relative effect of the lift forces in the distribution of solid particles.     

An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows

In this paper, an immersed boundary (IB) method is developed to simulate compressible turbulent flows governed by the Reynolds‐averaged Navier‐Stokes equations. The flow variables at the IB nodes (interior nodes in the immediate vicinity of the solid wall) are evaluated via linear interpolation in the normal direction to close the discrete form of the governing equations. An adaptive wall function and a 2‐layer wall model are introduced to reduce the near‐wall mesh density required by the high resolution of the turbulent boundary layers. The wall shear stress modified by the wall modeling technique and the no‐penetration condition are enforced to evaluate the velocity at an IB node. The pressure and temperature at an IB node are obtained via the local simplified momentum equation and the Crocco‐Busemann relation, respectively. The SST k ? ω and S‐A turbulence models are adopted in the framework of the present IB approach. For the Shear‐Stress Transport (SST) k ? ω model, analytical solutions in near‐wall region are utilized to enforce the boundary conditions of the turbulence equations and evaluate the turbulence variables at an IB node. For the S‐A model, the turbulence variable at an IB node is calculated by using the near‐wall profile of the eddy viscosity. In order to validate the present IB approach, numerical experiments for compressible turbulent flows over stationary and moving bodies have been performed. The predictions show good agreements with the referenced experimental data and numerical results.     

Regimes of vortex shedding from an in-line oscillating circular cylinder in the uniform flow

A finite volume method for the time dependent viscous incompressible flow around an in-line oscillating circular cylinder at Reynolds number of 200, 855 is presented in this paper. The Navier-Stokes equations in a finite volume form are solved with a moving grid system, based on a time dependent coordinate transformation. To investigate the vortex-shedding characteristics behind the circular cylinder and the effects of Reynolds number and other non-dimensional parameters such as reduced amplitude and reduced frequency， several numerical schemes have been tested with different amplitude and frequency close to Sto and a harmonic at each Reynolds number. Present numerical results indicate several types of vortex shedding mode which is known mainly depending on the reduced frequency and also the reduced amplitude, which is called synchronization or lock-on.     

动失速型非定常分离涡结构的控制

本文对动失速型非定常分离涡结构的控制方法，在低速风洞中应用相平均测压技术进行了实验研究。在二元平板模型中部安装一作俯仰振荡的扰流板产生动失速型分离涡，在其上游安装另一用作控制的小扰流板。实验结果表明，应用前置的振荡小扰流板可影响并改变动失速分离涡的强度和对流特性。在最有利的控制相位下，涡吸力峰可降低４８％，涡对流时间可以推迟０．１９周期。对于间歇式振荡扰流板，采用相位提前控制方式比相位滞后控制方式更有效。     

Modeling the flow past an oscillating airfoil by the method of viscous vortex domains

The Lagrangian vortex method for solving the Navier-Stokes equations is applied for numerically modeling the unsteady flow past a wing airfoil executing angular oscillations in a viscous incompressible flow. Formulas relating the unsteady forces on the airfoil and the vorticity field are derived. The calculated results are compared with the experimental data for the NACA-0012 airfoil executing harmonic oscillations in an air flow at the Reynolds number Re = 4.4 × 10⁴.     

Control of vortex shedding behind circular cylinder for flows at low Reynolds numbers

It has been observed by researchers in the past that vortex shedding behind circular cylinders can be altered, and in some cases suppressed, over a limited range of Reynolds numbers by proper placement of a second, much smaller, ‘control’ cylinder in the near wake of the main cylinder. Results are presented for numerical computations of some such situations. A stabilized finite element method is employed to solve the incompressible Navier–Stokes equations in the primitive variables formulation. At low Reynolds numbers, for certain relative positions of the main and control cylinder, the vortex shedding from the main cylinder is completely suppressed. Excellent agreement is observed between the present computations and experimental findings of other researchers. In an effort to explain the mechanism of control of vortex shedding, the streamwise variation of the pressure coefficient close to the shear layer of the main cylinder is compared for various cases, with and without the control cylinder. In the cases where the vortex shedding is suppressed, it is observed that the control cylinder provides a local favorable pressure gradient in the wake region, thereby stabilizing the shear layer locally. Copyright © 2001 John Wiley & Sons, Ltd.     

Flow past a transversely oscillating square cylinder in free stream at low Reynolds numbers

This paper reports simulation results for free‐stream flow past an oscillating square cylinder at Re=100 and 150, for oscillating‐to‐natural‐shedding frequency ratios of 0.5?f_r?3.0 at a fixed oscillation amplitude of 0.2 of the cylinder width. The transformed governing equations are solved in a non‐inertial frame of reference using the finite volume technique. The ‘lock‐in’ phenomena, where the vortex shedding becomes one with the oscillation frequency, is observed near the natural shedding frequency (f_r≈1). Beyond the synchronization band, downstream recovery of the wake to its stationary (natural) state (frequency) is observed in cross‐stream velocity spectra. At higher forcing frequencies, a phase lag between the immediate and the far wake results in a shear layer having multi‐polar vortices. A ‘Vortex‐switch’ accompanied by a change in the direction of energy transfer is identified at the ‘lock‐in’ boundaries. The variation of aerodynamic forces is noticed to be different in the lock‐in regime. The velocity phase portrait in the far wake revealed a chaotic state of flow at higher excitation though a single (natural) frequency appears in the spectra. Copyright © 2008 John Wiley & Sons, Ltd.     

Unsteady theory of the hydrofoil of finite span

A three-dimensional unsteady flow around the hydrofoil has been dealt with in this paper. Under the assumption of infinitesimal wave and thin hydrofoil, analytic solution of perturbation velocity potential induced by the unsteady flow of a hydrofoil which is simplified as a series of rectangular vortex rings is achieved and then the unsteady forces can be evaluated. The paper had been accepted by the X VI th International Congress of IUTAM, Lyngby, Denmark, Aug. 19–25, 1985.     

Hybrid Unsteady RANS and PDF Method for Turbulent Non-Reactive and Reactive Flows

A hybrid unsteady Reynolds-averaged numerical simulation (U-RANS) and probability density function (PDF) method is developed for turbulent non-reactive and reactive flows. The resulting modeled equations are solved by a consistent hybrid finite volume and Lagrangian Monte-Carlo particle method. Both turbulent non-reactive and reactive flows in a rectangular channel containing a triangular-shaped bluff-body are simulated. One-step and two-step mechanisms for propane/air combustion are used for the reactive case. The time-averaged results are compared with both experimental data and numerical results from the literature using large eddy simulation (LES) and steady RANS. The results of the present method are in good agreement with the experimental data, and they improve the numerical results available in the literature.     

串列双圆柱绕流问题的数值模拟

本文运用有限体积方法,对绕串列放置的双圆柱的二维不可压缩流动进行了数值计算。为研究两圆柱不同间距对圆柱相互作用和尾流特征的影响,选取间距比Ｌ／Ｄ（Ｌ为两圆柱中心间的距离,Ｄ为圆柱直径）在１．５～５．０之间每隔０．５共八个有代表性的间距进行了计算模拟。计算均在Ｒｅ＝２００条件下进行。计算结果表明：对该绕流问题,流动特征在很大程度上取决于间距的大小。且间距存在一临界值,间距比从小于临界值变化到大于临界     

Measurements of pulsating and oscillating flows of non-Newtonian fluids through concentric and eccentric cylinders

The unsteady flow of non-Newtonian fluids through concentric and eccentric cylinders was investigated experimentally. Two experiments were carried out; one was pulsating flow and the other was flow under a constant pressure gradient with the inner cylinder oscillating longitudinally. The flow enhancement was examined and its dependence on the frequency of the oscillations and the eccentricity of the apparatus was determined.     

Characteristic Regimes of Transitional Separation Bubbles in Unsteady Flow

This experimental investigation deals with transition phenomena of a separated boundary layer under unsteady inlet flow conditions. The main purpose of this investigation is to understand the influence of the rotor-stator interaction in turbomachinery on the subsequent, highly loaded boundary layer. The research project is divided into two phases. In the first phase, which has been completed recently, only the variation of mean velocity caused by upstream blades was simulated in the experiments while the free-stream turbulence intensity was retained at a constant low level. The experiments are carried out in an Eifel-type wind tunnel to investigate the laminar separated boundary layer of a flat plate under oscillating inlet conditions. The adverse pressure gradient, similar to that of turbomachines, is generated by the contoured upper wall. The unsteadiness is produced by a rotating flap located downstream of the test section. The reduced frequency, the amplitude and the mean Reynolds number are varied to simulate the conditions prevailing in turbomachines. In addition to the Kelvin–Helmholtz instability of the separated shear layer, a lower frequency instability was observed. This is frequently referred to as `free shear layer flapping' and results in two distinctly different ways of re-attachment, depending primarily on the Reynolds number. For low momentum thickness Reynolds numbers at the separation point, large-scale vortices locked to the frequency of the unsteady main flow are identified. They originate nearly at the top of the separation bubble and are ejected downstream. A fully turbulent boundary layer develops after these vortices mix out. For higher Reynolds numbers, transition is completed within a short length of the free shear layer and there-attachment region. The characteristic momentum thickness Reynolds number separating these two regimes in unsteady flow is about 125. The Strouhal number (reduced frequency) does not appear to have any significant effect. Based on the experimental results, this behaviour is discussed in some detail. This revised version was published online in July 2006 with corrections to the Cover Date.     

Periodic solution of turbulent oscillating channel flows

The time-dependent turbulent Navier–Stokes equations are solved numerically by a finite element method with an algebraic eddy viscosity model (Baldwin–Lomax formulation) for oscillating turbulent channel flows. The method of averaging is used to analyse the resulting periodic motion of the fluid. Numerical results are obtained for various Strouhal numbers and relative amplitudes. A comparison is made between the numerical and published experimental results. It appears that for low relative amplitudes in a certain range of frequencies the agreement is satisfactory.