• International Journal of Fluid Machinery and Systems
• /
• v.4 no.1
• /
• pp.199-208
• /
• 2011

This work presents numerical results, using OpenFOAM, of the flow in the swirl flow generator test rig developed at Politehnica University of Timisoara, Romania. The work shows results computed by solving the unsteady Reynolds Averaged Navier Stokes equations. The unsteady method couples the rotating and stationary parts using a sliding grid interface based on a GGI formulation. Turbulence is modeled using the standard k-${\varepsilon}$ model, and block structured wall function ICEM-Hexa meshes are used. The numerical results are validated against experimental LDV results, and against design velocity profiles. The investigation shows that OpenFOAM gives results that are comparable to the experimental and design profiles. The unsteady pressure fluctuations at four different positions in the draft tube is recorded. A Fourier analysis of the numerical results is compared whit that of the experimental values. The amplitude and frequency predicted by the numerical simulation are comparable to those given by the experimental results, though slightly over estimated.

• Journal of Hydrospace Technology
• /
• v.1 no.1
• /
• pp.29-40
• /
• 1995

This paper describes the procedure to predict steady and unsteady performances of a contrarotating propeller(CRP) by a mixed formulation of the boundary value problem(BVP) far the flow around a CRP. The blade BVP is treated by a classical vortex lattice method, whereas the hub BVP is solved by a potential-based panel method. Blades and trailing wakes are represented by a vortex and/or source lattice system, and hubs are represented by normal dipole and source distributions. Both forward and aft propellers are solved simultaneously, thus treating the interaction effect without iteration. The unsteady performance is computed directly in time domain. The new numerical procedure requires a large amount of storage and computing time, which is however no longer a limit in a modern computer system. Sample computations show that the steady performance compares very well with the experiments. The predicted unsteady behavior shows that the dominant harmonics of the total forces are multiples of not only the number of blades of the forward and aft propellers but also the product of both blade numbers. The magnitude of the latter harmonics, present also in uniform oncoming flow, may reach abort 50% of the mean torque for the aft propeller, which in turn may cause a serious vibration problem in the complicated contrarotating shafting system.

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

A high-order accurate Euler flow solver based on a discontinuous Galerkin finite-element method has been developed for the numerical simulations of blade-vortex interaction phenomena on unstructured meshes. A free vortex in freestream was investigated to assess the vortex-preserving property and the accuracy of the present flow solver. Blade-vortex interaction problems in subsonic and transonic freestreams were simulated by adopting a multi-level solution-adaptive dynamic mesh refinement/coarsening technique. The results were compared with those of other numerical and experimental methods. It was shown that the present discontinuous Galerkin flow solver can preserve the vortex structure for significantly longer vortex convection time and can accurately capture the complex unsteady blade-vortex interaction flows, including generation and propagation of acoustic waves.

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

A high-order accurate Euler flow solver based on a discontinuous Galerkin finite-element method has been developed for the numerical simulations of blade-vortex interaction phenomena on unstructured meshes. A free vortex in freestream was investigated to assess the vortex-preserving property and the accuracy of the present flow solver. Blade-vortex interaction problems in subsonic and transonic freestreams were simulated by adopting a multi-level solution-adaptive dynamic mesh refinement/coarsening technique. The results were compared with those of other numerical and experimental methods. It was shown that the present discontinuous Galerkin flow solver can preserve the vortex structure for significantly longer vortex convection time and can accurately capture the complex unsteady blade-vortex interaction flows, including generation and propagation of acoustic waves.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.10
• /
• pp.1401-1411
• /
• 2003

A two-dimensional direct numerical simulation was performed to investigate the flame structure of CH$_4$$N_2$-air counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman's detailed chemistry were adopted in this simulation. The characteristic vortex and chemical time scales were introduced to quantify and investigate the extinction phenomenon during a flame-vortex interaction. The results showed that fuel- and air-side vortex cause an unsteady extinction. In this case, the flame interacting with a vortex was extinguished at much larger scalar dissipation rate than steady flame. It was also found that the air-side vortex extinguished a flame more rapidly than the fuel-side vortex. Furthermore, it was noted that the degree of unsteady effect experienced by a flame can be investigated by comparing the above two characteristic time scales, and this analysis could give an appropriate reason for the results of the previously reported experiment.

• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.2
• /
• pp.110-120
• /
• 2006

There have been developed many structural and fluid rotorcraft analysis models in rotorcraft community, and also lots of investigations have been conducted to combine these two models. These investigations turn out to be good at predicting the airloads precisely, but they have not taken the blade nonlinear deflection into account. For this reason, the present paper adopts a sophisticated structural model which can describe three-dimensional nonlinear deflection of the blade. And it is combined with two types of aerodynamic model. First one is generalized Greenberg type of finite-time aerodynamic model, which is originally established for a fixed wing, but later modified to be suitable for coupled flap-lag-torsional aeroelastic analysis of the rotor blade. Second aerodynamic model is based on the unsteady source-doublet panel method coupled with a free wake model. The advantages of the present method are capabilities to consider thickness of the blade and more precise wake effects. Transient responses of the airloads and structural deflections in time domain are mainly analyzed in this paper.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2005.11a
• /
• pp.86-89
• /
• 2005

A computational work has been performed to investigate the aerodynamics of a projectile which is launched from the two-stage light gas gun. A moving coordinate method for a multi-domain technique is employed to simulate unsteady projectile flows with a moving boundary. The effect of a virtual mass is added to the axisymmetric unsteady Euler equation system. The computed results reasonably capture the major flow characteristics which we generated in launching the projectile supersonically, such as the interaction between the shock wave and the blast wave, the interaction between the vortical flow and the barrel shock, and the steady under-expanded jet. The present computational results properly predict the velocity, acceleration, and drag histories of the projectile.

• Wind and Structures
• /
• v.13 no.2
• /
• pp.169-189
• /
• 2010

This paper describes the use of coupled Computational Fluid Dynamics (CFD) and Rigid Body Dynamics (RBD) in modelling the aerodynamic behaviour of wind-borne plate type objects. Unsteady 2D and 3D Reynolds Averaged Navier-Stokes (RANS) CFD models are used to simulate the unsteady and non-uniform flow field surrounding static, forced rotating, auto-rotating and free-flying plates. The auto-rotation phenomenon itself is strongly influenced by vortex shedding, and the realisable k-epsilon turbulence modelling approach is used, with a second order implicit time advancement scheme and equal or higher order advection schemes for the flow variables. Sequentially coupling the CFD code with a RBD solver allows a more detailed modelling of the Fluid-Structure Interaction (FSI) behaviour of the plate and how this influences plate motion. The results are compared against wind tunnel experiments on auto-rotating plates and an existing 3D analytical model.

• International Journal of Aeronautical and Space Sciences
• /
• v.10 no.2
• /
• pp.95-105
• /
• 2009

The flow fields around the HARTII rotor were numerically investigated using a viscous flow solver on adaptive unstructured meshes. An overset mesh and a deforming mesh technique were used to handle the blade motion including blade deflection, which was obtain from the HARTII experimental data. A solution-adaptive mesh refinement technique was also used to capture the rotor wake effectively. Comparison of the sectional normal force and pitching moment at 87% radial station between the two cases, with and without the blade deflection, showed that the blade loading is significantly affected by blade torsion. It was found that as the mesh was refined, the strength of tip vortex is better preserved, and the magnitude of high frequency blade loading, caused by blade-vortex interaction (BVI), is further magnified. It was also found that a proper time step size, which corresponds to the cell size, should be used to predict unsteady solutions accurately. In general, the numerical results in terms of the unsteady blade loading and the rotor wake show good agreement with the experimental data.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2009.05a
• /
• pp.79-82
• /
• 2009

Unsteady numerical analysis of an entire supersonic propulsion system from intake to nozzle was performed to study dynamic interaction between terminal shock in the intake and flame in the combustor. Both acceleration and cruise flight-modes were considered. Acoustic mode of the entire engine for both flight-modes were investigated by detail analysis of pressure fluctuation at each location of engine.