Posts for the month of August 2009

Dissolution Simulations

The simulations in the table below are being run as I write up the paper, most are very close to production quality.

ID Notes diff Issues Status $N$ movie
119 hydrodynamic limit solid waves, wrong value of X1 close to equilibrium stopped 32 movie
120 removed pytables dependencies to run on rareton diff stopped
121 finer grid density diff better initial conditions required stopped 32 movie
122 $\rho_i = 0.1 \left( \rho_i^s - \rho_i^l \right) \rho_i^l$ and deeper solid diff will induce solidification stopped 6 movie
123 $\rho_i = 0.5 \left( \rho_i^s - \rho_i^l \right) \rho_i^l$ diff drop wets top surface stopped 11
124 $\rho_i = 0.5 \left( \rho_i^s - \rho_i^l \right) \rho_i^l$ and deeper vapor diff will induce solidification stopped 11 movie
125 $\mu_s^{\text{new}} = 100 \mu_s $, hydrodynamic limit diff less waves, but better initial conditions required stopped 12
126 $X_1^{l, \text{new}} = 0.5 X_1^l$ diff no melting stopped 6
127 $M_{\phi}^{\text{new}} = 10 M_{\phi}$ diff increasing mobility will not induce melting stopped 12
128 $R_{\text{tol}}^{\text{new}} = 0.1 R_{\text{tol}}$ diff better initial conditions required stopped 12
129 $\mu_s^{\text{new}} = \mu_l$ and $\bar{M}_s^{\text{new}} = \bar{M}_l$ diff stopped 12
130 $M_{\phi}^{\text{new}} = 10 M_{\phi}$ diff increasing mobility will not induce melting stopped 4
131 $M_{\phi}^{\text{new}} = 10 M_{\phi}$ diff increasing mobility will not induce melting stopped 4
132 $\bar{M}_f = 1 \times 10^{3} \bar{M}_f$ diff X1 has weird profile due to smoothing stopped 12
133 smoothing $\rho$ and $X_1$ diff worked, but fast diffusion not required sropped 4
134 1D to see eventual melting diff not required smaller nx works stopped 4
135 $\bar{M}_f = 1 \times 10^{5} \bar{M}_f$ diff duffusion fast enough, large dissolution impossible stopped 4
136 $\bar{M}_f = 1 \times 10^{5} \bar{M}_f$, 1D, nx=100 diff not required smaller nx works stopped 4
137 1D, nx=50, $X_1^{l, \text{new}} = 0.1 X_1^l$ diff wrong initial conditions stopped 4
138 hydrodynamic limit diff stopped 12
139 1D, nx=10, $X_1^{l, \text{new}} = 0.1 X_1^l$, shift=0.0 diff nx < 29 doesn't seem to work or different shifts stopped 4
140 $X_1^{l, \text{new}} = 0.1 X_1^l$ diff stopped 4
141 $\bar{M}_f = 1 \times 10^{-6} $ diff stopped 12
142 nx=2 diff calculation broke down with nx =2 stopped 4
143 same as 137 but corrected initial conditions diff stopped 4
144 same as 141 but $\bar{M}_f = 1 \times 10^{-5} $ diff stopped 4
145 143: $\bar{M}_s = 1 \times 10^{-2} $ diff stopped 4
146 63: with modifed unsegregatedSolver from 145 $X_1^{l, \text{new}} = 0.1 X_1^l$ and NC1, NC2 = 0 between sweeps diff stopped 1 movie
147 146: $\bar{M}_s = \bar{M}_f = 1 \times 10^{-2}$ diff stopped 1 movie
148 147: $\bar{M}_s = 1 \times 10^{-11}$ diff stopped 1 movie
149 147: $\mu_s = 2 \times 10^{6}$ diff stopped 1 movie
150 148: $\mu_s = 2 \times 10^{6}$ diff stopped 1 movie
151 150: $\bar{M}_s = 1 \times 10^{-1}$ diff stopped 1
152 150: $\bar{M}_s = 1 \times 10^{-3}$ diff stopped 1
153 150: $\bar{M}_s = 1 \times 10^{-4}$ diff stopped 1
154 150: $\bar{M}_s = 1 \times 10^{-5}$ diff stopped 1
155 150: $\bar{M}_s = 1 \times 10^{-6}$ diff stopped 1
156 150: $\bar{M}_s = 1 \times 10^{-7}$ diff stopped 1
157 150: $\bar{M} = \bar{M}_l$ when $\phi < 0.97$ diff stopped 1
158 157: $\bar{M} = \bar{M}_s^{\phi^2} \bar{M}_s^{1 - \phi^2} $ diff stopped 1
159 154: half interface depth and N=1600 diff stopped 1
160 154: $N=1600 $ diff stopped 1
161 159: $\epsilon_i^{\text{new}} = 4 \epsilon_i $ diff stopped 1
162 149: $N=1600 $ diff stopped 1
163 154: $\epsilon_i^{\text{new}} = \epsilon_i / 4 $ diff stopped 1
164 154: $\epsilon_i^{\text{new}} = \epsilon_i / 16 $ diff stopped 1
165 154: $\epsilon_i^{\text{new}} = \epsilon_i / 64 $ diff stopped 1
166 150: $\bar{M}_{\text{int}} = 1 \times 10^{-3}$ diff stopped 1
167 166: $l=1$ diff stopped 1
168 167: box interpolation diff stopped 1
169 167: $\bar{M}_f = 1 \times 10^{-6}$, $l=2$ diff stopped 1
170 167: $\bar{M}_f = 1 \times 10^{-5}$, $l=2$ diff stopped 1
171 170: $\bar{M}_f = 1 \times 10^{-6}$, $l=1$ diff stopped 1
172 170: $\bar{M}_f = 1 \times 10^{-5}$, $l=1$ diff stopped 1
173 167: $\bar{M}_f = 1 \times 10^{-7}$, $l=0.5$ diff stopped 1
174 173: $\bar{M}_f = 1 \times 10^{-6}$, $l=0.5$ diff stopped 1
175 174: $\bar{M}_f = 1 \times 10^{-5}$, $l=0.5$ diff stopped 1
176 175: $\bar{M}_f = 1 \times 10^{-7}$, sharp 0.95 and 0.05 diff stopped 1
177 176: $\bar{M}_f = 1 \times 10^{-6}$ diff stopped 1
178 177: $\bar{M}_f = 1 \times 10^{-5}$ diff stopped 1
179 176: $\bar{M}_i = 1 \times 10^{-2}$ diff stopped 1
180 179: $\bar{M}_i = 1 \times 10^{-4}$ diff stopped 1
181 179: $\bar{M}_i = 1 \times 10^{-5}$ diff stopped 1
182 179: $\bar{M}_i = 1 \times 10^{-6}$ diff stopped 1
184 138: 2D, $\bar{M}_i = 1 \times 10^{-7}$ diff stopped 12
185 184: 2D, $\bar{M}_i = 1 \times 10^{-5}$ diff stopped 12
186 184: 2D, $\bar{M}_i = 1 \times 10^{-3}$ diff stopped 4
187 184: 2D, $\bar{M}_i = 1 \times 10^{-1}$ diff stopped 12
188 186: 2D, $X_1 = 0.5 X_1^{equ}$ diff stopped 12
189 186: 2D, $X_1 = 0.1 X_1^{equ}$ diff calculations broke down stopped 4
190 186: 2D, $X_1 = 0.01 X_1^{equ}$ diff calculations broke down stopped 4
191 189: 2D, $\bar{M}_f = 1 \times 10^{-6}$ diff stopped 4
192 189: 2D, $\bar{M}_f = 1 \times 10^{-5}$ diff stopped 4
193 189: 2D, $\bar{M}_f = 1 \times 10^{-4}$ diff stopped 4
194 189: 2D, $\bar{M}_f = 1 \times 10^{-3}$ diff stopped 4
195 186: 2D, $X_1 = 0.3 X_1^{equ}$, $R_{tol} = 1 \times 10^{-2}$ diff stopped 12
196 195: 2D, $X_1 = 1.0 X_1^{equ}$ diff stopped 12
197 196: 2D, factor=2 diff stopped 32
198 185: 2D, $\mu = \mu_s^{\phi^4} \mu_l^{\left(1 - \phi^4\right)}$ diff stopped 12
199 185: 2D, $\mu = \mu_s^{\phi^6} \mu_l^{\left(1 - \phi^6\right)}$ diff stopped 12
200 185: 2D, $\mu = \mu_s^{\phi^{10}} \mu_l^{\left(1 - \phi^{10}\right)}$ diff stopped 12
201 185: 2D, $\mu = \mu_s^{\phi^8} \mu_l^{\left(1 - \phi^8\right)}$ diff stopped 12
202 199: 2D, $\mu = \mu_s^{\phi^2} \mu_l^{\left(1 - \phi^2\right)}$ diff stopped 4
203 199: 2D, $\mu = \mu_s^{\phi^3} \mu_l^{\left(1 - \phi^3\right)}$ diff stopped 4
204 198: 2D, $X_1 = 0.5 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-7}$ diff stopped 12
205 198: 2D, $X_1 = 0.1 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-7}$ $\Delta t < 1 \times 10^{-11} $ diff crashed benson, hobson, alice (stopped) 12
206 198: 2D, $X_1 = 1.0 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-6}$ diff crashed stopped 12
207 198: 2D, $X_1 = 0.5 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-6}$ diff crashed stopped 4
208 198: 2D, $X_1 = 0.1 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-6}$ diff stopped 4
209 198: 2D, $X_1 = 0.1 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-5}$ diff crashed cluster 411 (stopped) 4
210 198: 2D, $X_1 = 0.1 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-4}$ diff cluster 412 (stopped) 4
211 198: 2D, $X_1 = 0.1 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-3}$ diff cluster 413 (stopped) 4
212 198: 2D, $X_1 = 1.0 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-7}$ $\nu_f = 2 \times 10^{-2}$ $\Delta t < 1 \times 10^{-11} $ diff crashed renfield, sancho, kato (stopped) 12
213 212: 2D, $X_1 = 1.0 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-7}$ $\nu_f = 2 \times 10^{-1}$ 5000 iterations diff crashed cluster 418 (stopped) 4
214 205: 2D, $X_1 = 0.1 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-7}$ $\nu_f = 2 \times 10^{-3}$ $R_{\text{tol}} = 1 \times 10^{-2} $ diff crashed poole (stopped) 4
215 212: 2D, $X_1 = 1.0 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-7}$ $\nu_f = 2 \times 10^{-2}$ $R_{\text{tol}} = 1 \times 10^{-2} $ diff crashed luggage (stopped) 12
216 205: 2D, $X_1 = 0.1 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-7}$ $\nu_f = 2 \times 10^{-3}$ $R_{\text{tol}} = 1 \times 10^{-1} $ $\Delta_{\text{CFL}} = 5 \times 10^{-2} $ diff crashed benson, hobson, alice (stopped) 12
217 212: 2D, $X_1 = 1.0 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-7}$ $\nu_f = 2 \times 10^{-2}$ $R_{\text{tol}} = 1 \times 10^{-1} $ $\Delta_{\text{CFL}} = 5 \times 10^{-2} $ diff crashed renfield, sancho, kato (stopped) 12
218 198: $ \bar{M} = \bar{M}_f^{1 - \phi^m} \bar{M}_s^{\phi^m} $ diff crashed luggage (stopped) 12
219 218: $ m = 2 $ diff cluster 424 (stopped) 4
220 218: $ m = 3 $ diff cluster 425 (stopped) 4
221 218: $ m = 4 $ diff cluster 426 (stopped) 4

New base simulation

ID Notes diff Issues Status $N$ movie
222 base diff benson, hobson, alice 12
223 $X_1 = 0.5 X_1^{equ}$ diff None luggage (stopped) 12
224 $X_1 = 0.1 X_1^{equ}$ diff crashed luggage (stopped) 12
225 $\nu_f = 2 \times 10^{-2}$ diff luggage 12
226 $\nu_f = 2 \times 10^{-1}$ diff crashed renfield, sancho, kato (stopped) 12
227 224: $X_1 = 0.1 X_1^{equ}$, $R_{\text{solver}}=1 \times 10^{-9}$, $\sigma_{\text{CFL}} = 0.05$ diff crashed luggage (stopped) 12
228 225: $\nu_f = 2 \times 10^{-2}$, $R_{\text{solver}}=1 \times 10^{-9}$, $\sigma_{\text{CFL}} = 0.05$ diff crashed renfield, sancho, kato (stopped) 12
229 226: $\nu_f = 2 \times 10^{-1}$, $R_{\text{tol}}=1 \times 10^{-3}$, $\alpha_{\text{relax}} = 1.0$ diff crashed cluster (stopped) 4
230 222: $X_1 = 2.0 X_1^{equ}$ diff crashed cluster 12.0 (stopped) 4
231 222: $\nu_f = 1 \times 10^{-2}$ diff crashed oddjob (stopped) 4
232 222: $\nu_f = 2 \times 10^{-2}$ $\bar{M}_f = 1 \times 10^{-6}$ diff crashed oddjob (stopped) 4
233 222: $X_1 = 0.1 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-6}$ diff cluster 7.0 4
234 225: $\sigma_{\text{CFL}} = 0.02$ diff crashed oddjob (stopped) 3
235 225: $\nu_f = 4 \times 10^{-3}$ diff None cluster 2.0 (stopped) 4
236 225: $\alpha = 0.2$ diff crashed cluster 31.0 (stopped) 4
237 224: $\sigma_{\text{CFL}} = 0.02$ diff crashed renfield, sancho, kato (stopped) 12
238 224: $\alpha = 0.2$ diff crashed cluster 3.0 (stopped) 4
239 230: $X_1 = 1.5 X_1^{equ}$ diff crashed cluster 4.0 (stopped) 4
240 222: $X_1 = 1.0 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-6}$ diff crashed luggage (stopped) 12
241 223: $X_1 = 0.5 X_1^{equ}$ $\bar{M}_f = 1 \times 10^{-6}$ diff crashed cluster 5.0 (stopped) 4
242 235: $X_1 = 1.0 X_1^{equ}$ $\nu_f = 8 \times 10^{-3}$ diff crashed cluster 6.0 (stopped) 4
243 226: $\nu_f = 2 \times 10^{-1}$, shift=dx diff None renfield, sancho, kato (stopped) 12
244 226: $\nu_f = 2 \times 10^{-1}$, shift=10 * dx diff None luggage (stopped) 12
245 226: $\nu_f = 2 \times 10^{-1}$, shift=100 * dx diff None cluster 28 (stopped) 4
246 226: $\nu_f = 2 \times 10^{-1}$ diff None cluster 29 (stopped) 1
247 225: $\nu_f = 2 \times 10^{-2}$, shift=dx diff crashed luggage (stopped) 12
248 224: $X_1 = 0.1 X_1^{equ}$, shift=dx diff crashed cluster 30 (stopped) 4
249 247: $\nu_f = 2 \times 10^{-2}$, shift=2 * dx diff crashed renfield, sancho, kato (stopped) 12
250 247: $\nu_f = 2 \times 10^{-2}$, $\sigma_{\text{CFL}} = 0.05$, shift=dx diff crashed luggage (stopped) 12
251 240: $X_1 = 1.0 X_1^{equ}$, $\bar{M}_f = 1 \times 10^{-6}$, shift=dx diff crashed cluster 31 (stopped) 4
252 243: $\nu_f = 2 \times 10^{-1}$, shift=dx, $\sigma_{\text{CFL}} = 0.01$ diff cluster 32 4
253 243: $\nu_f = 2 \times 10^{-4}$, shift=dx diff crashed cluster 33 (stopped) 4
254 253: $\nu_f = 2 \times 10^{-4}$, shift=dx, $\nu_s = 2 \times 10^{5}$ diff crashed renfield, sancho, kato (stopped) 12
255 254: $\nu_f = 2 \times 10^{-4}$, shift=dx, $\nu_s = 2 \times 10^{4}$ diff crashed renfield, sancho, kato (stopped) 12
256 255: $\nu_f = 2 \times 10^{-4}$, shift=dx, $\nu_s = 2 \times 10^{7}$ diff crashed cluster 34 (stopped) 4
257 251: $X_1 = 1.0 X_1^{equ}$, $\bar{M}_f = 1 \times 10^{-6}$, shift=2*dx, $\sigma_{\text{CFL}}=0.05$ diff renfield, sancho, kato 12
258 256: $\nu_f = 2 \times 10^{-4}$, shift=dx, $\nu_s = 2 \times 10^{3}$ diff crashed cluster 35 (stopped) 4
259 248: $X_1 = 0.1 X_1^{equ}$, shift=2 * dx diff crashed cluster 36 (stopped) 4
260 250: $\nu_f = 2 \times 10^{-2}$, $\sigma_{\text{CFL}} = 0.05$, shift=2 * dx diff crashed renfield, sancho, kato (stopped) 12
261 259: $X_1 = 0.1 X_1^{equ}$, shift=10 * dx diff crashed luggage (stopped) 12
262 257: $X_1 = 1.0 X_1^{equ}$, $\bar{M}_f = 1 \times 10^{-6}$, shift=10*dx, $\sigma_{\text{CFL}}=0.05$ diff crashed bunter (stopped) 4
263 260: $\nu_f = 2 \times 10^{-2}$, $\sigma_{\text{CFL}} = 0.05$, shift=10 * dx diff crashed renfield, sancho, kato (stopped) 12
  • Posted: 2009-08-27 13:02 (Updated: 2009-12-15 15:24)
  • Author: wd15
  • Categories: (none)
  • Comments (0)

Fixing LAMMPI bug

Had a weird problem with the following error when using LAMMPI. 

-----------------------------------------------------------------------------
The selected RPI failed to initialize during MPI_INIT.  This is a
fatal error; I must abort.
This occurred on host poole (n0).
The PID of failed process was 29080 (MPI_COMM_WORLD rank: 0)
-----------------------------------------------------------------------------
-----------------------------------------------------------------------------
One of the processes started by mpirun has exited with a nonzero exit
code.  This typically indicates that the process finished in error.
If your process did not finish in error, be sure to include a "return
0" or "exit(0)" in your C code before exiting the application.
PID 29081 failed on node n0 (127.0.0.1) with exit status 1.
-----------------------------------------------------------------------------

The solution involved using the command "ipcs" and "ipcrm" and is explained here:

<http://www.riothamus.net/drupal/?q=node/12>

The above post has a perl script to run it use 

exec `./ipcrmAll.pl -s`
  • Posted: 2009-08-13 11:43 (Updated: 2009-08-13 11:55)
  • Author: wd15
  • Categories: (none)
  • Comments (0)