Commit 95bad7bc authored by Matthieu Constant's avatar Matthieu Constant
Browse files

remis ap et bp

parent 4bcac4a0
Pipeline #5512 failed with stage
in 2 minutes and 1 second
......@@ -319,18 +319,17 @@ static void compute_force_fields(FluidProblem *problem, double dt) {
g[1] = problem->g;
problem->gamma[ip] = gamma;
double gammastar = gamma/(1+dt/m*gamma);
/*
for (int id = 0; id < D; ++id) {
ap[id] = gammastar*(vp[id]+dt/m*contactp[id]+dt*g[id]);
}
double bp = -gammastar*V*dt/m;
double cp = -gammastar/c;
*/
/*
for (int id = 0; id < D; ++id) {
ap[id] = 0; //gamma*(vp[id]);
}
double bp = 0;
double cp = -gamma/c;
double cp = -gamma/c;*/
for (int inode = problem->particle_nodes_n[ip]; inode < problem->particle_nodes_n[ip+1]; ++inode) {
int node = problem->particle_nodes[inode];
double w = problem->particle_nodes_weight[inode];
......@@ -420,10 +419,10 @@ static void fluid_problem_f(const FluidProblem *problem, const double *sol, doub
rho*(u[i]-uold[i])/dt
+ c*dp[i]
+ G[i]*rhoreduced*c
+ problem->volume_drag*u[i]*mu- cp*u[i];;
/*- ap[i] - bp*dp[i] */
+ problem->volume_drag*u[i]*mu
- ap[i] - bp*dp[i] - cp*u[i];
f00[(U+i)*n_fields+U+i] = rho/dt + problem->volume_drag*mu -cp;
f01[(U+i)*n_fields*D+P*D+i] = c;//-bp;
f01[(U+i)*n_fields*D+P*D+i] = c-bp;
// advection :
// dj(uj ui /c) = uj/c dj(ui) + ui/c dj(uj) - uj ui /(c*c) dj(c)
// = ui/c * (dj(uj)- uj/c*dj(c)) + uj/c*dj(ui)
......
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