Docs @ 21a4ebe

Author: MFC Action

documentation/architecture.html

@@ -429,7 +429,7 @@ <h1 class="doxsection"><a class="anchor" id="autotoc_md13"></a>
 <li><b>Add the module to <span class="tt">docs/module_categories.json</span></b> so it appears in this page</li>
 </ol>
 <p>Follow the pattern of existing modules like <span class="tt">m_body_forces</span> (simple) or <span class="tt">m_viscous</span> (more involved) as a template.</p>
-<div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-20</div> </div></div><!-- contents -->
+<div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-21</div> </div></div><!-- contents -->
 </div><!-- PageDoc -->
 </div><!-- doc-content -->
 <div id="page-nav" class="page-nav-panel">

documentation/case_constraints.html

@@ -1500,7 +1500,7 @@ <h2 class="doxsection"><a class="anchor" id="physics-warnings"></a>
 </table>
 <hr  />
 <p>💡 <b>Tip:</b> If you encounter a validation error, check the relevant section above or review <a href="https://github.com/MFlowCode/MFC/blob/master/toolchain/mfc/case_validator.py"><span class="tt">case_validator.py</span></a> for complete validation logic.</p>
-<div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-20</div> </div></div><!-- contents -->
+<div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-21</div> </div></div><!-- contents -->
 </div><!-- PageDoc -->
 </div><!-- doc-content -->
 <div id="page-nav" class="page-nav-panel">

documentation/cli-reference.html

@@ -981,7 +981,7 @@ <h3 class="doxsection"><a class="anchor" id="autotoc_md80"></a>
 <h3 class="doxsection"><a class="anchor" id="autotoc_md81"></a>
 Debug Logging (<span class="tt">-d, --debug-log</span>)</h3>
 <p>Enables debug logging for the Python toolchain (mfc.sh internals). This is for troubleshooting the build system, not the MFC simulation code.</p>
-<div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-20</div> </div></div><!-- contents -->
+<div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-21</div> </div></div><!-- contents -->
 </div><!-- PageDoc -->
 </div><!-- doc-content -->
 <div id="page-nav" class="page-nav-panel">

documentation/equations.html

@@ -579,7 +579,7 @@ <h4 class="doxsection"><a class="anchor" id="autotoc_md159"></a>
 6.1.1 Method of Classes (Commander and Prosperetti <a class="el" href="citelist.html#CITEREF_commander89">[13]</a>; Ando et al. <a class="el" href="citelist.html#CITEREF_ando11">[2]</a>)</h4>
 <p><b>Modified mixture pressure:</b></p>
 <p class="formulaDsp">
-\[p = (1 - \alpha)\,p_l + \alpha\left(\frac{R^3\,p_{bw}}{\bar{R}^3} - \frac{\rho\,R^3\,\dot{R}^2}{\bar{R}^3}\right)\]
+\[p = (1 - \alpha)\,p_l + \alpha\left(\frac{R^3\,p_{bw}}{\bar{R}^3} + \frac{\rho\,R^3\,\dot{R}^2}{\bar{R}^3}\right)\]
 </p>
 <p><b>Modified stiffened gas for the liquid phase:</b></p>
 <p class="formulaDsp">

documentation/examples.html

@@ -630,7 +630,7 @@ <h3 class="doxsection"><a class="anchor" id="autotoc_md311"></a>
 <ul>
 <li>Papanastasiou, T. C. (1987). Flows of materials with yield. <em>J. Rheol.</em> 31, 385.</li>
 </ul>
-<div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-20</div> </div></div><!-- contents -->
+<div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-21</div> </div></div><!-- contents -->
 </div><!-- PageDoc -->
 </div><!-- doc-content -->
 <div id="page-nav" class="page-nav-panel">

documentation/parameters.html

@@ -1635,7 +1635,7 @@ <h2 class="doxsection"><a class="anchor" id="general"></a>
 <tr class="markdownTableRowOdd">
 <td class="markdownTableBodyNone"><span class="tt">relativity</span>  </td><td class="markdownTableBodyNone">Logical (T/F)  </td><td class="markdownTableBodyNone">Enable special relativity  </td><td class="markdownTableBodyNone"></td><td class="markdownTableBodyNone"><span class="tt">relativity</span> requires <span class="tt">mhd</span> to be enabled; HLLD is not available for RMHD (<span class="tt">relativity</span>)  </td></tr>
 <tr class="markdownTableRowEven">
-<td class="markdownTableBodyNone"><span class="tt">relax</span>  </td><td class="markdownTableBodyNone">Logical (T/F)  </td><td class="markdownTableBodyNone">Enable relaxation terms  </td><td class="markdownTableBodyNone"></td><td class="markdownTableBodyNone">phase change requires <span class="tt">model_eqns</span>==2 with <span class="tt">relax_model</span> in [5,6] or <span class="tt">model_eqns</span>==3 with <span class="tt">relax_model</span> in [1,4,5,6]; <span class="tt">palpha_eps</span> must be positive; <span class="tt">palpha_eps</span> must be less than 1  </td></tr>
+<td class="markdownTableBodyNone"><span class="tt">relax</span>  </td><td class="markdownTableBodyNone">Logical (T/F)  </td><td class="markdownTableBodyNone">Enable relaxation terms  </td><td class="markdownTableBodyNone"></td><td class="markdownTableBodyNone">phase change requires <span class="tt">num_fluids</span> &gt;= 2 (liquid = 1, vapor = 2); phase change requires <span class="tt">model_eqns</span>==2 with <span class="tt">relax_model</span> in [5,6] or <span class="tt">model_eqns</span>==3 with <span class="tt">relax_model</span> in [1,4,5,6]; <span class="tt">palpha_eps</span> must be positive  </td></tr>
 <tr class="markdownTableRowOdd">
 <td class="markdownTableBodyNone"><span class="tt">relax_model</span>  </td><td class="markdownTableBodyNone">Integer  </td><td class="markdownTableBodyNone">Relaxation model type  </td><td class="markdownTableBodyNone"></td><td class="markdownTableBodyNone">phase change requires <span class="tt">model_eqns</span>==2 with <span class="tt">relax_model</span> in [5,6] or <span class="tt">model_eqns</span>==3 with <span class="tt">relax_model</span> in [1,4,5,6]  </td></tr>
 <tr class="markdownTableRowEven">
@@ -1750,7 +1750,7 @@ <h2 class="doxsection"><a class="anchor" id="autotoc_md445"></a>
 <div class="line"> </div>
 <div class="line"># Filter by type</div>
 <div class="line">./mfc.sh params -t real weno</div>
-</div><!-- fragment --><div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-20</div> </div></div><!-- contents -->
+</div><!-- fragment --><div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-21</div> </div></div><!-- contents -->
 </div><!-- PageDoc -->
 </div><!-- doc-content -->
 <div id="page-nav" class="page-nav-panel">

documentation/physics_constraints.html

@@ -487,6 +487,7 @@ <h3 class="doxsection"><a class="anchor" id="autotoc_md481"></a>
 <p><b>Stage:</b> common, pre_process, simulation, post_process | <b>Severity:</b> error</p>
 <p><b>Enforced checks:</b></p>
 <ul>
+<li>phase change requires <span class="tt">num_fluids</span> &gt;= 2 (liquid = 1, vapor = 2)</li>
 <li>phase change requires <span class="tt">model_eqns</span>==2 with <span class="tt">relax_model</span> in [5,6] or <span class="tt">model_eqns</span>==3 with <span class="tt">relax_model</span> in [1,4,5,6]</li>
 <li><span class="tt">palpha_eps</span> must be positive</li>
 <li><span class="tt">palpha_eps</span> must be less than 1</li>
@@ -739,7 +740,7 @@ <h3 class="doxsection"><a class="anchor" id="autotoc_md504"></a>
 <li><code class="f90l">patch_icpp(i)%</code><code class="f90r">hcid</code>=hcid requires <span class="tt">files_dir</span> to be set</li>
 <li><code class="f90l">patch_icpp(i)%</code><code class="f90r">hcid</code>=hcid requires <span class="tt">file_extension</span> to be set</li>
 </ul>
-<div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-20</div> </div></div><!-- contents -->
+<div style="text-align:center; font-size:0.75rem; color:#888; padding:16px 0 0;">Page last updated: 2026-08-21</div> </div></div><!-- contents -->
 </div><!-- PageDoc -->
 </div><!-- doc-content -->
 <div id="page-nav" class="page-nav-panel">

simulation/doxygen_crawl.html

@@ -345,6 +345,7 @@
 <a href="namespacem__bubbles.html#aec95ed87199e71943f283ecbf199ca4c"/>
 <a href="namespacem__bubbles.html#afbe3b969f2151f7b5b615f12ee93901e"/>
 <a href="namespacem__bubbles__ee.html"/>
+<a href="namespacem__bubbles__ee.html#a0a2216eea31da1917887f93bbdcdc5e3"/>
 <a href="namespacem__bubbles__ee.html#a523f8112128cbb1d0ec8b2bab674475e"/>
 <a href="namespacem__bubbles__ee.html#a74139ae2462f5d45e1834e11d66f5f28"/>
 <a href="namespacem__bubbles__ee.html#a75e41b12f290decf36c8317a5df45799"/>

simulation/m__bubbles__EE_8fpp_8f90.html

@@ -178,6 +178,8 @@
 <tr class="memdesc:a95ca95b3f22f62885bf1dd09eda53f3a"><td class="mdescLeft">&#160;</td><td class="mdescRight">Compute the right-hand side for Euler-Euler bubble transport.  <br /></td></tr>
 <tr class="memitem:ab15ff78772125079eec3132ae9a5bea1" id="r_ab15ff78772125079eec3132ae9a5bea1"><td class="memItemLeft">impure subroutine&#160;</td><td class="memItemRight"><a class="el" href="namespacem__bubbles__ee.html#ab15ff78772125079eec3132ae9a5bea1">m_bubbles_ee::s_compute_bubble_ee_source</a> (<a class="el" href="m__phase__change_8fpp_8f90.html#a6d2af1366163bbea81771912d405abfa">q_cons_vf</a>, q_prim_vf, rhs_vf, divu_in)</td></tr>
 <tr class="memdesc:ab15ff78772125079eec3132ae9a5bea1"><td class="mdescLeft">&#160;</td><td class="mdescRight">Compute the Euler-Euler bubble source terms.  <br /></td></tr>
+<tr class="memitem:a0a2216eea31da1917887f93bbdcdc5e3" id="r_a0a2216eea31da1917887f93bbdcdc5e3"><td class="memItemLeft">impure subroutine&#160;</td><td class="memItemRight"><a class="el" href="namespacem__bubbles__ee.html#a0a2216eea31da1917887f93bbdcdc5e3">m_bubbles_ee::s_compute_ptilde</a> (<a class="el" href="m__phase__change_8fpp_8f90.html#a6d2af1366163bbea81771912d405abfa">q_cons_vf</a>, q_prim_vf)</td></tr>
+<tr class="memdesc:a0a2216eea31da1917887f93bbdcdc5e3"><td class="mdescLeft">&#160;</td><td class="mdescRight">Evaluate the Euler-Euler pressure correction at cell centers and store it in ptil. This mirrors the modified mixture pressure that the HLLC solver applies to the reconstructed face states in m_riemann_solver_hllc.fpp, following Ando (2010), Eq. (2.5). It is a diagnostic: ptil is read only by the probe output path and never feeds back into the solution.  <br /></td></tr>
 </table><table class="memberdecls">
 <tr class="heading"><td colspan="2"><h2 id="header-var-members" class="groupheader"><a id="var-members" name="var-members"></a>
 Variables</h2></td></tr>

simulation/m__bubbles__EE_8fpp_8f90.js

@@ -3,6 +3,7 @@ var m__bubbles__EE_8fpp_8f90 =
     [ "m_bubbles_ee::s_comp_alpha_from_n", "namespacem__bubbles__ee.html#a523f8112128cbb1d0ec8b2bab674475e", null ],
     [ "m_bubbles_ee::s_compute_bubble_ee_source", "namespacem__bubbles__ee.html#ab15ff78772125079eec3132ae9a5bea1", null ],
     [ "m_bubbles_ee::s_compute_bubbles_ee_rhs", "namespacem__bubbles__ee.html#a95ca95b3f22f62885bf1dd09eda53f3a", null ],
+    [ "m_bubbles_ee::s_compute_ptilde", "namespacem__bubbles__ee.html#a0a2216eea31da1917887f93bbdcdc5e3", null ],
     [ "m_bubbles_ee::s_initialize_bubbles_ee_module", "namespacem__bubbles__ee.html#acdf857ad6b9f389ee1bd4dd81bbc23a5", null ],
     [ "m_bubbles_ee::bub_adv_src", "namespacem__bubbles__ee.html#a75e41b12f290decf36c8317a5df45799", null ],
     [ "m_bubbles_ee::bub_m_src", "namespacem__bubbles__ee.html#ab3e7e32a515a15bd575b5d1a5583ab11", null ],