• The Journal of the Acoustical Society of Korea
• /
• v.23 no.8
• /
• pp.576-582
• /
• 2004

A computer code that is based on the Pseudo-spectral finite difference algorithm using staggered grid is developed for the wave propagation modeling in the time domain. The advantage of a finite difference approximation is that any geometrically complicated media can be modeled. Staggered grids are advantageous as it provides much more accuracy than using a regular grid. Pseudo-spectral methods are those that evaluate spatial derivatives by multiplying a wavenumber by the Fourier transform of a pressure wave-field and performing the inverse Fourier transform. This method is very stable and reduces memory and the number of computations. The synthetic results by this algorithm agree with the analytic solution in the infinite and half space. The time domain modeling was implemented in various models. such as half-space. Pekeris waveguide, and range dependent environment. The snapshots showing the total wave-field reveals the Propagation characteristic or the acoustic waves through the complex ocean environment.

• Proceedings of the Acoustical Society of Korea Conference
• /
• autumn
• /
• pp.425-428
• /
• 2004

본 논문에서는 시간영역에서 해양 음 전달 해석을 위해 엇갈림 격자에서 유사 스펙트럴 알고리듬을 기반으로 한 전산조직을 개발하였다. 유사 스펙트럴 방법은 파수 영역에서 파수에 음압의 푸리에 변환을 곱한 후 이를 역 푸리에 변환함으로서 공간 도함수를 구하는 방법이다. 유사 스펙트럴 방법은 빠른 푸리에 변환법의 사용으로 계산속도가 빠르며, 엇갈림 격자에서 이 방법을 사용하면 음 전달 현상을 정확하고 안정되게 모사할 수 있다. 무한 및 반무한 영역에서 이 알고리듬에 의한 결과가 해석해와 잘 일치함을 확인하였고, 다양한 해양환경에서 시간영역 모델링을 수행하여 스냅사진을 얻어내었으며 이 스냅사진을 통해 복잡한 해양환경에서 신호의 전파 현상을 파악할 수 있었다.

• 한국방재학회:학술대회논문집
• /
• 2009.02b
• /
• pp.74.1-74.1
• /
• 2009

• Journal of Energy Engineering
• /
• v.17 no.4
• /
• pp.218-226
• /
• 2008

A component-scale two-phase analysis code has been developed for a realistic simulation of two-phase flow transients in a light water nuclear reactor component. In the code, a two-fluid three-field model is adopted and the governing equations are solved on an unstructured mesh. For the numerical solution scheme, the semi-implicit method used in the RELAP5 code was selected, which has been proved to be very stable and accurate for most of practical applications. However, some modifications were needed for its application to an unstructured non-staggered grid. This paper presents the modified semi-implicit numerical method for unstructured grid and the preliminary results of the calculations.

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.1-8
• /
• 2007

A code is developed using the hybrid Cartesian/immersed boundary method and it is applied to simulate flows around a three-dimensional deforming body. A new criterion is suggested to distribute the immersed boundary nodes based on edges crossing a body boundary. Velocities are reconstructed at the immersed boundary nodes using the interpolation along a local normal line to the boundary. Reconstruction of the pressure at the immersed boundary node is avoided using the hybrid staggered/non-staggered grid method. The developed code is validated through comparisons with other experimental and numerical results for the velocity profiles around a circular cylinder under the forced in-line oscillation and the pressure coefficient distribution on a sphere. The code is applied to simulate the flow fields around a plate whose tail is periodically flapping under a translation. The effects of the velocity and acceleration due to the deformation on the periodic shedding of pairs of tip vortices are investigated.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.20 no.4
• /
• pp.77-86
• /
• 2016

A numerical analysis is performed to investigate the effect of a supply channel and jet hole arrangement on the heat flow characteristics of impingement jet. The jet holes in a supply channel are composed of a single or staggered array from the center of a leading edge channel. The software ICEMCFD is used to generate the structured grids for calculation domain and a CFD code CFX 15.0 to perform the simulation. The present solutions are validated by comparison with the experimental and numerical ones of others. A comparison of mass flow rates of impingement jets and Nusselt numbers on the impingement surface for the single or staggered arrays is made.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.37-45
• /
• 2007

The hybrid Cartesian/immersed boundary method is applied to simulate fluid-structure interaction of a two-dimensional orbiting flexible foil. The elastic deformation of the flexible foil is modelled based on the dynamic equation of a thin-plate. At each time step, the locations and velocities of the Lagrangian control points on the flexible foil are used to reconstruct the boundary conditions for the flow solver based on the hybrid staggered/non-staggered grid. To test the developed code, the flow fields around a flapping elliptical wing are calculated. The time history of the vertical force component and the evolution of the vorticity fields are compared with recent other computations and good agreement is achieved. For the orbiting flexible foil, the vorticity fields are compared with those of the case without the deformation. The combined effects of the angle of attack and the orbit on the deformation are investigated. The grid independency study is carried out for the computed time history of the deformation at the tip.

• Journal of computational fluids engineering
• /
• v.7 no.1
• /
• pp.10-19
• /
• 2002

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.93-93
• /
• 2016

최근 증가하고 있는 기후변화에 의해 설계빈도를 상회하는 강우의 발생빈도가 증가하고 있으며, 이로 인한 도시유역의 내수범람 피해가 증가하고 있다. 도시유역에서 발생하는 침수 피해의 경우 인적 물적 자원이 집중되어 있는 도시의 특성으로 인해 침수로 인한 직접적 피해 규모가 상당할 뿐만 아니라 침수 발생 후 세균 및 박테리아에 의해 발생하는 수인성 전염병의 유행 등과 같은 2차적 피해 또한 심각한 사회적 비용을 초래할 수 있어 도시유역의 침수 피해를 저감시키기 위한 대책이 절실히 요구되어지고 있다. 도시유역의 침수를 예방하기 위한 대책은 구조적 비구조적 대책으로 구분되어 질 수 있으며 구조적 대책의 경우 침수 피해 예방에 직접적인 효과를 낼 수 있다는 장점이 있으나 대규모 사업예산 및 사업 기간으로 인해 직접적 효과를 보기까지 상대적으로 긴 시간이 필요할 뿐만 아니라 사업 진행 중 대상지역 거주민들의 민원으로 인한 갈등 조정 등으로 인해 사업실행에 어려움을 겪고 있다. 이러한 측면에서 비구조적 대책의 일환인 수치해석을 통한 침수피해 재현 및 침수원인 파악을 통한 구조개선 제안은 구조적 대책의 단점을 보완할 수있는 좋은 대안이 될 수 있다. 도시유역의 경우 비도시유역과 대조적인 차이점으로는 높은 비율의 불투수층, 복잡한 지형, 다수의 인공 구조물 및 배수관망 시스템 등을 들 수 있으며, 침수해석 모형의 정확도를 높이기 위해서는 복잡한 지형의 효율적인 처리가 무엇보다 중요하다. 일반적으로 이용되는 2차원 침수해석 모형들은 직교구조 격자 또는 비구조 격자를 이용하여 지형을 묘사하고 있으며 DEM 자료를 직접 사용하는 직교구조 격자의 경우 지형 데이터 생성이 상대적으로 쉽다는 장점이 있으나 복잡한 지형을 표현하기 위해서는 불필요한 지역까지 높은 해상도를 이용해야 하며 이로 인하여 모의시간이 지나치게 길어지는 문제점이 발생한다. 비구조 격자의 경우 상대적으로 복잡한 도시 유역을 잘 묘사할 수 있다는 장점이 있으나 격자망 생성에 필요한 데이터가 많고 격자망 생성에 지나치게 많은 시간과 노력이 소요된다는 단점이 있다. 따라서 본 연구에서는 위에서 언급한 두 가지 방법의 장점만을 취할 수 있도록 메쉬 세분화 기법을 이용한 2차원 침수해석 모형을 개발 하여 복잡한 지형은 고해상도 메쉬를 이용하여 보다 자세히 묘사하고 상대적으로 복잡하지 않은 지형은 저해상도 메쉬를 이용하여 계산시간을 단축시킬 수 있도록 하였다. 수치해석 기법으로는 엇갈림 격자를 이용하는 Leap-Frog 기법과 유한차분 (Finite difference Method)기법을 이용하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.290-297
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second order scheme.