Introduction to CFD Questions & Answers  




Here xcut is the xstation you specified in the Post() module. The program can't find a proper xstation at 0.1 m. If you are solving Problems 1 or 2, you don't need to create xstations, so get rid of the code related to the xstations within Post(). If you are solving Problems 3 or 4, you need only one xstation fixed to the $x$ just slightly before the domain exit.




As I explained today, ximax probably won't go down to zero for your problems if you use the TVD Van Leer limiter or WENO. The residual will thus “hang”. So, you can't rely on ximax to determine convergence. You need to proceed differently as explained in class.




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



$\pi$ 