• Journal of the Korea Institute of Military Science and Technology
• /
• v.12 no.6
• /
• pp.704-710
• /
• 2009

In this study, a comparison study of flutter analysis for the AGARD 445.6 wing with wind turnnel test data has been conducted in the subsonic, transonic and supersonic flow regions. Nonlinear aeroelastic using FSIPRO3D which is a generalized user-friendly fluid-structure analyses have been conducted for a 3D wing configuration considering shockwave and turbulent viscosity effects. The developed fluid-structure coupled analysis system is applied for aeroelastic computations combining computational structure dynamics(CSD), finite element method(FEM) and computations fluid dynamics(CFD) in the time domain. MSC/NASTRAN is used for the vibration analysis of a wing model, and then the result is applied to the FSIPRO3D module. the results for dynamic aeroelastic response using different turbulent models are presented for several Mach numbers. Calculated flutter boundary are compared with the wind-tunnel experimental and the results show very good agreements.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition. The velocity field is obtained by using the Biot-Savart integral. The Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well established for potential flow analysis. The present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

• Structural Engineering and Mechanics
• /
• v.69 no.2
• /
• pp.177-191
• /
• 2019

The local topography has a significant effect on the characteristics of seismic ground motion. This paper investigates the influence of topographic effects on the seismic response of a train-bridge system. A 3-D finite element model with local absorbing boundary conditions is established for the local site. The time histories of seismic ground motion are converted into equivalent loads on the artificial boundary, to obtain the seismic input at the bridge supports. The analysis of the train-bridge system subjected to multi-support seismic excitations is performed, by applying the displacement time histories of the seismic ground motion to the bridge supports. In a case study considering a bridge with a span of 466 m crossing a valley, the seismic response of the train-bridge system is analyzed. The results show that the local topography and the incident angle of seismic waves have a significant effect on the seismic response of the train-bridge system. Leaving these effects out of consideration may lead to unsafe analysis results.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.388-397
• /
• 2010

A two-dimensional time-domain, potential-theory-based fully nonlinear numerical wave tank (NWT) was developed by using boundary element method and the mixed Eulerian-Lagrangian (MEL) approach for free-surface node treatment. The NWT was applied to prediction of primary wave energy conversion efficiency of a bottom-mounted oscillating water column (OWC) wave power device. The nonlinear free-surface condition inside the chamber was specially devised to represent the pneumatic pressure due to airflow velocity and viscous energy loss at the chamber entrance due to wave column motion. The newly developed NWT technique was verified through comparison with given experimental results. The maximum energy extraction was estimated with various chamber-air duct volume ratios.

• Journal of Ocean Engineering and Technology
• /
• v.28 no.6
• /
• pp.508-516
• /
• 2014

A numerical method to investigate the non-linear motion characteristics of a TLP is established. A time domain simulation that includes the memory effect using the convolution integral is used to consider the transient effect of TLP motion. The hydrodynamic coefficients and wave force are calculated using a potential flow model based on the HOBEM(higher order boundary element method). The viscous drag force acting on the platform and tendons is also considered by using Morison’s drag. The results of the present numerical method are compared with experimental data. The focus is the nonlinear effect due to the viscous drag force on the TLP motion. The ringing, springing, and drift motion are due to the drag force based on Morison's formula.

• Journal of Environmental Science International
• /
• v.13 no.7
• /
• pp.655-660
• /
• 2004

This paper demonstrates the numerical simulation of three dimensional flow pattern for vehicular exhaust dispersion in the street canyons. The wind flow around buildings in urban is computed by the SIMPLEST method. The convection-diffusion equation was used to compute the $NO_X$ concentration level near buildings. Details are given of important boundary conditions and turbulence quantities variations. The simple turbulence model was used for unisotropic viscous effect. A control-volume based finite-difference method with the upwind scheme is employed for discretization equation. The simple turbulence model applied in this study has been verified through comparison between predicted and measured data near buildings. By the predictive results, the updraft induced by the presence of high-rise buildings is important in the transport of street level pollutant out from the street canyons. Our suggestion for reducing ground level pollution is to have high-rise buildings constructed or to reduce the channelling effect of street canyons.

• Journal of Navigation and Port Research
• /
• v.27 no.5
• /
• pp.473-478
• /
• 2003

Numerical study was made on the flow characteristics around a circular pipeline between parallel walls. The incompressible Navier-Stokes equations were solved by using a third-order upwind differential scheme. When the distance near a wall is small enough, the vortex shedding is almost completely suppressed because of the interaction with the wall boundary layer separation. This study aims to clarify the characteristics of the vortex shedding regime as the body approaches a wall as Reynolds number varies. The feature of separated vorticity dynamics is analyzed at different conditions with particular attention to the interaction between the pipeline wake and the induced separation on the plane walls.

• Journal of Ship and Ocean Technology
• /
• v.1 no.2
• /
• pp.19-30
• /
• 1997

The accurate prediction of the flow and the pressure distribution near the tip of the blade is crucial in determining the tip vortex cavitation inception which usually occurs on the blade tip or inside the core of the tip vortex just downstream of the blade tip. An improved boundary element method is applied to the prediction of the flow around propeller blades, with emphasis at the tip region. In the method, the Blow adapted grid and a higher order panel method, which combines a hyperboloidal panel geometry with a hi-quadratic dipole distribution, are used in order to accurately model the trailing wake geometry and the highly rolled-up regions in the wake. The method is applied to several propeller geometries and the results have been found to agree well to the existing experimental data. Inviscid flow methods are able to predict the pressures at the tip as well as the shape of the trailing wake. On the other hand, they are unable to determine the flow inside the viscous core of the tip vortex, where cavitation inception often occurs. Thus, a method is presented that treats the flow inside the viscous core. The inner flow is treated with a 2-D Clavier-stokes solution without making any assumptions for axisymmetric flow and conicity of the flow along the tip trajectory. The method can thus allow the treatment of general propeller blade configurations. The velocity and pressure distributions inside the core are shown and compared to those from other numerical methods.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.407-413
• /
• 2011

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).