ONERA M6 Wing Study # 2

Figure:1 Plot of the preasue coefficients on the M6 ONERA wing surface.
Introduction
This study is an example demonstrating the use of PAB3D to compute three-dimensional, transonic, turbulent flow over the ONERA M6 wing.
Download tar File
Most of the files for this example study are available in the Unix compressed tar file Onera.tar . The files can then be accessed by the commands
Grid
This study uses the grid lines from the Plot3d grid file (SGI unformatted, 3D, multi-zone) m6wing.x, which is from the original NPARC validation ONERA M6 Wing example case.
It consists of four zones wrapped as a C-grid about the wing leading edge.
Boundary Conditions
The boundaries at the surface of the wing are specified as no slip wall boundary conditions. The inflow boundary is specified with a constant total pressure, total temperature and flow angle. The farfield boundaries are specified as characteristic boundaries. The outflow boundary is specified as an extrapolation boundary. The boundary conditions are set within the control file tpab3d.cont. Procedure for setting the case, USER can refer to Tutorial.
Initial Conditions
The initial flow conditions are simply the freestream flow values as presented in Table 2. These corresond to a Reynolds number of 11.72 million based on the mean aerodynamic chord length.
Table 2. Freestream conditions.
|
Mach |
Pressure (psia) |
Temperature (R) |
Angle-of-Attack (deg) |
Angle-of-Sideslip (deg) |
|
0.8395 |
45.82899 |
460.0 |
3.06 |
0.0 |
These conditions correspond match the Mach number, Reynolds number, and angle-of-attack of test 2308 from the paper by Schmitt and Charpin in AGARD Report AR-138.
Computation Strategy
The computation is performed using the time-marching with mesh sequencing capabilities of PAB3D to march to a steady-state (time asymptotic) solution. Local time stepping is used at each iteration. The time-marching is performed until a convergence criterion is achieved.
Input Files and Computation
Step 1: Preprocessing of grid file can be done two ways; either by running autog3d program or using graphical user interface (UPMS).
If user prefers to use autog3d, he must change the values in control files as needed by user modeling requirements.
If user prefers UPMS graphical user interface, after completing the all needed steps (Ref: Tutorial), he/she must save the project and export this project. This prompts user to choose the solver from three codes (TLNS3D, PAB3D, CFL3D) available. Once user selects PAB3D and clicks OK, it presents user with interface as shown below

User should enter names of control files. The User filename must be user.cont. User can give name to G3D filename. For this case, we name the file: g3d.cont
Step 2: Run Pre2 to generate the database information and tpab3d.cont control file. This can be done by issuing the following command as shown below
Pre2 g3d.cont
Step 3: Keep all the generated control files and grid file in one subdirectory. Now run PAB3D code by issuing the following command.
On single processor: Pab2
or
On multi processor: mpirun –p4pg hostfile Pab2
hostfile should contain the information about the cluster information. For example, if some one wants to run the code on a cluster machine with machine names e2, e3, e4, user should create a hostfile with the information as shown below.
E.g.: hostfile
e2 0
e3 1 ~hamid/bin/Pab2
e4 1 ~hamid/bin/Pab2
Here, executable Pab2 is located in sub-directory bin of hamid home directory.
Step 4: Once results are obtained, output files can be processed by using POST and tecplot. Refer POST tutorial for post processing.
Comparisons of the Results
The pressure coefficients at the several sections on the wing are examined and compared to experimental values. Comparisons were made at seven sections: 0.2, 0.44, 0.65, 0.8, 0.9, 0.95, 0.99.
The corresponding experimental profile data for these sections were taken from the file ONERAb114.tec.







Performance
The computations are performed on Pentium 3 with 800MHz clock speed and two-processor capability. Convergence of solution can be examined from residual history plots (Fig. 10 and 11).


Acknowledgments
The design of this web page and the grid for this test case were based on the ONERA M6 Wing validation Study # 1 or the Wind code by John W. Slater of NASA LeRC
References
Schmitt, V. and F. Charpin, "Pressure Distributions on the ONERA-M6-Wing at Transonic Mach Numbers," Experimental Data Base for Computer Program Assessment. Report of the Fluid Dynamics Panel Working Group 04, AGARD AR 138, May 1979.
John W. Slater, M6 ONERA Wing Study for Wind Validation WWW site.