|
22 | 22 | G_l = 1e06 |
23 | 23 |
|
24 | 24 |
|
25 | | -# primitive vartiables |
| 25 | +# primitive variables |
26 | 26 | patmos = 101325.0 # pa |
27 | 27 |
|
28 | 28 | # problem specific variable |
|
35 | 35 |
|
36 | 36 | # non-dim |
37 | 37 |
|
38 | | -# define characteristic density, length, time, stress material #make it liquid |
| 38 | +# define characteristic density, length, time, stress material |
| 39 | +#make it liquid |
39 | 40 | rho_char = rhol |
40 | 41 | length_char = lambda_wave |
41 | 42 | c_char = c_l # should be liquid |
|
86 | 87 | print( |
87 | 88 | json.dumps( |
88 | 89 | { |
89 | | - # Logistics ================================================================ |
90 | 90 | "run_time_info": "T", |
91 | | - #'sim_data' : 'T', |
92 | | - # ========================================================================== |
93 | | - # Computational Domain Parameters ========================================== |
| 91 | + # Computational Domain Parameters |
94 | 92 | "x_domain%beg": -dlengx / 2.0, |
95 | 93 | "x_domain%end": dlengx / 2.0, |
96 | 94 | "y_domain%beg": 0.0, |
|
109 | 107 | "t_step_start": tstart, |
110 | 108 | "t_step_stop": tstop, |
111 | 109 | "t_step_save": tsave, |
112 | | - # ========================================================================== |
113 | | - # Simulation Algorithm Parameters ========================================== |
| 110 | + # Simulation Algorithm Parameters |
114 | 111 | "num_patches": 2, |
115 | 112 | "model_eqns": 3, |
116 | 113 | "alt_soundspeed": "F", |
|
135 | 132 | "bc_y%end": -1, |
136 | 133 | "bc_z%beg": -1, |
137 | 134 | "bc_z%end": -1, |
138 | | - # ========================================================================== |
139 | | - # Turning on Hypoelasticity ================================================ |
140 | 135 | #'hypoelasticity' : 'T', |
141 | 136 | "hyperelasticity": "F", |
142 | | - # ========================================================================== |
143 | | - # Formatted Database Files Structure Parameters ============================ |
| 137 | + # Formatted Database Files Structure Parameters |
144 | 138 | "format": 1, |
145 | 139 | "precision": 2, |
146 | 140 | "prim_vars_wrt": "T", |
147 | 141 | "parallel_io": "T", |
148 | | - # ========================================================================== |
149 | | - # acoustics setting ========================================================= |
| 142 | + # acoustics setting |
150 | 143 | "acoustic_source": "T", |
151 | 144 | "num_source": 1, |
152 | 145 | "acoustic(1)%support": 3, |
|
160 | 153 | "acoustic(1)%length": dlengy, # length of the place ??? |
161 | 154 | "acoustic(1)%height": dlengz, |
162 | 155 | "acoustic(1)%dir": -math.pi, |
163 | | - # =========================================================================== |
164 | | - # Patch 1: Background ====================================================== |
| 156 | + # Patch 1: Background |
165 | 157 | "patch_icpp(1)%geometry": 9, |
166 | 158 | "patch_icpp(1)%x_centroid": 0.0, |
167 | 159 | "patch_icpp(1)%y_centroid": dlengy / 2.0, |
|
178 | 170 | "patch_icpp(1)%alpha(1)": alphal_back, |
179 | 171 | "patch_icpp(1)%alpha(2)": alphag_back, |
180 | 172 | #'patch_icpp(1)%tau_e(1)' : 0.0, |
181 | | - # ========================================================================== |
182 | | - # Patch 2: Lung ============================================================ |
| 173 | + # Patch 2: Lung |
183 | 174 | "patch_icpp(2)%geometry": 13, |
184 | 175 | "patch_icpp(2)%hcid": 301, |
185 | 176 | #'patch_icpp(2)%geometry' : 9, |
|
199 | 190 | "patch_icpp(2)%alpha_rho(2)": rhog_n * alphag_lung, |
200 | 191 | "patch_icpp(2)%alpha(1)": alphal_lung, |
201 | 192 | "patch_icpp(2)%alpha(2)": alphag_lung, |
202 | | - # ========================================================================== |
203 | | - # Fluids Physical Parameters =============================================== |
| 193 | + # Fluids Physical Parameters |
204 | 194 | "fluid_pp(1)%gamma": 1.0e00 / (gammal - 1.0e00), |
205 | 195 | "fluid_pp(1)%pi_inf": gammal * Bl_n / (gammal - 1.0e00), |
206 | 196 | "fluid_pp(1)%G": G_l_n, |
207 | 197 | "fluid_pp(2)%gamma": 1.0e00 / (gammag - 1.0e00), |
208 | 198 | "fluid_pp(2)%pi_inf": gammag * Bg_n / (gammag - 1.0e00), |
209 | 199 | "fluid_pp(2)%G": G_g_n, |
210 | | - # ============================================================================== |
211 | 200 | } |
212 | 201 | ) |
213 | 202 | ) |
214 | 203 |
|
215 | | -# ============================================================================== |
0 commit comments