Rheolef  7.1
an efficient C++ finite element environment
oldroyd_contraction.h
Go to the documentation of this file.
1 #include "contraction.h"
27  struct tau_upstream: base {
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 };
contraction.h
The contraction geometry: boundary conditions.
tensor
see the tensor page for the full documentation
oldroyd_contraction
Definition: oldroyd_contraction.h:26
oldroyd_contraction::tau_upstream::tau_upstream
tau_upstream(geo omega, Float We1, Float alpha1)
Definition: oldroyd_contraction.h:28
contraction
Definition: contraction.h:25
contraction::base::c
Float c
Definition: contraction.h:33
oldroyd_contraction::tau_upstream::operator()
tensor operator()(const point &x) const
Definition: oldroyd_contraction.h:30
oldroyd_contraction::tau_upstream::We
Float We
Definition: oldroyd_contraction.h:38
oldroyd_contraction::tau_upstream::alpha
Float alpha
Definition: oldroyd_contraction.h:38
Float
see the Float page for the full documentation
point
see the point page for the full documentation
oldroyd_contraction::tau_upstream
Definition: oldroyd_contraction.h:27
geo
see the geo page for the full documentation
contraction::base::umax
Float umax
Definition: contraction.h:33