본 논문은 2차원 천수방정식을 해석하기 위해 새롭고 간단한 MUSCL 재구성법을 제안하였다. 수정 MUSCL 기법은 보존변수의 재구성을 위해 계산격자와 인접격자의 보존변수 차에 대해 각 경계면에 균일하게 분배하는 기존 방법 대신 면적가중비를 사용하였으며, 이 방법은 정형격자 뿐만 아니라 비정형 격자의 사용에도 보존변수의 물리적 재구성이 가능하다. 또한, 본 연구에서는 비구조적 격자의 적용이 가능한 차원비분리 기법을 적용하였으며, 수정 MUSCL 기법의 사용으로 발생할 수 있는 수치진동을 제어하기 위해 TVD 기법의 경사제한자를 사용하였다. 하상경사항의 정확하고 효율적인 수치 처리를 위해 수정 MUSCL 기법을 수면경사법과 연계하였다. 제안한 기법을 적용한 유한체적모형을 건물의 영향을 고려한 댐 붕괴 해석 및 Bellos의 댐 붕괴 실험에 적용하고, 적용결과를 실험실 자료 및 기존 연구자의 계산결과와 비교하여 개발모형을 검증하였다.
충격파의 높이나 속도는 홍수제어조작이나 수로벽과 빠른 유속을 가지는 하천에서 교량의 설계에 중요한 자료가 된다. 따라서 광범위한 조건에서 흐름의 불연속면을 모의할 수 있는 수치모형이 요구된다. 본 연구에서는 천수방정식을 지배방정식으로 한 Godunov 형 유한체적법 모형을 개발하였다. Riemann 해법으로 Roe(1981)의 해법이 사용된다. 이 모형은 본 연구에서 비구조적격자(unstructured grids)를 사용하기 위해 개발된 수정 MUSCL을 도입하였다. 양해법을 쓰는 본 모형은 시간간격을 자동 계산한다. 개발된 모형을 전형적인 이차원 댐 붕괴파 모의, 수리모형 실험에서 행해진 붕괴파 모의, 그리고 수리모형 실험에서 행해진 만곡수로에서의 정상상태 모의 등에 적용하였다. 그 적용결과에 의해 다음과 같이 결론을 내었다. 1)유한체적법은, 충격파 모의를 위한 수치해석 기법인 Godunov 형 방법과 잘 결합될 수 있기 때문에 충격파를 모의하기에 적당한 방법이다. 2)수정 MUSCL과 결합된 유한체적법 모형이 충격파를 잘 포착함으로써 수정 MUSCL의 적용성이 입증되었다.
T. Buffard and S. Clain은 자료의 선형 재구축에서 계산 격자의 각 변에서 보존변수의 경사가 설정되는 다중 경사 MUSCL을 제안하였다. 이 연구에서는 천수방정식에 대한 수치모형의 개발을 위해 비구조 격자에서 공간에 대한 2차 정확도를 얻을 수 있는 다중 경사 MUSCL을 적용하였으며, 흐름률의 계산을 위해 근사 Riemann 해법 중에서 HLLL 기법을 이용하였다. 모형의 적용성 검토를 위해 유럽의 IMPACT 사업의 일환으로 실시된 '고립된 건물 시험'과 '모형 도시 홍수 실험'과 비교하였다. 건물의 전면에서 저항에 의한 갑작스런 수심 상승과 건물 사이 특정 위치에서 수심의 예측에는 한계가 있었으나, 이것은 같은 실험에 대한 여러 모형들의 적용 결과에서도 나타나는 문제인 것으로 확인되었다. 보다 세분된 계산 격자에서 '모형 도시 홍수 실험'에 대한 모의 결과가 측정값에 잘 부합되는 것으로 나타났다. 개발된 모형으로 댐 붕괴나 돌발 홍수에 의한 도시 침수와 같은 복잡한 현상을 잘 모의할 수 있음이 확인되었다.
로터 공기역학 모사를 위한 수치해석 기법이 제안되었다. 국부적인 유동특성변수를 4차 정확도 내삽의 MUSCL 접근방법을 유한체적 공식에 적용시킴으로써 로터 블레이드로 부터 멀리 떨어져 성긴 격자영역에서의 비 물리적인 수치 확산을 개선시켰다. 또한, 유동특성변수에 따라 각각 다른 제한자를 적용함으로써 수치적인 소산작용을 억제하고 수치적 안정성을 높였다. 대표적인 로터 공기역학 응용문제에 적용한 결과 후류포획성능이 크게 개선됨을 확인하였으나 공력소음계산 결과는 전통적인 MUSCL 기법에 비교하여 큰 개선점은 없었다.
천수방정식과 같은 쌍곡선형 미분방정식의 불연속 해에 대한 Riemann 해법은, 1950년대 말 공기동역학 분야에서 S. K. Godunov의 선구적인 시도 이후, 다양한 영역에서 성공적으로 적용되고 있다. 당초 제안된 해법은 공간에 대해 1차 정도였으나, 2차의 정도를 얻을 수 있는 기법이 1970년대 말 B. van Leer에 의해 제안되었으며, MUSCL로 불린다. 서로 인접한 격자의 보존변수가 고려된 경사가 도입되어 두 격자에 의해 공유되는 변의 좌 우에서 선형으로 보존변수가 재구축되는 MUSCL은 제한자와 함께 이용될 때, 구조 격자 체계에서 비교적 단순하면서도 효과적인 적용성이 입증되었다. 그런데, 이 기법을 2차원의 비구조 격자 체계에 적용하는 경우, 인접한 모든 격자의 보존변수를 고려한 평면의 경사를 결정해야 하는 어려움이 따른다. 특히, 삼각형 비구조 격자에 적용할 경우 최적의 평면을 결정하기 위해 Green-Gauss 적분식이나 최소-자승법 등을 이용하게 된다. 이에 비해, 2010년 T. Buffard와 S. Clain이 제안한 다중경사 기법은 격자의 각 변에서 경사가 각각 결정되는 방법으로 계산량이 많은 Green-Gauss 적분식이나 최소자승법을 피할 수 있는 장점이 있는 것으로 알려져 있다. 정확해가 알려진 두 경우에 대해 몇 가지 제한자를 적용한 결과를 1차 정도의 해와 함께 비교하였으며, superbee 제한자에 의한 결과가 우수하였으나, 희유파와 충격파가 맞닿는 곳에서 수치 분산이 나타났다. minmod 제한자의 결과가 대체로 무난하였으며, 이를 2차원 댐 붕괴 문제에 적용하여 1차 정도의 해와 비교하였다. 마찰이 없고 초기 수심이 댐 상류에서 10 m, 하류에서 5 m로서 물이 차 있는 경우, 1차 정도의 해에서 나타나는 수치 소산이 2차 정도에서는 발생되지 않았다. 댐 하류에서 초기에 수심이 영으로 바닥이 드러난 경우에서 마찰의 영향을 검토하였다. 마찰이 있는 경우, 마찰 경사 항의 Manning 계수를 0.04로 두었으며, 마찰에 의한 영향이 잘 드러났다. 수심이 50 mm 보다 작은 경우에는 마찰을 적용하지 않았다. 이 연구는 환경부 '차세대 핵심환경기술개발 사업'의 지원에 의한 것이다.
최근 개발된 FUNWAVE-TVD 파랑모형을 이용하여 적용되어 온 TVD 기법들의 수렴도와 수치적인 안정성을 비교하였다. Yamamoto and Daiguji(1993)의 minmod limiter를 사용하는 4차 정확도의 MUSCL-TVD 기법과 Erduran et al.(2005)의 van-Leer limiter를 사용하는 4차 정확도의 MUSCL-TVD 기법, Zhou et al.(2001)의 van-Leer limiter를 사용하는 2차 정확도의 MUSCL-TVD 기법을 비교하였으며, 수리실험 관측치가 제시되어 있는 Vincent and Briggs(1989)의 불규칙 파랑실험에 적용하였다. 불규칙 파랑의 비쇄파 실험 결과에서 minmod limiter를 사용하는 4차 정확도의 기법은 van-Leer limiter를 사용하는 기법이 요구하는 격자의 크기만큼 세밀한 격자를 요구하지는 않지만, 더 낮은 CFL을 사용해야 안정적인 모의가 가능하였다. 반면에 van-Leer limiter를 사용하는 기법에서는 numerical dissipation을 줄이기 위하여 보다 세밀한 격자를 필요로 하지만 비교적 높은 CFL을 사용할 수 있는 것으로 나타났다. 각 기법의 numerical dissipation의 영향을 최대한 줄이기 위하여 공간격자를 충분히 줄인 쇄파 모의 실험에서는 비쇄파 실험에 비하여 각 기법의 특성이 명확히 나타났다. Numerical dissipation이 상대적으로 작은 minmod limiter를 사용하는 기법으로 모의할 때는 격자를 충분히 줄이면 수치적인 불안정성이 나타나며 수치해가 발산하는 결과를 보였지만, van-Leer limiter를 사용하는 기법에서는 비교적 낮은 CFL을 사용하여 쇄파 모의가 완료되었으며, 관측치를 잘 재현하는 결과를 보였다.
본 연구에서는 2차원 부정류 해석을 위한 유한체적 수치모형을 개발하였다. 개발된 모형에서는 천수방정식의 수치해를 구하기 위하여 Lax-Friedrichs 해법에 근간을 둔 ENO (essentially non-oscillatory) 기법을 적용하였으며, 공간적으로 2차 정도 모의를 위하여 수정 MUSCL 기법을 적용하였다. 개발된 모형의 적용성을 평가하기 위하여 해석해가 있는 1차원 댐 붕괴파 모의에 적용한 결과, 해석해와 유사한 결과를 나타냈으며, 수정 MUSCL 기법의 모의결과를 minmod 경사제한자를 사용한 모의결과와 비교한 결과, 수정 MUSCL 기법을 사용하는 경우 해석해와 더 근접한 결과를 나타냈다. 댐붕괴파의 2차원 해석에 대한 모형의 적용성을 평가하기 위하여 대칭형, 비대칭형 댐 붕괴에 적용한 결과, 저수지에서 수로 구간으로 충격파 전파를 잘 모의했으며, 기존의 연구결과와 유사한 모의결과를 나타냈다. 개수로의 부정류 해석에 대한 모형의 적용성을 평가하기 위하여 말단이 폐쇄된 수로의 입구에서 정상파가 주어졌을 때 수위와 유속을 모의하여 해석해와 비교하였으며, 입구와 말단이 폐쇄된 수로에서 파의 진동을 모의하여 모형의 거동을 평가하였다.
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.
Journal of the Korean Society for Industrial and Applied Mathematics
/
제15권4호
/
pp.291-306
/
2011
We consider the second order Strang splitting method to approximate the solution to the KdV equation. The model equation is split into three sets of initial value problems containing convection and dispersal terms separately. TVD MUSCL or MUSCL scheme is applied to approximate the convection term and the second order centered difference method to approximate the dispersal term. In time stepping, explicit third order Runge-Kutta method is used to the equation containing convection term and implicit Crank-Nicolson method to the equation containing dispersal term to reduce the CFL restriction. Several numerical examples of weakly and strongly dispersive problems, which produce solitons or dispersive shock waves, or may show instabilities of the solution, are presented.
Performance of first, second and third order accurate methods for calculation of in viscid fluxes in fluid flow governing equations are investigated here. For the purpose, an upwind method based on Roe's scheme is used to solve 2-dimensional Euler equations. To increase the accuracy of the method two different schemes are applied. The first one is a second and third order upwind-based algorithm with the MUSCL extrapolation Van Leer (1979), based on primitive variables. The other one is an upwind-based algorithm with the Chakravarthy extrapolation to the fluxes of mass, momentum and energy. The results show that the thickness of shock layer in the third order accuracy is less than its value in second order. Moreover, applying limiter eliminates the oscillations near the shock while increases the thickness of shock layer especially in MUSCL method using Van Albada limiter.
본 웹사이트에 게시된 이메일 주소가 전자우편 수집 프로그램이나
그 밖의 기술적 장치를 이용하여 무단으로 수집되는 것을 거부하며,
이를 위반시 정보통신망법에 의해 형사 처벌됨을 유념하시기 바랍니다.
[게시일 2004년 10월 1일]
이용약관
제 1 장 총칙
제 1 조 (목적)
이 이용약관은 KoreaScience 홈페이지(이하 “당 사이트”)에서 제공하는 인터넷 서비스(이하 '서비스')의 가입조건 및 이용에 관한 제반 사항과 기타 필요한 사항을 구체적으로 규정함을 목적으로 합니다.
제 2 조 (용어의 정의)
① "이용자"라 함은 당 사이트에 접속하여 이 약관에 따라 당 사이트가 제공하는 서비스를 받는 회원 및 비회원을
말합니다.
② "회원"이라 함은 서비스를 이용하기 위하여 당 사이트에 개인정보를 제공하여 아이디(ID)와 비밀번호를 부여
받은 자를 말합니다.
③ "회원 아이디(ID)"라 함은 회원의 식별 및 서비스 이용을 위하여 자신이 선정한 문자 및 숫자의 조합을
말합니다.
④ "비밀번호(패스워드)"라 함은 회원이 자신의 비밀보호를 위하여 선정한 문자 및 숫자의 조합을 말합니다.
제 3 조 (이용약관의 효력 및 변경)
① 이 약관은 당 사이트에 게시하거나 기타의 방법으로 회원에게 공지함으로써 효력이 발생합니다.
② 당 사이트는 이 약관을 개정할 경우에 적용일자 및 개정사유를 명시하여 현행 약관과 함께 당 사이트의
초기화면에 그 적용일자 7일 이전부터 적용일자 전일까지 공지합니다. 다만, 회원에게 불리하게 약관내용을
변경하는 경우에는 최소한 30일 이상의 사전 유예기간을 두고 공지합니다. 이 경우 당 사이트는 개정 전
내용과 개정 후 내용을 명확하게 비교하여 이용자가 알기 쉽도록 표시합니다.
제 4 조(약관 외 준칙)
① 이 약관은 당 사이트가 제공하는 서비스에 관한 이용안내와 함께 적용됩니다.
② 이 약관에 명시되지 아니한 사항은 관계법령의 규정이 적용됩니다.
제 2 장 이용계약의 체결
제 5 조 (이용계약의 성립 등)
① 이용계약은 이용고객이 당 사이트가 정한 약관에 「동의합니다」를 선택하고, 당 사이트가 정한
온라인신청양식을 작성하여 서비스 이용을 신청한 후, 당 사이트가 이를 승낙함으로써 성립합니다.
② 제1항의 승낙은 당 사이트가 제공하는 과학기술정보검색, 맞춤정보, 서지정보 등 다른 서비스의 이용승낙을
포함합니다.
제 6 조 (회원가입)
서비스를 이용하고자 하는 고객은 당 사이트에서 정한 회원가입양식에 개인정보를 기재하여 가입을 하여야 합니다.
제 7 조 (개인정보의 보호 및 사용)
당 사이트는 관계법령이 정하는 바에 따라 회원 등록정보를 포함한 회원의 개인정보를 보호하기 위해 노력합니다. 회원 개인정보의 보호 및 사용에 대해서는 관련법령 및 당 사이트의 개인정보 보호정책이 적용됩니다.
제 8 조 (이용 신청의 승낙과 제한)
① 당 사이트는 제6조의 규정에 의한 이용신청고객에 대하여 서비스 이용을 승낙합니다.
② 당 사이트는 아래사항에 해당하는 경우에 대해서 승낙하지 아니 합니다.
- 이용계약 신청서의 내용을 허위로 기재한 경우
- 기타 규정한 제반사항을 위반하며 신청하는 경우
제 9 조 (회원 ID 부여 및 변경 등)
① 당 사이트는 이용고객에 대하여 약관에 정하는 바에 따라 자신이 선정한 회원 ID를 부여합니다.
② 회원 ID는 원칙적으로 변경이 불가하며 부득이한 사유로 인하여 변경 하고자 하는 경우에는 해당 ID를
해지하고 재가입해야 합니다.
③ 기타 회원 개인정보 관리 및 변경 등에 관한 사항은 서비스별 안내에 정하는 바에 의합니다.
제 3 장 계약 당사자의 의무
제 10 조 (KISTI의 의무)
① 당 사이트는 이용고객이 희망한 서비스 제공 개시일에 특별한 사정이 없는 한 서비스를 이용할 수 있도록
하여야 합니다.
② 당 사이트는 개인정보 보호를 위해 보안시스템을 구축하며 개인정보 보호정책을 공시하고 준수합니다.
③ 당 사이트는 회원으로부터 제기되는 의견이나 불만이 정당하다고 객관적으로 인정될 경우에는 적절한 절차를
거쳐 즉시 처리하여야 합니다. 다만, 즉시 처리가 곤란한 경우는 회원에게 그 사유와 처리일정을 통보하여야
합니다.
제 11 조 (회원의 의무)
① 이용자는 회원가입 신청 또는 회원정보 변경 시 실명으로 모든 사항을 사실에 근거하여 작성하여야 하며,
허위 또는 타인의 정보를 등록할 경우 일체의 권리를 주장할 수 없습니다.
② 당 사이트가 관계법령 및 개인정보 보호정책에 의거하여 그 책임을 지는 경우를 제외하고 회원에게 부여된
ID의 비밀번호 관리소홀, 부정사용에 의하여 발생하는 모든 결과에 대한 책임은 회원에게 있습니다.
③ 회원은 당 사이트 및 제 3자의 지적 재산권을 침해해서는 안 됩니다.
제 4 장 서비스의 이용
제 12 조 (서비스 이용 시간)
① 서비스 이용은 당 사이트의 업무상 또는 기술상 특별한 지장이 없는 한 연중무휴, 1일 24시간 운영을
원칙으로 합니다. 단, 당 사이트는 시스템 정기점검, 증설 및 교체를 위해 당 사이트가 정한 날이나 시간에
서비스를 일시 중단할 수 있으며, 예정되어 있는 작업으로 인한 서비스 일시중단은 당 사이트 홈페이지를
통해 사전에 공지합니다.
② 당 사이트는 서비스를 특정범위로 분할하여 각 범위별로 이용가능시간을 별도로 지정할 수 있습니다. 다만
이 경우 그 내용을 공지합니다.
제 13 조 (홈페이지 저작권)
① NDSL에서 제공하는 모든 저작물의 저작권은 원저작자에게 있으며, KISTI는 복제/배포/전송권을 확보하고
있습니다.
② NDSL에서 제공하는 콘텐츠를 상업적 및 기타 영리목적으로 복제/배포/전송할 경우 사전에 KISTI의 허락을
받아야 합니다.
③ NDSL에서 제공하는 콘텐츠를 보도, 비평, 교육, 연구 등을 위하여 정당한 범위 안에서 공정한 관행에
합치되게 인용할 수 있습니다.
④ NDSL에서 제공하는 콘텐츠를 무단 복제, 전송, 배포 기타 저작권법에 위반되는 방법으로 이용할 경우
저작권법 제136조에 따라 5년 이하의 징역 또는 5천만 원 이하의 벌금에 처해질 수 있습니다.
제 14 조 (유료서비스)
① 당 사이트 및 협력기관이 정한 유료서비스(원문복사 등)는 별도로 정해진 바에 따르며, 변경사항은 시행 전에
당 사이트 홈페이지를 통하여 회원에게 공지합니다.
② 유료서비스를 이용하려는 회원은 정해진 요금체계에 따라 요금을 납부해야 합니다.
제 5 장 계약 해지 및 이용 제한
제 15 조 (계약 해지)
회원이 이용계약을 해지하고자 하는 때에는 [가입해지] 메뉴를 이용해 직접 해지해야 합니다.
제 16 조 (서비스 이용제한)
① 당 사이트는 회원이 서비스 이용내용에 있어서 본 약관 제 11조 내용을 위반하거나, 다음 각 호에 해당하는
경우 서비스 이용을 제한할 수 있습니다.
- 2년 이상 서비스를 이용한 적이 없는 경우
- 기타 정상적인 서비스 운영에 방해가 될 경우
② 상기 이용제한 규정에 따라 서비스를 이용하는 회원에게 서비스 이용에 대하여 별도 공지 없이 서비스 이용의
일시정지, 이용계약 해지 할 수 있습니다.
제 17 조 (전자우편주소 수집 금지)
회원은 전자우편주소 추출기 등을 이용하여 전자우편주소를 수집 또는 제3자에게 제공할 수 없습니다.
제 6 장 손해배상 및 기타사항
제 18 조 (손해배상)
당 사이트는 무료로 제공되는 서비스와 관련하여 회원에게 어떠한 손해가 발생하더라도 당 사이트가 고의 또는 과실로 인한 손해발생을 제외하고는 이에 대하여 책임을 부담하지 아니합니다.
제 19 조 (관할 법원)
서비스 이용으로 발생한 분쟁에 대해 소송이 제기되는 경우 민사 소송법상의 관할 법원에 제기합니다.
[부 칙]
1. (시행일) 이 약관은 2016년 9월 5일부터 적용되며, 종전 약관은 본 약관으로 대체되며, 개정된 약관의 적용일 이전 가입자도 개정된 약관의 적용을 받습니다.