inertia_upw.icc

The inertia term of the Navier-Stokes equation with the discontinuous Galerkin method – upwinding variant

#include "sgn.icc"
form inertia_upw (field w, trial u, test v,
integrate_option iopt = integrate_option())
{
return integrate ("internal_sides",
0.5*abs(dot(average(w),normal()))*dot(jump(u),jump(v)));
}
form d_inertia_upw (field w, trial dw, field u, test v,
integrate_option iopt = integrate_option())
{
return integrate ("internal_sides",
0.5*compose (sgn, dot(average(w),normal()))
*dot(average(dw),normal())*dot(jump(u),jump(v)));
}
form
see the form page for the full documentation
field
see the field page for the full documentation
rheolef::normal
details::field_expr_v2_nonlinear_terminal_function< details::normal_pseudo_function< Float > > normal()
normal: see the expression page for the full documentation
Definition: field_expr_terminal.h:439
rheolef::integrate
std::enable_if< details::is_field_expr_v2_nonlinear_arg< Expr >::value &&! is_undeterminated< Result >::value, Result >::type integrate(const geo_basic< T, M > &omega, const Expr &expr, const integrate_option &iopt, Result dummy=Result())
see the integrate page for the full documentation
Definition: integrate.h:202
rheolef::details::compose
class rheolef::details::field_expr_v2_nonlinear_node_unary compose
sgn.icc
The sign function.
test
see the test page for the full documentation
u
Definition: leveque.h:25
inertia_upw
form inertia_upw(field w, trial u, test v, integrate_option iopt=integrate_option())
Definition: inertia_upw.icc:26
trial
see the test page for the full documentation
d_inertia_upw
form d_inertia_upw(field w, trial dw, field u, test v, integrate_option iopt=integrate_option())
Definition: inertia_upw.icc:32
sgn
Float sgn(Float x)
Definition: sgn.icc:25
rheolef::details::dot
rheolef::details::is_vec dot