- finish update

Author: rois1995

Common/include/geometry/dual_grid/CPoint.hpp

@@ -514,27 +514,6 @@ class CPoint {
   inline su2double& GetWall_Distance(unsigned long iPoint) { return Wall_Distance(iPoint); }
   inline const su2double& GetWall_Distance(unsigned long iPoint) const { return Wall_Distance(iPoint); }
 
-  /*!
-   * \brief Get the index of the closest wall element.
-   * \param[in] iPoint - Index of the point.
-   * \param[out] ClosestWall_Elem - ID of the closest element on a wall boundary.
-  */
-  inline unsigned long GetClosestWall_Elem(unsigned long iPoint) {return ClosestWall_Elem(iPoint);}
-
-  /*!
-   * \brief Get the marker of the closest wall marker.
-   * \param[in] iPoint - Index of the point.
-   * \param[out] ClosestWall_Marker - MarkerID of the closest wall boundary.
-   */
-  inline unsigned long GetClosestWall_Marker(unsigned long iPoint) {return ClosestWall_Marker(iPoint);}
-
-  /*!
-   * \brief Get the rank of the closest wall marker.
-   * \param[in] iPoint - Index of the point.
-   * \param[out] ClosestWall_Rank - RankID of the closest wall boundary.
-   */
-  inline unsigned long GetClosestWall_Rank(unsigned long iPoint) {return ClosestWall_Rank(iPoint);}
-
   /*!
    * \brief Set the value of the distance to the nearest wall.
    * \param[in] iPoint - Index of the point.
@@ -567,24 +546,6 @@ class CPoint {
    */
   inline su2double *GetNormal(unsigned long iPoint) { return Normals[iPoint]; }
 
-  /*!
-   * \brief Set the value of the normal of the nearest wall element.
-   * \param[in] iPoint - Index of the point.
-   * \param[in] normal - Value of the normal.
-   */
-  template<typename Normals_type>
-  inline void SetNormal(unsigned long iPoint, Normals_type const&normal) { 
-    for (unsigned long iDim = 0; iDim < nDim; iDim++)
-      Normals(iPoint,iDim) = normal[iDim];
-  }
-
-  /*!
-   * \brief Set the value of the normal of the nearest wall element.
-   * \param[in] iPoint - Index of the point.
-   * \return normal to the normal of the nearest wall element.
-   */
-  inline su2double *GetNormal(unsigned long iPoint) { return Normals[iPoint]; }
-
   /*!
    * \brief Set the value of the distance to a sharp edge.
    * \param[in] iPoint - Index of the point.

SU2_CFD/include/numerics/turbulent/transition/trans_diffusion.hpp

@@ -105,74 +105,6 @@ class CAvgGrad_TransLM final : public CAvgGrad_Scalar<FlowIndices> {
 };
 
 
-/*!
- * \class CAvgGrad_TransSLM
- * \brief Class for computing viscous term using average of gradient with correction (Simplified LM transition model).
- * \ingroup ViscDiscr
- * \author S. Kang.
- */
-template <class FlowIndices>
-class CAvgGrad_TransSLM final : public CAvgGrad_Scalar<FlowIndices> {
-private:
-  using Base = CAvgGrad_Scalar<FlowIndices>;
-  using Base::Laminar_Viscosity_i;
-  using Base::Laminar_Viscosity_j;
-  using Base::Eddy_Viscosity_i;
-  using Base::Eddy_Viscosity_j;
-  using Base::Density_i;
-  using Base::Density_j;
-  using Base::ScalarVar_i;
-  using Base::ScalarVar_j;
-  using Base::Proj_Mean_GradScalarVar;
-  using Base::proj_vector_ij;
-  using Base::Flux;
-  using Base::Jacobian_i;
-  using Base::Jacobian_j;
-
-  /*!
-   * \brief Adds any extra variables to AD
-   */
-  void ExtraADPreaccIn() override {}
-
-  /*!
-   * \brief LM transition model specific steps in the ComputeResidual method
-   * \param[in] config - Definition of the particular problem.
-   */
-  void FinishResidualCalc(const CConfig* config) override {
-    const bool implicit = config->GetKind_TimeIntScheme() == EULER_IMPLICIT;
-
-    /*--- Compute mean effective dynamic viscosity ---*/
-    const su2double diff_i_gamma = Laminar_Viscosity_i + Eddy_Viscosity_i;
-    const su2double diff_j_gamma = Laminar_Viscosity_j + Eddy_Viscosity_j;
-
-    const su2double diff_gamma = 0.5*(diff_i_gamma + diff_j_gamma);
-
-    Flux[0] = diff_gamma*Proj_Mean_GradScalarVar[0];
-
-    /*--- For Jacobians -> Use of TSL (Thin Shear Layer) approx. to compute derivatives of the gradients ---*/
-    if (implicit) {
-      const su2double proj_on_rho_i = proj_vector_ij/Density_i;
-      Jacobian_i[0][0] = -diff_gamma*proj_on_rho_i;
-
-      const su2double proj_on_rho_j = proj_vector_ij/Density_j;
-      Jacobian_j[0][0] = diff_gamma*proj_on_rho_j;
-    }
-  }
-
-public:
-  /*!
-   * \brief Constructor of the class.
-   * \param[in] val_nDim - Number of dimensions of the problem.
-   * \param[in] val_nVar - Number of variables of the problem.
-   * \param[in] correct_grad - Whether to correct gradient for skewness.
-   * \param[in] config - Definition of the particular problem.
-   */
-  CAvgGrad_TransSLM(unsigned short val_nDim, unsigned short val_nVar, bool correct_grad, const CConfig* config)
-    : CAvgGrad_Scalar<FlowIndices>(val_nDim, val_nVar, correct_grad, config){
-  }
-
-};
-
 /*!
  * \class CAvgGrad_TransSLM
  * \brief Class for computing viscous term using average of gradient with correction (Simplified LM transition model).

SU2_CFD/include/variables/CVariable.hpp

@@ -1699,51 +1699,6 @@ class CVariable {
    */
   inline virtual void SetIntermittencyEff(unsigned long iPoint, su2double val_Intermittency_eff) {}
 
-
-  /*!
-   * \brief Set Value of Transition Momentum Thickness Reynolds number from correlations.
-   */
-  inline virtual void SetCorr_Rec(unsigned long iPoint, su2double val_Corr_Rec) {};
-  inline virtual void SetTu(unsigned long iPoint, su2double val_Tu) {};
-  inline virtual void SetLambda_theta(unsigned long iPoint, su2double val_Lambda_theta) {};
-  inline virtual void Setduds(unsigned long iPoint, su2double val_duds) {};
-  inline virtual void SetRe_v(unsigned long iPoint, su2double val_Re_v) {};
-  inline virtual void SetProd(unsigned long iPoint, su2double val_Prod) {};
-  inline virtual void SetDestr(unsigned long iPoint, su2double val_Destr) {};
-  inline virtual void SetF_onset1(unsigned long iPoint, su2double val_F_onset1) {};
-  inline virtual void SetF_onset2(unsigned long iPoint, su2double val_F_onset2) {};
-  inline virtual void SetF_onset3(unsigned long iPoint, su2double val_F_onset3) {};
-  inline virtual void SetF_onset(unsigned long iPoint, su2double val_F_onset1) {};
-  inline virtual void SetNormal(unsigned long iPoint, su2double val_normal_x, su2double val_normal_y, su2double val_normal_z) {};
-
-  /*!
-   * \brief Get Value of Transition Momentum Thickness Reynolds number from correlations.
-   */
-  inline virtual su2double GetCorr_Rec(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetTu(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetLambda_theta(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double Getduds(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetRe_v(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetProd(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetDestr(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetF_onset1(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetF_onset2(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetF_onset3(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetF_onset(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetNormal_x(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetNormal_y(unsigned long iPoint) const { return 0.0; }
-  inline virtual su2double GetNormal_z(unsigned long iPoint) const { return 0.0; }
-
-  /*!
-   * \brief Set Value of Momentum Thickness Reynolds number from correlations (substitute to the second equation of original LM model).
-   */
-  inline virtual void SetRe_t(unsigned long iPoint, su2double val_Re_t) {};
-
-  /*!
-   * \brief Get Value of Momentum Thickness Reynolds number from correlations (substitute to the second equation of original LM model).
-   */
-  inline virtual su2double GetRe_t(unsigned long iPoint) const { return 0.0; }
-  
   /*!
    * \brief Set Value of Transition Momentum Thickness Reynolds number from correlations.
    */