Computational Aerodynamics Questions & Answers | |
|
|
|
|
Your grid has issues. You should make the grid so that the spacing is more or less uniform everywhere. Don't worry about making the grid parallel to the shocks. The shocks can go through the cells at an angle — this is fine. Simply make sure that your bottom wall has the right shape, that the cells are more or less of the same size everywhere (and are not distorted as in your mesh), and the shocks will appear correctly positioned.
|
|
|
|
The instructions are in the CFDWARP HOWTO here: |
|
|
|
Please attach a picture of the pressure contours around your cylinder. Maybe your flow is not yet at steady-state. Because you are doing the cylinder case, you don't need the xstation. Inside the Post() module, delete the part about the x-station and only keep the line Fpressure[dim]=_Fpressure(... Then, add a bit of code to find the drag coefficient using Fpressure[1] (the pressure force acting on the body along $x$).
|
|
|
|
Problem 4 is harder than the others. You need to find $\sigma$ and $\delta$ for each shock using the obliqueshock program. Once this is done, you need to find the bottom surface of the ramp that will make all shocks meet each other at one point (such is only function of $H$, and the 3 $\delta$s and $\sigma$s). Thus, this requires a bit of geometrical calculations on paper. I recommend another problem if this is taking too long.
|
|
|
|
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.
|
|
||
$\pi$ |