Influence of elastic deformation of the blade on the flow field in the blade cascade

Flutter of turbine and compressor blades represents a serious problem for designers of large turbomachines. In real machines, measurements of flutter conditions are hardly possible. Therefore, tests on linear blade cascades with movable blades play important role in investigations of flutter. In the current research, a simple blade cascade for controlled flutter testing was developed. The blades in this test rig undergo high-frequency oscillations which induce inertial forces. The influence of the elastic deformation of the blade on the flow field is studied in this work by means of experiments and numerical simulations. First, the computational and experimental modal analysis was done to obtain eigenfrequencies and modal damping. The deformed shape of the blade due to high-frequency oscillation was acquired by structural transient analysis. The influence of the elastic deformation on the flow field was then studied by CFD analysis for two incidence angles both for the deformed (elastic) and undeformed (idealized rigid) blade. Flow field was only very weakly distorted due to the blade elastic deformation. The total torque induced by aerodynamic and inertial forces was evaluated. The inertial loading is an order of magnitude larger than loading due to fluid flow.


Introduction
Blade flutter in turbines and compressors represents a serious problem for designers and operators of large turbomachines.Lack of experimental data constrains faster progress in this area.Flutter of turbine and compressor blades is defined as flow-induced vibration where the vibration frequency is not integer multiple of rotational frequency.A similar phenomenon, nonsynchronous vibration (NSV), occurs at higher incidence angles and higher reduced frequencies [1,2].
Flutter in turbomachines appears due to the interaction of unstable aerodynamics and bending or torsional vibrating mode [3].Growth of oscillation amplitudes that may lead to blade fatigue failure occurs when aerodynamic damping is negative, i.e. when the work per oscillation cycle is positive.When flutter occurs, the oscillation of blades in bladed disk assembly are not in phase but are shifted.This shift is expressed by the Inter Blade Phase Angle (IPBA) and the travelling flutter wave can be observed.The initiation of flutter travelling waves is not yet fully explained [4].
In the past years, a considerable number of studies that focused on wind turbine flutter have been published, e.g.[5][6][7][8].With the growing number of wind and solar power plants, turbines in standard power plants must be operated in wider range of conditions to compensate the fluctuation of supply by renewable sources.The turbines are then operated in off-design conditions and flutter in low-pressure stages, where long and slender blades are installed, can occur [9].The problem of flutter becomes even more severe due to modern design of compressor blisks with negligible mechanical damping and increased structural coupling [2].
As flutter is an aeroelastic phenomenon, the structural and aerodynamic properties of the blades can be tuned to avoid flutter.Structural parameter of high importance is the modeshape.The reduced frequency k = b∕u , where = 2 f is the circular frequency, b the semi-chord of the blade and u is flow velocity, is both a structural and aerodynamic parameter.An important aerodynamic property is the incidence angle -flow separation can occur when the incidence angle is deflected to off-design conditions which can lead to stall flutter [3].
Many attempts to provide "flutter free" design of blades were carried out.Improvement of the flutter stability of a fan blade was presented by Stapelfeldt and Vahdati [10].They suppressed flutter by altering the radial distribution of stagger angle which led to lower efficiency.As an alternative approach they reduced flutter by drawing bleed air from the casing above the tip of the blade.It is known that well-tuned symmetric designs contribute to the fluid-structure coupling [2].In 2000, Nowinski [11] demonstrated benefits of blade mistuning.Mistuning has beneficial effect on blade flutter but unwanted effect on forced response-due to mistuning the forced-response oscillations are present in wider range of frequencies.Therefore, it is required to perform optimization on bladed disk for both forced response and blade flutter.Practical optimization method was proposed by Kaneko [12] or Raeisi [13].Vanti et al. [14] considered numerical procedure to optimize both efficiency and flutter-free design of the compressor blade-row.In 2018, Corral et al. [15] experimentally validated the benefits of intentional alternate mistuning of blades by flutter free tests.They showed that mistuning can suppress flutter and also modulate vibration amplitude.Geometric optimization of the low-pressure turbine blade airfoil to enhance stability of the turbine was realized by Peeren and Vogeler [16] and the design recommendations were given.Tani et al. studied possibility of flutter suppression by change of torsion axis position [17].
Although the development of numerical methods and experimental testing was intensive in the past decade, the problem of blade flutter has not yet been fully resolved.There is a need to further improve aeroelastic design tools to prevent self-excited vibrations.CFD predictions of aeroelastic instabilities are possible but not always reliable, especially for transonic flow regimes [2].Numerical simulations have to be supplied by experimental testing.
Even with improvements in the measurement techniques it is not possible to carry out flutter research on real turbines [18].Simplified experimental test setups have to be employed.Experimental test setups for flutter measurements can be (a) linear cascades and (b) annular cascades.Annular cascades are more realistic, but linear cascades offer easier instrumentation of measuring techniques.The motion of the test object can be flow-induced (free flutter testing) or controlled (controlled flutter testing).In case of free flutter testing, the flow conditions are modified until flutter develops.Controlled flutter testing is often used in research.Oscillations of the blades are prescribed and the aerodynamic response is measured.From these data, aerodynamic damping can be evaluated [3].
A new test facility for controlled flutter testing was developed in the Institute of Thermomechanics of the Czech Academy of Sciences in cooperation with Technical University of Liberec, see Fig. 1.This advanced experimental setup was developed upon experience gained while working with the NASA Transonic Flutter Cascade [19].The cascade in a transonic suck-down wind tunnel consists of five blades with a simplified flat geometry (see Fig. 2).The rotational oscillation of the middle blade can be controlled by a mechanical drive.The oscillation of the middle blade is monitored using a laser triangulation sensor, and the pressure distribution over the blade is measured using miniature unsteady pressure transducers.Detailed description of the test facility can be found in [20].Using this facility, an experimental study investigating the limits of applicability of the quasisteady approximation was carried out in [21].
A similar transonic linear cascade test rig was constructed at the Royal Institute of Technology (KTH) by Glodic et al [22].The oscillations of blades in a traveling wave mode are induced by piezoelectric actuators.These actuators allow oscillation frequencies f = 1−2.5 kHz which corresponds to reduced frequencies k = 2 − 4 .Another transonic linear cas- cade was developed in Indian Institute of Science [23].The oscillations of central blade are driven by barrel cam mechanism up to frequency 250 Hz.
Subsonic linear cascade for controlled flutter measurement is present at University of West Bohemia in Pilsen [18].Annular test rigs for controlled flutter measurement were developed at KTH (subsonic annular section) [24] and at École Polytechnique Fédérale de Lausanne (annular, up to M = 1.4 ) [25]).Test rigs for controlled flutter testing that the author is aware of are summarized in Table 1.
Tall slender blades of modern turbomachines have first flap and torsional eigenfrequencies up to 500 Hz.Therefore, the forced oscillation of blades in the current test facility has to be in identical frequency range.This produces high inertial forces acting on the blade, and therefore it is necessary to take into account the elastic deformation of the blade.Consequently, the influence of the elastic blade deformation onto the flow field and its two-dimensionality has to be assessed.Knowledge of inertial and aerodynamic loading is also necessary for the structural design of the blades in this experimental facility.
This study utilizes experimental and computational modal analysis to determine the eigenfrequencies and damping of the blade, providing insights into its dynamic behaviour (Sect.2-Modal analysis).In Sect. 3 Transient structural analysis is performed to obtain the deformed shape of the blade under highfrequency oscillations.Section 4 CFD analysis compares the flow around the deformed blade with that around the ideally rigid (undeformed) blade, evaluating the influence of elastic deformation on the flow field.These combined techniques contribute to a better understanding of the structural and aerodynamic aspects of the blade cascade for controlled flutter testing, facilitating the development of robust test blades that accurately simulate real-world conditions and enhance our understanding of flow-induced vibrations in turbomachinery.

Modal analysis
In the current experimental setup, where the middle blade undergoes forced oscillation, the natural frequencies of the blade need to be known.When the oscillation frequency approaches resonance, the structural stresses could outreach the material strength, or the elastic deformations of the blade might exceed the geometric limits given by the dimensions of the sidewall slots.The eigenfrequencies of the blade can be quite easily obtained by computational modal analysis.However, the numerical simulation of stresses and deformations during forced oscillation is dependent on the damping, especially near resonance.The damping of the blade can be identified only experimentally, because the overall damping is dependent  EPFL [25] 1983 Annular [11] Up to M = 1.4 NASA [19] 1986 Linear Transonic KTH [24] 2005 Annular section Subsonic UWB [18] 2018 Linear Subsonic IISC [23] 2019 Linear Transonic KTH [22] 2020 Linear Transonic CAS [20] 2020 Linear Transonic on the internal damping of the material, geometry of the structure and friction damping.For that reason, free experimental modal analysis was performed.
The structure in free modal analysis is unconstrained, placed on the soft pad with damping an order of magnitude higher than measured structure.

Computational modal analysis
Two types of computational modal analyses are performed: free and fixed.The equation of motion for an undamped system, expressed in matrix notation using the linear elastic material model assumptions, is given by where [M] is the mass matrix, u is the displacement vector and [K] is the structural stiffness matrix.For a linear system, free vibrations will be harmonic of the form: where i represents the mode shape of the ith natural frequency (eigenvector), i is the ith natural circular frequency in radians per second, and t is time.This leads to an eigenvalue problem which may be solved for up to n values of and n eigenvectors satisfying equation ( 1), where n is the number of DOFs.The main dimensions of the blade are: length l = 341 mm , width w = 120 mm and thickness (as well as radius of leading and trailing edge) t = 5 mm .Material of the blade is high-strength stainless steel with Young's modulus E = 210 GPa , Poisson ratio = 0.3 and den- sity = 7850 kg m −3 .
Free modal analysis of single blade was done for comparison with free experimental modal analysis.The fixed modal analysis was performed to evaluate the eigenfrequencies of the blade in the test rig.The boundary conditions for fixed modal analysis for displacement u were u ⋅ n = 0 in the bearings, where n is unit outer normal, and u = 0 in the hexagonal shoulder.For the free modal analysis there are no boundary conditions, just the six rigid-body motion modes with zero natural frequencies have to be removed from further analysis.For both free and fixed modal analyses, no loads are applied.The governing equations were solved numerically in ANSYS Mechanical on a mesh consisting of 590,000 elements (see Fig. 3). (1) Oscillation frequency of the blades in the test rig should not be near the eigenfrequency of the first torsional eigenmode to avoid resonance.In the fixed computational modal analysis, the first eigenfrequency (torsional eigenmode) is 270 Hz and the eigenfrequency of the second (bending) eigenmode is 541 Hz.The results of computational modal analysis are summarized in Table 2 and corresponding eigenmodes are depicted in Fig. 4.

Experimental modal analysis
The experimental modal analysis was performed in the free configuration, which allows measurement of blade out of the test rig, see Fig. 5. Therefore, the constraints are not considered, which means that the frictional damping is neglected.This simplification keeps vertical displacements, predicted by the transient structural simulation, on the safe side.
The blade was arranged on two flexible elements and a mesh of 15 measuring points was prepared.Measurement apparatus consisted of impact hammer Brüel &Kjaer 8206-002 and the miniature accelerometer Brüel &Kjaer 4394 which was placed in point 14 (grey dot in Fig. 5).In every measuring point, five measurements (impacts) were performed and the average transfer function was computed.The modal parameters (eigenfrequencies, modal damping and eigenmodes) of the system were determined by regression calculation in the frequency range up to 1600 Hz.
Experimentally determined first torsional eigenfrequency (2nd overall) is 711 Hz compared to the computational 694 Hz, which is suitable agreement.Measured modal damping ratio for the first torsional eigenmode is 0.03%.Results are summarized in Table 3.
Experimental modal analysis provides modal damping ratio , which is used to compute constants and for the Rayleigh damping model used in the structural transient simulation.In this simulation, the damping matrix B is estimated as M and K are the mass and stiffness matrices and con- stants and are given by , where 1 = 2 f 1 , 2 = 2 f 2 are angular eigenfre- quencies, f 1 = 313 Hz and f 2 = 711 Hz are eigen- frequencies obtained from experimental free modal analysis.

Transient structural analysis
Structural transient analysis was performed to obtain the deformed shape of the blade due to inertial forces.The equation of motion for damped system expressed in matrix notation is given by: where [M] is the mass matrix, u is the displacement vector, [C] is the damping matrix and F(t) is load vector.
The boundary conditions for displacement u were in the bearings.In the hexagonal shoulder, harmonic rotational displacement was prescribed.The traction boundary condition for all other surfaces is where is the stress tensor.Same material constants as in Sect.2.1 are employed.Linear elastic material model was used.The damping of the structure was described by Rayleigh damping model (3).Computational mesh is identical as in Sect.2.1.The timestep was t = 0.0001 s .A computational time of T = 0.1 s , i.e. 20 periods of oscillation, was necessary in order to pass the initial transient and reach steady-state oscillations.
The governing equations are solved numerically using the Finite Element Method in ANSYS Mechanical by implicit solver based on Newmark time integration method on a mesh consisting of 590,000 elements (same as in Sect.2.1, see Fig. 3).Mesh convergence study was performed and the observed parameter was average local von Mises stress (Fig. 6 left).The results do not significantly vary if mesh element size is lower than h = 1.5 mm (Fig. 6 right). ( The element size h = 1.25 mm was used for calcula- tions.Anticipated frequency of forced oscillation f = 200 Hz and angular amplitude 0 = 1 • were selected.

Results of the transient analysis
Figure 7 shows the temporal evolution of the vertical displacement at the rear end of the leading edge (see point A in Fig. 2), evaluated both for the idealized rigid and real elastic blades.The oscillations of the elastic blade at time t = 0.08 s are assumed to be stabilized at the maximum value of vertical displacement circa 3 mm.Out of this value, 1 mm is due to rigid rotation, and 2 mm corresponds to elastic deformation caused by inertial forces.
The deformed geometry of the blade obtained by transient structural analysis, which is further used in the CFD simulation, was taken at two times.For the positive incidence angle of the blade the time is t 1 = 0.0937 s , which corresponds to rotation angle of the hexagonal shoulder = 1 • .And for negative incidence angle of the blade the time is t 2 = 0.096 s , corresponding to rotation angle of the hexagonal shoulder = −1 • , see Fig. 7.

CFD analysis
The goal of the CFD simulations is to evaluate the influence of the blade elastic deformations on the flow field.The difference of the geometry of the rigid and elastic blade twisted due to inertial forces is depicted in Fig. 8.The flow field is numerically computed for the blade motion frozen in four limit configurations which are summarized in Table 4. Configurations are combination of positive/negative incidence angle and rigid/elastic behaviour of blade.

Geometry and mesh
Similarly, as in [26], the cascade (see Fig. 10) has a slope of 31.5 • and pitch 74.52 mm .The span of the blade exposed to air flow is 160 mm .The computa- tional mesh is composed of hexahedral elements in the free stream, with a structured layer capturing the boundary layer of the blades.The values of y+ on the surface of the blades range between 3-6.The refinement is not realized towards the sidewalls, therefore the side boundary layers are not well resolved.Even with these compromises necessary to keep the computational cost of the simulation reasonable, the mesh consists of up to 16,000,000  elements.The mesh dependence was investigated, see Fig. 9. Based on this analysis, the "Fine mesh", shown in Fig. 10, was selected for all simulations.
The computation took 15 h on a workstation with 8 cores Intel Xeon W-2245 3.90 GHz.

Mathematical model
The airflow is modeled as a steady flow of a compressible ideal gas described by compressible Navier-Stokes equations: where is the density, u is the velocity vector, t is time, f is density of external volume forces, = −pI + (∇ ⋅ u)I + 2 D is the stress tensor, p is pressure, I is identity tensor, = − 2 3 is second viscosity coefficient, is dynamic viscosity coefficient, D is rate of deformation tensor, E is the total energy per unit volume, q is volume heat source and q is the heat flux vector.
Air is considered as an ideal gas with specific heat capacity c p = 1005 J.kg.K −1 , molar mass M = 28.966kg.kmol −1 and dynamic viscosity = 1.81 × 10 −5 Pa ⋅ s .Simulating a suction-type wind tunnel with atmospheric entrance conditions, the boundary conditions at the inlet are the total pressure p tot = 100, 000 Pa , total temperature T tot = 300 K , tur- bulence intensity T u = 2 % , turbulent viscosity ratio T = 10 , zero velocity gradient u n = 0 and zero heat flux

Numerical solution
The governing equations ( 9), ( 10) and ( 11) are solved numerically using the Finite Volume Method by software ANSYS Fluent 2021 R2.The computational setup was: implicit density based solver, k-omega SST turbulence model with log-law not employed.The solver turns on wall functions in wall-adjacent cells where y+ < 11.2.Zero initial conditions were set and full multigrid initialization was employed as a pivotal step in preparing the initial conditions for the simulation.Discretization schemes of the second order were used: second order upwind discretization for the convective term, turbulent kinetic energy and specific dissipation rate, and least squares scheme for the diffusion term.Euler implicit discretization in time was used and CFL (Courant-Friedrichs-Lewy) condition is used to compute the time step.The computations were stopped when residuals of continuity, energy, turbulent kinetic energy, specific dissipation rate and velocities were lower than 1 × 10 −5 .
Fig. 9 Influence of the computational mesh on the isentropic Mach number distribution over the chord of blade 3 in NEG-RIG configuration  The CFD model was validated using the experimental data obtained during wind tunnel measurements for outlet isentropic Mach number M 2is = 0.5 .Experi- ments were carried out in test section depicted in Fig. 1.For these measurements, all blades were kept stationary.The experimental results were compared to the CFD calculations with idealised rigid blade in plane P2 (see Fig. 12).In the experimental setup, the static pressure is measured in 10 positions along the blade chord.Further details about the measurement setup can be found in [21].Isentropic Mach numbers are computed from measured static pressures using (12).Comparison of CFD results for idealised rigid blade and experiment for outlet Mach number M 2is = 0.5 are presented in Fig. 11.The comparison between the CFD and experimental data revealed a very good agreement on upper side of the blade.There is slight yet acceptable difference between CFD and experiment on the lower side of the blade.

CFD results
The results are presented in terms of isentropic Mach number where p is static pressure, p tot is the total pressure at the inlet and = C p ∕C v = 1.4 is the specific heat ratio.Isentropic Mach numbers are evaluated in three ( 12) planes perpendicular to the blade cascade.Sections P1 a P3 are 10 mm from the side walls and section P2 is in the mid-span, see Fig. 12.The local incidence angles for the elastic blade in these planes are P3 = ± 2.37 • , P2 = ± 3.16 • and P1 = ± 3.65 • .The evolution of incidence angle along blade span of elastic blade is shown in Fig. 13.The incidence angle for the rigid blade is caused only by rotation and is constant, P1 = P2 = P3 = ± 1 • .

Positive incidence angle results
Figure 14 shows comparison of the flow fields between POS-EL and POS-RIG cases ( = 1 • ) in plane P2 for outlet isentropic Mach number M 2is = 0.5.
Elastic configuration with higher incidence angle of the middle blade induces faster airflow in the higher interblade channel and slower airflow in the In both cases POS-EL and POS-RIG, the flow is mostly subsonic.Only, in the POS-EL case, the higher incidence angle of the middle blade causes small supersonic area shortly downstream of the leading edge of fourth blade of cascade.
The difference between the elastic and rigid blade in positive incidence angle for outlet Mach number M 2is = 0.7 is significantly larger.The results in the form of isentropic Mach number are depicted in Fig. 15, sonic lines are highlighted in black.In the POS-RIG configuration, there is a large supersonic area in the interblade channel between blade 1 and 2, compared to smaller supersonic area in the case of elastic blade.Opposite phenomenon is visible in the interblade channel between blade 3 and 4.There is larger supersonic area in the POS-EL case due to the higher incidence angle of blade 3 and thus narrower interblade channel.

Negative incidence angle results
Flow fields in plane P2 for cases NEG-EL and NEG-RIG for outlet isentropic Mach number M 2is = 0.5 are depicted in Fig. 16.Higher incidence angle in the NEG-EL case induces faster airflow in the lower interblade channel.The difference of airflow velocity in the lower and higher interblade channel of the middle blade is more significant than in the cases of positive incidence angle (POS-EL, POS-RIG).
In both NEG-EL and NEG-RIG cases there are supersonic areas near the leading edge of the middle blade.In the elastic case the supersonic area is slightly larger.This is in contrast to the POS-EL case where supersonic area was present only on the leading edge of the fourth blade.
Figure 17 shows comparison of the flow fields between NEG-EL and NEG-RIG cases ( = −1 • ) in plane P2 for outlet Mach number M 2is = 0.7 .There is only very little difference in the field of the isentropic Mach number between elastic and rigid blade.
Profiles of the isentropic Mach number along the blade chord are plotted in Fig. 18.For outlet Mach number M 2is = 0.5 , there is considerable difference between the elastic and rigid blade in all planes P1, P2 and P3 due to the higher local incidence angle in the elastic case.The changes of the isentropic Mach number profiles in planes P1, P2 and P3 can be hardly observed.Bigger difference is visible in the case of negative incidence angle (compared to the positive incidence angle) in the isentropic Mach number behind the leading edge between rigid and elastic blade.
For outlet Mach number M 2is = 0.7 and positive incidence angle, the insentropic Mach number profiles are considerably different for the EL and RIG case (see Fig. 19).The same was also clearly visible in Fig. 15, which shows change of locations of supersonic areas between elastic and rigid blade.In case of negative incidence angle and outlet Mach number M 2is = 0.7 , the isentropic Mach number profiles of rigid and elastic case are the most similar.Clearly the

Aerodynamic moments
Total force and moment acting on the middle blade was computed from the flow variables provided by the CFD simulations.Inertial moments were evaluated by structural transient simulations for the elastic case.In case of rigid blade, the inertial moments were estimated by calculation for a flat plate, see Tab. 6.
The aerodynamic moments for POS-RIG, NEG-RIG and NEG-EL are 2-3 times higher for outlet isentropic Mach number M 2is = 0.7 than for M 2is = 0.5 .The faster flow M 2is = 0.7 led to approx.two times lower aerodynamic moment than for M 2is = 0.5 in the POS-EL case.The highest aerodynamic moment 10.26 Nm was observed in NEG-EL configuration.Nevertheless, for the frequency of oscillation f = 200 Hz , the aerodynamic moments are at least ten times lower than inertial moments due to the high frequency oscillation.Inertial moments increase quadratically with the oscillating frequency.The situation is more complicated for the aerodynamic moments.For subsonic flow, they increase with flow velocity and incidence angle.But the value of aerodynamic moment is non-trivial for transonic and supersonic flow regimes.
This difference between aerodynamic and inertial moments justifies the uncoupled approach of flow field computations along the "frozen" blade.Performing full fluid-structure simulations in this case would not help to obtain more reliable data.Even with this simplified model it can be assumed that the inertial forces are more severe from the point of stress induced in the middle blade.

Conclusion
A new blade cascade for controlled flutter testing at transonic speeds was developed.The middle blade in the cascade is exposed to high-frequency oscillations to simulate blade flutter which induces high inertial forces in the blade.For the measurements it is important to know how the elastic deformation of the blade due to inertial forces affects flow field along the blade span -if the flow along the span remains two dimensional.Evaluation of the influence of elastic deformation on the flow field in the blade cascade for controlled flutter measurement was performed by combination of measurements and numerical simulations.
Free experimental modal analysis was performed to estimate damping of the structure.Fixed computational modal analysis provided eigenfrequencies of the blade to avoid resonance during measurements.Free experimental and computational modal analysis was compared for validation which gave good agreement.The deformed shape of the blade was obtained by structural transient simulation at peak positive and negative rotation angle = ± 1 • after stabilization of the initial transient oscillations.The flow fields in the blade cascade for outlet isentropic Mach number M 2is = 0.5 and M 2is = 0.7 were computed for four configurations -positive and negative incidence angle for rigid (idealized) and elastic (real) blade.The flow fields in the elastic configurations were affected by higher incidence angles, but surprisingly only very weakly.Significant difference between elastic and rigid blade was observed only for positive incidence angle and outlet isentropic Mach number M 2is = 0.7 .For the given Mach numbers and given cases of either elastic or rigid blade, the flow fields remain nearly two-dimensional along the blade span.
The structural and aerodynamic moments acting on the blade were computed.Moments induced by inertial forces due to high-frequency oscillation are by an order of magnitude larger than moments caused by the fluid flow.In this case the stresses in the oscillating blade are mainly induced by inertial forces.
Experimental data for a lower isentropic Mach number M 2is = 0.5 was available and CFD model was validated, demonstrating good agreement, especially on the upper side of the blade.
The loading conditions of blades in the experimental setup for forced flutter testing were determined, providing valuable insights for blade design.The computed structural and aerodynamic moments clearly demonstrate that inertial forces due to highfrequency oscillation induce significantly larger moments than those caused by fluid flow for analysed flow regimes.This highlights the importance of considering inertial forces in blade design to ensure the blade's ability to withstand high-frequency oscillations and accurately simulate real-world conditions.These conclusions will guide the development of robust test blades that can cover a wider frequency range, effectively simulating the flow-induced vibrations experienced by real turbomachine blades.In particular, exploring the use of fiber composite materials in blade design offers the potential for improved strength-to-weight ratio and enhanced damping properties, further enhancing the test blade performance and accuracy.Technical Library in Prague.

Declaration
Confict of interest The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Fig. 1
Fig. 1 New test facility for forced flutter research

Fig. 6 Fig. 7 Fig. 8
Fig. 6 Area for evaluation of average local von Mises stress (left).Average local von Mises stress as a function of mesh size (right)

T n = 0 .
Two flow regimes were investigated.In order to reach the outlet isentropic Mach number M 2isA = 0.5 and M 2isB = 0.7 , the static pressure at the outlet boundary condition is set to p 2A = 84, 302 Pa and p 2B = 72, 092 Pa , respectively (see Tab. 5).The back- flow boundary conditions consist of: turbulent intensity T ub = 5% , turbulent viscosity ratio ( T ) b = 10 and total temperature T totb = 300 K .For the velocity, zero gradi- ent is enforced by extrapolating the velocity values from the interior nodes.On the top, bottom and side walls and blades, the no-slip boundary condition u = 0 is specified, further pressure gradient p n = 0 and heat flux T n = 0.

Fig. 10
Fig.10 Geometry and mesh of the 3D computational domain and detail of the mesh refinement near the blade

Fig. 13 Fig. 14
Fig.13 Evolution of the incidence angle along the blade span due to elastic deformation (see also Fig.8)

Fig. 15
Fig. 15 Isentropic Mach number distribution in the elastic (top) and rigid (bottom) case in section cut P2 in rotation = 1 for M 2is = 0.7 .Sonic lines highlighted in black

Table 1
Test rigs for controlled flutter measurements

Table 2
Results of computational free and fixed modal analysis

Table 3
Results of experimental and computational free modal analysis: first torsional (second overall) eigenfrequency and modal damping ratio

Table 4
Computed configurations of the blade

Table 5
Boundary conditions at the outlet