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> 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;
213 StokesCosine (
const double viscosity = 1.,
const double reaction = 0.);
218 void set_parameters (
const double viscosity,
const double reaction);
222 std::vector<std::vector<double> > &values)
const;
226 std::vector<std::vector<double> > &values)
const;
251 std::vector<std::vector<double> > &values)
const;
253 std::vector<std::vector<
Tensor<1,2> > > &gradients)
const;
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);
309 std::vector<std::vector<double> > &values)
const;
311 std::vector<std::vector<
Tensor<1,2> > > &gradients)
const;
313 std::vector<std::vector<double> > &values)
const;
316 double lambda ()
const;
318 const double Reynolds;
326 DEAL_II_NAMESPACE_CLOSE
std::vector< std::vector< double > > aux_values
const double lm
Auxiliary variable 1-lambda.
virtual void vector_value(const Point< dim > &points, Vector< double > &value) const
const double lp
Auxiliary variable 1+lambda.
static const double lambda
The exponent of the radius.
virtual void vector_gradients(const std::vector< Point< dim > > &points, std::vector< std::vector< Tensor< 1, dim > > > &gradients) const =0
double viscosity
The viscosity.
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
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.
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
virtual void vector_values(const std::vector< Point< dim > > &points, std::vector< std::vector< double > > &values) const =0
void pressure_adjustment(double p)
virtual double value(const Point< dim > &points, const unsigned int component) const