On the properties of discrete spatial filters for CFD

The spatial filtering of variables in the context of Computational Fluid Dynamics (CFD) is a common practice. Most of the discrete filters used in CFD simulations are locally accurate models of continuous operators. However, when filters are adaptative, i.e. the filter width is not constant, or meshes are irregular, discrete filters sometimes break relevant global properties of the continuous models they are based on. For example, the principle of maxima and minima reduction or conservation are eventually infringed. In this paper, we analyze the properties of analytic continuous convolution filters and extract those we consider to define filtering. Then, we impose the accomplishment of these properties on explicit discrete filters by means of constraints. Three filters satisfying the derived conditions are deduced and compared to common differential discrete CFD filters on synthetic fields. Tests on the developed discrete filters show the fulfillment of the imposed properties. In particular, the problem of maxima and minima generation is resolved for physically relevant cases. The tests are conducted on the basis of the eigenvectors of graph Laplacian matrices of meshes. Thus, insight into the relations between filtering and oscillation growth on general meshes is provided. Further tests on singularity fields and on isentropic vortices have also been conducted to evaluate the performance of filters on basic CFD fields. Results confirm that imposing the proposed conditions makes discrete filters properties consistent with those of the continuous ones.

[1]  M. Germano,et al.  Turbulence: the filtering approach , 1992, Journal of Fluid Mechanics.

[2]  Arthur Veldman,et al.  A symmetry-preserving discretisation and regularisation model for compressible flow with application to turbulent channel flow , 2014 .

[3]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[4]  P. Moin,et al.  A dynamic subgrid‐scale eddy viscosity model , 1990 .

[5]  O. Vasilyev,et al.  Commutative discrete filtering on unstructured grids based on least-squares techniques , 2001 .

[6]  F. X. Trias,et al.  On the construction of discrete filters for symmetry-preserving regularization models , 2011 .

[7]  Pierre Sagaut,et al.  Discrete filters for large eddy simulation , 1999 .

[8]  Alison L. Marsden,et al.  Construction of commutative filters for LES on unstructured meshes , 2000 .

[9]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[10]  Julia S. Mullen,et al.  Filtering techniques for complex geometry fluid flows , 1999 .

[11]  A. W. Vreman The adjoint filter operator in large-eddy simulation of turbulent flow , 2004 .

[12]  B. Geurts,et al.  Realizability conditions for the turbulent stress tensor in large-eddy simulation , 1994, Journal of Fluid Mechanics.

[13]  P. Moin,et al.  A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries , 2007 .

[14]  Roel Verstappen On restraining the production of small scales of motion in a turbulent channel flow , 2008 .

[15]  P. Moin,et al.  A General Class of Commutative Filters for LES in Complex Geometries , 1998 .

[16]  A. Gorobets,et al.  DNS and regularization modeling of a turbulent differentially heated cavity of aspect ratio 5 , 2013 .

[17]  F. Xavier Trias,et al.  Symmetry-preserving discretization of Navier-Stokes equations on collocated unstructured grids , 2014, J. Comput. Phys..

[18]  P. Sagaut BOOK REVIEW: Large Eddy Simulation for Incompressible Flows. An Introduction , 2001 .

[19]  Per Lötstedt,et al.  Nonlinear filters for efficient shock computation , 1989 .

[20]  C. Meneveau,et al.  Large eddy simulation study of fully developed wind-turbine array boundary layers , 2010 .

[21]  A. Veldman,et al.  Symmetry-preserving discretization of turbulent flow , 2003 .

[22]  P. Moin,et al.  Numerical Simulation of Turbulent Flows , 1984 .

[23]  T. Hughes,et al.  Large Eddy Simulation and the variational multiscale method , 2000 .

[24]  Jens Nørkær Sørensen,et al.  Actuator Line Simulation of Wake of Wind Turbine Operating in Turbulent Inflow , 2007 .

[25]  M. Germano Differential filters for the large eddy numerical simulation of turbulent flows , 1986 .

[26]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[27]  T. Lund The use of explicit filters in large eddy simulation , 2003 .

[28]  A. Oliva,et al.  Low-frequency unsteadiness in the vortex formation region of a circular cylinder , 2013 .

[29]  Jungil Lee,et al.  A dynamic subgrid-scale eddy viscosity model with a global model coefficient , 2006 .

[30]  Christophe Bailly,et al.  A shock-capturing methodology based on adaptative spatial filtering for high-order non-linear computations , 2009, J. Comput. Phys..

[31]  H. C. Yee,et al.  A class of high resolution explicit and implicit shock-capturing methods , 1989 .

[32]  P. Moin,et al.  The basic equations for the large eddy simulation of turbulent flows in complex geometry , 1995 .

[33]  B. Geurts,et al.  Large-eddy simulation of the turbulent mixing layer , 1997, Journal of Fluid Mechanics.

[34]  E. Koschmieder Taylor vortices between eccentric cylinders , 1976 .

[35]  J. Blair Perot,et al.  Discrete Conservation Properties of Unstructured Mesh Schemes , 2011 .

[36]  A. Harten High Resolution Schemes for Hyperbolic Conservation Laws , 2017 .