search

ECP-copa / ExaMPM

Material point method proxy application based on Cabana.
BSD 3-Clause "New" or "Revised" License
9 stars 11 forks source link

CFL condition #43

Closed kwitaechong closed 1 year ago

kwitaechong commented 1 year ago

add CFL condition based on wave speed c = sqrt( Bulk_modulus / rho ) and max velocity

streeve commented 1 year ago

@kwitaechong I updated quite a bit since the solver still assumed a fixed step size. Can you test before I approve and merge?

streeve commented 1 year ago

@abisner can you take a quick look as well?

kwitaechong commented 1 year ago

using delta_t = 0.001 is already sufficiently low and thus create output file in a uniform time interval. On the other hand, using delta_t = 0.1 violates cfl condition and the new computed dt = 0.00134849, 0.00129002.. and so on was used. Looks working fine