• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.410-415
• /
• 2001

An analytical and numerical examination of second-order fractional-step methods and boundary condition for the incompressible Navier-Stokes equations is presented. In this study, the compatibility condition for pressure Poisson equation and its boundary conditions, stability, and numerical accuracy of canonical fractional-step methods has been investigated. It has been found that satisfaction of compatibility condition depends on tentative velocity and pressure boundary condition, and that the compatible boundary conditions for type D method and approximately compatible boundary conditions for type P method are proper for divergence-free velocity for type D and approximately divergence-free for type P method. Instability of canonical fractional-step methods is induced by approximation of implicit viscous term with explicit terms, and the stability criteria have been founded with simple model problems and numerical experiments of cavity flow and Taylor vortex flow. The numerical accuracy of canonical fractional-step methods with its consistent boundary conditions shows second-order accuracy except $D_{MM}$ condition, which make approximately first-order accuracy due to weak coupling of boundary conditions.

• Wind and Structures
• /
• v.26 no.6
• /
• pp.355-367
• /
• 2018

This paper studies the aerodynamic stability of a tensioned, geometrically nonlinear orthotropic membrane structure with hyperbolic paraboloid in sag direction. Considering flow separation, the wind field around membrane structure is simulated as the superposition of a uniform flow and a continuous vortex layer. By the potential flow theory in fluid mechanics and the thin airfoil theory in aerodynamics, aerodynamic pressure acting on membrane surface can be determined. And based on the large amplitude theory of membrane and D'Alembert's principle, interaction governing equations of wind-structure are established. Then, under the circumstance of single-mode response, the Bubnov-Galerkin approximate method is applied to transform the complicated interaction governing equations into a system of second-order nonlinear differential equation with constant coefficients. Through judging the frequency characteristic of the system characteristic equation, the critical velocity of divergence instability is determined. Different parameter analysis shows that the orthotropy, geometrical nonlinearity and scantling of structure is significant for preventing destructive aerodynamic instability in membrane structures. Compared to the model without considering flow separation, it's basically consistent about the divergence instability regularities in the flow separation model.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.55-63
• /
• 1995

The flow field around a high-speed train including cross-wind effects has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations in the inertial frame using the iterative time marching scheme. The governing equations are differenced with 1st-order accurate backward difference scheme for the time derivatives, 3th-order accurate QUICK scheme for the convective terms and 2nd-order accurate central difference scheme for the viscous terms. The Marker-and-Cell concept was applied to efficiently solve continuity equation, which is differenced with 2nd-order accurate central difference scheme. The 4th-order artificial damping is added to the continuity equation for numerical stability. A C-H type of elliptic grid system is generated around a high-speed train including ground. The Baldwin-Lomax turbulent model was implemented to simulate the turbulent flows. To validate the present procedure, the flow around a high speed train at constant yaw angle of $45^{\circ}\;and\;90^{\circ}$ has been simulated. The simulation shows 3-D vortex generation in the lee corner. The flow separation is also observed around the rear of the train. It has concluded that the results of present study properly agree with physical flow phenomena.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.46-53
• /
• 1997

The flow field around a three dimensional minivan-like body has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using second-order accurate schemes for the time derivatives, and third/second-order scheme for the spatial derivatives. The Marker-and-Cell concept is applied to efficiently solve continuity equation. A H-H type of multi-block grid system is generated around a three dimensional minivan-like body. Turbulent flows have been modeled by the Baldwin-Lomax turbulent model. To validate present procedure, the flows around the Ahmed body with 12.5° of slant angle are simulated. A good agreement with other numerical results is achived. After code validation, the flows around a mimivan-like body are simulated. The simulation shows three dimensional vortex-pair just behind body. The flow separation is also observed on the rear of the body. It has concluded that the results of present study properly agreed with physical flow phenomena.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.2
• /
• pp.174-180
• /
• 2004

Failure of tube bundles due to excessive flow-induced vibrations continues to affect the performance of nuclear power plant Early experimental studies concentrated on rigid structures and later investigators dealt with elastic structures because of their importance in many engineering fields. On the other hand, much less numerical work has been carried out, because of the numerical complexity associated with the problem. Conventional approaches usually decoupled the flow solution from the structural problem. The present numerical study proposes the methodology in analyzing the fluidelastic instability occurring in tube bundles by coupling the Computational fluid Dynamics (C%) with the tube equation of motions. The motion of the structures is modeled by a spring-damper-mass system that allows transnational motion in two directions (a two-degree-of-freedom system). The fluid motion and the cylinder response are solved in an iterative way, so that the interaction between the fluid and the structure can be accounted for property. The aim of the present work is to predict the fluidelstic instability of tube bundles and the associated phenomena, such as the response of the cylinder, the unsteady lift and drag on the cylinder, the vortex shedding frequency.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.658-663
• /
• 2003

A process of 3-D stereo particle image velocimetry(PIV)was developed for the measurement of an illuminated sliced section field of 3-D complex flows. The present method includes modeling of camera by a calibrator based on the homogeneous coordinate system, transformation of the oblique-angled image to the right-angled image, identification of 2-D velocity vectors by 2-D cross-correlation equation, stereo matching of 2-D velocity vectors of two cameras, accurate calculation of 3-D velocity vectors by homogeneous coordinate system, removal of error vectors by a statistical method followed by a continuity equation criteria, and finally 3-D display as the post processing. An experimental system was also used for the application of the proposed method. Two high speed digital CCD cameras and an Argon-Ion Laser for the illumination were adopted to clarify the time-dependent characteristics of the leading edge extension(LEX) in a highly swept shape applied to a delta wing found in modern air-fighters.

• Journal of the Korean Society of Visualization
• /
• v.10 no.1
• /
• pp.27-33
• /
• 2012

This paper describes a new design of small-scale horizontal wind blade, called spiral wind turbine blade. Theoretical and numerical approaches on the prediction of aerodynamic performance of the blade have been conducted. A theoretical equation is successfully derived using the angular momentum equation to predict aerodynamic characteristics according to the design shape parameters of spiral blade. To be compared with the theoretical value, a numerical simulation using ANSYS CFX v12.1 is performed on the same design with the theoretical one. Large scale tip vortex is captured and graphically presented in this paper. The TSR-$C_p$ diagram shows a typical parabolic relation in which the maximum efficiency of the blade approximately 25% exists at TSR=2.5. The numerical simulation agrees well with that of the theoretical result except at the low rotational speed region of 0~20 rad/s.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.206-218
• /
• 2014

The wave attenuation by floating breakwaters in high amplitude waves, which can lead to wave overtopping and breaking, is examined by numerical simulations. The governing equations, the Navier-Stokes equations and the continuity equation, are calculated in a fixed Cartesian grid system. The body boundaries are defined by the line segment connecting the points where the grid line and body surface meet. No-slip and divergence free conditions are satisfied at the body boundary cell. The nonlinear waves near the moving body is defined using the modified marker-density method. To verify the present numerical method, vortex induced vibration on an elastically mounted cylinder and free roll decay are numerically simulated and the results are compared with those reported in the literature. Using the present numerical method, the wave attenuations by three kinds of floating breakwaters are simulated numerically in a regular wave to compare the performance.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.3
• /
• pp.342-351
• /
• 1999

A jet injected normally into a cross flow has been found to have the cross section of a horseshoe shape. It occurs by a twin vortex motion in the region downstream of the jet injection. Such a flow is inherently and highly three-dimensional and numerical calculations should play an important role. The three-dimensional momentum equations with buoyancy effect and energy equation are solved to obtain the velocity distributions, center-line trajectories, cross sectional shape and entrainment. The density difference is sufficiently small, so that the Boussinesq approximation is considered to be valid. The SIMPLE algorithm is applied in a staggered grid system of a calculational domain for the numerical method.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.754-759
• /
• 2001

This paper presents unsteady computational investigations and wind tunnel tests on the flow field around a square cylinder with a gap between the body and the ground plane. Two-dimensional unsteady, incompressible Navier-Stokes codes are developed for the computation of the viscous turbulent flows. By computing the flow around a square cylinder without ground effect, three two-equation turbulence models are evaluated and the developed code is validated. The results show a good agreement with experimental values and other computational results. Critical gap height at which the formation of Karman vortex streets is interrupted, is demonstrated and another transition regime is pointed out