Rheolef  7.2
an efficient C++ finite element environment
oldroyd_cavity.h
Go to the documentation of this file.
1 #include "cavity.h"
27  static Float u0_top (const point& x) { return 16*sqr(x[0])*sqr(1-x[0]); }
28  static field velocity_field (const space& Xh) {
29  geo omega = Xh.get_geo();
30  space Wh (omega["top"], Xh.get_approx());
31  field uh (Xh, 0.);
32  uh[0]["top"] = lazy_interpolate (Wh, u0_top);
33  return uh;
34  }
35  struct tau_upstream {
37  tensor operator() (const point& x) const { return tensor(); }
38  };
39 };
The driven cavity benchmark: 2D boundary conditions.
see the Float page for the full documentation
see the field page for the full documentation
see the geo page for the full documentation
see the point page for the full documentation
see the space page for the full documentation
see the tensor page for the full documentation
field_basic< T, M > lazy_interpolate(const space_basic< T, M > &X2h, const field_basic< T, M > &u1h)
see the interpolate page for the full documentation
Definition: field.h:871
Definition: cavity.h:25
tensor operator()(const point &x) const
tau_upstream(geo, Float, Float)
static Float u0_top(const point &x)
static field velocity_field(const space &Xh)