16 #ifndef dealii__flow_function_h 17 #define dealii__flow_function_h 20 #include <deal.II/base/config.h> 21 #include <deal.II/base/function.h> 22 #include <deal.II/base/point.h> 23 #include <deal.II/base/thread_management.h> 25 DEAL_II_NAMESPACE_OPEN
73 std::vector<std::vector<double> > &values)
const = 0;
80 std::vector<std::vector<
Tensor<1,dim> > > &gradients)
const = 0;
90 std::vector<std::vector<double> > &values)
const = 0;
93 virtual double value (
const Point<dim> &points,
const unsigned int component)
const;
96 virtual void vector_gradient_list (
const std::vector<
Point<dim> > &points,
104 std::size_t memory_consumption ()
const;
151 std::vector<std::vector<double> > &values)
const;
155 std::vector<std::vector<double> > &values)
const;
159 const double Reynolds;
185 StokesCosine (
const double viscosity = 1.,
const double reaction = 0.);
189 void set_parameters (
const double viscosity,
const double reaction);
193 std::vector<std::vector<double> > &values)
const;
197 std::vector<std::vector<double> > &values)
const;
222 std::vector<std::vector<double> > &values)
const;
224 std::vector<std::vector<
Tensor<1,2> > > &gradients)
const;
226 std::vector<std::vector<double> > &values)
const;
229 double Psi(
double phi)
const;
231 double Psi_1(
double phi)
const;
233 double Psi_2(
double phi)
const;
235 double Psi_3(
double phi)
const;
237 double Psi_4(
double phi)
const;
268 Kovasznay (
const double Re,
bool Stokes =
false);
272 std::vector<std::vector<double> > &values)
const;
274 std::vector<std::vector<
Tensor<1,2> > > &gradients)
const;
276 std::vector<std::vector<double> > &values)
const;
279 double lambda ()
const;
281 const double Reynolds;
289 DEAL_II_NAMESPACE_CLOSE
std::vector< std::vector< double > > aux_values
virtual double value(const Point< dim > &points, const unsigned int component) const
const double lm
Auxiliary variable 1-lambda.
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_value(const Point< dim > &points, Vector< double > &value) const
std::vector< std::vector< Tensor< 1, dim > > > aux_gradients
double reaction
The reaction parameter.
virtual void vector_laplacian_list(const std::vector< Point< dim > > &points, std::vector< Vector< double > > &values) const
virtual void vector_value_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)