Commit 6e1289b2 authored by Matthieu Constant's avatar Matthieu Constant
Browse files

all terms

parent c7b86890
......@@ -558,10 +558,10 @@ static void fluid_problem_f(FluidProblem *problem, double *f0, double *f1, doubl
double R = dudt + u[i]*drhodt/rho + dp[i]/rho + u[i]*divu/c + utau[i]/c + uugradrho[i]/rho + (i==1 ? -c*problem->g : 0.);
f0[U+i] = rho*dudt + u[i]*drhodt + (i==1 ? -c*rho*problem->g : 0.);
for (int j = 0; j < D; ++j) {
f1[(U+i)*D+j] = -rho*u[i]*u[j]/c + mu*(tau[i][j]+tau[j][i]) + (i==j ? -p : 0) + R*u[j]*rho*taup*0;
f1[(U+i)*D+j] = -rho*u[i]*u[j]/c + mu*(tau[i][j]+tau[j][i]) + (i==j ? -p : 0) + R*u[j]*rho*taup;
}
f1[(U+i)*D+i] += 0*divu*tauc*rho;
f1[P*D+i] = -u[i] + 1e-6*(dp[i]+(i==1 ? -rho*problem->g : 0.));
f1[(U+i)*D+i] += divu*tauc*rho;
f1[P*D+i] = -u[i] + taup*R;
if(problem->n_fluids == 2){
f1[P*D+i] += taua*(divu+(c-cold)/dt)*u[i];
f1[A*D+i] = -u[i]*a + taua*Ra*u[i] + taup*R*a + da[i]*(1e-7+extraa);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment