cs_les_mu_t.h

Go to the documentation of this file.

#ifndef CS_TURBULENCE_VIS_LES_H
#define CS_TURBULENCE_VIS_LES_H
 
/*============================================================================
 * Turbulent viscosity for LES models.
 *============================================================================*/
 
/*
  This file is part of code_saturne, a general-purpose CFD tool.
 
  Copyright (C) 1998-2026 EDF S.A.
 
  This program is free software; you can redistribute it and/or modify it under
  the terms of the GNU General Public License as published by the Free Software
  Foundation; either version 2 of the License, or (at your option) any later
  version.
 
  This program is distributed in the hope that it will be useful, but WITHOUT
  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
  FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
  details.
 
  You should have received a copy of the GNU General Public License along with
  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
  Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/
 
/*----------------------------------------------------------------------------*/
 
/*----------------------------------------------------------------------------
 * Local headers
 *----------------------------------------------------------------------------*/
 
#include "base/cs_base.h"
 
/*=============================================================================
 * Local Macro definitions
 *============================================================================*/
 
/*============================================================================
 * Type definition
 *============================================================================*/
 
/*=============================================================================
 * Public function prototypes
 *============================================================================*/
 
/*----------------------------------------------------------------------------*/
/*
 * \brief Calculation of  Mij:Mij and Mij:Lij for dynamic Smagorinsky model
 *
 * Please refer to the
 * <a href="../../theory.pdf#dynsmago"><b>dynamic Smagorinsky model</b></a>
 * section of the theory guide for more informations.
 *
 * \param[out]    s_n           strain rate (sqrt(2SijSij))
 * \param[out]    sf_n          filtered strain rate
 * \param[out]    f_vel         filtered velocity
 * \param[out]    mijmij        Mij:Mij
 * \param[out]    mijlij        Mij:Lij
 *
 !*/
/*----------------------------------------------------------------------------*/
 
extern "C" void
cs_les_mu_t_smago_dyn_prepare(cs_real_t s_n[],
                              cs_real_t sf_n[],
                              cs_real_3_t *f_vel,
                              cs_real_t mijmij[],
                              cs_real_t mijlij[]);
 
/*----------------------------------------------------------------------------*/
/*
 * \brief Calculation of turbulent viscosity for
 *        a dynamic Smagorinsky LES model
 *
 * \f[ smago = \dfrac{L_{ij}M_{ij}}{M_{ij}M_{ij}} \f]
 *
 * \f[ \mu_T = \rho smago L^2  \sqrt{2 S_{ij}S_{ij}} \f]
 * \f[ S_{ij} = \dfrac{\der{u_i}{x_j} + \der{u_j}{x_i}}{2}\f]
 *
 *
 !*/
/*----------------------------------------------------------------------------*/
 
extern "C" void
cs_les_mu_t_smago_dyn(void);
 
/*----------------------------------------------------------------------------*/
/*
 * \brief Calculation of the turbulent viscosity for a Smagorinsky LES model
 *
 * \f[ \mu_T = \rho (C_{S} l)^2  \sqrt{2 S_{ij}S_{ij}} \f]
 * \f[ S_{ij} = \dfrac{\der{u_i}{x_j} + \der{u_j}{x_i}}{2}\f]
 *
 !*/
/*----------------------------------------------------------------------------*/
 
extern "C" void
cs_les_mu_t_smago_const(void);
 
/*----------------------------------------------------------------------------*/
/*
 * \brief Compute the turbulent viscosity for the WALE LES model.
 *
 * The turbulent viscosity is:
 * \f$ \mu_T = \rho (C_{wale} L)^2 * \dfrac{(\tens{S}:\tens{Sd})^{3/2}}
 *                                         {(\tens{S} :\tens{S})^(5/2)
 *                                         +(\tens{Sd}:\tens{Sd})^(5/4)} \f$
 * with \f$ \tens{S}  = \frac{1}{2}(\gradt \vect{u} + \transpose{\gradt \vect{u}})\f$
 * and  \f$ \tens{Sd} = \deviator{(\symmetric{(\tens{S}^2)})}\f$
 *
 !*/
/*----------------------------------------------------------------------------*/
 
extern "C" void
cs_les_mu_t_wale(void);
 
/*----------------------------------------------------------------------------*/
/*
 * \brief Calculation of turbulent viscosity for
 *        a k-SGS LES model
 !*/
/*----------------------------------------------------------------------------*/
 
extern "C" void
cs_les_mu_t_ksgs(void);
 
/*----------------------------------------------------------------------------*/
/*
 * \brief Calculation of turbulent viscosity for
 *        a tau-SGS LES model
 !*/
/*----------------------------------------------------------------------------*/
 
extern "C" void
cs_les_mu_t_tausgs(void);
 
/*----------------------------------------------------------------------------*/
 
#endif /* CS_LES_MU_T_H */