MSES

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 1
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
MSES

MSES
Developed by Professor Mark Drela
at the Massachusetts Institute of Technology
gy ((MIT))

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 2
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
MSES

MSES is a coupled viscous/inviscid Euler method for airfoil design and

analysis featuring:

Subsonic, Transonic and Supersonic Single- and Multi-Element Airfoil Analysis

 Forced or free boundary layer transition
 Transitional separation bubble modeling
 Lift and drag predictions to just beyond CL max
 Blunt trailing-edge treatment

Polar and Mach Number Sweeps Mode

 Program steps through angles of attack
 Program steps through Mach numbers while holding lift or angle of attack constant

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 3
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
Input files

blade.xxx – airfoil geometry and grid domain

 Created by
y user NAME
XINL XOUT YBOT YTOP
 Used in MSET and AIRSET programs X(1,1) Y(1,1)
 „xxx“ – variable extension X(2,1) Y(2,1)
X(3,1) Y(3,1)
 XINL – x-location of the left grid inlet plane . .
. .
 XOUT – x-location of the right grid outlet plane X(i,1) Y(i,1)
 YBOT – y
y-location
location of the lowest grid streamline 999.0 999.0
X(1,2) Y(1,2)
 YTOP – y-location of the topmost grid streamline X(2,2) Y(2,2)
X(3,2) Y(3,2)
. .
. .
X(i,2) Y(i,2)
999.0 999.0
X(1,2) Y(1,2)
. .
. .

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 4
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
Input files

mses.xxx – runtime parameters

 Used in solver programs MSES, MSIS, MPOLAR GVAR(1)
( ) GVAR(2)
( ) … GVAR(N)
( )
and
d MPOLIS
GCON(1) GCON(2) … GCON(N)
 GVAR(1) … GVAR(N) – list of integers specifying the
MACHIN CLIFIN ALFAIN
global variables to be used
ISMOM IFFBC
 GCON(1) … GCON(N) – list of integers specifying
the global constraints to be used REYNIN ACRIT
 MACHIN – free stream Mach number XTRS1 XTRP1 XTRS2 XTRP2 …
 CLIFIN – specified CL (if CL option invoked in GCON) MCRIT MUCON

 ALFAIN – free stream angle of attack

 ISMOM – momentum/entropy conservation
 IFFBC – farfield boundary conditions
 XTRSn, XTRPn – top and bottom surface transition
strip x/chord location on element n
 MCRIT – critical Mach number above which artificial
dissipation is added
 MUCON – artificial dissipation coefficient

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 5
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
Executing MSES

mrun xxx
 With “xxx” extension to distinguish
g case being
g run
 For repetitive program execution
 Menu for executing all programs

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 6
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
AIRSET

Airfoil geometry manipulator

 Translate elements
 Rotate elements
 Scale elements
 Split off flap
 Modify contour
 Add/delete corner points

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 7
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
MSET

Grid and flow initializer

 Options
p 1,, 2,, 3,, 4 must be issued
as a minimum.
 1 – specifies angle of attack to
generate a panel solution to trace
stagnation
t ti lines
li and
d th
the upper and
d
lower farfield streamlines.
 2 – initial surface gridding.
Distributes grid nodes along the
streamlines on the element
surfaces. Intermediate streamline
nodes are generated in the
flowfield interior.
interior
The resulting grid is not yet
suitable for MSES calculation!

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 8
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
MSET

Grid smoothing
 3 – elliptic
p g grid smoothing
g of linearly
y
interpolated grid. Eliminates all kinks,
overlaps and makes all grid streamlines
correspond to streamlines of inviscid
incompressible
p flow.
 4 – writes out the initial solution to file
mdat.xxx, which is then ready for the MSES
solver.
 The grid can be viewed with option 5.
Options 2,3 or 1,2,3 can be repeated to
obtain acceptable grid.

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 9
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
MSES

Flow solver
 Solves the Euler equations.
 Has two input files: mdat.xxx and mses.xxx.
 User specifies “n” – maximum number of iterations to be executed.
 Possibility to coarse/refine the grid to speed up the convergence by typing “-” or “+” at the
iteration number prompt.
prompt
 After each iteration, MSES prints out a table of flow parameters change (r.m.s. and max)

 The number of iterations specified before is only a limit.

 MSES will terminate early if the r.m.s. drop below the convergence tolerance limits and the
maximum
i changes
h d
drop 10
10x ththese lilimits.
it
 Writes the solution back to mdat.xxx

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 10
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
MPLOT

Solution plotter
 Displays
p y solution in mdat.xxx whether it is converged
g
or not.

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 11
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
MPOLAR

α-sweep calculation driver

 Conveniently sweeps through a range of specified angles of attack.
 Takes advantage of the quadratic convergence of the Newton method.
 Much more efficient than running a sequence of independent cases with MSES.
 Requires an alfas.xxx file, which contains the sequence of angles of attack.
 Writes
W it outt the
th integrated
i t t d parameters
t to
t polar.xxx
l and
d the
th surface
f pressure distribution
di t ib ti and
d
boundary layer variables to polarx.xxx
 Results from polar.xxx can be plotted using program PPLOT.

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 12
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
MPOLARC, MPOLARM

CL-sweep calculation driver

 Essentially the same as MPOLAR.
 Converges on specified lift coefficients.
 Requires a cls.xxx file, which contains the sequence of lift coefficients.
 Output files polar.xxx and polarx.xxx have same format as for MPOLAR.

Mach-sweep calculation driver

 Generates a Mach sweep in a manner similar to a polar.
 Requires a file
f machs.xxx containing the Mach numbers to be executed.
 Results are written to files msweep.xxx and msweepx.xxx

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 13
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
Results Comparison

Comparison with VZLÚ

wind tunnel results:
 Cruise configuration
 MSES overpredicts CL max
and αcrit
 Overpredicts dCL/dα
 Good drag agreement

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 14
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
Results Comparison

Comparison with VZLÚ

wind tunnel results:
 Take-off configuration [0;0]

 Flap deflection 20°

+ykl
 CL max comparison +kl
 Disagreementt caused
Di db
by
+xkl
wrong αcrit prediction

Wind Tunnel MSES

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 15
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
Other Features

Computational Design of Multi-Element Airfoil Systems

 Mixed inverse design mode
 Modal inverse design capability
 Modal geometry perturbations

Multi-Point Optimization of Single- and Multi-Element Airfoils With Geometric

Thickness and Area Constraints

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 16
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic
MSES

THANK YOU FOR YOUR ATTENTION!

VZLU CESAR Training Workshop

Aeronautical Research and Test Institute March 18th - 19th, 2009 17
Beranových 130, 199 05 Prague, Czech Republic Prague, Czech Republic