Computational Aerodynamics Questions & Answers  




I need more information. Please show a plot of the pressure contours. Also, show your code in the Post() module. Further, how are you invoking warp with opm? I hope you read in your data file... Otherwise, this will simply postprocess the initial conditions..




Everything seems fine except your matlab values. Are you doing a simple averaging using matlab? This will of course be wrong. You need to do a mass flux averaging (do $\frac{1}{\dot{m}}\int P d\dot{m}$ in matlab not $\frac{1}{H}\int P dy$).




The problem may be that you are imposing free stream boundary conditions on all boundary nodes..? You should check your boundary conditions are well imposed using the on flag.




This looks fine: your boundary conditions are well imposed. The contours you showed previously seem to have not been iterated at all. After the iteration process, are you reading in the data file (obtained after 1000 iterations or so) before outputting the post file?




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



$\pi$ 