NUMERICAL SIMULATION OF UNSTEADY FLOW AROUND A CIRCULAR CYLINDER AT HIGH REYNOLDS NUMBERS

A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.

Determination of Strouhal Numbers at High Reynolds Numbers

Leeward of natural elevations, like mountains and hills, the air flow becomes turbulent and often times damaging and hazardous to aviation and downwind populations. There is currently a trend for massive construction projects, the result of which are megastructures that behave similarly to these natural elevations and create analogous turbulence conditions. Examples five mega projects were analyzed, and it was estimated that the Reynolds number variation in these buildings, is from 6.10g and 7.109, for winds between 10 m/s and 50 m/s. In this work, the authors present a first numerical approach to this phenomenon by calculating the Strouhal numbers induced by winds blowing against large-volume bodies, in the range of high Reynolds numbers. For this study, satellite images depicting von K^irm^n cloud streets leeward of isolated islands were used. The methodology employed was based on a satellite image where streets watch von K^rnfin vortex, from NOAA-ARL （National Oceanic and Atmospheric Administration-Air Resource Laboratory） prognosis was obtained for a grid point near the island, then determined the inversion layer and meteorological data （wind, temperature and pressure）, was measured from the satellite image the distances separating the vortices to calculate the period, the Reynolds number and Strouhal. The studied results of the cases are displayed graphically, where it is possible to observe a data dispersion as well as a rising trend of the Strouhal number as the Reynolds number increases.

Unsteady Numerical Simulation of Flow around 2-D Circular Cylinder for High Reynolds Numbers

In this paper, 2-D computational analyses were conducted for unsteady high Reynolds number flows around a smooth circular cylinder in the supercritical and upper-transition flow regimes, i.e. 8.21×104〈Re〈1.54×106. The calculations were performed by means of solving the 2-D Unsteady Reynolds-Averaged Navier-Stokes （URANS） equations with a k-ε turbulence model. The calculated results, produced flow structure drag and lift coefficients, as well as Strouhal numbers. The findings were in good agreement with previous published data, which also supplied us with a good understanding of the flow across cylinders of different high Reynolds numbers. Meanwhile, an effective measure was presented to control the lift force on a cylinder, which points the way to decrease the vortex induced vibration of marine structure in future.

Control of mean and fluctuating forces on a circular cylinder at high Reynolds numbers

A narrow strip is used to control mean and fluctuating forces on a circular cylinder at Reynolds numbers from 2.0 ×10^4 to 1.0 ×^ 10^5. The axes of the strip and cylinder are parallel. The control parameters are strip width ratio and strip position characterized by angle of attack and distance from the cylinder. Wind tunnel tests show that the vortex shedding from both sides of the cylinder can be suppressed, and mean drag and fluctuating lift on the cylinder can be reduced if the strip is installed in an effective zone downstream of the cylinder. A phenomenon of mono-side vortex shedding is found. The strip-induced local changes of velocity profiles in the near wake of the cylinder are measured, and the relation between base suction and peak value in the power spectrum of fluctuating lift is studied. The control mechanism is then discussed from different points of view.

Direct numerical simulation of turbulent flow over backward-facing at high Reynolds numbers

Direct numerical simulation(DNS) was performed for the first time to study the flow over a backward-facing step at a high Reynolds number on a coarse grid.The flow over backward-facing step is the typical turbulent flow controlled by large eddy,in which the effect of small eddy could be negligible as an approximation.The grid dimension could easily satisfy the resolution requirement to describe the characteristics of a large eddy flow.Therefore,direct numerical simulation of N-S equations to obtain the turbulent flow field on the coarse grid could be realized.Numerical simulation of a two-dimensional flow over a backward-facing step at a Reynolds number Re=44000 was conducted using Euler-Lagrange finite element scheme based on the efficient operator-splitting method(OSFEM).The flow field was descretized by triangle meshes with 16669 nodes.The overall computational time only took 150 min on a PC.Both the characteristics of time-averaged and instantaneous turbulent flow were simultaneously obtained.The analysis showed that the calculated results were in good agreement with the test data.Hence,the DNS approach could become the reality to solve the complex turbulent flow with high Reynolds numbers in practical engineering.

Wavelet solutions of Burgers' equation with high Reynolds numbers

A wavelet method is proposed to solve the Burgers’equation.Following this method,this nonlinear partial differential equation is first transformed into a system of ordinary differential equations using the modified wavelet Galerkin method recently developed by the authors.Then,the classical fourth-order explicit Runge–Kutta method is employed to solve the resulting system of ordinary differential equations.Such a wavelet-based solution procedure has been justified by solving two test examples:results demonstrate that the proposed method has a much better accuracy and efficiency than many other existing numerical methods,and whose order of convergence can go up to 5.Most importantly,our results also indicate that the present wavelet method can readily deal with those fluid dynamics problems with high Reynolds numbers.

Characteristics of settling of dilute suspension of particles with different density at high Reynolds numbers

Dilute suspension of particles with same density and size develops clusters when settle at high Reynolds number(≥250).It is due to particles entrapment in the wakes produced by upstream particles.In this work,this phenomenon is studied for suspension having particles with different densities by numerical simulations.The particle-fluid interactions are modelled using immersed boundary method and inter-particle collisions are modelled using discrete element method.In simulations,settling Reynolds number is always kept above 250 and the suspension solid volume fraction is nearly 0.1 percent.Two particle density ratios(i.e.density of heavy particles to lighter particles)equal to 4:1 and 2:1 and particles with same density are studied.For each density ratio,the percentage volume fraction of each particle density is nearly varied from 0.8 to 0.2.Settling characteristics such as microstructures of settling particle,average settling velocity and velocity fluctuations of settling particles are studied.Simulations show that for different density particles settling characteristics of suspension is largely dominated by heavy particles.At the end of paper,the underlying physics is explained for the anomalies observed in simulation.

Simulating high Reynolds number flow in two-dimensional lid-driven cavity by multi-relaxation-time lattice Boltzmann method

By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number （Re） flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.

The Simulation of High Reynolds Number Cavity Flow Based on Fractional Volumetric Lattice Boltzmann Method

The fractional volumetric lattice Boltzmann method with much better stability was used to simulate two-chmensional cavity flows. Because the effective viscosity was reduced by the fraction factor, it is very effective for simulating high Reynolds number flows. Simulations were carried out on a uniform grids system. The stream lines and the velocity profiles obtained from the simulations agree well with the standard lattice Boltzmann method simulations. Comparisons of detailed flow patterns with other studies via location of vortex centers are also satisfactory.

Study of poly-disperse aerosols deposition in turbulent flow with different Reynolds numbers

The turbulent deposition mechanism is one of the main mechanisms of aerosol deposition in nuclear power plant tubes.An experimental study of poly-disperse aerosol deposition in a horizontal tube is conducted,where the nominal Reynolds number(Re)is in a range of 3600–200,000.The aerosol deposition velocity first increases and then decreases with the increase of Res,and at high Re,particle rebound occurs during aerosol deposition in the tube.When the Re is low,the aerosol deposition velocity increases with the increase of aerosol diameter.When the Re is greater than 60,000,the deposition velocity first increases and then decreases with the increase of aerosol diameter due to particle surface rebound.A new aerosol deposition model has been developed by establishing the energy conservation equation of the rebounded particles in the viscous sublayer.The calculated results of the new model are in good agreement with these experimental results,and the error between the aerosol deposition velocity calculated by the model and experimental results is between−60%and 150%.

An artificial viscosity augmented physics-informed neural network for incompressible flow

Physics-informed neural networks(PINNs)are proved methods that are effective in solving some strongly nonlinear partial differential equations(PDEs),e.g.,Navier-Stokes equations,with a small amount of boundary or interior data.However,the feasibility of applying PINNs to the flow at moderate or high Reynolds numbers has rarely been reported.The present paper proposes an artificial viscosity(AV)-based PINN for solving the forward and inverse flow problems.Specifically,the AV used in PINNs is inspired by the entropy viscosity method developed in conventional computational fluid dynamics(CFD)to stabilize the simulation of flow at high Reynolds numbers.The newly developed PINN is used to solve the forward problem of the two-dimensional steady cavity flow at Re=1000 and the inverse problem derived from two-dimensional film boiling.The results show that the AV augmented PINN can solve both problems with good accuracy and substantially reduce the inference errors in the forward problem.

NUMERICAL STUDY OF PARTICLE DISTRIBUTION IN THE WAKE OF GAS-PARTICLE TWO-PHASE FLOWS PAST A CIRCULAR CYLINDER AT HIGH REYNOLDS NUMBER

Particle-laden gas flows past a circular cylinder at the Reynolds number of 2×10^(5) were numerically investigated. The Discrete Vortex Method (DVM) was employed to evaluate the unsteady gas flow fields and a Lagrangian approach was applied for tracking individual solid particles. The vortex patterns and the distributions of particles with different Stokes numbers were obtained. Numerical results show that: (1) at small Stokes number (St=0.01) the particles move with the fluid and could be found evenly throughout the flow, (2) the regions around the vortex cores, where few particles exist, become wider as the stokes number of particles increases from 0.01 to 1.0, (3) at middle Stokes number (St=1.0, 10) centrifugal forces throw the particles out of the wake vortices, (4) at high Stokes number (St=100, 1000) the particles are not affected by the vortices,and their motion is determined by their inertia effects.

Study on the dynamic characteristics of flow over building cluster at high Reynolds number by large eddy simulation

In this paper,the dynamic characteristics of building clusters are simulated by large eddy simulation at high Reynolds number for both homogeneous and heterogeneous building clusters.To save the computational cost a channel-like flow model is applied to the urban canopy with free slip condition at the upper boundary.The results show that the domain height is an important parameter for correct evaluation of the dynamic characteristics.The domain height must be greater than 8h(h is the average building height)in order to obtain correct roughness height while displacement height and roughness sublayer are less sensitive to the domain height.The Reynolds number effects on the dynamic characteristics and flow patterns are investigated.The turbulence intensity is stronger inside building cluster at high Reynolds number while turbulence intensity is almost unchanged with Reynolds number above the building cluster.Roughness height increases monotonously with Reynolds number by 20%from Re*=103 to Re*=105 but displacement height is almost unchanged.Within the canopy layer of heterogeneous building clusters,flow structures vary between buildings and turbulence is more active at high Reynolds number.

NUMERICAL STUDY OF PARTICLE DISTRIBUTION IN WAKE OF LIQUID-PARTICLE FLOWS PAST A CIRCULAR CYLINDER USING DISCRETE VORTEX METHOD

Particle-laden water flows past a circular cylinder were numerically investigated. The discrete vortex method （DVM） was employed to evaluate the unsteady water flow fields and a Lagrangian approach was applied for tracking individual solid particles. A dispersion function was defined to represent the dispersion scale of the particle. The wake vortex patterns, the distributions and the time series of dispersion functions of particles with different Stokes numbers were obtained. Numerical results show that the particle distribution in the wake of the circular cylinder is closely related to the particle＇s Stokes number and the structure of wake vortices： （1） the intermediate sized particles with Stokes numbers, St, of 0.25, 1.0 and 4.0 can not enter the vortex cores and concen- trate near the peripheries of the vortex structures, （2） in the circular cylinder wake, the dispersion intensity of particles decreases as St is increased from 0.25 to 4.0.

SUPG finite element method based on penalty function for lid-driven cavity flow up to Re = 27500

A streamline upwind/Petrov-Galerkin （SUPG） finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.

Finite difference streamline diffusion method using nonconforming space for incompressible time-dependent Navier-Stokes equations

This paper proposes a new nonconforming finite difference streamline diffusion method to solve incompressible time-dependent Navier-Stokes equations with a high Reynolds number. The backwards difference in time and the Crouzeix-Raviart （CR） element combined with the P0 element in space are used. The result shows that this scheme has good stabilities and error estimates independent of the viscosity coefficient.

A NEW STABILIZED FINITE ELEMENT METHOD FOR SOLVING TRANSIENT NAVIER-STOKES EQUATIONS WITH HIGH REYNOLDS NUMBER

In this paper,we present a new stabilized finite element method for transient Navier-Stokes equations with high Reynolds number based on the projection of the velocity and pressure.We use Taylor-Hood elements and the equal order elements in space and second order difference in time to get the fully discrete scheme.The scheme is proven to possess the absolute stability and the optimal error estimates.Numerical experiments show that our method is effective for transient Navier-Stokes equations with high Reynolds number and the results are in good agreement with the value of subgrid-scale eddy viscosity methods,Pet ro-Galerkin finite element method and st reamline diffusion method.

ACCURATE COMPUTATION AND INTERPOLATION TECHNIQUE OF FINITE ANALYTIC COEFFICIENTS

The finite analytic method (FA) developed in the last decade is an effective numerical method for solving fluid flow problems. However, because of the limitation in the present computer, large round-off errors are found in calculating FA coefficients when Reynolds number is large. This paper investigates the cause of this difficulty and presents a special programming technique in making an accurate computation of FA coefficients. Then a fundamental function known as 'Pe' is tabulated by the accurate computation. In practical application the interpolation technique is employed so that the FA coefficients can be obtained reliably and quickly.

Numerical Research of Flow past a Circular Cylinder with Splitter Plate at a Subcritical Reynolds Number Region

Numerical research of flow past a circular cylinder with a splitter at the subcritical Reynolds number region of 5 × 10~4—9 × 10~4 was researched based on Computational Fluid Dynamics(CFD) through solving twodimensional incompressible unsteady Reynolds-averaged Navier-Stokes(URANS) equations with the shear stress transport(SST) k-ω turbulence model. Three different grid resolutions were employed in the verification and validation study of the adopted turbulence model. Various fluid characteristics such as Strouhal number, lift coefficient of the cylinder and the splitter with respect to various splitter lengths and different Reynolds numbers were investigated. It was revealed that the lift coefficient ratio of the splitter over the cylinder remains near 1.6 when the splitter length is 1.5—4 times the cylinder's diameter. Vortex shedding is strongly inhibited when the splitter length is greater than a critical value of around four times the cylinder's diameter. The phase difference of the lift coefficient on the upper and lower surface of the splitter varies between-30?and 30?. The maximal lift coefficients are reached when the splitter length is about 2 times the cylinder's diameter. Besides, the splitter length has little influence on the separation angle around the cylinder.

A Finite Element Variational Multiscale Method for Stationary Incompressible Magnetohydrodynamics Equations

In this paper,we propose a variational multiscale method(VMM)for the stationary incompressible magnetohydrodynamics equations.This method is defined by large-scale spaces for the velocity field and the magnetic field,which aims to solve flows at high Reynolds numbers.We provide a new VMM formulation and prove its stability and convergence.Finally,some numerical experiments are presented to indicate the optimal convergence of our method.