| Line 36: | Line 36: | ||
|- | |- | ||
| <math>\boldsymbol{U}=\left\{ \begin{matrix}\rho \\ \rho {v}_{1} \\ \rho {v}_{2} \\ \rho {v}_{3} \\ \rho e \\ \rho \tilde{\nu }\end{matrix}\right\} ;{\boldsymbol{F}}_{i=}\left\{ \begin{matrix}\rho {v}_{i} \\ \rho {v}_{1}{v}_{i}+p{\delta }_{i1} \\ \rho {v}_{2}{v}_{i}+p{\delta }_{i2} \\ \rho {v}_{3}{v}_{i}+p{\delta }_{i3} \\ {v}_{i}\left( \rho e+p\right) \\ \rho \tilde{\nu }{v}_{i}\end{matrix}\right\} ;{\boldsymbol{G}}_{i}=</math><math>\left\{ \begin{matrix} 0 \\ -{\tau }_{i1} \\ -{\tau }_{i2} \\ -{\tau }_{i3} \\ -{\tau }_{ij}{v}_{j}-{q}_{i}-{\psi }_{i} \\ {\varphi }_{i}\end{matrix}\right\} ;\boldsymbol{Q}=</math><math>\left\{ \begin{matrix} 0 \\ 0 \\ 0 \\ 0 \\ 0 \\P+D+{D}_{T}\end{matrix}\right\}</math> | | <math>\boldsymbol{U}=\left\{ \begin{matrix}\rho \\ \rho {v}_{1} \\ \rho {v}_{2} \\ \rho {v}_{3} \\ \rho e \\ \rho \tilde{\nu }\end{matrix}\right\} ;{\boldsymbol{F}}_{i=}\left\{ \begin{matrix}\rho {v}_{i} \\ \rho {v}_{1}{v}_{i}+p{\delta }_{i1} \\ \rho {v}_{2}{v}_{i}+p{\delta }_{i2} \\ \rho {v}_{3}{v}_{i}+p{\delta }_{i3} \\ {v}_{i}\left( \rho e+p\right) \\ \rho \tilde{\nu }{v}_{i}\end{matrix}\right\} ;{\boldsymbol{G}}_{i}=</math><math>\left\{ \begin{matrix} 0 \\ -{\tau }_{i1} \\ -{\tau }_{i2} \\ -{\tau }_{i3} \\ -{\tau }_{ij}{v}_{j}-{q}_{i}-{\psi }_{i} \\ {\varphi }_{i}\end{matrix}\right\} ;\boldsymbol{Q}=</math><math>\left\{ \begin{matrix} 0 \\ 0 \\ 0 \\ 0 \\ 0 \\P+D+{D}_{T}\end{matrix}\right\}</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(2) |
|} | |} | ||
| Line 74: | Line 74: | ||
|- | |- | ||
| <math>{\tau }_{ij}=\left( \mu +{\mu }_{t}\right) \left[ \frac{\partial {v}_{i}}{\partial {x}_{j}}+\right. </math><math>\left. \frac{\partial {v}_{j}}{\partial {x}_{i}}\right] -\frac{2}{3}\left( \mu +\right. </math><math>\left. {\mu }_{t}\right) \frac{\partial {v}_{k}}{\partial {x}_{k}}{\delta }_{ij}-</math><math>\frac{2}{3}\rho k{\delta }_{ij}</math> | | <math>{\tau }_{ij}=\left( \mu +{\mu }_{t}\right) \left[ \frac{\partial {v}_{i}}{\partial {x}_{j}}+\right. </math><math>\left. \frac{\partial {v}_{j}}{\partial {x}_{i}}\right] -\frac{2}{3}\left( \mu +\right. </math><math>\left. {\mu }_{t}\right) \frac{\partial {v}_{k}}{\partial {x}_{k}}{\delta }_{ij}-</math><math>\frac{2}{3}\rho k{\delta }_{ij}</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(3) |
|} | |} | ||
| Line 84: | Line 84: | ||
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" | ||
| <math>{q}_{i}=\left( K+{K}_{t}\right) \frac{\partial u}{\partial {x}_{i}}</math> | | <math>{q}_{i}=\left( K+{K}_{t}\right) \frac{\partial u}{\partial {x}_{i}}</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(4) |
|} | |} | ||
| Line 94: | Line 94: | ||
|- | |- | ||
| <math>{\psi }_{i}=\left( \mu +\frac{{\mu }_{t}}{{\sigma }_{k}}\right) \frac{\partial k}{\partial {x}_{i}}</math> | | <math>{\psi }_{i}=\left( \mu +\frac{{\mu }_{t}}{{\sigma }_{k}}\right) \frac{\partial k}{\partial {x}_{i}}</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(5) |
|} | |} | ||
| Line 105: | Line 105: | ||
|- | |- | ||
| <math>{\varphi }_{i}=\frac{1}{\sigma }\left( \mu +\rho \tilde{\nu }\right) \frac{\partial \tilde{\nu }}{\partial {x}_{i}}</math> | | <math>{\varphi }_{i}=\frac{1}{\sigma }\left( \mu +\rho \tilde{\nu }\right) \frac{\partial \tilde{\nu }}{\partial {x}_{i}}</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(6) |
|} | |} | ||
| Line 114: | Line 114: | ||
|- | |- | ||
| <math>{D}_{T}=-\frac{1}{\sigma }\, \rho \, {cb}_{2}\frac{\partial \tilde{\nu }}{\partial {x}_{i}}\frac{\partial \tilde{\nu }}{\partial {x}_{i}}+</math><math>\, \frac{1}{\sigma }(\nu +\tilde{\nu })\, \frac{\partial \rho }{\partial {x}_{i}}\frac{\partial \tilde{\nu }}{\partial {x}_{i}}</math> | | <math>{D}_{T}=-\frac{1}{\sigma }\, \rho \, {cb}_{2}\frac{\partial \tilde{\nu }}{\partial {x}_{i}}\frac{\partial \tilde{\nu }}{\partial {x}_{i}}+</math><math>\, \frac{1}{\sigma }(\nu +\tilde{\nu })\, \frac{\partial \rho }{\partial {x}_{i}}\frac{\partial \tilde{\nu }}{\partial {x}_{i}}</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(7) |
|} | |} | ||
| Line 124: | Line 124: | ||
|- | |- | ||
| <math>P=-{cb}_{1}\, S\, \rho \tilde{\nu };\, \, D={cw}_{1}\, {f}_{w}\, \rho \, {\left( \frac{\tilde{\nu }}{d}\right) }^{2}</math> | | <math>P=-{cb}_{1}\, S\, \rho \tilde{\nu };\, \, D={cw}_{1}\, {f}_{w}\, \rho \, {\left( \frac{\tilde{\nu }}{d}\right) }^{2}</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(8) |
|} | |} | ||
| Line 136: | Line 136: | ||
<math display="inline">S=\Omega +\, {f}_{v2}\, \frac{\tilde{\nu }}{{{\Kappa }^{2}d}^{2}}\, ;\, \Kappa =</math><math>0.41</math> <math display="inline"></math> <math>\boldsymbol{\Omega }=\, \, \frac{1}{2}\, rot\mathit{\boldsymbol{V;\, }}\Omega =\sqrt{\boldsymbol{\Omega }\cdot \boldsymbol{\Omega }}</math> | <math display="inline">S=\Omega +\, {f}_{v2}\, \frac{\tilde{\nu }}{{{\Kappa }^{2}d}^{2}}\, ;\, \Kappa =</math><math>0.41</math> <math display="inline"></math> <math>\boldsymbol{\Omega }=\, \, \frac{1}{2}\, rot\mathit{\boldsymbol{V;\, }}\Omega =\sqrt{\boldsymbol{\Omega }\cdot \boldsymbol{\Omega }}</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(9) |
|} | |} | ||
| Line 148: | Line 148: | ||
<math>r=\frac{\tilde{\nu }}{{{S\, \Kappa }^{2}d}^{2}}\, ;\, {c}_{w1}=\frac{{\, cb}_{1}}{{\Kappa }^{2}}+</math><math>\, \frac{1+{\, cb}_{2}}{\sigma }</math> | <math>r=\frac{\tilde{\nu }}{{{S\, \Kappa }^{2}d}^{2}}\, ;\, {c}_{w1}=\frac{{\, cb}_{1}}{{\Kappa }^{2}}+</math><math>\, \frac{1+{\, cb}_{2}}{\sigma }</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(10) |
|} | |} | ||
| Line 156: | Line 156: | ||
|- | |- | ||
| <math>{\mu }_{t}=\rho \, \tilde{\nu }\, {f}_{v1};\, {K}_{t}={\mu }_{t}\frac{\gamma }{{Pr}_{t}};\, \, \gamma =</math><math>\frac{{c}_{p}}{{c}_{v}}=1.4</math> | | <math>{\mu }_{t}=\rho \, \tilde{\nu }\, {f}_{v1};\, {K}_{t}={\mu }_{t}\frac{\gamma }{{Pr}_{t}};\, \, \gamma =</math><math>\frac{{c}_{p}}{{c}_{v}}=1.4</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(11) |
|} | |} | ||
| Line 167: | Line 167: | ||
|- | |- | ||
| <math>p=\left( \gamma -1\right) \rho u;</math> | | <math>p=\left( \gamma -1\right) \rho u;</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(12) |
|} | |} | ||
| Line 176: | Line 176: | ||
|- | |- | ||
| <math>u={c}_{v}\, T=e-\frac{1}{2}{v}_{i}{v}_{i}-k</math> | | <math>u={c}_{v}\, T=e-\frac{1}{2}{v}_{i}{v}_{i}-k</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(13) |
|} | |} | ||
| Line 186: | Line 186: | ||
|- | |- | ||
| <math>\mu ={\mu }_{\infty }\frac{S+{u}_{\infty }}{S+u}{\left( \frac{u}{{u}_{\infty }}\right) }^{\frac{3}{2}};\quad \quad K=</math><math>{K}_{\infty }\, \frac{{S}_{k}+{u}_{\infty }}{{S}_{k}+u}{\left( \frac{u}{{u}_{\infty }}\right) }^{\frac{3}{2}}</math> | | <math>\mu ={\mu }_{\infty }\frac{S+{u}_{\infty }}{S+u}{\left( \frac{u}{{u}_{\infty }}\right) }^{\frac{3}{2}};\quad \quad K=</math><math>{K}_{\infty }\, \frac{{S}_{k}+{u}_{\infty }}{{S}_{k}+u}{\left( \frac{u}{{u}_{\infty }}\right) }^{\frac{3}{2}}</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(14) |
|} | |} | ||
| Line 197: | Line 197: | ||
|- | |- | ||
| <math>\begin{matrix}\begin{matrix}{v}_{i}=\overline{{v}_{i}}\quad on\quad {\Gamma }_{v}\\\rho =\overline{\rho }\quad \, \, on\quad {\Gamma }_{\rho }\\u=\overline{u}\quad \, \, on\quad {\Gamma }_{u}\end{matrix}\\\tilde{\nu }=\overline{\tilde{\nu }}\quad \quad on\quad {\Gamma }_{\tilde{\nu }}\end{matrix}</math> | | <math>\begin{matrix}\begin{matrix}{v}_{i}=\overline{{v}_{i}}\quad on\quad {\Gamma }_{v}\\\rho =\overline{\rho }\quad \, \, on\quad {\Gamma }_{\rho }\\u=\overline{u}\quad \, \, on\quad {\Gamma }_{u}\end{matrix}\\\tilde{\nu }=\overline{\tilde{\nu }}\quad \quad on\quad {\Gamma }_{\tilde{\nu }}\end{matrix}</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(15) |
|} | |} | ||
| Line 209: | Line 209: | ||
<math display="inline">{{q}_{i}\, n}_{i}\, =\hat{q}\quad on\quad {\Gamma }_{q}</math> | <math display="inline">{{q}_{i}\, n}_{i}\, =\hat{q}\quad on\quad {\Gamma }_{q}</math> | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(16) |
|} | |} | ||
| Line 264: | Line 264: | ||
| <math>\begin{matrix}\Delta {\left( \rho \right) }_{I+1}^{n+1}=\Delta t\left\{ -\frac{{\partial \left( {\rho v}_{i}\right) }^{n}}{{\partial x}_{i}}+\frac{\Delta t}{2}\left[ f\frac{\partial }{{\partial x}_{i}}\left( \frac{{\partial p}^{n}}{{\partial x}_{i}}\right) +\left( 1-f\right) \frac{\partial }{\partial {x}_{k}}\left( {{v}_{k}}^{n}\, \frac{{\partial \left( {\rho v}_{i}\right) }^{n}}{\partial {x}_{i}}\right) \right] \right\} +\\\quad \quad \quad \quad \quad \quad \quad \quad \quad +\frac{\Delta t}{2}\left\{ -\frac{\partial \Delta {\left( {\rho v}_{i}\right) }_{I}^{n+1}}{{\partial x}_{i}}+\frac{\Delta t}{2}\left[ f\frac{\partial }{{\partial x}_{i}}\left( \frac{{\partial \Delta p}_{I}^{n+1}}{{\partial x}_{i}}\right) +\left( 1-f\right) \frac{\partial }{\partial {x}_{k}}\left( {{v}_{k}}^{n}\, \frac{{\partial \left( \Delta {\rho v}_{i}\right) }_{I}^{n+1}}{\partial {x}_{i}}\right) \right] \right\} \end{matrix}\,</math> | | <math>\begin{matrix}\Delta {\left( \rho \right) }_{I+1}^{n+1}=\Delta t\left\{ -\frac{{\partial \left( {\rho v}_{i}\right) }^{n}}{{\partial x}_{i}}+\frac{\Delta t}{2}\left[ f\frac{\partial }{{\partial x}_{i}}\left( \frac{{\partial p}^{n}}{{\partial x}_{i}}\right) +\left( 1-f\right) \frac{\partial }{\partial {x}_{k}}\left( {{v}_{k}}^{n}\, \frac{{\partial \left( {\rho v}_{i}\right) }^{n}}{\partial {x}_{i}}\right) \right] \right\} +\\\quad \quad \quad \quad \quad \quad \quad \quad \quad +\frac{\Delta t}{2}\left\{ -\frac{\partial \Delta {\left( {\rho v}_{i}\right) }_{I}^{n+1}}{{\partial x}_{i}}+\frac{\Delta t}{2}\left[ f\frac{\partial }{{\partial x}_{i}}\left( \frac{{\partial \Delta p}_{I}^{n+1}}{{\partial x}_{i}}\right) +\left( 1-f\right) \frac{\partial }{\partial {x}_{k}}\left( {{v}_{k}}^{n}\, \frac{{\partial \left( \Delta {\rho v}_{i}\right) }_{I}^{n+1}}{\partial {x}_{i}}\right) \right] \right\} \end{matrix}\,</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(17) |
|} | |} | ||
| Line 276: | Line 276: | ||
| <math>f=\left\{ \begin{matrix}\, \, 1\quad \quad \quad \quad \, if\quad \quad \, \quad M\leq 1\, \\-M+2\quad \quad \, \, if\quad \, 1<M<2\quad \, \, \\\, 0\quad \quad \quad \, \, \, \, \, if\quad \quad \quad \, M\geq 2\end{matrix}\, \right.</math> | | <math>f=\left\{ \begin{matrix}\, \, 1\quad \quad \quad \quad \, if\quad \quad \, \quad M\leq 1\, \\-M+2\quad \quad \, \, if\quad \, 1<M<2\quad \, \, \\\, 0\quad \quad \quad \, \, \, \, \, if\quad \quad \quad \, M\geq 2\end{matrix}\, \right.</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(18) |
|} | |} | ||
| Line 295: | Line 295: | ||
|<math>\begin{matrix}{\left. \boldsymbol{\Delta \rho }\right| }_{I+1}^{n+1}=f\left[ {\Delta t}_{ext}\, {{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{B}}_{i}\, {\left. {\boldsymbol{f}}_{i}^{\rho }\right| }^{n}-\boldsymbol{P}{\boldsymbol{\, p}}^{n}-{\boldsymbol{r}}^{n}\right) +\frac{{\Delta t}_{ext}}{2}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -\, {\boldsymbol{B}}_{i}\, {\left. {\boldsymbol{\Delta f}}_{i}^{\rho }\right| }_{I}^{n+1}-\boldsymbol{P}\, {\boldsymbol{\Delta p}}_{I}^{n+1}\right) \right] +\\\quad \quad \quad +\left( 1-f\right) \left[ {\Delta t}_{ext}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{f}}_{i}^{\rho }\right| }^{n}\right) +\frac{{\Delta t}_{ext}}{2}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -\, {\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{\Delta f}}_{i}^{\rho }\right| }_{I}^{n+1}\right) \right] \end{matrix}</math> | |<math>\begin{matrix}{\left. \boldsymbol{\Delta \rho }\right| }_{I+1}^{n+1}=f\left[ {\Delta t}_{ext}\, {{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{B}}_{i}\, {\left. {\boldsymbol{f}}_{i}^{\rho }\right| }^{n}-\boldsymbol{P}{\boldsymbol{\, p}}^{n}-{\boldsymbol{r}}^{n}\right) +\frac{{\Delta t}_{ext}}{2}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -\, {\boldsymbol{B}}_{i}\, {\left. {\boldsymbol{\Delta f}}_{i}^{\rho }\right| }_{I}^{n+1}-\boldsymbol{P}\, {\boldsymbol{\Delta p}}_{I}^{n+1}\right) \right] +\\\quad \quad \quad +\left( 1-f\right) \left[ {\Delta t}_{ext}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{f}}_{i}^{\rho }\right| }^{n}\right) +\frac{{\Delta t}_{ext}}{2}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -\, {\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{\Delta f}}_{i}^{\rho }\right| }_{I}^{n+1}\right) \right] \end{matrix}</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(19a) |
|} | |} | ||
| Line 307: | Line 307: | ||
|<math>\quad \begin{matrix}{\left. \boldsymbol{\Delta \rho }{\boldsymbol{v}}_{j}\right| }_{I+1}^{n+1}={\Delta t}_{ext}\, {{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -\, {\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{f}}_{ij}^{\rho v}\right| }^{n}-{\boldsymbol{D}}_{ij}{{\boldsymbol{\, v}}_{i}}^{n}+{{\boldsymbol{t}}_{j}}^{n}\right) +\\\quad \quad \quad \quad \, +\frac{{\Delta t}_{ext}}{2}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, \boldsymbol{\Delta }{\left. {\boldsymbol{f}}_{ij}^{\rho v}\right| }_{I}^{n+1}-{\boldsymbol{D}}_{ij}{\left. {\boldsymbol{\, \Delta v}}_{i}\right| }_{I}^{n+1}\right) \end{matrix}</math> | |<math>\quad \begin{matrix}{\left. \boldsymbol{\Delta \rho }{\boldsymbol{v}}_{j}\right| }_{I+1}^{n+1}={\Delta t}_{ext}\, {{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -\, {\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{f}}_{ij}^{\rho v}\right| }^{n}-{\boldsymbol{D}}_{ij}{{\boldsymbol{\, v}}_{i}}^{n}+{{\boldsymbol{t}}_{j}}^{n}\right) +\\\quad \quad \quad \quad \, +\frac{{\Delta t}_{ext}}{2}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, \boldsymbol{\Delta }{\left. {\boldsymbol{f}}_{ij}^{\rho v}\right| }_{I}^{n+1}-{\boldsymbol{D}}_{ij}{\left. {\boldsymbol{\, \Delta v}}_{i}\right| }_{I}^{n+1}\right) \end{matrix}</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(19b) |
|} | |} | ||
| Line 319: | Line 319: | ||
|<math>\, \begin{matrix}{\left. \boldsymbol{\Delta \rho e}\right| }_{I+1}^{n+1}={\Delta t}_{ext}{{\boldsymbol{\, M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{f}}_{i}^{\rho e}\right| }^{n}-{\boldsymbol{E}}_{i}\, {{\boldsymbol{v}}_{i}}^{n}-\boldsymbol{K}\, {\boldsymbol{u}}^{n}+{\boldsymbol{h}}^{n}\right) +\\\quad \quad \quad \quad \quad \quad +\frac{{\Delta t}_{ext}}{2}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{\Delta f}}_{i}^{\rho e}\right| }_{I}^{n+1}-{\boldsymbol{E}}_{i}\, {\left. {\boldsymbol{\Delta v}}_{i}\right| }_{I}^{n+1}-\boldsymbol{K\, }{\left. \boldsymbol{\Delta u}\right| }_{I}^{n+1}\right) \end{matrix}</math> | |<math>\, \begin{matrix}{\left. \boldsymbol{\Delta \rho e}\right| }_{I+1}^{n+1}={\Delta t}_{ext}{{\boldsymbol{\, M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{f}}_{i}^{\rho e}\right| }^{n}-{\boldsymbol{E}}_{i}\, {{\boldsymbol{v}}_{i}}^{n}-\boldsymbol{K}\, {\boldsymbol{u}}^{n}+{\boldsymbol{h}}^{n}\right) +\\\quad \quad \quad \quad \quad \quad +\frac{{\Delta t}_{ext}}{2}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{\Delta f}}_{i}^{\rho e}\right| }_{I}^{n+1}-{\boldsymbol{E}}_{i}\, {\left. {\boldsymbol{\Delta v}}_{i}\right| }_{I}^{n+1}-\boldsymbol{K\, }{\left. \boldsymbol{\Delta u}\right| }_{I}^{n+1}\right) \end{matrix}</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(19c) |
|} | |} | ||
| Line 331: | Line 331: | ||
|<math>\begin{matrix}{\left. \boldsymbol{\Delta \rho }\tilde{\boldsymbol{\nu }}\right| }_{I+1}^{n+1}=\, {\Delta t}_{ext}{{\boldsymbol{\, M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{f}}_{i}^{\rho \tilde{\nu }}\right| }^{n}\boldsymbol{-M}\, {\boldsymbol{q}}^{n}-{\boldsymbol{A}}_{i}\, {\left. {\boldsymbol{q}}_{{\mathit{\boldsymbol{v}}}_{\mathit{\boldsymbol{i}}}}\right| }^{n}-{\boldsymbol{D}}_{\nu }\, {\tilde{\boldsymbol{\nu }}}^{n}\right) +\\\quad \quad \quad \quad \quad +\, \frac{{\Delta t}_{ext}}{2}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{\Delta f}}_{i}^{\rho \tilde{\nu }}\right| }_{I}^{n+1}-\boldsymbol{M}\, {\boldsymbol{\Delta q}}_{I}^{n+1}-\, {\boldsymbol{D}}_{\nu }\, \boldsymbol{\Delta }{\tilde{\boldsymbol{\nu }}}_{I}^{n+1}\right) \end{matrix}\quad \, \,</math> | |<math>\begin{matrix}{\left. \boldsymbol{\Delta \rho }\tilde{\boldsymbol{\nu }}\right| }_{I+1}^{n+1}=\, {\Delta t}_{ext}{{\boldsymbol{\, M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{f}}_{i}^{\rho \tilde{\nu }}\right| }^{n}\boldsymbol{-M}\, {\boldsymbol{q}}^{n}-{\boldsymbol{A}}_{i}\, {\left. {\boldsymbol{q}}_{{\mathit{\boldsymbol{v}}}_{\mathit{\boldsymbol{i}}}}\right| }^{n}-{\boldsymbol{D}}_{\nu }\, {\tilde{\boldsymbol{\nu }}}^{n}\right) +\\\quad \quad \quad \quad \quad +\, \frac{{\Delta t}_{ext}}{2}{{\boldsymbol{M}}_{\boldsymbol{D}}}^{-1}\left( -{\boldsymbol{H}}_{i}\, {\left. {\boldsymbol{\Delta f}}_{i}^{\rho \tilde{\nu }}\right| }_{I}^{n+1}-\boldsymbol{M}\, {\boldsymbol{\Delta q}}_{I}^{n+1}-\, {\boldsymbol{D}}_{\nu }\, \boldsymbol{\Delta }{\tilde{\boldsymbol{\nu }}}_{I}^{n+1}\right) \end{matrix}\quad \, \,</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(19d) |
|} | |} | ||
| Line 356: | Line 356: | ||
<math>\boldsymbol{M}=\int_{{\Omega }_{E}}^{}{\boldsymbol{\phi }}^{T}\boldsymbol{\phi }d\Omega ;\, \, {\boldsymbol{M}}_{\boldsymbol{D}}=</math><math>\left\{ \begin{matrix}\frac{{\Omega }_{E}}{8}\quad \quad \quad \, \, \hbox{for diagonal terms.}\, \\\\0\quad \quad \quad \, \hbox{off-diagonal terms.}\end{matrix}\right.</math> | <math>\boldsymbol{M}=\int_{{\Omega }_{E}}^{}{\boldsymbol{\phi }}^{T}\boldsymbol{\phi }d\Omega ;\, \, {\boldsymbol{M}}_{\boldsymbol{D}}=</math><math>\left\{ \begin{matrix}\frac{{\Omega }_{E}}{8}\quad \quad \quad \, \, \hbox{for diagonal terms.}\, \\\\0\quad \quad \quad \, \hbox{off-diagonal terms.}\end{matrix}\right.</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(20a) |
|} | |} | ||
| Line 373: | Line 373: | ||
<math>\begin{matrix}{\left. \displaystyle {\boldsymbol{q}}_{{\mathit{\boldsymbol{v}}}_{\mathit{\boldsymbol{i}}}}\right| }^{n}=-{cb}_{1}\, {{S}_{E}}^{n}\, {\left\{ \rho \tilde{\nu }{v}_{i}\right\} }^{n}+{cw}_{1}\, {{{f}_{w}}_{E}}^{n}\, {\left\{ \rho \, {\left( \frac{\tilde{\nu }}{d}\right) }^{2}{v}_{i}\right\} }^{n}-\frac{{cb}_{2}}{\sigma }\, {\left( {\rho }_{E}\, {\left. \frac{\partial \tilde{\nu }}{\partial {x}_{i}}\right| }_{E}{\left. \frac{\partial \tilde{\nu }}{\partial {x}_{i}}\right| }_{E}\right) }^{n}{{\boldsymbol{v}}_{i}}^{n}+\\ \displaystyle+\, \frac{1}{\sigma }{\left[ \left( {\nu }_{E}+{\tilde{\nu }}_{E}\right) {\left. \frac{\partial \rho }{\partial {x}_{i}}\right| }_{E}{\left. \frac{\partial \tilde{\nu }}{\partial {x}_{i}}\right| }_{E}\right] }^{n}{{\boldsymbol{v}}_{i}}^{n}\end{matrix}</math> | <math>\begin{matrix}{\left. \displaystyle {\boldsymbol{q}}_{{\mathit{\boldsymbol{v}}}_{\mathit{\boldsymbol{i}}}}\right| }^{n}=-{cb}_{1}\, {{S}_{E}}^{n}\, {\left\{ \rho \tilde{\nu }{v}_{i}\right\} }^{n}+{cw}_{1}\, {{{f}_{w}}_{E}}^{n}\, {\left\{ \rho \, {\left( \frac{\tilde{\nu }}{d}\right) }^{2}{v}_{i}\right\} }^{n}-\frac{{cb}_{2}}{\sigma }\, {\left( {\rho }_{E}\, {\left. \frac{\partial \tilde{\nu }}{\partial {x}_{i}}\right| }_{E}{\left. \frac{\partial \tilde{\nu }}{\partial {x}_{i}}\right| }_{E}\right) }^{n}{{\boldsymbol{v}}_{i}}^{n}+\\ \displaystyle+\, \frac{1}{\sigma }{\left[ \left( {\nu }_{E}+{\tilde{\nu }}_{E}\right) {\left. \frac{\partial \rho }{\partial {x}_{i}}\right| }_{E}{\left. \frac{\partial \tilde{\nu }}{\partial {x}_{i}}\right| }_{E}\right] }^{n}{{\boldsymbol{v}}_{i}}^{n}\end{matrix}</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(20b) |
|} | |} | ||
| Line 390: | Line 390: | ||
<math>\begin{matrix}\displaystyle{\boldsymbol{h}}^{n}={\displaystyle\int}_{{\Gamma }_{E}}^{}{{\boldsymbol{\phi }}^{\ast }}^{T}\, \left( \boldsymbol{\phi }{{\boldsymbol{\, v}}_{j}}^{n}\right) \left[ \eta \left( \frac{\partial \boldsymbol{\phi }}{\partial {x}_{i}}\, {{\boldsymbol{v}}_{j}}^{n}+\frac{\partial \boldsymbol{\phi }}{\partial {x}_{j}}{{\boldsymbol{v}}_{i}}^{n}\right) +\lambda \left( \frac{\partial \boldsymbol{\phi }}{\partial {x}_{k}}\, {{\boldsymbol{v}}_{k}}^{n}\right) {\delta }_{ij}\right] {n}_{i}\, \, d\Gamma +\\ \displaystyle +\, {\displaystyle\int}_{{\Gamma }_{E}}^{}{{\boldsymbol{\phi }}^{\ast }}^{T}\, \left( {K+K}_{t}\right) \, \left( \frac{\partial \boldsymbol{\phi }}{\partial {x}_{i}}{\boldsymbol{u}}^{n}\right) \, {n}_{i}\, \, d\Gamma \end{matrix}</math> | <math>\begin{matrix}\displaystyle{\boldsymbol{h}}^{n}={\displaystyle\int}_{{\Gamma }_{E}}^{}{{\boldsymbol{\phi }}^{\ast }}^{T}\, \left( \boldsymbol{\phi }{{\boldsymbol{\, v}}_{j}}^{n}\right) \left[ \eta \left( \frac{\partial \boldsymbol{\phi }}{\partial {x}_{i}}\, {{\boldsymbol{v}}_{j}}^{n}+\frac{\partial \boldsymbol{\phi }}{\partial {x}_{j}}{{\boldsymbol{v}}_{i}}^{n}\right) +\lambda \left( \frac{\partial \boldsymbol{\phi }}{\partial {x}_{k}}\, {{\boldsymbol{v}}_{k}}^{n}\right) {\delta }_{ij}\right] {n}_{i}\, \, d\Gamma +\\ \displaystyle +\, {\displaystyle\int}_{{\Gamma }_{E}}^{}{{\boldsymbol{\phi }}^{\ast }}^{T}\, \left( {K+K}_{t}\right) \, \left( \frac{\partial \boldsymbol{\phi }}{\partial {x}_{i}}{\boldsymbol{u}}^{n}\right) \, {n}_{i}\, \, d\Gamma \end{matrix}</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(20c) |
|} | |} | ||
| Line 423: | Line 423: | ||
<math>{r}_{I+1}^{\rho e}={\left[ \frac{\sum _{N}^{}{\left| {\rho e}_{I+1}-{\rho e}_{I}\right| }^{2}}{\sum _{N}^{}{{\rho e}_{I}}^{2}}\right] }^{\frac{1}{2}}\leq To</math> | <math>{r}_{I+1}^{\rho e}={\left[ \frac{\sum _{N}^{}{\left| {\rho e}_{I+1}-{\rho e}_{I}\right| }^{2}}{\sum _{N}^{}{{\rho e}_{I}}^{2}}\right] }^{\frac{1}{2}}\leq To</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(21) |
|} | |} | ||
| Line 435: | Line 435: | ||
| <math>{R}^{n+1}={\left[ \sum _{N}^{}{\left| {\rho }^{n+1}-{\rho }^{n}\right| }^{2}\right] }^{\frac{1}{2}}\leq TO{L}_{t}</math> | | <math>{R}^{n+1}={\left[ \sum _{N}^{}{\left| {\rho }^{n+1}-{\rho }^{n}\right| }^{2}\right] }^{\frac{1}{2}}\leq TO{L}_{t}</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(22) |
|} | |} | ||
| Line 451: | Line 451: | ||
| <math>{\, \Phi }_{N}=\frac{1}{8}\left( 1+{\xi }_{1N}{\xi }_{1}\right) \left( 1+{\xi }_{2N}{\xi }_{2}\right) \left( 1+\right. </math><math>\left. {\xi }_{3N}{\xi }_{3}\right)</math> | | <math>{\, \Phi }_{N}=\frac{1}{8}\left( 1+{\xi }_{1N}{\xi }_{1}\right) \left( 1+{\xi }_{2N}{\xi }_{2}\right) \left( 1+\right. </math><math>\left. {\xi }_{3N}{\xi }_{3}\right)</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(23) |
|} | |} | ||
| Line 464: | Line 464: | ||
| <math>\frac{\partial \boldsymbol{\phi }}{\partial {x}_{j}}={I}_{ij}^{-1}\frac{\partial \boldsymbol{\phi }}{\partial {\xi }_{i}}</math> | | <math>\frac{\partial \boldsymbol{\phi }}{\partial {x}_{j}}={I}_{ij}^{-1}\frac{\partial \boldsymbol{\phi }}{\partial {\xi }_{i}}</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(24) |
|} | |} | ||
| Line 479: | Line 479: | ||
\displaystyle\frac{\partial \boldsymbol{\phi }}{\partial {\xi }_{3}}{\boldsymbol{x}}_{1}\, \, \displaystyle\frac{\partial \boldsymbol{\phi }}{\partial {\xi }_{3}}{\boldsymbol{x}}_{2}\, \, \displaystyle\frac{\partial \boldsymbol{\phi }}{\partial {\xi }_{3}}{\boldsymbol{x}}_{3}\end{matrix}\right]</math> | \displaystyle\frac{\partial \boldsymbol{\phi }}{\partial {\xi }_{3}}{\boldsymbol{x}}_{1}\, \, \displaystyle\frac{\partial \boldsymbol{\phi }}{\partial {\xi }_{3}}{\boldsymbol{x}}_{2}\, \, \displaystyle\frac{\partial \boldsymbol{\phi }}{\partial {\xi }_{3}}{\boldsymbol{x}}_{3}\end{matrix}\right]</math> | ||
|} | |} | ||
| − | | style="text-align: right | + | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(25) |
|} | |} | ||
Revision as of 14:20, 25 September 2019
Abstract
As already demonstrated by different authors, the Taylor-Galerkin (TG) scheme, in the context of the Finite Element Method (FEM), is particularly suitable for the solution of supersonic flows. The TG scheme, using hexahedral finite elements with analytical evaluation of element matrices, is applied in this work. Tools to avoid locking and a shock capturing technique for the solution of supersonic viscous and non-viscous compressible flows are also employed. However, TG scheme usually presents instabilities in subsonic flows. Even in cases in which the free stream Mach number corresponds to supersonic flows, there are always flow regions, specifically near the walls of the immersed obstacles, where the speed is lower, and the local Mach number corresponds to a subsonic flow. The Characteristic Based Split (CBS) scheme was developed in order to obtain a single method to improve the behavior with respect to TG method in subsonic and supersonic regimes. In the last two decades some works have shown advantages in convergence rates of the CBS method when compared to the TG algorithm. However, simulation time increases in the CBS method since split operations, typical of this algorithm, imply in additional element loops. In this paper a hybrid algorithm called Modified-Taylor-Galerkin scheme (MTG) is proposed. This algorithm presents advantages with respect to TG and CBS schemes in terms of convergence properties and computational processing time. In order to get an efficient algorithm, the element matrices are analytically integrated. This is performed with two different approaches. In the first approach the inverse matrix and the determinant of the Jacobian matrix at element level are evaluated with a reduced integration form, using the point located in the center of the element for mass, convective, diffusive and stabilization element matrices; all these matrices are integrated analytically. In the second approach, mass and convective matrices are calculated by a complete integration scheme (including the inverse matrix and the determinant of the Jacobian matrix at element level in the analytical expression to be integrated) and the diffusive and stabilization matrices are calculated with a reduced integration form, using the point located at the center of the element to calculate the inverse matrix and the determinant of the Jacobian matrix at element level. Finally, this work incorporates the Spalart-Allmaras (S-A) turbulence model using a conservative version of the transport equation, as proposed by the authors of the original S-A model in a later paper. Algorithms are tested to determine convergence rate improvements in both laminar and turbulent cases and for different Mach numbers (supersonic, transonic and subsonic flows).
Keywords: Compressible turbulent, FEM, Taylor-Galerkin, CBS, Spalart-Allmaras.
1. Introduction
Different authors have shown that, in the context of FEM, the TG scheme is particularly suitable for the solution of compressible flows involving supersonic shock waves at high Mach numbers [1-2]. Burbridge and Awruch [3] applied this scheme using hexahedral finite elements with analytical evaluation of the element matrices but evaluating the inverse matrix and the determinant of the Jacobian matrix at the center of the element, together with a shock capturing technique for the solution of viscous and non-viscous compressible laminar flows [4]. However, for subsonic regimes (with M <0.9) or for incompressible flows, the TG scheme presents certain instabilities, as reported by Zienkiewicz and Taylor [5].
Even in supersonic flows there are always some regions of subsonic flow near the walls of immersed obstacles. It is therefore necessary to improve the convergence of the algorithms for low Mach numbers, even when they are applied to flows where Mach number at the freestream is supersonic.
To address this problem the CBS scheme was developed by some researchers, such as Zienkiewicz and Codina [6], to obtain a single algorithm with a good behavior, both in subsonic and supersonic regimes. During the last two and a half decades this algorithm was established as an acknowledged tool for the computation of a wide spectrum of incompressible or compressible flows at different Mach numbers [7-10]. The work presented by Zienkiewicz et al. in Reference [11] shows the explicit form of the algorithm and its performance for subsonic, transonic and supersonic flows.
However, in spite of convergence rates improvements, it is observed that computational processing time increases due to split operations involved in the CBS scheme. For this reason, an alternative hybrid algorithm is presented here, looking for improvements in convergence properties with respect to the TG algorithm without an increase of the computational processing time.
It is important to mention that, as in previous works [3,12], an analytical integration of the element matrices was used here, evaluating the inverse matrices and the determinants of the Jacobian matrices at the center of the element. But in this work, an alternative approach was also applied, integrating mass and convective matrices in an exact analytical form (including the inverse matrices and the determinants of the Jacobian matrices in the analytical expressions) using Maxima symbolic resolution program (http://maxima.sourceforge.net/). Stabilization and convective matrices were integrated in reduced form (evaluating the inverse matrices and the determinants of the Jacobian matrices at the center of the elements), as performed in [3,12].
Finally, the turbulence model of Spalart-Allmaras (S-A) was implemented here in order to test the behavior of these schemes, and their convergence properties, in the context of the FEM for turbulent flows. The form of the S-A model implemented here was proposed for the first time by Allmaras et al. [13] in 2012, using the transport equation for the turbulent modified kinematic viscosity variable Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \tilde{\nu }}
in conservative form.
TG, CBS and MTG schemes are presented here with some cases of laminar and turbulent flows, comparing their performances.
2. Governing equations with Spalart-Allmaras turbulence model
In an Eulerian description, the dimensionless system of partial differential equations governing the fluid dynamics problem in conservative form, and filtered with a Favre filter (FANS), may be written as follows:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \frac{\partial \boldsymbol{U}}{\partial t}+\frac{\partial {\boldsymbol{F}}_{i}}{\partial {x}_{i}}+
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \frac{\partial {\boldsymbol{G}}_{i}}{\partial {x}_{i}}+\boldsymbol{Q}=\mathit{\boldsymbol{0}} |
(1) |
where
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \boldsymbol{U}=\left\{ \begin{matrix}\rho \\ \rho {v}_{1} \\ \rho {v}_{2} \\ \rho {v}_{3} \\ \rho e \\ \rho \tilde{\nu }\end{matrix}\right\} ;{\boldsymbol{F}}_{i=}\left\{ \begin{matrix}\rho {v}_{i} \\ \rho {v}_{1}{v}_{i}+p{\delta }_{i1} \\ \rho {v}_{2}{v}_{i}+p{\delta }_{i2} \\ \rho {v}_{3}{v}_{i}+p{\delta }_{i3} \\ {v}_{i}\left( \rho e+p\right) \\ \rho \tilde{\nu }{v}_{i}\end{matrix}\right\} ;{\boldsymbol{G}}_{i}=
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left\{ \begin{matrix} 0 \\ -{\tau }_{i1} \\ -{\tau }_{i2} \\ -{\tau }_{i3} \\ -{\tau }_{ij}{v}_{j}-{q}_{i}-{\psi }_{i} \\ {\varphi }_{i}\end{matrix}\right\} ;\boldsymbol{Q}= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left\{ \begin{matrix} 0 \\ 0 \\ 0 \\ 0 \\ 0 \\P+D+{D}_{T}\end{matrix}\right\} |
(2) |
The dimensionless form has been obtained using a reference length, the freestream density and the freestream speed of sound. In Eq. (2):
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {v}_{i}\,}
are velocity vector components.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \rho}
is the specific mass.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle p}
is the thermodynamic pressure.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\tau }_{ij}\,}
are the viscous components of the stress tensor.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle e}
is the specific total energy.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \tilde{\nu }}
is the transport equation variable of the S-A conservative turbulence model.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {q}_{i}}
are the conduction heat flux vector components.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\delta }_{ij}\,}
is the Kronecker delta.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\varphi }_{i}}
is the diffusive term of the S-A transport equation.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\psi }_{i}}
is the diffusive term associated with turbulent kinetic energy ( Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k}
) in the energy equation.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle P}
is the production term and D is the destruction term of the S-A transport equation.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {D}_{T}\,}
contains additional terms of the S-A transport equation.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {x}_{i}\,}
and t are the spatial and temporal coordinates, respectively.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle i,j}
=1,2,3.
The viscous stress components of a Newtonian fluid are given by:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\tau }_{ij}=\left( \mu +{\mu }_{t}\right) \left[ \frac{\partial {v}_{i}}{\partial {x}_{j}}+\right.
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left. \frac{\partial {v}_{j}}{\partial {x}_{i}}\right] -\frac{2}{3}\left( \mu +\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left. {\mu }_{t}\right) \frac{\partial {v}_{k}}{\partial {x}_{k}}{\delta }_{ij}- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \frac{2}{3}\rho k{\delta }_{ij} |
(3) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mu}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\mu }_{t}}
are the molecular viscosity and the eddy viscosity, respectively, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k}
is the turbulent kinetic energy and the subscript Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k}
runs from 1 to 3.
The components of the conduction heat flux vector are:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {q}_{i}=\left( K+{K}_{t}\right) \frac{\partial u}{\partial {x}_{i}} | (4) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle K}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {K}_{t}}
are the thermal conductivity and the turbulent thermal conductivity respectively and u is the specific internal energy.
The terms Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\psi }_{i}}
in the energy equation are given by the following expression:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\psi }_{i}=\left( \mu +\frac{{\mu }_{t}}{{\sigma }_{k}}\right) \frac{\partial k}{\partial {x}_{i}} | (5) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\sigma }_{k}}
is a constant which depends on the turbulence model. Since the S-A turbulence model does not include the calculation of turbulent kinetic energy, in this work, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\psi }_{i}}
in Eq. (2) and the third term in Eq. (3) are neglected.
The terms Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\varphi }_{i}}
in the S-A transport equation are given by the following expression:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\varphi }_{i}=\frac{1}{\sigma }\left( \mu +\rho \tilde{\nu }\right) \frac{\partial \tilde{\nu }}{\partial {x}_{i}} | (6) |
The transport equation also includes the additional term Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {D}_{T}}
, expressed by
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {D}_{T}=-\frac{1}{\sigma }\, \rho \, {cb}_{2}\frac{\partial \tilde{\nu }}{\partial {x}_{i}}\frac{\partial \tilde{\nu }}{\partial {x}_{i}}+
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \, \frac{1}{\sigma }(\nu +\tilde{\nu })\, \frac{\partial \rho }{\partial {x}_{i}}\frac{\partial \tilde{\nu }}{\partial {x}_{i}} |
(7) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \sigma}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {cb}_{2}}
are constants of the S-A turbulence model.
Finally, the production and destruction terms, in the conservative transport equation of S-A are given by
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): P=-{cb}_{1}\, S\, \rho \tilde{\nu };\, \, D={cw}_{1}\, {f}_{w}\, \rho \, {\left( \frac{\tilde{\nu }}{d}\right) }^{2} | (8) |
where the first term Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle P} represents the production and the second term D the destruction, being Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {cb}_{1}} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle S} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {cw}_{1}} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {f}_{w}}
constants and variables of the S-A turbulence model and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle d}
is the smallest distance from a field point (or node) to the nearest solid boundary (wall). The transport equation used here is expressed in its conservative form (Allmaras et al. [13]). Here the trip-term and the term of laminar suppression, both included in the original expressions of the model, are not considered.
To close the system of equations it is still necessary to define all parameters associated with the S-A turbulence model; they are given by
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \sigma =\frac{2}{3};{\, cb}_{1}=0.1355;\, {\, cb}_{2}=0.622
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \boldsymbol{\Omega }=\, \, \frac{1}{2}\, rot\mathit{\boldsymbol{V;\, }}\Omega =\sqrt{\boldsymbol{\Omega }\cdot \boldsymbol{\Omega }}
|
(9) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Omega}
is the magnitude of the angular velocity vector, and
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {f}_{v2}=1-\, \, \frac{\chi }{1+\chi \, {f}_{v1}}\, ;\, \chi =\frac{\rho \tilde{\nu }}{\mu }\, ;\, {f}_{v1}=
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \frac{{\chi }^{3}}{{\chi }^{3}+\, {{c}_{v1}}^{3}};\, {c}_{v1}=7.1
|
(10) |
The turbulent viscosity and turbulent thermal conductivity are defined as follows:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\mu }_{t}=\rho \, \tilde{\nu }\, {f}_{v1};\, {K}_{t}={\mu }_{t}\frac{\gamma }{{Pr}_{t}};\, \, \gamma =
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \frac{{c}_{p}}{{c}_{v}}=1.4 |
(11) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {Pr}_{t}}
is the turbulent Prandtl number, which is adopted as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {Pr}_{t}=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 0.9
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {c}_{v}\,}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {c}_{p}\,}
are the specific heat coefficients at constant volume and at constant pressure, respectively.
The state equation for a perfect gas (air) is given by:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): p=\left( \gamma -1\right) \rho u; | (12) |
The internal energy u and the temperature are related to the independent field variables by the following expression:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): u={c}_{v}\, T=e-\frac{1}{2}{v}_{i}{v}_{i}-k | (13) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle T}
is the temperature, being Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k}
neglected in the S-A model.
The Sutherland law is used in this work to establish the dependence of viscosity and thermal conductivity with respect to temperature (dimensionless internal energy). This law can be expressed as follows:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \mu ={\mu }_{\infty }\frac{S+{u}_{\infty }}{S+u}{\left( \frac{u}{{u}_{\infty }}\right) }^{\frac{3}{2}};\quad \quad K=
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {K}_{\infty }\, \frac{{S}_{k}+{u}_{\infty }}{{S}_{k}+u}{\left( \frac{u}{{u}_{\infty }}\right) }^{\frac{3}{2}} |
(14) |
Values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle S}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {S}_{k}}
may be found in reference [14] for the Sutherland law based on temperature; they must be multiplied by Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {c}_{v}}
and divided by the square of the freestream sound speed to obtain the dimensionless form, which was used in this work.
Initial and boundary conditions must be applied to Eq. (1) in order to obtain a mathematical model correctly defined. The forced boundary conditions (or Dirichlet boundary conditions) are given by:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \begin{matrix}\begin{matrix}{v}_{i}=\overline{{v}_{i}}\quad on\quad {\Gamma }_{v}\\\rho =\overline{\rho }\quad \, \, on\quad {\Gamma }_{\rho }\\u=\overline{u}\quad \, \, on\quad {\Gamma }_{u}\end{matrix}\\\tilde{\nu }=\overline{\tilde{\nu }}\quad \quad on\quad {\Gamma }_{\tilde{\nu }}\end{matrix} | (15) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \overline{{v}_{i}}}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \overline{\rho }\, ,\, \overline{u}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \overline{\tilde{\nu }}}
are prescribed values of the velocity vector components, the density (specific mass), the specific internal energy and the variable Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \tilde{\nu }}
, at the boundary surfaces Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Gamma }_{v}}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Gamma }_{\rho }}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Gamma }_{u}}
and Г Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \tilde{v}\quad}
respectively.
The natural boundary conditions (or Neumann boundary conditions) are defined by the following expressions:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {{\sigma }_{ij}\, n}_{j}=\hat{{t}_{i}}\quad on\quad {\Gamma }_{\sigma }
|
(16) |
In Eq. (16)
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {n}_{i}\,}
are the cosines of the angles formed by the normal vector to the surface Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Gamma }_{\sigma }}
or Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Gamma }_{q}}
and the global reference axes Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {x}_{i}}
.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat{{t}_{i}}\,}
are the tractions acting on Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Gamma }_{\sigma }}
.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\sigma }_{ij}\, =\, -p{\delta }_{ij}+\, {\tau }_{ij}}
are the components of the Cauchy stress tensor.
- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat{q}\,}
is the heat flow acting on the perpendicular direction to Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, {\Gamma }_{q}}
.
- Effects of thermal radiation on boundary surfaces have not been considered.
For non-viscous flows the transport equation of the variable Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \tilde{\nu }}
and the diffusive terms Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol{G}}_{i}}
in Eq. (1) should be omitted. In this case, only normal velocity components to the solid boundaries are prescribed with null values and only the normal components of the surface tractions on the outlet boundaries are considered.
3. Finite element formulation
In this section the well-known TG, CBS and Single Step CBS algorithms are mentioned but the respective equations are not shown for the sake of brevity. The reader can look at the references [1,3,5,6,8] to get more details. Instead, the full matrix equations for the proposed MTG scheme are presented in the section 4.
3.1 Taylor-Galerkin (TG) scheme
To discretize Eq. (1) in time, the unknown vector Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{U}}
is expanded in Taylor series omitting higher order terms. Then the space discretization is done using the classic Bubnov-Galerkin technique in the context of the finite element method. The Green-Gauss theorem is applied to terms containing second order derivatives (terms with third or higher order derivatives are neglected) to weaken the demands with respect to continuity of the element interpolation functions and their derivatives, obtaining the matrix equations [1-3].
A conditionally stable explicit time integration version of the TG scheme [4,5,15], may be obtained working with a lumped mass matrix (where only elements of the main diagonal are different from zero).
Also, it is important to point out that the increments of the flux variables obtained for the current time step Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, n+} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 1\,
with the TG scheme depend of flow variables at the previous time step Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle n\,}
and at the current time step Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, n+}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 1 . So, an iterative scheme to solve terms involving flow variables at the new time step is required [1,3].
3.2 Characteristic Based Split (CBS) method
The CBS method consists of four steps [6]. In the first step, intermediate values of the conservative flow variables of the momentum equations are calculated, omitting the terms containing pressure gradients. In the second step, the continuity equation is solved to obtain the density increment. In the third step, the flow variables of the momentum equations are updated. In the fourth step, the energy equation and the S-A transport equation in its conservative form are solved. Finally, pressure is calculated using the state equation.
Once the temporal discretization is performed, spatial domain discretization is carried out using the classic Bubnov-Galerkin method. The Green-Gauss theorem is applied to terms containing second order derivatives to weaken the demands with respect to continuity of the element interpolation functions and their derivatives.
This method may be applied with implicit, semi-implicit and explicit time integration schemes [11]. An explicit scheme was implemented in this work in order to compare its behavior with respect to other explicit algorithms presented here.
Matrix expressions obtained with the CBS method may be found in references [5,6,8].
3.3 Single-Step CBS method
An alternative of the CBS method consists in the elimination of the split in the momentum equation in order to obtain a one-step algorithm. Some authors [5,6] suggest that, although the single-step CBS scheme saves computational processing time with respect to the classical CBS method, it is not suitable to simulate flows with high Mach numbers. The explicit version of the single-step CBS was implemented in this work. The common matrix expressions obtained with single-step CBS scheme can be found in reference [5].
4. The modified Taylor-Galerkin scheme
It may be observed that stabilization problems arise when TG scheme is used, and they are linked with density and pressure values, being these variables basically controlled by the continuity equation. The TG method is an efficient technique in terms of computational processing time, but schemes such as CBS have much better convergence properties for problems involving a wide range of Mach numbers. In order to link both features (computational efficiency and good convergence properties), a modification of the TG scheme is proposed. The new scheme is called Modified Taylor-Galerkin (MTG) scheme, which consists to change the stabilization term on the continuity equation by a similar expression, but defined in terms of pressure, as used in the CBS method. Since the TG scheme has good convergence rates for supersonic flows, the modification is only implemented in regions where the local Mach number is less than one. For this reason, a blending function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle f}
is proposed such that, the continuity equation discretized in time results:
|
(17) |
where
|
(18) |
The blending function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle f}
depends of the Mach number, which is evaluated locally. Momentum and energy equations as well as their iterative terms are not modified, remaining as in the original TG scheme.
Spatial discretization is performed using the classical Bubnov-Galerkin method. The Green-Gauss theorem is applied to terms containing second order derivatives to weaken the demands with respect to continuity of the element interpolation functions and their derivatives.
The following matrix expressions are obtained for the element vectors of independent unknowns’ increments:
- Mass conservation equation:
|
(19a) |
- Momentum conservation equations:
|
(19b) |
- Energy conservation equation:
|
(19c) |
- Transport equation of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \tilde{\nu }}
(being Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \tilde{\nu }}
a variable used in the S-A model):
|
(19d) |
An iteration procedure is required for those right-hand terms evaluated on the current time step Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle n+}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 1
, where the subscript Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I}
indicates the previous iteration and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I+}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 1
corresponds to the current one. Element matrices are indicated below (with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle i,j,k=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 1,2,3 ):
|
(20a) |
The element vectors are:
|
(20b) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, {\boldsymbol{v}}_{i}=\left\{ {v}_{i}\right\}}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{p}=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left\{ p\right\} \boldsymbol{,\, }\boldsymbol{u}=\left\{ u\right\}
are Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \,}
the velocity components, pressure and internal energy element vectors, containing their corresponding eight nodal values. The quantities between braces are element vectors built with the nodal values of the variables inside the braces. The boundary vectors are:
|
(20c) |
In these expressions Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Omega }_{E}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Gamma }_{E}}
are the element volume and its boundary surface respectively; expressions with the subscript Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle E}
are element averages; Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{\phi }}
= [Ф1, Ф2, … Ф8] is the row matrix containing the interpolation functions for each local node and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol{\phi }}^{\ast }}
is the row matrix containing the interpolation functions evaluated at the boundary surface, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{M}}
is the consistent mass matrix, while Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol{M}}_{\boldsymbol{D}}}
is the lumped mass matrix and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \left\{ 1\right\}}
is a vector where all its components are equal to one. Terms that involve turbulent kinetic energy have been neglected.
The variables with superscript Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle n}
are evaluated on previous time step and the variables with superscript Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle n+}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 1
are evaluated at current time step. Boundary vectors appearing in the S-A transport equation and in the iterative terms of all the equations were neglected.
Once these equations are assembled for all the elements of the mesh and the corresponding boundary conditions are applied, the nodal values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \rho}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \rho {v}_{j}}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \rho e} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \rho \tilde{\nu }}
can be calculated at each time step using an iterative scheme. Thus, the nodal values of the components of velocity, the specific total energy and the variable Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \tilde{\nu }}
can be obtained immediately. Then the nodal values of the specific internal energy are calculated and, using the state equation, the nodal pressure values are also obtained.
When viscous terms are not considered, terms containing matrices Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol{D}}_{ij}} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol{E}}_{i}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{K}}
as well as vectors Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {{\boldsymbol{t}}_{j}}^{n}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol{h}}^{n}}
and the transport equation corresponding to the turbulence model are omitted.
Stabilization matrices for TG, CBS and MTG schemes are evaluated using a local (internal) time step Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Delta t}_{in}} . For the TG scheme there are no appreciable differences in using a local internal time step or a single time step Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Delta t}_{ext}}
[16].
5. Stability condition, shock capture and convergence
Considering that these algorithms were implemented in their explicit form, they are conditionally stable and the stability condition of Courant-Friedrichs-Lewy (CFL) for each element, in dimensionless form, is used to ensure stability. A safety coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \beta}
is adopted with values Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 0.2\leq \beta \leq \, 0.5}
[3].
To capture the strong discontinuities and eliminate high frequency oscillations near shock waves, a well-known method is used to add artificial viscosity selectively in regions where high-pressure gradients are observed [17]. An artificial damping coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {C}_{AD}}
is used to control de amount of artificial viscosity.
The convergence of the iterative process is obtained when the following conditions are simultaneously satisfied:
|
(21) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N}
is an index of nodes ( Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 1,2,3,\ldots ,NNODES ), Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle NNODES}
is the total number of nodes of the finite element mesh and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle To}
is the tolerance (in this paper Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle To=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {10}^{-3}
was adopted). The solution process ends when the variable Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle t}
(time) reaches a value previously established by the user, or when the steady state is reached. This last condition is considered fulfilled when the following expression is satisfied:
|
(22) |
In this work Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle TO{L}_{t}={10}^{-5}} was adopted.
6. Explicit integration of element matrices
Eight nodes isoparametric hexahedral elements were used in this work. Figures 1 and 2 show the element in the physical and computational spaces, respectively. The interpolation functions of this element are given by the following expressions:
|
(23) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N=1,2,3,\ldots ,8}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\xi }_{jN}}
are the natural coordinates of the node Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N}
. The derivatives of the interpolation functions with respect to the global coordinates are given by:
|
(24) |
being Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {I}_{ij}^{-1}}
the components of the inverse of the Jacobian matrix Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathit{\boldsymbol{\, }}\mathit{\boldsymbol{I}}}
, which can be expressed as follows:
|
(25) |
|
| Figure 1. Physical space |
|
| Figure 2. Computational space |
Here Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol{x}}_{i}}
is the vector containing global coordinates of the eight nodes of each element, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, \left| \mathit{\boldsymbol{I}}\right|}
is the determinant of the Jacobian matrix and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle d\Omega =}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left| \mathit{\boldsymbol{I}}\right| \, d{\xi }_{1}\, d{\xi }_{2\, }\, d{\xi }_{3\, }
is the differential of the element volume.
By substituting Eqs. (23), (24) and (25) on the element matrices and vectors defined by the TG, CBS and MTG schemes, integrals in the computational domain are obtained. These integrals are commonly calculated by means of a numerical integration scheme to obtain components of the element matrices that will be later assembled to get the global system of equations. However, in this work, matrices were calculated by means of analytical integration using the Maxima symbolic resolution software. This implies an improvement in computational processing time efficiency.
In a first approach the inverse matrix elements and determinants of the Jacobian matrix were evaluated at the center of the element, where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\xi }_{1}=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \, {\xi }_{2\, }=\, {\xi }_{3\, }=0,
and these values were used to integrate matrices corresponding to variables time derivatives, convective, diffusive and stabilization terms [3,12]. Special attention must be given to elements distortion and an hourglass control technique to avoid spurious modes on diffusive terms are also necessary [12].
Alternatively, in a second approach, the mass matrix Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{\, }\boldsymbol{M}}
(representing variables time derivatives) and matrix Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol{B}}_{i}}
(representing convective terms) were integrated analytically, including the inverse matrix elements and determinant of the Jacobian matrix in the analytical expressions (without evaluating these terms at the center of the element). However, to integrate matrices corresponding to stabilization and diffusive terms, the procedure used in the first approach was followed (inverse matrix elements and determinants of the Jacobian matrix were evaluated at the center of the element). Both integration approaches were used for the implementation of TG, CBS and MTG schemes.
7. Numerical simulations and comparative results
7.1 Transonic flow around a NACA 0012 airfoil
The laminar transonic flow around a NACA 0012 airfoil with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 0.85\,
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle Re=2000\,}
is analyzed. Results obtained with version 17.1 of ANSYS/Fluent are compared with those obtained with the MTG scheme for the same problem. All simulations performed for this case were obtained with TG, MTG and CBS schemes with the second integration approach, that is, a complete analytical integration (including inverse matrix and determinant of the Jacobian matrix in the expressions to be evaluated analytically) for the mass and convective matrices, while the inverse matrix and determinant of the Jacobian matrix were evaluated at the center of the elements for the stabilization and diffusive matrices.
This problem is described with dimensionless quantities. The mesh used to solve this problem with 386 560 hexahedral elements and 775 684 nodes is shown in Figure 3. The coordinates axis are Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {x}_{1}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y={x}_{2}}
. Figures 4 and 5 show details of the finite element mesh. It is considered as being three-dimensional mesh with a single element in the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle z=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {x}_{3}
direction. The same mesh is used to run this problem with ANSYS/Fluent.
Considering that the mesh is three-dimensional with a single element in the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {x}_{3}}
direction, the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, {v}_{3}=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 0.0
condition is imposed on the lateral boundaries to avoid mass flow through those boundaries, ensuring the two-dimensional characteristic of the flow.
|
| Figure 3. Side view of whole domain and the finite element mesh |
|
| Figure 4. Side view of the mesh |
|
| Figure 5. Details of the airfoil leading edge |
The dimensionless inlet boundary conditions are:
|
On the solid boundaries the non-slip condition is applied. Then, on the airfoil wall:
|
together with the stagnation temperature, which is specified using the specific internal energy, according to the following expression:
|
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle r=0.89}
is the recovery factor (an empirical factor introduced because in practice the energy recovery is not perfect [14]). The conditions at the output boundaries are specified using the vectors Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol{t}}_{j}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{q}}
.
Initial conditions are given by:
|
These initial conditions are defined for all nodes of the finite element mesh, excepting those where Dirichlet boundary conditions are prescribed.
Using the safety coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \beta =0.2} , the effective dimensionless time step is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Delta t=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {0.5}^{\ast }{10}^{-4} . Artificial viscosity is not applied for this case Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle ({C}_{AD}=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 0.0) .
Figure 6 shows pressure contours obtained with the MTG scheme. Pressure contours obtained with CBS scheme are practically coincident with results obtained with MTG scheme and are not shown here. Figure 7 shows pressure contours obtained with ANSYS/Fluent 17.1 available at the GFC (Grupo de Fluidodinámica Computacional) of the UNLP (Universidad Nacional de La Plata). An excellent agreement was obtained between MTG and ANSYS/Fluent simulations. Figure 8 shows details of pressure contours at the leading edge of the airfoil obtained with the MTG scheme (CBS shows no appreciable difference). Figure 9 shows the same image, but obtained with the TG scheme, where pressure oscillations in zones of low air speeds can be observed. CBS and MTG schemes do not present these oscillations.
|
| Figure 6. MTG pressure contours |
|
| Figure 7. ANSYS / Fluent pressure contours |
|
| Figure 8. MTG pressure contours |
|
| Figure 9. TG pressure contours |
Figure 10 shows an image of the skin friction coefficient comparing the MTG scheme with values obtained with ANSYS/Fluent 17.1 showing a good agreement. Figure 11 shows the comparison of pressure coefficient obtained with ANSYS/Fluent 17.1 and the MTG scheme on the NACA 0012 airfoil extrados.
Here, the skin friction coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {c}_{f}}
and the pressure coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {c}_{p}}
were calculated using the following expressions:
|
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \tau}
is the wall tangent stress and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle p}
is the pressure, both acting on each point of the solid boundary surface.
Results obtained for these coefficients with MTG scheme applying an internal time step are very close to values obtained with ANSYS/Fluent 17.1.
|
| Figure 10. Skin friction coefficient over the airfoil extrados surface |
|
| Figure 11. Pressure coefficient over the airfoil extrados surface |
Figure 12 shows the convergence of the different schemes. The CBS method presents a clear improvement in density and pressure convergence when used with a local interior time step [16]. The use of an internal time step with the original TG method does not show comparative advantage with respect to the same method but taking a single time step. Notably, the performance of the MTG scheme with a local internal time step, is equivalent to that obtained with the CBS algorithm with a local internal time step.
|
| Figure 12. Residues showing convergence of the different algorithms |
Additional simulations were carried out with the single-step CBS method. Results are not shown here, but they do not show significant differences compared to those obtained with the CBS or MTG schemes.
Relative times for the different simulations are indicated in Table 1, assuming a reference time with a value equal to one for the fastest scheme (the single-step CBS).
| Algorithm | Relative wall clock time |
| Single-step CBS | 1.0 |
| TG | 1.01 |
| MTG | 1.04 |
| CBS | 1.37 |
It should be noted that the cases shown in Table 1 correspond to simulations of the same problem, with the same mesh and boundary conditions, on the same computer (Dell Precision 7610, 16 Cores, 64 Gb ram), and using the same number of threads (16) in parallel. The times that are compared are wall clock times between two recordings of results with the same number of time steps. This problem was also solved using the first integration approach, that is, analytical integration with evaluation of the inverse matrices and determinants of the Jacobian matrices at the center of the elements for all the mass, convective, stabilization and diffusive matrices. There are not significant differences between both integration approaches.
It should also be pointed out that the codes for the different schemes were developed starting from the same original TG MPI parallel code, including only those modifications required by CBS, single-step CBS and MTG schemes. This ensures that the code structure is basically the same, with smallest differences required by the different algorithms.
7.2 Turbulent supersonic flow on flat plate
The turbulent compressible flow on a flat plate in supersonic regime with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 2.0\,
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle Re=15\ast {10}^{6}\,}
is presented.
This problem is described at the NASA Langley Research Center Turbulence Modeling Resource on the web (https://turbmodels.larc.nasa.gov/), where results are shown for several turbulence models and they are compared with the semi-empirical theoretical correlations of Van Driest.
The geometry of this problem is shown in Figure 13. Considering that the mesh is three-dimensional with a single element in the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {x}_{2}}
direction, the condition Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {v}_{2}=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 0.0
is imposed on the lateral boundaries to avoid the mass flow through them, keeping the two-dimensional characteristic of the flow.
|
| Figure 13. 3D view and geometric characteristics |
Figure 14 shows the side view of the finite element mesh, which has 13 056 hexahedral elements and 26 578 nodes. The mesh is the same found in the NASA web repository for turbulence models (3-D with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 2\times 137 \times 97
nodes and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 2\times 113
points on the flat plate).
 
|
| Figure 14. Side view of the mesh |
The mesh was refined near the wall such that the estimate of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {y}^{+}}
is close to 0.3.
The boundary conditions at the inlet are:
|
On the solid boundaries (over the plate) the non-slip condition is applied:
|
together with the stagnation temperature, which is specified using the specific internal energy, according to the following expression:
|
with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle r=0.89} . The conditions in the output boundary are specified using the vectors Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol{t}}_{j}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{q}}
.
The initial conditions are given by:
|
These initial values are defined in all nodes of the mesh, except for the nodes where Dirichlet boundary conditions were prescribed.
Using the safety factor Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \beta =0.2} , the dimensionless time step is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Delta t=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {1.0}^{\ast }{10}^{-7} . The artificial viscosity coefficient used here is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {C}_{AD}=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 0.3 .
The pressure contours obtained with the TG, MTG and CBS schemes shown no appreciable differences. However, the researchers who developed the CBS scheme [6] point out that the single-step CBS scheme is not suitable for high Mach numbers. The simulations carried out in this work with this scheme for high Mach numbers indicate that the convergence of the single-step CBS presents some instabilities. Therefore, observations given in the reference [6] is confirmed here: single-step CBS scheme is not recommended for supersonic flows.
In Figure 15 an image of the skin friction coefficient, comparing with values provided by the semi-empirical correlation of Van Driest II (see the NASA Langley Research Center Turbulence Modeling Resource) is shown. Good agreement is observed. There are no appreciable differences between results obtained with different schemes (TG, CBS and MTG). All these schemes seem to be adequate in terms of quality of the results.
|
| Figure 15. Skin friction coefficient over the plate |
Finally, Figure 16 shows the dimensionless velocity profile, corresponding to the position of the plate where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {Re}_{\theta }=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 10\, 000
. As can be observed, results obtained in this work are very close to values shown at the NASA repository for the velocity profile, corresponding to the semi-empirical correlation of Van Driest I.
|
| Figure 16. Velocity profile at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {Re}_{\theta }=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 10\, 000 |
The convergence between the TG, MTG and CBS schemes does not differ significantly when the residuals are plotted (although they are not shown here); however, there are differences in the relative computational times, as shown in Table 2.
| Algorithm | Relative wall clock time |
| Single-step CBS | Not recommended |
| TG | 1.0 |
| MTG | 1.03 |
| CBS | 1.30 |
In this way, it is observed that, for supersonic Mach numbers, the three algorithms are suitable from the point of view of convergence rates, but there are appreciable differences in terms of wall clock computational time, being the CBS scheme the most expensive.
7.3 Axisymmetric Shock Wave Boundary Layer Interaction (ASWBLI). Turbulent cold hypersonic flow
The following case corresponds to the turbulent flow over an ogive-shaped obstacle in a cold hypersonic regime: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 7.11 , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \frac{Re}{L}=57\, 060\, \frac{1}{cm}} . This problem is described at the NASA Langley Research Center Turbulence Modeling Resource where results for several turbulence models are compared with experimental data. The purpose is to provide a validation case for turbulence models. Unlike verification, which seeks to establish whether a numerical model has been implemented correctly by comparison with other software or code, the validation of CFD results implies to compare with those obtained by experimental data, establishing the ability of a numerical model to reproduce the physics of the problem.
The experimental study conducted by Kussoy and Horstman [18] involved a cone cylinder with a taper angle of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 20}
degrees. The leading edge is a 10-degree cone that makes the transition through a circular arc to a constant-area cylinder with a radius of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 100\, mm}
. The cone cylinder with a taper angle of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 20}
degrees, designed to produce an oblique shock and an interaction zone between the boundary layer and the shock wave, is located Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 1390\, mm}
downstream from the leading edge, as can be seen in Figure 17. As described in reference [19], the leading edge can be excluded from the CFD calculations, provided that the conditions of the input limits are slightly adjusted (the Mach number is changed from Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 7.05}
to Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 7.11}
. The wall (cylinder and cone) has a constant temperature of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 311\, K} . The fluid domain and the boundary conditions are shown in the Figure 18.
|
| Figure 17. Solid body (from [18]) |
|
| Figure 18. Domain and boundaries |
Although this is an axisymmetric problem, its resolution is made here with a 3-D mesh, taking a quarter of the geometry by symmetry considerations. In addition, this problem was simulated by the dimensionless codes developed for this work. In the two lateral boundaries is imposed that the normal velocity component is zero. Inlet boundary conditions are:
|
The non-slip condition is applied on solid boundaries. Then
|
The temperature at the wall, which has a known value Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle (311\, K)} , is specified using the dimensionless specific internal energy:
|
and the outlet boundary conditions are specified using Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol{t}}_{j}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{q}}
.
The initial conditions are given by:
|
These initial values are defined in all nodes of the finite element mesh, except for the nodes where Dirichlet boundary conditions were prescribed.
Figure 19 shows the side view of the finite element mesh, which has 1 600 000 hexahedral elements. The mesh is the same as found in the NASA web repository for turbulence models (2-D with 2-zones of 81x101 and 81x101 nodes) but rotated to form a three-dimensional domain mesh corresponding to a quarter of the space surrounding the obstacle with 100 subdivisions in the quarter circumferential length. The total number of elements is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 80 \times (2\times 100) \times 100 . Figure 20 shows a detailed view of the mesh where the flare starts.
|
| Figure 19. Finite element mesh |
|
| Figure 20. Mask in the zone where the cone begins |
The mesh was refined near the wall such that the estimate of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): y+
is close to 0.5.
Using the safety factor Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \beta =0.2} , the dimensionless time step is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Delta t=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {1.0}^{\ast }{10}^{-5} . The artificial viscosity coefficient used here is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {C}_{AD}=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 0.3 .
Figure 21 compares the experimental velocity profiles with those obtained in this work and by the NASA WindUS program. The position in which these data were obtained corresponds to the plane Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): -6
(the origin of coordinates Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x}
is located where the flare starts). Since there are no appreciable differences between the results obtained with TG, CBS and MTG, only the curve corresponding to the MTG scheme is shown here.
|
| Figure 21. Comparison of velocity profiles |
Figure 22 compares the temperature profiles for the same location. Both figures show proximity between the values obtained with the WindUS program [19] and the MTG code implemented in this work. In turn, results of both codes show a good agreement with the experimental data obtained in reference [18].
|
| Figure 22. Comparison of temperature profiles |
Figure 23 shows a comparison of the pressure distribution along the wall of the ogive. The Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle s}
coordinate is measured along the surface of the solid with its origin Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle s=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 0\,
located when the flare starts.
|
| Figure 23. Pressure distribution on the wall |
Finally, Figure 24 shows a comparison of the heat transfer through the wall of the solid obstacle. The definition of wall heat transfer Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {q}_{w}}
is given by: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {q}_{w}=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): -K\, {(dT/dy)}_{w} , where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle K}
is the thermal conductivity of the fluid (but here it is calculated from the corresponding dimensionless quantities).
|
| Figure 24. Heat transfer through the wall |
A good agreement between results obtained with the MTG scheme and "WindUs" codes. As observed from Figure 23, the S-A model tends to slightly over-predict the pressure immediately downstream of the shock. Then the model starts to under-predict this quantity at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle s=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 7
approximately. The S-A model also under-predicts the heat transfer though the solid wall after position Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle s=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 7 . The differences between the results of both codes and the experimental data for positions above Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle s=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 7\,
can be attributed to the adopted RANS/FANS turbulence model (S-A), provided that all references show similar results [19,20].
Like the previous case, the convergence between the TG, MTG and CBS schemes does not differ significantly when the residuals are considered. However, there are differences in the relative wall clock computational times, as shown in Table 3.
| Algorithm | Relative wall clock time |
| Single-step CBS | Not recommended |
| TG | 1.0 |
| TGM | 1.03 |
| CBS | 1.30 |
It is observed that, for hypersonic flows, the three schemes (i.e. TG, MTG and CBS) are suitable from the point of view of convergence rates, but there are appreciable differences in terms of wall clock computational time, being the CBS scheme the most expensive.
The single-step CBS scheme was used for this case showing some numerical oscillations in the flare starting zone, showing again that this algorithm is not recommended for high Mach numbers.
8. Conclusions
Comparing density residuals, CBS and MTG methods show better convergence rate than the TG scheme. But the best convergence rate is obtained using CBS and MTG with local Δtin. The TG method exhibits an oscillatory behavior in transonic and subsonic flows. This behavior is not observed when CBS or MTG schemes are used.
Using a local internal time step with the original TG method, density convergence is not improved, and pressure oscillation remains. However, by making a change in the stabilization term of the continuity equation with a blending function and using a local internal time step, the MTG scheme presents a convergence rate comparable to that obtained with the CBS method.
Although the CBS method does not require an iterative process, as required by the TG and MTG schemes, split operations are performed to calculate first an auxiliary momentum variable, which is used to obtain the density to finally update the momentum values. Then the energy and the S-A variable are calculated. These split operations imply in additional computational time cost. In addition, iterative terms of TG and MTG schemes converge rapidly with very few iterative steps and split operations are not necessary. These features make TG and MTG schemes more efficient than CBS method.
In this way, the MTG scheme is more efficient in terms of computational processing time than the CBS method and with a comparable convergence rate for all ranges of Mach numbers.
The alternative of the single-step CBS method, although efficient and suitable for transonic regimes, is not recommended for high Mach numbers, such as supersonic or hypersonic regimes.
With the modifications proposed here for the TG scheme, a modified method was obtained (the MTG scheme), which is the best option when convergence properties and computational processing times are considered for a wide range of Mach numbers.
References
[5] Zienkiewicz O.C., Taylor R.L. The finite element method. Vol. 3: Fluid Dynamics. Fifth edition, Butterworth Heinemann, Oxford, 2000.
[14] White F. Viscous fluid flow. McGraw Hill Book Co, First edition, New York, 1974.
[15] Huebner K.H., Thornton E.A., Byrom T.G. The finite element method for engineers. Third Edition, John Wiley and Sons, Inc. New York, 1995.
[17] Argyris J., Doltsinis I.S., Friz H. Hermes space shuttle: exploration of reentry aerodynamics. Computer Methods in Applied Mechanics and Engineering, 73:1-51, 1989.
Document information
Published on 26/09/19
Accepted on 19/09/19
Submitted on 25/05/19
Volume 35, Issue 3, 2019
DOI: 10.23967/j.rimni.2019.09.006
Licence: CC BY-NC-SA license
Share this document
Keywords
claim authorship
Are you one of the authors of this document?