• 한국지하수토양환경학회지:지하수토양환경
• /
• 제7권4호
• /
• pp.10-16
• /
• 2002

본 연구의 목적은 누수대수층으로 분리된 다층-대수층에 대한 지하수 해석이다. Crank-Nicolson방법에 의한 유한차분법을 적용하여 1차원이며 정상상태인 2중 대수층 구조에 대해 해석해와 비교하였다. 수치해와 해석해는 거의 일치하였으므로 수치해를 2차원의 확장된 다층-대수층 구조에 적용하였다. 이는 한 개 또는 여러 개의 대수층에서 양수하는 경우에 각 대수층에서의 수두값을 계산할 수 있게 하였다. 본 연구는 지하수의 효율적인 운영에 도움이 될 것이다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제15권4호
• /
• pp.291-306
• /
• 2011

We consider the second order Strang splitting method to approximate the solution to the KdV equation. The model equation is split into three sets of initial value problems containing convection and dispersal terms separately. TVD MUSCL or MUSCL scheme is applied to approximate the convection term and the second order centered difference method to approximate the dispersal term. In time stepping, explicit third order Runge-Kutta method is used to the equation containing convection term and implicit Crank-Nicolson method to the equation containing dispersal term to reduce the CFL restriction. Several numerical examples of weakly and strongly dispersive problems, which produce solitons or dispersive shock waves, or may show instabilities of the solution, are presented.

• Journal of applied mathematics & informatics
• /
• 제16권1_2호
• /
• pp.53-68
• /
• 2004

Numerical solutions for the viscous Cahn-Hilliard equation are considered using the Crank-Nicolson type finite difference method which conserves the mass. The corresponding stability and error analysis of the scheme are shown. The decay speeds of the solution in $H^1-norm$ are shown. We also compare the evolution of the viscous Cahn-Hilliard equation with that of the Cahn-Hilliard equation numerically and computationally, which has been given as an open question in Novick-Cohen[13].

• Journal of applied mathematics & informatics
• /
• 제29권5_6호
• /
• pp.1477-1487
• /
• 2011

We propose a numerical method for solving Allen-Cahn equation, in both one-dimensional and two-dimensional cases. The new scheme that is explicit, stable, and easy to compute is obtained and the proposed method provides a straightforward and effective way for nonlinear evolution equations.

• 한국수학교육학회지시리즈B:순수및응용수학
• /
• 제13권1호
• /
• pp.19-38
• /
• 2006

We analyze the spectral collocation approximation for a parabolic partial integrodifferential equations(PIDE) with a weakly singular kernel. The space discretization is based on the spectral collocation method and the time discretization is based on Crank-Nicolson scheme with a graded mesh. We obtain the stability and second order convergence result for fully discrete scheme.

• Advances in Energy Research
• /
• 제5권3호
• /
• pp.255-275
• /
• 2017

The advantages of using natural circulation (NC) as a cooling system, has prompted the worldwide development to investigate this phenomenon more than before. The interesting application of the NC in low power experimental facilities and research reactors, highlights the obligation of study in these laminar flows. The inherent oscillations of NC between hot source and cold sink in low Grashof numbers necessitates stability analysis of cooling flow with experimental or numerical schemes. For this type of analysis, numerical methods could be implemented to desired mass, momentum and energy equations as an efficient instrument for predicting the behavior of the flow field. In this work, using the explicit, implicit and Crank-Nicolson methods, the fluid flow parameters in a natural circulation experimental test loop are obtained and the sensitivity of solving approaches are discussed. In this way, at first, the steady state and transient results from explicit are obtained and compared with experimental data. The implicit and crank-Nicolson scheme is investigated in next steps and in subsequent this research is focused on the numerical aspects of instability prediction for these schemes. In the following, the assessment of the flow behavior with coarse and fine mesh sizes and time-steps has been reported and the numerical schemes convergence are compared. For more detail research, the natural circulation of fluid was modeled by ANSYS-CFX software and results for the experimental loop are shown. Finally, the stability map for rectangular closed loop was obtained with employing the Nyquist criterion.

• 대한토목학회논문집
• /
• 제10권4호
• /
• pp.141-150
• /
• 1990

하구나 연안에서 해수의 순환형상을 모사(模寫)하게 위해 천수방정식(淺水方程式)을 여러 가지 경계조건 하에서 수치해석하였다. 공간영역은 Galerkin방법으로 이산화(離散化)하였으며 시간영역에 대해서는 유한차분법(Crank-Nicolson방법)을 사용하였다. 네 가지 검정실험이 해석적인 해가 있는 일차원 수로에서 행하여졌으며, 해석해를 구할 수 없는 이차원 모형에도 적용되었다. 해석해가 있는 경우 수치모사 결과가 이와 잘 일치하였으며, 이차원 모형에서의 결과도 매우 합당함을 알 수 있었다. 또 일차원 문제에서 4점 bilinear요소와 삼각형 요소를 사용한 결과를 각각 비교하였으며 시간적분도 2단계 Lax-Wendroff방법을 사용하여 결과를 비교하였다. 음해법을 사용할 경우 비교적 정확한 결과를 얻을 수 있으나 요소의 갯수가 많아지면 구성되는 대수방정식(代數方程式)이 커지기 때문에 각 시간마다의 계산량이 엄청나게 늘어나게 되며 양해법을 사용할 때는 원하는 만큼의 정확한 결과를 얻기 위하여 시간간격이나 공간격자 간격을 선정하는데 각별히 유의하여야 할 것이다.

• 대한수학회보
• /
• 제57권1호
• /
• pp.219-244
• /
• 2020

In this paper, we approximate the solution of the coupled nonlinear Schrödinger equations by using a fully discrete finite element scheme based on the standard Galerkin method in space and implicit midpoint discretization in time. The proposed scheme guarantees the conservation of the total mass and the energy. First, a priori error estimates for the fully discrete Galerkin method is derived. Second, the existence of the approximated solution is proved by virtue of the Brouwer fixed point theorem. Moreover, the uniqueness of the solution is shown. Finally, convergence orders of the fully discrete Crank-Nicolson scheme are discussed. The end of the paper is devoted to some numerical experiments.

• Journal of applied mathematics & informatics
• /
• 제39권5_6호
• /
• pp.623-641
• /
• 2021

In this paper, numerical treatment of singularly perturbed parabolic delay differential equations is considered. The considered problem have small delay on the spatial variable of the reaction term. To treat the delay term, Taylor series approximation is applied. The resulting singularly perturbed parabolic PDEs is solved using Crank Nicolson method in temporal direction with non-standard finite difference method in spatial direction. A detail stability and convergence analysis of the scheme is given. We proved the uniform convergence of the scheme with order of convergence O(N^-1 + (∆t)²), where N is the number of mesh points in spatial discretization and ∆t is mesh length in temporal discretization. Two test examples are used to validate the theoretical results of the scheme.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.392-397
• /
• 2001

A numerical study is performed to investigate the interaction between subsonic axial turbine blade boundary layer and periodically oncoming rotor induced wakes. An implicit scheme for solving the compressible Navier-Stokes equation is developed, which adopts a 4th-order compact difference for spatial discretiztion, a 2nd order Crank-Nicolson scheme for temporal discretization and the dynamic eddy viscosity model as the subgrid scale model. The efficiency and the accuracy of the proposed method are verified by applying to some benchmark problems such as laminar cylinder flow, laminar airfoil cascade flow and a transitional flat plate boundary layer flow. Computational results show good agreements with previous experimental and numerical results. Finally, flow through a stator cascade is simulated at $Re = 7.5{\times}10^5$ without free-stream turbulence intensity. The velocity fields and skin friction coefficients in the transitional region show similar trends with previous boundary layer natural transition.