Reference documentation for deal.II version 8.1.0
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Modules Pages
flow_function.h
1 // ---------------------------------------------------------------------
2 // @f$Id: flow_function.h 30036 2013-07-18 16:55:32Z maier @f$
3 //
4 // Copyright (C) 2007 - 2013 by the deal.II authors
5 //
6 // This file is part of the deal.II library.
7 //
8 // The deal.II library is free software; you can use it, redistribute
9 // it, and/or modify it under the terms of the GNU Lesser General
10 // Public License as published by the Free Software Foundation; either
11 // version 2.1 of the License, or (at your option) any later version.
12 // The full text of the license can be found in the file LICENSE at
13 // the top level of the deal.II distribution.
14 //
15 // ---------------------------------------------------------------------
16 
17 #ifndef __deal2__flow_function_h
18 #define __deal2__flow_function_h
19 
20 
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>
25 
27 
28 namespace Functions
29 {
49  template <int dim>
50  class FlowFunction : public Function<dim>
51  {
52  public:
57  FlowFunction();
58 
62  virtual ~FlowFunction();
63 
70  void pressure_adjustment(double p);
71 
80  virtual void vector_values (const std::vector<Point<dim> > &points,
81  std::vector<std::vector<double> > &values) const = 0;
90  virtual void vector_gradients (const std::vector<Point<dim> > &points,
91  std::vector<std::vector<Tensor<1,dim> > > &gradients) const = 0;
106  virtual void vector_laplacians (const std::vector<Point<dim> > &points,
107  std::vector<std::vector<double> > &values) const = 0;
108 
109  virtual void vector_value (const Point<dim> &points, Vector<double> &value) const;
110  virtual double value (const Point<dim> &points, const unsigned int component) const;
111  virtual void vector_value_list (const std::vector<Point<dim> > &points,
112  std::vector<Vector<double> > &values) const;
113  virtual void vector_gradient_list (const std::vector<Point<dim> > &points,
114  std::vector<std::vector<Tensor<1,dim> > > &gradients) const;
119  virtual void vector_laplacian_list (const std::vector<Point<dim> > &points,
120  std::vector<Vector<double> > &values) const;
121 
122  std::size_t memory_consumption () const;
123 
124  protected:
131 
132  private:
133 
139 
144  mutable std::vector<std::vector<double> > aux_values;
145 
150  mutable std::vector<std::vector<Tensor<1,dim> > > aux_gradients;
151  };
152 
162  template <int dim>
163  class PoisseuilleFlow : public FlowFunction<dim>
164  {
165  public:
172  PoisseuilleFlow<dim> (const double r,
173  const double Re);
174  virtual ~PoisseuilleFlow();
175 
176  virtual void vector_values (const std::vector<Point<dim> > &points,
177  std::vector<std::vector<double> > &values) const;
178  virtual void vector_gradients (const std::vector<Point<dim> > &points,
179  std::vector<std::vector<Tensor<1,dim> > > &gradients) const;
180  virtual void vector_laplacians (const std::vector<Point<dim> > &points,
181  std::vector<std::vector<double> > &values) const;
182 
183  private:
184  const double radius;
185  const double Reynolds;
186  };
187 
188 
202  template <int dim>
203  class StokesCosine :
204  public FlowFunction<dim>
205  {
206  public:
213  StokesCosine (const double viscosity = 1., const double reaction = 0.);
218  void set_parameters (const double viscosity, const double reaction);
219  virtual ~StokesCosine();
220 
221  virtual void vector_values (const std::vector<Point<dim> > &points,
222  std::vector<std::vector<double> > &values) const;
223  virtual void vector_gradients (const std::vector<Point<dim> > &points,
224  std::vector<std::vector<Tensor<1,dim> > > &gradients) const;
225  virtual void vector_laplacians (const std::vector<Point<dim> > &points,
226  std::vector<std::vector<double> > &values) const;
227 
228  private:
230  double viscosity;
232  double reaction;
233  };
234 
235 
245  {
246  public:
249 
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;
256  private:
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;
268  const double omega;
270  static const double lambda;
272  const double coslo;
274  const double lp;
276  const double lm;
277  };
278 
288  class Kovasznay : public FlowFunction<2>
289  {
290  public:
305  Kovasznay (const double Re, bool Stokes = false);
306  virtual ~Kovasznay();
307 
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;
314 
316  double lambda () const;
317  private:
318  const double Reynolds;
319  double lbda;
320  double p_average;
321  const bool stokes;
322  };
323 
324 }
325 
326 DEAL_II_NAMESPACE_CLOSE
327 
328 #endif
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