• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.2
• /
• pp.93-114
• /
• 2019

In this article, we introduce a finite difference method for solving the Navier-Stokes equations in rectangular domains. The method is proved to be energy stable and shown to be second-order accurate in several benchmark problems. Due to the guaranteed stability and the second order accuracy, the method can be a reliable tool in real-time simulations and physics-based animations with very dynamic fluid motion. We first discuss a simple convection equation, on which many standard explicit methods fail to be energy stable. Our method is an implicit Runge-Kutta method that preserves the energy for inviscid fluid and does not increase the energy for viscous fluid. Integration-by-parts in space is essential to achieve the energy stability, and we could achieve the integration-by-parts in discrete level by using the Marker-And-Cell configuration and central finite differences. The method, which is implicit and second-order accurate, extends our previous method [1] that was explicit and first-order accurate. It satisfies the energy stability and assumes rectangular domains. We acknowledge that the assumption on domains is restrictive, but the method is one of the few methods that are fully stable and second-order accurate.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.11 no.6
• /
• pp.18-25
• /
• 2007

In a plastic metal forming of thermally rate-sensitive material, the localized shear band stems from evolution of a narrow region in which intensive plastic flow occurs. And it give rise to fatal fracture with plastic instability. The objectives of this study are to investigate the localization behavior by using numerical method and predict the failure for WHA(Tungsten Heavy Alloy). In this work, the implicit finite difference scheme is used because of the advantage about convergence and the numerical stability. This study is based on an analysed material with hardening as well as thermally softening behavior which includes isotropic strain hardening and observed the extension of localization within shear band according to material properties.

• Proceedings of the IEEK Conference
• /
• 2002.07a
• /
• pp.137-140
• /
• 2002

In this paper, a new algorithm for two-dimensional (2-D) finite-difference time-domain (FDTD) method is presented. By this new method, the maximum time step size can be increased over the Courant-Friedrich-Levy (CFL) condition restraint. This new algorithm is adapted from an Alternating-Direction Implicit FDTD (ADI-FDTD) method. However, unlike the ADI-FDTD algorithm. the alternation is performed with respect to the blocks of fields rather than with respect to each respective coordinate direction. Moreover. this method can be efficiently simulated with parallel computation. and it is more efficient than the conventional FDTD method in terms of CPU time. Numerical formulations are shown and simulation results are presented to demonstrate the effectiveness and efficiency of our proposed method.

• Geomechanics and Engineering
• /
• v.6 no.6
• /
• pp.545-560
• /
• 2014

For slope stability analysis, an alternative to the classical limit equilibrium method (LEM) of slices is the shear strength reduction method (SRM), which can be integrated into finite element analysis or finite difference analysis. Recently, probabilistic analysis of earth slopes has been very attractive because it is capable to take the soil uncertainty into account. However, the SRM is less commonly extended to probabilistic framework compared to a variety of probabilistic LEM analysis of earth slopes. To overcome some limitations that hinder the development of probabilistic SRM stability analysis, a new procedure based on recursive algorithm FORM with sensitivity analysis in the space of original variables is proposed. It can be used to deal with correlated non-normal variables subjected to implicit limit state surface. Using the proposed approach, a probabilistic finite element analysis of the stability of an existing earth dam is carried out in this paper.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2000.03b
• /
• pp.631-638
• /
• 2000

Pioneering work by Terzaghi imparted scientific and mathematical bases to many aspects of this subject and many people use this theory to measure the consolidation settlement until now. In this paper, Finite Difference Methods for consolidation are considered. First, it is shown the stability criterion of Explicit scheme and the Crank-Nicolson scheme, although unconditionally stable in the mathematical sense, produces physically unrealistic solutions when the time step is large. it is also shown that The Fully Implicit scheme shows more satisfactory behavior, but is less accurate for small time steps. and then we need to decide what scheme is more proper to consolidation. The purpose of this paper is to suggest the pertinent scheme to consolidation.

• Structural Engineering and Mechanics
• /
• v.48 no.5
• /
• pp.711-736
• /
• 2013

In this paper new implicit time integration called N-IHOA is presented for dynamic analysis of high damping systems. Here, current displacement and velocity are assumed to be functions of the velocities and accelerations of several previous time steps, respectively. This definition causes that only one set of weighted factors is calculated from the Taylor series expansion which leads to a simple approach and reduce the computational efforts. Moreover a comprehensive study on stability of the proposed method i.e., N-IHOA compared with IHOA integration which is performed based on amplification matrices proves the ability of the N-IHOA in high damping vibrations such as control systems. Also, wide range of numerical examples which contains single/multi degrees of freedom, damped/un-damped, free/forced vibrations from finite element/finite difference demonstrate that the accuracy and efficiency of the proposed time integration is more than the common approaches such as the IHOA, the Wilson-${\theta}$ and the Newmark-${\beta}$.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.1
• /
• pp.58-69
• /
• 1995

The present paper is devoted to constructing a numerical algorithm for solving un steady problems on generation, propagation and interaction of nonlinear waves at a surface of ideal fluid, within the framework of the potential-flow model. The numerical scheme is implicit. with non-linearity iteration at every step of time. the finite-difference method with boundary-fitted coordinates are presented in favor for validity and high efficiency of the numerical model developed. Among these arguments, there are the results of calculations of two test problems-on stretching of a liquid ellipse and on wave generation by lifting a portion of a bottom.

• 한국방재학회:학술대회논문집
• /
• 2007.02a
• /
• pp.181-185
• /
• 2007

A free-surface correction(FSC) method is presented to solve the 3-D shallow water equations. Using the mode splitting process, FSC method can simulate shallow water flows under the hydrostatic assumption. For the hydrostatic pressure calculation, the momentum equations are firstly discretized using a semi-implicit scheme over the vertical direction leading to the tri-diagonal matrix systems. A semi-implicit scheme has been adopted to reduce the numerical instability caused by relatively small vertical length scale compare to horizontal one. and, as the free surface correction step the final horizontal velocity fields are corrected after the final surface elevations are obtained. Finally, the vertical final velocity fields can be calculated from the continuity equation. The numerical model is applied to the calculation of the simulation of flow fields in a rectangular open channel with the tidal influence. The comparisons with the analytical solutions show overall good agreements between the numerical results and analytical solutions.

• East Asian mathematical journal
• /
• v.40 no.1
• /
• pp.109-118
• /
• 2024

This paper proposes a numerical method based on domain decomposition to find approximate solutions for one-dimensional convection-diffusion equations with Neumann boundary conditions. First, the equations are transformed into convection-diffusion equations with Dirichlet conditions. Second, the author introduces the Prediction/Correction Domain Decomposition (PCDD) method and estimates errors for the interface prediction scheme, interior scheme, and correction scheme using known error estimations. Finally, the author compares the PCDD algorithm with the fully explicit scheme (FES) and the fully implicit scheme (FIS) using three examples. In comparison to FES and FIS, the proposed PCDD algorithm demonstrates good results.

• Proceedings of the Korean Geotechical Society Conference
• /
• 1999.03a
• /
• pp.349-356
• /
• 1999

In general, the term soft ground includes clayey soils, which have large compressibility and small shear resistance due to the external load. All process of consolidation in compressible soils can be explained in terms of a transfer of load from an incompressible pore-water to a compressible soil structure. Therefore, one of the most important subjects about the characteristics of the time-dependent consolidation of the clay foundation by the change of load may be the presumption of the final settlement caused by consolidation and the degree of consolidation according to the time. The problems of discontinuous layer interface are very important in the algorithm and programming for the analysis of multi-layered soils using a numerical analysis, finite difference method. Better results can be obtained by the Process for discontinuous layer interface, since it can help consolidation analysis to model the actual ground. The purpose of this paper Provides an efficient computer algorithm based on numerical analysis using finite difference method(F.D.M.) which account for multi-layered soils to determine the degree of consolidation and excess pore pressures relative to time and positions more realistically.