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.     

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

The Navier–Stokes equations, which are the governing equations for a steady, viscous, incompressible fluid rotating about the z-axis with angular velocity ω, are linearized using the Oseen approximation. Two parameters, namely the Reynolds number Re = Ua/v and Re_ω = 2ωa²/v (the Reynolds number w.r.t. rotation), enter the linearized equations. These equations are solved by the Peaceman–Rachford ADI method and the resulting algebraic equations are solved by the SOR method. Streamlines are plotted and compared with the Oseen solution for the non-rotating case. The magnitude of the vorticity vector with increasing θ is also plotted.     

NUMERICAL SOLUTION OF THE SLOW TRANSLATION OF A SPHERE MOVING ALONG THE AXIS OF A ROTATING VISCOUS FLUID

The motion of a sphere along the axis of rotation of an incompressible viscous fluid that is rotating as a solid mass is investigated by means of numerical methods for small values of Reynolds numbers and moderate values of Taylor numbers. The Navier-Stokes equations governing the steady, axisymmetric, viscous flow can be written as three coupled, nonlinear, elliptic partial differential equations for the stream function, vorticity and rotational velocity component. Finite difference method is used for solving the governing equations. Second order derivatives are approximated by central differences and nonlinear terms are approximated by upwind differences. Results are presented mostly in the form of graphs of the streamlines and vorticity lines. When 1/ Ro > 2.2, separation occurs and reverse flow is obtained.     

Numerical study of vortex shedding from a rotating cylinder immersed in a uniform flow field

A numerical study is made of the unsteady two‐dimensional, incompressible flow past an impulsively started translating and rotating circular cylinder. The Reynolds number (Re) and the rotating‐to‐translating speed ratio (α) are two controlled parameters, and the influence of their different combinations on vortex shedding from the cylinder is investigated by the numerical scheme sketched below. Associated with the streamfunction (ψ)–vorticity (ω) formulation of the Navier–Stokes equations, the Poisson equation for ψ is solved by a Fourier/finite‐analytic, separation of variable approach. This approach allows one to attenuate the artificial far‐field boundary, and also yields a global conditioning on the wall vorticity in response to the no‐slip condition. As for the vorticity transport equation, spatial discretization is done by means of finite difference in which the convection terms are handled with the aid of an ENO (essentially non‐oscillatory)‐like data reconstruction process. Finally, the interior vorticity is updated by an explicit, second‐order Runge–Kutta method. Present computations fall into two categories. One with Re=10³ and α≤3; the other with Re=10⁴ and α≤2. Comparisons with other numerical or physical experiments are included. Copyright © 2000 John Wiley & Sons, Ltd.     

Asymptotic Behavior of Solutions of the Compressible Navier-Stokes Equations on the Half Space

Asymptotic behavior of solutions of the compressible Navier-Stokes equations on the half space R ⁿ₊(n≥2) is considered around a given constant equilibrium. A solution formula for the linearized problem is derived, and L^p estimates for solutions of the linearized problem are obtained for 2≤p≤∞. It is shown that, as in the case of the Cauchy problem, the leading part of the solution of the linearized problem is decomposed into two parts, one that behaves like diffusion waves and the other one purely diffusively. There, however, are some aspects different from the Cauchy problem, especially in considering spatial derivatives. It is also shown that the solution of the linearized problem approaches for large times the solution of the nonstationary Stokes problem in some L^p spaces; and, as a result, a solution formula for the nonstationary Stokes problem is obtained. Large-time behavior of solutions of the nonlinear problem is then investigated in L^p spaces for 2≤p≤∞ by applying the results on the linearized analysis and the weighted energy method. The results indicate that there may be some nonlinear interaction phenomena not appearing in the Cauchy problem.     

On coupling the Reynolds‐averaged Navier–Stokes equations with two‐equation turbulence model equations

Two methods for coupling the Reynolds‐averaged Navier–Stokes equations with the q–ω turbulence model equations on structured grid systems have been studied; namely a loosely coupled method and a strongly coupled method. The loosely coupled method first solves the Navier–Stokes equations with the turbulent viscosity fixed. In a subsequent step, the turbulence model equations are solved with all flow quantities fixed. On the other hand, the strongly coupled method solves the Reynolds‐averaged Navier–Stokes equations and the turbulence model equations simultaneously. In this paper, numerical stabilities of both methods in conjunction with the approximated factorization‐alternative direction implicit method are analysed. The effect of the turbulent kinetic energy terms in the governing equations on the convergence characteristics is also studied. The performance of the two methods is compared for several two‐ and three‐dimensional problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Vorticity and Regularity for Viscous Incompressible Flows under the Dirichlet Boundary Condition. Results and Related Open Problems

In reference [7] it is proved that the solution of the evolution Navier–Stokes equations in the whole of R ³ must be smooth if the direction of the vorticity is Lipschitz continuous with respect to the space variables. In reference [5] the authors improve the above result by showing that Lipschitz continuity may be replaced by 1/2-H?lder continuity. A central point in the proofs is to estimate the integral of the term (ω · ∇)u · ω, where u is the velocity and ω = ∇ × u is the vorticity. In reference [4] we extend the main estimates on the above integral term to solutions under the slip boundary condition in the half-space R ₊³. This allows an immediate extension to this problem of the 1/2-H?lder sufficient condition. The aim of these notes is to show that under the non-slip boundary condition the above integral term may be estimated as well in a similar, even simpler, way. Nevertheless, without further hypotheses, we are not able now to extend to the non slip (or adherence) boundary condition the 1/2-H?lder sufficient condition. This is not due to the “nonlinear" term (ω · ∇)u · ω but to a boundary integral which is due to the combination of viscosity and adherence to the boundary. On the other hand, by appealing to the properties of Green functions, we are able to consider here a regular, arbitrary open set Ω.      

Numerical analysis of the rotating viscous flow approaching a solid sphere

A numerical simulation is performed to investigate the flow induced by a sphere moving along the axis of a rotating cylindrical container filled with the viscous fluid. Three‐dimensional incompressible Navier–Stokes equations are solved using a finite element method. The objective of this study is to examine the feature of waves generated by the Coriolis force at moderate Rossby numbers and that to what extent the Taylor–Proudman theorem is valid for the viscous rotating flow at small Rossby number and large Reynolds number. Calculations have been undertaken at the Rossby numbers (R_o) of 1 and 0.02 and the Reynolds numbers (R_e) of 200 and 500. When R_o=O(1), inertia waves are exhibited in the rotating flow past a sphere. The effects of the Reynolds number and the ratio of the radius of the sphere and that of the rotating cylinder on the flow structure are examined. When R_o ? 1, as predicted by the Taylor–Proudman theorem for inviscid flow, the so‐called ‘Taylor column’ is also generated in the viscous fluid flow after an evolutionary course of vortical flow structures. The initial evolution and final formation of the ‘Taylor column’ are exhibited. According to the present calculation, it has been verified that major theoretical statement about the rotating flow of the inviscid fluid may still approximately predict the rotating flow structure of the viscous fluid in a certain regime of the Reynolds number. Copyright © 2004 John Wiley & Sons, Ltd.     

Efficient preconditioning techniques for finite‐element quadratic discretization arising from linearized incompressible Navier–Stokes equations

We develop an efficient preconditioning techniques for the solution of large linearized stationary and non‐stationary incompressible Navier–Stokes equations. These equations are linearized by the Picard and Newton methods, and linear extrapolation schemes in the non‐stationary case. The time discretization procedure uses the Gear scheme and the second‐order Taylor–Hood element P₂?P₁ is used for the approximation of the velocity and the pressure. Our purpose is to develop an efficient preconditioner for saddle point systems. Our tools are the addition of stabilization (penalization) term r?(div(·)), and the use of triangular block matrix as global preconditioner. This preconditioner involves the solution of two subsystems associated, respectively, with the velocity and the pressure and have to be solved efficiently. Furthermore, we use the P₁?P₂ hierarchical preconditioner recently proposed by the authors, for the block matrix associated with the velocity and an additive approach for the Schur complement approximation. Finally, several numerical examples illustrating the good performance of the preconditioning techniques are presented. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

The Fundamental Solution of the Linearized Navier–Stokes Equations for Spinning Bodies in Three Spatial Dimensions – Time Dependent Case

Explicit formulae for the fundamental solution of the linearized time dependent Navier–Stokes equations in three spatial dimensions are obtained. The linear equations considered in this paper include those used to model rigid bodies that are translating and rotating at a constant velocity. Estimates extending those obtained by Solonnikov in [23] for the fundamental solution of the time dependent Stokes equations, corresponding to zero translational and angular velocity, are established. Existence and uniqueness of solutions of these linearized problems is obtained for a class of functions that includes the classical Lebesgue spaces L^p(R³), 1 < p < ∞. Finally, the asymptotic behavior and semigroup properties of the fundamental solution are established.     

Motion of a sphere in an electrically conducting rotating fluid

The combined effect of rotation and magnetic field is investigated for the axisymmetric flow due to the motion of a sphere in an inviscid, incompressible electrically conducting fluid having uniform rotation far upstream. The steady-state linearized equations contain a single parameter α=1/2βR _m, β being the magnetic pressure number and R _m the magnetic Reynolds number. The complete solution for the flow field and magnetic field is obtained and the distribution of vorticity and current density is found. The induced vorticity is O(α⁴) and the current density is O(R _m) on the sphere.     

Direct Numerical Simulation of Transitional Turbulent Flow in a Closed Rotor-Stator Cavity

The transitional turbulent regime in confined flow between a rotating and a stationary disc is studied using direct numerical simulation. Besides its fundamental importance as a three-dimensional prototype flow, such flows frequently arise in many industrial devices, especially in turbomachinary applications. The present contribution extends the DNS simulation into the turbulent flow regime, to a rotational Reynolds number Re =3 × 10⁵. An annular rotor-stator cavity of radial extension ΔR and height H, is considered with L = 4.72(L = ΔR/H) and Rm = 2.33 (Rm = (R ₁+ R ₀)/ΔR). The direct numerical simulation is performed by integrating the time-dependent Navier–Stokes equations until a statistically steady state is reached. A three-dimensional spectral method is used with the aim of providing both very accurate instantaneous fields and reliable statistical data. The instantaneous quantities are analysed in order to enhance our knowledge of the physics of turbulent rotating flows. Also, the results have been averaged so as to provide target turbulence data for any subsequent modelling attempts at reproducing the flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

Unsteady supersonic flow over an oscillatory aerofoil using a second-order TVD scheme

The unsteady flow over an oscillatory NACA0012 aerofoil has been simulated by the calculation with Euler equations. The equations are discretized by an implicit Euler in time, and a second-order space-accurate TVD scheme based on flux vector splitting with van Leer's limiter. Modified eigenvalues are proposed to overcome the slope discontinuities of split eigenvalues at Mach = 0·0 and ± 1·0, and to generate a bow shock in front of the aerofoil. A moving grid system around the aerofoil is generated by Sorenson's boundary fitted co-ordinates for each time step. The calculations have been done for two angles of attack θ = 5·0° sin (ωt) and θ = 3·0° + 3·0° sin (ωt) for the free-stream Mach numbers 2·0 and 3·0. The results show that pressure and Mach cells flow along characteristic lines. To examine unsteady effects, the responses of wall pressure and normal force coefficients are analysed by a Fourier series expansion.     

NEAR-WAKE FLOW STRUCTURE OF A CYLINDER WITH A HELICAL SURFACE PERTURBATION

The near-wake of a circular cylinder having a helical wire pattern about its surface is characterized using a technique of high-image-density velocimetry. Patterns of vorticity in three orthogonal planes show substantial influence of a wire having a diameter an order of magnitude smaller than the cylinder diameter. The distinctive patterns of vorticity in these three planes are associated with lack of formation of large-scale Kármán-like clusters of vorticity (ω_z) in the near-wake region of the cylinder. The instantaneous structure of the separating spanwise vorticity (ω_z) layers on either side of the cylinder involve small-scale concentrations of vorticity analogous to the well-known Kelvin–Helmholtz vortices from a smooth cylinder. Moreover, a dual vorticity layer, i.e., two adjacent layers of like vorticity (ω_z), can form from one side of the cylinder. Along the span of the cylinder, distributions of instantaneous velocity and transverse vorticity (ω_y) show a spatially periodic sequence of wake-like patterns, each of which has features in common with the very near-wake of a two-dimensional bluff body, including a large velocity defect bounded by vorticity layers with embedded small-scale vorticity concentrations. In the cross-flow plane of the wake, patterns of streamwise vorticity (ω_x) show small-scale, counter-rotating pairs of vorticity concentrations (ω_x) emanating from the inclined helical perturbation, rather than isolated concentrations of vorticity of like sign, which would indicate single streamwise vortices. All of the aforementioned patterns of small-scale vorticity concentrations are scaled according to the local wake width/local pitch of the helical wire pattern in the respective plane of observation.     

On the Navier�CStokes equations with rotating effect and prescribed outflow velocity

We consider the equations of Navier–Stokes modeling viscous fluid flow past a moving or rotating obstacle in \mathbb R^d{\mathbb R^d} subject to a prescribed velocity condition at infinity. In contrast to previously known results, where the prescribed velocity vector is assumed to be parallel to the axis of rotation, in this paper we are interested in a general outflow velocity. In order to use L ^p -techniques we introduce a new coordinate system, in which we obtain a non-autonomous partial differential equation with an unbounded drift term. We prove that the linearized problem in \mathbb R^d{\mathbb R^d} is solved by an evolution system on L^p_s(\mathbb R^d){L^p_{\sigma}(\mathbb R^d)} for 1 < p < ∞. For this we use results about time-dependent Ornstein–Uhlenbeck operators. Finally, we prove, for p ≥ d and initial data u₀ ? L^p_s(\mathbb R^d){u_0\in L^p_{\sigma}(\mathbb R^d)}, the existence of a unique mild solution to the full Navier–Stokes system.     

Direct numerical simulation of the amplification of a 2D temporal disturbance in plane Poiseuille flow

A direct numerical scheme is developed to study the temporal amplification of a 2D disturbance in plane Poiseuille flow. The transient non-linear Navier–Stokes equations are applied in a region of wavelength moving with the wave propagation speed. The complex amplitude involved in the perturbation functions is considered as the initial input of the non-linear stability equations. In this study a fully implicit finite difference scheme with five points in the flow direction and three points in the normal direction is developed so that numerical simulation of the amplification of a two-dimensional temporal disturbance in plane Poiseuille flow can be investigated. The growth and decay of the disturbance with time are presented and neutral stability curves which are in good agreement with existing solutions can be determined. The critical conditions as a function of the amplitude A₀ of the disturbance are presented. Fixing the wavelength, the Navier–Stokes equations are solved up to Re=10,000 a friction factor increasing with Reynolds number is observed. The 2D non-linear behaviour of the streamfunction, vorticity and velocity components at Re=10,000 are also exhibited. © 1998 John Wiley & Sons, Ltd.     

Translation of a sphere in a rotating viscous fluid: A numerical study

The translation of a sphere moving along the axis of a rotating viscous fluid is studied by the finite difference method at moderate Reynolds (up to R = 500) and Taylor (up to T = 100) numbers. Suppression of the separation is observed with increasing rotation parameter T. The drag coefficient is also presented. It is observed that the drag coefficient is less than that with no rotation in the range 0<N<0·7, where N = 2T/R is the inverse Rossby number. The same phenomenon was observed experimentally by Maxworthy in the range 0<N<0·75±0·03.     

Turbulence in a skewed three‐dimensional wall‐bounded shear flow: effect of mean vorticity on structure modification

Mean‐flow three‐dimensionalities affect both the turbulence level and the coherent flow structures in wall‐bounded shear flows. A tailor‐made flow configuration was designed to enable a thorough investigation of moderately and severely skewed channel flows. A unidirectional shear‐driven plane Couette flow was skewed by means of an imposed spanwise pressure gradient. Three different cases with 8°, 34°and 52°skewing were simulated numerically and the results compared with data from a purely two‐dimensional plane Couette flow. The resulting three‐dimensional flow field became statistically stationary and homogeneous in the streamwise and spanwise directions while the mean velocity vector V and the mean vorticity vector Ω remained parallel with the walls. Mean flow profiles were presented together with all components of the Reynolds stress tensor. The mean shear rate in the core region gradually increased with increasing skewing whereas the velocity fluctuations were enhanced in the spanwise direction and reduced in the streamwise direction. The Reynolds shear stress is known to be closely related to the coherent flow structures in the near‐wall region. The instantaneous and ensemble‐averaged flow structures were turned by the skewed mean flow. We demonstrated for the medium‐skewed case that the coherent structures should be examined in a coordinate system aligned with V to enable a sound interpretation of 3D effects. The conventional symmetry between Case 1 and Case 2 vortices was broken and Case 1 vortices turned out to be stronger than Case 2. This observation is in conflict with the common understanding on the basis of the spanwise (secondary) mean shear rate. A refined model was proposed to interpret the structure modifications in three‐dimensional wall‐flows. What matters is the orientation of the mean vorticity vector Ω relative to the vortex vorticity vector ω _v, that is, the sign of Ω · ω _v. In the present situation, Ω · ω _v > 0 for the Case 1 vortices causing a strengthening relative to the Case 2 vortices. Copyright © 2011 John Wiley & Sons, Ltd.     

Finite-element analysis of fluid flow and heat transfer for staggered bundles of cylinders in cross flow

The two-dimensional Navier-Stokes equations and the energy equation governing steady laminar incompressible flow are solved by a penalty finite-element model for flow across finite depth, five-row deep, staggered bundles of cylinders. Pitch to diameter ratios of 1·5 and 2·0 are considered for cylinders in equilateral triangular and square arrangements. Reynolds numbers studied range from 100 to 400, and a Prandtl number of 0·7 is used. Velocity vector fields, streamline patterns, vorticity, pressure and temperature contours, local and average Nusselt numbers, pressure and shear stress distributions around the cylinder walls and drag coefficients are presented. The results obtained agree well with available experimental and numerical data.