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

The Oldroyd problem on the contraction benchmark – boundary conditions.

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.
tensor operator()(const point &x) const