Computational Aerodynamics Questions & Answers  


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.






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$.




Use a 1D extrapolation polynomial. 2D is too time consuming to compute.




No, $A$ is a $2 \times 2$ matrix determined from an average state function of $Z_L$ and $Z_R$.




I updated the answers to make them more clear.




Hm, I see a problem in your approach. You shouldn't be finding a polynomial when determining the flux with a TVD minmod2 limiter.




I don't understand why using reconstructionevolution would give you a different answer. You need to find $u_{\rm R}$ by interpolating $u$. Once $u$ is interpolated and $u_{\rm L}$ and $u_{\rm R}$ are found, then apply reconstruction evolution.




For $f$ and $u$ given, there is only one wave speed not two.



$\pi$ 