Commit 8d611dca authored by Jonathan Lambrechts's avatar Jonathan Lambrechts
Browse files

no V dp in force on fluid

parent da553552
......@@ -151,7 +151,7 @@ static void particle_force_f(FluidProblem *problem, double *f0, double *sol, dou
for (int j = 0; j < D; ++j){
if(!in_derivative)
problem->particle_force[ip*D+j] = (-gamma*Du[j]-vol*dp[j]);
f0[U+j] = -gamma*Du[j]-vol*dp[j];
f0[U+j] = -gamma*Du[j]; //-vol*dp[j];
}
if(!in_derivative)
problem->particle_force[ip*D+1] += g*mass;
......
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