• 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.

• 한국전산유체공학회지
• /
• 제12권1호
• /
• pp.43-52
• /
• 2007

A comparative study on flux functions for the 2-dimensional Euler equations has been conducted. Explicit 4-stage Runge-Kutta method is used to integrate the equations. Flux functions used in the study are Steger-Warming's, van Leer's, Godunov's, Osher's(physical order and natural order), Roe's, HLLE, AUSM, AUSM+, AUSMPW+ and M-AUSMPW+. The performance of MUSCL limiters and MLP limiters in conjunction with flux functions are compared extensively for steady and unsteady problems.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 추계 학술대회논문집
• /
• pp.36-40
• /
• 2006

A comparative study on flux functions for the 2-dimensional Euler equations has been conducted. Explicit 4-stage Runge-Kutta method is used to integrate the equations. Flux functions used in the study are Steger-Warming's, van Leer's. Godunov's, Osher's(physical order and natural order), Roe's, HILE, AUSM, AUSM+ and AUSMPW+. The performance of MUSCL limiters and MLP limiters in conjunction with flux functions are compared extensively for steady and unsteady problems.

• 대한수학회보
• /
• 제60권4호
• /
• pp.985-1001
• /
• 2023

Recently, Alzer and Choi [2] introduced and studied a set of the four linear Euler sums with parameters. These sums are parametric extensions of Flajolet and Salvy's four kinds of linear Euler sums [9]. In this paper, by using the method of residue computations, we will establish two explicit combined formulas involving two parametric linear Euler sums S⁺⁺_p,q (a, b) and S^+-_p,q (a, b) defined by Alzer and Choi, which can be expressed in terms of a linear combinations of products of trigonometric functions, digamma functions and Hurwitz zeta functions.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제26권1호
• /
• pp.49-66
• /
• 2022

The G-Euler process has been proposed to overcome the difficulties of the calculation of the exponential function of the Jacobian. It is an explicit method that uses the exponential function of the scalar skew-symmetric matrix. We define the moving shapes of true solutions and the moving shapes of numerical solutions. It is discussed whether the moving shape of the numerical solution matches the moving shape of the true solution. The match rates of these two kinds of moving shapes are sequentially calculated by the G-Euler process without using the true solution. It is shown that the closer the minimum match rate is to 100%, the more closely the numerical solutions follow the true solutions to the end. The minimum match rate indicates the reliability of the numerical solution calculated by the G-Euler process. The graphs of the Lorenz system in Perko [1] are different from those drawn by the G-Euler process. By the way, there is no basis for claiming that the Perko's graphs are reliable.

• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.66-72
• /
• 1997

Optimum nozzle design exploiting the method of characteristic(M.O.C) has been in application as an efficient design methodology targeting a less weighted and short expansion nozzle. This paper treats the optimum nozzle design and the analysis of the inviscid compressible flow inside. Based on traditional Rao's method, the optimum nozzle design is coded with minor modifications for the identification of the control surface across which the mass flux should be conserved. Internal flow field is simulated numerically by M.O.C and implicit/explicit Taylor-Galerkin finite element method(F.E.M) with the aid of adaptive remeshing to capture the shock wave, hence improve the accuracy. Designed and calculated flow fields due to the separate analyses show that the mass flux predicted by optimum nozzle design with M.O.C is not conserved across the control surface and the sonic line should be located upstream of the nozzle throat. Rao's optimum nozzle design methodology exaggerates the momentum thrust and tends to overemphasize the engine performance loss.

• 융합보안논문지
• /
• 제5권3호
• /
• pp.93-97
• /
• 2005

MacCormack 방법은 hyperbolic 편미분 방정식의 근을 구하는데 많이 쓰이는 방법으로 그 정확도가 2차 오더가 된다. 하지만 이 방법으로 편미분방정식을 풀 경우 불연속인 점에서는 엔트로피를 만족하지 않는 경우가 있어 우리는 임의의 항을 첨가하여 근을 구해야한다. 이 임의의 항을 첨가하지 않고 직접 방정식으로부터 구하는 방법을 생각하는데 있어서 기존의 MacCormack 방법에 새 central scheme의 개념을 이용하면 전형적인 MacCormack 방법의 정확도와 장점을 보존할 수 있다. 이 새로운 방법을 이용하여 1D Burgers' 방정식과 1D Euler gas dynamic 방정식에 활용하여 그 결과를 살펴본다.

• 대한수학회보
• /
• 제61권3호
• /
• pp.671-697
• /
• 2024

In this paper, we formally introduce the notion of a general parametric digamma function Ψ(−s; A, a) and we find the Laurent expansion of Ψ(−s; A, a) at the integers and poles. Considering the contour integrations involving Ψ(−s; A, a), we present some new identities for infinite series involving Dirichlet type parametric harmonic numbers by using the method of residue computation. Then applying these formulas obtained, we establish some explicit relations of parametric linear Euler sums and some special functions (e.g. trigonometric functions, digamma functions, Hurwitz zeta functions etc.). Moreover, some illustrative special cases as well as immediate consequences of the main results are also considered.

• 한국콘텐츠학회논문지
• /
• 제11권9호
• /
• pp.1-8
• /
• 2011

본 연구는 질량-스프링 모델 기반의 직물 모델에서 질점의 움직임을 분석하여 실시간 직물 애니메이션이 가능한 퍼지 추론 기법을 제안한다. 지금까지 직물과 같은 탄성체를 표현하기 위한 많은 기법들은 질량-스프링 모델을 사용하였다. 직물은 다수의 질량과 스프링의 조합으로 구성되어 변형 가능한 면을 이루게되고, 면의 움직임은 운동법칙을 기반으로 수치적분을 통해 계산될 수 있다. 제안된 방법과 동일한 직물구조에서 Explicit 오일러 방법은 ${\Delta}t$ > 0.01 일 경우 불안정성 문제가 나타났으며, Implicit 오일러 방법은 ${\Delta}t$ = 0.03 에서도 애니메이션이 생성되지만 많은 양의 선형 시스템을 계산해야 하는 단점을 가지고 있어서 실시간 처리에 부적합하다. 본 연구는 질량-스프링 모델에서 질점의 움직임을 계산하기 위하여 ${\Delta}t$ = 0.03을 가지면서도 실시간 처리가 가능한 방법을 제안한다.

• 에너지공학
• /
• 제6권1호
• /
• pp.52-57
• /
• 1997

유동방향으로 초원형 형상을 갖는 잘록한 관내의 압축성 유동에 관한 동특성을 수치적으로 연구하였다. MacCormack의 양해법, 즉 예측자/보정자 단계를 거친 시간 진행법을 이용하여 Euler 방정식의 해를 구하였는데, 관내 유동은 이차원, 비점성, 압축성 유동이라 가정하였다. 관의 직경비와 형상비가 압력분포에 미치는 영향을 광범위하게 고찰하였으며, 본 연구에서 개발한 전산프로그램을 이용한 수치 결과는 상용코드인 FLUENT를 이용한 결과와 비교하여 일치된 결과를 얻을 수 있었다.