brinkman-solutions.hpp

Go to the documentation of this file.

#ifndef BRINKMAN_SOLUTIONS_HPP
#define BRINKMAN_SOLUTIONS_HPP
 
#include <Eigen/Dense>
 
#include "MeshBuilder2D.hpp"
 
#include "ns-solutions.hpp"
 
namespace HArDCore2D
{
  namespace BrinkmanSolutions
  {   
    //------------------------------------------------------------------------------
    // Interface for exact solutions
 
    struct IBrinkmanExactSolution : public NavierStokesSolutions::IExactSolution {
      IBrinkmanExactSolution(std::string name, double viscosity, double permeability);
      virtual ~IBrinkmanExactSolution();
      
      inline std::string name() const {
        return m_name;
      }
 
      inline double viscosity() const {
        return m_viscosity;
      }
 
      inline double permeability() const {
    return m_permeability;
      }
    
      virtual Eigen::Vector2d momentumForcingTerm(const Eigen::Vector2d & x, const double & t = 0) const = 0;
      virtual double compressibilityForcingTerm(const Eigen::Vector2d & x, const double & t = 0) const = 0;
      virtual Eigen::Vector2d velocity(const Eigen::Vector2d & x, const double & t = 0) const = 0;
      virtual double pressure(const Eigen::Vector2d & x, const double & t = 0) const = 0;
      virtual Eigen::Matrix2d velocityGradient(const Eigen::Vector2d & x, const double & t = 0) const = 0;
      virtual Eigen::Vector2d pressureGradient(const Eigen::Vector2d & x, const double & t = 0) const = 0;
 
    protected:
      double m_permeability;
    };
  
    //------------------------------------------------------------------------------
    // Linear velocity
 
    struct LinearVelocity : public IBrinkmanExactSolution {
      LinearVelocity();
      
      Eigen::Vector2d momentumForcingTerm(const Eigen::Vector2d & x, const double & t = 0) const;
      double compressibilityForcingTerm(const Eigen::Vector2d & x, const double & t = 0) const;    
      Eigen::Vector2d velocity(const Eigen::Vector2d & x, const double & t = 0) const;   
      double pressure(const Eigen::Vector2d & x, const double & t = 0) const;
      Eigen::Matrix2d velocityGradient(const Eigen::Vector2d & x, const double & t = 0) const;    
      Eigen::Vector2d pressureGradient(const Eigen::Vector2d & x, const double & t = 0) const;
    };
    
    //------------------------------------------------------------------------------
    // Linear pressure
 
    struct LinearPressure : public IBrinkmanExactSolution {
      LinearPressure();
      
      Eigen::Vector2d momentumForcingTerm(const Eigen::Vector2d & x, const double & t = 0) const;
      double compressibilityForcingTerm(const Eigen::Vector2d & x, const double & t = 0) const;    
      Eigen::Vector2d velocity(const Eigen::Vector2d & x, const double & t = 0) const;   
      double pressure(const Eigen::Vector2d & x, const double & t = 0) const;
      Eigen::Matrix2d velocityGradient(const Eigen::Vector2d & x, const double & t = 0) const;    
      Eigen::Vector2d pressureGradient(const Eigen::Vector2d & x, const double & t = 0) const;
    };
 
    //------------------------------------------------------------------------------
    // Quadratic velocity
 
    struct QuadraticVelocity : public IBrinkmanExactSolution {
      QuadraticVelocity();
      
      Eigen::Vector2d momentumForcingTerm(const Eigen::Vector2d & x, const double & t = 0) const;
      double compressibilityForcingTerm(const Eigen::Vector2d & x, const double & t = 0) const;    
      Eigen::Vector2d velocity(const Eigen::Vector2d & x, const double & t = 0) const;   
      double pressure(const Eigen::Vector2d & x, const double & t = 0) const;
      Eigen::Matrix2d velocityGradient(const Eigen::Vector2d & x, const double & t = 0) const;    
      Eigen::Vector2d pressureGradient(const Eigen::Vector2d & x, const double & t = 0) const;
    };
 
    //------------------------------------------------------------------------------
    // Trigonometric
 
    struct Trigonometric : public IBrinkmanExactSolution {
      Trigonometric(double viscosity = 1., double permeability = 1.);
    
      Eigen::Vector2d momentumForcingTerm(const Eigen::Vector2d & x, const double & t = 0) const;
      double compressibilityForcingTerm(const Eigen::Vector2d & x, const double & t = 0) const;    
      Eigen::Vector2d velocity(const Eigen::Vector2d & x, const double & t = 0) const;   
      double pressure(const Eigen::Vector2d & x, const double & t = 0) const;
      Eigen::Matrix2d velocityGradient(const Eigen::Vector2d & x, const double & t = 0) const;    
      Eigen::Vector2d pressureGradient(const Eigen::Vector2d & x, const double & t = 0) const;
    };
 
  } // namespace BrinkmanSolutions
 
} // namespace HArDCore2D
#endif