• Korea-Australia Rheology Journal
• /
• v.18 no.2
• /
• pp.109-118
• /
• 2006

A finite-element method was employed to analyze axisymmetric unsteady motion of a deformable bubble near the wall. In the present study a deformable bubble in a Newtonian medium near the wall was considered. In solving the governing equations a structured mesh generator was used to describe the collapse of highly deformed bubbles with the Arbitrary Lagrangian Eulerian (ALE) method being employed in order to capture the transient bubble boundary effectively. In order to check the accuracy of the present FE analysis we compared the results of our FE solutions with the result of the collapse of spherical bubbles in a large body of fluid in which solutions can be obtained using a 1D FE analysis. It has been found that 1D and 2D bubble deformations are in good agreement for spherically symmetric problems confirming the validity of the numerical code. Non-spherically symmetric problems were also solved for the collapse of bubble located near a plane solid wall. We have shown that a microjet develops at the bubble boundary away from the wall as already observed experimentally. We have discussed the effect of Reynolds number and distance of the bubble center from the wall on the transient collapse pattern of bubble.

• The KSFM Journal of Fluid Machinery
• /
• v.14 no.3
• /
• pp.39-44
• /
• 2011

In this study, performance analyses have been conducted for a 5MW class wind turbine blade model. Advanced computational analysis system based on computational fluid dynamics(CFD) and computational structural dynamics(CSD) has been developed in order to investigate detailed dynamic responsed of wind turbine blade. Reynolds-averaged Navier-Stokes (RANS) equations with K-${\epsilon}$ turbulence model are solved for unsteady flow problems of the rotating turbine blade model. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of the 3D turbine blade for fluid-structure interaction (FSI) problems. Predicted aerodynamic performance considering structural deformation effect of the blade show different results compared to the case of rigid blade model.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2011.04a
• /
• pp.663-670
• /
• 2011

In this study, the design and verification of 6kW class lift-type vertical-axis wind turbine (VAWT) has been conducted using advanced CAE technique based on computational fluid dynamics (CFD), finite element method (FEM), and computational structural dynamics (CSD). Designed aerodynamic performance of the VAWT model is tested using unsteady CFD method. Designed structural safety is also tested through the evaluation of maximum induced stress level and resonance characteristics using FEM and CSD methods. It is importantly shown that the effect of master eccentricity due to rotational inertia needs to be carefully considered to additionally investigate dynamic stress and deformation level of the designed VAWT system.

• 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.

• Structural Engineering and Mechanics
• /
• v.85 no.1
• /
• pp.105-118
• /
• 2023

The present work is an attempt to develop a simple and accurate finite element formulation for the assessment of thermal shock/thermally induced vibrations in pretwisted and tapered functionally graded material thin (FGM) blades obtained from Voigt and local representative volume elements (LRVE) homogenization models, based on neutral surface approach. The neutral surface of the FGM blade does not coincide with its mid-surface. A finite element model (FEM) is developed using first-order shear deformation theory (FSDT) and the FGM turbine blade is modelled according to the shallow shell theory. The top and the bottom layers of the FGM blade are made of pure ceramic and pure metal, respectively and temperature-dependent material properties are functionally graded in the thickness direction, the position of the neutral surface also depends on the temperature. The material properties are estimated according to two different homogenization models viz., Voigt or LRVE. The top layer of the FGM blade is subjected to high temperature and the bottom surface is either thermally insulated or kept at room temperature. The solution of the nonlinear profile of the temperature in the thickness direction is obtained from the Fourier law of heat conduction in the unsteady state. The results obtained from the present FEM are compared with the benchmark examples. Next, the effect of angle of twist, intensity of thermal shock, variable chord and span and volume fraction index on the transient response due to thermal shock obtained from the two homogenization models viz., Voigt and LRVE scheme is investigated. It is shown that there can be a significant difference in the transient response calculated by the two homogenization models for a particular set of material and geometric parameters.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.8
• /
• pp.1067-1074
• /
• 2010

In this study, a loosely coupled fluid-structure interaction (FSI) analysis was conducted for a low-pressure (LP) final-stage rotor blade. Preliminary FSI analyses of a $15^{\circ}$ sweptback wing and a NASA Rotor 37 compressor blade were performed for verifying the boundary conditions. The results were compared with the established literatures for each model. The FSI analysis of the $15^{\circ}$ sweptback wing was carried out under both stable and unstable conditions. The excessive deformation of the wing was observed within 0.05 s under the unstable condition which is higher than the divergence speed of a wing compared with the stable condition. On the basis of the results of a steady-state study, an unsteady state FSI analysis was conducted for a NASA Rotor 37. Different deformations were observed at trailing edge of the blade in the static FSI and dynamic FSI analysis. A 3D FE model of a LP rotor was generated from the span-wise section data. In order to develop a reasonable model, an impact test was performed and compared to the FE model. Using this FE model, the steady-state FSI analysis was performed successfully.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.41 no.5
• /
• pp.361-366
• /
• 2017

As the recent development of computing architecture and application software technology, real world simulation, which is the ultimate destination of computer simulation, is emerging as a practical issue in several research sectors. In this paper, metal plate motion in a square shock tube for small time interval was calculated using a supercomputing-based fluid-structure-combustion multi-physics simulation tool called Illinois Rocstar, developed in a US national R amp; D program at the University of Illinois. Afterwards, the simulation results were compared with those from experiments. The coupled solvers for unsteady compressible fluid dynamics and for structural analysis were based on the finite volume structured grid system and the large deformation linear elastic model, respectively. In addition, a strong correlation between calculation and experiment was shown, probably because of the predictor-corrector time-integration scheme framework. In the future, additional validation studies and code improvements for higher accuracy will be conducted to obtain a reliable open-source software research tool.

• Journal of Korea Water Resources Association
• /
• v.54 no.8
• /
• pp.629-641
• /
• 2021

In this study, a two-dimensional numerical model (Nays2DH) was applied to analyze the process of morphological changes in the river channel bed depending on the changes in the amount of flooding after fully opening the Sejong weir, which was constructed upstream of the Geum River. For this, numerical simulations were performed by assuming the flow conditions, such as a non-uniform flow (NF), unsteady flows (single flood event, SF), and a continuous flood event (CF). Here, in the cases of the SF and CF, the normalized hydrograph was calculated from real flood events, and then the hydrograph was reconfigured by the peak flow discharge according to the scenario, and then it was employed as the flow discharge at the upstream boundary condition. In this study, to quantitatively evaluate the morphological changes, we analyzed the time changes in the bed deformation the bed relief index (BRI), and we compared the aerial photographs of the study area and the numerical simulation results. As simulation results of the NF, when the steady flow discharge increases, the ratio of lower width to depth decreases and the speed of bar migration increases. The BRI initially increases, but the amount of change decreased with time. In addition, when the steady flow discharge increases, the BRI increased. In the case of SF, the speed of bar migration decreased with the change of the flow discharge. In terms of the morphological response to the peak flood discharge, the time lag also indicated. In other words, in the SF, the change of channel bed indicates a phase lag with respect to the hydraulic condition. In the result of numerical simulation of CF, the speed of bar migration depending on the peak flood discharges decreased exponentially despite the repeated flood occurrences. In addition, as in the result of SF, the phase lag indicated, and the speed of bar migration decreased exponentially. The BRI increased with time changes, but the rate of increase in the BRI was modest despite the continuous peak flooding. Through this study, the morphological changes based on the hydrological characteristics of the river were analyzed numerically, and the methodology suggested that a quantitative prediction for the river bed change according to the flow characteristic can be applied to the field.