Dynamic nonlinear algebraic models with scale-similarity dynamic procedure for large-eddy simulation of turbulence

A dynamic nonlinear algebraic model with scale-similarity dynamic procedure (DNAM-SSD) is proposed for subgrid-scale (SGS) stress in large-eddy simulation of turbulence. The model coefficients of the DNAM-SSD model are adaptively calculated through the scale-similarity relation, which greatly simplifies the conventional Germano-identity based dynamic procedure (GID). The a priori study shows that the DNAM-SSD model predicts the SGS stress considerably better than the conventional velocity gradient model (VGM), dynamic Smagorinsky model (DSM), dynamic mixed model (DMM) and DNAM-GID model at a variety of filter widths ranging from inertial to viscous ranges. The correlation coefficients of the SGS stress predicted by the DNAM-SSD model can be larger than 95% with the relative errors lower than 30%. In the a posteriori testings of LES, the DNAM-SSD model outperforms the implicit LES (ILES), DSM, DMM and DNAM-GID models without increasing computational costs, which only takes up half the time of the DNAM-GID model. The DNAM-SSD model accurately predicts plenty of turbulent statistics and instantaneous spatial structures in reasonable agreement with the filtered DNS data. These results indicate that the current DNAM-SSD model is attractive for the development of highly accurate SGS models for LES of turbulence.

whose model coefficient for the original version is statically adjusted by the experimental and DNS data in the early stage. Germano et al. [26] and Lilly [27] pioneered the development of a dynamical procedure based on the Germano identity through the leastsquares algorithm, which makes the parameter of the dynamic Smagorinsky model (DSM) dynamically determined as the flow changes. Subsequently, the dynamical versions of some conventional SGS models with the Germano-identity based dynamic procedure (GID) were proposed [4][5][6][7], including the dynamic mixed model (DMM) [28][29][30][31], the dynamic Clark model [32], the dynamic localization model [33], etc.
The Smagorinsky model [8,9,26,27] constructs the SGS stress with the linear constitutive relation based on the Boussinesq hypothesis, which requires the alignment between the SGS stress and the filtered strain-rate tensor. Pope [34] derived the general expression between the Reynolds stress and the averaged strain-rate and rotation-rate tensors with eleven integrity basis tensors based on the theory of invariants. Due to the expensive calculations of the high-order basis tensors in the general expression, the numerical verification of Pope's general viscous hypothesis was only limited to the two-dimensional turbulence [34]. Lund and Novikov [35] showed that the sixth invariant can be expressed as the ratio of the other five invariants, and reduced the original eleven polynomial basis tensors to five, which greatly simplified the computational complexity of the nonlinear algebraic SGS model in LES calculations. Especially, the anisotropic part of the SGS stress can be expressed as the general expression of the resolved strain-rate and rotationrate tensors with five model coefficients [35]. Speziale et al. [36,37] further simplified Lund's general expression to a quadratic constitutive relation for the Reynolds stress. The model coefficients of the nonlinear algebraic model were mostly determined by the DNS data in the early research work. Wong [38] proposed a two-parameter dynamic nonlinear algebraic model (DNAM) using the quadratic constitutive relation with the Germano-identity based dynamic procedure. Kosović [39] applied the nonlinear constitutive relation to the shear-driven boundary layers at high Reynolds number. Wang et al. [40,41] proposed a dynamic SGS model based on the quadratic nonlinear constitutive expression with local stability. Marstorp et al. [42] proposed an explicit algebraic SGS stress model with the equilibrium assumptions made on the partial-differential equations of SGS stress, and successfully applied to the rotational channel flow. Recently, a stochastic extension of the explicit algebraic SGS models has been developed by Rasam et al. [43] In our previous research work, a nonlinear algebraic model based on the artificial neural network (ANN-NAM) was proposed [44], whose model coefficients are predicted by the invariant-input ANN with embedded invariance. The ANN-NAM model [44] reconstructs the SGS stress and statistics of velocity with high accuracy both in the a priori and a posteriori analyses of LES. Wang et al. [45] proposed an ANN-based semi-explicit spatial gradient model with embedded invariance. A dynamic version of the spatial gradient model (DSGM) [46] was proposed for the parameter determination strategy. Yuan et al. developed deconvolutional artificial neural network (DANN) [21] and dynamic iterative approximate deconvolution (DIAD) models [22] to recover the local unfiltered velocity with the neighboring spatial stencils of the filtered velocity. A scale-similarity-based dynamic procedure (SSD) was proposed to adaptively calculate the weights of the spatial stencil [22]. The DIAD model with the SSD procedure is superior to the other conventional dynamic SGS models in the reconstruction of the statistics of velocity and transient coherent structures of turbulence [22].
The SSD procedure is inspired by the scale-similarity hypothesis [12], assuming that the original unfiltered velocity shares the consistent constitutive relation with the resolved filtered velocity in the inertial region. Bardina et al. [12] developed the scale-similarity model which is formally consistent with the Leonard stress derived by the Germano decomposition [26]. Liu et al. [13] extended the scale-similarity model by introducing the test-level filter. He et al. [47] developed a universal form of Lagrangian velocity correlations at different space separations based on the scale-similarity hypothesis. Stallcup et al. [48,49] proposed an adaptive scale-similar closure using the generalized representations of the SGS terms by solving a local system identification problem at the test-filter scale. The adaptive scale-similar model can accurately represent the SGS terms near the smallest resolved scales with the minimal tensor representation [48,49].
In the current research, a novel dynamic nonlinear algebraic model with scale-similarity dynamic procedure (DNAM-SSD) is developed for reconstructing the unclosed SGS stress in LES of incompressible turbulence. The performance of the DNAM-SSD model is examined by comparing with those of some classical SGS models both in the a priori and a posteriori testings of LES at two filter widths 16h DNS and 32h DNS with the corresponding grid resolutions of N 128 3 and 64 3 . The computational accuracy and costs of the newly-proposed scale-similarity based dynamic (SSD) procedure are compared to the conventional Germano-identity based dynamic (GID) procedure. The remainder of the paper is organized as follows. The governing equations of LES will be described in Section 2. The introductions of the conventional SGS models and DNAM models are respectively illustrated in Sections 3 and 4. Section 5 will conduct the numerical simulation of incompressible isotropic turbulence. The a priori and a posteriori studies are correspondingly provided in Sections 6 and 7. Conclusions are drawn in Section 8.

Governing equations of large-eddy simulation
The incompressible turbulence is governed by the Navier-Stokes equations, whose dimensionless conservation form is written as [1] u i where u i denotes the i-th velocity component (i 1, 2, 3 represents the three directions of the Cartesian coordinate system, respectively.), p is the pressure, Re is the Reynolds number, and i is the i-th large-scale force component. [21,50,51] For brevity and simplicity, we adopt the summation convection for the repeated indices by default in this paper. Besides, the governing dimensionless parameter for the incompressible turbulence, namely, the Taylor microscale Reynolds number Re is given by [1] Re u rms 3 , where denotes the kinematic viscosity and u rms u i u i is the root-mean-square (rms) value of the velocity magnitude. Here, " " represents a spatial average over the entire computational domain. In addition, the Taylor microscale is expressed as [1] u rms 5 , where 2 S ij S ij denotes the dissipation rate and S ij 1 2 u i x j u j x i is the strain-rate tensor.
For the large-eddy simulation, the resolved large scales are separated from the subgrid small scales by the spatial filtering operation, which is introduced as [ where f x represents the arbitrary physical variable, and an overbar stands for the low-pass spatial filtering. Here, denotes the entire physical domain, with G and respectively being the spatial filter function and filter width. The governing equations for the LES can be obtained by applying the spatial filtering on the Eqs. (1) and (2), correspondingly, which can be derived as [1,2] u i x i 0, Here, the unclosed SGS stress ij in the Eq. (7) is defined by [4][5][6] ij u i u j u i u j .
The SGS stress involves the nonlinear interactions between the resolved large scales and under-solved small scales, therefore additional SGS stress modeling is required to close the governing equations of LES. In the following two sections, the conventional SGS models and the proposed dynamic nonlinear algebraic models with scale-similarity dynamic procedure (SSD) are respectively described for the LES computations.

Conventional SGS models
The explicit modeling for the unclosed SGS stress can be divided into the functional modeling and structural modeling. The functional models mimic the forward energy transfer from the resolved large scales to the residual small scales by constructing the explicit eddy-viscosity forms, while the structural modeling is established by the truncated series expansions or the hypothesis of scale similarity to correctly recover the SGS stress with high accuracy. A typical functional model is the dynamic Smagorinsky model (DSM), whose constitutive relation for the deviatoric SGS stress is given by [26,27] where S 2S ij S ij 1 2 is the characteristic filtered strain rate, and S ij 1 2 u i x j u j x i is the filtered strain-rate tensor. The superscript "A" represents the trace-free part of the arbitrary variables, namely, A ij ij kk ij 3. Here, the isotropic SGS stress kk is absorbed into the pressure term. C 2 S is the Smagorinsky coefficient, which can be determined by the Germano identity dynamic procedure (GID). The test-filter level SGS stress with the double-filtering scale 2 is expressed as [26,27] ij u i u j u i u j , where a tilde stands for the test filtering operation at the filter scale . The deviatoric part of ij can be modeled based on the scale-invariance hypothesis, defined by [26,27] These two SGS stresses with different filter scales, namely, ij and ij satisfy the Germano identity, expressed as [26] ij ij ij where the Leonard stress ij can be calculated using the resolved filtered field for LES calculations. Therefore, the optimal Smagorinsky coefficient C 2 S can be further determined by the least-squares algorithm, namely [27] where A ij ij 1 3 ij kk , and ij ij ij . Here ij 2 2 S S ij , ij 2 2 S S ij . A typical structural model is the velocity gradient model (VGM) based on the truncated Taylor series expansions, given by [11] ij 2 12 The dynamic mixed model (DMM) combines the scale-similarity model with the dissipative Smagorinsky term, which can overcome the deficiency of numerical instability in the structural modeling of the SGS stress. The SGS stresses constructed by the DMM model at scales and are expressed, respectively, as [12,28,52] ij where h 1,ij 2 2 S S ij , h 2,ij u i u j u i u j , H 1,ij 2 2 S S ij , and H 2,ij u i u j u i u j .
Here, the hat stands for the test filtering at scale 4 . Similar to the DSM model, the model coefficients C 1 and C 2 are calculated by the Germano identity dynamic procedure, namely [21,22] where M ij H 1,ij h 1,ij , and N ij H 2,ij h 2,ij .

Dynamic nonlinear algebraic models (DNAM)
In the SGS stress modeling, the constitutive relation of the unclosed SGS stress can be regarded as the function of the local filtered quantities, i.e., the filtered strain-rate tensor S ij and filtered rotation-rate tensor ij , namely [34,35] where the filtered rotation-rate tensor ij 1 2 u i x j u j x i . For brevity and simplicity of the tensorial polynomials, the matrix multiplications for the tensor contractions are expressed as [34,35,40] A general expression of the modeled SGS stress [Eq. (19)] can be expanded to the sum of an infinite number of tensorial polynomials with the form S m 1 n 1 S m 2 n 2 , where m i and n i are positive integers. The infinite tensorial polynomials can be reduced to a finite number by the Cayley-Hamilton theorem [34,35,40], thus the modeled SGS stress can be expressed as the linear combination of the basis tensors formed by the product of S and , namely [34] ij 11 n 1 g n 1 , 2 , where T n ij is the n-th basis tensor and the model coefficients g n are functions of the six integrity invariants m (m 1, 2, , 6). Here, the eleven basis tensors T n ij and six independent invariants m are respectively expressed as [34] If the model coefficients g n are relaxed as the ratios of polynomials of these integrity invariants, the number of the above basis tensors can be reduced from eleven to five. In accordance with the dimensional consistency, the anisotropic part of the modeled SGS stress can be given by [35] A ij 2 where the characteristic filtered strain rate S 2S ij S ij 1 2 , and C n are five dimensionless model coefficients. The corresponding basis tensors n ij that satisfy the consistent dimension with the square of the velocity gradient are defined by In the paper, two dynamic procedures are adopted to determine the model coefficients C n of the dynamic nonlinear algebraic models (DNAM). One is the Germano identity dynamic (GID) procedure based on the scale-invariance assumption, and the other is the newly proposed scale-similarity dynamic (SSD) procedure in accordance with the scale-similarity relation. The rest of this section will be divided into two subsections to respectively introduce these two different modeling approaches.

DNAM models with Germano identity dynamic procedure (DNAM-GID)
Similar to the conventional dynamic SGS models (e.g. DSM and DMM models), the DNAM model with Germano identity dynamic procedure, abbreviated as DNAM-GID, introduces the test-filter level SGS stress ij at the double-filtering scale 2 , modeled by Here, n ij is the n-th basis tensor at the test filter scale 2 , expressed as Consistent with Eq. (12), the modeled SGS stresses ij and ij at different filter scales satisfy the Germano identity, namely where n ij The model coefficients C n can be further calculated by the least-squares algorithm, derived by For the DNAM-GID model, the optimal model coefficients C n can be obtained by solving the system of five linear equations in Eq. (29).

DNAM models with scale-similarity dynamic procedure (DNAM-SSD)
In this paper, we propose a novel scale-similarity dynamic procedure for the DNAM model to determine the optimal model coefficients dynamically. The real SGS stress can be regarded as the nonlinear function of the velocity u i and the filter kernel at scale , whereas the SGS stress modeled by the DNAM model has the nonlinear constitutive relation with the local filtered physical quantities (e.g. the filtered strain-rate and rotation-rate tensors S and ), Based on the scale-similarity hypothesis, the modeled SGS stress at the filter scale shares the consistent model coefficients C n with that at the filter scale , namely The constitutive equation of the SGS stress is assumed to be invariant to the physical field, therefore we can replace the unfiltered velocity u i with the filtered velocity u i in Eq. (31) and obtain It is worth noting that the DNAM-SSD model only calculates in the DNAM-GID model, therefore the scale-similarity dynamic procedure simplifies the conventional dynamic procedure based on the Germano identity. Besides, in the following sections, we can show that the DNAM model with the proposed scale-similarity dynamic procedure performs better than that with the conventional GID procedure both in the a priori and the a posteriori testings of incompressible turbulence.

Numerical simulation of incompressible isotropic turbulence
In order to validate the performance of the proposed DNAM-SSD model, the numerical simulation of incompressible isotropic turbulence is performed in a cubic box of 2 3 with periodic boundary conditions at the Taylor Reynolds number Re 250.
The pseudo-spectral approach with the two-thirds dealiasing rule is adopted for the spatial discretization of the governing equation. A second-order explicit Adams-Bashforth scheme [53] is applied to the time advancement. The large-scale forcing is implemented on the two lowest wavenumber shells [21,50,51] to keep the turbulence in equilibrium.
In the paper, we use N 1024 3 uniform grids in the DNS calculation with the grid spacing h DNS 2 1024. The kinematic viscosity is set to 1 Re 0.001. The detailed one-point statistics of DNS calculation are summarized in Table 1. Here, k max 2 3h DNS represents the largest effective wavenumber after the fully two-thirds dealiasing, and rms i i stands for the root-mean-square value of the vorticity magnitude , where u represents the vorticity which is the curl of the velocity field. The Kolmogorov length scale and the integral length scale L I represent the smallest resolved scale and the largest characteristic scale, which are defined, respectively, by 3 1 4 , where 2 S ij S ij is the dissipation rate. The total turbulent kinetic energy E k 1 2 u i u i 0 E k dk, and E k stands for the velocity spectrum. The resolution parameter k max 2.1 is found to be sufficient enough for the convergence of turbulent kinetic energy at all scales [54,55].
In the paper, the filtered physical quantities and the real SGS stress ij are calculated using a Gaussian filter, which is expressed as [1,2] We select two filter scales ( 16h DNS and 32h DNS ) for model verification, and the corresponding cutoff wavenumbers are k c 32 and 16, respectively. Figure 1 shows the velocity spectra of the DNS and filtered DNS at both filter widths ( 16h DNS and 32h DNS ). The filtered velocity spectra almost overlap with the DNS data at the lowwavenumber region satisfying the Kolmogorov k 5 3 scaling, while generally diminish with the increasing of wavenumbers, and drop rapidly at the region larger than the truncated wavenumber k c . More kinetic energy is filtered out at a larger filter scale, therefore the filtered velocity spectrum at 32h DNS is lower than that at 16h DNS . Overall 95% and 88% of the turbulent kinetic energy is retained in the filtered velocity field at the filter widths 16h DNS and 32h DNS , respectively.

A priori study of the DNAM models
In the a priori analysis, twenty snapshots of DNS data at equal temporal intervals during two large-eddy turnover periods ( L I u rms ) are adopted to examine the model accuracy of the DNAM-GID and DNAM-SSD models with several filter scales ranging from 4h DNS to 64h DNS . Two evaluation metrics are used to quantify the distinction between the real value (Q real ) and the modeled value (Q model ) for targeted variable Q, namely the correlation coefficient C Q and the relative error E r Q , respectively defined by [21,22] where represents the ensemble average of total samples. In the a priori study, we first investigate the impact of integrity basis tensors n ,A ij on the SGS stress ij by calculating the correlation coefficients at different filter widths, shown in S 2 A , therefore we keep all five basis tensors without any simplification in the paper. With the increasing of the filter widths, the correlation coefficients between the first basis tensor 1 ,A ij and the SGS stress constantly increase, while those of the other four terms gradually drop but are still higher than those of the first term. These results indicate that the classical Smagorinsky model (linear relation with only the first basis tensor 1 ,A ij ) cannot fully reconstruct the SGS stress. Figures 2 and 3 respectively illustrate the correlation coefficients and relative errors of the normal and shear components of the SGS stress for different SGS models at a number of filter scales ranging from the inertial region to the dissipation range.
Here, the VGM model is the velocity gradient model (see Eq. (14)) which has a high a priori accuracy among the classical SGS models. The DNAM-LS model is a DNAM model with a priori knowledge of DNS data, whose model coefficients are calculated by the least-squares method using the real SGS stress, namely  29)) and the newly proposed scale-similarity dynamic procedure (cf. Eq. (33)), respectively.  The DNAM-LS model has the highest correlation coefficients and the lowest relative errors with the SGS stress, since the DNS data are used to determine the model coefficients. The correlation coefficients and relative errors predicted by the proposed DNAM-SSD model are very close to the DNAM-LS model at all filter scales, which are much better than the DNAM-GID and VGM models. The DNAM-SSD model predicts the SGS stress accurately with the correlation coefficients overall higher than 92% and the relative errors less than 40% ranging from the viscous region to the inertial region. In contrast, the DNAM-GID model gives the worst prediction among these SGS models with the relative errors approximately over 40%. It is worth noting that the DNAM-SSD model performs better than the conventional VGM model at all filter widths, indicating that the basis tensors of the DNAM model are more complete than the velocity gradient in reconstructing the SGS stress.
In order to further quantify the model accuracy of different SGS models in the a priori analysis, we compare the correlation coefficients and relative errors of the SGS stress at  Tables 3 and 4, respectively. The DSM and DMM models use the Germano-identity dynamic procedure (see Eqs. (13), (17) and (18)) to dynamically determine model coefficients. The DNAM-SSD and DNAM-LS models give the best prediction of the SGS stress with correlation coefficients higher than 95% and 92% as well as relative errors lower than 30% and 38% at corresponding filter widths 16h DNS and 32h DNS among these SGS models. In contrast, the DSM model performs the worst compared to other SGS models at both filter scales, whose correlation coefficients are lower than 30% and relative errors are nearly 100%. The DNAM-GID model predicts the SGS stress tangibly worse than the DNAM-SSD and DNAM-LS models, but is still much better than the classical DMM model with the consistent GID dynamic procedure. Besides, the performance of the VGM model in the a priori study is between the DNAM-GID model and DNAM-SSD model at both filter widths. We also compare the DNAM-SSD model with the artificial neural network-based spatial gradient models (ANN-SGM) proposed by Wang et al. [45] The ANNSGM-7-49 model consists of four fully-connected layers of neurons (4:20:20:49) and takes the integrity invariants as input to learn the model coefficients of the velocity gradient products for the neighboring seven-point stencil [45]. The DNAM-SSD model can accurately reconstruct the SGS stress with the similar accuracy to the ANN-based SGS model (ANNSGM-7-49) in the a priori analysis. These results demonstrate that the basis tensors of the DNAM model are more complete in modeling the SGS stress compared to those of the DMM model and VGM model. In the a priori analysis, the proposed scale-similarity dynamic procedure (SSD) shows distinct advantages over the conventional Germano-identity dynamic procedure (GID) in determining the model coefficients of SGS models.
Finally, we evaluate the SGS energy transfer for different SGS models by comparing the normalized SGS energy flux DNS shown in Fig. 4, where ij S ij represents the SGS energy flux and DNS denotes the dissipation rate calculated using the DNS data. The PDFs of the SGS energy flux reconstructed by the DNAM-SSD and DNAM-LS models coincide with the filtered DNS data at both filter scales 16h DNS and 32h DNS , which are obviously better than the VGM and DNAM-GID models. In comparison, the DNAM-GID model fails to predict the SGS energy transfer, indicating that the model-coefficient determination by the proposed SSD procedure is superior to that of the conventional GID procedure, and is well approximated to the DNAM-LS model in the a priori study.

A posteriori study of the DNAM models
The a posteriori testing of LES is important to illustrate the practical performance of the SGS models. In LES computations, the kinematic viscosity 0.001 is consistent with that of DNS. It is worth noting that an explicit filtering operation at the test-level filter scale is only introduced to determine the model coefficients, and no additional explicit filtering at the filter width is performed on the primary variables u i in the computations of LES. Two filter scales 16h DNS and 32h DNS are selected to study the impact of filter widths on the SGS modeling. The newly proposed DNAM-SSD model is compared to the classical SGS models, including the implicit LES (ILES), the dynamic Smagorinsky model (DSM), the dynamic mixed model (DMM), the dynamic nonlinear algebraic model with conventional Germano-identity dynamic procedure (DNAM-GID) and the ANN-based SGS model (ANNSGM-7-49) [45]. It has been found that the filter-to-grid ratio FGR h LES 2 can effectively reduce the influence of the spatial discretization errors on the SGS stress modeling [56][57][58]. Therefore, we fix the FGR value to 2 and the corresponding grid points of LES are N 128 3 and 64 3 for the selected filter widths 16h DNS and 32h DNS . It is worth noting that both filter scales lie in the inertial range (cf. Fig. 1) and the scale-invariance assumption still holds, which are essential for the conventional dynamic models based on the Germano identity. The time step of LES is selected as t LES 10 t DNS for all SGS models. The coarse-grained LES computations without any SGS models suffer from the numerical instability due to the nonlinear interactions between resolved large scales and residual scales. Implicit LES methods adopt the artificial dissipation to mimic the forward kinetic energy transfer from large scales to small scales. In this paper, we use the six-order compact-difference filtering scheme to provide necessary artificial dissipation for the numerical stability of coarse-grained LES computations, expressed as [59,60] where a check " " represents the dissipative filtering operation applied to the velocity field after each calculation of ILES, the parameter 0 f 0.5, and four coefficients a n are determined by the Taylor-series expansion, namely [59,60] In the paper, the control parameter f 0.495 is chosen for a commonly-used sixorder low-dissipation filtering scheme [59,60]. The pure scale-similarity type SGS model (DNAM-SSD) is numerically unstable due to lack of sufficient dissipation. Therefore, we couple the consistent compact-difference dissipative filtering scheme with the DNAM-SSD model to ensure the numerical stability of the LES calculations. For the ANNSGM-7-49 model, an artificial dissipation with the fourth-order hyperviscosity is introduced to maintain the numerical stability of the LES calculations [45]. The same instantaneous snapshot of the filtered DNS data is used as the initialization of the LES calculations for different SGS models. Table 5 summarizes the average computational time for the SGS stress modeling at both filter scales 16h DNS and 32h DNS . Compared to the classical DSM and DMM models, the DNAM model would not particularly increase too much computational cost: the modeling time of the DNAM-GID model is about 1.1 times that of the DMM model. The proposed SSD procedure simplifies the calculation process of the conventional GID method, which greatly reduces the computation costs and only takes about half the time of the DNAM-GID model (0.6 times that of the DMM model) at both filter widths 16h DNS and 32h DNS . In the a posteriori testings of LES computations, we first compare the velocity spectra of different SGS models with those of the DNS and filtered DNS (fDNS) data at both filter scales 16h DNS and 32h DNS shown in Figs. 5 and 6. LES of the incompressible turbulence is governed by the filtered Navier-Stokes equations (Eqs. (6) and (7)), therefore, the statistics of an ideal LES would be close to that of the fDNS data. The velocity spectra predicted by the ILES model are insufficiently dissipated and obviously higher than the fDNS data, while those modeled by the DSM and DMM models exhibit the tilde distribution due to excessive dissipation. The ILES generally reconstructs the inter-scale interactions and subgrid-scale effects by the artificial dissipation or dissipative spatial discretization schemes. Therefore, the ILES results are generally obtained without any explicit filtering operations [1][2][3]. Moreover, there is no filter width in ILES, and ILES results are very sensitive to the change of grid number. Statistics of ILES could converge to those of DNS when the grid is refined to the DNS level. For the DSM and DMM models, velocity spectra near the truncated wavenumbers are diminished by the model dissipation, resulting in the blockage of the kinetic energy cascade from large scales to small scales. Therefore the kinetic energy accumulates in the region of medium wavenumbers. The predicted velocity spectra of the DNAM-GID model are very similar to the DMM model, which has a significant deviation from the fDNS data. In contrast, the DNAM-SSD model outperforms the other SGS models and accurately reconstructs the velocity spectra at both filter scales 16h DNS and 32h DNS , which is very close to the ANNSGM-7-49 model. It is  Fig. 7. In order to guarantee the numerical stability of the LES with the DSM model, the model coefficient C S is restricted to be non-negative. The DSM model only reconstructs the forward SGS energy transfer whose modeled SGS energy flux is always non-negative. The PDFs of the SGS energy flux for the DMM and DNAM-GID models are obviously narrower than the fDNS data. In comparison, the DNAM-SSD model can accurately mimic both the forward SGS energy transfer and the energy backscatter, which is very similar to the ANN-based model (ANNSGM-7-49).
We further compare the normalized strain-rate tensor reconstructed by different SGS models, whose PDFs of the normal and shear components at both filter widths 16h DNS and 32h DNS are illustrated in Figs. 8 and 9, respectively. In the figures, S rms ij,fDNS S fDNS ij 2 represents the root-mean-square value of the strain-rate tensor calculated using the filtered DNS data shown in Table 6. The strain-rate tensor represents the local  stands for the root-mean-square value of the characteristic strain rate given from the filtered DNS data (see Table 6). ILES obviously fails to predict the characteristic strain rate, while the DSM, DMM and DNAM-GID models cannot well capture the peak of the PDFs in reconstructing the characteristic strain rate. In comparison, the PDFs of the characteristic strain rate predicted by the DNAM-SSD model are very close to the fDNS data. We finally examine the reconstruction ability of the turbulent coherent structure by comparing the transient contours of the normalized vorticity, shown in Fig. 11. Here, the vorticity magnitude is normalized by the root-mean-square values calculated using the filtered DNS data: rms fDNS fDNS i fDNS i (see Table 6), where u denotes the resolved vorticity which is the curl of the resolved velocity field. The snapshots of LES calculations for different SGS models are selected on an arbitrary XY slice at the consistent Fig. 9 PDFs of the normalized strain-rate tensor in a posteriori study at filter scale 32h DNS time with approximately two large-eddy turnover periods. For the DSM and DNAM-GID models, some small-scale flow structures are excessively dissipated and only large scales are maintained. Compared to the other SGS models, the vortex structures reconstructed by the DNAM-SSD model exhibit more similar spatial distribution to the fDNS data, and more multiple-scale flow structures are accurately recovered by the proposed DNAM-SSD model. In order to test the impact of the explicit filters on the accuracy of the DNAM models, we choose two different types of explicit filters (the top-hat and differential Helmholtz filters). The top-hat filter in one dimension is expressed as [57,61,62] where n h DNS . The explicit form of differential Helmholtz filter is expressed as [63] u i 1 2 2 1 u i 1 2 2 x k x k where 2 2 24 [63]. The velocity spectra for different SGS models (ILES, DSM, DMM, DNAM-GID and DNAM-SSD) at the filter width 16h DNS with the top-hat and Helmholtz filters are displayed in Figs. 12 and 13, respectively. The LES results using both the top-hat and Helmholtz filters are very similar to those with the Gaussian filter (see Fig. 5). The ILES model exhibits insufficient dissipation, while the DSM, DMM and DNAM-GID models show bump distributions with excessive dissipation. In comparison, the DNAM-SSD model accurately predicts the velocity spectrum, which is almost coincident with the fDNS data. These results demonstrate that the accuracy of the DNAM models is not significantly affected by the type of the explicit filters.

Conclusions
In the current work, we develop a dynamic nonlinear algebraic model with the newly proposed scale-similarity dynamic procedure (DNAM-SSD) for the large-eddy simulation of turbulence. In the DNAM-SSD model, the model coefficients are dynamically determined based on the scale-similarity relation, which greatly simplifies the conventional dynamic procedure based on the Germano identity (GID). The a priori analysis demonstrates that the proposed DNAM-SSD model outperforms the conventional velocity gradient model (VGM) and DNAM-GID model at a number of filter scales ranging from the inertial to dissipation ranges. The DNAM-SSD model gives the best prediction of the SGS stress with correlation coefficients higher than 95% and 92% as well as relative errors lower than 30% and 38% at corresponding filter scales 16h DNS and 32h DNS in comparison with the dynamic Smagorinsky model (DSM), dynamic mixed model (DMM), VGM model and the DNAM-GID model, respectively. The proposed SSD procedure shows significant advantages over the conventional GID approach in determining the model coefficients of SGS models.
In the a posteriori testings of LES, the performance of the proposed DNAM-SSD model is examined at both filter widths 16h DNS and 32h DNS . The classical implicit-LES (ILES), DSM, DMM and DNAM-GID models are used for comparisons of the a posteriori model accuracy. ILES fails to predict the statistics of turbulence with insufficient dissipation, while the DSM and DMM models are over-dissipative, leading to the fact that small scales are diminished by the excessive dissipation. The results predicted by the DNAM-GID model are very similar to those of the DMM model and have obvious deviations from the filtered DNS data. In contrast, the predictions of the DNAM-SSD model are very close to the filtered DNS data in the velocity spectra, the statistics of SGS energy flux and strain rate, as well as the instantaneous spatial structures of the vorticity magnitude at both filter scales without increasing the computational cost. The modeling time of the DNAM-SSD model is only half the time of the DNAM-GID model (0.6 times that of the DMM model) at both filter widths. These results demonstrate that the current DNAM-SSD model is an effective framework for enhancing the advanced SGS stress modeling of LES.