Open VOGEL/Tutorial 1 Part 1

From Wikibooks, open books for an open world
< Open VOGEL
Jump to navigation Jump to search

Modeling wings[edit]

In this tutorial we are going to analyze a flying wing. We are going to start by converting our design into a calculation model and then we are going to take a look at the results to see how we can use them. In parallel, we are going to make some comments about how the software works, and how it was meant to be used.

In Tucan the job is divided in three main steps, namely:

  1. Modeling
  2. Simulation set up and calculation
  3. Post processing of results

For each phase there is a corresponding tab in the main ribbon, so while you evolve in the analysis, you will jump from left to right into each one of them. That will be the usual way of working with the program. As you might have noticed, the HMI has no icons. Instead, we prefer giving each component a visible name describing its functionality. Therefore, in most situations, the HMI is self describing.

Creating and editing a component[edit]

In the modeling phase, we are going to translate our design into a meshed model for numerical analysis. And here is something you have to keep in mind. We are now focusing our attention in the aerodynamic analysis of the model, not in the constructive details, so we will have to leave superfluous details behind. The situation is like this: we have this numerical method based in inviscid fluid theory to test the model, so, how do we adapt our design to make it useful? During the analysis phase, the method is conditioning what we can model. Unfortunately, it is not the other way around. So we have to be clever enough to use the method in our favor.

Before you start up the program, you will need to gather basic information about your design: the dimensions of the model, the kind of airfoils you are going to test along with their polar curves, and so on.

In this tutorial we are going to analyze a tapered wing with the following characteristics:

Property Value
Wingspan 10m
Root chord 1.5m
Tip chord 0.5m
Leading edge sweep-back 20°
Flap length 2.5m
Flap chord ratio 0.3
Airfoil NACA 2411

Adding a lifting surface[edit]

Tucan normally starts up with an empty work-space. To introduce a wing (or any other component) go to the Model flap of the main ribbon, click on Add and then pick Lifting surface in the pop-up window. You will see that a new wing will immediately appear in the model.

The wing we have just added contains default parameters, so we will have to reshape it to our needs. To do that, you first have to select it. You might have already noticed that in Tucan there are two ways of selecting a component:

  1. By clicking on it
  2. By selecting it from the combo box located in the ribbon

The combo box will always contain all of the current models, so if you loose one of sight (something that could happen when you work with fuselages), you can always find it back there.

Editing a lifting surface[edit]

Now let's see how to edit our wing. You will see that this is actually very simple. Just select the wing and then click on Edit (under the name entry box) and you will see that a panel will emerge from the left side of the screen. In general, when you edit a component you will follow the same procedure, with the difference that for some of them you will get an edition window in place of a panel.

Lifting surface edition panel

The lifting surface edition panel contains all of the parameters that you can modify to create a wing model. A lifting surface in Tucan is a collection of consecutive regions that go from the root to the wingtip in span-wise direction. The program will let you add as many of these regions as you want, at the end or in between existent ones. If you click on Add a new region is appended at the tip of the wing. If you click on Insert, a new region is added next to the current one and closer to the wingtip. You can navigate through the different regions using the numeric up-down located next to the label Region. Keep in mind that regions are needed when the geometrical or elastic parameters of the wing change abruptly in span-wise direction, or when you need a variation different than linear. By default, when we add a new wing, we get one region of 1m by 1m. For our specific case, because we are going to model flaps and ailerons, we need two regions. So we are going to add a second one. When you do that, this second regions will be focused, so move the numeric up-down to go back to the first one. Then fill in the properties of the first region. To recreate the model you can see in the pictures, follow the next instructions:

  • Under the flap Geometry, choose 8 chord-wise panels, 1.5m as root chord, 2.5m length, 1.0m tip chord, 20° sweepback and 3° dihedral.
  • Under the flap Flaps, choose flapped, 5° of flap deflection, 0.3 flap to chord ratio, and 3 flapped panels.
  • Then change to region 2, and do the same: 2.5m length, 0.5m tip chord, 20° sweepback and 3° dihedral.

Note that you can play with the different geometrical parameters and you will immediately see how the model is affected. This feature can be very useful if you are designing a model rather than trying to recreate something. Ok, now the wing is almost ready. It is time to define the airfoil properties!

Camber lines[edit]

Our wing cross section is a NACA 2411 airfoil. How do we model this? Tucan will not model the thickness of the wing, so the last two digits are irrelevant for us. We will only focus on the camber line given by the two first digits.

Every time you create a project, the program will internally make room for a small database of profiles, and you can access this database from the lifting surface edition panel, by clicking in the Camber lines button located in the Profile flap of the Region panel. The next picture should make this clear.

Camber line edition dialog

By default, wing regions are created using a symmetric camber line. This type is always there at the top of the list and you cannot remove it or edit. If you want a different one, then you have to add it to the local database and then assign it to the target region by clicking Ok. You have two options to generate a camber line:

  1. You draw it yourself point-by-point
  2. You generate one using the generator

The current version of the software only counts with a NACA 2 digits camber line generator, so if you have a different one, then you have to either load it manually, or develop a new generator in Visual Studio and share it with us.

So for now click Add, and a new item will appear on the list. Select the item and name it NACA 24XX. Then generate a camber line using the NACA generator with e/C=0.02 and x/C=0.4.

Now click Ok, and you will see that the current region has been updated in the model to reflect this new change. Do now the same for the other region. This time, when you open the camber line database form, you will notice that the previously made camber line is already there ready to be used. The idea is that you can actually reuse every camber line you make in any of the regions of any of the lifting surfaces you add. You don't need to redefine camber lines for different regions or wings.

Now we can proudly say that our geometrical model is complete. We will leave the elastic properties for a different tutorial (they are only required for the aero-elastic analysis), and see now how we can take into account the skin drag of our wing using airfoil polar curves.

Polar curves[edit]

Inside the project there is a small database of polar curve families. A family of polar curves is a set of curves that represent the skin drag coefficient of an airfoil as a function of the local lift coefficient and the Reynolds number. Be careful with this. Unlike camber lines, polar curves do depend on the thickness of the airfoil, and therefore you will have to introduce different families for a family of airfoils having the same camber line.

Polar curve edition dialog

You can add as many families as you want, and for each family, you can add as many curves as you need for different Reynolds numbers. When the actual computation of the skin drag takes place, what will happen is that the local chord of the wing will be computed for each chord-wise strip of panels, and based on that length and the provided viscosity, the local Reynolds number will be computed. If the Reynolds number of a given strip is in-between two of the provided polar curves (which will typically be the case) the program will deduce a local skin drag coefficient by linear interpolation of the polar curves.

Tucan lets you introduce parabolic polars and custom polars. With this later type you can use any external source of information, such as XFoil or your own wind tunnel measurements, to prepare the data. You can manually add the points one by one, but I can assure you it will take a lot of time. To avoid this, what you have to do is to write a table in a simple text file editor (such as Kate) with as first column the CL values, and as second column the CD values. These two columns must be separated by a tab character. Then you just copy the data (Crt+A and Ctr+C) and go back to the form, right click on the table and select Get from clipboard.

To try this out I have prepared the next data table for you:

-0.8 0.0175
-0.6 0.0136
-0.5 0.0122
-0.4 0.0112
-0.3 0.0105
-0.2 0.0101
-0.1 0.0099
 0 0.0098
 0.1 0.01
 0.2 0.0102
 0.3 0.0106
 0.4 0.0115
 0.5 0.0124
 0.6 0.0138
 0.7 0.0159
 0.8 0.018
 0.9 0.0205

Final model

Note that you can navigate with the mouse through the curve and every time the cursor hovers over a node, the local Cd/Cl values will be displayed on a label. You can also drag these nodes with the mouse if you need to tune-up the curve manually.

After you have loaded the polar, click Ok and it will be assigned to the current wing section. Don't forget to do the same for the remaining section!

With all this, our model is ready for the analysis. In the next section we will se how to setup the simulation.

Running the simulation[edit]

Now that our model is ready, let's see how to perform some calculations. In this first tutorial we will limit our study to the steady state analysis of the wing, that is to say, at constant stream velocity and after the wake nodes have converged to a final position. We will take a cruise speed, and then analyze the equilibrium and static stability using the flap deflection and a reference weight. What we can do in this case is to run four simulations using different incidence angles and flap deflections at the target speed.

After doing this, we will have sufficient information to write the three longitudinal aerodynamic coefficients as linear functions of the incidence and deflection angles.

Go to the Simulation panel of the main ribbon and set the next simulation parameters:

Steps Interval
30m/s 1.225kg/m³ 0.00001780kg/ms 60 0.02s

The air properties that we have declared here are required to compute the dynamic pressure and the Reynolds number.

Next, you will have to go to the wing editor and under the flap Primitives mark the Convect wake and Trailing edge only check boxes and set the cutting step to 40. This later parameter is required to trim the wake and to get rid of the starting vortex. In general, the total number of time steps will have to be around 10 or 20 times larger than the largest trimming step. In this case you can see that we opted for 20 extra steps after the trimming step. This should provide sufficient steps for the wake and the air-loads to converge to the steady state values.

In order to produce the simulations in the table above, we will have to adjust the incidence and the flap deflection angles. We have already seen how the flap deflection can be adjusted, but we have not mentioned yet how to change the incidence angle. For this there are actually two options, and both should produce the same results. The first one is to change the direction of the stream velocity by adding a vertical component (along the Z axis) and slightly modifying the X component. The second option is to rotate the wing model around the Y axis using the theta angle. For our case now, because we only have one component, this second option is simpler. However, when many components are present, you will see that it is easier to just re-orientate the stream.

So now we are finally ready to hit the Start button! When you do it, the settings dialog will popup. You can see that many of the simulation parameters we have already set before, so here we can take a second look at them to be sure everything is alright. This form also presents some extra options that you might want to use. In our specific case, we wish to turn on the wake extension feature, so that the effect of trimming the wake is mitigated.

The results of the four simulations are given in the next table.

0.16147 0.02466 -0.22074
0.58605 0.03286 -0.67111
0.32856 0.02931 -0.40366
0.75137 0.03684 -0.85224

Using an intel Core i7 processor, each simulation should take around 18 seconds.

(tutorial to be continued)