• Nuclear Engineering and Technology
• /
• 제53권7호
• /
• pp.2084-2094
• /
• 2021

Delta-form-based methods for solving high order spatial discretization schemes are introduced into the reactor S_N transport equation. Due to the nature of the delta-form, the final numerical accuracy only depends on the residuals on the right side of the discrete equations and have nothing to do with the parts on the left side. Therefore, various high order spatial discretization methods can be easily adopted for only the transport term on the right side of the discrete equations. Then the simplest step or other robust schemes can be adopted to discretize the increment on the left hand side to ensure the good iterative convergence. The delta-form framework makes the sweeping and iterative strategies of various high order spatial discretization methods be completely the same with those of the traditional S_N codes, only by adding the residuals into the source terms. In this paper, the flux limiter method and weighted essentially non-oscillatory scheme are used for the verification purpose to only show the advantages of the introduction of delta-form-based solving methods and other high order spatial discretization methods can be also easily extended to solve the S_N transport equations. Numerical solutions indicate the correctness and effectiveness of delta-form-based solving method.

• Structural Engineering and Mechanics
• /
• 제3권4호
• /
• pp.391-410
• /
• 1995

The application of adaptive finite element method to dynamic problems is investigated. Both the kinetic and strain energy errors induced by space and time discretization were estimated in a consistent manner and controlled by the simultaneous use of the adaptive mesh generation and the automatic time stepping. Also an optimal ratio of spatial discretization error to temporal discretization error was discussed. In this study it was found that the best performance can be obtained when the specified spatial and temporal discretization errors have the same value. Numerical examples are carried out to verify the performance of the procedure.

• 한국수학교육학회지시리즈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.

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

In order to reduce the excessive numerical dissipation, the new spatial discretization scheme is introduced. The present method in this paper has the formula that has an additional procedure of defining transferred properties at a cell-interface, based on AUSMPW+. The newly defined transferred property could eliminate numerical dissipation effectively in non-flow aligned grid system. In addition, the present method guarantees the monotonic characteristic in capturing a discontinuity. Through a stationary or moving contact discontinuity and a stationary or moving shock discontinuity, a vortex discontinuity and shock wave/ boundary layer interaction, it is verified that the accuracy of the present method is improved.

• 대한기계학회논문집B
• /
• 제29권6호
• /
• pp.703-710
• /
• 2005

A finite element discretization of the advection and redistancing equations of level set method has been studied. It has been shown that Galerkin spatial discretization combined with Crank-Nicolson temporal discretization of the advection equation of level set yields a good result and that consistent streamline upwind Petrov-Galerkin(CSUPG) discretization of the redistancing equation gives satisfactory solutions for two test problems while the solutions of streamline upwind Petrov-Galerkin(SUPG) discretization are dissipated by the numerical diffusion added for the stability of a hyperbolic system. Furthermore, it has been found that the solutions obtained by CSUPG method are comparable to those by second order ENO method.

• Geomechanics and Engineering
• /
• 제17권3호
• /
• pp.227-236
• /
• 2019

In practice, the natural slopes are frequently with soils of spatial properties and irregular features. The traditional limit analysis method meets an inherent difficulty to deal with the stability problem for such slopes due to the normal condition in the associated flow rule. To overcome the problem, a novel technique based on the upper bound limit analysis, which is called the discretization technique, is employed for the stability evaluation of complex slopes. In this paper, the discretization mechanism for complex slopes was presented, and the safety factors of several examples were calculated. The good agreement between the discretization-based and previous results shows the accuracy of the proposed mechanism, proving that it can be an alternative and reliable approach for complex slope stability analysis.

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

• 한국항공우주학회지
• /
• 제33권9호
• /
• pp.56-65
• /
• 2005

시간유한요소법은 시간영역을 고정시키고 행렬 미분방정식 형태의 공간전파 관계식을 풂으로써 시간과 공간에 대한 동적 해석을 수행하는 방법이다. 이 방법은 공간이산화 유한요소법이나 시/공간 동시이산화 유한요소법에 비해 공간에 관한 자유도가 발생하는 것이 두드러진 특징으로, 이를 이용하여 분포형 구동기의 공간에 따른 특성을 최적화하는 데에 효율적으로 사용될 수 있다. 본 논문에서는 임의의 초기조건을 반영할 수 있도록 구성된 상태변수 벡터를 이용하여 구조물을 시간영역에서 이산화하고, 공간영역에서 전파관계식 및 경계조건을 이용하여 공간전파 관계식을 형성하였다. 이 때 구동기의 공간에 따른 형상 분포는 설계되어야 할 변수의 함수이고, 시간반응은 형상함수를 이용하여 이산화 하였다. 포텐셜 에너지 및 운동에너지를 구조물의 변위제어에 적절한 최적의 성능지수로 설정하고, 이를 최소화하도록 미지의 함수인 구동기의 분포형상을 구하였다. 일반적으로 구조물은 임의의 초기조건에서 외란을 받게 되나, 본 연구에서는 구현가능한 제어법칙을 이용하여 최종시간에서 안정화(rest) 조건을 만족한다고 가정하였다. 구동기 분포형상 최적화를 위해 상태/준상태 방정식을 유도하였다. 서브행렬 재형상화와 시/공간 경계조건을 통해 상태변수와 준상태변수에 대한 Ricatti 미분방정식을 유도하였다. 이를 통해 구동기 분포형상 최적화를 구현하였으며, 수치 시뮬레이션을 통해 적절한 구동기의 분포형상 최적화를 수행할 수 있음을 보였다.

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

An automated procedure which allows adaptation of spatial and time discretization simultaneously in finite element analysis of linear elastodynamic problems is presented. For dynamic problems having responses dominated by high frequency modes, such as those with impact, explosive, traveling and earthquake loads high gradient stress regions change their locations from time to time. And the time step size may need to vary in order to deal with whole process ranging from transient phase to steady state phase. As the sizes of elements in space vary in different regions, the procedure also permits different time stepping. In such a way, the best performance attainable by the finite element method can be achieved. In this study, we estimate both of the kinetic energy error and stran energy error induced by spatial and time discretization in a consistent manner. Numerical examples are used to demonstrate the performance of the procedure.

• Journal of applied mathematics & informatics
• /
• 제4권1호
• /
• pp.17-30
• /
• 1997

In this work we approximate the solution of initialboun-dary value problem using a Galerkin-finite element method for the spatial discretization and Implicit Runge-Kutta method for the spatial discretization and implicit Runge-Kutta methods for the time stepping. To deal with the nonlinear term f(x, t, u), we introduce the well-known extrapolation sheme which was used widely to prove the convergence in $L_2$-norm. We present computational results showing that the optimal order of convergence arising under $L_2$-norm will be preserved in $L_{\infty}$-norm.