cngf-pf

continuum model for granular flows with pore-pressure dynamics (renamed from 1d_fd_simple_shear)
git clone git://src.adamsgaard.dk/cngf-pf # fast
git clone https://src.adamsgaard.dk/cngf-pf.git # slow
Log | Files | Refs | README | LICENSE Back to index

commit 838d87537c8aae08d8ce35f4ec6db0dfddfb67f7
parent 3e2725c9388be5990e38846b37e38322d66a0e17
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 16 Feb 2021 15:52:14 +0100

simulation.c: rename timestep variable in calculation function

Diffstat:
Msimulation.c | 12++++++------
1 file changed, 6 insertions(+), 6 deletions(-)

diff --git a/simulation.c b/simulation.c @@ -908,7 +908,7 @@ find_flux(const struct simulation *sim) void set_coupled_fluid_transient_timestep(struct simulation *sim, const double safety) { - double max_gamma_dot, mu, xi, max_I, t_val; + double max_gamma_dot, mu, xi, max_I, dt; /* max expected strain rate */ max_gamma_dot = 1.0/sim->d; @@ -936,10 +936,10 @@ set_coupled_fluid_transient_timestep(struct simulation *sim, const double safety sim->sigma_n[sim->nz - 1] - sim->p_f_mod_ampl, sim->rho_s); - t_val = xi * (sim->alpha + sim->phi[sim->nz - 1] * sim->beta_f) * sim->mu_f - / (sim->phi[sim->nz - 1] * sim->phi[sim->nz - 1] - * sim->phi[sim->nz - 1] * sim->L_z * max_I); + dt = xi * (sim->alpha + sim->phi[sim->nz - 1] * sim->beta_f) * sim->mu_f + / (sim->phi[sim->nz - 1] * sim->phi[sim->nz - 1] + * sim->phi[sim->nz - 1] * sim->L_z * max_I); - if (sim->dt > safety * t_val) - sim->dt = safety * t_val; + if (sim->dt > safety * dt) + sim->dt = safety * dt; }