• Bulletin of the Korean Mathematical Society
• /
• v.45 no.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.

• Journal of computational fluids engineering
• /
• v.13 no.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.10a
• /
• 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
• /
• v.27 no.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.

• Communications of the Korean Mathematical Society
• /
• v.24 no.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
• /
• v.12 no.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.

• Journal of Korea Water Resources Association
• /
• v.36 no.1
• /
• pp.1-11
• /
• 2003

This paper describes a finite volume, two-dimensional model. It adopts a recently developed essentially non-oscillatory(ENO) schemes based on the Lax-Friedrichs solver, which was modified for a finite volume grid, and employs a modified MUSCL(Monotonic Upstream centered Scheme for Conservation Law) for second-order accuracy in space. To demonstrate the applications of the model, it is applied to solve the 1-D and 2-D dam-break problems. The model in conjunction with the modified MUSCL showed a better agreement with analytical solutions than the minmod function in 1-D dam-break problems and is satisfactorily validated with documented published data in 2-D dam-break problems. The model was applied to tidal wane entering channel at one end, and the results showed a good agreement with analytical solutions. In the channel with reflective boundary conditions specified at the extremities, the model was capable of accurately simulating the wave propagation.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.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.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• 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.05a
• /
• 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.