Kernel.hpp

Go to the documentation of this file.

//
// Author(s)    : - Jean-luc CHARLES I2M-DuMAS/ENSAM Talence France
//                  <jean-luc.charles@ensam.eu>
//                - Damien ANDRE     SPCTS/ENS Ceramique industrielle, Limoges France
//                  <damien.andre@unilim.fr>
//                - Jeremie GIRARDOT I2M-DuMAS/ENSAM Talence France
//                  <jeremie.girardot@ensam.eu>
//                - Cedric Hubert    LAMIH/UVHC Valenciennes France
//                  <cedric.hubert@univ-valenciennes.fr>
//                - Ivan IORDANOFF   I2M-DuMAS-MPI/ENSAM Talence France
//                  <ivan.iordanoff@ensam.eu>
//
// Copyright (C) 2008-2016 JL. CHARLES, D. ANDRE, I. IORDANOFF, J. GIRARDOT
//
//
//
//
//
// 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 3 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, see <http://www.gnu.org/licenses/>.
 
 
 
#ifndef GranOO_SPH_Kernel_hpp
#define GranOO_SPH_Kernel_hpp
 
 
#include "GranOO3/Geom/Vector.hpp"
 
namespace GranOO3
{
  namespace SPH
  {
    namespace Kernel
    {
 
      struct K0 // The default main kernel
      {
        static double k_def(const double &h, const Geom::Vector &d);
        static Geom::Vector gk_def(const double &h, const Geom::Vector &d);
        static double lk_def(const double &h, const Geom::Vector &d);
        static Geom::Vector gk_pres(const double &h, const Geom::Vector &d);
        static double lk_visc(const double &h, const Geom::Vector &d);
      };
 
 
 
      inline double
      K0::k_def(const double &h, const Geom::Vector &d) {
        const double r = d.norm();
        if (r <= h)
          return (315./ (64*M_PI*pow(h,9.)) ) * pow( (h*h) - (r*r), 3.);
        else
          return 0.;
      }
 
      inline Geom::Vector
      K0::gk_def(const double &h, const Geom::Vector &d) {
        double Norm = d.norm();
 
        if (Norm<=h)
          return -945*d*pow((h*h-pow(Norm,2)),2)/(32*M_PI*pow(h,9));
 
        return Geom::Vector(0., 0., 0);
      }
 
      inline double
      K0::lk_def(const double &h, const Geom::Vector &d) {
        double Norm = d.norm();
        if (Norm<=h)
          return -945*(h*h-pow(Norm,2))*(3*h*h-7*pow(Norm,2))/(32*M_PI*pow(h,9));
        else
          return 0.;
      }
 
      inline Geom::Vector
      K0::gk_pres(const double &h, const Geom::Vector &d) {
        const double r = d.norm();
        if (r > 0. && r <= h)
          return (-45. / (M_PI*pow(h,6.))) * (d/r) * pow(h-r,2.);
 
        return Geom::Vector(0., 0., 0.);
      }
 
      inline double
      K0::lk_visc(const double &h, const Geom::Vector &d) {
        double r = d.norm();
        if (r <= h)
          return (45. / (M_PI*pow(h,6)) ) * (h -r);
        else
          return 0.;
      }
 
    }
  }
}
 
 
 
 
 
 
#endif