|
144 | 144 | k_6 = \mathbf{u}\!\bigl(\mathbf{x}_n + \Delta t(-\tfrac{3}{2}k_1 + 2k_2 - \tfrac{1}{2}k_3 + k_4),\;t_n + \Delta t\bigr),\\ |
145 | 145 | \mathbf{x}_{n+1} = \mathbf{x}_n + \Delta t\bigl(\tfrac{1}{20}k_1 + \tfrac{1}{4}k_4 + \tfrac{1}{5}k_5 + \tfrac{1}{2}k_6\bigr).\end{split}\]</div> |
146 | 146 | <p>This scheme incurs a global error of order <span class="math notranslate nohighlight">\(O(\Delta t^6)\)</span> with six velocity evaluations.</p> |
| 147 | +<table class="docutils align-default"> |
| 148 | +<thead> |
| 149 | +<tr class="row-odd"><th class="head"><p><code class="docutils literal notranslate"><span class="pre">output_dir</span></code></p></th> |
| 150 | +<th class="head"><p>Path to directory for saving results</p></th> |
| 151 | +<th class="head"><p><code class="docutils literal notranslate"><span class="pre"><string></span></code></p></th> |
| 152 | +<th class="head"></th> |
| 153 | +</tr> |
| 154 | +<tr class="row-even"><th class="head"><p><code class="docutils literal notranslate"><span class="pre">auto_wall</span></code></p></th> |
| 155 | +<th class="head"><p>Flag for automatic wall detection</p></th> |
| 156 | +<th class="head"><p><code class="docutils literal notranslate"><span class="pre">0</span></code></p></th> |
| 157 | +<th class="head"><p>Manually indicate computation bounds.</p></th> |
| 158 | +</tr> |
| 159 | +<tr class="row-odd"><th class="head"></th> |
| 160 | +<th class="head"></th> |
| 161 | +<th class="head"><p>(<code class="docutils literal notranslate"><span class="pre">1</span></code>)</p></th> |
| 162 | +<th class="head"><p>Auto detect bounds from input mesh.</p></th> |
| 163 | +</tr> |
| 164 | +<tr class="row-even"><th class="head"><p><code class="docutils literal notranslate"><span class="pre">wall_xrange</span></code></p></th> |
| 165 | +<th class="head"><p>X-direction wall bounds [xmin, xmax]</p></th> |
| 166 | +<th class="head"><p><code class="docutils literal notranslate"><span class="pre">[<float>,</span> <span class="pre"><float>]</span></code></p></th> |
| 167 | +<th class="head"><p>Dimensionalized coordinates bounds.</p></th> |
| 168 | +</tr> |
| 169 | +<tr class="row-odd"><th class="head"><p><code class="docutils literal notranslate"><span class="pre">wall_yrange</span></code></p></th> |
| 170 | +<th class="head"><p>Y-direction wall bounds [ymin, ymax]</p></th> |
| 171 | +<th class="head"><p><code class="docutils literal notranslate"><span class="pre">[<float>,</span> <span class="pre"><float>]</span></code></p></th> |
| 172 | +<th class="head"></th> |
| 173 | +</tr> |
| 174 | +<tr class="row-even"><th class="head"><p><code class="docutils literal notranslate"><span class="pre">wall_zrange</span></code></p></th> |
| 175 | +<th class="head"><p>Z-direction wall bounds [zmin, zmax]</p></th> |
| 176 | +<th class="head"><p><code class="docutils literal notranslate"><span class="pre">[<float>,</span> <span class="pre"><float>]</span></code></p></th> |
| 177 | +<th class="head"></th> |
| 178 | +</tr> |
| 179 | +<tr class="row-odd"><th class="head"><p><code class="docutils literal notranslate"><span class="pre">x</span></code></p></th> |
| 180 | +<th class="head"><p>Number of Lagrangian seed points in X direction</p></th> |
| 181 | +<th class="head"><p><code class="docutils literal notranslate"><span class="pre"><int></span></code></p></th> |
| 182 | +<th class="head"></th> |
| 183 | +</tr> |
| 184 | +<tr class="row-even"><th class="head"><p><code class="docutils literal notranslate"><span class="pre">y</span></code></p></th> |
| 185 | +<th class="head"><p>Number of Lagrangian seed points in Y direction</p></th> |
| 186 | +<th class="head"><p><code class="docutils literal notranslate"><span class="pre"><int></span></code></p></th> |
| 187 | +<th class="head"></th> |
| 188 | +</tr> |
| 189 | +<tr class="row-odd"><th class="head"><p><code class="docutils literal notranslate"><span class="pre">z</span></code></p></th> |
| 190 | +<th class="head"><p>Number of Lagrangian seed points in Z direction</p></th> |
| 191 | +<th class="head"><p><code class="docutils literal notranslate"><span class="pre"><int></span></code></p></th> |
| 192 | +<th class="head"></th> |
| 193 | +</tr> |
| 194 | +<tr class="row-even"><th class="head"><p><code class="docutils literal notranslate"><span class="pre">T_range</span></code></p></th> |
| 195 | +<th class="head"><p>Time range [<code class="docutils literal notranslate"><span class="pre">T_start</span></code>, <code class="docutils literal notranslate"><span class="pre">T_end</span></code>] in frames</p></th> |
| 196 | +<th class="head"><p><code class="docutils literal notranslate"><span class="pre">[<float>,</span> <span class="pre"><float>]</span></code></p></th> |
| 197 | +<th class="head"><p><code class="docutils literal notranslate"><span class="pre">T_end</span></code> no larger than <code class="docutils literal notranslate"><span class="pre">N_frames</span></code>.</p></th> |
| 198 | +</tr> |
| 199 | +<tr class="row-odd"><th class="head"><p><code class="docutils literal notranslate"><span class="pre">dt</span></code></p></th> |
| 200 | +<th class="head"><p>Time step for particle advection in frames</p></th> |
| 201 | +<th class="head"><p><code class="docutils literal notranslate"><span class="pre"><float></span></code></p></th> |
| 202 | +<th class="head"></th> |
| 203 | +</tr> |
| 204 | +</thead> |
| 205 | +<tbody> |
| 206 | +<tr class="row-even"><td><p><code class="docutils literal notranslate"><span class="pre">physical_dt</span></code></p></td> |
| 207 | +<td><p>Physical time step between two input frames</p></td> |
| 208 | +<td><p><code class="docutils literal notranslate"><span class="pre"><float></span></code></p></td> |
| 209 | +<td><p>Perform physical advection.</p></td> |
| 210 | +</tr> |
| 211 | +</tbody> |
| 212 | +</table> |
147 | 213 | </section> |
148 | 214 | </section> |
149 | 215 | <section id="ftle-computation"> |
|
0 commit comments