• Journal of Mechanical Science and Technology
• /
• 제16권10호
• /
• pp.1264-1275
• /
• 2002

This paper is concerned with the investigation on the stability and convergence characteristics of the Crank-Nicolson-Galerkin scheme that is widely being employed for the numerical approximation of parabolic-type partial differential equations. Here, we present the theoretical analysis on its consistency and convergence, and we carry out the numerical experiments to examine the effect of the time-step size △t on the h- and P-convergence rates for various mesh sizes h and approximation orders P. We observed that the optimal convergence rates are achieved only when △t, h and P are chosen such that the total error is not affected by the oscillation behavior. In such case, △t is in linear relation with DOF, and furthermore its size depends on the singularity intensity of problems.

• 대한수학회지
• /
• 제61권4호
• /
• pp.621-657
• /
• 2024

In this study, the numerical solution of a space-fractional Burgers' equation with initial and boundary conditions is considered. This equation is the simplest nonlinear model for diffusive waves in fluid dynamics. It occurs in a variety of physical phenomena, including viscous sound waves, waves in fluid-filled viscous elastic pipes, magneto-hydrodynamic waves in a medium with finite electrical conductivity, and one-dimensional turbulence. The proposed QBS/CNG technique consists of the Galerkin method with a function basis of quadratic B-splines for the spatial discretization of the space-fractional Burgers' equation. This is then followed by the Crank-Nicolson approach for time-stepping. A linearized scheme is fully constructed to reduce computational costs. Stability analysis, error estimates, and convergence rates are studied. Finally, some test problems are used to confirm the theoretical results and the proposed method's effectiveness, with the results displayed in tables, 2D, and 3D graphs.

• 대한수학회보
• /
• 제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.

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

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

• 대한기계학회논문집A
• /
• 제24권8호
• /
• pp.2042-2049
• /
• 2000

This paper is concerned with the comparison of forging characteristics between forward and backward processes, through the three-dimensional finite element simulation, for the aluminum powder forging of engine pistons. Starting from the theoretical formulation of velocity and temperature fields in the sintered preform during the process, we examine the comparative distributions of relative density, effective stress and temperature as well as the variations of total forging load and total volume reduction. Through the comparative results, we find out that the forward method provides better forging characteristics than the backward method.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회 논문집
• /
• pp.414-418
• /
• 2005

A finite element model is developed for the numerical simulation of bed elevation change in a curved channel. The SU/PG (Streamline-Upwind/Petrov-Galerkin) method is used to solve 2D shallow water equations and the BG (Bubnov-Galerkin) method is used for the Exner equation. For the time derivative terms, the Crank-Nicolson scheme is used. The developed model is a decoupled model in a sense that the bed elevation does not change simultaneously with the flow during the computational time step. The total load formula with is used for the sediment transport model. The slip conditions are described along the lateral boundaries. The effects of gravity force due to geometry change and the secondary flows in a curved channel are considered in the model. For the verification, the model is applied to two laboratory experiments. The first is $140^{\circ}$ bended channel data at Delft Hydraulics Laboratory and the second is $140^{\circ}$ bended channel data at Laboratory of Fluid Mechanics of the Delft University of Technology. The finite element grid is constructed with linear quadrilateral elements. It is found that the computed results are in good agreement with measured data, showing a point bar at the inner bank and a pool at the outer bank.

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

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.