圆柱管道内运动的小球受到的黏滞阻力的数值分析

为了更好地运用落球法测量研究流体的黏滞系数,研究小球在黏性流体中下落的受力情况,本文对小球在充满黏性流体的圆柱管道的下落过程进行分析。利用COMSOL4.4仿真模拟,建立了合理的仿真模型,并分析了小球受到的黏滞阻力与小球的大小、下落位置的关系。结果表明:选择速度项二阶近似、压强项一阶近似的离散化方法,可以得到和理论值非常相符的仿真结果;当下落过程中小球球心始终在圆柱轴线上时,小球受到的黏滞阻力相对于Stokes力的修正系数,是小球半径与圆柱管道半径的比例函数,本文得到了更大范围的符合理论解的修正系数;当下落过程中小球的球心偏离圆柱轴线时,对于同样大小的小球,黏滞阻力、压强力、黏性力均随着球心到轴线的距离先减小后增大,且具有不同的极小值点。     

Simulation of the Action of a Permeable Barrier on an Ideal Incompressible Channel Flow

Plane ideal incompressible flow in a rectangular channel partitioned by a thin permeable barrier (lattice) is considered. In flowing through the lattice the stream suddenly (jumpwise) changes direction and loses energy. The flow is assumed to be vortical; the vorticity is discontinuous on the lattice. A mathematical formulation of the problem for the stream function is proposed in the form of a nonlinear elliptic equation with coefficients discontinuous on the lattice line. A numerical solution is constructed using the finite-element iteration method. The results of the numerical simulation show how the flow velocity profile in the channel can be controlled by means of permeable barriers.     

Features of the supercritical transition in the wake of a circular cylinder

The features of the wake behind a uniform circular cylinder atRe=200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations by solving 3-D incompressible Navier-Stokes equations using mixed spectral-spectral-element method. The high-order splitting algorithm based on the mixed stiffly stable scheme is employed in the time discretization. Due to the nonlinear evolution of the secondary instability of the wake, the spanwise modes with different wavelengths emerge. The spanwise characteristic length determines the transition features and global properties of the wake. The existence of the spanwise phase difference of the primary vortices shedding is confirmed by Fourier analysis of the time series of the spanwise vorticity and attributed to the dominant spanwise mode. The spatial energy distributions of various modes and the velocity profiles in the near wake are obtained. The numerical results indicate that the near wake is in 3-D quasi-periodic laminar state with transitional behaviors at this supercritical Reynolds number. The project supported by the State Key Fundamental Research Project of “Large Scale Scientific Computation Research” (G199903281)     

A numerical study of the wake behind a blunt‐trailing‐edged body. Part 2: the topology of the flow

This project focuses on the wake behind a two‐dimensional blunt‐trailing‐edged body. The body has an elliptical nose followed by a straight section that ends in a blunt base. For this geometry, the boundary layer separation points are defined and fixed. This permits a better assessment of the vital influence of the boundary layers upon the wake. It is hoped that the study of this special geometry will contribute to a better understanding of the wake flow. The main aim of this article is to study the topology of the wake by means of direct numerical simulation (DNS). The idea is to mimic some experimental work techniques. In an experiment, by means of proper synchronization, one is able to arrive at the mean‐in‐the‐phase fields. In the present work, similar data are obtained through the use of numerical tools. The topology of the wake flow in terms of the movement of the coherent structures, fields of Reynolds stresses, position of saddles and centers, is obtained for certain values of the body's aspect ratio. A comparison with the circular cylinder experimental data is also provided. An important result is the finding of some specific stress distributions around the main structures. The stress level lines are a hint of where the energy transfer of the mean flow to the turbulence is more intense. Copyright © 2011 John Wiley & Sons, Ltd.     

Controlling the Turbulent Flow Past a Thick Airfoil by Means of Flow Enhancement in Vortex Cells Using Suction from Central Body Surfaces

On the basis of a numerical simulation of the turbulent steady-state flow past a thick airfoil with vortex cells built into the body contour, an unconventional technique for controlling flow separation by means of distributed suction from central bodies embedded in the cells is analyzed over a wide range of Reynolds numbers and suction velocities.     

A three‐dimensional vortex particle‐in‐cell method for vortex motions in the vicinity of a wall

A new vortex particle‐in‐cell method for the simulation of three‐dimensional unsteady incompressible viscous flow is presented. The projection of the vortex strengths onto the mesh is based on volume interpolation. The convection of vorticity is treated as a Lagrangian move operation but one where the velocity of each particle is interpolated from an Eulerian mesh solution of velocity–Poisson equations. The change in vorticity due to diffusion is also computed on the Eulerian mesh and projected back to the particles. Where diffusive fluxes cause vorticity to enter a cell not already containing any particles new particles are created. The surface vorticity and the cancellation of tangential velocity at the plate are related by the Neumann conditions. The basic framework for implementation of the procedure is also introduced where the solution update comprises a sequence of two fractional steps. The method is applied to a problem where an unsteady boundary layer develops under the impact of a vortex ring and comparison is made with the experimental and numerical literature. Copyright © 2001 John Wiley & Sons, Ltd.     

An attempt to improve accuracy of higher‐order statistics and spectra in direct numerical simulation of incompressible wall turbulence by using the compact schemes for viscous terms

We attempt to improve accuracy in the high‐wavenumber region in DNS of incompressible wall turbulence such as found in fully developed turbulent channel flow. In particular, it is shown that the improvement of accuracy of viscous terms in the Navier–Stokes equations leads to the improvement of accuracy of higher‐order statistics and various spectra. It is emphasized that increase in required computational cost will not be crucial when incompressible flow is simulated, because the introduction of a higher‐order scheme into the viscous terms does not increase computational cost for solving the Poisson equation. We introduced fourth‐order and eighth‐order central compact schemes for discretizing the viscous terms in DNS of a fully developed turbulent channel flow. The results are compared with those using second‐order and fourth‐order central‐difference schemes applied to the viscous terms and those obtained by the spectral method. The results show that accuracy improvement of the viscous terms improve accuracy of higher‐order statistics (i.e., skewness and flatness factors of streamwise velocity fluctuation) and various spectra of velocity and pressure fluctuations in the high‐wavenumber region. Copyright © 2013 John Wiley & Sons, Ltd.     

O(2)‐symmetry breaking bifurcation: with application to the flow past a sphere in a pipe

The steady, axisymmetric laminar flow of a Newtonian fluid past a centrally‐located sphere in a pipe first loses stability with increasing flow rate at a steady O(2)‐symmetry breaking bifurcation point. Using group theoretic results, a number of authors have suggested techniques for locating singularities in branches of solutions that are invariant with respect to the symmetries of an arbitrary group. These arguments are presented for the O(2)‐symmetry encountered here and their implementation for O(2)‐symmetric problems is discussed. In particular, how a bifurcation point may first be detected and then accurately located using an ‘extended system’ is described. Also shown is how to decide numerically if the bifurcating branch is subcritical or supercritical. The numerical solutions were obtained using the finite element code ENTWIFE. This has enabled the computation of the symmetry breaking bifurcation point for a range of sphere‐to‐pipe diameter ratios. A wire along the centerline of the pipe downstream of the sphere is also introduced, and its effect on the critical Reynolds number is shown to be small. Copyright © 2000 John Wiley & Sons, Ltd.     

Flow Visualization of Superbursts and of the Log-Layer in a DNS at $$ \operatorname{Re} _{\tau } = 950 $$

This paper uses direct numerical simulations (DNS) of turbulent flow in a channel at (Del álamo, Jiménez, Zandonade, Moser J Fluid Mech 500:135–144, 2004) to provide a picture of the turbulent structures making large contributions to the Reynolds shear stress. Considerable work of this type has been done for the viscous wall region at smaller , for which a log-layer does not exist. Recent PIV measurements of turbulent velocity fluctuations in a plane parallel to the direction of flow have emphasized the dominant contribution of large scale structures in the outer flow. This prompted Hanratty and Papavassiliou (The role of wall vortices in producing turbulence. In: Panton, R.L. (ed) Self-sustaining Mechanism of Wall Turbulence. Computational Mechanics Publications, Southampton, pp. 83–108, 1997) to use DNS at to examine these structures in a plane perpendicular to the direction of flow. They identified plumes which extend from the wall to the center of a channel. The data at are used to explore these results further, to examine the structure of the log-layer, and to test present notions about the viscous wall layer.     

On the distributed Lagrange multiplier/fictitious domain method for rigid‐particle‐laden flows: a proposition for an alternative formulation of the Lagrange multipliers

The distributed Lagrange multiplier/fictitious domain method proposed for the direct numerical simulation of particle‐laden flows is considered in this work. First, it is demonstrated that improved accuracy is obtained with a coupled numerical scheme, whereby the pressure and the Lagrange multiplier fields enforcing incompressibility and rigid body motion, respectively, are calculated and applied together. However, the convergence characteristics of the iterative solution of the coupled scheme are poor because symmetric but indefinite and poorly conditioned matrices are produced. An analysis is then presented, which suggests that the cause for the matrix pathologies lies in the interaction of the respective matrix operators enforcing incompressibility and rigid body motion. On the basis of this analysis, an alternative formulation is developed for the Lagrange multipliers, being now composed of a set of forces distributed only on the particle boundary together with a set of couples distributed within the particle core. The new formulation is tested with several types of flows with stationary or moving particles under creeping or finite Reynolds number conditions and it is demonstrated that it produces correct results and better conditioned matrices, thus enabling faster and more reliable convergence of the conjugate gradient method. The analysis and tests, therefore, support the expectation that the proposed formulation is promising and worthy of further study and improvement. Copyright © 2011 John Wiley & Sons, Ltd.     

Verification testing in computational fluid dynamics: an example using Reynolds‐averaged Navier–Stokes methods for two‐dimensional flow in the near wake of a circular cylinder

Verification testing was performed for various Reynolds‐averaged Navier–Stokes methods for uniform flow past a circular cylinder at Re= 5232. The standard and renormalized group (RNG) versions of the k–ε method were examined, along with the Boussinesq, Speziale and Launder constitutive relationships. Wind tunnel experiments for flow past a circular cylinder were also performed to obtain a comparative data set. Preliminary studies demonstrate poor convergence for the Speziale relationship. Verification testing with the standard and RNG k–ε models suggests that the simulations exhibit global monotonic convergence for the Boussinesq models. However, the global order of accuracy of the methods was much lower than the expected order of accuracy of 2. For this reason, pointwise convergence ratios and orders of accuracy were computed to show that not all sampling locations had converged (standard k–ε model: 19% failed to converge; RNG k–ε model: 14% failed to converge). When the non‐convergent points were removed from consideration, the average orders of accuracy are closer to the expected value (standard k–ε model: 1.41; RNG k–ε model: 1.27). Poor iterative and global grid convergence was found for the RNG k–ε/Launder model. The standard and RNG k–ε models with the Boussinesq relationship were compared with experimental data and yielded results significantly different from the experiments. Copyright © 2003 John Wiley & Sons, Ltd.     

Three‐dimensional transient Navier–Stokes solvers in cylindrical coordinate system based on a spectral collocation method using explicit treatment of the pressure

A spectral collocation method is developed for solving the three‐dimensional transient Navier–Stokes equations in cylindrical coordinate system. The Chebyshev–Fourier spectral collocation method is used for spatial approximation. A second‐order semi‐implicit scheme with explicit treatment of the pressure and implicit treatment of the viscous term is used for the time discretization. The pressure Poisson equation enforces the incompressibility constraint for the velocity field, and the pressure is solved through the pressure Poisson equation with a Neumann boundary condition. We demonstrate by numerical results that this scheme is stable under the standard Courant–Friedrichs–Lewy (CFL) condition, and is second‐order accurate in time for the velocity, pressure, and divergence. Further, we develop three accurate, stable, and efficient solvers based on this algorithm by selecting different collocation points in r‐, ? ‐, and z‐directions. Additionally, we compare two sets of collocation points used to avoid the axis, and the numerical results indicate that using the Chebyshev Gauss–Radau points in radial direction to avoid the axis is more practical for solving our problem, and its main advantage is to save the CPU time compared with using the Chebyshev Gauss–Lobatto points in radial direction to avoid the axis. Copyright © 2010 John Wiley & Sons, Ltd.