• Journal of applied mathematics & informatics
• /
• 제40권3_4호
• /
• pp.603-618
• /
• 2022

In this paper, a numerical solution of the generalized diffusion equation with a delay has been obtained by a numerical technique based on the Galerkin finite element method by applying the cubic B-spline basis functions. The time discretization process is carried out using the forward Euler method. The numerical scheme is required to preserve the delay-independent asymptotic stability with an additional restriction on time and spatial step sizes. Both the theoretical and computational rates of convergence of the numerical method have been examined and found to be in agreement. As it can be observed from the numerical results given in tables and graphs, the proposed method approximates the exact solution very well. The accuracy of the numerical scheme is confirmed by computing L₂ and L_∞ error norms.

• 한국대기환경학회지
• /
• 제21권5호
• /
• pp.547-555
• /
• 2005

In this study, an analytical expression for aerosol mass transfer at spherical collector in the low Knudsen number region was obtained. Happel's zero shear stress cell model was extended in the low Knudsen number region and the result was compared with numerical solution results. The zero vorticity model based on the Kuwabara's cell model was also extended in the low Knudsen number region and compared with Happel's results. The results showed that both analytic and numerical solution agree very well with each other in low Knudsen number region. Happel's zero shear stress model also agrees with Kuwabara's zero vorticity model without significant loss of accuracy. The obtained solution converges to the original solution of Lee et al. (1999) when Knudsen number approaches to zero. Subsequently, this study derived most general type of analytic solution for aerosol mass transfer of spherical collector including the finite Knudsen number region.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 추계학술대회논문집A
• /
• pp.458-463
• /
• 2000

The mechanical structures generally have discontinuous parts such as the cracks, notches and holes owing to various reasons. In this paper, in order to analyze effectively these singularity problems using the finite element method, a mixed analysis method which an analytical solution and finite element solutions are simultaneously used is newly proposed. As the analytical solution is used in the singularity region and the finite element solutions are used in the remaining regions except this singular zone, this analysis method reasonably provides for the numerical solution of a singularity problem. Through various numerical examples, it is shown that the proposed analysis method is very convenient and gives comparatively accurate solution.

• 한국지하수토양환경학회지:지하수토양환경
• /
• 제13권2호
• /
• pp.54-61
• /
• 2008

수평 혹은 경사 형태 특수 정호 양수량에 대한 시공간적 수위 강하를 지하수 수치 모델링을 활용하여, 평가하였다. 지하수 수치 모델링은 상용 프로그램인 FEFLOW(version 5.1)의 1차원 선형 불연속 특징 요소를 활용하여 수행되었으며, 수치해의 검증을 위해 Zhan과 Zlotnik(2002)이 제안한 연속된 점 형태 배출원 배열 방식 준 해석해와 비교하였다. 비교 검증 결과, 수치해와 준해석해는 최대 수위 강하가 나타나는 양수 최인접부를 제외하고는 거의 일치한 형태를 보여주었다. 검증된 수치적 방법을 이용하여, 강변여과 방식 취수가 검토되는 현장에 대한 수위강하를 정량적으로 평가할 수 있었다.

• 한국항공우주학회지
• /
• 제30권3호
• /
• pp.27-36
• /
• 2002

AUSM계열 수치기법의 수치적 불안정성에 대한 원인과 해결방안에 대한 연구를 수행하였다. Euler 유동에서 수치적 불안정성은 제어면에 수직한 방향의 유동속도가 영인 영역에서 발생하며 이 영역에서 Eule r 방정식은 근본적으로 부정해를 가지게 되어 무수히 많은 해를 가지게 된다. 지배방정식 자체로는 유일해를 찾는 것이 불가능하고 주위의 유동조건이나 외부교란에 의해 유일해를 결정하게 된다. 이러한 특징은 충격파 영역에서 교란이 존재할 경우 초기 상태에 대한 정보를 상실하게 되어 충격파 불안정성을 유발하게 된다. slip유동을 정확히 계산할 수 있는, 즉 유일해를 결정할 수 없는, 수치기법은 충격파 불안정성을 근본적으로 제거할 수 없다.

• 한국지하수토양환경학회지:지하수토양환경
• /
• 제20권2호
• /
• pp.10-21
• /
• 2015

The present study introduces a novel numerical approach for solving dispersion dominated problems with Cauchy boundary condition in an Eulerian-Lagrangian scheme. The study reveals the incapability of traditional Neuman approach to address the dispersion dominated problems with Cauchy boundary condition, even though it can produce reliable solution in the advection dominated regime. Also, the proposed numerical approach is applied to a real field problem of radioactive contaminant migration from radioactive waste repository which is a major current waste management issue. The performance of the proposed numerical approach is evaluated by comparing the results with numerical solutions of traditional FDM (Finite Difference Method), Neuman approach, and the analytical solution. The results show that the proposed numerical approach yields better and reliable solution for dispersion dominated regime, specifically for Peclet Numbers of less than 0.1. The proposed numerical approach is validated by applying to a real field problem of radioactive contaminant migration from radioactive waste repository of varying Peclet Number from 0.003 to 34.5. The numerical results of Neuman approach overestimates the concentration value with an order of 100 than the proposed approach during the assessment of radioactive contaminant transport from nuclear waste repository. The overestimation of concentration value could be due to the assumption that dispersion is negligible. Also our application problem confirms the existence of real field situation with advection dominated condition and dispersion dominated condition simultaneously as well as the significance or advantage of the proposed approach in the real field problem.

• 한국콘텐츠학회논문지
• /
• 제13권2호
• /
• pp.505-511
• /
• 2013

경사진 지형을 갖는 축대칭 지형에 적용이 가능한 경계처리기법을 개발하였다. 섬 지형의 경우 복잡한 지형으로 인하여 유한요소모형을 사용하여 파의 변형을 해석하는 것이 좋지만 해수와 접하는 섬의 단면이 연직이 아닌 경우에는 수심이 0이 되어 경계면을 적절하게 처리하기 어렵다는 단점이 있다. 본 연구에서는 장파에 대한 해석해를 활용하여 임의의 경사진 경계면에 적용가능한 경계처리기법을 개발하였다. 이를 위해 지배방정식으로 완경사 방정식을 사용하였으며 계산 영역을 해석해 영역과 수치해 영역으로 구분하여 해석해 영역에 기존의 해석해를 적용한 후 수치해와 결합하여 모델을 완성하였다. 유도된 해는 기존의 해석해와 비교하여 그 타당성을 검증하였다.

• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2004년도 학술대회지
• /
• pp.179-182
• /
• 2004

This paper deals with a new mathematical formulation of nonlinear wave profile based on Banach fixed point theorem. As application of the formulation and its solution procedure, some numerical solutions was presented in this paper and nonlinear equation was derived. Also we introduce a new operator for iteration and getting solution. A numerical study was accomplished with Stokes' first-order solution and iteration scheme, and then we can know the nonlinear characteristic of Stokes' high-order solution. That is, using only Stokes' first-oder(linear) velocity potential and an initial guess of wave profile, it is possible to realize the corresponding high-oder Stokian wave profile with tile new numerical scheme which is the method of iteration. We proved the mathematical convergence of tile proposed scheme. The nonlinear strategy of iterations has very fast convergence rate, that is, only about 6-10 iterations arc required to obtain a numerically converged solution.

• Journal of Advanced Marine Engineering and Technology
• /
• 제30권4호
• /
• pp.494-504
• /
• 2006

It is well known that the two-parameter $J-A_2$ solution can well characterize the crack-tip fields and quantify the crack-tip constraint for different flawed geometries in variety of loading conditions. However, this solution fails to do so for bending dominated specimens or geometries at large deformation because of the influence of significant global bending stress on the crack-tip field. To solve this issue, a modified $J-A_2$ solution is developed in this paper by introducing an additional term to address the global bending influence. Using the $J_2$ flow theory of plasticity and within the small-strain framework detailed finite element analyses are carried out for the single edge notched bend (SENB) specimen with a deep crack in A533B steel at different deformation levels ranging from small-scale Yielding to large-scale Yielding conditions. The numerical results of the crack-tip stress field are then compared with those determined from the $J-A_2$ solution and from the modified $J-A_2$ solution at the same level of applied loading Results indicate that the modified $J-A_2$ solution largely improves the $J-A_2$ solution, and match very well with the numerical results in the region of interest at all deformation levels. Therefore, the proposed solution can effectively describe the crack-tip field and the constraint for bending dominated specimens or geometries.

• 한국방재학회 논문집
• /
• 제11권3호
• /
• pp.111-116
• /
• 2011

본 연구에서는 횡방향 하중을 받는 말뚝의 해석 해를 유도하고 그로부터 말뚝의 거동분석과 더불어 긴 말뚝으로의 거동기준이 되는 임계 말뚝길이(critical pile length)를 지반조건을 달리하여 구하고 비교 분석하였다. 또한 해석 해에 의한 말뚝의 거동과 p-y곡선을 적용한 수치해석을 통해 말뚝의 거동을 비교분석하였다. 해석 해에 따르면 밀도를 달리한 지반조건에 대해 무차원 임계 말뚝길이는 해석에서 고려한 세 가지 말뚝머리 경계조건에 있어 2.3~3.2 사이였다. 해석 해에 의한 결과와 수치해석에 의한 결과를 비교하면 말뚝길이에 따른 말뚝의 변형과 모멘트 분포양상은 유사하였다. 말뚝머리 변형량은 해석에 의한 경우가 수치해석에 의한 경우보다 보수적인 값을 보여주었으며 휨모멘트의 값은 해석에 의한 값과 수치해석에 의한 값 사이에 큰 차이를 보이지 않았다.