• Journal of electromagnetic engineering and science
• /
• v.11 no.1
• /
• pp.11-15
• /
• 2011

Preserving high-order accuracy in high-order FDTD solutions across dielectric interfaces is very important for practical time-domain electromagnetic simulations. This paper presents a fourth-order accurate numerical boundary scheme for the planar dielectric interface to be used in the fourth-order FDTD method proposed earlier by the author. The interface scheme for the two-dimensional (2-D) transverse magnetic (TM) polarization case is derived and validated by monitoring the $L_2$ norm errors in the numerical solutions of a partially-filled cavity demonstrating its fourth-order convergence and long-time numerical stability in the presence of the planar dielectric interface.

• 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 applied mathematics & informatics
• /
• v.19 no.1_2
• /
• pp.527-535
• /
• 2005

Numerical differentiation is one of the main topics which have been studied by many researchers. If we use the forward difference scheme or the centered difference scheme, the convergence rates to the derivative are O(h) and O($h^2$), respectively. In this paper, using the Lagrange Interpolation, we construct a simple high order numerical differentiation scheme which has the convergence rate O($h^{2k}$) if we have 2k+1 equally spaced nodes. Our scheme is constructive.

• Journal of the Korean Mathematical Society
• /
• v.60 no.2
• /
• pp.273-302
• /
• 2023

In this paper, a fully discrete numerical scheme for the viscoelastic Oldroyd flow is considered with an introduced auxiliary variable. Our scheme is based on the finite element approximation for the spatial discretization and the backward Euler scheme for the time discretization. The integral term is discretized by the right trapezoidal rule. Firstly, we present the corresponding equivalent form of the considered model, and show the relationship between the origin problem and its equivalent system in finite element discretization. Secondly, unconditional stability and optimal error estimates of fully discrete numerical solutions in various norms are established. Finally, some numerical results are provided to confirm the established theoretical analysis and show the performances of the considered numerical scheme.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.9
• /
• pp.1148-1156
• /
• 2000

The QUICK scheme for convection terms is modified for unstructured finite volume method by using linear reconstruction technique and validated through the computation of two well defined laminar flows. It uses two upstream grid points and one downstream grid point in approximating the convection terms. The most upstream grid point is generated by considering both the direction of flow and local grid line. Its value is calculated from surrounding grid points by using a linear construction method. Numerical error by the modified QUICK scheme is shown to decrease about 2.5 times faster than first order upwind scheme as grid size decreases. Computations are also carried out to see effects of the skewness and irregularity of grid on numerical solution. All numerical solutions show that the modified QUICK scheme is insensitive to both the skewness and irregularity of grid in terms of the accuracy of solution.

• Journal of applied mathematics & informatics
• /
• v.36 no.3_4
• /
• pp.181-203
• /
• 2018

In this paper, we have constructed an overlapping Schwarz method for singularly perturbed second order convection-diffusion equations. The method splits the original domain into two overlapping subdomains. A hybrid difference scheme is proposed in which on the boundary layer region we use the central finite difference scheme on a uniform mesh while on the non-layer region we use the mid-point difference scheme on a uniform mesh. It is shown that the numerical approximations which converge in the maximum norm to the exact solution. When appropriate subdomains are used, the numerical approximations generated from the method are shown to be first order convergent. Furthermore it is shown that, two iterations are sufficient to achieve the expected accuracy. Numerical examples are presented to support the theoretical results. The main advantages of this method used with the proposed scheme is it reduces iteration counts very much and easily identifies in which iteration the Schwarz iterate terminates.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.5 no.1
• /
• pp.1-9
• /
• 2001

BGK schemes are developed by improving the standard BGK numerical method. Shceme 1 uses the Osher's Godunov type solution and scheme 2 are developed to overcome the problems of scheme 1.The improved schemes show many unique properties such as entropy condition, positivity condition, higher order gas evolution model, which lead to an high degree of robustness and accuracy. The scheme 2 especially overcomes the shortcomings of the scheme 1 and posseses many superior properties that cannot be found ohter numerical schemes, and is expected to apply various problems with high accuracy and robustness.

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.471-476
• /
• 2011

In this study, a modified dispersion-correction scheme describing tsunami propagation on variable water depths is proposed by introducing additional terms to the previous numerical scheme. The governing equations used in previous tsunami propagation models are slightly modified to consider the effects of a bottom slope. The numerical dispersion of the proposed model replaces the physical dispersion of the governing equations. Then, the modified scheme is employed to simulate tsunami propagation on variable water depths and numerical results are compared with those of the previous tsunami propagation model.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.8-18
• /
• 1999

In devising a numerical approximation for the convective spatial transport of a fluid mechanical quantity, it is noted that the convective motion of a scalar quantity occurs in one-way, or from upstream to downstream. This consideration leads to a new scheme termed a pure upwind difference scheme (PUDS) in which an estimated value for a fluid mechanical quantity at a control surface is not influenced from downstream values. The formal accuracy of the proposed scheme is third order accurate. Two typical benchmark problems of a wall-driven fluid flow in a square cavity and a buoyancy-driven natural convection in a tall cavity are computed to evaluate performance of the proposed method. for comparison, the widely used simple upwind scheme, power-law scheme, and QUICK methods are also considered. Computation results are encouraging: the proposed PUDS sensitized to the convection direction produces the least numerical diffusion among tested convection schemes, and, notable improvements in representing recirculation of fluid stream and spatial change of a scalar. Although the formal accuracy of PUDS and QUICK are the same, the accuracy difference of approximately a single order is observed from the revealed results.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1263-1267
• /
• 2007

A more complete two-dimensional vertical numerical model has been developed to describe the saltwater intrusion in an estuary. The model is based on the previous studies in order to obtain a better accuracy. The non-linear terms of the governing equations are analyzed and the $\sigma$-coordinate system is employed in the vertical direction with full transformation which is recently issued in several studies because numerical errors can be generated during the coordinate transformation of the diffusion term. The advection terms of the governing equations are discretized by an upwind scheme in second-order of accuracy. By employing an explicit scheme for the longitudinal direction and an implicit scheme for the vertical direction, the numerical model is free from the restriction of temporal step size caused by a relatively small grid ratio. In previous researches, some terms induced from the transformation have been intentionally excluded since they are asked the complicate discretization of the numerical model. However, the lack of these terms introduces significant errors during the numerical simulation of scalar transport problems, such as saltwater intrusion and sediment transport in an estuary. The numerical accuracy attributable to the full transformation is verified by comparing results with a previous model in a simply sloped topography. The numerical model is applied to the Han River estuary. Very reasonable agreements for salinity intrusion are observed.