47 uh = (1-theta)/theta*uh1 - (1-2*theta)/theta*uh0;
50 field Duh = inv_mt*integrate(ddot(D(uh),xi));
54 update_transport_stress (uh);
56 typename P::tau_upstream tau_up (Th.get_geo(), We, alpha);
57 field lh = integrate (ddot(c4*tau_h1 + 2*c5*D(uh1),xi))
58 + integrate (
"boundary",
59 max(0, -dot(uh,normal()))*ddot(tau_up,xi));
61 transport.solve (
lh, tau_h);
66 typename P::tau_upstream tau_up (Th.get_geo(), We, alpha);
68 auto ma = 0.5*((1-a)*grad(uh) - (1+a)*trans(grad(uh)));
69 auto beta_a = tau*ma + trans(ma)*tau;
70 th = integrate (ddot(grad_h(tau)*uh + beta_a,xi))
71 + integrate (
"boundary", max(0, -dot(uh,normal()))*ddot(tau,xi))
72 + integrate (
"internal_sides",
73 - dot(uh,normal())*ddot(jump(tau),average(xi))
74 + 0.5*abs(dot(uh,normal()))*ddot(jump(tau),jump(xi)));
75 thb = integrate (
"boundary", max(0, -dot(uh,normal()))*ddot(tau_up,xi));