The question that you have is fundamental in LES, since many confusion still exists....
 implicit filtering: The governing filtered equations in LES are never acted upon an explicit filtering in the code (apart from the use of SGS dynamic procedure). Therefore, the filtering is provided implicitly by two causes, the computational grid (induces a cutoff at Nyquist frequency) and the discretization used in the code (that can change the shape of the implicit transfer function). The SGS model is some way imply the type of filtering. If you see the mathematical LES equation, you could confuse them for URANS equations! The difference is in the SGS term and in the fact that in LES you generally use high order accurate discretizations. The type of filtering (sharp cutoff, tophat, gaussian, ecc.) is due to the discretization.  Explicit filtering: you explicitly apply a filtering procedure on the nonlinear terms. This way you want to explicitly control the shape of the filtering (that you really choose) and eliminate numerical errors at high resolved wavenumbers Much more details about this issue were described in: http://www.sciencedirect.com/science...21999111000933 Good work 
Thank you very much FMDenaro!!
I am much clear now.So Basically, I just solve it like the normal NS eqns, but with the subgrid term and a higher order discretization... As for explicit filtering, I come across various discrete explicit filter equations in journals, esp. the 3  point filter like this, Now my question is: Is the way it works that the implicit LES equations are solved first, at every cell, and then the explicit discrete filter is applied to the implicit LES solution to obtain the "explicit filtered field"? Is it just like an averaging operation to obtain the explicit filtered value at a cell by taking weighted averages of the values in the cells in the neighboring region? So I assume that, when discrete explicit Gaussian filters are used, the "weights" are decided by the cutoff length of the filter and the mesh spacing at every cell. i.e. the ratio (where is the cutoff length and is the mesh spacing) is used in the definition of the Gaussian discrete explicit filter. This way, wont it make the explicit LES a two step process?  increasing computational time and cost? What is the biggest benefit of this explicit filtering, except for the fact that the LES solution is no longer dependent on the grid spacing? I dont understand the idea as to why "averaging" an already implicit filtered field would improve accuracy.It would be nice if you could shed some light on this. Once again, I thank you for your time !! P.S. : I am going through your paper, and I must say that it clears a lot of things I had questions about.Its really good. 
I'm happy to be of some help.
About explicit filtering, I can resume the two general way you can find in literature (if you see the references in my paper you will find 3  4 papers discussing explicit filtering)  Apply an explicit filtering to the implicit LES solution at each timestep. It was demonstrated by some authors (e.g. Lund) that this is not consistent, you will get a [V_bar]_bar solution that does not satisfy a filtered equation.  Apply an explicit filtering to the nonlinear flux, that is [U_bar* U_bar]_bar. This is more congruent to the governing filtered equation, it controls the numerical errors and generally you need a filter lenght from 2 to 4 times the mesh size. Of course you will get a reduced resolution. Explicit filtering requires some more computational effort. It is still controversial if use it or not ... 
That makes things more clear...
Thank you very much for your help! 
Just to add a comment, at best of my knowledge explicit filtering is only an acadamic issue, no commercial LES code (FLuent, TransAT, etc) or opensource (OpenFOAM, Code_Saturne) use that, the filter is always implicitly defined by the discretization.
That makes very often confusing to distinguish in the manuals the LES from URAN equations, they are written formally in the same way.. 
I am very interested in filtering issues for LES. As a relatively new researcher in the field, I really agree with Filippo when he says that there is much confusion about filtering; this is probably due to the fact that there is not a unique and accepted classification regarding implicit / explicit filtering techniques and so on. I still don't understand why many authors, when using LES, do not even mention which type of filtering has been used, or do not link the grid resolution to the Nyquist cutoff wavenumber and the dissipation scales.
Furthermore, there is often confusion between the mathematical model used to derive the LES equations (i.e. the convolution filter model, at foremost) and the practical LES, in which grid and discretization filters are always present. Even in a very good book like that by Pope such concepts are pretty misleading, in my opinion. OpenFOAM offers a number of dynamic models, and therefore some discrete filters are already implemented; one interesting thing which I would like to try is to exploit those builtin routines to create a solver for explicitly filtered LES. I wonder if someone has already done it... 
Quote:
you are right and it is very difficult for a newcomer in LES to focus on filtering issue... Many often the authors say simply "tophat" when they use second order FD/FV and spectral cutoff when using spectral methods... but the real shape of the transfer function can be affected by many factors... For years I tried to understand among several papers the correct way to address such an issue... I tried in my paper on JCP (cited above) to clarify what I now understand about filtering, I hope that reading it you can find some useful topic :) But I know that many things remain to be addressed ... 
Quote:
I suggest starting with a very simple test: use a velocity field coming from a well established simulation (with known spectra) and apply the discrete filters implemented in OF. Then analyze the resulting spectra in the directions to understand what the filter real does... This is the first thing to understand 
Quote:
Quote:

Quote:
I completely agree on this. The practical form of LES filtered equations and the theoretical form (convolution filter mathematical derivations) etc are not related to clearly by many authors. I'm glad Prof. Denaro cleared the confusion here. I have a question  when we apply the filtered equations, we solve for the subgrid stresses in the equation at every grid/cell, don't we? So this means that we *assume* there are scales smaller than the subgrid scale in every region (every gridpoint) of the flow. What if this is not the case in some regions? Won't the code give unphysical results in these regions? Is there a way in which the subgrid stresses can be evaluated in only in some regions of the grid, while we don't in others? Thanks for your time. 
Quote:
that happens automatically in a wellbuilt SGS model for LES.... for example, the dynamic Smagorinsky procedure leads to a vanishing eddy viscosity if the flow is laminar. You can see many papers where the eddy visosity profile is reported, for example in channel flow, and goes to zero in the viscous sublayer. This is a great difference with RANS formulation where the model acts along all scales, affecting also laminar region. 
Quote:
Thanks again. I was actually thinking of implementing the ordinary Smagorinsky model with a Cs =0.1 ( which the Fluent documentation says is a good guess to avoid excessive diffusion in large scale structures). In dynamic Smagorinsky, I know that we have to apply a test filter with a width twice that of the grid filter and the difference in the resolved scales between the two filters would aid in computation of the constant C_s. Now, Doesn't this actually look like explicit filtering? If my understanding is correct, we have to resolve the grid filter , and the apply some sort of discrete explicit filter (like the 3 point filter in my earlier post) to account for the test filter, and use the resolved fields obtained from both the filters, to calculate the smagorinsky constant using the Germano identity. is this correct? or is method of applying the test filter different ? Thanks. 
Quote:
The explicit filtering in LES involves the application of a filter on the convective flux in the momentum equation whilst (as you stated) the dynamic procedure is based on the application of a testfilter on all the terms of the equation, and the Germano identity is derived for the SGS tensors. Both procedures require an algorithm for doing an explicit filtering. What kind of discrete filtering works the best is a further complex issue.... 
Quote:
Its more clear now.....I guess I will implement the ordinary Smagorinsky model as of now, since this is going to be my first LES code, and then move on to the dynamic model if required. Thank you for you help all along, Prof Denaro! 
Quote:
Now, my doubt is: which part of a (let's say dynamic) SGS model is responsible for this? Is it the dynamic calculation of the constant? I guess the standard Smagorinsky is not capable of doing that, as the constant is fixed and is calculated assuming, a priori, a certain position of the filter width along the spectrum, is that right? 
Quote:
Have a good work 
Quote:
well, you can observe that the filtered velocity can be written in the Fourier space as u_bar(k)=G(k)*u(k). When apply a test filtering (its transfer function I call T), you have [u_bar]_test(k)=T(k)*G(k)*u(k). In laminar flows you have only very low wavenumber components where T and G become practically =1. So, the testfiltered velocity is equal to the gridfiltered velocity and the SGS term will vanish ... see the Germano identity and modelling. For the static Smagorinsky model, when the flow has some direction in which is not homogeneous (for example in channel flow), it is advised to use the Van Driest damping function that depends on the local position. Note that convergence of an LES solution towards DNS for refined grids is true only if you do not use explicit filtering. 
Ok, clear, thank you.
I was also trying to better understand the relationship between the filter width and the model constant. Considering a static Smagorinsky, I guess we cannot use the same constant if the filter width is in the dissipation range (i.e. very fine grid) or if it is in the inertial subrange. For example, if I remember well, the Lilly's constant for the Smagorinsky model is calculated assuming the filter width in the inertial subrange. However, this is maybe unuseful, because for small filter widths the viscous stress should be much higher than the residual stress. 
Quote:
Delta is someway related to the computational grid measure, therefore for non uniform grids also Delta=Delta(x). Therefore (Cs*Delta)^2 depends on the position. You can see the great advantage of the dynamic procedure that allows us to compute locally all the term (Cs*Delta)^2. The static Smagorinsky requires some ad hoc adjustement. As I wrote before, the example of the channel is good. Assume a grid very refined near the wall, so that you have some points before y+=1. This way practically the filter size is at DNS resolution in the boundary layer. The static Smagorinsky viscosity must be damped by using the Van Driest dumping. Conversely, the dynamic Smagorinsky works automatically! ;) 
All times are GMT 4. The time now is 12:59. 