6.3. MMS04 manufactured solutions
The manufactured solutions for the No-slip wall BC, single phase, 3D, curl-based functions are listed below.
Gas pressure:
(6.30)\[\begin{split}p_{g} = p_{g0} &+ p_{\text{gx}}\cos\left( A_{p_{\text{gx}}}\text{πx} \right) + p_{\text{gy}}\cos\left( A_{p_{\text{gy}}}\text{πy} \right) + p_{\text{gxy}}\cos\left( A_{p_{\text{gxy}}}\text{πxy} \right) \\
&+ p_{\text{gz}}\sin\left( A_{p_{\text{gz}}}\text{πz} \right) + p_{\text{gyz}}\sin\left( A_{p_{\text{gyz}}}\text{πyz} \right) + p_{\text{gzx}}\cos\left( A_{p_{\text{gzx}}}\text{πzx} \right)\end{split}\]
Gas velocity components:
(6.31)\[\begin{split}u_{g} = x^{2}\big[& A_{w_{\text{gy}}}\pi w_{\text{gy}}\cos\left( A_{w_{\text{gy}}}\text{πy} \right) + A_{w_{\text{gxy}}}\pi w_{\text{gxy}}x\cos\left( A_{w_{\text{gxy}}}\text{πxy} \right) \\
&- A_{v_{\text{gyz}}}\pi v_{\text{gyz}}y\cos\left( A_{v_{\text{gyz}}}\text{πyz} \right) + A_{w_{\text{gyz}}}\pi w_{\text{gyz}}z\cos\left( A_{w_{\text{gyz}}}\text{πyz} \right) \\
&+ A_{v_{\text{gz}}}\pi v_{\text{gz}}\sin\left( A_{v_{\text{gz}}}\text{πz} \right) + A_{v_{\text{gzx}}}\pi v_{\text{gzx}}x\sin\left( A_{v_{\text{gzx}}}\text{πzx} \right) \big]\end{split}\]
(6.32)\[\begin{split}v_{g} = x^{2}\big[& - A_{w_{\text{gxy}}}\pi w_{\text{gxy}}y\cos\left( A_{w_{\text{gxy}}}\text{πxy} \right) + A_{u_{\text{gyz}}}\pi u_{\text{gyz}}y\cos\left( A_{u_{\text{gyz}}}\text{πyz} \right) \\
&+ A_{w_{\text{gx}}}\pi w_{\text{gx}}\sin\left( A_{w_{\text{gx}}}\text{πx} \right) - A_{u_{\text{gz}}}\pi u_{\text{gz}}\sin\left( A_{u_{\text{gz}}}\text{πz} \right) \\
&- A_{u_{\text{gzx}}}\pi u_{\text{gzx}}x\sin\left( A_{u_{\text{gzx}}}\text{πzx} \right) + A_{w_{\text{gzx}}}\pi w_{\text{gzx}}z\sin\left( A_{w_{\text{gzx}}}\text{πzx} \right) \big] \\
&+ 2x\big[ - w_{g0} - w_{\text{gx}}\cos\left( A_{w_{\text{gx}}}\text{πx} \right) - w_{\text{gz}}\cos\left( A_{w_{\text{gz}}}\text{πz} \right) \\
&- w_{\text{gzx}}\cos\left( A_{w_{\text{gzx}}}\text{πzx} \right) - w_{\text{gy}}\sin\left( A_{w_{\text{gy}}}\text{πy} \right) \\
&- w_{\text{gxy}}\sin\left( A_{w_{\text{gxy}}}\text{πxy} \right) - w_{\text{gyz}}\sin\left( A_{w_{\text{gyz}}}\text{πyz} \right)\big]\end{split}\]
(6.33)\[\begin{split}w_{g} = x^{2}\big[& A_{v_{\text{gx}}}\pi v_{\text{gx}}\cos\left( A_{v_{\text{gx}}}\text{πx} \right) - A_{u_{\text{gyz}}}\pi u_{\text{gyz}}z\cos\left( A_{u_{\text{gyz}}}\text{πyz} \right) \\
&+ A_{u_{\text{gy}}}\pi u_{\text{gy}}\sin\left( A_{u_{\text{gy}}}\text{πy} \right) + A_{u_{\text{gxy}}}\pi u_{\text{gxy}}x\sin\left( A_{u_{\text{gxy}}}\text{πxy} \right) \\
&- A_{v_{\text{gxy}}}\pi v_{\text{gxy}}y\sin\left( A_{v_{\text{gxy}}}\text{πxy} \right) - A_{v_{\text{gzx}}}\pi v_{\text{gzx}}z\sin\left( A_{v_{\text{gzx}}}\text{πzx} \right) \big] \\
&+ 2x\big[ v_{g0} + v_{\text{gy}}\cos\left( A_{v_{\text{gy}}}\text{πy} \right) + v_{\text{gxy}}\cos\left( A_{v_{\text{gxy}}}\text{πxy} \right) \\
&+ v_{\text{gz}}\cos\left( A_{v_{\text{gz}}}\text{πz} \right) + v_{\text{gzx}}\cos\left( A_{v_{\text{gzx}}}\text{πzx} \right) \\
&+ v_{\text{gx}}\sin\left( A_{v_{\text{gx}}}\text{πx} \right) + v_{\text{gyz}}\sin\left( A_{v_{\text{gyz}}}\text{πyz} \right) \big]\end{split}\]
Gas volume fraction:
(6.34)\[\varepsilon_{g} = 1.0\]
The parameters appearing in the manufactured solutions are as follows:
Table 6.6 Parameters in MMS04 manufactured solutions.
\(p_{g0}\) |
100.0 |
\(v_{g0}\) |
9.0 |
\(u_{g0}\) |
7.0 |
\(w_{g0}\) |
8.0 |
\(p_{\text{gx}}\) |
20.0 |
\(v_{\text{gx}}\) |
-5.0 |
\(u_{\text{gx}}\) |
3.0 |
\(w_{\text{gx}}\) |
-4.0 |
\(p_{\text{gy}}\) |
-50.0 |
\(v_{\text{gy}}\) |
4.0 |
\(u_{\text{gy}}\) |
-4.0 |
\(w_{\text{gy}}\) |
3.5 |
\(p_{\text{gz}}\) |
20.0 |
\(v_{\text{gz}}\) |
5.0 |
\(u_{\text{gz}}\) |
-3.0 |
\(w_{\text{gz}}\) |
4.2 |
\(p_{\text{gxy}}\) |
-25.0 |
\(v_{\text{gxy}}\) |
-3.0 |
\(u_{\text{gxy}}\) |
2.0 |
\(w_{\text{gxy}}\) |
-2.2 |
\(p_{\text{gyz}}\) |
-10.0 |
\(v_{\text{gyz}}\) |
2.5 |
\(u_{\text{gyz}}\) |
1.5 |
\(w_{\text{gyz}}\) |
2.1 |
\(p_{\text{gzx}}\) |
10.0 |
\(v_{\text{gzx}}\) |
3.5 |
\(u_{\text{gzx}}\) |
-2.0 |
\(w_{\text{gzx}}\) |
2.5 |
\(A_{p_{\text{gx}}}\) |
0.4 |
\(A_{v_{\text{gx}}}\) |
0.8 |
\(A_{u_{\text{gx}}}\) |
0.5 |
\(A_{w_{\text{gx}}}\) |
0.85 |
\(A_{p_{\text{gy}}}\) |
0.45 |
\(A_{v_{\text{gy}}}\) |
0.8 |
\(A_{u_{\text{gy}}}\) |
0.85 |
\(A_{w_{\text{gy}}}\) |
0.9 |
\(A_{p_{\text{gz}}}\) |
0.85 |
\(A_{v_{\text{gz}}}\) |
0.5 |
\(A_{u_{\text{gz}}}\) |
0.4 |
\(A_{w_{\text{gz}}}\) |
0.5 |
\(A_{p_{\text{gxy}}}\) |
0.75 |
\(A_{v_{\text{gxy}}}\) |
0.9 |
\(A_{u_{\text{gxy}}}\) |
0.6 |
\(A_{w_{\text{gxy}}}\) |
0.4 |
\(A_{p_{\text{gyz}}}\) |
0.7 |
\(A_{v_{\text{gyz}}}\) |
0.4 |
\(A_{u_{\text{gyz}}}\) |
0.8 |
\(A_{w_{\text{gyz}}}\) |
0.8 |
\(A_{p_{\text{gzx}}}\) |
0.8 |
\(A_{v_{\text{gzx}}}\) |
0.6 |
\(A_{u_{\text{gzx}}}\) |
0.9 |
\(A_{w_{\text{gzx}}}\) |
0.75 |