• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 추계 학술대회논문집
• /
• pp.11-17
• /
• 2003

The new multi-dimensional higher order interpolation scheme called MHIS is developed. Firstly, multi-dimensional TVD condition is derived based on one-dimensional TVD condition. Using multi-dimensional TVD condition, 2nd, 3rd and 5th order MHIS are presented. By help of multi-dimensional TVD condition, it is possible to captured a discontinuity monotonically even in a multi-dimensional flow. It is verified through several test cases that the accuracy and the robustness of MHIS are enhanced in regions of shock discontinuities as well as boundary-layers.

• 한국해안·해양공학회논문집
• /
• 제30권4호
• /
• pp.143-152
• /
• 2018

최근 개발된 FUNWAVE-TVD 파랑모형을 이용하여 적용되어 온 TVD 기법들의 수렴도와 수치적인 안정성을 비교하였다. Yamamoto and Daiguji(1993)의 minmod limiter를 사용하는 4차 정확도의 MUSCL-TVD 기법과 Erduran et al.(2005)의 van-Leer limiter를 사용하는 4차 정확도의 MUSCL-TVD 기법, Zhou et al.(2001)의 van-Leer limiter를 사용하는 2차 정확도의 MUSCL-TVD 기법을 비교하였으며, 수리실험 관측치가 제시되어 있는 Vincent and Briggs(1989)의 불규칙 파랑실험에 적용하였다. 불규칙 파랑의 비쇄파 실험 결과에서 minmod limiter를 사용하는 4차 정확도의 기법은 van-Leer limiter를 사용하는 기법이 요구하는 격자의 크기만큼 세밀한 격자를 요구하지는 않지만, 더 낮은 CFL을 사용해야 안정적인 모의가 가능하였다. 반면에 van-Leer limiter를 사용하는 기법에서는 numerical dissipation을 줄이기 위하여 보다 세밀한 격자를 필요로 하지만 비교적 높은 CFL을 사용할 수 있는 것으로 나타났다. 각 기법의 numerical dissipation의 영향을 최대한 줄이기 위하여 공간격자를 충분히 줄인 쇄파 모의 실험에서는 비쇄파 실험에 비하여 각 기법의 특성이 명확히 나타났다. Numerical dissipation이 상대적으로 작은 minmod limiter를 사용하는 기법으로 모의할 때는 격자를 충분히 줄이면 수치적인 불안정성이 나타나며 수치해가 발산하는 결과를 보였지만, van-Leer limiter를 사용하는 기법에서는 비교적 낮은 CFL을 사용하여 쇄파 모의가 완료되었으며, 관측치를 잘 재현하는 결과를 보였다.

• 한국수자원학회논문집
• /
• 제34권2호
• /
• pp.177-186
• /
• 2001

2차 이상의 정확도를 가지며 불연속면에서 수치진동이 발생하지 않는 수치모형의 개발을 위해 풍상차분기법에 기초한 TVD기법이 소개되었다. 이 수치모형을 불연속면 이 존재하는 경우와 존재하지 않는 경우에 대해 적용하였으며, 이 결과 1차 정확도의 풍상차분기법은 시간이 지나면서 수치점성의 영향이 커졌으며 2차 정확도의 Lax-Wendroff기법의 경우에는 불연속면에서 진동이 발생하였다. 그러나 TVD기법은 모든 경우에서 만족스러운 결과를 예측하였다.

• 한국수자원학회논문집
• /
• 제36권4호
• /
• pp.597-608
• /
• 2003

본 연구에서는 천수방정식을 이용하여 2차원 개수로에서 하상과 하폭이 급격히 변화하는 경우에 발생하는 천이류를 해석하였다. 불연속점 근처에서 발생하는 수치진동을 제어하면서 시간과 공간에 대한 2차 정확도를 확보하기 위하여 WAF 기법에 TVD 조건을 갖는 흐름률 제한자를 도입하였으며, Riemann해를 계산하기 위하여 3개의 전파속도를 고려하는 HLLC 방법을 이용하였다. 개수로에서 단면변화를 고려한 2차원 해석을 할 경우, 격자 구성과 경계 처리에서 어려움이 발생한다. 이와 같은 어려움을 해결하기 위하여 일반좌표계를 도입하여 하폭이 변화하는 구간에 발생하는 천이류를 수치모의하였다.

• 한국항공우주학회지
• /
• 제32권9호
• /
• pp.12-26
• /
• 2004

수치점성을 효과적으로 줄이기 위한 새로운 AUSM-type 수치기법을 개발하였다. TVD 제한자 분석을 통하여 제어면에서의 불성치를 보다 정확하게 예측할 수 있는 기준을 찾아내고 이를 이용하여 제어면의 물성치를 적절히 정의하는 M-AUSMPW+를 개발하였다. M-AUSMPW+의 장점은 다차원 유동을 해석하는데 있어 분명히 나타난다. M-AUSMPW+ 는 유동과 격자계가 일치하지 않는 격자계에서 수치점성을 효과적으로 제거할 수 있기 때문에 다차원 유동을 계산하는데 있어 기존의 어떠한 수치기법보다 정확하고 효과적인 계산이 가능하다. 접촉불연속면, 와류 유동, 충격파 경계층 상호작용, 정성 충격파관 문제를 통해 이를 확인하였다.

• 한국전산유체공학회지
• /
• 제12권2호
• /
• pp.53-62
• /
• 2007

A high resolution numerical method aimed at solving cavitating flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media at isothermal condition and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• 한국해안·해양공학회논문집
• /
• 제36권1호
• /
• pp.11-19
• /
• 2024

본 연구에서는 완전비선형 Boussinesq 방정식을 사용한 FUNWAVE-TVD 모델에 경사구조물의 월파 계산이 가능하도록 EurOtop(2018)의 경험식을 추가하여 프로그램을 수정하였다. 수정된 수치 모형은 경사제에 대한 CLASH data 및 수리모형실험 자료와 비교하여 그 타당성이 검증되었다. 본 모형은 피복 블록의 조도계수 차이에 따른 월파량의 변화를 정확히 재현하였고, 상대여유고 증가에 따른 월파량의 감소율을 잘 재현하였다. 경사제의 월파량은 피복 블록의 배치 상태에 따라서 큰 차이를 보였다. Tetrapods를 정적(regular positioning)으로 배치한 경우는 난적(random positioning)으로 배치한 경우보다 월파량이 크게 증가하였고 Rocks로 배치했을 때와 일치하였다. 한편, Rocks를 1열로 배치한 경우는 2열로 배치한 경우보다 월파량이 커졌으며 이는 구조물 표면의 조도와 투수성의 영향 때문으로 판단된다.

• 한국전산유체공학회지
• /
• 제19권3호
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.49-54
• /
• 1998

3-dimensional Euler solver is parallelized. The spatial discretization method is the 2nd order TVD scheme and DADI method with multigrid is used as a time integration. In order to parallelize this solver, the domain decomposition method with overlapped grid and message passing techniques are used. The informations on the each inter-processor bound-aries are communicated with MPI library. Finally, the parallel performance repsented by calculating the ONERA M6 wing at transonic flow condition using CRAY T3E and C90.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 춘계 학술대회논문집
• /
• pp.128-134
• /
• 2005

In the CFD area, the numerical analysis of high Mach number flow over a blunt-body poses many issues. Various numerical schemes have been developed to cover the issues, but the traditional schemes are still used widely due to the complexities of new schemes and intricacy of modifying the established codes. In the present study, the well-known Roe's FDS based on TVD-MUSCL scheme is used for the solution of very high Mach number three-dimensional flows posing carbuncle and non-physical phenomena in numerical analysis. A parametric study was carried out to account for the effects of the entropy fixing, grid configurations and initial condition. The carbuncle phenomena could be easily overcome by the entropy fixing, and the non-physical solution could be eliminated by the use of the modified initial condition regardless of entropy fixing and grid configurations.