Rheolef  7.2
an efficient C++ finite element environment
oldroyd_contraction.h

The Oldroyd problem on the contraction benchmark – boundary conditions

#include "contraction.h"
struct tau_upstream: base {
tau_upstream (geo omega, Float We1, Float alpha1)
: base(omega), We(We1), alpha(alpha1) {}
tensor operator() (const point& x) const {
tensor tau;
Float dot_gamma = - 2*base::umax*x[1]/sqr(base::c);
tau(0,0) = 2*alpha*We*sqr(dot_gamma);
tau(0,1) = tau(1,0) = alpha*dot_gamma;
tau(1,1) = 0;
return tau;
}
};
};
see the Float page for the full documentation
see the geo page for the full documentation
see the point page for the full documentation
see the tensor page for the full documentation
The contraction geometry: boundary conditions.
base(geo omega)
Definition: contraction.h:27
tensor operator()(const point &x) const
tau_upstream(geo omega, Float We1, Float alpha1)