• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.12
• /
• pp.1153-1161
• /
• 2010

An efficient solution method for harmonic balance techniques with Fourier transform is presented for periodic unsteady flow problems. The present partially-implicit harmonic balance treats the flux terms implicitly and the harmonic source term is solved explicitly. The convergence of the partially Implicit method is much faster than the explicit Runge-Kutta harmonic balance method. The method does not need to compute the additional flux Jacobian matrix from the implicit harmonic source term. Compared with fully implicit harmonic balance method, this partial approach turns out to have good convergence property. Oscillating flows over NACA0012 airfoil are considered to verify the method and to compare with results of explicit R-K(Runge-Kutta) and dual time stepping methods.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.21 no.1
• /
• pp.45-50
• /
• 2013

In this study, lumped-vortex element method and thin airfoil theory were used to analyze aerodynamic characteristics of airfoils with relative motion that had camber lines of NACA $44{\times}{\times}$ airfoil in 2-dimensional unsteady incompressible potential flow. Velocity disturbance due to airfoil was calculated by lumped-vortex element model and force distribution on airfoil by unsteady Bernoulli's equation. Variables in relative motion were considered the period p, the amplitude of flapping $A_f$ and pitching $A_p$, and the phase difference between flapping and pitching ${\phi}_p$ and the angle of attack ${\alpha}$. Due to movement of an airfoil, dag was induced in 2-dimensional unsteady incompressible potential flow. The numerical results show that the aerodynamic characteristics of the airfoil with flapping and pitching at the same time are illustrated. Especially the mean lift coefficient became smaller, but drag coefficient became larger.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.56 no.2
• /
• pp.1-9
• /
• 2014

A modeling system is constructed by integrating an one-dimensional unsteady flow simulation model and a hydrologic model to simulate flood flows in drainage channel networks of paddy field districts. The modeling system's applicability is validated by simulating flood discharges from a paddy field district, which consists of nine paddy fields and one drainage channel. The simulation results are in good agreement with the observed. Particularly, in the verification stage, the relative errors of peak flows and peak depths between the observed and simulated hydrographs range 8.96 to 10.26 % and -10.26 to 2.97 %, respectively. The modeling system's capability is compared with that of a water balance equation-based model; it is revealed that the modeling system's accuracy is superior to the other model. In addition, the simulations of flood discharges from large-sized paddy fields through drainage channels show that the flood discharge patterns are affected by drainage outlet management for paddy fields and physical characteristics of the drainage channels. Finally, it is concluded that to efficiently design drainage channel networks, it is necessary to analyze the results from simulating flood discharges of the drainage channel networks according to their physical characteristics and connectivities.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.4
• /
• pp.48-59
• /
• 1996

The numerical step in the unsteady viscous flow analysis can be divided in the space analysis step satisfying continuity equation and the time marching step. In this study the spectral method is applied to solve the pressure Poisson equation in the space analysis step. If the highest order differential term of the pressure Poisson equation is transformed by Fourier series, pressure arid its first derivatives can be expressed by the integral form of Fourier series. So Gibb's phenomena can be eliminated and the spectral method can be applied to non-periodic problems. The numerical analysis of unsteady viscous flow around 2-dimensional circular cylinder and wing is carried out and compared for verification.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.95-104
• /
• 2002

Numerical solution of the full Navier-Stokes equation as well as the energy equation has been obtained for the unsteady natural convection in a rectangular enclosure. One side wall was maintained at very high temperature simulating fires. Especially the effect of surface radiation was taken into account. While the enclosed air was assumed to be transparent, the internal walls directly interacted one another through the surface radiation. Due to a significant temperature difference in the flow field, the equation of state was used instead of the Boussinesq approximation. It was found that the rapid heating of the adiabatic ceiling and floor by the incoming radiation from the hot wall made the evolution at thermo-fluid field highly unstable in the initial period. Therefore, the secondary cells brought about at the floor region greatly affected the heat transfer mechanism inside the enclosure. The heat transfer rate was augmented by the radiation, resulting in requiring less time for the flow to reach the steady state. At the steady state neglecting radiation two internal hydraulic jumps were clearly observed in upper/left as well as in lower/right comer. However, the hydraulic jump in the lower/right comer could not be observed for the case including radiation due to its high momentum flow over the bottom wall. Radiation resulted in a faster establishment of the steady state phenomena.

• Journal of computational fluids engineering
• /
• v.6 no.2
• /
• pp.22-31
• /
• 2001

The two-dimensional unsteady incompressible viscous flow behind rectangular bluff bodies between two parallel plates was analyzed numerically. The steady state flow and the vortex flow behind rectangular bluff bodies in the channel were investigated for two regimes i.e., the laminar(Re = 100, 300, 500) and the turbulent flows(Re = 10⁴∼10/sup 6/). The vortex shedding was generated by a physical disturbance(6%) numerically imposed at the rear of the bluff bodies for a short time. It was observed that the perturbed flow became periodic after a transient period. And in the case of unsteady inflow, the sinusoidal pulsatile flow was applied as the inlet condition in the turbulent flow of Reynolds number of 1.0×10/sup 5/. FLUENT code was employed to solve the problems. The power-law scheme was used to get stable linearized equations and the PISO algorithm was applied to finding the solution of them.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.30-37
• /
• 2005

To evaluate LBM we performed the simulation of the unsteady two dimensional flow over a square cylinder in a channel in moderate Reynolds number range, $100\~500$ by using LBM and Fractional-Step method. Frist of all we compared LBM solution of Poiseuille flow applied Farout and periodic boundary conditions with the analytical solution to verify the applicability of the boundary conditions. For LBM simulation the calculation domain was formed by structured 500x100 grids. Prescribed maximum velocity and density inlet and Farout boundary conditions were imposed on the in-out boundaries. Bounceback boundary condition was applied to the channel and the cylinder waifs. The flow patterns and vortex shedding strouhal numbers were compared with previous research results. The flow patterns by LBM were in agreement with the flow pattern by fractional step method. Furthermore the strouhal number computed by LBM simulation result was more accurate than that of fractional step method through the comparison of the previous research results.

• Wind and Structures
• /
• v.39 no.1
• /
• pp.31-45
• /
• 2024

This paper examines the flow characteristics around an inclined prism with a U-shaped cross-section ("U-profile") and investigates the connection between the flow and flow-induced vibrations. The study employs a combined approach that involves wind tunnel experiments and computational fluid dynamics (CFD) using an unsteady Reynolds-averaged Navier-Stokes (RANS) turbulence model. Distinct vortex formation patterns are observed in the flow field surrounding the stationary inclined profile. When the cavity of the profile faces away from the incoming flow, large vortices develop behind the profile. Conversely, when the cavity is oriented towards the oncoming flow, these vortices form within the cavity. Notably, due to the slow movement of these large vortices through the cavity, the frequency at which vortices are shed in the negative inclination case is lower compared to the positive inclination, where they form in the wake. Wind tunnel experiments reveal an intermittent transition between the two vortex formation patterns at zero inclination. Large vortices sporadically emerge both in the cavity and behind the profile. The simulation results demonstrate that when these large vortices occur at a frequency close to the structure's natural frequency, they induce prominent pitch vibrations. This phenomenon is also sought after and presented in coupled vibration experiments. Additionally, the simulations indicate that when the natural frequency of the structure is considerably lower than the vortex shedding frequency, this type of vibration can be observed.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.7-13
• /
• 2010

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 5.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}10^6$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental datum in the low aspect ratio cavity (L/D = ~4.5). In the high aspect ratio cavity, however, there are other low dominant frequencies of the leading edge shear layer with the dominant frequencies of the feedback mechanism.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.13-18
• /
• 2009

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought about by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 2.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}106$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental data in the low aspect ratio cavity (L/D = ~ 4.5). In the large aspect ratio cavity, however, there are other low dominant frequencies due to the leading edge shear layer with the dominant frequencies of the feedback mechanism. The characteristics of the acoustic wave propagation are analyzed using the Correlation of Pressure Distribution (CPD).