Rheolef
7.1
an efficient C++ finite element environment
oldroyd_contraction.h
Go to the documentation of this file.
1
#include "
contraction.h
"
26
struct
oldroyd_contraction
:
contraction
{
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
}
38
Float
We
,
alpha
;
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
Definition:
contraction.h:26
contraction::base::umax
Float umax
Definition:
contraction.h:33