• Advances in aircraft and spacecraft science
• /
• v.4 no.2
• /
• pp.203-218
• /
• 2017

This paper focuses on the computational study of nonlinear effects of unsteady aerodynamics for non-classical aileron buzz. It aims at a comprehensive investigation of the aileron buzz phenomenon under varying flow parameters using the describing function technique with multiple inputs. The limit cycle oscillatory behavior of an asymmetrical airfoil is studied initially using a CFD-based numerical model and direct time marching. Sharp increases in limit cycle amplitude for varying Mach numbers and angles of attack are investigated. An aerodynamic describing function is developed in order to estimate the variation of limit cycle amplitude and frequency with Mach number and angle of attack directly, without time marching. The describing function results are compared to the amplitudes and frequencies predicted by the CFD calculations for validation purposes. Furthermore, a limited sensitivity analysis is presented to demonstrate the potential of the approach for aeroelastic design.

• International Journal of Fluid Machinery and Systems
• /
• v.6 no.2
• /
• pp.87-93
• /
• 2013

In order to study the effect of the fluid-structure interaction (FSI) on the simulation results, the external characteristics and internal flow features of a diffuser pump were analyzed with a two-way flow solid coupling method. And the static and dynamic structure analysis of the blade was also caculated with the FEA method. The steady flow field is based on Reynolds Averaged N-S equations with standard $k-{\varepsilon}$ turbulent model, the unsteady flow field is based on the large eddy simulation, and the structure response is based on elastic transient structural dynamic equation. The results showed that the effect of FSI on the head prediction based on CFD really exists. At the same radius, the van mises stress on the nodes closed shroud and hub was larger than other nodes. A large deformation region existed near inlet side at the middle of blades. The strength of impeller satisfied the strength requirement with static stress analysis based on the fourth strength theory. The dynamic stress varied periodically with the impeller rotating. It was also found that the fundamental frequency of the dynamic stress is the rotating frequency and its harmonic frequency. The frequency of maximum stress amplitude at node 1626 was 7 times of the rotating frequency. The frequency of maximum stress amplitude at node 2328 was 14 times of the rotating frequency. No matter strength failure or fatigue failure, the root of blades near shroud is the key region to analyse.

• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.361-379
• /
• 2012

This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.42 no.3
• /
• pp.191-198
• /
• 2014

During multi-physics or multi-phase simulations accompanying fluid- structure- thermal interaction, data transfer problems always arise along non- matching interfaces caused by different computational meshes for each physical domain. Common- refinement scheme, among many available methods, is attractive since it is known to yield conservative and accurate data transfer for non- matching interface cases. This is particularly important in simulating compressible unsteady fluid- structure- thermal interaction inside solid propellant rockets, where grid size along solid- fluid interfaces is substantially different. From this perspective, we examine performances of common- refinement- based data transfer scheme between structured quadrilateral (structure part) and unstructured triangular (fluid part) meshes by comparing computed results with other data transfer methods.

• Nuclear Engineering and Technology
• /
• v.46 no.5
• /
• pp.633-640
• /
• 2014

In this research, the surface roughness affecting the pressure drop in a pipe used as the steam generator of a PWR was studied. Based on the CFD (Computational Fluid Dynamics) technique using a commercial code named ANSYS-FLUENT, a straight pipe was modeled to obtain the Darcy frictional coefficient, changed with a range of various surface roughness ratios as well as Reynolds numbers. The result is validated by the comparison with a Moody chart to set the appropriate size of grids at the wall for the correct consideration of surface roughness. The pressure drop in a full-scale U-shaped pipe is measured with the same code, correlated with the surface roughness ratio. In the next stage, we studied a reduced scale model of a U-shaped heat pipe with experiment and analysis of the investigation into fluid-structure interaction (FSI). The material of the pipe was cut from the real heat pipe of a material named Inconel 690 alloy, now used in steam generators. The accelerations at the fixed stations on the outer surface of the pipe model are measured in the series of time history, and Fourier transformed to the frequency domain. The natural frequency of three leading modes were traced from the FFT data, and compared with the result of a numerical analysis for unsteady, incompressible flow. The corresponding mode shapes and maximum displacement are obtained numerically from the FSI simulation with the coupling of the commercial codes, ANSYS-FLUENT and TRANSIENT_STRUCTURAL. The primary frequencies for the model system consist of three parts: structural vibration, BPF(blade pass frequency) of pump, and fluid-structure interaction.

• International Journal of Fluid Machinery and Systems
• /
• v.3 no.4
• /
• pp.315-323
• /
• 2010

This work investigates the influence of water compressibility on pressure pulsations induced by rotor-stator interaction (RSI) in hydraulic machinery, using the commercial CFD solver ANSYS-CFX. A pipe flow example with harmonic velocity excitation at the inlet plane is simulated using different grid densities and time step sizes. Results are compared with a validated code for hydraulic networks (SIMSEN). Subsequently, the solution procedure is applied to a simplified 2.5-dimensional pump-turbine configuration in prototype with different speeds of sound as well as in model scale with an adapted speed of sound. Pressure fluctuations are compared with numerical and experimental data based on prototype scale. The good agreement indicates that the scaling of acoustic effects with an adapted speed of sound works well. With respect to pressure fluctuation amplitudes along the centerline of runner channels, incompressible solutions exhibit a linear decrease while compressible solutions exhibit sinusoidal distributions with maximum values at half the channel length, coinciding with analytical solutions of one-dimensional acoustics. Furthermore, in compressible simulation the amplification of pressure fluctuations is observed from the inlet of stay vane channels to the spiral case wall. Finally, the procedure is applied to a three-dimensional pump configuration in model scale with adapted speed of sound. Normalized Pressure fluctuations are compared with results from prototype measurements. Compared to incompressible computations, compressible simulations provide similar pressure fluctuations in vaneless space, but pressure fluctuations in spiral case and penstock may be much higher.

• The KSFM Journal of Fluid Machinery
• /
• v.2 no.3 s.4
• /
• pp.45-51
• /
• 1999

Centrifugal fans are widely used and the noise generated by these machines causes one of the most serious problems. In general, the centrifugal fan noise is often dominated by tones at BPF(blade passage frequency) and its higher harmonics. This is a consequence of the strong interaction between the flow discharged from the impeller and the cutoff in the casing. However, only a few researches have been carried out on predicting the noise due to the difficulty in obtaining detailed information about the flow field and casing effects on noise radiation. The objective of this study is to develop a prediction method for the unsteady flow field and the acoustic pressure field of a centrifugal fan and to calculate the effects of rotating velocity, flow rate, cut-off distance and the number of blades and its effects on the noise of the fan. We assume that the impeller rotates with a constant angular velocity and the flow field around the impeller is incompressible and inviscid. So, a discrete vortex method (DVM) is used to model the centrifugal fan and to calculate the flow field. The force of each element on the blade is calculated with the unsteady Bernoulli equation. Lowson's method is used to predict the acoustic source. The cut-off distance is the most important factor effecting the noise generation. Acoustic pressure is proportional to 2.8, which shows the same scaling index as the experimental result. In this paper, the cut-off distance is found to be the dominant parameter offecting the acoustic pressure.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.7
• /
• pp.647-654
• /
• 2010

In this study, an efficient ornithopter aerodynamic model, which is applicable to ornithopter wing design considering fluid-structure interaction or ornithopter flight dynamics and control simulation, was proposed and experimentally validated through the wind tunnel experiments. Due to the ornithopter aerodynamics governed by unsteady low Reynolds number flow, an experimental device was specially designed and developed. A part of the experimental device, 2-axis loadcell, was situated in the non-inertial frame; the dynamic calibration method was established to compensate the inertial load for pure aerodynamic load measurements. The characteristics of proposed aerodynamic model were compared with the experimental data in terms of mean and root-mean-square values of lift and drag coefficients with respect to the flow speed, flapping frequency, and fixed angle of attack.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.367-373
• /
• 2008

Spinning detonations propagating in a circular tube were numerically investigated with a one-step irreversible reaction model governed by Arrhenius kinetics. Activation energy is used as parameter as 10, 20, 27 and 35, and the specific heat ratio and the heat release are fixed as 1.2 and 50. The time evolution of the simulation results was utilized to reveal the propagation mechanism of single-headed spinning detonation. The track angle of soot record on the tube wall was numerically reproduced with various levels of activation energy, and the simulated unique angle was the same as that of the previous reports. The maximum pressure histories of the shock front on the tube wall showed stable pitch at Ea=10, periodical unstable pitch at Ea=20 and 27 and unstable pitch consisting of stable, periodical unstable and weak modes at Ea=35, respectively. In the weak mode, there is no Mach leg on the shock front, where the pressure level is much lower than the other modes. The shock front shapes and the pressure profiles on the tube wall clarified the mechanisms of these stable and unstable modes. In the stable pitch at Ea=10, the maximum pressure history on the tube wall remained nearly constant, and the steady single Mach leg on the shock front rotated at a constant speed. The high and low frequency pressure oscillations appeared in the periodical unstable pitch at Ea=20 and 27 of the maximum pressure history. The high frequency was one cycle of a self-induced oscillation by generation and decay in complex Mach interaction due to the variation in intensity of the transverse wave behind the shock front. Eventually, sequential high frequency oscillations formed the low frequency behavior because the frequency behavior was not always the same for each cycle. In unstable pitch at Ea=35, there are stable, periodical unstable and weak modes in one cycle of the low frequency oscillation in the maximum pressure history, and the pressure amplitude of low frequency was much larger than the others. The pressure peak appeared after weak mode, and the stable, periodical unstable and weak modes were sequentially observed with pressure decay. A series of simulations of spinning detonations clarified that the unsteady mechanism behind the shock front depending on the activation energy.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.367-373
• /
• 2008

Spinning detonations propagating in a circular tube were numerically investigated with a one-step irreversible reaction model governed by Arrhenius kinetics. Activation energy is used as parameter as 10, 20, 27 and 35, and the specific heat ratio and the heat release are fixed as 1.2 and 50. The time evolution of the simulation results was utilized to reveal the propagation mechanism of single-headed spinning detonation. The track angle of soot record on the tube wall was numerically reproduced with various levels of activation energy, and the simulated unique angle was the same as that of the previous reports. The maximum pressure histories of the shock front on the tube wall showed stable pitch at Ea=10, periodical unstable pitch at Ea=20 and 27 and unstable pitch consisting of stable, periodical unstable and weak modes at Ea=35, respectively. In the weak mode, there is no Mach leg on the shock front, where the pressure level is much lower than the other modes. The shock front shapes and the pressure profiles on the tube wall clarified the mechanisms of these stable and unstable modes. In the stable pitch at Ea=10, the maximum pressure history on the tube wall remained nearly constant, and the steady single Mach leg on the shock front rotated at a constant speed. The high and low frequency pressure oscillations appeared in the periodical unstable pitch at Ea=20 and 27 of the maximum pressure history. The high frequency was one cycle of a self-induced oscillation by generation and decay in complex Mach interaction due to the variation in intensity of the transverse wave behind the shock front. Eventually, sequential high frequency oscillations formed the low frequency behavior because the frequency behavior was not always the same for each cycle. In unstable pitch at Ea=35, there are stable, periodical unstable and weak modes in one cycle of the low frequency oscillation in the maximum pressure history, and the pressure amplitude of low frequency was much larger than the others. The pressure peak appeared after weak mode, and the stable, periodical unstable and weak modes were sequentially observed with pressure decay. A series of simulations of spinning detonations clarified that the unsteady mechanism behind the shock front depending on the activation energy.