17 #ifndef __deal2__flow_function_h
18 #define __deal2__flow_function_h
21 #include <deal.II/base/config.h>
22 #include <deal.II/base/function.h>
23 #include <deal.II/base/point.h>
24 #include <deal.II/base/thread_management.h>
26 DEAL_II_NAMESPACE_OPEN
81 std::vector<std::vector<double> > &values)
const = 0;
91 std::vector<std::vector<
Tensor<1,dim> > > &gradients)
const = 0;
107 std::vector<std::vector<double> > &values)
const = 0;
110 virtual double value (
const Point<dim> &points,
const unsigned int component)
const;
122 std::size_t memory_consumption ()
const;
177 std::vector<std::vector<double> > &values)
const;
181 std::vector<std::vector<double> > &values)
const;
185 const double Reynolds;
222 std::vector<std::vector<double> > &values)
const;
226 std::vector<std::vector<double> > &values)
const;
250 virtual void vector_values (
const std::vector<
Point<2> > &points,
251 std::vector<std::vector<double> > &values)
const;
252 virtual void vector_gradients (
const std::vector<
Point<2> > &points,
253 std::vector<std::vector<
Tensor<1,2> > > &gradients)
const;
254 virtual void vector_laplacians (
const std::vector<
Point<2> > &points,
255 std::vector<std::vector<double> > &values)
const;
258 double Psi(
double phi)
const;
260 double Psi_1(
double phi)
const;
262 double Psi_2(
double phi)
const;
264 double Psi_3(
double phi)
const;
266 double Psi_4(
double phi)
const;
305 Kovasznay (
const double Re,
bool Stokes =
false);
308 virtual void vector_values (
const std::vector<
Point<2> > &points,
309 std::vector<std::vector<double> > &values)
const;
310 virtual void vector_gradients (
const std::vector<
Point<2> > &points,
311 std::vector<std::vector<
Tensor<1,2> > > &gradients)
const;
312 virtual void vector_laplacians (
const std::vector<
Point<2> > &points,
313 std::vector<std::vector<double> > &values)
const;
318 const double Reynolds;
326 DEAL_II_NAMESPACE_CLOSE
std::vector< std::vector< double > > aux_values
const double lm
Auxiliary variable 1-lambda.
double lambda() const
The value of lambda.
virtual void vector_values(const std::vector< Point< dim > > &points, std::vector< std::vector< double > > &values) const
virtual void vector_value(const Point< dim > &points, Vector< double > &value) const
const double lp
Auxiliary variable 1+lambda.
virtual void vector_values(const std::vector< Point< dim > > &points, std::vector< std::vector< double > > &values) const
static const double lambda
The exponent of the radius.
StokesCosine(const double viscosity=1., const double reaction=0.)
virtual void vector_gradients(const std::vector< Point< dim > > &points, std::vector< std::vector< Tensor< 1, dim > > > &gradients) const =0
void set_parameters(const double viscosity, const double reaction)
double viscosity
The viscosity.
virtual void vector_laplacians(const std::vector< Point< dim > > &points, std::vector< std::vector< double > > &values) const
double Psi_3(double phi) const
The 3rd derivative of Psi()
virtual void vector_gradient_list(const std::vector< Point< dim > > &points, std::vector< std::vector< Tensor< 1, dim > > > &gradients) const
std::vector< std::vector< Tensor< 1, dim > > > aux_gradients
virtual void vector_gradients(const std::vector< Point< dim > > &points, std::vector< std::vector< Tensor< 1, dim > > > &gradients) const
double reaction
The reaction parameter.
virtual void vector_value_list(const std::vector< Point< dim > > &points, std::vector< Vector< double > > &values) const
virtual void vector_laplacian_list(const std::vector< Point< dim > > &points, std::vector< Vector< double > > &values) const
const double omega
The angle of the reentrant corner.
double Psi_2(double phi) const
The 2nd derivative of Psi()
virtual void vector_gradients(const std::vector< Point< dim > > &points, std::vector< std::vector< Tensor< 1, dim > > > &gradients) const
const double coslo
Cosine of lambda times omega.
virtual void vector_laplacians(const std::vector< Point< dim > > &points, std::vector< std::vector< double > > &values) const =0
StokesLSingularity()
Constructor setting upsome data.
virtual void vector_values(const std::vector< Point< dim > > &points, std::vector< std::vector< double > > &values) const =0
double Psi(double phi) const
The auxiliary function Psi.
Kovasznay(const double Re, bool Stokes=false)
void pressure_adjustment(double p)
double Psi_1(double phi) const
The derivative of Psi()
virtual void vector_laplacians(const std::vector< Point< dim > > &points, std::vector< std::vector< double > > &values) const
double Psi_4(double phi) const
The 4th derivative of Psi()
virtual double value(const Point< dim > &points, const unsigned int component) const