SU2/doxygen/sgs__model_8inl_source.html

/*!

 * \file sgs_model.inl

 * \brief In-Line subroutines of the <i>sgs_model.hpp</i> file.

 * \author E. van der Weide, T. Economon, P. Urbanczyk

 * \version 8.5.0 "Harrier"

 *

 * SU2 Project Website: https://su2code.github.io

 *

 * The SU2 Project is maintained by the SU2 Foundation

 * (http://su2foundation.org)

 *

 * Copyright 2012-2026, SU2 Contributors (cf. AUTHORS.md)

 *

 * SU2 is free software; you can redistribute it and/or

 * modify it under the terms of the GNU Lesser General Public

 * License as published by the Free Software Foundation; either

 * version 2.1 of the License, or (at your option) any later version.

 *

 * SU2 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

 * Lesser General Public License for more details.

 *

 * You should have received a copy of the GNU Lesser General Public

 * License along with SU2. If not, see <http://www.gnu.org/licenses/>.

 */


#pragma once


inline CSGSModel::CSGSModel(void){}

inline CSGSModel::~CSGSModel(void){}


inline su2double CSGSModel::ComputeEddyViscosity_2D(const su2double rho,

                                                    const su2double dudx,

                                                    const su2double dudy,

                                                    const su2double dvdx,

                                                    const su2double dvdy,

                                                    const su2double lenScale,

                                                    const su2double distToWall) {

  return 0.0;

}


inline su2double CSGSModel::ComputeEddyViscosity_3D(const su2double rho,

                                                    const su2double dudx,

                                                    const su2double dudy,

                                                    const su2double dudz,

                                                    const su2double dvdx,

                                                    const su2double dvdy,

                                                    const su2double dvdz,

                                                    const su2double dwdx,

                                                    const su2double dwdy,

                                                    const su2double dwdz,

                                                    const su2double lenScale,

                                                    const su2double distToWall) {

  return 0.0;

}


inline void CSGSModel::ComputeGradEddyViscosity_2D(const su2double rho,

                                                   const su2double drhodx,

                                                   const su2double drhody,

                                                   const su2double dudx,

                                                   const su2double dudy,

                                                   const su2double dvdx,

                                                   const su2double dvdy,

                                                   const su2double d2udx2,

                                                   const su2double d2udy2,

                                                   const su2double d2udxdy,

                                                   const su2double d2vdx2,

                                                   const su2double d2vdy2,

                                                   const su2double d2vdxdy,

                                                   const su2double lenScale,

                                                   const su2double distToWall,

                                                         su2double &dMuTdx,

                                                         su2double &dMuTdy) {

  dMuTdx = dMuTdy = 0.0;

}


inline void CSGSModel::ComputeGradEddyViscosity_3D(const su2double rho,

                                                   const su2double drhodx,

                                                   const su2double drhody,

                                                   const su2double drhodz,

                                                   const su2double dudx,

                                                   const su2double dudy,

                                                   const su2double dudz,

                                                   const su2double dvdx,

                                                   const su2double dvdy,

                                                   const su2double dvdz,

                                                   const su2double dwdx,

                                                   const su2double dwdy,

                                                   const su2double dwdz,

                                                   const su2double d2udx2,

                                                   const su2double d2udy2,

                                                   const su2double d2udz2,

                                                   const su2double d2udxdy,

                                                   const su2double d2udxdz,

                                                   const su2double d2udydz,

                                                   const su2double d2vdx2,

                                                   const su2double d2vdy2,

                                                   const su2double d2vdz2,

                                                   const su2double d2vdxdy,

                                                   const su2double d2vdxdz,

                                                   const su2double d2vdydz,

                                                   const su2double d2wdx2,

                                                   const su2double d2wdy2,

                                                   const su2double d2wdz2,

                                                   const su2double d2wdxdy,

                                                   const su2double d2wdxdz,

                                                   const su2double d2wdydz,

                                                   const su2double lenScale,

                                                   const su2double distToWall,

                                                         su2double &dMuTdx,

                                                         su2double &dMuTdy,

                                                         su2double &dMuTdz) {

  dMuTdx = dMuTdy = dMuTdz = 0.0;

}


inline CSmagorinskyModel::CSmagorinskyModel(void) : CSGSModel() {

  const_smag  = 0.1;

  filter_mult = 2.0;

}


inline CSmagorinskyModel::~CSmagorinskyModel(void){}


inline su2double CSmagorinskyModel::ComputeEddyViscosity_2D(const su2double rho,

                                                            const su2double dudx,

                                                            const su2double dudy,

                                                            const su2double dvdx,

                                                            const su2double dvdy,

                                                            const su2double lenScale,

                                                            const su2double distToWall) {

  /* Constant coefficient Smagorinsky SGS is calculated:

   * ( C_s * L_c )^2 * |S(x,t)|

   * C_s = Smagorinsky constant

   * L_c = Filter width

   * S(x,t) = Rate of Strain Tensor ( 1/2 [ du_i/dx_j + du_j/dx_i] )

   */

  const su2double C_s_filter_width = const_smag*filter_mult*lenScale;


  const su2double S12          = 0.5*(dudy + dvdx);

  const su2double strain_rate2 = 2.0*(dudx*dudx + dvdy*dvdy + 2.0*S12*S12);


  /* Return the SGS dynamic viscosity. */

  return rho*C_s_filter_width*C_s_filter_width*sqrt(strain_rate2);

}


inline su2double CSmagorinskyModel::ComputeEddyViscosity_3D(const su2double rho,

                                                            const su2double dudx,

                                                            const su2double dudy,

                                                            const su2double dudz,

                                                            const su2double dvdx,

                                                            const su2double dvdy,

                                                            const su2double dvdz,

                                                            const su2double dwdx,

                                                            const su2double dwdy,

                                                            const su2double dwdz,

                                                            const su2double lenScale,

                                                            const su2double distToWall) {

  /* Constant coefficient Smagorinsky SGS is calculated:

   * ( C_s * L_c )^2 * |S(x,t)|

   * C_s = Smagorinsky constant

   * L_c = Filter width

   * S(x,t) = Rate of Strain Tensor ( 1/2 [ du_i/dx_j + du_j/dx_i] )

   */

  const su2double C_s_filter_width = const_smag*filter_mult*lenScale;


  const su2double S12 = 0.5*(dudy + dvdx);

  const su2double S13 = 0.5*(dudz + dwdx);

  const su2double S23 = 0.5*(dvdz + dwdy);


  const su2double strain_rate2 = 2.0*(dudx*dudx + dvdy*dvdy + dwdz*dwdz

                               +      2.0*(S12*S12 + S13*S13 + S23*S23));


  /* Return the SGS dynamic viscosity. */

  return rho*C_s_filter_width*C_s_filter_width*sqrt(strain_rate2);

}


inline void CSmagorinskyModel::ComputeGradEddyViscosity_2D(const su2double rho,

                                                           const su2double drhodx,

                                                           const su2double drhody,

                                                           const su2double dudx,

                                                           const su2double dudy,

                                                           const su2double dvdx,

                                                           const su2double dvdy,

                                                           const su2double d2udx2,

                                                           const su2double d2udy2,

                                                           const su2double d2udxdy,

                                                           const su2double d2vdx2,

                                                           const su2double d2vdy2,

                                                           const su2double d2vdxdy,

                                                           const su2double lenScale,

                                                           const su2double distToWall,

                                                                 su2double &dMuTdx,

                                                                 su2double &dMuTdy) {

  SU2_MPI::Error("Not implemented yet", CURRENT_FUNCTION);

}


inline void CSmagorinskyModel::ComputeGradEddyViscosity_3D(const su2double rho,

                                                           const su2double drhodx,

                                                           const su2double drhody,

                                                           const su2double drhodz,

                                                           const su2double dudx,

                                                           const su2double dudy,

                                                           const su2double dudz,

                                                           const su2double dvdx,

                                                           const su2double dvdy,

                                                           const su2double dvdz,

                                                           const su2double dwdx,

                                                           const su2double dwdy,

                                                           const su2double dwdz,

                                                           const su2double d2udx2,

                                                           const su2double d2udy2,

                                                           const su2double d2udz2,

                                                           const su2double d2udxdy,

                                                           const su2double d2udxdz,

                                                           const su2double d2udydz,

                                                           const su2double d2vdx2,

                                                           const su2double d2vdy2,

                                                           const su2double d2vdz2,

                                                           const su2double d2vdxdy,

                                                           const su2double d2vdxdz,

                                                           const su2double d2vdydz,

                                                           const su2double d2wdx2,

                                                           const su2double d2wdy2,

                                                           const su2double d2wdz2,

                                                           const su2double d2wdxdy,

                                                           const su2double d2wdxdz,

                                                           const su2double d2wdydz,

                                                           const su2double lenScale,

                                                           const su2double distToWall,

                                                                 su2double &dMuTdx,

                                                                 su2double &dMuTdy,

                                                                 su2double &dMuTdz) {

  SU2_MPI::Error("Not implemented yet", CURRENT_FUNCTION);

}


inline CWALEModel::CWALEModel(void) : CSGSModel() {

  const_WALE = 0.325;

}


inline CWALEModel::~CWALEModel(void){}


inline su2double CWALEModel::ComputeEddyViscosity_2D(const su2double rho,

                                                     const su2double dudx,

                                                     const su2double dudy,

                                                     const su2double dvdx,

                                                     const su2double dvdy,

                                                     const su2double lenScale,

                                                     const su2double distToWall) {

  /* Compute the length scale in WALE. */

  const su2double lenScaleWale = const_WALE*lenScale;


  /* Compute the strain rate tensor, which is symmetric. */

  const su2double S11 = dudx, S22 = dvdy;

  const su2double S12 = 0.5*(dudy + dvdx);


  /* Compute the values of the Sd tensor. First without the trace

     correction of the diagonal terms. */

  su2double Sd11 = dudx*dudx + dudy*dvdx;

  su2double Sd22 = dvdx*dudy + dvdy*dvdy;


  const su2double Sd12 = 0.5*(dudx*dudy + dudy*dvdy + dvdx*dudx + dvdy*dvdx);


  /* Correct the diagonal elements, such that the trace of the Sd tensor is zero

     Note that this comes from the 3D formulation. */

  const su2double thirdTrace = (Sd11 + Sd22)/3.0;


  Sd11 -= thirdTrace;

  Sd22 -= thirdTrace;


  /* Compute the summation of both tensors. */

  const su2double sumS  = S11 *S11  + S22 *S22  + 2.0*S12 *S12;

  const su2double sumSd = Sd11*Sd11 + Sd22*Sd22 + 2.0*Sd12*Sd12;


  /* Compute the kinematic eddy viscosity. */

  const su2double sumSdPow3_2 = sumSd*sqrt(sumSd);

  const su2double sumSdPow5_4 = sqrt(sumSdPow3_2*sumSd);

  const su2double sumSPow5_2  = sumS*sumS*sqrt(sumS);

  const su2double denom       = sumSPow5_2 + sumSdPow5_4;


  const su2double nuEddy = lenScaleWale*lenScaleWale*sumSdPow3_2

                         / max(denom, 1.e-20);


  /* Return the SGS dynamic viscosity. */

  return rho*nuEddy;

}


inline su2double CWALEModel::ComputeEddyViscosity_3D(const su2double rho,

                                                     const su2double dudx,

                                                     const su2double dudy,

                                                     const su2double dudz,

                                                     const su2double dvdx,

                                                     const su2double dvdy,

                                                     const su2double dvdz,

                                                     const su2double dwdx,

                                                     const su2double dwdy,

                                                     const su2double dwdz,

                                                     const su2double lenScale,

                                                     const su2double distToWall) {

  /* Compute the length scale in WALE. */

  const su2double lenScaleWale = const_WALE*lenScale;


  /* Compute the strain rate tensor, which is symmetric. */

  const su2double S11 = dudx, S22 = dvdy, S33 = dwdz;

  const su2double S12 = 0.5*(dudy + dvdx);

  const su2double S13 = 0.5*(dudz + dwdx);

  const su2double S23 = 0.5*(dvdz + dwdy);


  /* Compute the values of the Sd tensor. First without the trace

     correction of the diagonal terms. */

  su2double Sd11 = dudx*dudx + dudy*dvdx + dudz*dwdx;

  su2double Sd22 = dvdx*dudy + dvdy*dvdy + dvdz*dwdy;

  su2double Sd33 = dwdx*dudz + dwdy*dvdz + dwdz*dwdz;


  const su2double Sd12 = 0.5*(dudx*dudy + dudy*dvdy + dudz*dwdy

                       +      dvdx*dudx + dvdy*dvdx + dvdz*dwdx);

  const su2double Sd13 = 0.5*(dudx*dudz + dudy*dvdz + dudz*dwdz

                       +      dwdx*dudx + dwdy*dvdx + dwdz*dwdx);

  const su2double Sd23 = 0.5*(dvdx*dudz + dvdy*dvdz + dvdz*dwdz

                       +      dwdx*dudy + dwdy*dvdy + dwdz*dwdy);


  /* Correct the diagonal elements, such that the trace of the Sd tensor is zero. */

  const su2double thirdTrace = (Sd11 + Sd22 + Sd33)/3.0;


  Sd11 -= thirdTrace;

  Sd22 -= thirdTrace;

  Sd33 -= thirdTrace;


  /* Compute the summation of both tensors. */

  const su2double sumS  = S11*S11 + S22*S22 + S33*S33

                        + 2.0*(S12*S12 + S13*S13 + S23*S23);

  const su2double sumSd = Sd11*Sd11 + Sd22*Sd22 + Sd33*Sd33

                        + 2.0*(Sd12*Sd12 + Sd13*Sd13 + Sd23*Sd23);


  /* Compute the kinematic eddy viscosity. */

  const su2double sumSdPow3_2 = sumSd*sqrt(sumSd);

  const su2double sumSdPow5_4 = sqrt(sumSdPow3_2*sumSd);

  const su2double sumSPow5_2  = sumS*sumS*sqrt(sumS);

  const su2double denom       = sumSPow5_2 + sumSdPow5_4;


  const su2double nuEddy = lenScaleWale*lenScaleWale*sumSdPow3_2

                         / max(denom, 1.e-20);


  /* Return the SGS dynamic viscosity. */

  return rho*nuEddy;

}


inline void CWALEModel::ComputeGradEddyViscosity_2D(const su2double rho,

                                                    const su2double drhodx,

                                                    const su2double drhody,

                                                    const su2double dudx,

                                                    const su2double dudy,

                                                    const su2double dvdx,

                                                    const su2double dvdy,

                                                    const su2double d2udx2,

                                                    const su2double d2udy2,

                                                    const su2double d2udxdy,

                                                    const su2double d2vdx2,

                                                    const su2double d2vdy2,

                                                    const su2double d2vdxdy,

                                                    const su2double lenScale,

                                                    const su2double distToWall,

                                                          su2double &dMuTdx,

                                                          su2double &dMuTdy) {

  SU2_MPI::Error("Not implemented yet", CURRENT_FUNCTION);

}


inline void CWALEModel::ComputeGradEddyViscosity_3D(const su2double rho,

                                                    const su2double drhodx,

                                                    const su2double drhody,

                                                    const su2double drhodz,

                                                    const su2double dudx,

                                                    const su2double dudy,

                                                    const su2double dudz,

                                                    const su2double dvdx,

                                                    const su2double dvdy,

                                                    const su2double dvdz,

                                                    const su2double dwdx,

                                                    const su2double dwdy,

                                                    const su2double dwdz,

                                                    const su2double d2udx2,

                                                    const su2double d2udy2,

                                                    const su2double d2udz2,

                                                    const su2double d2udxdy,

                                                    const su2double d2udxdz,

                                                    const su2double d2udydz,

                                                    const su2double d2vdx2,

                                                    const su2double d2vdy2,

                                                    const su2double d2vdz2,

                                                    const su2double d2vdxdy,

                                                    const su2double d2vdxdz,

                                                    const su2double d2vdydz,

                                                    const su2double d2wdx2,

                                                    const su2double d2wdy2,

                                                    const su2double d2wdz2,

                                                    const su2double d2wdxdy,

                                                    const su2double d2wdxdz,

                                                    const su2double d2wdydz,

                                                    const su2double lenScale,

                                                    const su2double distToWall,

                                                          su2double &dMuTdx,

                                                          su2double &dMuTdy,

                                                          su2double &dMuTdz) {

  SU2_MPI::Error("Not implemented yet", CURRENT_FUNCTION);

}


inline CVremanModel::CVremanModel(void) : CSGSModel() {


  /* const_Vreman = 2.5*Cs*Cs where Cs is the Smagorinsky constant */

  const_Vreman = 0.07;

}


inline CVremanModel::~CVremanModel(void){}


inline su2double CVremanModel::ComputeEddyViscosity_2D(const su2double rho,

                                                     const su2double dudx,

                                                     const su2double dudy,

                                                     const su2double dvdx,

                                                     const su2double dvdy,

                                                     const su2double lenScale,

                                                     const su2double distToWall) {

  SU2_MPI::Error("Not implemented yet", CURRENT_FUNCTION);

  return 0;

}


inline su2double CVremanModel::ComputeEddyViscosity_3D(const su2double rho,

                                                     const su2double dudx,

                                                     const su2double dudy,

                                                     const su2double dudz,

                                                     const su2double dvdx,

                                                     const su2double dvdy,

                                                     const su2double dvdz,

                                                     const su2double dwdx,

                                                     const su2double dwdy,

                                                     const su2double dwdz,

                                                     const su2double lenScale,

                                                     const su2double distToWall) {


  su2double alpha11 = dudx;

  su2double alpha22 = dvdy;

  su2double alpha33 = dwdz;


  //Check if it is necessary to remove the trace.

  const su2double tmp = (alpha11 + alpha22 + alpha33)/3.0;

  alpha11 -= tmp;

  alpha22 -= tmp;

  alpha33 -= tmp;


  const su2double lenScale2 = lenScale * lenScale;

  const su2double alpha12 = dudy;

  const su2double alpha13 = dudz;

  const su2double alpha23 = dvdz;


  const su2double alpha21 = dvdx;

  const su2double alpha31 = dwdx;

  const su2double alpha32 = dwdy;


  const su2double beta11  = lenScale2*alpha11*alpha11 + lenScale2*alpha12*alpha12 + lenScale2*alpha13*alpha13 ;

  const su2double beta12  = lenScale2*alpha11*alpha21 + lenScale2*alpha12*alpha22 + lenScale2*alpha13*alpha23 ;

  const su2double beta13  = lenScale2*alpha11*alpha31 + lenScale2*alpha12*alpha32 + lenScale2*alpha13*alpha33 ;

  const su2double beta22  = lenScale2*alpha21*alpha21 + lenScale2*alpha22*alpha22 + lenScale2*alpha23*alpha23 ;

  const su2double beta23  = lenScale2*alpha21*alpha31 + lenScale2*alpha22*alpha32 + lenScale2*alpha23*alpha33 ;

  const su2double beta33  = lenScale2*alpha31*alpha31 + lenScale2*alpha32*alpha32 + lenScale2*alpha33*alpha33 ;


  su2double B = beta11*beta22-beta12*beta12+beta11*beta33-beta13*beta13+beta22*beta33-beta23*beta23;

            B = (B + fabs(B))*0.5;

  const su2double denon    = alpha11*alpha11+alpha22*alpha22+alpha33*alpha33 +

                            alpha12*alpha12+alpha13*alpha13+alpha23*alpha23 +

                            alpha21*alpha21+alpha31*alpha31+alpha32*alpha32;


  const su2double nuEddy_Vreman = sqrt(B/(denon+1.0E-20));


 /* Return the SGS dynamic viscosity. */

 return rho*const_Vreman*nuEddy_Vreman;


}


inline void CVremanModel::ComputeGradEddyViscosity_2D(const su2double rho,

                                                    const su2double drhodx,

                                                    const su2double drhody,

                                                    const su2double dudx,

                                                    const su2double dudy,

                                                    const su2double dvdx,

                                                    const su2double dvdy,

                                                    const su2double d2udx2,

                                                    const su2double d2udy2,

                                                    const su2double d2udxdy,

                                                    const su2double d2vdx2,

                                                    const su2double d2vdy2,

                                                    const su2double d2vdxdy,

                                                    const su2double lenScale,

                                                    const su2double distToWall,

                                                    su2double &dMuTdx,

                                                    su2double &dMuTdy) {

  SU2_MPI::Error("Not implemented yet", CURRENT_FUNCTION);

}


inline void CVremanModel::ComputeGradEddyViscosity_3D(const su2double rho,

                                                    const su2double drhodx,

                                                    const su2double drhody,

                                                    const su2double drhodz,

                                                    const su2double dudx,

                                                    const su2double dudy,

                                                    const su2double dudz,

                                                    const su2double dvdx,

                                                    const su2double dvdy,

                                                    const su2double dvdz,

                                                    const su2double dwdx,

                                                    const su2double dwdy,

                                                    const su2double dwdz,

                                                    const su2double d2udx2,

                                                    const su2double d2udy2,

                                                    const su2double d2udz2,

                                                    const su2double d2udxdy,

                                                    const su2double d2udxdz,

                                                    const su2double d2udydz,

                                                    const su2double d2vdx2,

                                                    const su2double d2vdy2,

                                                    const su2double d2vdz2,

                                                    const su2double d2vdxdy,

                                                    const su2double d2vdxdz,

                                                    const su2double d2vdydz,

                                                    const su2double d2wdx2,

                                                    const su2double d2wdy2,

                                                    const su2double d2wdz2,

                                                    const su2double d2wdxdy,

                                                    const su2double d2wdxdz,

                                                    const su2double d2wdydz,

                                                    const su2double lenScale,

                                                    const su2double distToWall,

                                                    su2double &dMuTdx,

                                                    su2double &dMuTdy,

                                                    su2double &dMuTdz) {

  SU2_MPI::Error("Not implemented yet", CURRENT_FUNCTION);

}


B
#define B(i, j)
Definition blas_structure.cpp:132

CBaseMPIWrapper::Error
static void Error(const std::string &ErrorMsg, const std::string &FunctionName)
Definition mpi_structure.cpp:131

CSGSModel
Base class for defining the LES subgrid scale model.
Definition sgs_model.hpp:43

CSGSModel::CSGSModel
CSGSModel(void)
Constructor of the class.
Definition sgs_model.inl:30

CSGSModel::ComputeGradEddyViscosity_3D
virtual void ComputeGradEddyViscosity_3D(const su2double rho, const su2double drhodx, const su2double drhody, const su2double drhodz, const su2double dudx, const su2double dudy, const su2double dudz, const su2double dvdx, const su2double dvdy, const su2double dvdz, const su2double dwdx, const su2double dwdy, const su2double dwdz, const su2double d2udx2, const su2double d2udy2, const su2double d2udz2, const su2double d2udxdy, const su2double d2udxdz, const su2double d2udydz, const su2double d2vdx2, const su2double d2vdy2, const su2double d2vdz2, const su2double d2vdxdy, const su2double d2vdxdz, const su2double d2vdydz, const su2double d2wdx2, const su2double d2wdy2, const su2double d2wdz2, const su2double d2wdxdy, const su2double d2wdxdz, const su2double d2wdydz, const su2double lenScale, const su2double distToWall, su2double &dMuTdx, su2double &dMuTdy, su2double &dMuTdz)
Virtual function to determine the gradients of the eddy viscosity for the given function arguments fo...
Definition sgs_model.inl:78

CSGSModel::ComputeGradEddyViscosity_2D
virtual void ComputeGradEddyViscosity_2D(const su2double rho, const su2double drhodx, const su2double drhody, const su2double dudx, const su2double dudy, const su2double dvdx, const su2double dvdy, const su2double d2udx2, const su2double d2udy2, const su2double d2udxdy, const su2double d2vdx2, const su2double d2vdy2, const su2double d2vdxdy, const su2double lenScale, const su2double distToWall, su2double &dMuTdx, su2double &dMuTdy)
Virtual function to determine the gradients of the eddy viscosity for the given function arguments fo...
Definition sgs_model.inl:58

CSGSModel::ComputeEddyViscosity_3D
virtual su2double ComputeEddyViscosity_3D(const su2double rho, const su2double dudx, const su2double dudy, const su2double dudz, const su2double dvdx, const su2double dvdy, const su2double dvdz, const su2double dwdx, const su2double dwdy, const su2double dwdz, const su2double lenScale, const su2double distToWall)
Virtual function to determine the eddy viscosity for the given function arguments for a 3D simulation...
Definition sgs_model.inl:43

CSGSModel::ComputeEddyViscosity_2D
virtual su2double ComputeEddyViscosity_2D(const su2double rho, const su2double dudx, const su2double dudy, const su2double dvdx, const su2double dvdy, const su2double lenScale, const su2double distToWall)
Virtual function to determine the eddy viscosity for the given function arguments for a 2D simulation...
Definition sgs_model.inl:33

CSGSModel::~CSGSModel
virtual ~CSGSModel(void)
Destructor of the class.
Definition sgs_model.inl:31

CSmagorinskyModel::filter_mult
su2double filter_mult
Multiplier to get filter width from grid length scale.
Definition sgs_model.hpp:234

CSmagorinskyModel::~CSmagorinskyModel
~CSmagorinskyModel(void) override
Destructor of the class.
Definition sgs_model.inl:122

CSmagorinskyModel::CSmagorinskyModel
CSmagorinskyModel(void)
Constructor of the class.
Definition sgs_model.inl:117

CSmagorinskyModel::ComputeGradEddyViscosity_3D
void ComputeGradEddyViscosity_3D(const su2double rho, const su2double drhodx, const su2double drhody, const su2double drhodz, const su2double dudx, const su2double dudy, const su2double dudz, const su2double dvdx, const su2double dvdy, const su2double dvdz, const su2double dwdx, const su2double dwdy, const su2double dwdz, const su2double d2udx2, const su2double d2udy2, const su2double d2udz2, const su2double d2udxdy, const su2double d2udxdz, const su2double d2udydz, const su2double d2vdx2, const su2double d2vdy2, const su2double d2vdz2, const su2double d2vdxdy, const su2double d2vdxdz, const su2double d2vdydz, const su2double d2wdx2, const su2double d2wdy2, const su2double d2wdz2, const su2double d2wdxdy, const su2double d2wdxdz, const su2double d2wdydz, const su2double lenScale, const su2double distToWall, su2double &dMuTdx, su2double &dMuTdy, su2double &dMuTdz) override
function to determine the gradients of the eddy viscosity for the given function arguments for a 3D s...
Definition sgs_model.inl:197

CSmagorinskyModel::ComputeEddyViscosity_2D
su2double ComputeEddyViscosity_2D(const su2double rho, const su2double dudx, const su2double dudy, const su2double dvdx, const su2double dvdy, const su2double lenScale, const su2double distToWall) override
Function to determine the eddy viscosity for the given function arguments for a 2D simulation.
Definition sgs_model.inl:124

CSmagorinskyModel::const_smag
su2double const_smag
Smagorinsky Constant C_s.
Definition sgs_model.hpp:233

CSmagorinskyModel::ComputeEddyViscosity_3D
su2double ComputeEddyViscosity_3D(const su2double rho, const su2double dudx, const su2double dudy, const su2double dudz, const su2double dvdx, const su2double dvdy, const su2double dvdz, const su2double dwdx, const su2double dwdy, const su2double dwdz, const su2double lenScale, const su2double distToWall) override
Function to determine the eddy viscosity for the given function arguments for a 3D simulation.
Definition sgs_model.inl:146

CSmagorinskyModel::ComputeGradEddyViscosity_2D
void ComputeGradEddyViscosity_2D(const su2double rho, const su2double drhodx, const su2double drhody, const su2double dudx, const su2double dudy, const su2double dvdx, const su2double dvdy, const su2double d2udx2, const su2double d2udy2, const su2double d2udxdy, const su2double d2vdx2, const su2double d2vdy2, const su2double d2vdxdy, const su2double lenScale, const su2double distToWall, su2double &dMuTdx, su2double &dMuTdy) override
Function to determine the gradients of the eddy viscosity for the given function arguments for a 2D s...
Definition sgs_model.inl:177

CVremanModel::ComputeGradEddyViscosity_3D
void ComputeGradEddyViscosity_3D(const su2double rho, const su2double drhodx, const su2double drhody, const su2double drhodz, const su2double dudx, const su2double dudy, const su2double dudz, const su2double dvdx, const su2double dvdy, const su2double dvdz, const su2double dwdx, const su2double dwdy, const su2double dwdz, const su2double d2udx2, const su2double d2udy2, const su2double d2udz2, const su2double d2udxdy, const su2double d2udxdz, const su2double d2udydz, const su2double d2vdx2, const su2double d2vdy2, const su2double d2vdz2, const su2double d2vdxdy, const su2double d2vdxdz, const su2double d2vdydz, const su2double d2wdx2, const su2double d2wdy2, const su2double d2wdz2, const su2double d2wdxdy, const su2double d2wdxdz, const su2double d2wdydz, const su2double lenScale, const su2double distToWall, su2double &dMuTdx, su2double &dMuTdy, su2double &dMuTdz) override
function to determine the gradients of the eddy viscosity for the given function arguments for a 3D s...
Definition sgs_model.inl:497

CVremanModel::ComputeGradEddyViscosity_2D
void ComputeGradEddyViscosity_2D(const su2double rho, const su2double drhodx, const su2double drhody, const su2double dudx, const su2double dudy, const su2double dvdx, const su2double dvdy, const su2double d2udx2, const su2double d2udy2, const su2double d2udxdy, const su2double d2vdx2, const su2double d2vdy2, const su2double d2vdxdy, const su2double lenScale, const su2double distToWall, su2double &dMuTdx, su2double &dMuTdy) override
Function to determine the gradients of the eddy viscosity for the given function arguments for a 2D s...
Definition sgs_model.inl:477

CVremanModel::ComputeEddyViscosity_3D
su2double ComputeEddyViscosity_3D(const su2double rho, const su2double dudx, const su2double dudy, const su2double dudz, const su2double dvdx, const su2double dvdy, const su2double dvdz, const su2double dwdx, const su2double dwdy, const su2double dwdz, const su2double lenScale, const su2double distToWall) override
Function to determine the eddy viscosity for the given function arguments for a 3D simulation.
Definition sgs_model.inl:425

CVremanModel::ComputeEddyViscosity_2D
su2double ComputeEddyViscosity_2D(const su2double rho, const su2double dudx, const su2double dudy, const su2double dvdx, const su2double dvdy, const su2double lenScale, const su2double distToWall) override
Function to determine the eddy viscosity for the given function arguments for a 2D simulation.
Definition sgs_model.inl:414

CVremanModel::CVremanModel
CVremanModel(void)
Constructor of the class.
Definition sgs_model.inl:406

CVremanModel::const_Vreman
su2double const_Vreman
Vreman Constant c=2.5*Cs*Cs.
Definition sgs_model.hpp:609

CVremanModel::~CVremanModel
~CVremanModel(void) override
Destructor of the class.
Definition sgs_model.inl:412

CWALEModel::const_WALE
su2double const_WALE
WALE Constant Cw.
Definition sgs_model.hpp:421

CWALEModel::ComputeEddyViscosity_3D
su2double ComputeEddyViscosity_3D(const su2double rho, const su2double dudx, const su2double dudy, const su2double dudz, const su2double dvdx, const su2double dvdy, const su2double dvdz, const su2double dwdx, const su2double dwdy, const su2double dwdz, const su2double lenScale, const su2double distToWall) override
Function to determine the eddy viscosity for the given function arguments for a 3D simulation.
Definition sgs_model.inl:287

CWALEModel::CWALEModel
CWALEModel(void)
Constructor of the class.
Definition sgs_model.inl:236

CWALEModel::ComputeEddyViscosity_2D
su2double ComputeEddyViscosity_2D(const su2double rho, const su2double dudx, const su2double dudy, const su2double dvdx, const su2double dvdy, const su2double lenScale, const su2double distToWall) override
Function to determine the eddy viscosity for the given function arguments for a 2D simulation.
Definition sgs_model.inl:242

CWALEModel::ComputeGradEddyViscosity_3D
void ComputeGradEddyViscosity_3D(const su2double rho, const su2double drhodx, const su2double drhody, const su2double drhodz, const su2double dudx, const su2double dudy, const su2double dudz, const su2double dvdx, const su2double dvdy, const su2double dvdz, const su2double dwdx, const su2double dwdy, const su2double dwdz, const su2double d2udx2, const su2double d2udy2, const su2double d2udz2, const su2double d2udxdy, const su2double d2udxdz, const su2double d2udydz, const su2double d2vdx2, const su2double d2vdy2, const su2double d2vdz2, const su2double d2vdxdy, const su2double d2vdxdz, const su2double d2vdydz, const su2double d2wdx2, const su2double d2wdy2, const su2double d2wdz2, const su2double d2wdxdy, const su2double d2wdxdz, const su2double d2wdydz, const su2double lenScale, const su2double distToWall, su2double &dMuTdx, su2double &dMuTdy, su2double &dMuTdz) override
function to determine the gradients of the eddy viscosity for the given function arguments for a 3D s...
Definition sgs_model.inl:367

CWALEModel::ComputeGradEddyViscosity_2D
void ComputeGradEddyViscosity_2D(const su2double rho, const su2double drhodx, const su2double drhody, const su2double dudx, const su2double dudy, const su2double dvdx, const su2double dvdy, const su2double d2udx2, const su2double d2udy2, const su2double d2udxdy, const su2double d2vdx2, const su2double d2vdy2, const su2double d2vdxdy, const su2double lenScale, const su2double distToWall, su2double &dMuTdx, su2double &dMuTdy) override
Function to determine the gradients of the eddy viscosity for the given function arguments for a 2D s...
Definition sgs_model.inl:347

CWALEModel::~CWALEModel
~CWALEModel(void) override
Destructor of the class.
Definition sgs_model.inl:240

su2double
passivedouble su2double
Definition code_config.hpp:181

sqrt
FORCEINLINE ARRAY_T sqrt(const ARRAY_T &x)
Definition special_vectorization.hpp:116

CURRENT_FUNCTION
#define CURRENT_FUNCTION
Definition mpi_structure.hpp:59