• Journal of Ocean Engineering and Technology
• /
• v.12 no.1
• /
• pp.85-98
• /
• 1998

본 논문의 수치해법은 경계치문제를 풀기 위하여 코시이론(Cauchy's theorem)을 사용하였다. 경계치문제는 완전한 물체표면조건과 자유표면조건을 만족시키는 초기치문제로 귀결된다. 현 수치해법에서 무한영역은 수치계산 영역인 비선형 영역과 선형 자유표면조건을 만족하는 선형영역으로 나누어진다. 선형영역의 해는 과도 그린(Green)함수를 사용하여 정합조건을 부과함으로써, 수치계산은 비선형 영역에서만 수행된다. 본 논문에서 저자는 수치계산 영역에서 코시이론을 사용하여 적분방정식을 도출하였고, 무한영역의 해는 정합면에서 과도 그린함수를 사용하여 표현하였다. 본 수치계산에서 자유표면에 요소 재분배법을 적용함으로써 쇄파현상에 대해서도 안정적인 수치해석을 할 수 있었다. 본 논문에서 개발된 수치방법을 적용한 문제는 다음과 같다. 첫째는 자유표면에서 실린더가 강제동요하는 경우에 자유표면형상과 힘을 계산하여 이전의 실험치 및 계산치와 비교하였다. 두번째로는 실린더가 자유수면하에서 일정한 속도로 항주하는 경우에는 조파저항과 양력을 계산하여 고차 스펙트럴법과 비교하였다.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.6
• /
• pp.1250-1256
• /
• 1989

본 연구에서는 수치계산의 불안정성 때문에 과거 고차 정확도 유한차분법 들로 수치계산하기 어려웠던 날카로운 돌출부 위를 흐르는 난류유동현상을 ASQUICK 차분법으로 수치해석하고 그 수렴해를 얻음으로써 수치계산의 안정성을 확인하며 수치계산 결과를 HYBRID차분법에 의한 수치해석치 및 재순환난류유동의 특성치라 할 수 있는 재순환영역길이 실험치와 상호비교함으로써 수치계산의 정확도 향상효과를 확인하고자 한다. 더불어 수치계산의 수렴과정에 큰 영향을 미치는 것으로 잘 알려 져 있는 속도 압력교정법으로써의 PISO법을 ASQUICK차분법과 함께 사용해 PISO-ASQU- ICK조합의 유용성을 확인하고자 한다.

• Journal of Advanced Marine Engineering and Technology
• /
• v.21 no.4
• /
• pp.331-355
• /
• 1997

본 연구에서는 에어챔버가 설치된 펌프관로계에서 수격현상에 대한 현장실험을 실시하고, 에어챔버의 입력변수 및 설계인자가 수격현상에 어떠한 영향을 미치는지 수치계산과 실험을 통하여 조사한 후 다음과 같은 결론을 얻었다. (1) 수치계산에 사용된 입력변수중 폴리트로프지수, 오리피스 유량계수, 압력파 전파속도순으로 수격현상에 대한 영향이 크게 나타났으며, 수치계산은 이 변수값이 각각 1.3, 0.7, 1,050m/s일 때 측정된 압력변동 및 공기체적변화를 가장 잘 모사하였다. (2) 에어챔버내 초기공기체적이 증가하면 압축공기의 압력변화율이 작아지기 때문에 수격현상의 주기는 길어지고, 관로내 압력변화는 감소하였다. 측관에 설치한 오리피스 내경을 줄임으로써 상승압은 상당히 완화되었으며, 시간이 지남에 따라 수격작용도 보다 빨리 감쇄되었다. 수치계산결과는 수격현상의 1-2주기까지 측정값과 비교적 잘 일치하였으며, 펌프관로계에 작용하는 최대.최소압력 및 그 발생시간도 합리적으로 예측되었다. (3) 수치계산은 측정된 과도특성 값들을 약 5%이내의 오차로 예측하였으며, 개발된 전산프로그램은 실제 펌프관로계에 대한 수격완화장치를 설계하는데 매우 유용하게 사용될 것으로 판단된다.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.54-62
• /
• 1998

본 논문은 2차원 해면효과의 수치계산 결과를 정리하였다. 지면으로부터의 높이변화에 따른 점성유동장을 계산하기 위하여 지배방정식으로는 비압축성 RANS 방정식을, 시간에 대하여서는 음해법으로 프로그램을 구성하였다. 압력항은 가상압축성과 4차 수치확산항을 추가하는 것에 의해 계산하였으며, 높은 레이놀즈 수에서의 효과적인 계산을 위해 Baldwin- Lomax 난류모델을 도입하였다. 해면효과가 없는 무한유중에서의 NACA-0012 단면 계산결과를 실험 데이터와 비교하는 것에 의해 프로그램의 타당성을 확인하였다. NACA-6409와 두께 비 4.6%의 날개에 대하여 해면효과를 고려한 계산을 수행하였다. 계산결과, 높이의 변화에 따라 계산된 무차원계수, 압력 및 속도분포는 해면효과의 특성을 잘 보여주고 있다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.88-88
• /
• 2019

컴퓨터 성능향상과 수치해석기법의 발달로 인해 Navier-Stokes 방정식에 기초한 수치모델을 활용한 3차원 유동/파동장 해석이 증가하고 있는 추세이다. 그러나 아직까지 Navier-Stokes 방정식 모델의 계산부하를 PC에서 소화하기에는 무리가 따른다. 게다가 실험실 스케일을 벗어나, 실제 현장을 계산영역으로 설정할 경우에는 계산량이 엄청나게 증가하게 된다. 이것을 극복하기 위해서는 반듯이 병렬계산을 수행하여야 한다. 본 연구에서는 계산부하가 큰 Navier-Stokes 방정식 기반의 3차원 수치모델 LES-WASS-3D를 활용한 대용량 병렬계산체계를 구축한다. 나아가 3차원 정밀 또는 광역의 유동/파동장 해석에 있어서 병렬계산체계의 성능과 적용성을 검토한다. 현재 보급되고 있는 PC들은 모두 멀티프로세서가 장착됨으로 손쉽게 병렬계산을 수행할 수 있다. 그러나 정밀 또는 광역해석을 위해서는 대용량 병렬계산 컴퓨터가 요구된다. 따라서 본 연구에서는 보조프로세서를 장착한 공유메모리 환경의 고성능 병렬계산체계를 구축한다. 나아가 포트란 기반의 순차코드로 구축된 기존 3차원 Navier-Stokes 방정식 모델 LES-WASS- 3D를 병렬코드로 변환한다. 병렬계산 성능 및 적용성을 검토하기 위한 수치해석을 수행한다. 이상의 과정을 통해 본 연구에서 구축한 병렬계산체계의 성능 및 적용성을 확인할 수 있었다. 그리고 3차원 유동/파동장 해석에 있어서 정확도 향상뿐 아니라, 계산영역을 확장할 수 있는 계기가 마련되었다. 또한 유동/파동 해석보다 많은 계산시간이 필요한 지형변동 해석에도 충분히 적용될 수 있다고 판단된다.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.90-98
• /
• 1998

본 논문은 2차원 해면효과의 수치계산 결과를 정리하였다. 지면으로부터의 높이변화에 따른 점성유동장을 계산하기 위하여 지배방정식으로는 비압축성 RANS방정식을, 시간에 대하여 서는 음해법으로 프로그램을 구성하였다. 압력항은 가상압축성을 도입 4차 수치확산항을 추가하는 것에 의해 계산하였으며, 높은 레이놀즈수에서의 효과적인 계산을 위해 Baldwin-Lomax 난류 모델을 도입하였다. 해면효과가 없는 무한유중에서의 NACA-0012단면 계산결과를 실험데이터와 비교하는 것에 의해 프로그램의 타당성을 확인하였다. NACA-6409와 두께비 $4.6\%$의 날개에 대하여 해면효과를 고려한 계산을 수행하였다. 높이의 변화에 따라 계산된 무차원계수, 압력 및 속도분포는 해면효과의 특성을 잘 보여주고 있다.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 1998.09a
• /
• pp.12-16
• /
• 1998

현재 쯔나미의 계산에 있어서 수치계산의 안정성이나, 파형의 재현을 고려하기 위한 기준 등은 제시된 바 있다. 대표적인 예로 CFL 안정조건과 같이 계산의 불안정을 방지하기 위한 기준을 비롯하여 한 파장내에 적어도 20개의 공간격자를 확보하여야 한다는 등의 기준(Shuto,1991)이 그것이다. 또한, 쯔나미의 거동특성 중 하나인 분산현상에 대해서도 Imamura 등(1988)은 수치분산의 영향을 추정하였고, 저차방정식으로부터 고차의 근사해를 구하는 등의 수법을 제안하기도 하였다. (중략)

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.105-105
• /
• 2021

본 연구에서는 2차원 비선형 천수모형과 수심평균된 스칼라 이송모형을 해석하는 수치모형에 대해 기술하였다. 수치모형의 정확성을 보장함과 동시에 안정성을 높이기 위해 유한체적법, 플럭스 재구성 및 minmod 제한자를 사용하였다. 비선형 천수방정식의 이송항과 바닥 경사항은 계산된 수심의 양수 보존과 흐름의 정상 상태를 보장하기 위한 second order well-balanced positivity preserving central-upwind method를 이용하여 수치적으로 이산화되었다. 마찬가지로, 이송-확산 방정식 내 이송항은 동일한 2차 풍상차분법을 통해 수치적으로 풀이하였다. 격자점 경계면에서의 불연속으로 인한 수치진동을 방지하기 위해 이송항의 계산에 포함된 보존항의 차이로 인해 발생하는 스칼라의 수치확산을 최소화하기 위해 무차원의 비소산함수를 도입하였다. 또한, 확산항은 유한차분법을 이용하여 이산화하였다. 제안된 수치모형은 시간미분항의 계산을 위해 오일러 기법을 적용하여 계산된 수심 및 스칼라의 양수 보존여부와 함께 정지된 흐름의 정상 상태의 보존여부를 확인하였다. 제안된 수치모형의 해석 정확성을 평가하기 위해 1, 2차원 공간 내 다양한 흐름 조건에서의 해석해를 이용한 3개의 벤치마크 테스트를 수행하였다. 평균 제곱근 오차(Root Mean Squared Error, RMSE)를 산정하여 수치모형의 성능을 정량적으로 평가하였으며, 비소산함수를 적용함에 따라 스칼라의 수치확산이 감소하게 되었음을 확인하였다. 또한, 세 차례의 벤치마크 테스트 결과는 공통적으로 수치모형에 의해 계산된 결과값이 비소산함수를 고려함에 따라 해석해와 잘 일치함을 확인하였다.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.3 s.33
• /
• pp.114-122
• /
• 1999

준쇄파(쇄파발생정 리플형태의 불안정파)를 수치적으로 해석하고 예측하기 위하여 수중날개를 대상으로 점성유동장의 계산을 수행하였다. Navier-Stokes 방정식을 사용하여 자유수면을 계산하였고 정도향상을 위하여 Euler 형태의 자유표면조건에 고차의 유한차문법을 적용하였다. 이산화 과정에서 자유표면 격자에 3차 풍상미분항을 적용시켜 수치계산을 수행하였고, 계산결과를 사용하여 준쇄파의 생성조건을 규명하였다.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.1
• /
• pp.65-72
• /
• 1993

A numerical method for nonlinear free-surface-wave problem is developed in this paper. The final goal of this study is to simulate the towing tank experiment of a ship model and to partially replace the experiment by the numerical model. The exact problem in the scope of potential flow theory is formulated by a variational principle based on the classical Hamilton's principle. A localized finite element method is used in the present numerical computations which made use of the following two notable steps. The first step is an efficient treatment of the numerical radiation condition by using the intermediate nonlinear-to-linear transition buffer subdomain between the fully nonlinear and linear subdomains. The second is the use of a modal analysis in the final stage of the solution procedures, which enables us to reduce the computation time drastically. With these improvements the present method can treat a much larger computational domain than that was possible previously. A pressure patch on the free surface was chosen as an example. From the present computed results we could investigate the effect of nonlinearity on the down-stream wave pattern more clearly than others, because much larger computational domain was treated. We found, specifically, the widening of the Kelvin angle and the increase of the wave numbers and the magnitude of wave profiles.