비정렬 격자계에서 형상 최적설계 기법을 이용한 유체기계의 캐비테이션 저감 연구

Abstract

본 연구에서는 유체기계에서 발생하는 캐비테이션 현상의 저감을 위한 형상 최적설계 기법에 대한 연구를 수행하였다. 개발된 설계 기법을 이용하여 2차원 익형과 3차원의 프로펠러 및 워터제트 펌프에 대한 최적설계를 수행하였으며 설계 기법의 유용성을 검증하였다. 본 연구는 크게 캐비테이션 유동해석을 위한 다상유동 해석코드 개발, 최적설계 기법 및 코드 개발 및 개발된 코드를 이용한 최적설계 응용연구의 세 가지로 요약할 수 있다. 본 연구에서는 캐비테이션 유동해석을 위한 다상유동 모사 방법으로써 균질혼합 모델을 사용하였다. 이 방법은 액체상과 기체상 간의 경계면을 무시하는 방법으로써 각 상간의 상호작용을 포착하는데에는 한계가 있으나 계산 방법이나 계산 시간의 측면에서 매우 효율적인 방법이다. 또한 매우 복잡한 형상에 대해서도 적용이 가능하고 유체기계들의 추진 성능을 비교적 정확히 예측하는 것으로 알려져 있다. 액체상과 기체상 간의 상변화는 증기압과 국소압력의 차이에 비례하여 응축율과 기화율이 결정되는 Merkle의 캐비테이션 모델을 사용하여 산출하였다. 계산 영역 내의 모든 유동장은 비압축성 유동으로 가정하여 유동해석 코드의 강건성을 확보하였다. 다상 유동장에 대한 난류 점성계수를 계산하기 위하여 혼합물의 밀도와 점성계수가 고려된 Spalart-Allmaras의 1-방정식 난류모형을 사용하였다. 액체와 기체의 혼합물이 존재하는 영역에서 난류 점성계수가 크게 예측되는 문제를 해결하기 위하여 Reboud가 제안한 보정법을 사용하였다. 개발된 유동해석 코드는 2차원에서 NACA66 수중익형에 대하여 유동해석을 수행하였고 해석결과를 실험치와 비교하여 검증하였다. 익형 표면의 압력계수가 실험치와 매우 잘 일치함을 확인하였고 증기압 영역에 대한 검증을 통해 캐비테이션의 발생 길이도 정확하게 예측하는 것으로 판단하였다. 3차원에서는 P4381 프로펠러 형상에 대하여 다양한 전진비와 캐비테이션 수 조건에 대하여 유동해석을 수행하였다. 블레이드 표면의 캐비테이션 발생 면적과 추력계수 및 토크계수를 실험치와 비교하여 검증하였으며 실험결과와 매우 잘 일치하는 것을 확인 할 수 있었다. 로터-스테이터 간의 상대운동이 있는 AxWJ-2 펌프 형상에 대한 유동해석을 추가적으로 수행하였고 그 결과를 실험치와 비교하여 검증하였다. 캐비테이션 계수에 따라 로터 블레이드에서 발생하는 캐비테이션 발생면적과 캐비테이션에 의한 전압력수두, 추력 및 토크의 성능 저하율 등이 실험결과와 잘 일치함을 확인하였다. 최적점의 탐색을 위하여 유전 알고리즘을 사용하였다. 유전 알고리즘에서 반복되는 세대 진화과정의 적합도 판정에서 과도한 계산시간이 요구되는 전산유체역학의 사용을 회피할 목적으로 대리 모델(surrogate model)을 사용하였다. 대리 모델로는 방대한 데이터베이스에 대한 함수 근사화 문제에서 우수한 성능을 갖는 것으로 알려진 방사형 기저함수 기반의 인공신경망을 사용하였다. 설계 변수들의 표본 공간에서 표본들의 고른 분포도를 확보하기 위하여 improved distributed hypercube sampling 방법을 사용하였고 선택된 표본들에 대한 계산 격자는 초기 형상에 대한 계산 격자를 격자변형 기법을 사용하여 생성하였다. 근사 함수에 대하여 탐색된 최적점에 대하여 전산유체역학 해석을 통해 정확도를 검증하여 설계완료 조건 이하의 오차를 가질 때에 설계를 완료하였다. 오차가 설계완료 조건 이상일 때에는 검증값을 데이터베이스에 추가하여 함수근사화, 최적점 탐색 및 검증의 과정을 반복하였다. 개발된 유동해석 코드와 최적설계 코드를 이용하여 유체기계의 캐비테이션 발생 저감을 위한 최적설계 응용연구를 수행하였다...