brinkman/parameters.h

Go to the documentation of this file.

/**********************************************************************
 *  Copyright (C) 2011 - 2014 by the authors
 *  Distributed under the MIT License
 *
 * See the files AUTHORS and LICENSE in the project root directory
 *
 **********************************************************************/

#ifndef __brinkman_parameters_h
#define __brinkman_parameters_h

#include <deal.II/base/logstream.h>
#include <deal.II/base/subscriptor.h>

namespace Brinkman
{
class Parameters : public dealii::Subscriptor
{
public:
  Parameters();

  Parameters(double resistance);
  std::vector<double> viscosity;
  std::vector<double> resistance;
  std::vector<double> graddiv_stabilization;

  double saffman;
};

inline Parameters::Parameters()
  : saffman(.1)
{
}

inline Parameters::Parameters(double res)
  : viscosity(2, 0.)
  , resistance(2, 0.)
  , graddiv_stabilization(2, 1.)
  , saffman(.1)
{
  dealii::deallog << "Brinkman " << res << std::endl;
  viscosity[0] = 1.;
  resistance[1] = res;
}
}

#endif