• 한국해양공학회지
• /
• 제12권1호
• /
• pp.85-98
• /
• 1998

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

• 대한기계학회논문집
• /
• 제13권6호
• /
• pp.1250-1256
• /
• 1989

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

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2019년도 학술발표회
• /
• 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차원 유동/파동장 해석에 있어서 정확도 향상뿐 아니라, 계산영역을 확장할 수 있는 계기가 마련되었다. 또한 유동/파동 해석보다 많은 계산시간이 필요한 지형변동 해석에도 충분히 적용될 수 있다고 판단된다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2021년도 학술발표회
• /
• pp.105-105
• /
• 2021

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

• 대한조선학회논문집
• /
• 제30권1호
• /
• pp.65-72
• /
• 1993

본 논문에서는 비선형 자유표면파 문제에 대한 수치해법을 개발하였다. 궁극적인 물리적 모델은 수조시험에서의 조파저항 실험으로서 이를 모사하는 수치실험을 위한 수치해법을 개발하는데 그 목적이 있다. 수치해법으로서는 선형 문제에서 이미 수치방사조건의 효율성이 입증된 국소유한요소법을 비선형 문제에 응용하였다. 수치 계산과정에서는 교란원인 압력 분포면으로부터 멀리 떨어진 곳에서는 선형해를 이용하고 교란원 근처에서는 엄밀한 비선형 조건을 만족시켰다. 비선형 영역과 선형 영역 사이에는 비선형-선형 천이 완충영역을 수치 계산영역 내에 도입하여 적절히 수치 정합을 하였다. 수치 계산 과정중 각 축차 단계에서의 모드 해석을 이용하여 계산 시간을 월등히 줄일 수 있었다. 수치 계산 모형으로는 수조 실험을 모사한 자유표면 압력 분포면이 균일한 속도로 전진하는 문제를 택하였다. 본 수치계산 방법의 효율성으로 인하여 기존의 수치 계산 방법에 비해 월등이 큰 계산영역을 사용 수치계산을 수행할 수 있었다. 교란원에서 먼 하류에서 채택된 비선형 자유표면파를 선형 계산 결과와 비교하여 그동안 부분적으로 연구된 비선형파의 특성을 규명할 수 있었다. 비선형 효과에 의해 캘빈각이 증가하며, 캘빈각 근처의 파고와 파수가 증가하는 것을 확인하였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2017년도 학술발표회
• /
• pp.322-322
• /
• 2017

천수 흐름에 대한 수치해석에서 매우 작은 수심의 발생은 해가 불안정해지는 주요 원인 중 하나이며, 경사면이 잠기고 드러나는 그 전선에서 그 현상은 더욱 두드러질 수 있다. 특히, 지배 방정식이 보존형으로 기술되는 경우, 흐름률이나 생성항의 계산에서 수심에 의한 나눗셈이 불가피하므로 보존변수를 정확하게 계산하는 것이 해의 안정성을 도모하기 위한 관건이 된다. 이러한 기대에 부응할 수 있는 수치해법으로 흐름률을 정확한 계산할 수 있는 Riemann 해법을 들 수 있다. 또한, 생성항을 정확하게 계산할 수 있도록 계산 격자를 적절하게 구성하고 그 격자가 완전히 잠기지 않을 경우에 대해 물리적으로 타당하게 처리할 필요가 있다. 이 연구에서는 흐름률의 계산에 근사 Riemann 해법을 적용하여 포물면 지형을 지나는 천수 흐름에 대해 모의하였다. 1981년에 W. C. Thacker는 회전 포물면 위의 천수 문제에 대해 천수방정식의 정확해를 처음으로 유도하였다. 이 문제는 지형의 잠김과 드러남이 다수의 계산 격자에서 지속적으로 이루어지기 때문에 천수흐름의 수치 모의에서 극도로 혹독한 조건의 시험으로 알려져 있다. 회전 포물면 위의 천수 문제에 대해 근사 Riemann 해법에 따른 자료의 재구축 방법, 잠김과 드러남의 처리 등에 대해 검토하였다.

• 한국해양공학회지
• /
• 제13권3호통권33호
• /
• pp.114-122
• /
• 1999

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

• 한국콘텐츠학회:학술대회논문집
• /
• 한국콘텐츠학회 2014년도 추계 종합학술대회 논문집
• /
• pp.41-42
• /
• 2014

대형 안테나 해석에서 큰 화두는 안테나 수치 계산 속도를 높이는 것이다. 현재 인기를 얻고 있는 병렬 처리 방식인 OpenMP, MPI, CUDA를 이용하여 안테나 수치 계산을 병렬화할 경우 발생하는 단점을 제시하고, 각 병렬 처리법의 장점도 소개한다.

• Journal of Advanced Marine Engineering and Technology
• /
• 제21권4호
• /
• pp.331-355
• /
• 1997

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

• 동력기계공학회지
• /
• 제16권5호
• /
• pp.40-46
• /
• 2012

본 연구에서는 자유표면 유동문제를 효율적으로 계산하기 위한 방법으로서, 수치 파수조를 구현하여 잠수체에 의한 조파현상을 시간영역에서 다룰 수 있는 수치해법에 대하여 소개하였다. 그리고 이를 이용하여 수행된 연구내용들을 검토하고 양력물체의 경우를 포함하기 위하여 개선된 수학적 정식화 및 수치해법의 개요를 다루었다. 임의의 운동을 하는 양력물체에 의한 조파현상을 전산기로 구현하는 수치 파수조는 과중한 계산시간이 문제가 되는데, 이는 수치 Kutta 조건의 구현과 양력표면 후방의 wake 영역을 계산에서 고려해야 하기 때문이다. 따라서 한층 더 수치계산의 효율성이 중요하다고 판단되므로, 본 연구에서 소개된 3차원 고차 스펙트럴/경계요소법(High-Order Spectral/Boundary Element Method)은 자유표면 요소수를 N이라 할 때 그 계산량이 NlogN에 비례(N이 클 때는 거의 선형적으로 비례)하여 증가하므로 기존의 방법들 보다 매우 효율적인 수치해법이라 할 수 있다. 향후, 본 연구의 타당성 검증을 위한 수치코드의 개선과 여러 가지 수치계산결과 비교 등의 노력이 더 필요하다고 생각된다.