• 제목/요약/키워드: Crank-Nicholson

검색결과 21건 처리시간 0.032초

종확산 방정식에 대한 Eulerian-Lagrangian 연산자 분리방법 (Eulerian-Lagrangian Split-Operator Method for the Longitudinal Dispersion Equation)

  • 전경수;이길성
    • 대한토목학회논문집
    • /
    • 제14권1호
    • /
    • pp.131-141
    • /
    • 1994
  • 특성곡선을 고려한 세가지 연산자 분리방법을 오염원의 종확산 문제에 적용하여, 그 결과를 Eulerian 기법들의 계산결과와 비교하였다. 연산자 분리방법의 이송방정식에 대한 수치 기법들로는 generalized upwind, two-point fourth-order 및 sixth-order Holly-Preissmann 기법들을 각각 적용하였으며, 확산 방정식에 대한 수치기법으로는 Crank-Nicholson 기법을 적용하였다. Holly-Preissmann 기법을 사용하는 연산자 분리방법들이 Eulerian 기법들에 비하여 매우 정확한 계산결과를 나타내었다. Eulerian 기법들의 경우에는 이송항의 근사방법으로서 중앙차분을 취하는 기법들은 수치진동을, 후방차분을 취하는 기법들은 수치분산을 각각 보였으며, 이러한 현상들은 종확산계수의 값이 작을수록 더욱 뚜렷하게 나타났다.

  • PDF

Eulerian-Lagrangian 방법을 이용한 1차원 종확산방정식의 수치모형 (Numerical Modeling of One-Dimensional Longitudinal Dispersion Equation using Eulerian-Lagrangian Method)

  • 서일원;김대근
    • 물과 미래
    • /
    • 제27권2호
    • /
    • pp.155-166
    • /
    • 1994
  • Eulerian-Lagrangian 방법을 이용하여 1차원 종확산방정식의 수치모형을 비교·분석하였다. 본 연구에서서 비교·분석한 모형은 지배방정식을 연산자 분리방법에 의해서 이송만을 지배하는 이송방정식과 확산만을 지배하는 확산방정식으로 분리한다. 이송방정식은 특성곡선을 따라서 유체입자를 추적하는 특성곡선법을 사용하여 해를 구하고, 그 결과를 고정된 Eulerian 격자상에 보간하였고, 확산방정식은 상기 고정격자상에서 Crank-Nicholson 유한차분법을 사용하여 해를 구하였다. 이송방정식의 풀이에서 다양한 보간방법이 적용되었는데, 일반적으로 Hermite 보간다항식을 사용한 경우가 Lagrange 보간다항식을 사용한 경우보다 수치확산 및 수치진동 등의 오차를 최소화할 수 있어서 더욱 우수한 것으로 밝혀졌다.

  • PDF

난류전단 흐름에서의 비정상 수평 선오염원의 확산에 관한 수치해법 (A Numerical Method for Dispersion of Unsteady Horizontal Line Source in Turbulent Shear Flow)

  • 전경수
    • 물과 미래
    • /
    • 제29권4호
    • /
    • pp.187-198
    • /
    • 1996
  • 난류전단 흐름에서의 비정상 수평 선오염원의 확산에 관한 수치모형을 개발하였다. 계산기법으로는 비정상 이송확산 방정식을 종방향 이송 및 연직방향 확산으로 분리하고, 이들 방정식을 방시간 간격에 대하여 번갈아 계산하는 단계분리 유한차분기법을 사용하였다. 종방향 이송방정식에 대해서는 Holly-Preissmann 기법을, 연직방향 확산방정식에 대해서는 Crank-Nicholson 기법을 각각 적용하였다. 개발된 모형을 난류전단 흐름에서의 정상 수평 선오염원의 확산 문제에 적용하여 계산결과를 반해석해와 비교함으로써 모형을 검증하였다. 또한 난류전단 흐름내로 순간적으로 방류된 면오염원의 확산문제에 계산모형을 적용하였다. 마찰계수에 대한 민감도 분석 결과, 동일한 무차원 시각에서의 혼합 정도는 마찰계수에 관계없이 거의 일정한 것으로 나타났다. 또한 동일한 정도의 혼합상태에 도달하는 데 소요되는 유하거리는 마찰계수의 제곱근에 반비례함을 알 수 있었다.

  • PDF

Eulerian-Lagrangian 혼합모형에 의한 종확산 방정식의 수치해법 (An Eulerian-Lagrangian Hybrid Numerical Method for the Longitudinal Dispersion Equation)

  • 전경수;이길성
    • 물과 미래
    • /
    • 제26권3호
    • /
    • pp.137-148
    • /
    • 1993
  • 종확산 방정식에 대한 유한차분 모형으로서, 5차의 보간다항식을 사용한 Holly-Preissmann 기법과 Generalized Crank-Nicholson 기법을 결합한 혼합모형을 개발하였다. 순간적으로 부하된 오염원의 종확산문제에 본 모형 및 특성곡선을 고려한 다른 수치기법들을 적용하여 정확해와 비교하였다. 보 모형에 의한 계산결과, Courant 수에 관계없이 수치진동이 전혀 발생하지 않았으며, 최대농도 발생지점도 정확해와 일치하였다. 모형의 적용에 있어서 시간가중치 $\theta$의 값이 작을수록 계산의 정확성이 전반적으로 향상되는 것으로 나타났으며, $\theta$의 값을 크게 할수록 최대농도값을 과대평가하는 경향을 보였다. 전반적으로 Courant 수가 작을수록 정확한 계산결과를 나타내고 있으나 그 민감도는, 특히 $\theta$의 값이 작을수록, 매우 작게 나타났다. 3차의 보간다항식을 사용하는 혼합모형 및 연산자 분리방법들과의 비교결과, 이송항이 지배적일수록 본 모형이 정확해와 가장 근사한 계산결과를 보임을 알 수 있었다.

  • PDF

유한차분 모형에 의한 일차원 이송-확산방정식 계산결과의 비교 (Comparison of the Results of Finite Difference Method in One-Dimensional Advection-Dispersion Equation)

  • 이희영;이재철
    • 물과 미래
    • /
    • 제28권4호
    • /
    • pp.125-136
    • /
    • 1995
  • 특정곡선을 고려한 ELM을 이송-확산방정식에 적용하여 그 결과를 Eulerian 기법(Stone-Brian, QUICKEST)과 비교하였다. 이송항의 계산을 위해서는 Lagrangian 보간법과 Cubic spline 보간법을 이용하였고 확산항의 계산에 있어서는 Crank-Nicholson 방법을 이용하였다. 수치모형의 적용결과는 다음과 같다. (1) Gaussian hill에의 적용:Lagrangian 보간법을 사용하여 계산한 경우가 가장 정확한 결과를 보였다. Cubic spline 보간법을 사용한 경우와 QUICKEST 방법의 경우에는 Peclet수가 50인 경우에 감쇠현상을 보였다. Stone-Brian방법은 Peclet수 10,50에서 위상오차가 발생하였다. (2) Advanced front에의 적용: 모든 방법이 Peclet수 1,4에서 정확한 결과를 얻었다. Peclet수가 50인 경우에 Lagrangian 보간법을 사용하여 계산한 경우와 Stone-Brian 방법은 증폭오차가 발생하였고 Cubic spline 보간법을 사용한 경우와 QUICKEST 방법의 경우는 수치진동 현상을 보였다.

  • PDF

비보존성 오염물질에 대한 종확산 방정식의 수치해법 (A Numerical Method for Longitudinal Dispersion Equation for Nonconservative Contaminants)

  • 유명관;전경수
    • 한국수자원학회논문집
    • /
    • 제32권6호
    • /
    • pp.607-616
    • /
    • 1999
  • 비보존성 오염물질의 종확산에 관한 수치모형을 개발하였다. 계산기법으로는 종확산 방정식을 이송, 감쇠 및 확산 방정식으로 분리하고, 이들 방정식을 1/3 시간 간격에 대하여 번갈아 계산하는 단계분리 유한차분기법을 사용하였다. 이송방정식에 대해서는 Holly-Preissmann 기법을, 감쇠방정식에 대해서는 해석적 방법을, 확산방정식에 대해서는 Crank-Nicholson 기법을 각각 사용하였다. 오염물질이 불균일 흐름 내로 연속적으로 유입되는 경우 및 균일 흐름 내로 순간적으로 부하되는 경우에 대한 종확산 문제에 모형을 적용하여 계산결과를 정확해와 비교함으로써 모형을 검증하였다. 또한 감쇠방정식의 수치해법으로써 Euler 방법을 사용하는 기존의 모형에 계산결과를 비교하였다. 감쇠계수가 커질수록 본 모형이 기존의 모형에 비하여 더욱 정확한 계산결과를 나타내었다.

  • PDF

상용 여객기의 유동 및 공력 특성 해석 (Flow and Aerodynamic Characteristics Analyses of A Commercial Passenger Airplane)

  • 김양균;김성초;김정수;이기만;진학수
    • 대한기계학회:학술대회논문집
    • /
    • 대한기계학회 2007년도 춘계학술대회B
    • /
    • pp.2857-2861
    • /
    • 2007
  • Flow and aerodynamic characteristics were analyzed numerically for a commercial passenger airplane, Boeing 747-400, flying in the cruising condition. The model geometry with 100:1 in scale was obtained by the photo scanning measurement with the maximum error of 1.4% comparing with the real airplane dimension. The three-dimensional inviscid steady compressible governing equations were solved by the finite volume method in the unstructured grid system. The convective terms were treated by the Crank-Nicholson and first-order upwind schemes. In the computational results, the strong wing-tip vortices were clearly observed and the pressure contours on the airplane surface were suggested. The lift and drag forces in the wing with engines increase by 1.49% and 3.9%, respectively compared with the case without engines. The aerodynamic forces were estimated quantitatively for each element which consists of the airplane.

  • PDF

Numerical Modeling of One-Dimensional Longitudinal Dispersion Equation using Eulerian Method

  • Seo, Il-Won;Kim, Dae-Geun
    • Korean Journal of Hydrosciences
    • /
    • 제6권
    • /
    • pp.51-66
    • /
    • 1995
  • Various Eulerian-Lagerangian numerical models for the one-dimensional longtudinal dispersion equation are studied comparatively. In the models studied, the transport equation is decoupled into two component parts by the operator-splitting approach ; one part governing advection and the other dispersion. The advection equation has been solved using the method of characteristics following flud particles along the characteristic line and the result are interpolated onto an Eulerian grid on which the dispersion equation is solved by Crank-Nicholson type finite difference method. In solving the advection equation, various interpolation schemes are tested. Among those, Hermite interpo;ation po;ynomials are superor to Lagrange interpolation polynomials in reducing both dissipation and dispersion errors.

  • PDF

A study on temporal accuracy of OpenFOAM

  • Lee, Sang Bong
    • International Journal of Naval Architecture and Ocean Engineering
    • /
    • 제9권4호
    • /
    • pp.429-438
    • /
    • 2017
  • Cranke-Nicolson scheme in native OpenFOAM source libraries was not able to provide 2nd order temporal accuracy of velocity and pressure since the volume flux of convective nonlinear terms was 1st accurate in time. In the present study the simplest way of getting the volume flux with 2nd order accuracy was proposed by using old fluxes. A possible numerical instability originated from an explicit estimation of volume fluxes could be handled by introducing a weighting factor which was determined by observing the ratio of the finally corrected volume flux to the intermediate volume flux at the previous step. The new calculation of volume fluxes was able to provide temporally accurate velocity and pressure with 2nd order. The improvement of temporal accuracy was validated by performing numerical simulations of 2D Taylor-Green vortex of which an exact solution was known and 2D vortex shedding from a circular cylinder.

비압축성 Navier-Stokes 방정식에 대한 내재적 속도 분리 방법 (An implicit velocity decoupling procedure for the incompressible Navier-Stokes equations)

  • 김경연;백승진;성형진
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 2000년도 추계 학술대회논문집
    • /
    • pp.129-134
    • /
    • 2000
  • An efficient numerical method to solve the unsteady incompressible Navier-Stokes equations is developed. A fully implicit time advancement is employed to avoid the CFL(Courant-Friedrichs-Lewy) restriction, where the Crank-Nicholson discretization is used for both the diffusion and convection terms. Based on a block LU decomposition, velocity-pressure decoupling is achieved in conjunction with the approximate factorization. Main emphasis is placed on the additional decoupling of the intermediate velocity components with only n th time step velocity The temporal second-order accuracy is Preserved with the approximate factorization without any modification of boundary conditions. Since the decoupled momentum equations are solved without iteration, the computational time is reduced significantly. The present decoupling method is validated by solving the turbulent minimal channel flow unit.

  • PDF