parameter_defaults.h (2660B)
1 #ifndef ONED_FD_SIMPLE_SHEAR_ 2 #define ONED_FD_SIMPLE_SHEAR_ 3 4 #include <math.h> 5 #include <stdio.h> 6 #include "arrays.h" 7 #include "simulation.h" 8 9 #define DEFAULT_SIMULATION_NAME "unnamed_simulation" 10 11 /* Simulation settings */ 12 void 13 init_sim() 14 { 15 snprintf(sim.name, sizeof(sim.name), DEFAULT_SIMULATION_NAME); 16 17 sim.G = 9.81; 18 19 sim.P_wall = 120e3; 20 sim.mu_wall = 0.45; 21 sim.v_x_bot = 0.0; 22 sim.v_x_fix = NAN; 23 sim.v_x_limit = NAN; 24 25 sim.nz = -1; /* cell size equals grain size if negative */ 26 27 /* lower values of A mean that the velocity curve can have sharper curves, 28 * e.g. at the transition from μ ≈ μ_s */ 29 sim.A = 0.40; /* Loose fit to Damsgaard et al 2013 */ 30 31 /* lower values of b mean larger shear velocity for a given stress ratio 32 * above yield */ 33 sim.b = 0.9377; /* Henann and Kamrin 2016 */ 34 35 /* Henann and Kamrin 2016 */ 36 /* sim.mu_s = 0.3819; */ 37 /* sim.C = 0.0; */ 38 39 /* Damsgaard et al 2013 */ 40 sim.mu_s = tan(DEG2RAD(22.0)); 41 sim.C = 0.0; 42 sim.phi = initval(0.25, 1); 43 sim.d = 0.04; /* Damsgaard et al 2013 */ 44 45 sim.phi_min = 0.2; 46 sim.phi_max = 0.8; 47 sim.dilatancy_angle = 1.0; 48 49 /* Iverson et al 1997, 1998: Storglaciaren till */ 50 /* sim.mu_s = tan(DEG2RAD(26.3)); */ 51 /* sim.C = 5.0e3; */ 52 /* sim.phi = initval(0.22, 1); */ 53 /* sim.d = ??; */ 54 55 /* Iverson et al 1997, 1998: Two Rivers till */ 56 /* sim.mu_s = tan(DEG2RAD(17.8)); */ 57 /* sim.C = 14.0e3; */ 58 /* sim.phi = initval(0.37, 1); */ 59 /* sim.d = ??; */ 60 61 /* Tulaczyk et al 2000a: Upstream B till */ 62 /* sim.mu_s = tan(DEG2RAD(23.9)); */ 63 /* sim.C = 3.0e3; */ 64 /* sim.phi = initval(0.35, 1); */ 65 /* sim.d = ??; */ 66 67 /* grain material density [kg/m^3] */ 68 sim.rho_s = 2.6e3; /* Damsgaard et al 2013 */ 69 70 /* spatial settings */ 71 sim.origo_z = 0.0; 72 sim.L_z = 1.0; 73 74 /* temporal settings */ 75 sim.t = 0.0; 76 sim.dt = 2.0; 77 sim.t_end = 1.0; 78 sim.file_dt = 0.1; 79 sim.n_file = 0; 80 81 /* fluid settings */ 82 sim.fluid = 0; 83 84 sim.rho_f = 1e3; 85 86 /* Water at 20 deg C, Goren et al 2011 */ 87 /* sim.beta_f = 4.5e-10; */ 88 /* sim.mu_f = 1.0-3; */ 89 90 /* Water at 0 deg C */ 91 sim.beta_f = 3.9e-10; /* doi:10.1063/1.1679903 */ 92 sim.mu_f = 1.787e-3; /* Cuffey and Paterson 2010 */ 93 94 /* Damsgaard et al 2015 */ 95 sim.k = initval(1.9e-15, 1); 96 97 /* Iverson et al 1994: Storglaciaren */ 98 /* sim.k = initval(1.3e-14, 1); */ 99 100 /* Engelhardt et al 1990: Upstream B */ 101 /* sim.k = initval(2.0e-16, 1); */ 102 103 /* Leeman et al 2016: Upstream B till */ 104 /* sim.k = initval(4.9e-17, 1); */ 105 106 /* no fluid-pressure variations */ 107 sim.p_f_top = 0.0; 108 sim.p_f_mod_ampl = 0.0; 109 sim.p_f_mod_freq = 1.0; 110 sim.p_f_mod_phase = 0.0; 111 sim.p_f_mod_pulse_time = NAN; 112 sim.p_f_mod_pulse_shape = 0; 113 } 114 115 #endif