Required field
Required field
Required field

# Tutorial: Harmonic Analysis of an Airfoil (2/2)

This article provides a step-by-step tutorial for a harmonic analysis simulation of an airfoil, including post-processing. It is a continuation of the Frequency Analysis of an Airfoil (1/2) tutorial. Animation: Displacement of the airfoil due to harmonic excitation at 25 Hz. Displacements are magnified by a factor of 10.

This tutorial teaches how to:

• Set up and run a harmonic simulation;
• Assign boundary conditions, material and other properties to the simulation;
• Mesh with the SimScale standard meshing algorithm.

We are following the typical SimScale workflow:

1. Prepare the CAD model for the simulation;
2. Set up the simulation;
3. Create the mesh;
4. Run the simulation and analyze the results.

## 1. Prepare the CAD Model and Select the Analysis Type

In case you already have a project from the Frequency Analysis of an Airfoil Tutorial, feel free to continue using that workbench. The geometry and mesh are the same as the frequency analysis tutorial.

Otherwise, you can click the button below. It will copy the tutorial project containing the geometry into your workbench.

The following picture demonstrates what should be visible after importing the tutorial project.

### 1.1 Create the Simulation

As a first step, click on the geometry name.

Hitting the ‘create simulation’ button leads to the following options:

Choose ‘harmonic‘ as the analysis type and ‘create the simulation‘.

## 2. Assigning the Material and Boundary Conditions

To have an overview, the following picture shows the boundary conditions applied for this simulation:

We will assign a pressure load at different eigenfrequencies based on the results of the frequency analysis tutorial as well as fixed support. Furthermore, we are going to use aluminum as a material for the wing.

### 2.1 Define a Material

In the left-hand side panel, click on the ‘+ button‘ next to the materials. In doing so, a list of pre-defined material appears.

Choose ‘aluminium‘ from the library and hit ‘apply‘.

Now you successfully created the material. The blade is assigned to it automatically. All you need to do is confirm it by hitting the checkbox next to Aluminium.

### 2.2 Assign the Boundary Conditions

In the next step, boundary conditions need to be assigned, using Figure 4 as a reference. For an overview of the boundary conditions available, please have a look at this documentation page.

a. Fixed Support

The first boundary condition is a fixed support. It represents zero displacement in all directions. Proceed as below:

After hitting the ‘+ button’ next to boundary conditions, a drop-down menu opens. From it, one can choose between different boundary conditions. After choosing ‘fixed support‘, assign it to one of the airfoil’s tips:

b. Pressure

Still based on Figure 5, create a second boundary condition and choose ‘pressure‘. This will be used for the harmonic excitation. Let’s apply a pressure of ‘500‘ $$Pa$$ on the bottom face of the airfoil. Phase angle should be kept at 0 degrees in case of a single force. Figure 8: Setting up the pressure boundary condition. It is applied to the bottom face of the airfoil.

For applications where more than one force is applied, oftentimes a non-zero phase angle will be specified. For example, in a car engine analysis, we have several pistons with unsynchronized forces. In this case, you can set one of the forces with a phase angle of 0 and define the phase of the other forces accordingly.

Important

For harmonic analysis, at least one load-related boundary condition needs to be defined. Otherwise, there is no prescribed load for the harmonic excitation and the simulation will fail.

### 2.3 Numerics and Simulation Control

Under numerics, change the linear system relative residual to ‘1e-4‘. This value is the maximum relative residual tolerated by the solver.

Under simulation control, you can specify a list of excitation frequencies, at which the loads will be excited. To specify meaningful excitations for the harmonic analysis, it is important to perform a frequency analysis on the geometry, to gain more insight about the natural frequencies.

Therefore, the results from the frequency analysis tutorial are important in this step. These are the first 10 modes for the airfoil geometry, and their respective natural frequencies. You can find the following figure in the results of the frequency analysis tutorial here:

These 10 frequencies span from 4 to 255 $$Hz$$. Let’s focus on the first 3 modes for the harmonic analysis.

In simulation control, change the settings to the following:

• Set excitation frequencies to ‘frequency list‘;
• Change start frequency to ‘2.5‘ $$1/s$$;
• End frequency is ‘30‘ $$1/s$$;
• Frequency stepping will be ‘2.5‘ $$1/s$$. With these settings, a total of 12 frequencies will be analyzed.

### 2.4 Result Control

With result controls, users can monitor parameters in points of interest within the geometry.

To set up a point control, proceed as below:

1. Click on the ‘+ button‘ next to point data;
2. Make sure field selection is set to ‘displacement‘. Component selection is ‘all‘ and complex number is ‘magnitude and phase‘;
3. Lastly, click on the ‘+ button‘ next to geometry primitives to specify a point.

For the coordinates, specify ‘0.1‘ in the x-direction, ‘0‘ in the y-direction and ‘-0.8‘ in the z-direction: Figure 13: Coordinates for the data point. More points can be created, to monitor areas of interest.

## 3. Mesh

To get the mesh, we recommend using the standard algorithm, which is a good choice in general as it is quite automated and delivers good results for the most geometries.

The mesh is the same used for the frequency analysis test. Therefore, change the global fineness to ‘2‘. 2nd order elements will also be used here, to provide more accurate results. Lastly, the gap refinement factor should be set to ‘0.5‘. Finally, ‘generate‘ the mesh.

After about 4 minutes, the 303.3k nodes mesh is ready and looks like this: Figure 15: Mesh of the airfoil. This is the same mesh used for the frequency analysis tutorial.

## 4. Start the Simulation

Now you can ‘start’ the simulation and, after about 35 minutes, the results should be available:

## 5. Post-Processing

First, let’s have a look at the result controls that we specified. Recalling the natural frequencies, we can see that peaks in displacement occur around those frequencies: Figure 18: Displacement on the result control point. For 25 Hertz, the displacement was the biggest.

The phase angles observed for this simulation are either 0 or 180 degrees. This indicates that the peaks in both displacement and forces occur at the same time.

Important

If we were modeling damping, a phase shift would occur, and phase angles between 0 and 180 would also be observed. This is caused by the energy being dissipated out of the system. In SimScale, damping is modeled under materials.

Phase angles between 0 and 180 are also observed when two or more forces are present, with a phase angle between them.

Plotting displacement contours can also provide valuable insight into the behavior of the parts. Below, we have activated the visualization of displacements. For 25 Hz, a bending effect occurs: