1d_fd_simple_shear_transient

commit 9d72e8471f5ca61b81ac0f9951116c3057b01467
parent 05b6c6a8eb96032d53f60ae53158a3870dc7cc5d
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Mon, 26 Aug 2019 13:48:23 +0200

Fix transient stress behavior when limiting shear velocity

Diffstat:
M
main.c
 | 
4
+++-

1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/main.c b/main.c
@@ -107,7 +107,7 @@ main(int argc, char* argv[])
 	struct simulation sim;
 	const char* optstring;
 	unsigned long iter, stressiter;
-	double new_phi, new_k, filetimeclock, res_norm;
+	double new_phi, new_k, filetimeclock, res_norm, mu_wall_orig;
 
 	/* load with default values */
 	sim = init_sim();
@@ -291,6 +291,7 @@ main(int argc, char* argv[])
 
 	filetimeclock = 0.0;
 	iter = 0;
+	mu_wall_orig = sim.mu_wall;
 	while (sim.t <= sim.t_end) {
 
 		stressiter = 0;
@@ -337,6 +338,7 @@ main(int argc, char* argv[])
 		} while ((!isnan(sim.v_x_fix) || !isnan(sim.v_x_limit))
 		         && fabs(res_norm) > RTOL);
 
+		sim.mu_wall = mu_wall_orig;
 		sim.t += sim.dt;
 		filetimeclock += sim.dt;
 		iter++;