Pipe Flow and Wall Turbulence Using a Modified Navier-Stokes Equation

We use a derived incompressible modified Navier-Stokes equation to model pipe flow and wall turbulence.We reproduce the observed flattened paraboloid velocity profiles of turbulence that cannot be obtained directly using standard incompressible Navier-Stokes equation.The solutions found are in harmony with multi-valued velocity fields as a definition of turbulence.Repeating the procedure for the flow of turbulent fluid between two parallel flat plates we find similar flattened velocity profiles.We extend the analysis to the turbulent flow along a single wall and compare the results with experimental data and the established controversial von Karman logarithmic law of the wall.     

Lattice Boltzmann model for three-dimensional decaying homogeneous isotropic turbulence

We implement a lattice Boltzmann method (LBM) for decaying homogeneous isotropic turbulence based on an analogous Galerkin filter and focus on the fundamental statistical isotropic property. This regularized method is constructed based on orthogonal Hermite polynomial space. For decaying homogeneous isotropic turbulence, this regularized method can simulate the isotropic property very well. Numerical studies demonstrate that the novel regularized LBM is a promising approximation of turbulent fluid flows, which paves the way for coupling various turbulent models with LBM.     

New Sedov-Type Solution of Isotropic Turbulence

The starting point lies in the results obtained by Sedov （1944） for isotropic turbulence with a self-preserving hypothesis. A careful consideration of the mathematical structure of the Karman-Howarth equation leads to an exact analysis of all cases possible and to all admissible solutions of the problem. I study this interesting problem from a new point of view. New solutions are obtained. Based on these exact solutions, some physical significant consequences of recent advances in the theory of self-preserved homogeneous statistical solution of the Navier-Stokes equations are presented.     

Curvature Effect on Compressible Turbulent Flow over a Wavy Wall

A fully developed compressible turbulent flow in a channel with a lower wavy wall and a upper plane wall is studied using large eddy simulation. We mainly attempt to deal with the curvature effect on compressible turbulent flow over the wavy wall. Some typical quantities including the mean turbulence statistics, dilatation and baroclinic terms in the enstrophy equation, turbulent kinetic energy budgets and the near-wall turbulent structures are analysed. The results obtained in this study provide physical insight into the understanding of the effects of curvature and compressibility on wall-bounded compressible turbulent flow.     

Controlling Flow Turbulence Using Local Pinning Feedback

Flow turbulence control in two-dimensional Navier-Stokes equation is considered. By applying local pinning control only to a single component of flow velocity field, the flow turbulence can be controlled to desirable targets. It is found that with certain number of controllers there exist an optimal control strength at which control error takes minimum value, and larger and smaller control strengths give worse control efficiency. The physical mechanism underlying these strange control results is analysed based on the interactions between different types of modes.     

Nonperturbative renormalization group approach to turbulence

We suggest a new, renormalization group (RG) based, nonperturbative method for treating the intermittency problem of fully developed turbulence which also includes the effects of a finite boundary of the turbulent flow. The key idea is not to try to construct an elimination procedure based on some assumed statistical distribution, but to make an ansatz for possible RG transformations and to pose constraints upon those, which guarantee the invariance of the nonlinear term in the Navier-Stokes equation, the invariance of the energy dissipation, and other basic properties of the velocity field. The role of length scales is taken to be inverse to that in the theory of critical phenomena; thus possible intermittency corrections are connected with the outer length scale. Depending on the specific type of flow, we find different sets of admissible transformations with distinct scaling behaviour: for the often considered infinite, isotropic, and homogeneous system K41 scaling is enforced, but for the more realistic plane Couette geometry no restrictions on intermittency exponents were obtained so far. Received: 28 December 1997 / Accepted: 6 August 1998     

Quasi-stationary states in a circular geometry

The quasi-stationary final states of decaying two-dimensional turbulence on a circular domain, with an initial flow field containing either no or a substantial amount of angular momentum, have been investigated numerically. The production of angular momentum is almost absent for these flows on a circular domain with a no-slip wall. Its presence or absence essentially determines the character of the quasi-stationary final state. Based on a minimum-enstrophy principle a diagram is constructed that provides insight into the development of the typical late-time flow patterns on a circular domain with a no-slip wall. The quasi-stationary final states found in the present numerical study can be understood based on the predictions from the minimum-enstrophy principle.     

Three-dimensional lattice Boltzmann method for simulating blood flow in aortic arch

The three-dimensional （3D） lattice Boltzmann models, 3DQ15, 3DQ19 and 3DQ27, under different wall boundary conditions and lattice resolutions have been investigated by simulating Poiseuille flow in a circular cylinder for a wide range of Reynolds numbers. The 3DQ19 model with improved Fillippova and Hanel （FH） curved boundary condition represents a good compromise between computational efficiency and reliability. Blood flow in an aortic arch is then simulated as a typical haemodynamic application. Axial and secondary fluid velocity and effective wall shear stress profiles in a 180° bend are obtained, and the results also demonstrate that the lattice Boltzmann method is suitable for simulating the flow in 3D large-curved vessels.     

Development of Lattice Boltzmann Flux Solver for Simulation of Incompressible Flows

A lattice Boltzmann flux solver (LBFS) is presented in this work for simulation of incompressible viscous and inviscid flows. The new solver is based on Chapman-Enskog expansion analysis, which is the bridge to link Navier-Stokes (N-S) equations and lattice Boltzmann equation (LBE). The macroscopic differential equations are discretized by the finite volume method, where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers. The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. LBFS is validated by its application to simulate the viscous decaying vortex flow, the driven cavity flow, the viscous flow past a circular cylinder, and the inviscid flow past a circular cylinder. The obtained numerical results compare very well with available data in the literature, which show that LBFS has the second order of accuracy in space, and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.     

Upper bounds on the convective heat transport in a rotating fluid layer of infinite Prandtl number: case of large Taylor numbers

By means of the Howard-Busse method of the optimum theory of turbulence we obtain upper bounds on the convective heat transport in a horizontal fluid layer heated from below and rotating about a vertical axis. We consider the interval of large Taylor numbers where the intermediate layers of the optimum fields expand in the direction of the corresponding internal layers. We consider the 1 - α-solution of the arising variational problem for the cases of rigid-stress-free, stress-free, and rigid boundary conditions. For each kind of boundary condition we discuss four cases: two cases where the boundary layers are thinner than the Ekman layers of the optimum field and two cases where the boundary layers are thicker than the Ekman layers. In most cases we use an improved solution of the Euler-Lagrange equations of the variational problem for the intermediate layers of the optimum fields. This solution leads to corrections of the thicknesses of the boundary layers of the optimum fields and to lower upper bounds on the convective heat transport in comparison to the bounds obtained by Chan [J. Fluid Mech. 64, 477 (1974)] and Hunter and Riahi [J. Fluid Mech. 72, 433 (1975)]. Compared to the existing experimental data for the case of a fluid layer with rigid boundaries the corresponding upper bounds on the convective heat transport is less than two times larger than the experimental results, the corresponding upper bound on the convective heat transport, obtained by Hunter and Riahi is about 10% higher than the bound obtained in this article. When Rayleigh number and Taylor number are high enough the upper bound on the convective heat transport ceases to depend on the boundary conditions. Received 30 January 2001 and Received in final form 28 May 2001     

Minimum description of the onset of pipe turbulence

To demonstrate essentials of the mechanism for the onset of turbulence in a pipe at Re=2000, 48 degrees of freedom are enough. The derivation from the Navier-Stokes equation uses a novel type of modes which guarantee linear stability. For the reduction of the nonlinear interactions, the modes are grouped in 3 blocks. Facilitated by these simplifications the interdependence between linear and nonlinear processes is analysed, however, just for a special example. A phenomenon resembling backflow is identified. Received 20 October 1997     

Laser Doppler velocimetry analysis of transitional pipe flow

The objective of this study has been to experimentally analyze the correlation structure of the strong temporal intermittency which characterizes pipe flow close to the transition to turbulence. In doing so transitional pipe flow has been analyzed by Laser Doppler velocimetry and the Reynolds number dependence of the covariance function has been studied. The range which has been analyzed covers the transition to turbulence and moderately developed turbulence (Reynolds number from 1 500 to 5 000). The correlation structure which has been evidenced is generally in agreement with the deterministic, dynamical, interpretation of temporal intermittency which explains the intermittent behavior as a result of a saddle node bifurcation. However, the analysis has evidenced fluctuations even before the onset of turbulence. The structure of these fluctuations is perfectly autoregressive which leads us to conclude that the transition to turbulence can be viewed as a transition from linear randomness to (non-linear) homogeneity. Received 29 March 1999 and Received in final form 6 September 1999     

Two-dimensional, self-consistent, three-moment simulation of RFglow discharge

A two-dimensional self-consistent nonequilibrium fluid model is used to simulate radio frequency (RF) glow discharges to evaluate the quantitative effects of the radial and axial flow dynamics inside a cylindrically symmetric parallel-plate geometry. This model is based on the three moments of the Boltzmann equation and on Poisson's equation. Radial/axial flow dynamics of plasma in low-pressure parallel plate RF glow discharges are investigated. Instead of uniform profiles along the radial direction, which are assumed in one-dimensional models, nonplate profiles are obtained from the two-dimensional simulations. Ionization rate and three moment distributions of plasma density, average velocity, and mean energy are presented in a two-dimensional configuration. The maximum ionization rate occurs in the radial sheath region and agrees with experimental results. Variations in ion density distributions at different positions, various gas pressures frequencies, and applied fields are discussed     

Immersed boundary method for the simulation of 2D viscous flow based on vorticity–velocity formulations

A new immersed boundary method based on vorticity–velocity formulations for the simulation of 2D incompressible viscous flow is proposed in present paper. The velocity and vorticity are respectively divided into two parts: one is the velocity and vorticity without the influence of the immersed boundary, and the other is the corrected velocity and the corrected vorticity derived from the influence of the immersed boundary. The corrected velocity is obtained from the multi-direct forcing to ensure the well satisfaction of the no-slip boundary condition at the immersed boundary. The corrected vorticity is derived from the vorticity transport equation. The third-order Runge–Kutta for time stepping, the fourth-order finite difference scheme for spatial derivatives and the fourth-order discretized Poisson for solving velocity are applied in present flow solver. Three cases including decaying vortices, flow past a stationary circular cylinder and an in-line oscillating cylinder in a fluid at rest are conducted to validate the method proposed in this paper. And the results of the simulations show good agreements with previous numerical and experimental results. This indicates the validity and the accuracy of present immersed boundary method based on vorticity–velocity formulations.     

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier–Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.     

Nonlinear transport in the Boltzmann limit

Formal expressions for the irreversible fluxes of a simple fluid are obtained as functionals of the thermodynamic forces and local equilibrium time correlation functions. The Boltzmann limit of the correlation functions is shown to yield expressions for the irreversible fluxes equivalent to those obtained from the nonlinear Boltzmann kinetic equation. Specifically, for states near equilibrium, the fluxes may be formally expanded in powers of the thermodynamic gradients and the associated transport coefficients identified as integrals of time correlation functions. It is proved explicitly through nonlinear Burnett order that the time correlation function expressions for these transport coefficients agree with those of the Chapman-Enskog expansion of the nonlinear Boltzmann equation. For states far from equilibrium the local equilibrium time correlation functions are determined in the Boltzmann limit and a similar equivalence to the Boltzmann equation solution is established. Other formal representations of the fluxes are indicated; in particular, a projection operator form and its Boltzmann limit are discussed. As an example, the nonequilibrium correlation functions for steady shear flow are calculated exactly in the Boltzmann limit for Maxwell molecules.Research supported in part by NSF grant PHY 76-21453.     

Turbulence in diffusion replicator equation

Dynamical behaviors in the diffusion replicator equation of three species are numerically studied. We point out the significant role of the heteroclinic cycle in the equation, and analyze the details of the turbulent solution that appeared in this system. Firstly, the bifurcation diagram for a certain parameter setting is drawn. Then it is shown that the turbulence appears with the supercritical Hopf bifurcation of a stationary uniform solution and it disappears under a subcritical-type bifurcation. Secondly, the statistical property of the turbulence near the supercritical Hopf onset point is analyzed precisely. Further, the correlation lengths and correlation times obey some characteristic scaling laws.     

Rotation effect on near-wall turbulence statistics and flow structures

Turbulent flow in an axially rotating pipe, involving complicated physical mecha- nism of turbulence, is a typical problem for the study of rotating turbulent flow. The pipe rotation induces two effects on the flow. One is the stabilizing effect due to the centrifu- gal and Coriolis forces, which accounts for the relaminarization of the turbulence[1—3] and the reduction of the friction coefficient at the pipe wall. The behavior is also related to the wall streaks inclining to the azimuthal di…     

Large Eddy Simulation of Turbulent Channel Flow with 3D Roughness Using a Roughness Element Model

Large eddy simulation of turbulent channel flow with dense and small 3D roughness elements is carried out using a roughness element model Profiles of mean Reynolds stress, mean velocity and rms velocity as well as turbulent structures near the wall are obtained. The shear stress in the rough wall is larger than that in the smooth wall side and the rough wall has a larger influence on the channel flow. Profiles of mean streamwise velocity near the wall have logarithmic velocity distributions for both smooth and roughness walls, while there is a velocity decrease for the rough wall due to larger fractional drag. All the three components of rms velocities in the rough wall region are larger than that in the smooth wall region, and the roughness elements on the wall increase turbulent intensity in all directions. The s~reak spacing and average diameter of near wall quasi-s~reamwise vortices increase with the presence of roughness elements on the wall and it is shown that the rough wall induces complex and strong streamwise vortices. Results of dense and small 3D roughness elements in both turbulent statistics and structure, obtained with a relatively simple method, are found to be comparable to related experiments.