A high‐order element based adaptive mesh refinement strategy for three‐dimensional unstructured grid

Adaptive mesh refinement (AMR) shows attractive properties in automatically refining the flow region of interest, and with AMR, better prediction can be obtained with much less labor work and cost compared to manually remeshing or the global mesh refinement. Cartesian AMR is well established; however, AMR on hybrid unstructured mesh, which is heavily used in the high‐Reynolds number flow simulation, is less matured and existing methods may result in degraded mesh quality, which mostly happens in the boundary layer or near the sharp geometric features. User intervention or additional constraints, such as freezing all boundary layer elements or refining the whole boundary layer, are required to assist the refinement process. In this work, a novel AMR strategy is developed to handle existing difficulties. In the new method, high‐order unstructured elements are first generated based on the baseline mesh; then the refinement is conducted in the parametric space; at last, the mesh suitable for the solver is output. Generating refined elements in the parametric space with high‐order elements is the key of this method and this helps to guarantee both the accuracy and robustness. With the current method, 3‐dimensional hybrid unstructured mesh of huge size and complex geometry can be automatically refined, without user intervention nor additional constraints. With test cases including the 2‐dimensional airfoil and 3‐dimensional full aircraft, the current AMR method proves to be accurate, simple, and robust.     

Accurate and robust adaptive mesh refinement for aerodynamic simulation with multi‐block structured curvilinear mesh

A multi‐block curvilinear mesh‐based adaptive mesh refinement (AMR) method is developed to satisfy the competing objectives of improving accuracy and reducing cost. Body‐fitted curvilinear mesh‐based AMR is used to capture flow details of various length scales. A series of efforts are made to guarantee the accuracy and robustness of the AMR system. A physics‐based refinement function is proposed, which is proved to be able to detect both shock wave and vortical flow. The curvilinear mesh is refined with cubic interpolation, which guarantees the aspect ratio and smoothness. Furthermore, to enable its application in complex configurations, a sub‐block‐based refinement strategy is developed to avoid generating invalid mesh, which is the consequence of non‐smooth mesh lines or singular geometry features. A newfound problem of smaller wall distance, which negatively affects the stability and is never reported in the literature, is also discussed in detail, and an improved strategy is proposed. Together with the high‐accuracy numerical scheme, a multi‐block curvilinear mesh‐based AMR system is developed. With a series of test cases, the current method is verified to be accurate and robust and be able to automatically capture the flow details at great cost saving compared with the global refinement. Copyright © 2015 John Wiley & Sons, Ltd.     

Efficient computation of mean drag for the subcritical flow past a circular cylinder using general Galerkin G2

General Galerkin (G2) is a new computational method for turbulent flow, where a stabilized Galerkin finite element method is used to compute approximate weak solutions to the Navier–Stokes equations directly, without any filtering of the equations as in a standard approach to turbulence simulation, such as large eddy simulation, and thus no Reynolds stresses are introduced, which need modelling. In this paper, G2 is used to compute the drag coefficient c_D for the flow past a circular cylinder at Reynolds number Re=3900, for which the flow is turbulent. It is found that it is possible to approximate c_D to an accuracy of a few percent, corresponding to the accuracy in experimental results for this problem, using less than 10⁵ mesh points, which makes the simulations possible using a standard PC. The mesh is adaptively refined until a stopping criterion is reached with respect to the error in a chosen output of interest, which in this paper is c_D. Both the stopping criterion and the mesh‐refinement strategy are based on a posteriori error estimates, in the form of a space–time integral of residuals times derivatives of the solution of a dual problem, linearized at the approximate solution, and with data coupling to the output of interest. Copyright © 2008 John Wiley & Sons, Ltd.     

Mesh adaptation framework for embedded boundary methods for computational fluid dynamics and fluid-structure interaction

Embedded Boundary Methods (EBMs) are often preferred for the solution of Fluid-Structure Interaction (FSI) problems because they are reliable for large structural motions/deformations and topological changes. For viscous flow problems, however, they do not track the boundary layers that form around embedded obstacles and therefore do not maintain them resolved. Hence, an Adaptive Mesh Refinement (AMR) framework for EBMs is proposed in this paper. It is based on computing the distance from an edge of the embedding computational fluid dynamics mesh to the nearest embedded discrete surface and on satisfying the y⁺ requirements. It is also equipped with a Hessian-based criterion for resolving flow features such as shocks, vortices, and wakes and with load balancing for achieving parallel efficiency. It performs mesh refinement using a parallel version of the newest vertex bisection method to maintain mesh conformity. Hence, while it is sufficiently comprehensive to support many discretization methods, it is particularly attractive for vertex-centered finite volume schemes where dual cells tend to complicate the mesh adaptation process. Using the EBM known as FIVER, this AMR framework is verified for several academic FSI problems. Its potential for realistic FSI applications is also demonstrated with the simulation of a challenging supersonic parachute inflation dynamics problem.     

ADAPTIVE GRID REFINEMENT USING CELL-LEVEL AND GLOBAL IMBALANCES

A methodology for local solution-adaptive mesh refinement in computational fluid dynamics (CFD) using cell-level and global kinetic energy balances is formulated and tested. Results are presented for two two-dimensional steady incompressible laminar benchmark problems: a lid-driven cavity (Reynolds number Re=1000) and a backward-facing step (Re=400). It is demonstrated that local kinetic energy imbalance correlates with local solution accuracy, that normalized global imbalance is an appropriate criterion for halting mesh refinement and that a specified level of accuracy is realized at lower computational effort using local refinement compared with a uniform finer mesh. © 1997 by John Wiley and Sons, Ltd.     

Mesh adaptation for large‐eddy simulations in complex geometries

Large‐eddy simulation (LES) consists in explicitly simulating the large scales of the fluid motion and in modeling the influence of the smallest scales. Thanks to the steady growth of computational resources, LES can now be used to simulate realistic systems with complex geometries. However, when LES is used in such complex geometries, an adequate mesh has to be determined to perform valid LES. In this work, a strategy is proposed to assess the quality of a given mesh and to adapt it locally. Two different criteria are used as mesh adaptation criteria. The first criterion is defined to ensure a correct discretization of the mean field, whereas the second criterion is defined to ensure enough explicit resolution of turbulent scales motions. The use of both criteria is shown in canonical flow cases. As a second part of this work, a numerical strategy for mesh adaptation in high‐performance computing context is proposed by coupling the flow solver, YALES2, and the remeshing library, MMG3D, for massively parallel computations. This coupling enables an efficient and parallel remeshing of grids alleviating any memory or performance issues encountered in sequential tools. This strategy is finally applied to the simulation of the isothermal flow in a complex meso‐combustor to demonstrate the applicability of the adaptation methodology to complex turbulent flows. Copyright © 2015 John Wiley & Sons, Ltd.     

Wall modeling for implicit large-eddy simulation and immersed-interface methods

We propose and analyze a wall model based on the turbulent boundary layer equations (TBLE) for implicit large-eddy simulation (LES) of high Reynolds number wall-bounded flows in conjunction with a conservative immersed-interface method for mapping complex boundaries onto Cartesian meshes. Both implicit subgrid-scale model and immersed-interface treatment of boundaries offer high computational efficiency for complex flow configurations. The wall model operates directly on the Cartesian computational mesh without the need for a dual boundary-conforming mesh. The combination of wall model and implicit LES is investigated in detail for turbulent channel flow at friction Reynolds numbers from Re _τ  = 395 up to Re _τ =100,000 on very coarse meshes. The TBLE wall model with implicit LES gives results of better quality than current explicit LES based on eddy viscosity subgrid-scale models with similar wall models. A straightforward formulation of the wall model performs well at moderately large Reynolds numbers. A logarithmic-layer mismatch, observed only at very large Reynolds numbers, is removed by introducing a new structure-based damping function. The performance of the overall approach is assessed for two generic configurations with flow separation: the backward-facing step at Re _h = 5,000 and the periodic hill at Re _H = 10,595 and Re _H = 37,000 on very coarse meshes. The results confirm the observations made for the channel flow with respect to the good prediction quality and indicate that the combination of implicit LES, immersed-interface method, and TBLE-based wall modeling is a viable approach for simulating complex aerodynamic flows at high Reynolds numbers. They also reflect the limitations of TBLE-based wall models.     

A multiscale approach to hybrid RANS/LES wall modeling within a high-order discontinuous Galerkin scheme using function enrichment

We present a novel approach to hybrid Reynolds-averaged Navier-Stokes (RANS)/ large eddy simulation (LES) wall modeling based on function enrichment, which overcomes the common problem of the RANS-LES transition and enables coarse meshes near the boundary. While the concept of function enrichment as an efficient discretization technique for turbulent boundary layers has been proposed in an earlier article by Krank & Wall (A new approach to wall modeling in LES of incompressible flow via function enrichment. J Comput Phys. 2016;316:94-116), the contribution of this work is a rigorous derivation of a new multiscale turbulence modeling approach and a corresponding discontinuous Galerkin discretization scheme. In the near-wall area, the Navier-Stokes equations are explicitly solved for an LES and a RANS component in one single equation. This is done by providing the Galerkin method with an independent set of shape functions for each of these two methods; the standard high-order polynomial basis resolves turbulent eddies, where the mesh is sufficiently fine and the enrichment automatically computes the ensemble-averaged flow if the LES mesh is too coarse. As a result of the derivation, the RANS model is applied solely to the RANS degrees of freedom, which effectively prevents the typical issue of a log-layer mismatch in attached boundary layers. As the full Navier-Stokes equations are solved in the boundary layer, spatial refinement gradually yields wall-resolved LES with exact boundary conditions. Numerical tests show the outstanding characteristics of the wall model regarding grid independence, superiority compared to equilibrium wall models in separated flows, and achieve a speed-up by two orders of magnitude compared to wall-resolved LES.     

Influence of mesh deformation on the quality of large eddy simulations

The influence of mesh motion on the quality of large eddy simulation (LES) was studied in the present article. A three‐dimensional, turbulent pipe flow (Re_τ=360) was considered as a test case. Simulations with both stretching and static meshes were carried out in order to understand how mesh motion affects the turbulence statistics. The spatial filtering of static and moving mesh direct numerical simulation (DNS) data showed how an ideal LES would perform, while the comparison of DNS cases with static and moving meshes revealed that no significant numerical errors arise from the mesh motion when the simulation is fully resolved. The comparison of the filtered fields of the DNS with a moving mesh with the corresponding LES fields revealed different responses to mesh motion from different numerical approaches. A straightforward test was applied in order to verify that the moving mesh works consistently in LES: when the mesh is stretched in the streamwise direction, the moving mesh results should be in between the two extremal resolutions between which the mesh is stretched. Numerical investigations using four different LES approaches were carried out. In addition to the Smagorinsky model, three implicit LES approaches were used: linear interpolation (non‐dissipative), the Gamma limiter (dissipative), and the scale‐selective discretisation (slightly dissipative). The results indicate that while the Smagorinsky and the scale‐selective discretisation approaches produce results consistent with the resolution of the non‐static mesh, the implicit LES with linear interpolation or the Gamma scheme do not. Copyright © 2016 John Wiley & Sons, Ltd.     

Adaptive mesh refinement and load balancing based on multi-level block-structured Cartesian mesh

We developed a framework for a distributed-memory parallel computer that enables dynamic data management for adaptive mesh refinement and load balancing. We employed simple data structure of the building cube method (BCM) where a computational domain is divided into multi-level cubic domains and each cube has the same number of grid points inside, realising a multi-level block-structured Cartesian mesh. Solution adaptive mesh refinement, which works efficiently with the help of the dynamic load balancing, was implemented by dividing cubes based on mesh refinement criteria. The framework was investigated with the Laplace equation in terms of adaptive mesh refinement, load balancing and the parallel efficiency. It was then applied to the incompressible Navier–Stokes equations to simulate a turbulent flow around a sphere. We considered wall-adaptive cube refinement where a non-dimensional wall distance y⁺ near the sphere is used for a criterion of mesh refinement. The result showed the load imbalance due to y⁺ adaptive mesh refinement was corrected by the present approach. To utilise the BCM framework more effectively, we also tested a cube-wise algorithm switching where an explicit and implicit time integration schemes are switched depending on the local Courant-Friedrichs-Lewy (CFL) condition in each cube.     

Implicit LES of free and wall‐bounded turbulent flows based on the discontinuous Galerkin/symmetric interior penalty method

This paper presents the second validation step of a compressible discontinuous Galerkin solver with symmetric interior penalty (DGM/SIP) for the direct numerical simulation (DNS) and the large eddy simulation (LES) of complex flows. The method has already been successfully validated for DNS of an academic flow and has been applied to flows around complex geometries (e.g. airfoils and turbomachinery blades). During these studies, the advantages of the dissipation properties of the method have been highlighted, showing a natural tendency to dissipate only the under‐resolved scales (i.e the smallest scales present on the mesh), leaving the larger scales unaffected. This phenomenon is further enhanced as the polynomial order is increased. Indeed, the order increases the dissipation at the largest wave numbers, while its range of impact is reduced. These properties are spectrally compatible with a subgrid‐scale model, and hence DGM may be well suited to be used for an implicit LES (ILES) approach. A validation of this DGM/ILES approach is here investigated on canonical flows, allowing to study the impact of the discretisation on the turbulence for under‐resolved computations. The first test case is the LES of decaying homogeneous isotropic turbulence (HIT) at very high Reynolds number. This benchmark allows to assess the spectral behaviour of the method for implicit LES. The results are in agreement with theory and are even slightly more accurate than other numerical results from literature, obtained using a pseudo‐spectral (PS) method with a state‐of‐the‐art subgrid‐scale model. The second benchmark is the LES of the channel flow. Three Reynolds numbers are considered: Re_τ=395, 590 and 950. The results are compared with DNS of Moser et al. and Hoyas et al., also using PS methods. Both averaged velocity and fluctuations are globally in good agreement with the reference, showing the ability of the method to predict equilibrium wall‐bounded flow turbulence. To show that the method is able to perform accurate DNS, a DNS of HIT at Re_λ=64 and a DNS of the channel flow at Re_τ=180 are also performed. The effects of the grid refinement are investigated on the channel flow at Re_τ=395, highlighting the improvement of the results when refining the mesh in the spanwise direction. Finally, the modification of the ILES parameters, that is the Riemann solver and of the SIP coefficient, is studied on both cases, showing a significant influence on the choice of the Riemann solver. Copyright © 2015 John Wiley & Sons, Ltd.     

Adaptive mesh refinement techniques for high‐order shock capturing schemes for multi‐dimensional hydrodynamic simulations

The numerical simulation of physical phenomena represented by non‐linear hyperbolic systems of conservation laws presents specific difficulties mainly due to the presence of discontinuities in the solution. State of the art methods for the solution of such equations involve high resolution shock capturing schemes, which are able to produce sharp profiles at the discontinuities and high accuracy in smooth regions, together with some kind of grid adaption, which reduces the computational cost by using finer grids near the discontinuities and coarser grids in smooth regions. The combination of both techniques presents intrinsic numerical and programming difficulties. In this work we present a method obtained by the combination of a high‐order shock capturing scheme, built from Shu–Osher's conservative formulation (J. Comput. Phys. 1988; 77 :439–471; 1989; 83 :32–78), a fifth‐order weighted essentially non‐oscillatory (WENO) interpolatory technique (J. Comput. Phys. 1996; 126 :202–228) and Donat–Marquina's flux‐splitting method (J. Comput. Phys. 1996; 125 :42–58), with the adaptive mesh refinement (AMR) technique of Berger and collaborators (Adaptive mesh refinement for hyperbolic partial differential equations. Ph.D. Thesis, Computer Science Department, Stanford University, 1982; J. Comput. Phys. 1989; 82 :64–84; 1984; 53 :484–512). Copyright © 2006 John Wiley & Sons, Ltd.     

The direct simulation Monte Carlo method using unstructured adaptive mesh and its application

The implementation of an adaptive mesh‐embedding (h‐refinement) scheme using unstructured grid in two‐dimensional direct simulation Monte Carlo (DSMC) method is reported. In this technique, local isotropic refinement is used to introduce new mesh where the local cell Knudsen number is less than some preset value. This simple scheme, however, has several severe consequences affecting the performance of the DSMC method. Thus, we have applied a technique to remove the hanging node, by introducing the an‐isotropic refinement in the interfacial cells between refined and non‐refined cells. Not only does this remedy increase a negligible amount of work, but it also removes all the difficulties presented in the originals scheme. We have tested the proposed scheme for argon gas in a high‐speed driven cavity flow. The results show an improved flow resolution as compared with that of un‐adaptive mesh. Finally, we have used triangular adaptive mesh to compute a near‐continuum gas flow, a hypersonic flow over a cylinder. The results show fairly good agreement with previous studies. In summary, the proposed simple mesh adaptation is very useful in computing rarefied gas flows, which involve both complicated geometry and highly non‐uniform density variations throughout the flow field. Copyright © 2002 John Wiley & Sons, Ltd.     

A high‐order flux reconstruction adaptive mesh refinement method for magnetohydrodynamics on unstructured grids

We report our recent development of the high‐order flux reconstruction adaptive mesh refinement (AMR) method for magnetohydrodynamics (MHD). The resulted framework features a shock‐capturing duo of AMR and artificial resistivity (AR), which can robustly capture shocks and rotational and contact discontinuities with a fraction of the cell counts that are usually required. In our previous paper, 36 we have presented a shock‐capturing framework on hydrodynamic problems with artificial diffusivity and AMR. Our AMR approach features a tree‐free, direct‐addressing approach in retrieving data across multiple levels of refinement. In this article, we report an extension to MHD systems that retains the flexibility of using unstructured grids. The challenges due to complex shock structures and divergence‐free constraint of magnetic field are more difficult to deal with than those of hydrodynamic systems. The accuracy of our solver hinges on 2 properties to achieve high‐order accuracy on MHD systems: removing the divergence error thoroughly and resolving discontinuities accurately. A hyperbolic divergence cleaning method with multiple subiterations is used for the first task. This method drives away the divergence error and preserves conservative forms of the governing equations. The subiteration can be accelerated by absorbing a pseudo time step into the wave speed coefficient, therefore enjoys a relaxed CFL condition. The AMR method rallies multiple levels of refined cells around various shock discontinuities, and it coordinates with the AR method to obtain sharp shock profiles. The physically consistent AR method localizes discontinuities and damps the spurious oscillation arising in the curl of the magnetic field. The effectiveness of the AMR and AR combination is demonstrated to be much more powerful than simply adding AR on finer and finer mesh, since the AMR steeply reduces the required amount of AR and confines the added artificial diffusivity and resistivity to a narrower and narrower region. We are able to verify the designed high‐order accuracy in space by using smooth flow test problems on unstructured grids. The efficiency and robustness of this framework are fully demonstrated through a number of two‐dimensional nonsmooth ideal MHD tests. We also successfully demonstrate that the AMR method can help significantly save computational cost for the Orszag‐Tang vortex problem.     

Application of adaptive mesh refinement method to complex flows in clean rooms

The adaptive mesh refinement (AMR) method is developed for three-dimensional turbulent complex flows in clean rooms using the finite volume method with a collocated grid arrangement. Clean rooms have many interesting and complex flow characteristics especially the secondary flows and the recirculation regions. The accurate numerical solution of the flows is important for the efficient design of clean rooms. The use of the conventional uniform grid requires such a high computational time and data storage capacity that they make computational fluid dynamics (CFD) less attractive for the design optimization. The AMR method is, therefore, applied by using the fine grid only in the required regions and using the coarse grid in the other regions. The velocity is chosen as the main parameter for the grid refinement because it is the most influential parameter in clean rooms. The results show that the present AMR method can reduce the computational time by eight times and the data storage requirement is only 37% of that using the conventional method, while the same order of accuracy can be maintained. The present AMR method is, therefore, proved to be a promising technique for solving three-dimensional turbulent complex flows in clean rooms.     

Solution adaptive grids applied to low Reynolds number flow

A numerical study has been undertaken to investigate the use of a solution adaptive grid for flow around a cylinder in the laminar flow regime. The main purpose of this work is twofold. The first aim is to investigate the suitability of a grid adaptation algorithm and the reduction in mesh size that can be obtained. Secondly, the uniform asymmetric flow structures are ideal to validate the mesh structures due to mesh refinement and consequently the selected refinement criteria. The refinement variable used in this work is a product of the rate of strain and the mesh cell size, and contains two variables C_m and C_str which determine the order of each term. By altering the order of either one of these terms the refinement behaviour can be modified. Copyright © 2003 John Wiley & Sons, Ltd.     

A novel Cartesian cut cell solver for incompressible viscous flow in irregular domains

A Cartesian cut cell solver with solution‐based adaptive mesh refinement is developed for simulating viscous, incompressible flows with arbitrary complex geometries. The cut cells are automatically generated using Volume CAD (VCAD), a framework for storing geometric and material attribute data. Unlike earlier cut cell methods, this solver organizes the cutting patterns into only six categories and further subdivides the resulting pentagon into two quadrilaterals, such that mesh data can be stored by uniform data structure and the post‐processing of flow data can be handled conveniently. A novel method is proposed to treat minuscule cut cells without the process of cell merging. A collocated finite volume method, which can be used even when multiple cell shapes and orthogonal and non‐orthogonal grids exist in the decomposition, is employed to discretize the Navier–Stokes equations. A modified SIMPLE‐based smoothing pressure correction scheme is applied in this cut cell method to suppress checkerboard pressure oscillations caused by collocated arrangement. The solver is first used to simulate a channel flow to demonstrate its calculation accuracy expressed with L₁ and L_∞ norm errors and then the method is utilized to solve three benchmark problems of flow and heat transfer within irregular domains to verify its feasibility, efficiency, accuracy and potential in engineering applications. Copyright © 2010 John Wiley & Sons, Ltd.     

The use of dynamic grid adaptation algorithms for the modelling of flow around a circular cylinder in sub‐critical flow regime

In the present study a dynamic grid adaptation (DGA) algorithm is utilized for predicting flow around a circular cylinder in sub‐critical flow regime at a Reynolds number of 1.4×10⁵. The reason for adopting a DGA algorithm is the unsteadiness of the flow field which makes a conventional mesh inefficient. The concept being adopted is to concentrate mesh refinement in regions with high gradients and high turbulent viscosity, while in the region further downstream where the flow is fully developed a coarser mesh will develop and turbulence is modelled with the large eddy simulation (LES) turbulence model. The aim of the study is to present an appropriate variable for mesh refinement, which accomplishes a high rate of mesh refinement in the region with high gradients. The new variable is a product of the local mesh cell size and the rate of strain and includes two additional variables to allow control over the refinement behaviour. The results are compared with experimental data at the corresponding Reynolds number and also with numerical results obtained with conventional mesh. It is demonstrated that DGA algorithms can give results of a very high quality for a mesh that is significantly smaller than for a conventional mesh. Copyright © 2003 John Wiley & Sons, Ltd.     

A Mesh Adaptation Strategy to Predict Pressure Losses in LES of Swirled Flows

Large-Eddy Simulation (LES) has become a potent tool to investigate instabilities in swirl flows even for complex, industrial geometries. However, the accurate prediction of pressure losses on these complex flows remains difficult. The paper identifies localised near-wall resolution issues as an important factor to improve accuracy and proposes a solution with an adaptive mesh h-refinement strategy relying on the tetrahedral fully automatic MMG3D library of Dapogny et al. (J. Comput. Phys. 262, 358-378, 2014) using a novel sensor based on the dissipation of kinetic energy. Using a joint experimental and numerical LES study, the methodology is first validated on a simple diaphragm flow before to be applied on a swirler with two counter-rotating passages. The results demonstrate that the new sensor and adaptation approach can effectively produce the desired local mesh refinement to match the target losses, measured experimentally. Results shows that the accuracy of pressure losses prediction is mainly controlled by the mesh quality and density in the swirler passages. The refinement also improves the computed velocity and turbulence profiles at the swirler outlet, compared to PIV results. The significant improvement of results confirms that the sensor is able to identify the relevant physics of turbulent flows that is essential for the overall accuracy of LES. Finally, in the appendix, an additional comparison of the sensor fields on tetrahedral and hexahedral meshes demonstrates that the methodology is broadly applicable to all mesh types.     

Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.