• 한국수자원학회논문집
• /
• 제43권1호
• /
• pp.105-117
• /
• 2010

본 논문은 2차원 천수방정식을 해석하기 위해 새롭고 간단한 MUSCL 재구성법을 제안하였다. 수정 MUSCL 기법은 보존변수의 재구성을 위해 계산격자와 인접격자의 보존변수 차에 대해 각 경계면에 균일하게 분배하는 기존 방법 대신 면적가중비를 사용하였으며, 이 방법은 정형격자 뿐만 아니라 비정형 격자의 사용에도 보존변수의 물리적 재구성이 가능하다. 또한, 본 연구에서는 비구조적 격자의 적용이 가능한 차원비분리 기법을 적용하였으며, 수정 MUSCL 기법의 사용으로 발생할 수 있는 수치진동을 제어하기 위해 TVD 기법의 경사제한자를 사용하였다. 하상경사항의 정확하고 효율적인 수치 처리를 위해 수정 MUSCL 기법을 수면경사법과 연계하였다. 제안한 기법을 적용한 유한체적모형을 건물의 영향을 고려한 댐 붕괴 해석 및 Bellos의 댐 붕괴 실험에 적용하고, 적용결과를 실험실 자료 및 기존 연구자의 계산결과와 비교하여 개발모형을 검증하였다.

• 한국전산유체공학회지
• /
• 제4권2호
• /
• pp.57-66
• /
• 1999

Numerical Analysis has been done for the supersonic off-design jet flow due to the pressure difference between the jet and the ambient fluid. The difference of pressure generates an oblique shock or an expansion wave at the nozzle exit. The waves reflect repeatedly on the center axis and the sonic surface in the shear layer. The pressure difference is resolved across these reflected waves. In this paper, the axi-symmetric Navier-Stokes equation has been used with the κ-ε turbulence model. The second order TVD scheme with flux limiters, based on the flux vector split with the smooth eigenvalue split, has been used to capture internal shocks and other discontinuities. Numerical calculations have been done to analyze the off-design jet flow due to the pressure difference. The variation of pressure along the flow axis is compared with an experimental result and other numerical result. The characteristics of the interaction between the shock cell and the turbulence mixing layer have been analyzed.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.723-728
• /
• 2008

We performed numerical analysis of base drag phenomenon, when a projectile with backward step flies into atmosphere at supersonic speed. We compared with other researchers. From our previous studies that were 2-dimensional simulation, we found out from sophisticated simulations that need dense mesh points to compare base pressure and velocity profile after from base with experimental data. Therefore, we focus on high order spatial disceretization over 3rd order with TVD such as MUSCL TVD 3rd, 5th, and WENO 5th order, and Limiters such as minmod, Triad. Moreover, we enforce to flux averaging schemes such as Roe, RoeM, HLLE, AUSMDV. In present, one dimensional result of Euler tests, there are Sod, Lax, Shu-Osher and interacting blast wave problems. AUSMDV as a flux averaging scheme with MUSCL TVD 5th order as spatial resolution is good agreement with exact solutions than other combinations. We are carrying out the same approaches into 3-dimensional base flow only candidate flux schemes that are Roe, AUSMDV. Additionally, turbulence models are used in 3-dimensional flow, one is Menter s SST DES model and another is Sparlat-Allmaras DES/DDES model in Navier-Stokes equations.

• 한국전산유체공학회지
• /
• 제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.

• 천문학회보
• /
• 제37권1호
• /
• pp.66.2-66.2
• /
• 2012

Building a relativistic magnetohydrodynamic (RMHD) code based on upwind scheme is a challenging project, because eigenvalues and eigenvectors are not yet analytically given. Here, we present analytic expressions for eigenvalues and eigenvectors in isothermal flows. And then we show tests performed with a code based on the total variation diminishing (TVD) scheme.

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

In this paper, the preconditioned multistage time stepping methods which are popular multigrid smoothers is implemented for the compressible Navier-Stokes calculation with full-coarsening multigrid method. The convergence characteristic of the point-Jacobi and Alternating direction line Jacobi(DDADI) preconditioners are studied. The performance of 2nd order upwind numerical fluxes such as 2nd order upwind TVD scheme and MUSCL-type linear reconstruction scheme are compared in the inviscid and viscous turbulent flow caculations.

• 천문학회보
• /
• 제37권2호
• /
• pp.101.1-101.1
• /
• 2012

Building a relativistic magnetohydrodynamic (RMHD) code based on upwind schemes has been a challenging project, because of the absence of analytic expressions of eigenvalues and eigenvectors. We found analytic expressions of eigenvalues and eigenvectors for adiabatic RMHD flows which are relatively simple and manageable. Especially, our eigenvectors can handle all degenerate points. Using these analytic forms, we built a code based on the total variation diminishing (TVD) scheme, and successfully performed one-dimensional shock tube tests.

• 한국추진공학회지
• /
• 제10권2호
• /
• pp.1-14
• /
• 2006

본 논문은 약한 불안정 데토네이션 영역부터 강한 불안정 데토네이션 영역까지 여러 영역에 걸친 데토네이션 파 셀 구조 모사에 대한 수치적 문제점들을 살펴보았다. 비열 비 값이 변하는 점성 유체 역학 방정식 및 1단계 Arrhenius 반응 모델 해석을 위하여 MUSCL-type TVD 기법을 이용한 공간 차분과 4차 정확도의 Runge-Kutta 시간 적분을 이용하였다. 일련의 수치해석 연구는 여러 반응 상수 및 격자 해상도에 따른 데토네이션 셀 구조를 해석하기 위하여 요구되는 계산 조건을 구하기 위하여 다양한 데토네이션 현상 영역에서 수행되었다. 다른 영역의 데토네이션 현상에서 셀 구조를 포착하기 위한 계산 영역의 크기와 최소 격자 해상도를 찾아내기 위하여 정상 1차원 ZND 해석 결과와 전산 해석 결과를 비교 검토하였다.

• 한국해안·해양공학회논문집
• /
• 제32권6호
• /
• pp.569-574
• /
• 2020

천수방정식에 대한 불연속 갤러킨 기법 (DG)은 주로 양해법 기반으로 개발되어 적용되어 왔으나, 바닥마찰항의 처리, 과도한 CFL 조건 등의 불리한 점이 지적되어 왔다. 이에 대한 대안으로써, 본 연구에서는 음해법 기반의 모형을 개발하고 이를 적용하여 향후 가능성을 입증하였다. 본 논문에서 연구한 사례에서는 선형 삼각형 요소를 사용하였고, 수치흐름률로서 Roe 흐름률을 이용하였으며, TVD 특성 보존을 위한 기울기 제한자를 적용하였다. 적용 사례로서 실린더 주변의 흐름과 댐 붕괴류 문제 등에 대하여 적용하고, 기존의 실험치, 수치해와 비교하여 잘 일치함을 확인하였다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제15권4호
• /
• pp.291-306
• /
• 2011

We consider the second order Strang splitting method to approximate the solution to the KdV equation. The model equation is split into three sets of initial value problems containing convection and dispersal terms separately. TVD MUSCL or MUSCL scheme is applied to approximate the convection term and the second order centered difference method to approximate the dispersal term. In time stepping, explicit third order Runge-Kutta method is used to the equation containing convection term and implicit Crank-Nicolson method to the equation containing dispersal term to reduce the CFL restriction. Several numerical examples of weakly and strongly dispersive problems, which produce solitons or dispersive shock waves, or may show instabilities of the solution, are presented.