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

Parabolized stability equations for compressible flows in general curvilinear coordinate system are derived to deal with a broad range of transition prediction problems on complex geometry. A highly accurate finite difference PSE code has been developed using an implicit marching procedure. Compressible and incompressible flat plate flow stability under two-dimensional and three¬dimensional disturbances has been investigated to test the present code. Results of the present computation are found to be in good agreement with the multiple scale analysis and DNS data. Stability calculation results by the present PSE code for compressible boundary layer at Mach numbers ranging from 0.02 to 1.5 are also presented and are again seen to be as accurate as the spectral method.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.369-378
• /
• 2002

The application of Large Eddy Simulation (LES) in a curvilinear coordinate system to the flow around a square cylinder is presented. In order to obtain sufficient resolution near the side of the cylinder, we use an O-type grid. Even with a curvilinear coordinate system, it is difficult to avoid the numerical oscillation arising in high-Reynolds-number flows past a bluff body, without using an extremely fine grid used. An upwind scheme has the effect of removing the numerical oscillations, but, it is accompanied by numerical dissipation that is a kind of an additional sub-grid scale effect. Firstly, we investigate the effect of numerical dissipation on the computational results in a case where turbulent dissipation is removed in order to clarify the differences between the effect of numerical dissipation. Next, the applicability and the limitations of the present method, which combine the dynamic SGS model with acceptable numerical dissipation, are discussed.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.134-140
• /
• 2007

Nonlinear parabolized stability equations for compressible flow in general curvilinear coordinate system are derived to deal with a broad range of transition prediction problems on complex geometry. A highly accurate finite difference PSE code has been developed using an implicit marching procedure. Blasius flow is tested. The results of the present computation show good agreement with DNS data. Nonlinear interaction can make the T-S fundamental wave more unstable and the onset of its amplitude decay is shifted downstream relative to linear case. For nonlinear calculations, rather small difference in initial amplitude can produce large change during nonlinear region. Compressible secondary instability at Mach number 1.6 is also simulated and showed that 1.1% initial amplitude for primary mode is enough to trigger the secondary growth.

• Journal of the Society of Naval Architects of Korea
• /
• v.58 no.3
• /
• pp.191-197
• /
• 2021

This study addresses autonomous ship collision avoidance in narrow channels using curvilinear coordinates. Navigation in narrow channels or fairways is known to be much more difficult and challenging compared with navigation in the open sea. It is not straightforward to apply the existing collision avoidance framework designed for use in the open sea to collision avoidance in narrow channels due to the complexity of the problem. In this study, to generalize the autonomous navigation procedure for collision avoidance in narrow channels, we introduce a curvilinear coordinate system for collision-free path planning using a parametric curve, B-spline. To demonstrate the feasibility of the proposed algorithm, ship traffic simulations were performed and the results are presented.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.2
• /
• pp.198-209
• /
• 1999

This study reports the selection of dependent variables for momentum equations in general curvilinear coordinates. Catesian, covariant and contravariant velocity components were examined for the dependent variable. The focus of present study is confined to staggered grid system Each dependent variable selected for momentum equations are tested for several flow fields. Results show that the selection of Cartesian and covariant velocity components intrinsically can not satisfy mass conservation of control volume unless additional converting processes ore used. Also, Cartesian component can only be used for the flow field in which main-flow direction does not change significantly. Convergence rate for the selection of covariant velocity component decreases quickly as with the increase of non-orthogonality of grid system. But the selection of contravariant velocity component reduces the total mass residual of discretized equations rapidly to the limit of machine accuracy and the solutions are insensitive to the main-flow direction.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.11
• /
• pp.1500-1508
• /
• 2001

This paper represents the importance of dependent variables in non-orthogonal curvilinear coordinates just as the importance of those variables of convective scheme and turbulence model in computational fluid dynamics. Each of Cartesian, physical covariant and physical contravariant velocity components was tested as the dependent variables of momentum equations in the staggered grid system. In the flow past a circular cylinder, the results were computed to use each of three variables and compared to experimental data. In the skewed driven cavity flow, the results were computed to check the grid dependency of the variables. The results used in Cartesian and physical contravariant components of velocity in cylinder flow show the nearly same accuracy. In the case of Cartesian and contravariant component, the same number of vortex was predicted in the skewed driven cavity flow. Vortex strength of Cartesian component case has about 30% lower value than that of the other two cases.

• Journal of applied mathematics & informatics
• /
• v.32 no.5_6
• /
• pp.783-790
• /
• 2014

Most distance functions, including taxicab distance, are defined on Cartesian plane, and recent studies on distance functions have been mainly focused on Cartesian plane. However, most streets in cities include not only straight lines but also curves. Therefore, there is a significant need for a distance function to be defined on a curvilinear coordinate system. In this paper, we define a new function named polar taxicab distance, using polar coordinates. We prove that this function satisfies the conditions of distance function. We also investigate the geometric properties and classifications of circles in the plane with polar taxicab distance.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.402-407
• /
• 2006

본 연구는 곡선좌표계에서 유한차분기법(finite difference method)을 이용하여 2차원 흐름이 모의가능한 수치모형을 개발하는 것이다. 기존의 연구는 대부분 직교좌표계(cartesian coordinate system)에서의 격자망을 대상으로 개발되고 적용되었기 때문에 불규칙한 흐름의 경계 및 형상을 올바로 표현하기 어려웠다. 유한요소법이나 유한체적법같은 수치모의기법들이 개발되어 비구조격자체계를 구성하고 자연현상에 가까운 경계 표현할 수 있도록 개발되었다. 하지만 위의 기법들은 질량과 운동량과 같은 물리량을 보존하기 위해서 매우 조밀한 격자체계를 가져야만 한다. 이에 본 연구에서는 기존의 문제점들을 해결하기 위하여 곡선좌표계(curvilinear coordinate system)를 이용하여 지배방정식을 표현하고 2차원 흐름을 모의할 수 있는 모형을 구축한다. 수치모형은 leap-frog기법과 1차 정확도의 풍상차분기법(upwind scheme)을 사용하여 구성하였다. 본 연구에서 개발된 모형을 사각수조 및 만곡수로흐름에 적용하여 모의결과를 해석해 및 실험관측값과 비교하였다. 이로부터 본 수치모형이 해석해 및 실측치와 잘 일치하고 있음을 알 수 있었다.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• v.26 no.1
• /
• pp.1-14
• /
• 1998

Various numerical methods for the two dimensional shallow water equations have been applied to the problems of flood routing, tidal circulation, storm surges, and atmospheric circulation. These methods are often based on the Alternating Direction Implicity(ADI) method. However, the ADI method results in inaccuracies for large time steps when dealing with a complex geometry or bathymetry. Since this method reduces the performance considerably, a fully implicit method developed by Wilders et al. (1998) is used to improve the accuracy for a large time step. Finite Difference Methods are defined on a rectangular grid. Two drawbacks of this type of grid are that grid refinement is not possibile locally and that the physical boundary is sometimes poorly represented by the numerical model boundary. Because of the second deficiency several purely numerical boundary effects can be involved. A boundary fitted curvilinear coordinate transformation is used to reduce these difficulties. It the curvilinear coordinate transformation is used to reduce these difficulties. If the coordinate transformation is orthogonal then the transformed shallow water equations are similar to the original equations. Therefore, an orthogonal coorinate transformation is used for defining coordinate system. A multigrid (MG) method is widely used to accelerate the convergence in the numerical methods. In this study, a technique using a MG method is proposed to reduce the computing time and to improve the accuracy for the orthogonal to reduce the computing time and to improve the accuracy for the orthogonal grid generation and the solutions of the shallow water equations.

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

The development of unsteady three-dimensional incompressible viscous solver based on unsteady physical curvilinear coordinate system is presented. A 12-point finite analytic scheme based on local uniform grid spacing is extended for nonuniform grid spacing. The formulation of a condition is suggested to avoid the oscillation of the series summations produced by the application of the method of separation of variables. SIMPLER and pressure Poisson equation techniques are used for solving a velocity-pressure coupled problem. The matrix is solved using the Generalized Minimal RESidual (GMRES) method to enhance the convergence rate of unsteady flow solver and the Kinematic boundary condition of a free surface flow. It is demonstrated that the numerical solutions of these equations are less mesh sensitive.