• East Asian mathematical journal
• /
• 제40권1호
• /
• pp.109-118
• /
• 2024

This paper proposes a numerical method based on domain decomposition to find approximate solutions for one-dimensional convection-diffusion equations with Neumann boundary conditions. First, the equations are transformed into convection-diffusion equations with Dirichlet conditions. Second, the author introduces the Prediction/Correction Domain Decomposition (PCDD) method and estimates errors for the interface prediction scheme, interior scheme, and correction scheme using known error estimations. Finally, the author compares the PCDD algorithm with the fully explicit scheme (FES) and the fully implicit scheme (FIS) using three examples. In comparison to FES and FIS, the proposed PCDD algorithm demonstrates good results.

• 한국수학교육학회지시리즈B:순수및응용수학
• /
• 제17권4호
• /
• pp.289-298
• /
• 2010

In this paper, we propose a second-order prediction/correction (SPC) domain decomposition method for solving one dimensional linear hyperbolic partial differential equation $u_{tt}+a(x,t)u_t+b(x,t)u=c(x,t)u_{xx}+{\int}(x,t)$. The method can be applied to variable coefficients problems and singular problems. Unconditional stability and error analysis of the method have been carried out. Numerical results support stability and efficiency of the method.

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

A continuum-based design sensitivity analysis (DSA) method fur non-shape problems is developed for geometrically nonlinear elastic structures. The non-shape problem is characterized by the design variables that are not associated with the domain of system like sizing, material property, loading, and so on. Total Lagrangian formulation with the Green-Lagrange strain and the second Piola-Kirchhoff stress is employed to describe the geometrically nonlinear structures. The spatial domain is discretized using the 4-node isoparametric plane stress/strain elements. The resulting nonlinear system is solved using the Newton-Raphson iterative method. To take advantage of the derived analytical sensitivity In topology optimization, a fast and efficient design sensitivity analysis method, adjoint variable method, is employed and the material property of each element is selected as non-shape design variable. Combining the design sensitivity analysis method and a gradient-based design optimization algorithm, an automated design optimization method is developed. The comparison of the analytical sensitivity with the finite difference results shows excellent agreement. Also application to the topology design optimization problem suggests a very good insight for the layout design.

• 대한수학회논문집
• /
• 제9권2호
• /
• pp.439-448
• /
• 1994

It is known that the steepest descent(SD) method and the conjugate gradient(CG) method [1, 2, 5, 6] converge when these methods are applied to solve linear systems of the form Ax = b, where A is symmetric and positive definite. For some finite difference discretizations of elliptic problems, one gets positive definite matrices that are almost symmetric. Practically, the SD method and the CG method work for these matrices. However, the convergence of these methods is not guaranteed theoretically. The SD method is also called Orthores(1) in iterative method papers. Elman [4] states that the convergence proof for Orthores($\kappa$), with $\kappa$ a positive integer, is not heard. In this paper, we prove that the SD method and the CG method converge when the $\iota$$^2$ matrix norm of the non-symmetric part of a positive definite matrix is less than some value related to the smallest and the largest eigenvalues of the symmetric part of the given matrix.(omitted)

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

Many partial differential equations defined on a rectangular domain can be solved numerically by using a domain decomposition method. The most commonly used decompositions are the domain being decomposed in stripwise and rectangular way. Theories for non-overlapping domain decomposition(in which two adjacent subdomains share an interface) were often focused on the stripwise decomposition and claimed that extensions could be made to the rectangular decomposition without further discussions. In this paper we focus on the comparisons of the two ways of decompositions. We consider the unconditionally stable scheme, the MIP algorithm, for solving parabolic partial differential equations. The SOR iterative method is used in the MIP algorithm. Even though the theories are the same but the performances are different. We found out that the stripwise decomposition has better performance.

• Coupled systems mechanics
• /
• 제11권2호
• /
• pp.107-119
• /
• 2022

The aim of the work presented in this paper is development of numerical model for prediction of temperature distribution in pavement according to the measured meteorological parameters, with introduction of non-linear heat transfer coefficient which is a function of temerature difference between the air and the pavement. Developed model calculates heat radiated from the pavement back in the air, which is an important part of the heat trasfer process in the open air surfaces. Temperature of the pavement surface, heat radiation together with many meteorological parameters were measured in series during two years in order to validate the model and calibrate model parameters. Special finite element method for temperature heat transfer towards the soil together with the time integration scheme are used to solve the governing equation. It is proved that non-linear heat transfer coefficient, which is a function of time and temperature difference between the air and the pavement, is required to decribe this phenomena. Proposed model includes heat tranfer coefficient callibration for specific climate region, through the iterative inverse procedure.

• 한국해안해양공학회지
• /
• 제10권2호
• /
• pp.93-99
• /
• 1998

타원형 유한차분모형에서 개방 경계조건으로 포물선 근사식과 스폰지층 경계를 사용하여 모형의 개량을 도모하였다. 수치기법은 GCG(Generalized conjugate gradient)기법을 사용하였고 구형해저실험에서 포물형 근사식을 사용하여 부적절한 반사파를 상당 부분 제거할 수 있었다. 스폰지층 경계의 경우 2파장 이상의 스폰지층을 사용할 때 포물형 근사식과 유사한 결과를 얻을 수있었다. 직사각형 항만에 대한 실험을 통하여 임의 형상의 대상 해역에도 쉽게 모형을 적용할 수 있음을 확인하였다.

• Geomechanics and Engineering
• /
• 제4권4호
• /
• pp.263-279
• /
• 2012

The present paper deals with the analysis of combined footings resting on geosynthetic reinforced granular fill overlying stone column improved poor soil. An attempt has been made to study the influence of inclusion of geosynthetic layer on the deflection of the footing. The footing has been idealized as a beam having finite flexural rigidity. Granular fill layer has been represented by Pasternak shear layer and stone columns and poor soil have been represented by nonlinear Winkler springs. Nonlinear behavior of granular fill layer, stone columns and the poor soil has been considered by means of hyperbolic stress strain relationships. Governing differential equations for the soil-foundation system have been derived and solution has been obtained employing finite difference scheme by means of iterative Gauss Elimination method. Results of a detailed parametric study have been presented, for a footing supporting typically five columns, in non-dimensional form in respect of deflection with and without geosynthetic inclusion. Geosynthetic layer has been found to significantly reduce the deflection of the footing which has been quantified by means of parametric study.

• 대한기계학회논문집
• /
• 제11권2호
• /
• pp.262-270
• /
• 1987

본 논문에서는 Eulerian 좌표축에서 온도 경계층과 혼합된 비정상 경계층 유 동을 효과적으로 다룰 수 있는 좌표변환 방법과 수치해석 방법을 도입하였다. 이전 의 수치적 방법들은 축차적이며, 뒷정체점 부근의 경계층 두께가 시간과 더불어 지수 적으로 증가하므로, 격자점의 수를 경계층 두께의 증가에 따라 상당히 늘려 주어야 된 다. 그러나 여기에서는, 격자점의 수를 무리하게 늘릴 필요가 없으며 일반적인 비선 형 차분방정식을 정확도를 떨어뜨리지 않고 선형화시킴으로써 비축차적인 수치해들이 허용된다. 이런 선형화 방법은 Beam & Warming에 의해 최초로 압축성 Navier-Stokes 방정식에 사용되었고, Orlandi & Ferziger, Kim & Chang에 의해 경계층 유동에 확장되 었다.또한 뒷정체점 부근에서 경계층 두께의 증가로 인한, 격자점 증가의 필요를 피하기 위하여, 몇 가지 서로 다른 종류의 변환변수들을 시간과 공간에 따라 선별적으 로 사용하여 수치적인 경계층 두께가 거의 일정하도록 만들었다. 이와 같은 변환변 수들은 또한 쌍-포물선형인 현 지배방정식의 초기조건들을 쉽게 구할 수 있도록 허용 해 준다.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2008년도 공동학술대회
• /
• pp.51-56
• /
• 2008

등방성 매질에서의 파형역산에 대한 연구는 1980년대부터 꾸준히 이루어져 왔으나 이방성 매질에 대한 연구는 그렇지 못하다. 본 연구에서는 이방성 매질에 대한 시간영역 셀기반 유한 차분 모델링 기법을 이용해 2차원 TI 구조에서의 파형역산 알고리듬을 개발하였다. 반복적인 비선형 역산에서 최대 급경사 방향은 역시간 구조보정의 역전파 방법을 이용하여 간접적으로 계산하였고, 이를 정규화 시키기 위해 슈도-헤시안 행렬을 이용하였다. 본 연구에서 제시된 시간영역 파형역산 기법을 이방성 매질을 포함한 2층 구조와 이방성 Marmousi 모형 자료에 적용하고 이를 등방성 매질만을 고려한 기존의 파형역산 결과와 비교하였다. 본 연구의 결과를 통해 이방성 매질을 등방성 매질로 가정하고 파형역산을 수행할 경우 정확한 영상을 얻을 수 없기 때문에, 실제 탐사 자료의 파형역산을 수행할 경우 이방성 매질을 고려해야 좀 더 정확한 지하 구조를 파악할 수 있음을 확인하였다.