Hi all,
for maxadvec we find in the code
dustvelocity.f90: maxadvec=maxadvec+p%advec_uud(:,k) ! k runs over the species!
heatflux.f90: maxadvec = maxadvec + c_spitzer/dxmin_pencil
hydro.f90: maxadvec=maxadvec+p%advec_uu
neutralvelocity.f90: maxadvec=maxadvec+p%advec_uun
I can't see a reason why there is summing-up: shouldn't the maximum of these advective velocities matter for the timestep instead of the sum as they all act on different quantities?
In particular for dust with a high number of variables, the result could be absurd.
Best,
Matthias
Hi Axel,
I'm not questioning either of the plus signs in Uadvect + sqrt(cs^2+vA^2). What I refer to is that Uadvect itself is already a sum, e.g. in a two fluid setup with dust:
velocity(fluid1) + velocity(fluid2) + \sum_i velocity{dust species i)
This can't be meaningful as each of these velocities acts on a different field.
Best,
Matthias
I think it is less about problems, but more that due to the sum, the timestep becomes unnecessarily small.
Indeed, I could see this for samples where Uadvec matters for the timestep.
M.