round_cells.hpp

Go to the documentation of this file.

#ifndef ROUND_CELLS_HPP
#define ROUND_CELLS_HPP 1

#include "../point_motion.hpp"
#include "../../common/equation_of_state.hpp"
#include "../OuterBoundary.hpp"
#include "../geometric_outer_boundaries/PeriodicBox.hpp"

class RoundCells: public PointMotion
{
public:

        RoundCells(const PointMotion& pm, const EquationOfState& eos, OuterBoundary const& outer,
                double chi = 0.15, double eta = 0.02, bool cold = false, double cold_speed = 0.15);


        RoundCells(const PointMotion& pm, const EquationOfState& eos, double chi = 0.15, double eta = 0.02,bool cold = false
                , double cold_speed = 0.15);

  vector<Vector2D> operator()(const Tessellation& tess,const vector<ComputationalCell>& cells,double time,
          TracerStickerNames const& tracerstickernames) const;

  vector<Vector2D> ApplyFix(Tessellation const& tess, vector<ComputationalCell> const& cells, double time,
          double dt, vector<Vector2D> const& velocities, TracerStickerNames const& tracerstickernames)const;
private:

  Vector2D calc_dw(size_t i, const Tessellation& tess,const vector<ComputationalCell>& cells,
          TracerStickerNames const& tracerstickernames) const;

  Vector2D calc_dw(size_t i, const Tessellation& tess, double dt,vector<ComputationalCell> const& cells,
          TracerStickerNames const& tracerstickernames)const;

  const PointMotion& pm_;
  const EquationOfState& eos_;
  PeriodicBox pouter_;
  OuterBoundary const& outer_;
  const double chi_;
  const double eta_;
  const bool cold_;
  const double cold_speed_;
};

#endif // ROUND_CELLS_HPP