Rheolef  7.1
an efficient C++ finite element environment
oldroyd_cavity.h

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"] = interpolate (Wh, u0_top);
return uh;
}
struct tau_upstream {
tensor operator() (const point& x) const { return tensor(); }
};
};
oldroyd_cavity::tau_upstream::operator()
tensor operator()(const point &x) const
Definition: oldroyd_cavity.h:37
tensor
see the tensor page for the full documentation
field
see the field page for the full documentation
space
see the space page for the full documentation
oldroyd_cavity
Definition: oldroyd_cavity.h:26
rheolef::space_constant::tensor
@ tensor
Definition: space_constant.h:138
rheolef::interpolate
field_basic< T, M > interpolate(const space_basic< T, M > &V2h, const field_basic< T, M > &u1h)
see the interpolate page for the full documentation
Definition: interpolate.cc:233
oldroyd_cavity::u0_top
static Float u0_top(const point &x)
Definition: oldroyd_cavity.h:27
cavity
Definition: cavity.h:25
Float
see the Float page for the full documentation
cavity.h
The driven cavity benchmark: 2D boundary conditions.
point
see the point page for the full documentation
oldroyd_cavity::velocity_field
static field velocity_field(const space &Xh)
Definition: oldroyd_cavity.h:28
geo
see the geo page for the full documentation
oldroyd_cavity::tau_upstream::tau_upstream
tau_upstream(geo, Float, Float)
Definition: oldroyd_cavity.h:36