scholarly journals Comparison of Two Numerical Algorithms for Computing the Effects of Mistuning of Fan Flutter

Author(s):  
L. Salles ◽  
M. Vahdati

The aim of this paper is to study the effects of mistuning on fan flutter and to compare the prediction of two numerical models of different fidelity. The high fidelity model used here is a three-dimensional, whole assembly, time-accurate, viscous, finite-volume compressible flow solver. The Code used for this purpose is AU3D, written in Imperial College and validated for flutter computations over many years. To the best knowledge of authors, this is the first time such computations have been attempted. This is due to the fact that, such non-linear aeroelastic computations with mistuning require large amount of CPU time and cannot be performed routinely and consequently, faster (low fidelity) models are required for this task. Therefore, the second model used here is the aeroelastic fundamental mistuning model (FMM) and it based on an eigenvalue analysis of the linearized modal aeroelastic system with the aerodynamic matrix calculated from the aerodynamic influence coefficients. The influence coefficients required for this algorithm are obtained from the time domain non-linear Code by shaking one blade in the datum (tuned) frequency and mode. Once the influence coefficients have been obtained, the computations of aero damping require minimal amount of CPU time and many different mistuning patterns can be studied. The objectives of this work are to: 1. Compare the results between the two models and establish the capabilities/limitations of aeroelastic FMM, 2. Check if the introduction of mistuning would bring the experimental and computed flutter boundaries closer, 3. Establish a relationship between mistuning and damping. A rig wide-chord fan blade, typical of modern civil designs, was used as the benchmark geometry for this study. All the flutter analyses carried out in this paper are with frequency mistuning, but the possible consequences of mistuned mode shapes are briefly discussed at the end of this paper. Only the first family of modes (1F, first flap) is considered in this work. For the frequency mistuning analysis, the 1F frequency is varied around the annulus but the 1F mode shapes remain the same for all the blades. For the mode shape mistuning computations, an FE analysis of the whole assembly different mass blades is performed. The results of this work clearly show the importance of mistuning on flutter. It also demonstrates that when using rig test data for aeroelastic validation of CFD codes, the amount mistuning present must be known. Finally, it should be noted that the aim of this paper is the study of mistuning and not steady/unsteady validation of a CFD code and therefore minimal aerodynamic data are presented.

Author(s):  
Dimitri Franz ◽  
Loic Salles ◽  
Sina Stapelfeldt

In this paper the effects of mistuning on the flutter stability of a turbine blade are analysed. Two types of mistuning are considered, frequency mistuning and aerodynamic mistuning. The study concentrates on the the first family of modes (1F, first flap) as the blade fluttered in this mode during test. For the frequency mistuning analysis, the 1F frequency is varied around the annulus but the 1F mode shapes remain the same for all the blades. The mistuning analyses are performed by using a reduced order model (ROM) based on an eigenvalue analysis of the linearized modal aeroelastic system with the aerodynamic matrix calculated from the aerodynamic influence coefficients. The influence coefficients required for this algorithm are obtained from a three-dimensional, non-linear aeroelastic solver (AU3D) by shaking one blade in the datum (tuned) frequency and mode and recording aerodynamic forces on the other blades in the assembly. After the ROM is validated against the non-linear method for the tuned case, it is used for the mistuning and mis-staggering study as time-domain computations of such cases are very time consuming. The results of this paper indicate that, frequency mistuning is always stabilizing but aerodynamic mistuning can be destabilizing under certain conditions. Moreover, it is shown that the effect of frequency mistuning is much higher than the one of aerodynamic asymmetries and that structural coupling limits the effects of mistuning.


Author(s):  
Yoshiyuki Inoue ◽  
Md. Kamruzzaman

The LNG-FPSO concept is receiving much attention in recent years, due to its active usage to exploit oil and gas resources. The FPSO offloads LNG to an LNG carrier that is located close to the FPSO, and during this transfer process two large vessels are in close proximity to each other for daylong periods of time. Due to the presence of neighboring vessel, the motion response of both the vessels will be affected significantly. Hydrodynamic interactions related to wave effects may result in unfavorable responses or the risk of collisions in a multi-body floating system. Not only the motion behavior but also the second order drift forces are influenced by the neighboring structures due to interactions of the waves among the structures. A study is made on the time domain analysis to assess the behavior and the operational capability of the FPSO system moored in the sea having an LNG carrier alongside under environmental conditions such as waves, wind and currents. This paper presents an analysis tool to predict the dynamic motion response and non-linear connecting and mooring forces on a parallel-connected LNG-FPSO system due to non-linear exciting forces of wave, wind and current. Simulation for the mooring performance is also investigated. The three-dimensional source-sink technique has been applied to obtain the radiation forces and the transfer function of wave exciting forces on floating multi-bodies. The hydrodynamic interaction effect between the FPSO and the LNG carrier is included to calculate the hydrodynamic forces. For the simulation of a random sea and also for the generation of time depended wind velocity, a fully probabilistic simulation technique has been applied. Wind and current loads are estimated according to OCIMF. The effects of variations in wave, wind and current loads and direction on the slowly varying oscillations of the LNG and FPSO are also investigated in this paper. Finally, some conclusions are drawn based on the numerical results obtained from the present time domain simulations.


Author(s):  
Dora Foti ◽  
Mariella Diaferio ◽  
Nicola Ivan Giannoccaro ◽  
Salvador Ivorra

In the present chapter the theoretical basis of different methods developed for the calibration of FEMs are discussed. In general, Model Updating techniques are based on the use of appropriate functions that iteratively update selected physical properties (characteristics of the materials, stiffness of a link, etc.). In this way the correlation between the simulated response and the target value could improve if compared to an initial value. The FE model thus obtained can be used for a detailed structural analysis with a great confidence. The technique described in the first part of the chapter is applied to the evaluation of the structural properties of the tower of the Provincial Administration Building in Bari (Italy).The final purpose is to predict the performance of the tower to different combinations of static and dynamic loads, i.e. earthquakes or other induced vibrations. Ambient vibration tests have been performed on the above mentioned tower with the aim of determining its dynamic response and developing a procedure for modeling this building (Foti et al., 2012a). The Operation Modal Analysis (OMA) has been carried out both in the frequency domain and in the time domain to extract the dominant frequencies and mode shapes of the tower.


2000 ◽  
Vol 123 (1) ◽  
pp. 17-26 ◽  
Author(s):  
Shuangbiao Liu ◽  
Qian Wang

A necessary step in understanding failure problems of tribological elements is to investigate the contact performance of rough surfaces subjected to frictional heating. It is essential that the interfacial variables are obtained through solving the interactive thermomechanical contact problem. This paper studies the three dimensional thermomechanical contact of non-conforming rough surfaces, the model of which includes the normal surface displacements caused by the contact pressure, frictional shear, and frictional heating. Influence coefficients and frequency response functions for elastic and thermoelastic displacements, as well as those for temperature rises, are investigated for model construction. In order to develop an accurate and efficient solver, the numerical algorithms with the discrete convolution and fast Fourier transform techniques and the single-loop conjugated gradient method are used. The model modules are numerically verified and the thermomechanical performance of the rough surfaces in a point contact is studied.


Author(s):  
Dora Foti ◽  
Mariella Diaferio ◽  
Nicola Ivan Giannoccaro ◽  
Salvador Ivorra

In the present chapter the theoretical basis of different methods developed for the calibration of FEMs are discussed. In general, Model Updating techniques are based on the use of appropriate functions that iteratively update selected physical properties (characteristics of the materials, stiffness of a link, etc.). In this way the correlation between the simulated response and the target value could improve if compared to an initial value. The FE model thus obtained can be used for a detailed structural analysis with a great confidence. The technique described in the first part of the chapter is applied to the evaluation of the structural properties of the tower of the Provincial Administration Building in Bari (Italy).The final purpose is to predict the performance of the tower to different combinations of static and dynamic loads, i.e. earthquakes or other induced vibrations. Ambient vibration tests have been performed on the above mentioned tower with the aim of determining its dynamic response and developing a procedure for modeling this building (Foti et al., 2012a). The Operation Modal Analysis (OMA) has been carried out both in the frequency domain and in the time domain to extract the dominant frequencies and mode shapes of the tower.


Author(s):  
Y-S Wu ◽  
W-C Cui

The study of hydroelasticity of ships first gained momentum in the late 1970s with the work of Bishop and Price, who established the two-dimensional (2D) hydroelasticity theory of ships. The concept and basic principle presented in their work to embody the structure and the surrounding fluid as a coupled entirety was further employed and extended in the creation of the general linear three-dimensional (3D) theory of hydroelasticity for an arbitrary shaped flexible marine structure travelling with a forward speed in a seaway in the middle of 1980s (Wu, 1984; Price and Wu, 1985a; Bishop et al., 1986). Since then, great progress has been achieved in the development and application of 3D hydroelasticity theories. These include the more rigorous methods of frequency-domain linear analysis accounting for the forward speed effect and the steady flow effect, the time-domain linear 3D theory, the non-linear 3D theory and the numerical methods for a floating structure travelling in rough seas with large motions, experimental techniques of 3D flexible ship models, the hydroelasticity-based design and safety assessment, etc. This paper presents an overview of these developments and achievements of linear and non-linear 3D hydroelasticity theories of ships, and the corresponding numerical and experimental techniques.


Author(s):  
G. Jacobs ◽  
F. Theunissen

In order to understand how the algorithms underlying neural computation are implemented within any neural system, it is necessary to understand details of the anatomy, physiology and global organization of the neurons from which the system is constructed. Information is represented in neural systems by patterns of activity that vary in both their spatial extent and in the time domain. One of the great challenges to microscopists is to devise methods for imaging these patterns of activity and to correlate them with the underlying neuroanatomy and physiology. We have addressed this problem by using a combination of three dimensional reconstruction techniques, quantitative analysis and computer visualization techniques to build a probabilistic atlas of a neural map in an insect sensory system. The principal goal of this study was to derive a quantitative representation of the map, based on a uniform sample of afferents that was of sufficient size to allow statistically meaningful analyses of the relationships between structure and function.


Algorithms ◽  
2021 ◽  
Vol 14 (5) ◽  
pp. 129
Author(s):  
Yuan Li ◽  
Ni Zhang ◽  
Yuejiao Gong ◽  
Wentao Mao ◽  
Shiguang Zhang

Compared with continuous elements, discontinuous elements advance in processing the discontinuity of physical variables at corner points and discretized models with complex boundaries. However, the computational accuracy of discontinuous elements is sensitive to the positions of element nodes. To reduce the side effect of the node position on the results, this paper proposes employing partially discontinuous elements to compute the time-domain boundary integral equation of 3D elastodynamics. Using the partially discontinuous element, the nodes located at the corner points will be shrunk into the element, whereas the nodes at the non-corner points remain unchanged. As such, a discrete model that is continuous on surfaces and discontinuous between adjacent surfaces can be generated. First, we present a numerical integration scheme of the partially discontinuous element. For the singular integral, an improved element subdivision method is proposed to reduce the side effect of the time step on the integral accuracy. Then, the effectiveness of the proposed method is verified by two numerical examples. Meanwhile, we study the influence of the positions of the nodes on the stability and accuracy of the computation results by cases. Finally, the recommended value range of the inward shrink ratio of the element nodes is provided.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Souma Jinno ◽  
Shuji Kitora ◽  
Hiroshi Toki ◽  
Masayuki Abe

AbstractWe formulate a numerical method on the transmission and radiation theory of three-dimensional conductors starting from the Maxwell equations in the time domain. We include the delay effect in the integral equations for the scalar and vector potentials rigorously, which is vital to obtain numerically stable solutions for transmission and radiation phenomena in conductors. We provide a formalism to connect the conductors to any passive lumped-parameter circuits. We show one example of numerical calculations, demonstrating that the new formalism provides stable solutions to the transmission and radiation phenomena.


Sign in / Sign up

Export Citation Format

Share Document