Computational Aerodynamics Questions & Answers  




Very good question. We made this assumption only when deriving the GCI. So, within the $\rm GCI_{f}$ equation, you should change the term to $\left\vert\left(\frac{\triangle x_c}{\triangle x_f }\right)^p1\right\vert$. That is, the GCI should always be positive. However, when determining order of accuracy $p$, we did not use the GCI and we did not assume that $p$ should be greater than 1. So you should not change any of the equations used to determine $p$.




The userdefined constant $\epsilon$ is included to prevent a division by zero. Set it to a very small value. As for $\gamma_0$ and $\gamma_1$, they are fixed constants: don't change them.




For the inflow BC, the stagnation pressure can be assumed equal to the one in the freestream because the flow along a streamline is isentropic. But such is not the case for the outflow BC. What if there is a shock somewhere within the domain? Then, the entropy rises and the stagnation pressure will go down and not be equal to the one in the freestream. However, for external flows around a body, the pressure will eventually become equal to the freestream pressure even if shocks are present (as long as the BC is far away from the body). Hence why it's better at the outflow BC to choose to fix pressure rather than stagnation pressure.




OK, will look into this now.




You don't necessarily need to express $\phi$ as a function of $U$ to determine $\partial \phi/\partial U$ in the same way as you don't need to express $F$ as a function of $U$ to obtain $\partial F/\partial U$.



$\pi$ 