Loading...
Searching...
No Matches
oldroyd_theta_scheme1.h
Go to the documentation of this file.
template<class P>
oldroyd_theta_scheme<P>::oldroyd_theta_scheme()
: We(0), alpha(8./9), a(1), Re(1), delta_t(0.025), tol(1e-6), max_iter(500),
Th(), Xh(), Qh(), b(), c(), d(), mt(), inv_mt(), mu(), mp(), th(), thb(),
theta(), lambda(), eta(), nu(), c1(), c2(), c3(), c4(), c5(), stokes() {}
template<class P>
void oldroyd_theta_scheme<P>::reset(const geo& omega) {
Th = space (omega, "P1d", "tensor");
Xh = P::velocity_space (omega, "P2");
Qh = space (omega, "P1d");
theta = 1-1/sqrt(2.);
lambda = Re/(theta*delta_t);
eta = ((1 - alpha)*We + theta*delta_t)/(We + theta*delta_t);
nu = 1/((1-2*theta)*delta_t);
c1 = We/(We + theta*delta_t);
c2 = - We*theta*delta_t/(We + theta*delta_t);
c3 = alpha*theta*delta_t/(We + theta*delta_t);
c4 = 1/((1-2*theta)*delta_t) - 1/We;
c5 = alpha/We;
trial u (Xh), tau(Th), p (Qh);
test v (Xh), xi (Th), q (Qh);
mt = integrate (ddot(tau,xi));
mu = integrate (dot(u,v));
mp = integrate (p*q);
integrate_option iopt;
iopt.invert = true;
inv_mt = integrate (ddot(tau,xi), iopt);
b = integrate (-ddot(tau,D(v)));
c = integrate (lambda*dot(u,v) + 2*eta*ddot(D(u),D(v)));
d = integrate (-div(u)*q);
stokes = problem_mixed (c, d);
stokes.set_metric (mp);
}
geo
see the geo page for the full documentation
problem_mixed
see the problem_mixed page for the full documentation
space
see the space page for the full documentation
test
see the test page for the full documentation
trial
see the test page for the full documentation
u
point u(const point &x)
Definition diffusion_transport_tensor_exact.icc:25
oldroyd_theta_scheme::c4
Float c4
Definition oldroyd_theta_scheme.h:49
oldroyd_theta_scheme::inv_mt
form inv_mt
Definition oldroyd_theta_scheme.h:46
oldroyd_theta_scheme::th
form th
Definition oldroyd_theta_scheme.h:47
oldroyd_theta_scheme::theta
Float theta
Definition oldroyd_theta_scheme.h:49
oldroyd_theta_scheme::a
Float a
Definition oldroyd_theta_scheme.h:42
oldroyd_theta_scheme::mp
form mp
Definition oldroyd_theta_scheme.h:46
oldroyd_theta_scheme::c3
Float c3
Definition oldroyd_theta_scheme.h:49
oldroyd_theta_scheme::We
Float We
Definition oldroyd_theta_scheme.h:42
oldroyd_theta_scheme::mt
form mt
Definition oldroyd_theta_scheme.h:46
oldroyd_theta_scheme::delta_t
Float delta_t
Definition oldroyd_theta_scheme.h:42
oldroyd_theta_scheme::Xh
space Xh
Definition oldroyd_theta_scheme.h:45
oldroyd_theta_scheme::c2
Float c2
Definition oldroyd_theta_scheme.h:49
oldroyd_theta_scheme::eta
Float eta
Definition oldroyd_theta_scheme.h:49
oldroyd_theta_scheme::oldroyd_theta_scheme
oldroyd_theta_scheme()
Definition oldroyd_theta_scheme1.h:26
oldroyd_theta_scheme::stokes
problem_mixed stokes
Definition oldroyd_theta_scheme.h:50
oldroyd_theta_scheme::thb
field thb
Definition oldroyd_theta_scheme.h:48
oldroyd_theta_scheme::lambda
Float lambda
Definition oldroyd_theta_scheme.h:49
oldroyd_theta_scheme::alpha
Float alpha
Definition oldroyd_theta_scheme.h:42
oldroyd_theta_scheme::mu
form mu
Definition oldroyd_theta_scheme.h:46
oldroyd_theta_scheme::Re
Float Re
Definition oldroyd_theta_scheme.h:42
oldroyd_theta_scheme::d
form d
Definition oldroyd_theta_scheme.h:46
oldroyd_theta_scheme::tol
Float tol
Definition oldroyd_theta_scheme.h:42
oldroyd_theta_scheme::nu
Float nu
Definition oldroyd_theta_scheme.h:49
oldroyd_theta_scheme::c1
Float c1
Definition oldroyd_theta_scheme.h:49
oldroyd_theta_scheme::b
form b
Definition oldroyd_theta_scheme.h:46
oldroyd_theta_scheme::Th
space Th
Definition oldroyd_theta_scheme.h:45
oldroyd_theta_scheme::c
form c
Definition oldroyd_theta_scheme.h:46
oldroyd_theta_scheme::max_iter
size_t max_iter
Definition oldroyd_theta_scheme.h:43
oldroyd_theta_scheme::reset
void reset(const geo &omega)
Definition oldroyd_theta_scheme1.h:31
oldroyd_theta_scheme::Qh
space Qh
Definition oldroyd_theta_scheme.h:45
oldroyd_theta_scheme::c5
Float c5
Definition oldroyd_theta_scheme.h:49
p
Definition sphere.icc:25
u
Definition leveque.h:25