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
"
struct
oldroyd_cavity
:
cavity
{
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 {
tau_upstream
(
geo
,
Float
,
Float
) {}
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
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