un petit progres

doc/2fluids.tex

@@ -104,8 +104,8 @@ Balance equations resulting from applying Navier-Stokes equations to a control v
 \begin{align*}
 N_f\times\left\lbrace\begin{aligned}
 \rho_k\dpartial{\ub}{t}+\rho_k\nb\cdot\left(\dfrac{\ub\ub}{q_k}\right)\\
-\dpartial{q}{t}+\nb\cdot\ub+\varepsilon\nb\cdot \bm{r}_k\end{aligned}\right.&\begin{aligned}
-&=-q_k\nb p+\nb\cdot\taub +q_k\rho_k \gb\phantom{\dpartial{\ub}{t}}\\
+\dpartial{q_k}{t}+\nb\cdot\ub+\varepsilon\nb\cdot \bm{r}_k\end{aligned}\right.&\begin{aligned}
+&=-q_k\nb p+\nb\cdot(q_k\taub) +q_k\rho_k \gb\phantom{\dpartial{\ub}{t}}\\
 &=0\end{aligned}\\
 1-c-\sum_{k=1}^{N_f}q_k&\begin{aligned}&=0\phantom{\dpartial{q}{t}}\end{aligned}
 \end{align*}

src/fluid_problem.c

@@ -351,7 +351,7 @@ static void fluid_problem_f(FluidProblem *problem, double *f0, double *f1, doubl
         utau[i] += u[j]*tau[i][j];
       }
     }
-    f0[Q] = (q-qold)/dt;
+    f0[Q] = (q-qold)/dt+epsilonp*dq[1]*rho*problem->g;
     f0[P] -= q;
     for (int i = 0; i < D; ++i) {
       //f0[Q] += epsilonp*dq[i]*dp[i];

testcases/separation/rest.py

@@ -65,7 +65,7 @@ p = 6
 strong_boundaries = [
         #("Top",u0,u0,0),
         #("Top",u1,u1,0),
-        ("PtFix",q0,p,0),
+        ("PtFix",q1,p,0),
         #("Bottom",u0,u0,0),
         #("Bottom",u1,u1,0),
         ("Left",u0,u0,0),