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