Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
oldroyd_cavity.h

The Oldroyd problem on the driven cavity benchmark – boundary conditions.

The Oldroyd problem on the driven cavity benchmark – boundary conditions

#include "cavity.h"
static Float u0_top (const point& x) { return 16*sqr(x[0])*sqr(1-x[0]); }
static field velocity_field (const space& Xh) {
geo omega = Xh.get_geo();
space Wh (omega["top"], Xh.get_approx());
field uh (Xh, 0.);
uh[0]["top"] = lazy_interpolate (Wh, u0_top);
return uh;
}
struct tau_upstream {
tau_upstream (geo, Float, Float) {}
tensor operator() (const point& x) const { return tensor(); }
};
};
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
tensor operator()(const point &x) const
static Float u0_top(const point &x)
static field velocity_field(const space &Xh)