Commit aabd1688 authored by Jonathan Lambrechts's avatar Jonathan Lambrechts
Browse files

Merge branch 'drag_in_stab' into 'master'

add complete drag_in_stab factor

See merge request !7
parents 58ffcefd c6071167
Pipeline #7349 passed with stage
in 2 minutes and 21 seconds
......@@ -437,8 +437,8 @@ static void fluid_problem_f(const FluidProblem *problem, const double *sol, doub
f10[(P*D+i)*n_fields+U+i] += -1;
// PSPG
double pspg = taup/rho;
f1[P*D+i] += pspg*f0[U+i] + problem->stab_param*dp[i];
f11[(P*D+i)*(n_fields*D)+P*D+i] += pspg*f01[(U+i)*n_fields*D+P*D+i]+problem->stab_param;
f1[P*D+i] += pspg*(f0[U+i]+(problem->drag_in_stab?0:(1-c)*dp[i])) + problem->stab_param*dp[i];
f11[(P*D+i)*(n_fields*D)+P*D+i] += pspg*(f01[(U+i)*n_fields*D+P*D+i]+(problem->drag_in_stab?0:1-c))+problem->stab_param;
f10[(P*D+i)*n_fields+U+i] += pspg*f00[(U+i)*n_fields+(U+i)];
for (int j = 0; j < D; ++j) {
f11[(P*D+i)*n_fields*D+(U+i)*D+j] = pspg*f01[(U+i)*n_fields*D+(U+i)*D+j];
......
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