• 제목/요약/키워드: Runge-Kutta 해법

검색결과 23건 처리시간 0.019초

낮은 Mach수유동 해석을 위한 Preconditioning 가중계수의 의존성 (DEPENDENCE OF WEIGHTING PARAMETER IN PRECONDITIONING METHOD FOR SOLVING LOW MACH NUMBER FLOW)

  • 안영준;신병록
    • 한국전산유체공학회지
    • /
    • 제15권2호
    • /
    • pp.55-61
    • /
    • 2010
  • A dependence of weighting parameter in preconditioning method for solving low Mach number flow with incompressible flow nature is investigated. The present preconditioning method employs a finite-difference method applied Roe‘s flux difference splitting approximation with the MUSCL-TVD scheme and 4th-order Runge-Kutta method in curvilinear coordinates. From the computational results of benchmark flows through a 2-D backward-facing step duct it is confirmed that there exists a suitable value of the weighting parameter for accurate and stable computation. A useful method to determine the weighting parameter is introduced. With this method, high accuracy and stable computational results were obtained for the flow with low Mach number in the range of Mach number less than 0.3.

캐비테이션 유동해석을 위한 기-액 2상 국소균질 모델 (GAS-LIQUID TWO-PHASE HOMOGENEOUS MODEL FOR CAVITATING FLOW)

  • 신병록
    • 한국전산유체공학회지
    • /
    • 제12권2호
    • /
    • pp.53-62
    • /
    • 2007
  • A high resolution numerical method aimed at solving cavitating flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media at isothermal condition and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

압축성 기-액 이상매체중의 고속 유동현상 (HIGH-SPEED FLOW PHENOMENA IN COMPRESSIBLE GAS-LIQUID TWO-PHASE MEDIA)

  • 신병록
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 2007년도 추계 학술대회논문집
    • /
    • pp.249-257
    • /
    • 2007
  • A high resolution numerical method aimed at solving gas-liquid two-phase flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

  • PDF

부분 내재적 조화 균형법을 이용한 주기적인 2차원 비정상 유동 해석 (2-D Periodic Unsteady Flow Analysis Using a Partially Implicit Harmonic Balance Method)

  • 임동균;박수형;권장혁
    • 한국항공우주학회지
    • /
    • 제38권12호
    • /
    • pp.1153-1161
    • /
    • 2010
  • 본 연구에서는 주기적 비정상 유동 해석을 위해 푸리에 변환을 이용하는 조화 균형법의 효율적인 해법을 제안한다. 내재적으로 유속항을 처리하고 외재적으로 조화 원천항을 처리하였다. 외재적 조화 균형법 보다 더 빠르게 수렴 시킬 수 있으며 내재적 조화 균형법을 적용할 때 추가되는 자코비안 행렬을 처리할 필요가 없다. 또한 완전 내재적 기법에 상응하는 수준의 수렴안정성을 확인할 수 있었다. 2차원 비정상 유동 문제로 피칭하는 NACA0012 익형에 적용하였으며 이중 시간 적분법 및 외재적 Runge-Kutta기법의 해와 매우 일치하는 결과를 얻었다.

3차원 데토네이션 파의 수평 및 대각선 모드 파면 구조 (Transverse and Diagonal Mode Structures of Three-dimensional Detonation Wave)

  • 조덕래;최정열
    • 한국추진공학회:학술대회논문집
    • /
    • 한국추진공학회 2005년도 제24회 춘계학술대회논문집
    • /
    • pp.343-346
    • /
    • 2005
  • 전산 해석과 병렬처리를 이용하여 정사각 관 내부를 전파하는 데토네이션파의 삼차원 구조를 살펴보았다. 가변 비열비 공식과 간략화 된 일단계 Arrhenius 반응 모델과 연계된 비점성 유체 방정식을 MUSCL 기반 TVD 해법과 4단계 Runge-Kutta 적분 방법을 이용하여 해석하였다. 삼차원에서의 비정상 해석 결과로부터 그을음 막 기록(smoked-foil record)에서 같은 길이와 다른 폭을 가지는 수평 및 대각 방향 불안정에 의한 상세한 파면 구조를 파악할 수 있었다.

  • PDF

원호형 곡선보의 면외 자유진동에 관한 수치해석적 연구 (Out of Plane Free Vibrations of Circular Curved Beams)

  • 이병구;오상진
    • 전산구조공학
    • /
    • 제9권1호
    • /
    • pp.133-139
    • /
    • 1996
  • 이 논문은 원호형 곡선보의 면외 자유진동에 관한 연구이다. 곡선보 요소의 동적 평형방정식에 Timoshenko 이론을 적용하여 원호형 곡선보의 자유진동을 지배하는 상미분방정식을 유도하고 이를 수치해석하여 고유진동수를 산출할 수 있는 개략해법 중 하나인 수치해석기법을 개발하였다. 수치해석기법에서 미분방정식의 수치적분은 Runge-Kutta method를 이용하였고, 고유진동수의 결정은 Regular-Falsi method를 이용하였다. 실제 수치해석예에서는 회전-회전보, 고정-고정보에 대하여 시행하고 고유진동수에 미치는 무차원 변수들의 영향을 고찰하였다.

  • PDF

불연속 갤러킨 유한요소법을 이용한 1차원 천수방정식의 댐 붕괴류 및 천이류 해석 (Dam-Break and Transcritical Flow Simulation of 1D Shallow Water Equations with Discontinuous Galerkin Finite Element Method)

  • 윤광희;이해균;이남주
    • 대한토목학회논문집
    • /
    • 제34권5호
    • /
    • pp.1383-1393
    • /
    • 2014
  • 최근, 급속한 컴퓨터 하드웨어의 성능 향상과 전산유체역학 분야의 이론적 발전으로, 고차 정확도의 수치기법들이 계산수리학 분야에 적용되어 왔다. 본 연구에서는 1차원 천수방정식에 대한 수치 해법으로 TVD Runge-Kutta 불연속 갤러킨(RKDG) 유한요소법을 적용하였다. 대표적인 천이류(transcritical flow)의 예로 순간적인 댐 붕괴에 의한 댐 붕괴류(dam-break flow) 흐름과 지형변화에 의한 천이류를 모의하였다. 리만(Riemann) 근사해법으로 로컬 Lax-Friedrichs (LLF), Roe, HLL 흐름률(flux) 기법을 사용하였고, 불필요한 진동을 제거하기 위하여, 기울기 제한자로서 MUSCL 제한자를 사용하였다. 개발된 모델은 1차원 댐 붕괴류와 천이류에 적용하였다. 수치해석 결과는 해석해, 수리실험 결과와 비교하였다.

캐비테이션 유동해석을 위한 기- 2상 국소균질 모델 -제2보: 기-액 2상 매체중의 고속유동현상 (GAS-LIQUID TWO-PHASE HOMOGENEOUS MODEL FOR CAVITATING FLOW -Part II. HIGH SPEED FLOW PHENOMENA IN GAS-LIQUID TWO-PHASE MEDIA)

  • 신병록;박선호;이신형
    • 한국전산유체공학회지
    • /
    • 제19권3호
    • /
    • pp.91-97
    • /
    • 2014
  • A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

CERN이 제안한 미임계 평형 원자로에 대한 소고

  • 노태완;이지영;이재중
    • 한국원자력학회:학술대회논문집
    • /
    • 한국원자력학회 1997년도 춘계학술발표회논문집(1)
    • /
    • pp.96-101
    • /
    • 1997
  • 양자가속기를 외부 중성자 공급원으로 하여 미임계 운전을 가능하게 하고 토륨을 핵연료로 사용하므로 장주기 핵종과 핵무기 재료물질의 발생량을 현저히 줄일 수 있는 새로운 노형인 energy amplifier에 대한 연구가 CERN을 중심으로 활발히 진행되고 있다. 본 연구에서는 토륨주기에 대하여 고정 중성자속 조사에 의한 핵분열 및 방사붕괴에 관한 모델을 정립하여 다수의 연립선형 미분방정식으로 구성하여 Runge Kutta 5-6차 자동시간 간격 수치해법을 이용하여 계산하였다. 결과는 1014의 고정 중성자속에 대하여 충분한 U233의 생산이 평형상태에 도달하고 장주기 핵종도 우라늄 주기에 비하여 현저히 줄어듬을 보이므로 가속기를 이용한 토륨 핵연료 주기의 타당성을 확인하였다.

  • PDF

직교이방성 원통쉘의 동적좌굴 (Dynamic bucking of orthotropic cylindrical shells)

  • ;;김천욱
    • 대한기계학회논문집
    • /
    • 제5권4호
    • /
    • pp.266-273
    • /
    • 1981
  • 갑자기 작용하는 외압을 받는 직교이방성원통쉘의 동적좌굴을 해석하였다. Donnell-Karnam 형의 비선형방정식을 유도하였으며 쉘의 초기불완전성도 고려하였다. Galerkin의 방법을 사용하여 운동방정식을 구하고 Runge-Kutta 수직해법으로 비선형방정식을 풀었다. 쉘의 직교이방성특성이 처짐-하중 관계식의 비선형성에 미치는 영향을 검토하였으며 동적산출하중의 판별법을 정의하였다. 본 연구의 결과, 직교이방성원통쉘은 쉘의 초기불완전성에 그리 민감하지 않음을 보여주었다.