Unsteady Incompressible Flow Research Articles

Influence of internal flow structure on flow energy losses in a centrifugal pump with splitter blades using entropy generation theory

Centrifugal pumps are essential in various industrial applications, including renewable energy. To maximize the pump’s overall performance, estimating flow energy losses (FEL) is crucial. However, due to internal flow complexity, it is necessary to employ new methods and approaches to better understand FEL mechanisms. This study uses entropy generation theory to investigate the relationship between FEL and various flow patterns in all pump hydraulic components. Numerical simulations were conducted using a 3-dimensional incompressible unsteady flow at four different flow coefficients. The delayed detached eddy simulation (DDES) model was used, and the results showed good agreement with experimental data compared to the SST k-w model. Emphasis was given to the flow energy losses and the relative and absolute velocity distribution in the impeller and diffuser components at various flow coefficients. Flow energy losses primarily occur in the impeller (70%) followed by the diffuser (23%). Impeller losses are concentrated at the outlet region due to the wake-jet phenomena (28.6%), splitter blades region (15.3%), and impeller’s leading edge (LE) (10.7%). Vaned diffuser losses occur in the vanless zone (stator-rotor interaction) and near leading/trailing edges.Moreover, wall shear stress and the significant relative velocity gradient near the walls of the impeller and diffuser blades are the main contributors to the FEL in this region. This study provides insights into a better understanding of FEL mechanisms and highlights areas for improving pump performance.

Study on the integrated self-priming process of a vehicle-mounted self-priming pump

In order to explore the self-priming characteristics of the self-priming pump at the mobile pump truck, this paper established a complete three-dimensional circulatory piping system including the self-priming pump, tank, valves, inlet pipe and outlet pipe. The UDF(User Defined Functions) was used to realize the acceleration-constant speed operation process of the impeller, thus reflecting the actual changing state of the rotational speed. Based on the VOF(Volume Of Fluid) multiphase flow model and the Realizable k-ε turbulence model, a coupled numerical calculation of unsteady incompressible viscous flow was conducted for its self-priming process. The results show that the self-priming process of the pump can be roughly divided into four stages: the rapid suction stage, the shock exhaust stage, the rapid exhaust period and the pump residual gas discharge stage. The proportion of each stage in the total self-priming time showed an increasing trend. During the rapid suction stage, the water level in the vertical section of the inlet pipe showed a slow and then fast-rising pattern. During the shock exhaust stage, the average gas-phase volume fraction in the volute is lower than that of the impeller, and the gas content at the volute outlet is lower than that of the impeller inlet. The region at the inlet and outer edge of the impeller consistently experience significant energy losses.

A high‐order pure streamfunction method in general curvilinear coordinates for unsteady incompressible viscous flow with complex geometry

AbstractIn this paper, a high‐order compact finite difference method in general curvilinear coordinates is proposed for solving unsteady incompressible Navier‐Stokes equations. By constructing the fourth‐order spatial discretization schemes for all partial derivative terms of the pure streamfunction formulation in general curvilinear coordinates, especially for the fourth‐order mixed derivative terms, and applying a Crank‐Nicolson scheme for the second‐order temporal discretization, we extend the unsteady high‐order pure streamfunction algorithm to flow problems with more general non‐conformal grids. Furthermore, the stability of the newly proposed method for the linear model is validated by von‐Neumann linear stability analysis. Five numerical experiments are conducted to verify the accuracy and robustness of the proposed method. The results show that our method not only effectively solves problems with non‐conformal grids, but also allows grid generation and local refinement using commercial software. The solutions are in good agreement with the established numerical and experimental results.

Chaotic Vortex-Induced Vibrations of Rigid Cylinders with Nonlinear Snapping Support

Vortex-Induced Vibration (VIV) is a complex fluid–structure interaction in offshore structures. Traditionally, this phenomenon is considered periodic; however, many of its signals show chaotic behavior. The basic model already employed by other researchers is a rigid circular cylinder with linear springs and dampers. In this work, nonlinear snapping support is used to model nonlinearity in the system. To numerically simulate the flow, Reynolds-Averaged Navier–Stokes (RANS) equations for two-dimensional incompressible unsteady flows are applied. The degree of nonlinearity of the system can be changed by manipulating [Formula: see text], which is one of the geometric properties of the spring and takes values between 0 and 1. The 0–1 test, Poincaré section, and Fast Fourier Transform are used to analyze the cylinder and lift force behavior. Also, the Hilbert transform is applied to the signals, and the phase shift between displacement and lift force is obtained. The results show that the system behavior consists of branches: branch I and branch II. The large amplitudes occur in branch II. It is found that chaos emerges at the beginning of branch II, regardless of the value of [Formula: see text]. By raising the [Formula: see text] value, the span of branch II becomes more expansive, and its first point is placed at lower reduced velocities. Also, the wake dynamics becomes more regular at the end of branch I and more irregular at the beginning of branch II with the increase in [Formula: see text]. When the cylinder displacement signal is chaotic, the lift force behavior is also chaotic, but not vice versa.

Development of an accurate central finite-difference scheme with a compact stencil for the simulation of unsteady incompressible flows on staggered orthogonal grids

In scale-resolving simulations such as Large Eddy Simulations (LES), the spatial discretization scheme of the convective term plays a crucial role in avoiding interference between the numerical errors and the subgrid-scale model. Accurate schemes lead to lower truncation errors and better predictions of turbulent flows without the need for an excessively refined grid. To this end, a new second-order finite-difference scheme (HCDS6) has been developed for incompressible flows and orthogonal staggered grids. Compared to the standard second-order scheme, the new scheme has significantly lower dispersion errors. Compared to existing high-order schemes, the numerical stencil of HCDS6 is more compact, which makes it easier to implement, especially considering boundary conditions around complex geometries using Immersed Boundary Methods (IBM). The HCDS6 scheme conserves the discrete momentum with limited production or dissipation of discrete kinetic energy, which guarantees its numerical stability. Its performance is evaluated using an open-source CFD package called REEF3D. Three benchmarks demonstrate the key properties and performance of the scheme: the convection of an isentropic vortex, the Taylor-Green vortex flow, and turbulent channel flow. Its relatively low dispersion errors, combined with ease of implementation, make the HCDS6 scheme a promising candidate for efficient scale-resolving simulations of turbulent flows.

Investigation of MHD fractionalized viscous fluid and thermal memory with slip and Newtonian heating effect: a fractional model based on Mittag-Leffler kernel

This paper introduces an innovative approach for modelling unsteady incompressible natural convection flow over an inclined oscillating plate with an inclined magnetic effect that employs the Atangana-Baleanu time-fractional derivative (having a non-singular and non-local kernel) and the Mittag-Leffler function. The fractional model, which includes Fourier and Fick's equations, investigates memory effects and is solved using the Laplace transform. The Mittag-Leffler function captures power-law relaxation dynamics, which improves our understanding of thermal and fluid behaviour. Graphical examination shows the influence of fractional and physically involved parameters, leading to the conclusion that concentration, temperature, and velocity profiles initially grow and then decrease asymptotically with time. Moreover, the study emphasizes the impact of effective Prandtl and Schmidt numbers on temperature, concentration, and velocity levels in the fluid.

Soret and Dufour effects on an unsteady MHD flow about a permeable rotating vertical cone with variable fluid properties

The objective of the present work is to examine the characteristics of unsteady incompressible magnetohydrodynamic fluid flow around a permeable rotating vertical cone. The effects of thermal radiation, viscous dissipation, and the Soret and Dufour effects are investigated in the analysis of heat and mass transfer. The viscosity of the fluid is considered inversely proportional to the temperature, and the thermal conductivity of the fluid is considered directly proportional to the temper-ature. The governing equations are converted into ordinary differential equations using suitable similarity transformations, which are then solved numerically using bvp4c from MATLAB. Results obtained in this study are in excellent correlation with previously conducted studies. The results demonstrate that the Dufour and Soret effects subsequently reduce the heat transit rate (by –3.3%) and mass transit rate (by –1.2%) of the system. It is also detected that fluids with higher viscosity tend to increase tangential skin friction (+8.9%) and azimuthal skin friction (+8.3%). The heat transit rate of the system is found to be more efficient for fluids with higher viscosity and lower thermal conductivity and Eckert numbers. Further-more, the thickness of the momentum, thermal, and concentration boundary layers significantly reduces while the heat and mass transit rates (+17.8% and +18.3%, respectively) of the system become more efficient for greater values of the un-steadiness parameter.

Series solution of time-fractional mhd viscoelastic model through non-local kernel approach

The study of ramped condition in the context of unsteady incompressible magnetohydrodynamic Casson fluid flow over a moving vertical plate is a complex and important topic in fluid dynamics and heat transfer. This scenario combines several physical phenomena and has practical applications in various engineering and scientific fields. In this study, Casson fluid is considered unsteady under the influence of magnetic field. The fractional mathematical model is proposed by considering the effect of chemical reaction parameter of the flowing fluid. The governing equations are transformed into the dimensionless form and developed fractional models like Caputo-Fabrizio and Atangana-Baleanu Derivative. We used the Laplace transform technique to find the solution of the dimensionless governing equation analytically. The transformed solutions for velocity, energy and momentum balances developed in terms of series. MATHCAD software is being used for numerical computations and the physical attributes of material and fractional parameters are discussed. To analyze their behavior clearly, two-dimensional graphical results are plotted for velocity profile and temperature as well. It has been concluded that the fluid’s velocity are reduced for larger values of the fractional parameter and Prandtl number and is maximum for small values of both parameters. Further, the velocity behavior becomes larger for isothermal condition as compared to ramped conditions.

Optimization of mass and heat flux of MHD viscous fluid flow with constant proportional Caputo derivative by using response surface methodology: Sensitivity analysis

The study of viscous dissipation in the context of unsteady incompressible magnetohydrodynamic viscous fluid flow over a moving plate is a complex and important topic in fluid dynamics and heat transfer. This scenario combines several physical phenomena and has practical applications in various engineering and scientific fields. The fractional mathematical model is proposed by considering the effect of viscous dissipation force of the flowing fluid. For this purpose, fluid is assumed to be lying over a vertical plate that is moving in its own plane. The governing equations are transformed into the dimensionless form and developed fractional model with a new operator Constant Proportional Caputo Derivative. We used the Laplace transform technique to find the solution of the dimensionless governing equation analytically. The transformed solutions for energy and momentum balances developed in terms of series. The validation of the present model comparison graphs is plotted for temperature. The sensitivity of different input parameters on output response are analyzed and displayed the sensitivity results in tabular and graphical form and found that the permeable parameter is most sensitive to skin friction coefficient, and the Eckert number is most sensitive parameter among others to Nusselt Number. The novelty of present work is to develop a correlation between input parameters and output responses by using Response Surface Methodology (RSM). As no such correlation has been developed for optimization of mass and heat flux of MHD viscous fluid flow with Constant Proportional Caputo derivative by using RSM. By using this correlation, we have performed sensitivity analysis for output responses.

Numerical investigations on the entropy generation and heat transfer of an unsteady micropolar hybrid nanofluid flow over an inclined stretching surface

Studies on stagnation point flow on inclined surfaces have gained momentum in recent times due to its application in a variety of fields. This study aims at the heat transfer properties of an incompressible unsteady stagnation point flow of a non-Newtonian micropolar hybrid nanofluid over an inclined stretching sheet with thermal radiation, viscous dissipation, and Joule heating. The governing equations of the problem are transformed into a system of ordinary differential equations by employing similarity variables. They are then solved numerically using the Runge-Kutta 4 th order method along with the shooting approach. Numerical solutions on velocity, microrotation, and temperature distribution along with entropy generation and Bejan number are consolidated graphically and those on heat transfer rate and skin friction coefficient are tabulated. The results indicate that the heat transfer rate is higher in hybrid nanofluid as compared to that of mononanofluid, the angle of inclination of the sheet is unfavorable for the velocity profile, and the temperature of the hybrid nanofluid escalates with a rise in viscous dissipation. Entropy generation rate is reduced when the Hartmann number is increased and the Bejan number diminishes with elevating viscous dissipation. The results were validated using existing studies and comparing heat transfer rate values between the present results and the results obtained from NDSolve.

Machine learning mesh-adaptation for laminar and turbulent flows: applications to high-order discontinuous Galerkin solvers

We present a machine learning-based mesh refinement technique for steady and unsteady incompressible flows. The clustering technique proposed by Otmani et al. (Phys Fluids 35(2):027112, 2023) is used to mark the viscous and turbulent regions for the flow past a cylinder at Re=40\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Re=40$$\\end{document} (steady laminar flow), at Re=100\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Re=100$$\\end{document} (unsteady laminar flow), and at Re=3900\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Re=3900$$\\end{document} (unsteady turbulent flow). Within this clustered region, we use high mesh resolution, while downgrading the resolution outside, to show that it is possible to obtain levels of accuracy similar to those obtained when using a uniformly refined mesh. The mesh adaptation is effective, as the clustering successfully identifies the two flow regions, a viscous/turbulent dominated region (including the boundary layer and wake) that requires high resolution and an inviscid/irrotational region, which only requires low resolution. The new clustering sensor is compared with traditional feature-based sensors (Q-criterion and vorticity based) commonly used for mesh adaptation. Unlike traditional sensors that rely on problem-dependent thresholds, our novel approach eliminates the need for such thresholds and locates the regions that require adaptation. After the initial validation using flows past cylinders, the clustering technique is applied in an engineering context to study the flow around a horizontal axis wind turbine configuration which has been tested experimentally at the Norwegian University of Science and Technology. The data used within this framework are generated using a high-order discontinuous Galerkin solver, allowing to locally refine the polynomial order (p-refinement) in each element of the clustered region. For the laminar test cases, we can reduce the computational cost by 32% (steady Re=40\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Re=40$$\\end{document} case) and 20% (unsteady Re=100\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Re=100$$\\end{document} case), while we get a reduction of 33% for the Re=3900\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Re=3900$$\\end{document} turbulent case. In the context of the wind turbine, a reduction of 43% in computational cost is observed, while maintaining the accuracy.

Unsteady Compressed Williamson Fluid Flow Behavior under the Influence of a Fixed Magnetic Field (Numerical Study)

A numerical investigation is conducted into a two-dimensional mathematical model of magnetized unsteady incompressible Williamson fluid flow over a sensor surface with fixed thermal conductivity and external squeezing accompanied by viscous dissipation effect. Based on the flow geometry under consideration, the current flow model was created. The momentum equation takes into consideration the magnetic field when describing the impact of Lorentz forces on flow behavior. The energy equation takes varying thermal conductivity into account while calculating heat transmission. The extremely complex nonlinear, unstable governing flow equations for the now under investigation are coupled in nature. Due to the inability of analytical or direct methods, the Runge-Kutta scheme (RK-4) via similarity transformations approach is used to tackle the physical problem under consideration. The physical behavior of various control factors on the flow phenomena is described using graphs and tables. For increasing values of the Weissenberg parameter and the permeable velocity parameter, the temperature boundary layer thickens. As the permeable velocity parameter and squeezed flow index increased, the velocity profile shrank. The velocity profile grows as the magnetic number rises. Squeezed flow magnifying increases the Nusselt number's magnitude. Furthermore, the extremely complex nonlinear complex equations that arise in fluid flow issues are quickly solved by RK-4. The current findings in this article closely align with the findings that have been reported in the literature.

Enhancing the performance of vertical axis hydrokinetic Savonius turbines using a novel cambered hydrofoil profile for rotor blades

A novel rotor blade configuration using a cambered hydrofoil profile is developed to enhance the output power. Several design variables are considered including camber angle, overlap ratio, and camber position, in addition to integrating endplates. To assess the performance parameters at different design conditions, a comprehensive three-dimensional unsteady incompressible turbulent flow model based on Reynolds-Averaged Navier-Stokes equations combined with k-ω model is developed. The predicted numerical results is validated using the available numerical and experimental data. The findings indicate that using Savonius rotor with hydrofoil profile results in a substantial increase in the maximum output power for all investigated camber angles compared to the conventional semi-circular rotor. At a flow velocity of 0.4 m/s and an overlap ratio of 0.15, the Savonius rotor with hydrofoil profile attains a maximum power coefficient of 0.21 at a camber angle of 105°. By increasing the camber angle to 140°, the maximum power coefficient reaches 0.26 whereas it is 0.13 for the conventional design. The camber position at 50 % of the chord length achieves the highest power coefficient for all studied cases. In addition, adding the endplates significantly enhances the power coefficient. The current findings provide a promising approach in developing Savonius rotor blades configuration by utilizing hydrofoil profile.

Multiple solutions of unsteady flow of CNTs nanofluid over permeable shrinking surface with effects of dissipation and slip conditions

The objective of this paper is to analyze the unsteady incompressible flow of the viscous nanofluid on a contracting surface with viscous dissipation effects. Presented and contrasted are analyses of both multi-wall carbon nanotubes (MWNTs) and single-wall carbon nanotubes (SWNTs). As the common (or base) fluids, kerosene oil and water are utilized. In the existence of first-order thermal and velocity slip conditions, mathematical modeling and analysis are performed. Using the MATLAB software’s bvp4c solver tool, numerical solutions to the governing nonlinear modeled problems were obtained. This technique is particularly effective for developing many solutions to highly nonlinear differential equations. In addition, a comparison is done between this study and previously published works. The temperature, velocity, skin friction coefficient and heat-transfer rate have been explored for various significant factors included in the problem statements. In the unsteadiness parameter regime, dual solutions can be found. As the velocity slip parameter is increased, the flow slows down. In comparison to SWCNTs kerosene, MWCNTs kerosene oil has a greater velocity curve for the nanoparticles volume fraction. Increases in volume fraction decrease skin friction, whereas increases in the unsteadiness parameter speed up the drag force. Furthermore, as the Eckert number intensity increases, so do the temperature profiles in both solutions. Finally, the stability study revealed that the initial solution is robust, whereas the breakage in the second solution in the Nusselt number shows singularity, and thus the second solution is considered unstable.

Unsteady Casson fluid flow on MHD with an internal heat source

The present paper investigates the heat transfer features of a viscous fluid in an unsteady incompressible Casson fluid flow through a porous medium with an internal heat source under the influence of a magnetic field. A Casson fluid passes through the given porous medium using Darcy model. The study involves simulating the flow using mathematical models based on conservation laws and similarity transformations. The resulting highly nonlinear coupled ordinary differential equations are solved using the Bernoulli wavelet numerical approach. The study utilizes graphical tools to illustrate the flow behaviour under different physical factors, including the effect of the internal heat source. The normalization of the flow behaviour by the magnetic field suggests that it could potentially be used to control various flows. The study finds that various physical parameters which affects the velocity and temperature profiles of the fluid.

Solute Dispersion in an Unsteady Herschel-Bulkley Flow through an Inclined Stenosed Arter

Motivated by the concept of blood flow in a stenosed artery, this present research investigates the influence of stenosis shape in terms of height and arterial inclination on the blood flow and solute dispersion behaviour through an inclined stenosed artery. The blood rheology is depicted using the Herschel-Bulkley model in a laminar, axisymmetric and incompressible unsteady flow through the stenosed artery. The effect of stenosis is focused on the stenosis height for both sine and cosine stenosis. Parameters of arterial inclination are also investigated to observe the effect of inclination on the blood velocity and dispersion function. Perturbation method is adopted in solving for the blood flow velocity under the effect of stenosis height and arterial inclination. The dispersion function of solute dispersion is solved using the obtained blood velocity by adopting the Generalized Dispersion Model (GDM) in obtaining steady dispersion functions. This present study shows that the increase in stenosis height decreases both blood velocity and dispersion function. Meanwhile, the increase in arterial inclination increases the blood velocity and dispersion function. The effect of stenosis height also affects blood velocity and dispersion function for the sine stenosis more than the cosine stenosis.

Non-relaxed finite volume fractional step schemes for unsteady incompressible flows

Non-relaxed finite volume fractional step schemes for unsteady incompressible flows

CFD Simulation of Low Speed Turbulent Flow Problems Using Unsteady Rans and Large Eddy Simulation Approach

The present paper focuses on the recent development of an implicit pressure-based finite volume algorithm for numerical solution of Navier Stokes equation in an inertial frame of reference for prediction of unsteady incompressible flow problems. The algorithm uses boundary-conforming, multiblock structured grid with moving boundaries, collocated variable arrangement with momentum equations resolved along cartesian directions, second order accurate spatial and temporal discretisation schemes for the convective fluxes and a pressurevelocity solution strategy. Effect of turbulence is simulated using one of the two different approaches. In the Unsteady Reynolds Averaged (URANS) approach coupled to appropriate eddy viscosity based turbulence models, the Navier Stokes (NS) equations, averaged over the whole range of turbulent length scales of the flow, are solved numerically. On the other hand, in the Large Eddy Simulation (LES) approach, the model filtered 3D NS equations are directly solved for the flow variables to resolve the large scale turbulent motions whereas the transport processes at the fine subgrid scale level only are simulated using simple algebraic turbulence model. The capabilities and limitations of both the cost-effective URANS approach and the relatively expensive but rich in physics LES approach have been demonstrated for few application problems of engineering interest.

Unsteady incompressible flow over a backward-facing step using FEniCS

This study demonstrates unsteady incompressible flow over a backward-facing step and a cylinder section using a finite element scheme for solving the unsteady Navier-Stokes equation called the incremental pressure correction scheme (IPCS). The speed and accuracy of the study were boosted by using an open-source framework. The formulation was done using the C++/Python-based computing platform called FEniCS. Gmsh is the mesh generator used to make the meshes. Krylov solvers, parallel processing, and higher-order function spaces of order above two are available in FEniCS and these enable faster and more accurate analyses. Most commercial software packages do not include these functionalities. The analysis of the backward-facing step was performed using legacy FEniCS. The DFG benchmark 2D-3 unsteady analysis in Featflow was conducted using FEniCSx. The pressure difference between the front and back of the cylinder and the drag and lift coefficients were calculated and compared to the benchmark. The graphs of these characteristics were plotted with respect to time along with their mean squared errors. The outcomes demonstrate how efficient the open-source FEM platform on display is at handling fluid flow simulations. Therefore, any turbulent incompressible flow may be simulated using the same open-source structure. [copyright information to be updated in the production process]

A three-dimensional high-order flux reconstruction lattice boltzmann flux solver for incompressible laminar and turbulent flows

A three-dimensional (3D) high-order flux reconstruction lattice Boltzmann flux solver (FR-LBFS) is developed in this paper for accurate and efficient simulations of incompressible laminar and turbulent flows. The original lattice Boltzmann flux solver (LBFS) is a second-order finite volume method (FVM) which combines the advantages of conventional Navier-Stokes (N-S) solvers and lattice Boltzmann equation (LBE) solvers. But it is not convenient to construct high-order FVM, and the compactness is always a bottleneck. The present work separates LBFS from FVM and combines LBFS with high-order FR scheme. Thus, the FR-LBFS inherits the superiority of LBFS and FR and shares the following attractive features: (i) a fully-explicit arbitrary order compact method, (ii) weakly compressible (WC) model for unsteady incompressible flows and (iii) evaluating inviscid and viscous fluxes simultaneously and no particular technique requirement, such as local discontinuous Galerkin method (LDG), for viscous term discretization. To validate the current method, various 3D incompressible laminar isothermal and thermal numerical experiments are presented first. Good agreements are achieved between the current results and the previous experimental and numerical data whilst the present high-order solver adopt coarse meshes. Furthermore, the ability of the proposed method to perform accurate and stable computations of incompressible decaying homogeneous isotropic turbulent flows and turbulent channel flows with different mesh resolutions and accuracy orders are investigated. The results show that the present 3D FR-LBFS is a promising tool for under-resolved or implicit large eddy simulation (ILES) of incompressible turbulent flows.