case study xv

XV-15 Rotor Blade Analysis using the Blade Element Disk Method

6.9. xv-15 rotor blade analysis using the blade element disk method #, 6.9.1. highlights #.

In this case study we are going to show you how to:

Simulate an XV-15 tiltrotor airplane in both hovering helicopter mode and forward flight airplane mode using the steady BET Disk solver in Flow360

Convert dimensional input quantities to nondimensional quantities in the Flow360 configuration json file

Convert nondimensional quantities in the output to dimensional quantities

6.9.2. Introduction #

The development of the XV-15 tiltrotor was launched by the NASA Ames Research Center in collaboration with Bell Helicopters in the late 1960s and early 1970s. It has become the foundation of many VTOL aircrafts today. Unlike typical helicopter rotors, tiltrotor blades demand a compromised design to operate efficiently in both helicopter and propeller modes. Hence, the blades of the XV-15 tiltrotor have high twist, high solidity and small rotor radius. In this case study, we will apply our Blade Element Theory (BET) Model based solver in Flow360 to evaluate the performance of the XV-15 rotor in both helicopter hovering mode and airplane mode.

Fig. 6.9.1 The Bell Textron XV-15 Tiltrotor aircraft prepares to land (source: United States Coast Guard ) #

The definition of the cross sections of the rotor blade is shown in Table 6.9.1 . The main geometric characteristics of the XV-15 rotor blades are listed in Table 6.9.2 .

The radial chord and twist distribution is shown in Fig. 6.9.2 .

Fig. 6.9.2 XV-15 rotor blade’s twist and chord radial distribution #

More information about the geometry and a high-fidelity detached eddy simulation done using Flow360 can be found at

Feilin Jia, John Moore and Qiqi Wang, 2022, Assessment of Detached Eddy Simulation and Sliding Mesh Interface in Predicting Tiltrotor Performance in Helicopter and Airplane Modes

In this case study, the steady Navier-Stokes solver coupled with the Blade Element Theory Disk method is used to simulate the XV-15 rotor blades under two flight conditions: helicopter hovering mode and airplane mode. In order run cases in a few minutes, compared to the high fidelity unsteady simulations with full rotor geometry in the above paper, the mesh for this BET Disk example is much coarser. The mesh used in the current case study contains only 592K nodes. An overview of a slice of the mesh is shown below

Overview of volume mesh used in BET Disk simulations of XV-15 rotor, successively zooming in from top left->top right->bottom left->bottom right.

The freestream quantities are shown below. These quantities are needed to set up the nondimensional variables in the flow360 case configuration and translate the nondimensional output variables into dimensional values. More information on nondimensionalization in Flow360 can be found at Nondimensional Inputs and Nondimensional Outputs . In the above mesh, the grid unit \(L_{gridUnit}=1\,\text{inch}=0.0254\,\text{meter}\) . In this case study, the freestream conditions are set to standard sea level values as shown in Table 6.9.3

In some simulations, the freestream is not set to be at standard sea level on purpose. For example, in the “case 1a” from 3rd AIAA CFD High Lift Prediction Workshop , the viscosity of the freestream is adjusted to analyze the full-scale geometry at wind tunnel conditions, so please set the freestream properties based on your requirements.

6.9.3. Helicopter Hovering Mode #

In helicopter hovering mode, the freestream velocity is zero. Five blade collective angles are considered in the current study: \(0^o, 3^o, 5^o, 10^o, 13^o\) at \(r/R=0.75\) , corresponding to low, medium and high disk loadings. The flow conditions are:

Tip Mach Number, defined as \(U_\text{tip}/C_{\infty}\) , is 0.69, so \(U_\text{tip}=0.69C_\infty\) .

Reynolds Number (based on reference chord (14 inch) and blade tip speed) = \(4.95\times 10^6\) .

Reference Temperature = 288.15 K.

Here are some points to set up the case configuration file:

freestream->Mach is set to 0, because freestream velocity is zero.

freestream->MachRef has to be a non-zero number because the above “Mach” is 0. Theoretically, this value could be arbitrary, but we set it equal to the tip Mach number (0.69) for convenience.

Either freestream->Reynolds or freestream->muRef should be given to set a proper freestream viscosity. Both options are explained below:

Option 1 : set freestream->Reynolds . The freestream->Reynolds is based on grid unit as the reference length, thus it is mesh dependent. It’s definition is \(\rho_\infty U_\text{ref} L_\text{gridUnit}/{\mu_\infty}\) . In the case description, we know the Reynolds number based on tip speed and reference chord is \(4.95\times 10^6\) , so we must convert the value from the reference chord used to the \(L_{gridUnit}\) in the mesh.

thus the freestream->Reynolds is calculated by:

(6.9.2) # \[\begin{split}\text{freestream->Reynolds}&=\frac{\rho_\infty U_\text{ref} L_\text{gridUnit}}{\mu_\infty} = \frac{\rho_\infty\cdot \text{MachRef} \cdot C_\infty L_\text{gridUnit}}{\mu_\infty} \\ &=\frac{\rho_\infty \left(0.69 C_\infty\right) \text{chord}_\text{ref}}{\mu_\infty}\times\frac{\text{MachRef}}{0.69}\times\frac{L_\text{gridUnit}}{\text{chord}_\text{ref}} \\ &= 4.95\times 10^6 \times \frac{0.69}{0.69}\times \frac{1\, \text{inch}}{14\, \text{inch}} = 3.3536\times10^5\end{split}\]

The freestream->Reynolds is \(3.3536\times10^5\)

Option 2 : set freestream->muRef . The “muRef” here is a nondimensional quantity, it is mesh dependent because it depends on \(L_{gridUnit}\) . Its definition is (6.9.3) # \[\text{freestream->muRef} = \frac{\mu_\infty}{\rho_\infty C_\infty L_\text{gridUnit}}\]

If you already know the Reynolds number then you could use the following equation to get the freestream->muRef :

(6.9.4) # \[\begin{split}\text{freestream->muRef} &= \frac{\rho_\infty U_\text{tip}\text{Chord}_\text{ref}/Re}{\rho_\infty C_\infty L_\text{gridUnit}} =\frac{\rho_\infty \left(0.69 C_\infty\right) \text{14 inch}}{\rho_\infty C_\infty \text{1 inch} \cdot 4.95\times10^6} \\ &=\frac{0.69\cdot \text{14 inch}}{\text{1 inch} \cdot 4.95\times 10^6} = 1.95152\times10^{-6}\end{split}\]

Therefore, the freestream->muRef is \(1.95152\times10^{-6}\)

BETDisks->omega is the nondimensional angular speed of the rotor blades. Its definition is

Because we know \(\Omega=\frac{U_\text{tip}}{\text{Radius}}\) ,

(6.9.6) # \[\begin{split}\text{BETDisks->omega} &= \frac{U_\text{tip}}{\text{Radius}}\cdot\frac{L_\text{gridUnit}}{C_\infty} = \frac{0.69C_\infty\cdot \text{1 inch}}{\text{150 inch} \cdot C_\infty}\\ &= 0.69/150 = 0.0046\end{split}\]

Therefore, BETDisks->omega is 0.0046

After setting up the case configuration file, the case is ready to submit. The 592K-node mesh file, and its corresponding mesh and case configuration files can be downloaded via the following links:

Mesh file: XV15_BETDisk_R150_592K.lb8.ugrid

Mesh configuration file: Flow360Mesh.json

Case configuration file for \(\theta_{75}=0^o\) : Flow360_hover_pitch0.json

Case configuration file for \(\theta_{75}=3^o\) : Flow360_hover_pitch3.json

Case configuration file for \(\theta_{75}=5^o\) : Flow360_hover_pitch5.json

Case configuration file for \(\theta_{75}=10^o\) : Flow360_hover_pitch10.json

Case configuration file for \(\theta_{75}=13^o\) : Flow360_hover_pitch13.json

The only difference between the configuration files of the above 5 cases is the BETDisks/twists distribution.

For detailed instructions on how to upload a mesh, run a case and download the results, please refer to the Quick Start section of this documentation - these details will not be covered in this case study.

Some tips on setting the input quantities related to BET can be found in the BETDisks section of the knowledge base . Please note that there is no hub in the high fidelity detached eddy simulation of the XV-15 rotor blades in our paper mentioned above, so to match the model used in the high fidelity simulations, the chord length in \(r<0.09R\) should be set to 0 in the BET simulation configuration files. I.e. the chord length is set to 0 right before the first cross section ( \(r=0.09R\) ). This setting leads to a radial distribution similar to “chords_distribution_1” shown in BETDisks .

The forces and moments related to the BET Disk are saved in the “bet_forces_v2.csv” file. A detailed description can be found at BET Loading Output . Here we will convert those non-dimensional values into dimensional values:

thrust and torque.

Because the axial direction of the rotor is in the Z axis, the thrust is saved as “Disk0_Force_z” and the torque is saved as “Disk0_Moment_z” in the .csv file. The dimensional thrust and torque can be calculated by Eq.(8.3.2) and Eq.(8.3.3) :

(6.9.7) # \[\begin{split}\text{Thrust} &= \text{Disk0_Force_z}\cdot \rho_\infty C^2_\infty L^2_\text{gridUnit} \\ &= \text{Disk0_Force_z}\cdot 1.225 kg/m^3 \times 340.3^2 m^2/s^2 \times 0.0254^2 m^2 \\ &= \text{Disk0_Force_z}\cdot 91.5224 N\end{split}\] (6.9.8) # \[\begin{split}\text{Torque} &= \text{Disk0_Moment_z}\cdot \rho_\infty C^2_\infty L^3_\text{gridUnit} \\ &= \text{Disk0_Moment_z}\cdot 1.225 kg/m^3 \times 340.3^2 m^2/s^2 \times 0.0254^3 m^3 \\ &= \text{Disk0_Moment_z}\cdot 2.324669 N\cdot m\end{split}\]

The convergence history of the dimensional thrust and torque using steady BET Disk solver is shown in following figures:

Loading convergence of BET Disk simulation in hovering helicopter mode at various pitch angles.

sectional thrust and torque.

In the “bet_forces_v2.csv” file, the sectional thrust coefficients are provided. The process to convert the nondimensional Ct into its physical dimension, i.e. sectional thrust per unit span in SI is very similar to the one shown above. We need the dimensional quantities mentioned in Eq.(8.3.4) and Eq.(8.3.5) to compute the dimensional sectional thrust:

Radius of the rotor disk \(R = 150\times L_{gridUnit} = 3.81\,\text{meter}\)

rotating angular speed \(\Omega = V_{tip}/R = Mach_{tip}*C_{\infty}/R = 61.6237\, \text{rad/second}\)

reference chord \(\text{chord}_{\text ref}=14\times L_{gridUnit} = 0.3556\, \text{meter}\)

If we assume that we want to calculate the sectional thrust and torque at the first disk’s first blade’s second radial location, i.e. Disk0_Blade0_R1:

As a comparison example between the high fidelity full-rotor unsteady simulation and the BET Disk steady simulation, the physical sectional thrust on a blade per unit span, for the \(\theta_{75}=10^o\) case, is shown below

Sectional thrust distribution and history of total thrust in hovering mode, \(\theta_{75}=10^o\) .

The biggest differences between these high fidelity simulations and these BET Disk simulations are near the tip regions where blade-vortex interactions are strong. The flow around the tip can be highly three dimensional, making a BET Disk approximation locally inaccurate. This effect is proportional to the blade disk loading, thus it is larger at hovering or near-hovering conditions when the thrust is higher then in forward flight. Even then, the total thrust of the three blades in hover, compared to the Flow360 high fidelity unsteady simulation, is ~8% different. This level of accuracy makes the BET Disk a useful tool in the preliminary design stages.

To provide another overview of the BET disk accuracy; If we look at the propeller efficiency in hovering mode, the thrust coefficient, the torque coefficient and the figure of merit defined in Eq.(6.9.10) and compare them with several experimental data and numerical prediction of high-fidelity DES simulations, we get the graphs below:

where \(R\) is rotor disk radius and \(A\) is rotor disk area, i.e. \(\pi R^2\)

Comparison on thrust and torque coefficient and figure of merit in hovering mode at various pitch angles.

6.9.4. Airplane Mode #

In airplane mode, four blade collective angles are considered: \(26^o, 27^o, 28^o, 28.8^o\) at \(r/R=0.75\) . The flow conditions are:

Tip Mach Number = 0.54.

Reynolds Number (based on reference chord and tip speed, with no account for the inflow velocity) = \(4.5\times 10^6\) .

Reference chord = 14 inch.

Reference Temperature= 288.15 K.

Advance ratio (defined as inflow speed over tip speed) = 0.337

The mesh file and mesh configuration file are the same as the files used in the helicopter hovering example above. The case configuration files, however, are different:

Case configuration file for \(\theta_{75}=26^o\) : Flow360_airplane_pitch26.json

Case configuration file for \(\theta_{75}=27^o\) : Flow360_airplane_pitch27.json

Case configuration file for \(\theta_{75}=28^o\) : Flow360_airplane_pitch28.json

Case configuration file for \(\theta_{75}=28.8^o\) : Flow360_airplane_pitch28.8.json

The convergence history of the thrust coefficients and torque coefficients using the steady BET Disk solver are shown below

Convergence history of thrust coefficient and torque coefficient in airplane mode at various pitch angles.

As an example, similarly to the helicopter hovering case shown above, the physical sectional thrust for the \(\theta_{75}=26^o\) case is shown below. We have comparisons to the high fidelity DES simulation both on a blade per unit span basis and for the total force.

Please note the non-zero origin of the Y axis on the right plot of total force.

Sectional thrust distribution and thrust history in airplane mode at \(\theta_{75}=26^o\) .

To provide an overview of power efficiency in airplane mode at various pitching angles, the figure below shows the comparison of the thrust coefficient, torque coefficient and propeller propulsive efficiency, defined in Eq.(6.9.11) , between the BET Disk and high fidelity simulations.

Comparison of thrust coefficient, torque coefficient and propulsive efficiency in airplane mode at various pitch angles.

You may also be interested in

NISM- SERIES - XV: Research Analyst

Become a Master of Financial Analysis. Get NISM Series XV Certified

Instructor: Rohit Singh Language: English

About the course

The NISM-Series-XV: Research Analyst Certification Examination is designed to establish a common minimum knowledge benchmark for all individuals associated with research analyst activities under SEBI (Research Analyst) Regulations, 2014. This includes registered research analysts, employed research analysts, and partners involved in the preparation and/or publication of research reports or analysis.

Candidates are required to pass this examination to meet the requirements specified under Regulation 7(2) of the SEBI (Research Analysts) Regulations, 2014, aiming to elevate the quality of services offered by research analysts in the financial services industry.

Examination Objectives:

Upon successful completion of the examination, candidates should demonstrate proficiency in the following areas:

• Understanding the basics of Indian Securities Markets and the terminology used in equity and debt markets.
• Familiarity with both top-down and bottom-up approaches to fundamental research.
• Grasping the basic principles of micro and macro-economic analysis, including sources of information and the various macroeconomic variables influencing analysis.
• Identifying key industry drivers and understanding the sources of information for industry analysis.
• Comprehending the qualitative and quantitative dimensions of company analysis.
• Understanding the fundamentals of Risk and Return, Valuation Principles, and the rationale behind various Corporate Actions.
• Recognizing the characteristics of a well-constructed research report.

Assessment Structure:

The examination comprises 92 multiple-choice questions, each carrying 1 mark, and 2 case-based questions (each case containing 4 questions, also worth 1 mark each), totaling 100 marks. The assessment structure is as follows:

• Multiple Choice Questions: 92 questions (1 mark each) = 92 marks
• Case-based Questions: 2 cases (4 questions each, 1 mark each) = 8 marks

Candidates have 2 hours to complete the examination, with a passing score of 60. Negative marking applies, with a deduction of 25% of the marks assigned to each question for incorrect answers.

Test Details:

Module Name: NISM-Series-XV: Research Analyst Certification Examination

*Note: Payment inclusive of Service Tax. Additional Payment Gateway Charges apply.

*: Negative marking – 25% of the marks assigned to the question

~: 92 independent MCQs of 1 mark each; 2 caselets with 4 MCQs in each caselet of 1 mark each (total of 8 case-based MCQs)

What do we offer

Live learning.

Learn live with top educators, chat with teachers and other attendees, and get your doubts cleared.

Structured learning

Our curriculum is designed by experts to make sure you get the best learning experience.

Community & Networking

Interact and network with like-minded folks from various backgrounds in exclusive chat groups.

Learn with the best

Stuck on something? Discuss it with your peers and the instructors in the inbuilt chat groups.

Practice tests

With the quizzes and live tests practice what you learned, and track your class performance.

Get certified

Flaunt your skills with course certificates. You can showcase the certificates on LinkedIn with a click.

Reviews and Testimonials