You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
<tdclass="markdownTableBodyRight"><code>vel(i)</code> * </td><tdclass="markdownTableBodyCenter">Real </td><tdclass="markdownTableBodyCenter">Supported </td><tdclass="markdownTableBodyLeft">Velocity in direction $i$. </td></tr>
<tdclass="markdownTableBodyRight"><code>cf_val</code> * </td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyCenter">Supported</td><tdclass="markdownTableBodyLeft">Surface tension color function value</td></tr>
<tdclass="markdownTableBodyRight"><code>model_filepath</code></td><tdclass="markdownTableBodyCenter">String</td><tdclass="markdownTableBodyCenter">Not Supported </td><tdclass="markdownTableBodyLeft">Path to an STL or OBJ file (not all OBJs are supported).</td></tr>
287
+
<tdclass="markdownTableBodyRight"><code>cf_val</code>* </td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyCenter">Supported </td><tdclass="markdownTableBodyLeft">Surface tension color function value</td></tr>
288
288
<trclass="markdownTableRowEven">
289
-
<tdclass="markdownTableBodyRight"><code>model_scale(i)</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyCenter">Not Supported </td><tdclass="markdownTableBodyLeft">Model's (applied) scaling factor for component $i$. </td></tr>
289
+
<tdclass="markdownTableBodyRight"><code>model_filepath</code></td><tdclass="markdownTableBodyCenter">String</td><tdclass="markdownTableBodyCenter">Not Supported </td><tdclass="markdownTableBodyLeft">Path to an STL or OBJ file (not all OBJs are supported). </td></tr>
290
290
<trclass="markdownTableRowOdd">
291
-
<tdclass="markdownTableBodyRight"><code>model_rotate(i)</code></td><tdclass="markdownTableBodyCenter">Real </td><tdclass="markdownTableBodyCenter">Not Supported </td><tdclass="markdownTableBodyLeft">Model's (applied) angle of rotation about axis $i$. </td></tr>
<tdclass="markdownTableBodyRight"><code>model_translate(i)</code></td><tdclass="markdownTableBodyCenter">Real </td><tdclass="markdownTableBodyCenter">Not Supported </td><tdclass="markdownTableBodyLeft">Model's $i$-th component of (applied) translation. </td></tr>
293
+
<tdclass="markdownTableBodyRight"><code>model_rotate(i)</code></td><tdclass="markdownTableBodyCenter">Real </td><tdclass="markdownTableBodyCenter">Not Supported </td><tdclass="markdownTableBodyLeft">Model's (applied) angle of rotation about axis $i$. </td></tr>
294
294
<trclass="markdownTableRowOdd">
295
-
<tdclass="markdownTableBodyRight"><code>model_spc</code></td><tdclass="markdownTableBodyCenter">Integer</td><tdclass="markdownTableBodyCenter">Not Supported </td><tdclass="markdownTableBodyLeft">Number of samples per cell when discretizing the model into the grid. </td></tr>
295
+
<tdclass="markdownTableBodyRight"><code>model_translate(i)</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyCenter">Not Supported </td><tdclass="markdownTableBodyLeft">Model's $i$-th component of (applied) translation. </td></tr>
296
296
<trclass="markdownTableRowEven">
297
+
<tdclass="markdownTableBodyRight"><code>model_spc</code></td><tdclass="markdownTableBodyCenter">Integer </td><tdclass="markdownTableBodyCenter">Not Supported </td><tdclass="markdownTableBodyLeft">Number of samples per cell when discretizing the model into the grid. </td></tr>
298
+
<trclass="markdownTableRowOdd">
297
299
<tdclass="markdownTableBodyRight"><code>model_threshold</code></td><tdclass="markdownTableBodyCenter">Real </td><tdclass="markdownTableBodyCenter">Not Supported </td><tdclass="markdownTableBodyLeft">Ray fraction inside the model patch above which the fraction is set to one. </td></tr>
298
300
</table>
299
301
<p>*: These parameters should be prepended with <code>patch_icpp(j)%</code> where $j$ is the patch index.</p>
300
302
<p>The Table lists the patch parameters. The parameters define the geometries and physical parameters of fluid components (patch) in the domain at initial condition. Note that the domain must be fully filled with patche(s). The code outputs error messages when an empty region is left in the domain.</p>
303
+
<ul>
304
+
<li><code>tau_e(i)</code> is the <code>i</code>-th component of the elastic stress tensor, ordered as <code>tau_xx</code>, <code>tau_xy</code>, <code>tau_yy</code>, <code>tau_xz</code>, <code>tau_yz</code>, and <code>tau_zz</code>. 1D simulation requires <code>tau(1)</code>, 2D <code>tau(1:3)</code>, and 3D <code>tau(1:6)</code>.</li>
305
+
</ul>
301
306
<h3><aclass="anchor" id="autotoc_md8"></a>
302
307
Analytical Definition of Primitive Variables</h3>
303
308
<p>Some parameters, as described above, can be defined by analytical functions in the input file. For example, one can define the following patch:</p>
<p>When these parameters are undefined, fluids are treated as inviscid. Details of implementation of viscosity in MFC can be found in <ahref="references.md#Coralic15">Coralic (2015)</a>.</p>
442
449
<ul>
443
450
<li><code>fluid_pp(i)%cv</code>, <code>fluid_pp(i)%qv</code>, and <code>fluid_pp(i)%qvp</code> define $c_v$, $q$, and $q'$ as parameters of $i$-th fluid that are used in stiffened gas equation of state.</li>
451
+
<li><code>fluid_pp(i)%G</code> is required for <code>hypoelasticity</code>.</li>
<li><code>weno_avg</code> it activates the arithmetic average of the left and right, WENO-reconstructed, cell-boundary values. This option requires <code>weno_Re_flux</code> to be true because cell boundary values are only utilized when employing the scalar divergence method in the computation of velocity gradients.</li>
559
569
<li><code>surface_tension</code> activates surface tension when set to <code>'T'</code>. Requires <code>sigma</code> to be set and <code>num_fluids</code>. The color function in each patch should be assigned such that <code>patch_icpp(i)cf_val = 1</code> in patches where <code>patch_icpp(i)alpha = 1 - eps</code> and <code>patch_icpp(i)cf_val = 0</code> in patches where <code>patch_icpp(i)alpha = eps</code>.</li>
560
570
<li><code>viscous</code> activates viscosity when set to <code>'T'</code>. Requires <code>Re(1)</code> and <code>Re(2)</code> to be set.</li>
571
+
<li><code>hypoelasticity</code> activates elastic stress calculations for fluid-solid interactions. Requires <code>G</code> to be set in the fluid material's parameters.</li>
<tdclass="markdownTableBodyRight"><code>qm_wrt</code></td><tdclass="markdownTableBodyCenter">Logical </td><tdclass="markdownTableBodyLeft">Add the Q-criterion to the database </td></tr>
627
638
<trclass="markdownTableRowEven">
628
-
<tdclass="markdownTableBodyRight"><code>fd_order</code></td><tdclass="markdownTableBodyCenter">Integer</td><tdclass="markdownTableBodyLeft">Order of finite differences for computing the vorticity and the numerical Schlieren function [1,2,4]</td></tr>
639
+
<tdclass="markdownTableBodyRight"><code>tau_wrt</code></td><tdclass="markdownTableBodyCenter">Logical</td><tdclass="markdownTableBodyLeft">Add the elastic stress components to the database</td></tr>
629
640
<trclass="markdownTableRowOdd">
630
-
<tdclass="markdownTableBodyRight"><code>schlieren_alpha(i)</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">Intensity of the numerical Schlieren computed via <code>alpha(i)</code></td></tr>
641
+
<tdclass="markdownTableBodyRight"><code>fd_order</code></td><tdclass="markdownTableBodyCenter">Integer</td><tdclass="markdownTableBodyLeft">Order of finite differences for computing the vorticity and the numerical Schlieren function [1,2,4]</td></tr>
631
642
<trclass="markdownTableRowEven">
632
-
<tdclass="markdownTableBodyRight"><code>probe_wrt</code></td><tdclass="markdownTableBodyCenter">Logical</td><tdclass="markdownTableBodyLeft">Write the flow chosen probes data files for each time step</td></tr>
643
+
<tdclass="markdownTableBodyRight"><code>schlieren_alpha(i)</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">Intensity of the numerical Schlieren computed via <code>alpha(i)</code></td></tr>
633
644
<trclass="markdownTableRowOdd">
634
-
<tdclass="markdownTableBodyRight"><code>num_probes</code></td><tdclass="markdownTableBodyCenter">Integer</td><tdclass="markdownTableBodyLeft">Number of probes </td></tr>
645
+
<tdclass="markdownTableBodyRight"><code>probe_wrt</code></td><tdclass="markdownTableBodyCenter">Logical</td><tdclass="markdownTableBodyLeft">Write the flow chosen probes data files for each time step</td></tr>
635
646
<trclass="markdownTableRowEven">
636
-
<tdclass="markdownTableBodyRight"><code>probe(i)%[x,y,z]</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">Coordinates of probe $i$</td></tr>
647
+
<tdclass="markdownTableBodyRight"><code>num_probes</code></td><tdclass="markdownTableBodyCenter">Integer</td><tdclass="markdownTableBodyLeft">Number of probes</td></tr>
637
648
<trclass="markdownTableRowOdd">
638
-
<tdclass="markdownTableBodyRight"><code>output_partial_domain</code></td><tdclass="markdownTableBodyCenter">Logical</td><tdclass="markdownTableBodyLeft">Output part of the domain</td></tr>
649
+
<tdclass="markdownTableBodyRight"><code>probe(i)%[x,y,z]</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">Coordinates of probe $i$</td></tr>
639
650
<trclass="markdownTableRowEven">
640
-
<tdclass="markdownTableBodyRight"><code>[x,y,z]_outputbeg</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">Beginning of the output domain in the [x,y,z]-direction</td></tr>
651
+
<tdclass="markdownTableBodyRight"><code>output_partial_domain</code></td><tdclass="markdownTableBodyCenter">Logical</td><tdclass="markdownTableBodyLeft">Output part of the domain </td></tr>
641
652
<trclass="markdownTableRowOdd">
653
+
<tdclass="markdownTableBodyRight"><code>[x,y,z]_outputbeg</code></td><tdclass="markdownTableBodyCenter">Real </td><tdclass="markdownTableBodyLeft">Beginning of the output domain in the [x,y,z]-direction </td></tr>
654
+
<trclass="markdownTableRowEven">
642
655
<tdclass="markdownTableBodyRight"><code>[x,y,z]_outputend</code></td><tdclass="markdownTableBodyCenter">Real </td><tdclass="markdownTableBodyLeft">End of the output domain in the [x,y,z]-direction </td></tr>
643
656
</table>
644
657
<p>The table lists formatted database output parameters. The parameters define variables that are outputted from simulation and file types and formats of data as well as options for post-processing.</p>
0 commit comments