• 대한수학회보
• /
• 제45권4호
• /
• pp.781-795
• /
• 2008

In this paper, we describe the application procedure of the adaptive mesh refinement (AMR) for the weighted essentially non-oscillatory schemes (WENO), and observe the effects of the derived algorithm when problems have piecewise smooth solutions containing discontinuities. We find numerically that the dissipation of the WENO scheme can be lessened by the implementation of AMR while the accuracy is maintained. We deduce from the experiments that the AMR-implemented WENO scheme captures shocks more efficiently than the WENO method using uniform grids.

• 한국전산유체공학회지
• /
• 제13권3호
• /
• pp.35-43
• /
• 2008

A methodology for the simulation of compressible high Reynolds number flow over rigid and moving bodies on a structured Cartesian grid is described in this paper. The approach is based on a modified version of the Brinkman Penalization method. To avoid oscillations in the vicinity of the body and to simulate shcok-containing flows, a Weighted Essentially Non-Oscillatory scheme is used to discretize the spatial flux derivatives. For high Reynolds number viscous flow, two turbulence models of the two-equation Menter's SST URANS model and a two-equation Detached Eddy Simulation are implemented. Some simple flow examples are given to assess the accuracy of the technique. Finally, a moving grid capability is demonstrated.

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

A methodology for the simulation of compressible high Reynolds number flow over rigid and moving bodies on a structured Cartesian grid is described in this paper. The approach is based on a modified version of the Brinkman Penalization method. To avoid oscillations in the vicinity of the body and to simulate shcok-containing flows, a Weighted Essentially Non-Oscillatory scheme is used to discretize the spatial flux derivatives. For high Reynolds number viscous flow, two turbulence models of the two-equation Menter's SST URANS model and a two-equation Detached Eddy Simulation are implemented. Some simple flow examples are given to assess the accuracy of the technique. Finally, a moving grid capability is demonstrated.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제27권4호
• /
• pp.207-231
• /
• 2023

This paper aims to improve the alternative formulation of the fifth- and sixth-order accurate weighted essentially non-oscillatory (AWENO) finite difference schemes. The first is to derive the AWENO scheme with sixth-order accuracy in the smooth region of the solution. Second, a new weighted polynomial functions combining the perturbed forms with conserved variable to the AWENO is constructed; the new form of tunable functions are invented to maintain non-oscillatory property. Detailed numerical experiments are presented to illustrate the behavior of the new perturbational AWENO schemes. The performance of the present scheme is evaluated in terms of accuracy and resolution of discontinuities using a variety of one and two-dimensional test cases. We show that the resulted perturbational AWENO schemes can achieve fifth- and sixth-order accuracy in smooth regions while reducing numerical dissipation significantly near singularities.

• 대한수학회논문집
• /
• 제24권3호
• /
• pp.459-466
• /
• 2009

In this paper, a numerical scheme is introduced to solve the Cauchy problem for one-dimensional hyperbolic equations. The mesh points of the proposed scheme are distributed along characteristics so that the solution on the stencil can be easily and accurately computed. This is very important in reducing errors of the scheme because many numerical errors are generated when the solution is estimated over grid points. In addition, when characteristics intersect, the proposed scheme combines corresponding grid points into one and assigns new characteristic to the point in order to improve computational efficiency. Numerical experiments on the inviscid Burgers' equation have been presented.

• Ocean Systems Engineering
• /
• 제12권3호
• /
• pp.359-384
• /
• 2022

We develop in this work a new well-balanced preserving-positivity path-conservative central-upwind scheme for Saint-Venant-Exner (SVE) model. The SVE system (SVEs) under some considerations, is a nonconservative hyperbolic system of nonlinear partial differential equations. This model is widely used in coastal engineering to simulate the interaction of fluid flow with sediment beds. It is well known that SVEs requires a robust treatment of nonconservative terms. Some efficient numerical schemes have been proposed to overcome the difficulties related to these terms. However, the main drawbacks of these schemes are what follows: (i) Lack of robustness, (ii) Generation of non-physical diffusions, (iii) Presence of instabilities within numerical solutions. This collection of drawbacks weakens the efficiency of most numerical methods proposed in the literature. To overcome these drawbacks a reformulation of the central-upwind scheme for SVEs (CU-SVEs for short) in a path-conservative version is presented in this work. We first develop a finite-volume method of the first order and then extend it to the second order via the averaging essentially non oscillatory (AENO) framework. Our numerical approach is shown to be well-balanced positivity-preserving and shock-capturing. The resulting scheme could be seen as a predictor-corrector method. The accuracy and robustness of the proposed scheme are assessed through a carefully selected suite of tests.

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

본 연구에서는 2차원 부정류 해석을 위한 유한체적 수치모형을 개발하였다. 개발된 모형에서는 천수방정식의 수치해를 구하기 위하여 Lax-Friedrichs 해법에 근간을 둔 ENO (essentially non-oscillatory) 기법을 적용하였으며, 공간적으로 2차 정도 모의를 위하여 수정 MUSCL 기법을 적용하였다. 개발된 모형의 적용성을 평가하기 위하여 해석해가 있는 1차원 댐 붕괴파 모의에 적용한 결과, 해석해와 유사한 결과를 나타냈으며, 수정 MUSCL 기법의 모의결과를 minmod 경사제한자를 사용한 모의결과와 비교한 결과, 수정 MUSCL 기법을 사용하는 경우 해석해와 더 근접한 결과를 나타냈다. 댐붕괴파의 2차원 해석에 대한 모형의 적용성을 평가하기 위하여 대칭형, 비대칭형 댐 붕괴에 적용한 결과, 저수지에서 수로 구간으로 충격파 전파를 잘 모의했으며, 기존의 연구결과와 유사한 모의결과를 나타냈다. 개수로의 부정류 해석에 대한 모형의 적용성을 평가하기 위하여 말단이 폐쇄된 수로의 입구에서 정상파가 주어졌을 때 수위와 유속을 모의하여 해석해와 비교하였으며, 입구와 말단이 폐쇄된 수로에서 파의 진동을 모의하여 모형의 거동을 평가하였다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제20권3호
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.

• 한국전산유체공학회지
• /
• 제11권4호
• /
• pp.39-47
• /
• 2006

In this paper, the authors developed a multidimensional interpolation method inside a finite volume cell in the computation of high-order accurate numerical flux such as the fifth order WEND (weighted essentially non-oscillatory) scheme. This numerical method starts from a simple Taylor series expansion in a proper spatial order of accuracy, and the WEND filter is used for the reconstruction of sharp nonlinear waves like shocks in the compressible flow. Two kinds of interpolations are developed: one is for the cell-averaged values of conservative variables divided in one mother cell (Type 1), and the other is for the vertex values in the individual cells (Type 2). The result of the present study can be directly used to the cell refinement as well as the convective flux between finer and coarser cells in the Cartesian adaptive grid system (Type 1) and to the post-processing as well as the viscous flux in the Navier-Stokes equations on any types of structured and unstructured grids (Type 2).

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

The BGK equation, which is the kinetic model equation of Boltzmann equation, is solved using FDDO(finite difference with the discrete-ordinate method) to compute the rarefied flow of monatomic gas. Using reduced velocity distribution and discrete ordinate method, the scalar equation is transformed into a system of hyperbolic equations. High resolution ENO(Essentially Non-Oscillatory) scheme based on Harten-Yee's MFA(Modified Flux Approach) method with Strang-type explicit time integration is applied to solve the system equations. The calculated results are well compared with the experimental density field of NACA0012 airfoil, validating the developed computer code. Next. the computed results of circular cylinder flow for various Knudsen numbers are compared with the DSMC(Direct Simulation Monte Carlo) results by Vogenitz et al. The present scheme is found to be useful and efficient far the analysis of two-dimensional rarefied gas flows, especially in the transitional flow regime, when compared with the DSMC method.