• 천문학회보
• /
• 제39권2호
• /
• pp.119.2-119.2
• /
• 2014

Many astrophysical phenomena involve processes of magnetohydrodynamics (MHD) and relativistic magnetohydrodynamics (RMHD). A number of numerical schemes have been developed to solve the equations of ideal MHD and RMHD. Recent codes are based on upwind schemes which solve hyperbolic systems of equations following the characteristics of the systems. Upwind schemes stand out by their robustness, clarity of the underlying physical model, and ability of achieving high resolution. We present MHD and RMHD codes based on the total variation diminishing (TVD) and weighted essentially non-oscillatory (WENO) schemes, which are second and higher order accurate extensions of upwind schemes. We demonstrate the ability and limitation of codes based on upwind schemes through a series of tests.

• 대한기계학회논문집B
• /
• 제27권7호
• /
• pp.995-1006
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes for large eddy simulation is evaluated by a dynamic analysis. Large eddy simulation of isotropic turbulence is performed with various dissipative and non-dissipative schemes to investigate the effect of numerical dissipation on the resolved solutions. It is shown by the present dynamic analysis that upwind schemes reduce the aliasing error and increase the finite differencing error. The existence of optimal upwind scheme that minimizes total numerical error is verified. It is also shown that the finite differencing error from numerical dissipation is the leading source of numerical errors by upwind schemes. Simulations of a turbulent channel flow are conducted to show the existence of the optimal upwind scheme.

• 대한기계학회논문집B
• /
• 제27권7호
• /
• pp.984-994
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes for large eddy simulation is evaluated by a spectral, static error analysis. To investigate the effect of numerical dissipation on LES solutions, power spectra of discretization errors are evaluated for isotropic turbulence models in both continuous and discrete wavevector spaces. Contrary to the common belief, the aliasing errors from upwind-biased schemes are larger than those from comparable non-dissipative schemes. However, this result is the direct consequence of the definition of the power spectral density of the aliasing error, which poses the limitation of the static error analysis for upwind schemes.

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

본 연구에서는 다중 격지 Navier-Stokes 방정식 해석의 수렴성 향상을 위하여 상류 차분 기법들의 주파수 영역에서의 특생 해석을 수행하였다. 1차 상류 차분 기법에 기반한 내재적 연산자의 예조건화 특성 향상을 위하여 다차원적인 효과를 갖는 2차 상류 차분 기법이 전통적인 2차 상류 차분 기법에 비하여 우수한 예조건화 특성을 가지는 이유를 제시하였다. 주파수 영역에서의 해석 결과에 대한 검증을 위하여 완전 성김 다중 격자 기법과 예 조건화된 다단계 시간 전진 기법을 적용하였다. 비 점성 유동장 및 Spalart-Allmaras 난류 모델을 이용한 난류 유동장 해석을 수행하였으며, 주파수 영역해석의 결과와 일치하는 우수한 수렴 특성을 가짐을 확인하였다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.807-810
• /
• 2002

The suitability of high-order accurate, central and upwind-biased compact difference schemes is evaluated for the large-eddy simulations of flows in complex geometry. Two flow geometries are considered: channel and circular cylinder. The effects of numerical dissipation and aliasing error on the evaluation of subgrid scale stress are investigated by extending the analysis by Ghosal (1) to centered and upwind compact schemes. It is shown that the failure of upwind schemes mainly comes from the aliasing error.

• 한국수자원학회논문집
• /
• 제38권2호
• /
• pp.155-166
• /
• 2005

상류이송기법은 불연속 흐름을 해석할 수 있기 때문에 댐붕괴류, 천이류 등의 해석에 많이 이용되고 있다. 그러나 상류이송기법은 생성항 처리과정에서 발생하는 오차로 인해 불균일한 단면을 가진 자연하천에는 거의 적용되지 못하고 단순화된 하도에만 주로 적용되어 왔다. 본 논문에서는 생성항의 차분화를 위해서 정규화된 Jacobian을 사용하는 상류이송형 생성항 처리기법을 개발하였다. 적용 결과 본 연구에서 제안된 생성항 처리기법이 정확하면서 효율적인 것으로 나타났다. 본 연구에서 제안한 방법은 단순한 형태를 지니고 있으며 다른 상류이송기법에도 다양하게 적용될 수 있을 것으로 판단된다.

• 대한기계학회논문집B
• /
• 제24권9호
• /
• pp.1227-1238
• /
• 2000

This study examines the effects of the discretization schemes on numerical solutions of viscoelastic fluid flows. For this purpose, a temporally evolving mixing layer, a two-dimensional vortex pair interacting with a wall, and a turbulent channel flow are selected as the test cases. We adopt a fourth-order compact scheme (COM4) for polymeric stress derivatives in the momentum equations. For convective derivatives in the constitutive equations, the first-order upwind difference scheme (UD) and artificial diffusion scheme (AD), which are commonly used in the literature, show most stable and smooth solutions even for highly extensional flows. However, the stress fields are smeared too much and the flow fields are quite different from those obtained by higher-order upwind difference schemes for the same flow parameters. Among higher-order upwind difference schemes, a third-order compact upwind difference scheme (CUD3) shows most stable and accurate solutions. Therefore, a combination of CUD3 for the convective derivatives in the constitutive equations and COM4 for the polymeric stress derivatives in the momentum equations is recommended to be used for numerical simulation of highly extensional flows.

• 대한기계학회논문집B
• /
• 제27권7호
• /
• pp.973-983
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes is evaluated for large eddy simulation of turbulent flow. Two turbulent flows are considered: turbulent channel flow at Re = 23000 and flow over a circular cylinder at Re = 3900. The effects of numerical dissipation on the finite differencing and aliasing errors and the subgrid-scale stress are investigated. It is shown through the simulations that compact upwind schemes are not suitable for LES, whereas the fourth order-compact centered scheme is a good candidate for LES provided that proper dealiasing of nonlinear terms is performed. The classical issue on the aliasing error and the treatment of nonlinear terms is revisited with compact difference schemes.

• 한국전산유체공학회지
• /
• 제10권1호
• /
• pp.99-106
• /
• 2005

Performance of first, second and third order accurate methods for calculation of in viscid fluxes in fluid flow governing equations are investigated here. For the purpose, an upwind method based on Roe's scheme is used to solve 2-dimensional Euler equations. To increase the accuracy of the method two different schemes are applied. The first one is a second and third order upwind-based algorithm with the MUSCL extrapolation Van Leer (1979), based on primitive variables. The other one is an upwind-based algorithm with the Chakravarthy extrapolation to the fluxes of mass, momentum and energy. The results show that the thickness of shock layer in the third order accuracy is less than its value in second order. Moreover, applying limiter eliminates the oscillations near the shock while increases the thickness of shock layer especially in MUSCL method using Van Albada limiter.

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

Studies of the numerical characteristics of implicit upwind schemes, such as upwind ADI, Line Gauss-Seidel(LGS) and Point Gauss-Seidel(LU) algorithms, for preconditioned Navier-Stokes equations ate performed. All the algorithms are expressed in approximate factorization form and Von Neumann stability analysis and convergence studies are made. Preconditioning is applied for efficient convergence at low Mach numbers and low Reynolds numbers. For high aspect ratio computations, the ADI and LGS algorithms show efficient and uniform convergence up to moderate aspect ratio if we adopt viscous preconditioning based on min- CFL/max- VNN time-step definition. The LU algorithm, on the other hand, shows serious deterioration in convergence rate as the grid aspect ratio increases. Computations for practical applications also verify these results.