• 대한수학회보
• /
• 제54권1호
• /
• pp.145-158
• /
• 2017

We present a finite difference method for solving the Ohta-Kawasaki model, representing a model of mesoscopic phase separation for the block copolymer. The numerical methods for solving the Ohta-Kawasaki model need to inherit the mass conservation and energy dissipation properties. We prove these characteristic properties and solvability and unconditionally gradient stability of the scheme by using Hessian matrices of a discrete functional. We present numerical results that validate the mass conservation, and energy dissipation, and unconditional stability of the method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제23권1호
• /
• pp.1-18
• /
• 2019

In contrast to the well-developed convex splitting schemes for gradient flows of two-component system, there were few efforts on applying the convex splitting idea to gradient flows of multi-component system, such as the vector-valued Cahn-Hilliard (vCH) equation. In the case of the vCH equation, one need to consider not only the convex splitting idea but also a specific method to manage the partition of unity constraint to design an unconditionally energy stable scheme. In this paper, we propose a constrained Convex Splitting (cCS) scheme for the vCH equation, which is based on a convex splitting of the energy functional for the vCH equation under the constraint. We show analytically that the cCS scheme is mass conserving and unconditionally uniquely solvable. And it satisfies the constraint at the next time level for any time step thus is unconditionally energy stable. Numerical experiments are presented demonstrating the accuracy, energy stability, and efficiency of the proposed cCS scheme.

• Journal of electromagnetic engineering and science
• /
• 제3권1호
• /
• pp.39-44
• /
• 2003

An unconditionally stable compact 2D Alternating-Direction Implicit (ADI) FDTD method for calculating dispersion characteristics of waveguide structures is proposed. The numerical stability and numerical dispersion relation of the proposed method are also presented and discussed. Numerical wavelengths for the dominant and higher order modes in a hollow waveguide are obtained from numerical simulations and compared with those from the analytical dispersion relation. The numerical results show that the proposed scheme has the potential to successfully analyze a class of waveguides having locally fine geometry with reduced numerical costs.

• Structural Engineering and Mechanics
• /
• 제64권4호
• /
• pp.495-504
• /
• 2017

This paper proposes a new direct numerical integration algorithm for solving equation of motion in structural dynamics problems with nonlinear stiffness. The new implicit method's degree of accuracy is higher than that of existing methods due to the higher order of the acceleration. Two parameters are defined, leading to a new family of unconditionally stable methods, which helps to take greater time steps in integration and eliminate concerns about the duration of solving. The method developed can be utilized for a number of solid plane finite elements, examples of which are given to compare the proposed method with existing ones. The results indicate the superiority of the proposed method.

• Korean Journal of Mathematics
• /
• 제27권1호
• /
• pp.193-219
• /
• 2019

The Navier-Stokes equations with variable density are challenging problems in numerical analysis community. We recently built the 2nd order stabilized Gauge-Uzawa method [SGUM] to solve the Navier-Stokes equations with constant density and have estimated theoretically optimal accuracy. Also we proved that SGUM is unconditionally stable. In this paper, we apply SGUM to the Navier-Stokes equations with nonconstant variable density and find out the stability condition of the algorithms. Because the condition is rather strong to apply to real problems, we consider Allen-Cahn scheme to construct unconditionally stable scheme.

• 한국초전도ㆍ저온공학회논문지
• /
• 제1권2호
• /
• pp.37-42
• /
• 1999

The stability of high-temperature superconducting current leads for cryogenic devices are investigated. By assuming full transition from superconducting state to normal state at a transition temperature, the HTS current at a transition temperature, the HTS current lead shows bifurcation phenomenon. There is a bifurcation shape-factor, HTS leads have three steady state. Below the bifurcation shape-factor, the superconducting current lead is unconditionally stable, because there exists only one steady-factor HTS current lead is conditionally stable depending on the shape and intensity of disturbance.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 2000년도 춘계 학술발표회 논문집 Proceedings of EESK Conference-Spring
• /
• pp.151-158
• /
• 2000

The use of unconditionally stable implicit time integration techniques for pseudo-dynamic tests has been recently proposed and advanced by several researchers such as Thewalt and Mahin Nakashima and Shing. The developed implicit algorithms are based on a-Method of Hugest et al. In this paper a concise summary and explanation of implicit method for Pseudo dynamic test is presented. Especially The a-C method developed by shing at al. has been in-depth evaluated for this study. Important parameters of the a-C method have been analyzed by the simulation test.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2004년도 봄 학술발표회 논문집
• /
• pp.9-16
• /
• 2004

This paper deals with the development of computational schemes for the dynamic analysis of flexible and nonlinear multibody systems. Different from the existing method, this paper introduces the quaternion algebra to develop the equation of the conservation of energy. Simultaneously, Rodrigues parameters are used to express the finite rotation for the proposed scheme. The proposed energy scheme is derived such that it provides unconditionally stable conditions for the nonlinear problems. Several examples of dynamic systems are presented which illustrate the efficiency and accuracy of the developed energy schemes.

• 한국콘크리트학회:학술대회논문집
• /
• 한국콘크리트학회 2000년도 봄 학술발표회 논문집
• /
• pp.617-622
• /
• 2000

The use of unconditionally stable implicit time integration techniques for pseudo-dynamic test has been recently proposed and advanced by several researchers such as Thewalt and Mahin, Nakashima and Shing, etc. The developed implicit algorithms are based on the $\alpha$-Method of Huges et al. In this paper, a concise summary and explanation of implicit method for Pseudo dynamic tese is presented. Especially, The $\alpha$-C method developed by shing et al. has been in-depth evaluated for this study. Important parameters of the $\alpha$-C method have been analyzed by the simulation test.

• 한국진공학회:학술대회논문집
• /
• 한국진공학회 2000년도 제18회 학술발표회 논문개요집
• /
• pp.194-194
• /
• 2000

There are difficulties on transport modelling on high density plasma discharge, because of severe restrictions on space grid size and time step size. We present a new unconditionally stable algorithm for fluid simulation of high density process plasma. The origin of the restriction is investigated and a new method to solve the problem is suggested, The simulation result is compared with the other methods previously developed.