• 호남수학학술지
• /
• 제42권1호
• /
• pp.17-35
• /
• 2020

In this paper we introduce and consider a new class of variational inequalities with four operators. This class is called the extended general mixed quasi variational inequality. We show that the extended general mixed quasi variational inequality is equivalent to the fixed point problem. We use this alternative equivalent formulation to discuss the existence of a solution of extended general mixed quasi variational inequality and also develop several iterative methods for solving extended general mixed quasi variational inequality and its variant forms. We consider the convergence analysis of the proposed iterative methods under appropriate conditions. We also introduce a new class of resolvent equation, which is called the extended general implicit resolvent equation and establish an equivalent relation between the extended general implicit resolvent equation and the extended general mixed quasi variational inequality. Some special cases are also discussed.

• Structural Engineering and Mechanics
• /
• 제75권5호
• /
• pp.541-557
• /
• 2020

A novel family of controllable, dissipative structure-dependent integration methods is derived from an eigen-based theory, where the concept of the eigenmode can give a solid theoretical basis for the feasibility of this type of integration methods. In fact, the concepts of eigen-decomposition and modal superposition are involved in solving a multiple degree of freedom system. The total solution of a coupled equation of motion consists of each modal solution of the uncoupled equation of motion. Hence, an eigen-dependent integration method is proposed to solve each modal equation of motion and an approximate solution can be yielded via modal superposition with only the first few modes of interest for inertial problems. All the eigen-dependent integration methods combine to form a structure-dependent integration method. Some key assumptions and new techniques are combined to successfully develop this family of integration methods. In addition, this family of integration methods can be either explicitly or implicitly implemented. Except for stability property, both explicit and implicit implementations have almost the same numerical properties. An explicit implementation is more computationally efficient than for an implicit implementation since it can combine unconditional stability and explicit formulation simultaneously. As a result, an explicit implementation is preferred over an implicit implementation. This family of integration methods can have the same numerical properties as those of the WBZ-α method for linear elastic systems. Besides, its stability and accuracy performance for solving nonlinear systems is also almost the same as those of the WBZ-α method. It is evident from numerical experiments that an explicit implementation of this family of integration methods can save many computational efforts when compared to conventional implicit methods, such as the WBZ-α method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제17권3호
• /
• pp.197-207
• /
• 2013

The Cahn-Hilliard equation was proposed as a phenomenological model for describing the process of phase separation of a binary alloy. The equation has been applied to many physical applications such as amorphological instability caused by elastic non-equilibrium, image inpainting, two- and three-phase fluid flow, phase separation, flow visualization and the formation of the quantum dots. To solve the Cahn-Hillard equation, many numerical methods have been proposed such as the explicit Euler's, the implicit Euler's, the Crank-Nicolson, the semi-implicit Euler's, the linearly stabilized splitting and the non-linearly stabilized splitting schemes. In this paper, we investigate each scheme in finite-difference schemes by comparing their performances, especially stability and efficiency. Except the explicit Euler's method, we use the fast solver which is called a multigrid method. Our numerical investigation shows that the linearly stabilized stabilized splitting scheme is not unconditionally gradient stable in time unlike the known result. And the Crank-Nicolson scheme is accurate but unstable in time, whereas the non-linearly stabilized splitting scheme has advantage over other schemes on the time step restriction.

• 한국전산유체공학회지
• /
• 제16권4호
• /
• pp.72-83
• /
• 2011

A high-order discontinuous Galerkin method for the two-dimensional compressible Navier-Stokes equations was developed on unstructured triangular meshes. For this purpose, the BR2 methd(the second Bassi and Rebay discretization) was adopted for space discretization and an implicit Euler backward method was used for time integration. Numerical tests were conducted to estimate the convergence order of the numerical solutions of the Poiseuille flow for which analytic solutions are available for comparison. Also, the flows around a flat plate, a 2-D circular cylinder, and an NACA0012 airfoil were numerically simulated. The numerical results showed that the present implicit discontinuous Galerkin method is an efficient method to obtain very accurate numerical solutions of the compressible Navier-Stokes equations on unstructured meshes.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제25권4호
• /
• pp.224-261
• /
• 2021

The interest in implicit-explicit (IMEX) integration methods has emerged as an alternative for dealing in a computationally cost-effective way with stiff ordinary differential equations arising from practical modeling problems. In this paper, we introduce implicit-explicit second derivative linear multi-step methods (IMEX SDLMM) with error control. The proposed IMEX SDLMM is based on second derivative backward differentiation formulas (SDBDF) and recursive SDBDF. The IMEX second derivative schemes are constructed with order p ranging from p = 1 to 8. The methods are numerically validated on well-known stiff equations.

• 한국정보과학회논문지:시스템및이론
• /
• 제32권9호
• /
• pp.445-456
• /
• 2005

본 논문은 물리적인 힘을 기반으로 유체의 흐름을 실시간으로 시뮬레이션하기 위하여 유체 의 흐름을 지배하는 Wavier-Stokes 방정식에 대한 빠르고 정확한 풀이 기법을 제안한다 본 논문에서는 Navier-Stokes 방정식에 있는 비선형 항의 속도에 대한 초기값을 Stokes 방정식의 해로써 추정한다. 주어진 비선형 미분방정식의 해에 근사하게 초기값을 추정함으로써 정확하고 안정적인 풀이 기법을 만들 수 있었다. 또한 유한차분법(finite difference method)의 암시적(implicit) 방법 중에서 방대한 계산량을 피할 수 있는 ADI(Alternating Direction Implicit) 방법을 사용함으로써 큰 시간 간격(time-step)에 대해서 시스템이 안정적이며 계산속도 또한 빠르다. 실험 결과들은 특히 연기, 구름과 같이 큰 레이놀드 수(Reynolds number)를 가지는 유체에 대해서 탁월한 성능을 보여주었다.

• 한국전산구조공학회논문집
• /
• 제34권3호
• /
• pp.151-157
• /
• 2021

본 연구에서는 비국부 적분 연산기로 표현되는 페리다이나믹 방정식의 수렴성을 검토한다. 정적/준정적 손상 해석 문제를 효율적으로 해석하기 위해 페리다이나믹 방정식의 implicit 정식화가 필요하다. 이 과정에서 페리다이나믹 비국부 적분 방정식으로부터 대수방정식 형태가 나타나게 되어 시스템 행렬 계산을 위해 많은 시간이 소요되기 때문에, 효율적인 계산을 위해 수렴성이 중요한 요소가 된다. 특히 radial influence 함수를 적분 kernel로 사용하는 경우 fractional Laplacian 적분 방정식이 유도된다. 비국부 적분 연산기의 교윳값 성질에 의해 대수방정식의 condition number가 radial influence 함수의 차수 및 비국부 영역의 크기에 영향을 받는 것이 수학적으로 확인되었다. 본 연구에서는 이를 토대로 균열이 있는 페리다이나믹 정적 해석 문제를 Newton-Raphson 방법으로 해석할 때 적분 커널의 차수, 비국부 영역의 크기 등이 대수방정식의 condition number와 preconditioned conjugate gradient (PCG) 방법으로 계산 시 수렴성 및 계산 시간에 미치는 영향을 수치적으로 분석한다.

• 한국해안해양공학회지
• /
• 제12권2호
• /
• pp.70-80
• /
• 2000

Princeton Ocean Model(POM)에서 지배방정식을 반음해법으로 차분화하여 mode splitting을 제거하고 조간대 처리기법을 도입한 3차원 semi-implicit 모델을 수립하였다 운동방정식의 순압경도력항과 수직방향 와동점 성향, 그리고 수심적분된 연속방정식의 속도발산항을 음해법으로 처리하여 빠르게 움직이는 표면중력파가 야기하는 수치안정조건을 제거하여 더 큰 time step을 사용할 수 있는 semi-implicit 모델을 수립하였다 수립한 3차원 semi-implicit 모델의 유효성과 계산휴율을 확인하기 위하여 이상적인 3차원 수로에 적용한 결과 semi-implicit 모델이 POM과 같은 결과를 주었으나 POM보다 약 4,.4배 정도 빠르게 수행되어 향상된 계산효율을 보여주었다. mode splitting 기법을 사용하는 POM의 유속 결과는 조간대에서 발생한 noise가주 수로에까지 전파되어 불안정 한 결과를 준 반면에 semi-implicit 모델결과는 더 큰 time step을 사용함에도 불구하고 조간대와 주 수로 모두에서 noiserk 없는 안정된 결과를 주었다 3차원 semi-implicit 모델의 현장 적용성을 확인하기 위하여 경기만에 적용한 결과 semi-implicit 모델이 모델 영역 전반에 걸쳐 관측된 조석 및 조류의 크기 및 전파양상을 잘 재현하였다.

• 대한수학회논문집
• /
• 제14권4호
• /
• pp.833-842
• /
• 1999

In this paper, we construct a fully discrete solution of the incompressible Navier Stokes equations using implicit Runge kutta method. We prove the existence of the fully discrete solution.

• East Asian mathematical journal
• /
• 제29권1호
• /
• pp.69-82
• /
• 2013

In this paper we consider a doubly nonlinear Volterra equation related to the Leray-Lions with a nonsmooth kernel. By exploiting a suitable implicit time-discretization technique we obtain the existence of global strong solution.