Identification of local extinction topology in axisymmetric bluff‐body diffusion flames with a reactedness‐mixture fraction presumed probability density function model

The effects of finite‐rate chemistry, such as partial extinctions and re‐ignitions, are investigated in turbulent non‐pre‐mixed reacting flows stabilized in the wake of an axisymmetric bluff‐body burner. A two‐dimensional large‐eddy simulation procedure is employed that uses a partial equilibrium/two‐scalar reactedness mixture fraction probability density function (PDF) combustion sub‐model, which is applied at the sub‐grid scale (SGS) level. An anisotropic sub‐grid eddy–viscosity and two equations for the SGS turbulence kinetic and scalar energies complete the SGS closure model. The scalar covariances required in the joint PDF formulation are obtained from an extended scale‐similarity assumption between the resolved and the sub‐grid fluctuations. Extinction due to strong turbulence/chemistry interactions is recognized with the help of a ‘critical’, locally variable, turbulent Damkohler number criterion, while transient localized extinctions and re‐ignitions are treated with a Lagrangian transport equation for a reactedness progress variable. Comparisons with available experimental data suggested that the formulated approach was capable of identifying the effects of large‐scale vortex structure activity, which were inherent in the reacting wake and dominant in the counterpart isothermal flows that otherwise would have been obscured if a standard time‐averaged procedure had been used. Additionally, the post‐extinction and re‐ignition behaviour and its time‐varying interaction with the large‐scale structure dynamics were more appropriately addressed within the context of the present time‐dependent method. Copyright © 2001 John Wiley & Sons, Ltd.     

3D Numerical Simulation of a Laminar Experimental SWQ Burner with Tabulated Chemistry

Flame-wall interaction (FWI) plays an important role in enclosed combustion systems. For avoiding the complexity of close to reality combustors, in this study an atmospheric premixed V-shaped flame interacting with an isothermal cold wall in a side wall quenching (SWQ) configuration is investigated. A stoichiometric methane/air mixture is used as fuel. A three-dimensional (3D) numerical simulation, which resolves all flow structures is combined with a tabulated chemistry approach (flamelet generated manifold, FGM). Results are compared with experimental data and two-dimensional simulations. The FGM approach is a suitable trade-off between computationally expensive detailed chemistry simulations and over simplified single step mechanisms. 2D simulations are used to investigate the influence of the uncertainty of the wall temperature, to show that the resolution in 3D is sufficient and that the influence of the flame thickening on the wall heat fluxes can be determined. Our results show that the 3D FGM approach is in close agreement to experimentally obtained flow and temperature fields. The dimensionless wall heat flux and Péclet number matches the expected values of 0.16 and 7, respectively. However, during FWI the measured CO mole fractions are not reproduced accurately showing that the transported variables in the present approach of tabulated chemistry do not recover premixed flame structures near walls.     

Detached eddy simulation of cyclic large scale fluctuations in a simplified engine setup

Computational Fluid Dynamics using RANS-based modelling approaches have become an important tool in the internal combustion engine development and optimization process. However, these models cannot resolve cycle to cycle variations, which are an important aspect in the design of new combustion systems. In this study the feasibility of using a Detached Eddy Simulation (DES) SST model, which is a hybrid RANS/LES model, to predict cycle to cycle variations is investigated. In the near wall region or in regions where the grid resolution is not sufficiently fine to resolve smaller structures, the two-equation RANS SST model is used. In the other regions with higher grid resolution an LES model is applied. The case considered is a geometrically simplified engine, for which detailed experimental data for the ensemble averaged and single cycle velocity field are available from Boreé et al. [Boreé, J., Maurel, S., Bazile, R., 2002. Disruption of a compressed vortex, Physics of Fluids 14 (7), 2543–2556]. The fluid flow shows a strong tumbling motion, which is a major characteristic for modern turbo-charged, direct-injection gasoline engines. The general flow structure is analyzed first and the extent of the LES region and the amount of resolved fluctuations are discussed. Multiple consecutive cycles are computed and turbulent statistics of DES SST, URANS and the measured velocity field are compared for different piston positions. Cycle to cycle variations of the velocity field are analyzed for both computation and experiment with a special emphasis on the useability of the DES SST model to predict cyclic variations.     

An improved anisotropy-resolving subgrid-scale model for flows in laminar–turbulent transition region

Some types of mixed subgrid-scale (SGS) models combining an isotropic eddy-viscosity model and a scale-similarity model can be used to effectively improve the accuracy of large eddy simulation (LES) in predicting wall turbulence. Abe (2013) has recently proposed a stabilized mixed model that maintains its computational stability through a unique procedure that prevents the energy transfer between the grid-scale (GS) and SGS components induced by the scale-similarity term. At the same time, since this model can successfully predict the anisotropy of the SGS stress, the predictive performance, particularly at coarse grid resolutions, is remarkably improved in comparison with other mixed models. However, since the stabilized anisotropy-resolving SGS model includes a transport equation of the SGS turbulence energy, k_SGS, containing a production term proportional to the square root of k_SGS, its applicability to flows with both laminar and turbulent regions is not so high. This is because such a production term causes k_SGS to self-reproduce. Consequently, the laminar–turbulent transition region predicted by this model depends on the inflow or initial condition of k_SGS. To resolve these issues, in the present study, the mixed-timescale (MTS) SGS model proposed by Inagaki et al. (2005) is introduced into the stabilized mixed model as the isotropic eddy-viscosity part and the production term in the k_SGS transport equation. In the MTS model, the SGS turbulence energy, k_es, estimated by filtering the instantaneous flow field is used. Since the k_es approaches zero by itself in the laminar flow region, the self-reproduction property brought about by using the conventional k_SGS transport equation model is eliminated in this modified model. Therefore, this modification is expected to enhance the applicability of the model to flows with both laminar and turbulent regions. The model performance is tested in plane channel flows with different Reynolds numbers and in a backward-facing step flow. The results demonstrate that the proposed model successfully predicts a parabolic velocity profile under laminar flow conditions and reduces the dependence on the grid resolution to the same degree as the unmodified model by Abe (2013) for turbulent flow conditions. Moreover, it is shown that the present model is effective at transitional Reynolds numbers. Furthermore, the present model successfully provides accurate results for the backward-facing step flow with various grid resolutions. Thus, the proposed model is considered to be a refined anisotropy-resolving SGS model applicable to laminar, transitional, and turbulent flows.     

Large eddy simulation of turbulent premixed piloted flame using artificial thickened flame model coupled with tabulated chemistry

A sub-grid scale (SGS) combustion model, which combines the artificial thickened flame (ATF) model with the flamelet generated manifold (FGM) tabulation method, is proposed. Based on the analysis of laminar flame structures, two self-contained flame sensors are used to track the diffusion and reaction processes with different spatial scales in the flame front, respectively. The dynamic formulation for the proposed SGS combustion model is also performed. Large eddy simulations (LESs) of Bunsen flame F3 are used to evaluate the different SGS combustion models. The results show that the proposed SGS model has the ability in predicting the distributions of temperature and velocity reasonably, while the predictions for the distributions of some species need further improvement. The snapshots of instantaneous normalized progress variables reveal that the flame is more remarkably and severely wrinkled at the flame tip for flame F3. More satisfactory results obtained by the dynamic model indicate that it can preserve the premixed flame propagation characteristics better.     

Large Eddy Simulation of Turbulent Premixed Swirling Flames Using Dynamic Thickened Flame with Tabulated Detailed Chemistry

A sub-grid scale (SGS) combustion model by combining dynamic thickened flame (DTF) with flamelet generated manifolds (FGM) tabulation approach (i.e. DTF-FGM) is developed for investigating turbulent premixed combustion. In contrast to the thickened flame model, the dynamic thickening factor of the DTF model is determined from the flame sensor, which is obtained from the normalized gradient of the reaction progress variable from the one-dimensional freely propagating premixed flame simulations. Therewith the DTF model can ensure that the thickening of the flame is limited to the regions where it is numerically necessary. To describe the thermo-chemistry states, large eddy simulation (LES) transport equations for two characteristic scalars (the mixture fraction and the reaction progress variable) and relevant sub-grid variances in the DTF-FGM model are presented. As to the evaluation of different SGS combustion models, another model by utilizing the combination of presumed probability density function (PPDF) and FGM (i.e. PPDF-FGM) is also described. LES of two cases with or without swirl in premixed regime of the Cambridge swirl burner flames are performed to evaluate the developed SGS combustion model. The predicted results are compared with the experimental data in terms of the influence of different LES grids, model sensitivities to the thickening factor, the wrinkling factor, and the PPDF of characteristic scalars, the evaluation of different modelling approaches for the sub-grid variances of characteristic scalars, and the predictive capability of different SGS combustion models. It is shown that the LES results with the DTF-FGM model are in reasonable agreement with the experimental data, and better than the results with the PPDF-FGM approach due to its ability to predict better in regions where flame is not resolved.     

Large eddy simulation of turbulent premixed and stratified combustion using flame surface density model coupled with tabulation method

Large eddy simulations (LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density (FSD) model incorporated with two different wrinkling factor models, i.e., the Muppala and Charlette2 wrinkling factor models, is used to describe combustion/turbulence interaction, and the flamelet generated manifolds (FGM) method is employed to determine major scalars. This coupled sub-grid scale (SGS) combustion model is named as the FSD-FGM model. The FGM method can provide the detailed species in the flame which cannot be obtained from the origin FSD model. The LES results show that the FSD-FGM model has the ability of describing flame propagation, especially for stratified flames. The Charlette2 wrinkling factor model performs better than the Muppala wrinkling factor model in predicting the flame surface area change by the turbulence. The combustion characteristics are analyzed in detail by the flame index and probability distributions of the equivalence ratio and the orientation angle, which confirms that for the investigated stratified flame, the dominant combustion modes in the upstream and downstream regions are the premixed mode and the back-supported mode, respectively.     

Evaluations of SGS Combustion,Scalar Flux and Stress Models in a Turbulent Jet Premixed Flame

A newly developed fractal dynamic SGS (FDSGS) combustion model and a scale self-recognition mixed (SSRM) SGS stress model are evaluated along with other SGS combustion, scalar flux and stress models in a priori and a posteriori manners using DNS data of a hydrogen-air turbulent plane jet premixed flame. A posteriori tests reveal that the LES using the FDSGS combustion model can predict the combustion field well in terms of mean temperature distributions and peak positions in the transverse distributions of filtered reaction progress variable fluctuations. A priori and a posteriori tests of the scalar flux models show that a model proposed by Clark et al. accurately predicts the counter-gradient transport as well as the gradient diffusion, and introduction of the model of Clark et al. into the LES yields slightly better predictions of the filtered progress variable fluctuations than that of a gradient diffusion model. Evaluations of the stress models reveal that the LES with the SSRM model predicts the velocity fluctuations well compared to that with the Smagorinsky model.     

A 5-D Implementation of FGM for the Large Eddy Simulation of a Stratified Swirled Flame with Heat Loss in a Gas Turbine Combustor

Numerical simulations are foreseen to provide a tremendous increase in gas-turbine burners efficiency in the near future. Modern developments in numerical schemes, turbulence models and the consistent increase of computing power allow Large Eddy Simulation (LES) to be applied to real cold flow industrial applications. However, the detailed simulation of the gas-turbine combustion process remains still prohibited because of its enormous computational cost. Several numerical models have been developed in order to reduce the costs of flame simulations for engineering applications. In this paper, the Flamelet-Generated Manifold (FGM) chemistry reduction technique is implemented and progressively extended for the inclusion of all the combustion features that are typically observed in stationary gas-turbine combustion. These consist of stratification effects, heat loss and turbulence. Three control variables are included for the chemistry representation: the reaction evolution is described by the reaction progress variable, the heat loss is described by the enthalpy and the stratification effect is expressed by the mixture fraction. The interaction between chemistry and turbulence is considered through a presumed beta-shaped probability density function (PDF) approach, which is considered for progress variable and mixture fraction, finally attaining a 5-D manifold. The application of FGM in combination with heat loss, fuel stratification and turbulence has never been studied in literature. To this aim, a highly turbulent and swirling flame in a gas turbine combustor is computed by means of the present 5-D FGM implementation coupled to an LES turbulence model, and the results are compared with experimental data. In general, the model gives a rather good agreement with experimental data. It is shown that the inclusion of heat loss strongly enhances the temperature predictions in the whole burner and leads to greatly improved NO predictions. The use of FGM as a combustion model shows that combustion features at gas turbine conditions can be satisfactorily reproduced with a reasonable computational effort. The implemented combustion model retains most of the physical accuracy of a detailed simulation while drastically reducing its computational time, paving the way for new developments of alternative fuel usage in a cleaner and more efficient combustion.     

Flow Physics of a Bluff-Body Swirl Stabilized Flame and their Prediction by Means of a Joint Eulerian Stochastic Field and Tabulated Chemistry Approach

In the frame of this work a transported joint scalar probability density function (PDF) method is combined with the flamelet generated manifolds (FGM) tabulated chemistry approach for large eddy simulation (LES) modeling of a three-dimensional turbulent premixed swirl burner. This strategy accounts for the turbulence-chemistry interaction at reasonable computational costs. At the same time, it allows the usage of detailed chemistry mechanisms for the creation of the chemical database. The simulation results obtained are comparatively assessed along with complementary measurements. Furthermore, transient and time-averaged data are used to provide insight into the flow physics of the bluff-body swirl stabilized flame considered. The sensitivity of the results to different modeling approaches regarding the predicted flame shape and its dynamics is also investigated, where the implemented approach is compared with the well-established artificially thickened flame (ATF) combustion model. Consequently, the investigation conducted in this work aims to provide a complete picture on the ability of the proposed combustion model to reproduce the flow conditions within complex bluff-body swirl stabilized flames.     

Speckle optical measurement of a turbulent scalar field with high fluctuation amplitudes

A turbulent temperature field is produced in the gas flow downstream of a turbulence grid by passing the flow through a plane combustion front at the grid that serves as a flame holder. RMS values of the temperature fluctuations of up to 33 K are generated thereby. Properties of the turbulent scalar are measured with an optical speckle technique. The experimental results are compared with theoretical models described by Rotta (1972) and Driscoll and Kennedy (1985). From this comparison it is evident that, even at these high fluctuation amplitudes, the turbulent temperature field still behaves as a passive scalar.     

基于2D网格的轴对称浸没边界法

浸没边界法是处理颗粒两相流中运动边界问题的一种常用数值模拟方法. 当研究的物理问题的无量纲参数满足一定要求时, 该流场结构呈现轴对称状态. 为此本文提出了一种基于2D笛卡尔网格和柱坐标系的轴对称浸没边界法. 该算法采用有限体积法(FVM)对动量方程进行空间离散, 并通过阶梯状锐利界面替代真实的固体浸没边界来封闭控制方程. 为了提高计算效率, 本文采用自适应网格加密技术提高浸没边界附近网格分辨率. 由于柱坐标系的使用, 使得动量方程中的黏性项产生多余的源项, 我们对其作隐式处理. 此外, 在对小球匀速近壁运动进行直接数值模拟时, 由于球壁间隙很小, 间隙内的压力变化比较剧烈. 因此想要精确地解析流场需要很高的网格分辨率. 此时, 需要在一个时间步内多次实施投影步来保证计算的稳定性. 而在小球自由碰壁运动中, 我们通过引入一个润滑力模型使得低网格分辨率下也能模拟小球近壁处的运动. 最后通过小球和圆盘绕流、Stokes流小球近壁运动以及小球自由下落碰壁弹跳算例验证本算法对于轴对称流的静边界和动边界问题均是适用和准确的.      

A priori tests on numerical errors in large eddy simulation using finite differences and explicit filtering

When low‐order finite‐difference methods are applied in large eddy simulation (LES), the magnitude of the numerical error may be larger than that of the subgrid‐scale (SGS) term. In this paper, the effect of explicit filtering on the numerical error related to the spatial discretization of the convection term and the exact SGS term is studied a priori in the turbulent fully developed channel flow. As the filter width is increased the grid resolution is kept constant. Also filtering in the inhomogeneous wall‐normal direction is discussed. The main conclusions are related to two approaches to explicit filtering. In the traditional approach, the whole velocity field is filtered explicitly while in the alternative approach, only the non‐linear convection term of the Navier–Stokes equations is filtered explicitly. Based on the results presented in the paper it seems that the first approach leads to an unphysical situation. However, the later approach works in the desired way, and the numerical error becomes clearly smaller than the SGS term. The main difference between the two approaches seems to be the interpretation of the resolved non‐linear term in the filtered Navier–Stokes equations. Copyright © 2005 John Wiley & Sons, Ltd.     

Analysis of high-order velocity moments in a strained channel flow

In the current study, model expressions for fifth-order velocity moments obtained from the truncated Gram-Charlier series expansions model for a turbulent flow field probability density function are validated using data from direct numerical simulation (DNS) of a planar turbulent flow in a strained channel. Simplicity of the model expressions, the lack of unknown coefficients, and their applicability to non-Gaussian turbulent flows make this approach attractive to use for closing turbulent models based on the Reynolds-averaged Navier-Stokes equations. The study confirms validity of the model expressions. It also shows that the imposed flow deformation improves an agreement between the model and DNS profiles for the fifth-order moments in the flow buffer zone including when the flow reverses its direction. The study reveals sensitivity of particularly odd velocity moments to the grid resolution. A new length scale is proposed as a criterion for the grid generation near the wall and in the other flow areas dominated by high mean velocity gradients when higher-order statistics have to be collected from DNS.     

Parallel FEM LES with one-equation subgrid-scale model for incompressible flows

This article develops a parallel large-eddy simulation (LES) with a one-equation subgrid-scale (SGS) model based on the Galerkin finite element method and three-dimensional (3D) brick elements. The governing filtered Navier–Stokes equations were solved by a second-order accurate fractional-step method, which decomposed the implicit velocity–pressure coupling in incompressible flow and segregated the solution to the advection and diffusion terms. The transport equation for the SGS turbulent kinetic energy was solved to calculate the SGS processes. This FEM LES model was applied to study the turbulence of the benchmark open channel flow at a Reynolds number Re_τ = 180 (based on the friction velocity and channel height) using different model constants and grid resolutions. By comparing the turbulence statistics calculated by the current model with those obtained from direct numerical simulation (DNS) and experiments in literature, an optimum set of model constants for the current FEM LES model was established. The budgets of turbulent kinetic energy and vertical Reynolds stress were then analysed for the open channel flow. Finally, the flow structures were visualised to further reveal some important characteristics. It was demonstrated that the current model with the optimum model constants can predict well the organised structure near the wall and free surface, and can be further applied to other fundamental and engineering applications.     

The Effect of Partial Premixing and Heat Loss on the Reacting Flow Field Prediction of a Swirl Stabilized Gas Turbine Model Combustor

This work addresses the prediction of the reacting flow field in a swirl stabilized gas turbine model combustor using large-eddy simulation. The modeling of the combustion chemistry is based on laminar premixed flamelets and the effect of turbulence-chemistry interaction is considered by a presumed shape probability density function. The prediction capabilities of the presented combustion model for perfectly premixed and partially premixed conditions are demonstrated. The effect of partial premixing for the prediction of the reacting flow field is assessed by comparison of a perfectly premixed and partially premixed simulation. Even though significant mixture fraction fluctuations are observed, only small impact of the non-perfect premixing is found on the flow field and flame dynamics. Subsequently, the effect of heat loss to the walls is assessed assuming perfectly premixing. The adiabatic baseline case is compared to heat loss simulations with adiabatic and non-adiabatic chemistry tabulation. The results highlight the importance of considering the effect of heat loss on the chemical kinetics for an accurate prediction of the flow features. Both heat loss simulations significantly improve the temperature prediction, but the non-adiabatic chemistry tabulation is required to accurately capture the chemical composition in the reacting layers.     

Some Aspects of Scalar Dissipation

Scalar dissipation is of great importance in the theory and modelling of combustion and other reacting turbulent flows. Measurements of scalar dissipation are found to lack the quality assurance of checks available from the conservation equations. Conditional averages of the scalar dissipation, so important in turbulent reacting flow theory and modelling, have qualitative and quantitative dependences that are very dependent on the details of the flow and mixing conditions. Accordingly, effort needs to focus on viable means of modelling it. Fluctuations of the scalar dissipation about the conditional mean are also important. Research results in this area need to be made more accessible to the combustion scientist. Heat release effects, so important in turbulent premixed combustion, are found to be much less important in non-premixed combustion.     

On Fokker–Planck Equations for Turbulent Reacting Flows. Part 2. Filter Density Function for Large Eddy Simulation

The application of large eddy simulation (LES) to turbulent reacting flow calculations is faced with several closure problems. Suitable parametrizations for filtered reaction rates for instance are hardly available in general. A way to overcome these problems is investigated here. This is done by extending LES equations for filtered velocities and scalars (mass fractions of species and temperature) to equations that involve subgrid scale (SGS) fluctuations. Such equations are called filter density function (FDF) methods because they determine the FDF, which is essentially the probability density function of SGS variables. The FDF model considered involves only three parameters: C ₀ that controls the generation of velocity fluctuations and two parameters which determine the relaxation of velocity and scalar fluctuations. The consideration of this model may be seen as the analysis of a limiting case: the implications of the most simple equations for the dynamics of SGS fluctuations are investigated in this way. These equations were proved recently by various simulations. Here, the FDF model is used analytically to improve simpler methods. Existing models for the SGS stress tensor in velocity LES equations and the diffusion coefficient in scalar FDF equations are generalized in this way. The advantages of these models compared to existing ones are pointed out. These investigations provide further evidence for the suitability of the FDF model considered and they provide its parameters. A theoretical value C ₀ = 19/12 is derived, which agrees very well with the results of direct numerical simulation. This estimate implies the same value for the universal Kolmogorov constant of the energy spectrum, which is consistent with the results of many measurements. The other two model parameters can be obtained then by dynamic procedures. Therefore, the closure problems of LES equations are overcome in this way such that adjustable parameters are not involved.     

Large-eddy Simulation of Triangular-stabilized Lean Premixed Turbulent Flames: Quality and Error Assessment

In this numerical study, an algebraic flame surface wrinkling (AFSW) reaction submodel based on the progress variable approach is implemented in the large-eddy simulation (LES) context and validated against the triangular stabilized bluff body flame configuration measurements i.e. in VOLVO test rig. The quantitative predictability of the AFSW model is analyzed in comparison with another well validated turbulent flame speed closure (TFC) combustion model in order to help assess the behaviour of the present model and to further help improve the understanding of the flow and flame dynamics. Characterization of non-reacting (or cold) and reacting flows are performed using various subgrid scale models for consistent grid size variation with 300,000 (coarse), 1.2 million (intermediate) and 2.4 million (fine) grid cells. For non-reacting flows at inlet velocity of 17?m/s and inlet temperature 288?K, coarse grid leads to over prediction of turbulence quantities due to low dissipation at the early stage of flow development behind the bluff body that convects downstream eventually polluting the resulting solution. The simulated results with the intermediate (and fine) grid for mean flow and turbulence quantities, and the vortex shedding frequency (f_s) closely match experimental data. For combusting flows for lean propane/air mixtures at 35?m/s and 600?K, the vortex shedding frequency increase threefold compared with cold scenario. The predicted results of mean, rms velocities and reaction progress variable are generally in good agreement with experimental data. For the coarse grid the combustion predictions show a shorter recirculation region due to higher turbulent burning rate. Finally, both cold and reacting LES data are analyzed for uncertainty in the solution using two quality assessment techniques: two-grid estimator by Celik, and model and grid variation by Klein. For both approaches, the resolved turbulent kinetic energy is used to estimate the grid quality and error assessment. The quality assessment reveals that the cold flows are well resolved even on the intermediate mesh, while for the reacting flows even the fine mesh is locally not sufficient in the flamelet region. The Klein approach estimates that depending on the recirculation region in cold scenario both numerical and model errors rise near the bluff-body region, while in combusting flows these errors are significant behind the stabilizing point due to preheating of unburned mixture and reaction heat release. The total error mainly depends on the numerical error and the influence of model error is low for this configuration.     

Extension of the local domain-free discretization method to large eddy simulation of turbulent flows

In this work, an immersed boundary method, called the local domain-free discretization (DFD) method, is extended to large eddy simulation (LES) of turbulent flows. The discrete form of partial differential equations at an interior node may involve some nodes outside the solution domain. The flow variables at these exterior dependent nodes are evaluated via linear extrapolation along the direction normal to the wall. To alleviate the requirement of mesh resolution in the near-wall region, a wall model based on the turbulence boundary layer equations is introduced. The wall shear stress yielded by the wall model and the no-penetration condition are enforced at the immersed boundary to evaluate the velocity components at an exterior dependent node. For turbulence closure, a dynamic subgrid scale (SGS) model is adopted and the Lagrangian averaging procedure is used to compute the model coefficient. The SGS eddy viscosity at an exterior dependent node is set to be equal to that at the outer layer. To maintain the mass conservation near the immersed boundary, a mass source/sink term is added into the continuity equation. Numerical experiments on relatively coarse meshes with stationary or moving solid boundaries have been conducted to verify the ability of the present LES-DFD method. The predicted results agree well with the published experimental or numerical data.