
[Sponsors] 
LES  spatial, time and spectral filtering 

LinkBack  Thread Tools  Search this Thread  Display Modes 
November 2, 2018, 04:33 
LES  spatial, time and spectral filtering

#1 
New Member
Martin
Join Date: Nov 2016
Posts: 6
Rep Power: 9 
Hello,
I am reading Pierre Sagaut's book about LES for incompressible flow. I understood that in LES we solve the big turbulence scales and modelise the small's. The difference between resolved and modelised scales is fixed thanks to a filter. This filter can be spatial, spectral or time relative. My questions are the following :  What is the difference between filtering equation with a time filter and applying a time step to solve the equations ? Do we have to apply both ? Using time step is not it a way of time filtering ?  Can we superpose two filters (i.e. spatial and spectral or spectral and time) ?  I was thinking that by filtering space we implicitly filter time. Is that true ? Thank you in advance for your answers, Sincerely 

November 2, 2018, 05:26 

#2  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,738
Rep Power: 71 
Quote:
A lot of wellassessed questions on issues that are often not clear in LES  the time step defines only the cutoff filter implied by the Nyquist frequency pi/dt. This is an implicit filtering appearing when we discretize the time axis, no matter about the practical numerical method is used. Formally, the filter in time can be applied to the equations independently from the time step, for example if you want to resolve only large timescales (see the analogy to the URANS formulation) the timefilter width can be larger than the time step.  the spectral filter is also a spatial filter. You can filter in the spatial domain or in the spectral domain but you have always a correspondence between the two spaces. And yes, you can apply a filter any times you want on a variable. Just be careful that the spectral filter is idempotent, therefore you get the same original filtered field, at least if the same cutoff is used. This does not happen for filters such as the tophat.  Not exactly, you must think that a practical numerical discretization in time and space introduces implicitly a timespace filtering. That is differently induced by the chosen method of discretization. On the other hand, if you see the continuous equations, filtering in space does not correspond to a filter in time. I suggest to try develop some simple example using the 1D Burgers equation. Have a good reading 

November 2, 2018, 05:32 

#3 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,738
Rep Power: 71 
Forget to say that some of your questions are explained in my lecture notes on LES here https://www.researchgate.net/publica...CMJ6ACE7W9UnDA


November 2, 2018, 10:58 

#4 
Senior Member

But, suppressing some spatial scales wouldn't, in theory, suppress their associated time scales as well? I think Sagaut states this as well somewhere in his book...


November 2, 2018, 12:36 

#5  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,738
Rep Power: 71 
Quote:
Paolo, that could be a good argument to discuss. Here is just my opinion. Forget for now we are focusing on simulation of turbulence and let us start a general discussion to help also understanding of other people. If you discretize both space and time you introduce the cutoff frequencies pi/h and pi/dt. This is the Nyquist filter acting on any type of numerical solution (filtered or unfiltered governing equations). The relation between these two frequencies is governed by the CFL number, involving a characteristic local velocity. When CFL=O(1) it is clear the relation between the two cutoff frequencies. Now let us filter the equation using a spatial filter of width Delta. As you know, in explicit filtering Delta can be chosen independently from the computational spatial step, that is Delta=N*h. That defines the new spatial frequency pi/Delta<pi/h. The key of the discussion is if such new frequency induces a new filter in time such that pi/Delta_time<pi/dt. That would imply that a new CFL value exists to link them. Is the Delta_time = Delta/U somewhere introduced in a practical numerical LES solution? Considering now the turbulence field, this is somehow relating to the existence of a frozen velocity in the Taylor assumption. But we know that such velocity can be assumed only in specific cases and has a physical meaning. While Delta is a parameter we fix arbitrarily. Should we invoke the ergodicity hypothesis ? I am not sure about what Sagaut exactly considered in his book. Have you any idea? 

November 3, 2018, 07:02 

#6 
Senior Member

Admittedly, the demonstration on Sagaut book is not among those strongly convincing for me. Yet, I haven't the book with me right now and I can't add any detail.
But, if I picture a DNS snapshot of a velocity field and the same snapshot for the filtered DNS velocity field and a probe point where I register the time evolution, the only way the filtered field can give an unaltered history in the probe point, with respect to DNS, is by altering the bulk velocity in a very peculiar way... 

November 3, 2018, 08:14 

#7 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,738
Rep Power: 71 
My opinion is that the time filtering effect is a practical consequence of the discrete time integration that induces a smooth transfer function even when we assume that no timefiltering is introduced.


November 3, 2018, 18:21 

#8  
Senior Member

Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 13 
Quote:
If you apply a spatial filter to different snapshots of DNS solutions, there is no reason why the history of a probe point should be the same as for the DNS. Nevertheless you can state that the "change" of the history of a probe point is identical. This should hold for the smooth and discrete case. 

November 3, 2018, 18:46 

#9  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,738
Rep Power: 71 
Quote:
Any discrete time integration produces that... Just as example, if you integrate exactly in time an equation like df/dt + dF(f)/dx=0 you get f(x,t)f(x,t0)=  Int [t0,t] dF(f)/dx dt' = (tt0) <dF(f)/dx> where <*> is a timefiltering of width (tt0). But when you discretize it to solve numerically, a smoothing in the transfer function is introduced implicitly. Again, that does not require to think about turbulence or any other physics one want to solve. 

November 3, 2018, 19:06 

#10 
Senior Member

Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 13 
Dear FMDenaro,
i did not mention a time or spacial integration scheme. I only wanted to point out: No matter what integration scheme you are using, if you have the exact spacetime DNS solution and apply a spatial filter to different snapshots, i think you should compare the relative change at a given probe point and not the state. Your last explanations sound reasonable. Regards
__________________
Check out my side project: A multiphysics discontinuous Galerkin framework: Youtube, Gitlab. Last edited by Eifoehn4; November 4, 2018 at 06:01. 

November 11, 2018, 05:39 

#11 
New Member
Martin
Join Date: Nov 2016
Posts: 6
Rep Power: 9 
Thank you for the answer. FMDenaro your publication really helped me !
Sincerely, Martin 

November 11, 2018, 05:58 

#12  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,738
Rep Power: 71 
Quote:
Good, I suggest not to stop to the reading of my notes, follow the references to the publications of different authors, give yourself enough time to search and understand different ideas on the topic and, finally get your idea. There is no one exact idea and the LES field is still reach of several formulations and interpretation about the role of filtering. And if you see conflicting ideas about the role of the spatial filtering, while this is quite normal in the LES field, adding the timefiltering issue is stille somehow an obscure topic. Have a good study in LES! Regards Filippo 

November 11, 2018, 06:02 

#13  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,738
Rep Power: 71 
Quote:
The postfiltering of DNS fields is always a very good experience to do. Unfortunately, this operation get into theoretical troubles about its meaning when comparing to a real LES simulation. Are we able to filter DNS with a filter mimicking at the best the filter in the real LES? Are fields obtained by applying the static filtering on the DNS at each time correlated as happens in LES? The effect of the SGS model is dynamically different. But I always strongly suggest to make own experience in doing such exercise! 

Tags 
les, les filter, spatial, spectral, time 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Filtering DNS solutions Vs projecting filtered DNS solution  juliom  Main CFD Forum  5  May 19, 2016 17:06 
Computing the grids velocities with time spectral method.  jwang1_09  SU2  0  November 23, 2015 03:57 
Spectral LES solver, 2h waves ( wiggles )  kaya  Main CFD Forum  3  September 9, 2015 09:14 
Time Spectral Method in SU2  Rachit  SU2  2  July 14, 2015 17:24 
spectral methods for DNS and LES  ali  Main CFD Forum  1  January 16, 2006 10:47 