Commit 3b64a758 authored by Johannes Blaschke's avatar Johannes Blaschke
Browse files

tuning squirmer/pressure coupling

parent 60430edf
......@@ -453,8 +453,7 @@ void particle_hit_squirmer(particle* p, squirmer* s, const double dt, const int
#ifdef ENABLE_PFIELD
if(config.pressure == 1) {
// vmuladd_ip(& s->pressure_f, -config.pressure_k * p->p / sconfig.radius, & rs);
vmuladd_ip(& s->pressure_f, -p->p / sconfig.radius, & rs);
vmuladd_ip(& s->pressure_f, -config.pressure_k / (config.dt*config.ppercell) * p->p / sconfig.radius, & rs);
}
#endif
}
......
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