首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到10条相似文献,搜索用时 156 毫秒
1.
The pulsating cross-flow over a single circular cylinder at the subcritical Reynolds number ReD = 2580 is studied with the large eddy simulation (LES) technique using the standard Smagorinsky model as well as a dynamic model in which the test filtered quantities are evaluated through a truncated Taylor series expansion. The filtered equations are discretised using the finite volume method in an unstructured, collocated grid arrangement with a second-order accurate method, in both space and time. The predictions are compared against very detailed experiments for mean velocities and Reynolds stresses that were performed in a duct of cross-section 72 mm × 72 mm using the PIV technique. The effects of mesh refinement close to the cylinder as well as of subgrid scale model are also examined. The numerical predictions are in very good agreement with the measurements in terms of mean as well as turbulence quantities. The instantaneous flow patterns of the flow field are examined and the effect of the external flow pulsation on the wake characteristics such as vortex formation length, vortex strength, Strouhal number as well as the lift and drag coefficients is quantified. The vortex formation length is decreased while the mean drag, as well as the rms values of the drag and lift coefficients increase significantly under pulsating flow conditions. The performance of the LES technique is analysed in the light of the wake characteristics.  相似文献   

2.
An efficient solver integrating the restarted simpler generalized minimal residual method (SGMRES(m)) with finite volume method (FVM) on triangular grid is developed to simulate the viscoelastic fluid flows. In particular, the SGMRES(m) solver is used to solve the large-scale sparse linear systems, which arise from the course of FVM on triangular grid for modeling the Newtonian and the viscoelastic fluid flows. To examine the performance of the solver for the nonlinear flow equations of viscoelastic fluids, we consider two types of numerical tests: the Newtonian flow past a circular cylinder, and the Oldroyd-B fluid flow in a planar channel and past a circular cylinder. It is shown that the numerical results obtained by the SGMRES(m) are consistent with the analytical solutions or empirical values. By comparing CPU time of different solvers, we find our solver is a highly efficient one for solving the flow equations of viscoelastic fluids.  相似文献   

3.
Variational multiscale large-eddy simulations (VMS–LES) of the flow around a circular cylinder are carried out at different Reynolds numbers in the subcritical regime, viz. Re = 3900, 10,000 and 20,000, based on the cylinder diameter. A mixed finite-element/finite-volume discretization on unstructured grids is used. The separation between the largest and the smallest resolved scales is obtained through a variational projection operator and finite-volume cell agglomeration. The WALE subgrid scale model is used to account for the effects of the unresolved scales; in the VMS approach, it is only added to the smallest resolved ones. The capability of this methodology to accurately predict the aerodynamic forces acting on the cylinder and in capturing the flow features are evaluated for the different Reynolds numbers considered. The sensitivity of the results to different simulation parameters, viz. agglomeration level and numerical viscosity, is also investigated at Re = 20,000.  相似文献   

4.
Two-dimensional fluid flow around an oscillating circular cylinder is studied numerically at different values of oscillation frequency and amplitude. A novel finite element method which uses discretization along the characteristic line is used for simulation. The solver is coupled to a mesh movement scheme using the Arbitrary Lagrangian-Eulerian (ALE) formulation to account for body motion in the flow field. Two cases of cylinder motion have been studied, cross flow and inline oscillation. In both cases, occurrence of lock on is investigated and the bounds of the lock on region are determined. A comparison of the numerical results with the experimental data indicates that 2D simulation is valid up to Re = 300. Beyond that, 3D effects appear. By using flow visualization, effect of a cylinder oscillation on the flow field and wake pattern has been studied. Also, variation of the mean drag coefficient against the oscillation parameters is discussed. The numerical results are in good agreement with the experimental data available in the literature.  相似文献   

5.
A three-dimensional parallel unstructured non-nested multigrid solver for solutions of unsteady incompressible viscous flow is developed and validated. The finite-volume Navier–Stokes solver is based on the artificial compressibility approach with a high-resolution method of characteristics-based scheme for handling convection terms. The unsteady flow is calculated with a matrix-free implicit dual time stepping scheme. The parallelization of the multigrid solver is achieved by multigrid domain decomposition approach (MG-DD), using single program multiple data (SPMD) and multiple instruction multiple data (MIMD) programming paradigm. There are two parallelization strategies proposed in this work, first strategy is a one-level parallelization strategy using geometric domain decomposition technique alone, second strategy is a two-level parallelization strategy that consists of a hybrid of both geometric domain decomposition and data decomposition techniques. Message-passing interface (MPI) and OpenMP standard are used to communicate data between processors and decompose loop iterations arrays, respectively. The parallel-multigrid code is used to simulate both steady and unsteady incompressible viscous flows over a circular cylinder and a lid-driven cavity flow. A maximum speedup of 22.5 could be achieved on 32 processors, for instance, the lid-driven cavity flow of Re = 1000. The results obtained agree well with numerical solutions obtained by other researchers as well as experimental measurements. A detailed study of the time step size and number of pseudo-sub-iterations per time step required for simulating unsteady flow are presented in this paper.  相似文献   

6.
We present a validation strategy for enhancement of an unstructured industrial finite-volume solver designed for steady RANS problems for large-eddy-type simulation with near-wall modelling of incompressible high Reynolds number flow. Different parts of the projection-based discretisation are investigated to ensure LES capability of the numerical method. Turbulence model parameters are calibrated by using a minimisation of least-squares functionals for first and second order statistics of the basic benchmark problems decaying homogeneous turbulence and turbulent channel flow. Then the method is applied to the flow over a backward facing step at Reh = 37,500. Of special interest is the role of the spatial and temporal discretisation error for low order schemes. For wall-bounded flows, present results confirm existing best practice guidelines for mesh design. For free-shear layers, a sensor to quantify the resolution quality of the LES based on the resolved turbulent kinetic energy is presented and applied to the flow over a backward facing step at Reh = 37,500.  相似文献   

7.
An efficient spectral vanishing viscosity method for the large-eddy simulation of incompressible flows is proposed, both for standard spectral and spectral element approximations. The approach is integrated in a collocation spectral Chebyshev-Fourier solver and then used to compute the turbulent wake of a cylinder in a crossflow confined geometry (Reynolds number Re=3900)  相似文献   

8.
Incompressible turbulent flow over a backward facing step at Reh=5100 is investigated by large eddy simulations (LES). The ratio of the oncoming boundary layer thickness δ to the step height h was set to 1.2. Additionally channel flows at various Reτ numbers are presented for the validation of the numerical code. The results are compared with existing DNS and experimental databases. The present study focuses on different procedures for LES of engineering problems in complex geometries using structured rectangular grids. Two different methods that are able to treat complex geometrical configurations are implemented, examined and compared; namely the domain decomposition approach based on Schur’s complement and the immersed boundary method. In the present study both methods make use of a fast direct Poisson’s pressure solver based on a heavily modified version of the public domain package FISHPAK. The latter was optimised and fully parallelised for shared memory architectures, for solutions on rectangular grids stretched in one or two directions. The resulting code reaches performances of 1.0 μs/node/iter, allowing low cost computations on grids of the order of million points. The main objective of the present study was to investigate the potential of different methods for LES in complex geometrical configurations like bluff body flows and wakes. One of the main findings is that careful selection of numerical methods and implementation techniques can lead to accurate and very efficient codes, where the geometric complexity does not lead to algorithmic or numerical complexity.  相似文献   

9.
Neural networks based subgrid scale modeling in large eddy simulations   总被引:1,自引:0,他引:1  
In this paper a multilayer feed-forward neural network (NN) is used as subgrid scale (SGS) model in a large eddy simulation (LES). The NN was previously off-line trained using numerical data generated by a LES of a channel flow at Reτ=180 with Bardina's scale similar (BFR) SGS model. Results show the ability of NNs to identify and reproduce the highly nonlinear behavior of the turbulent flows, and therefore the possibility of using NN techniques in numerical simulations of turbulent flows.  相似文献   

10.
The effects of numerical viscosity, subgrid scale (SGS) viscosity and grid resolution are investigated in LES and VMS-LES simulations of the flow around a circular cylinder at Re=3900 on unstructured grids. The separation between the largest and the smallest resolved scales in the VMS formulation is obtained through a variational projection operator and finite-volume cell agglomeration. Three different non-dynamic eddy-viscosity SGS models are used both in classical and in VMS-LES. The so-called small-small formulation is used in VMS-LES, i.e. the SGS viscosity is computed as a function of the smallest resolved scales. Two different grid resolutions are considered. It is found that, for each considered SGS model, the amount of SGS viscosity introduced in the VMS-LES formulation is significantly lower than in classical LES. This, together with the fact that in the VMS formulation the SGS viscosity only acts on the smallest resolved scales, has a strong impact on the results. However, a significant sensitivity of the results to the considered SGS model remains also in the VMS-LES formulation. Moreover, passing from classical LES to VMS-LES does not systematically lead to an improvement of the quality of the numerical predictions.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号