Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
oldroyd_contraction.h
Go to the documentation of this file.
1
25#include "contraction.h"
28 tau_upstream (geo omega, Float We1, Float alpha1)
29 : base(omega), We(We1), alpha(alpha1) {}
30 tensor operator() (const point& x) const {
31 tensor tau;
32 Float dot_gamma = - 2*base::umax*x[1]/sqr(base::c);
33 tau(0,0) = 2*alpha*We*sqr(dot_gamma);
34 tau(0,1) = tau(1,0) = alpha*dot_gamma;
35 tau(1,1) = 0;
36 return tau;
37 }
39 };
40};
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
tau_upstream(geo omega, Float We1, Float alpha1)