Introduction to CFD Questions & Answers  
Question by Student 201327103
professor, i have some problem to run 'warp -r control.wrp -i data.01 -opm'. Terminal say 'could not interpolate at j=0 k=0 for xcut=1.000E-01 m..' what this mean?
05.29.18
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.
Question by Student 201327133
Professor, i have a question about ximax. You showed us that ximax being similar as zero when state is steady. I run the code about Question #3, and i've got the result about time,iter,etc... but there is no more less than 1.0E+3 at ximax. but i can't believe that it could be zero... is it correct result? or are there more exact solution?
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.
Question by Student 201227111
Professor, I don't understand what the value from -opm command mean exactly. At first, I thought it was the average of values at x=x_station and y=(node1 to node 60). So, I extract data at i=2,10,20,30,40,50,60 and j=1 to 60 using "-pr" command. And averaging each value, And compare to value obtained from "-opm command" with varying x_station. (compare density, Temperature, velocity) When i is small, "-opm" value is similar to my average value. But the bigger the i, the bigger the difference. When i=60, "-opm command" value is similar to average value of i=60, j=1 to 20(only 1/3 nodes). So, I'm stuck here, I don't what to do after. Sorry about my poor explanation.
05.30.18
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 post-process the initial conditions..
Question by Student 201227111
Professor! I bring more explain about previous question. I wanted to know the meaning of Post() modules command to compute mass flux average stagnation pressure.
plot.png  ./download/file.php?id=4214&sid=00431da0329bd07f22f1686717f1935d  ./download/file.php?id=4214&t=1&sid=00431da0329bd07f22f1686717f1935d
I simply gridding like this.(ie=60, je=60, dx is constant) Result is like this.
And My post code is same as original, only changing x_station[1] set to L-1e-5 (L is length of inlet)) and for loop range set 1 to 1.
node60_value.PNG
And I use "warp -r design3.wrp -i data.1500 -op post60.1 -pt gnuplot -pr 60 1 60 60" this command, I get properties that x is fixed in one line and y is vary.(like above picture) And, use matlab code, I compute average of each properties($\rho$, $T$, $q$). Then, I compare this average value to value that obtained from "warp -r design3.wrp -i data.1500 -opm". (Because properties vary along y axis but using "warp -r design3.wrp -i data.1500 -opm" is give only one value)
opm_value.PNG
then my computed average value is $\rho=0.2774$, $T=467.5837$, $q=2403.5$, but value obtained by "-opm command" is $\rho=0.6730$, $T=789.1027$, $q=2297.65$.
I think value obtained by Post() module is average of value on fixed x line, but it isn't. It only work when i is very small(only 2 or 3..). So I want to understand about meaning of value that obtained by Post() module. or I think something wrong?
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$).
Question by Student 201327103
professor this is my counter plot for cylinder.
9.png
I can't get shock wave. I Think my boundary condition is wrong, but i don't know where is wrong. here is my code
05.31.18
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.
Previous   1  ...  13 ,  14 ,  15    Next  •  PDF 1✕1 2✕1 2✕2
$\pi$