Docs @ 465bd4f

Author: MFC Action

documentation/md_case.html

@@ -762,9 +762,13 @@ <h2><a class="anchor" id="autotoc_md17"></a>
 <tr class="markdownTableRowOdd">
 <td class="markdownTableBodyRight"><code>perturb_sph_fluid</code>   </td><td class="markdownTableBodyCenter">Integer   </td><td class="markdownTableBodyLeft">Fluid component whose partial density is to be perturbed    </td></tr>
 <tr class="markdownTableRowEven">
-<td class="markdownTableBodyRight"><code>vel_profile</code>   </td><td class="markdownTableBodyCenter">Logical   </td><td class="markdownTableBodyLeft">Set the mean streamwise velocity to hyperbolic tangent profile    </td></tr>
+<td class="markdownTableBodyRight"><code>mixlayer_vel_profile</code>   </td><td class="markdownTableBodyCenter">Logical   </td><td class="markdownTableBodyLeft">Set the mean streamwise velocity to hyperbolic tangent profile    </td></tr>
 <tr class="markdownTableRowOdd">
-<td class="markdownTableBodyRight"><code>instability_wave</code>   </td><td class="markdownTableBodyCenter">Logical   </td><td class="markdownTableBodyLeft">Perturb the initial velocity field by instability waves   </td></tr>
+<td class="markdownTableBodyRight"><code>mixlayer_vel_coef</code>   </td><td class="markdownTableBodyCenter">Real   </td><td class="markdownTableBodyLeft">Coefficient for the hyperbolic tangent profile of a mixing layer    </td></tr>
+<tr class="markdownTableRowEven">
+<td class="markdownTableBodyRight"><code>mixlayer_perturb</code>   </td><td class="markdownTableBodyCenter">Logical   </td><td class="markdownTableBodyLeft">Perturb the initial velocity field by instability waves    </td></tr>
+<tr class="markdownTableRowOdd">
+<td class="markdownTableBodyRight"><code>mixlayer_domain</code>   </td><td class="markdownTableBodyCenter">Real   </td><td class="markdownTableBodyLeft">Domain size of a mixing layer for the linear stability analysis   </td></tr>
 </table>
 <p>The table lists velocity field parameters. The parameters are optionally used to define initial velocity profiles and perturbations.</p>
 <ul>
@@ -773,8 +777,13 @@ <h2><a class="anchor" id="autotoc_md17"></a>
 <li><code>perturb_flow</code> activates the perturbation of initial velocity by random noise.</li>
 <li><code>perturb_sph</code> activates the perturbation of initial partial density by random noise.</li>
 <li><code>perturb_sph_fluid</code> specifies the fluid component whose partial density is to be perturbed.</li>
-<li><code>vel_profile</code> activates setting the mean streamwise velocity to hyperbolic tangent profile. This option works only for 2D and 3D cases.</li>
-<li><code>instability_wave</code> activates the perturbation of initial velocity by instability waves obtained from linear stability analysis for a mixing layer with hyperbolic tangent mean streamwise velocity profile. This option only works for <code>n &gt; 0</code>, <code>bc_y%[beg,end] = -5</code>, and <code>vel_profile = 'T'</code>.</li>
+<li><code>mixlayer_vel_profile</code> activates setting of the mean streamwise velocity to hyperbolic tangent profile. This option works only for 2D and 3D cases.</li>
+<li><code>mixlayer_vel_coef</code> is a parameter for the hyperbolic tangent profile of a mixing layer when &lsquo;mixlayer_vel_profile = 'T&rsquo;`. The mean streamwise velocity profile is given as:</li>
+</ul>
+<p>$$ u = patch_icpp(1)%vel(1) * tanh(y_cc * mixlayer_vel_profile) $$</p>
+<ul>
+<li><code>mixlayer_perturb</code> activates the perturbation of initial velocity by instability waves obtained from linear stability analysis for a mixing layer with hyperbolic tangent mean streamwise velocity profile. This option only works for <code>n &gt; 0</code>, <code>bc_y%[beg,end] = -6</code>, <code>num_fluids = 1</code>, <code>model_eqns = 2</code> and &lsquo;mixlayer_vel_profile = 'T&rsquo;`.</li>
+<li><code>mixlayer_domain</code> defines the domain size to compute spatial eigenvalues of the linear instability analysis when &lsquo;mixlayer_perturb = 'T&rsquo;<code>. For example, the spatial eigenvalue in</code>x` direction in 2D problem will be $2 \pi \alpha / (mixlayer_domain*patch_icpp(1)%length_y)$ for $\alpha = 1$, $2$ and $4$.</li>
 </ul>
 <h2><a class="anchor" id="autotoc_md18"></a>
 11. Phase Change Model</h2>

post_process/m__eigen__solver_8f90.html

@@ -150,14 +150,19 @@
 <tr class="heading"><td colspan="2"><h2 class="groupheader"><a id="func-members" name="func-members"></a>
 Functions/Subroutines</h2></td></tr>
 <tr class="memitem:ae98c4a659de6724591ff45c1c72080b6" id="r_ae98c4a659de6724591ff45c1c72080b6"><td class="memItemLeft" align="right" valign="top">subroutine, public&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namespacem__eigen__solver.html#ae98c4a659de6724591ff45c1c72080b6">m_eigen_solver::cg</a> (nm, nl, ar, ai, wr, wi, zr, zi, fv1, fv2, fv3, ierr)</td></tr>
+<tr class="memdesc:ae98c4a659de6724591ff45c1c72080b6"><td class="mdescLeft">&#160;</td><td class="mdescRight">This subroutine calls the recommended sequence of subroutines from the eigensystem subroutine package (eispack) to find the eigenvalues and eigenvectors (if desired) of a complex general matrix.  <br /></td></tr>
 <tr class="separator:ae98c4a659de6724591ff45c1c72080b6"><td class="memSeparator" colspan="2">&#160;</td></tr>
 <tr class="memitem:a79d1a338e00b679ba5fcf166c529d77c" id="r_a79d1a338e00b679ba5fcf166c529d77c"><td class="memItemLeft" align="right" valign="top">subroutine, public&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namespacem__eigen__solver.html#a79d1a338e00b679ba5fcf166c529d77c">m_eigen_solver::cbal</a> (nm, nl, ar, ai, low, igh, scale)</td></tr>
+<tr class="memdesc:a79d1a338e00b679ba5fcf166c529d77c"><td class="mdescLeft">&#160;</td><td class="mdescRight">This subroutine is a translation of the algol procedure cbalance, which is a complex version of balance, num. math. 13, 293-304(1969) by parlett and reinsch. handbook for auto. comp., vol.ii-linear algebra, 315-326(1971). This subroutine balances a complex matrix and isolates eigenvalues whenever possible.  <br /></td></tr>
 <tr class="separator:a79d1a338e00b679ba5fcf166c529d77c"><td class="memSeparator" colspan="2">&#160;</td></tr>
 <tr class="memitem:ab54bfcf38d40c65fcd0887727d3b00f0" id="r_ab54bfcf38d40c65fcd0887727d3b00f0"><td class="memItemLeft" align="right" valign="top">subroutine, public&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namespacem__eigen__solver.html#ab54bfcf38d40c65fcd0887727d3b00f0">m_eigen_solver::corth</a> (nm, nl, low, igh, ar, ai, ortr, orti)</td></tr>
+<tr class="memdesc:ab54bfcf38d40c65fcd0887727d3b00f0"><td class="mdescLeft">&#160;</td><td class="mdescRight">This subroutine is a translation of a complex analogue of the algol procedure orthes, num. math. 12, 349-368(1968) by martin and wilkinson. handbook for auto. comp., vol.ii-linear algebra, 339-358(1971). Given a complex general matrix, this subroutine reduces a submatrix situated in rows and columns low through igh to upper hessenberg form by unitary similarity transformations.  <br /></td></tr>
 <tr class="separator:ab54bfcf38d40c65fcd0887727d3b00f0"><td class="memSeparator" colspan="2">&#160;</td></tr>
 <tr class="memitem:a41c2d326d4ab6e1c50656217cbc84b51" id="r_a41c2d326d4ab6e1c50656217cbc84b51"><td class="memItemLeft" align="right" valign="top">subroutine, public&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namespacem__eigen__solver.html#a41c2d326d4ab6e1c50656217cbc84b51">m_eigen_solver::comqr2</a> (nm, nl, low, igh, ortr, orti, hr, hi, wr, wi, zr, zi, ierr)</td></tr>
+<tr class="memdesc:a41c2d326d4ab6e1c50656217cbc84b51"><td class="mdescLeft">&#160;</td><td class="mdescRight">This subroutine is a translation of a unitary analogue of the algol procedure comlr2, num. math. 16, 181-204(1970) by peters and wilkinson. handbook for auto. comp., vol.ii-linear algebra, 372-395(1971). The unitary analogue substitutes the qr algorithm of francis (comp. jour. 4, 332-345(1962)) for the lr algorithm. This subroutine finds the eigenvalues and eigenvectors of a complex upper hessenberg matrix by the qr method. The eigenvectors of a complex general matrix can also be found if corth has been used to reduce this general matrix to hessenberg form.  <br /></td></tr>
 <tr class="separator:a41c2d326d4ab6e1c50656217cbc84b51"><td class="memSeparator" colspan="2">&#160;</td></tr>
 <tr class="memitem:a6d35cffb5e64987dbbb6211356e9136f" id="r_a6d35cffb5e64987dbbb6211356e9136f"><td class="memItemLeft" align="right" valign="top">subroutine&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namespacem__eigen__solver.html#a6d35cffb5e64987dbbb6211356e9136f">m_eigen_solver::cbabk2</a> (nm, nl, low, igh, scale, ml, zr, zi)</td></tr>
+<tr class="memdesc:a6d35cffb5e64987dbbb6211356e9136f"><td class="mdescLeft">&#160;</td><td class="mdescRight">This subroutine is a translation of the algol procedure cbabk2, which is a complex version of balbak, num. math. 13, 293-304(1969) by parlett and reinsch. handbook for auto. comp., vol.ii-linear algebra, 315-326(1971). This subroutine forms the eigenvectors of a complex general matrix by back transforming those of the correspondingbalanced matrix determined by cbal.  <br /></td></tr>
 <tr class="separator:a6d35cffb5e64987dbbb6211356e9136f"><td class="memSeparator" colspan="2">&#160;</td></tr>
 <tr class="memitem:a26981f13d48d92b739dfb0c30189417f" id="r_a26981f13d48d92b739dfb0c30189417f"><td class="memItemLeft" align="right" valign="top">subroutine, public&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namespacem__eigen__solver.html#a26981f13d48d92b739dfb0c30189417f">m_eigen_solver::csroot</a> (xr, xi, yr, yi)</td></tr>
 <tr class="separator:a26981f13d48d92b739dfb0c30189417f"><td class="memSeparator" colspan="2">&#160;</td></tr>