In this paper, a numerical analysis is carried out to study the drag of conical cavitators, supercavity generation devices for the high-speed underwater vehicle. The realizable k-∊ turbulence model and the Schnerr-Sauer cavitation model are applied to calculate steady-state supercavitating flows around cones of various cone angles. The calculated drags of the cones are decomposed of the pressure and the friction parts and their dependency on the geometry and the flow conditions have been analyzed. It is confirmed that the pressure drag coefficients of the cones can be estimated by a simple function of both the cone angle and the cavitation number while the friction drag coefficients approximately by well-known empirical formulas, e.g., Schults-Grunow's for the drag of the flat plate. Finally a practical method for estimating the total drags of supercavitating cones is suggested, which can be useful consequently for the design of conical cavitaors.
초공동(supercavity)을 이용한 고속화 기술은 수중운동체의 머리에 장착한 캐비테이터(cavitator)에서 발생시킨 초공동으로 몸체를 감싸서 물과 접촉하는 면적을 최소화시켜 항력을 줄임으로써 수중에서 고속으로 운항하는 기술이다. 러시아, 미국, 독일, 중국을 중심으로 초공동을 활용하여 수중에서 고속 운항하는 기술이 국방기술 분야에서 연구, 개발되어왔으며, 대표적인 사례로 초공동 어뢰(supercavitating torpedo)인 러시아의 Shkval (Kim & Nah, 2011)과 독일의 Barracuda를 들 수 있다. Shkval은 원판(disk) 형상, Barracuda는 원뿔(cone) 형상의 캐비테이터를 어뢰 머리에 장착하고 수중에서 200노트 이상으로 운항이 가능한 것으로 알려져 있는데, 일반 어뢰의 주행속도가 40노트 전후임을 감안하면 초공동 어뢰는 매우 빠른 속도로 항주함을 알 수 있다.
캐비테이터는 수중운동체의 항력을 최소화하도록 초공동을 발생시키고 유지하는 장치로서 수중운동체의 속도 성능을 결정하는 초공동 고속운항 기술의 핵심이다. 초공동 수중운동체의 개념 설계 또는 기본설계의 초기 단계에서는 목표 속도를 달성하기 위한 캐비테이터의 기본 형상과 제원을 선정해야하며 이를 위해서는 캐비테이터의 초공동 및 항력 특성을 추정할 수 있어야 한다.
현재까지 캐비테이터의 기본 형상으로는 원판 또는 원뿔이 채택되거나 제안되었는데, 이들 축대칭 형상의 캐비테이터에 대한 초공동 발생 및 항력에 대한 실험 데이터는 매우 부족하고 공개된 연구결과 또한 제한적이다. 그나마 원판과 원뿔 모형에서 발생하는 초공동 및 항력에 대한 실험 데이터 (Reichardt, 1946; Eisenberg & Pond, 1948; Self & Ripken, 1955)는 1960년대 이전에 공개된 것이며, 그 이후에는 러시아의 연구결과 (Epshtein, 1973; Logvinovich, 1975; Semenenko, 2001; Savchenko, 2002) 그리고 미국의 기초연구결과 (Brennen, 1969; Kirschner, et al.. 1995; Kunz, et al., 2001; Alyanak, et al., 2004)가 고작이다.
최근 국내에서도 초공동 수중운동체에 대한 연구 (Petitpas, et al., 2011; Ahn, et al., 2012; Kim, et al., 2013; Kim, et al., 2013; Park, 2013; Kim & Lee, 2014)가 활발히 진행되고 있지만 캐비테이터의 형상과 제원의 선정에 필요한 초공동 및 유체력 정보는 아직 부족한 실정이므로 캐비테이터의 기본 설계에 활용할 수 있는 초공동 및 항력 추정 기술의 확보가 요구된다. 특히, 캐비테이터의 항력은 캐비테이터에서 발생하는 초공동의 크기를 추정하는 점근 근사식과 준 경험식에 사용되는 핵심 정보이기 때문에 캐비테이터의 항력 추정 기술은 캐비테이터 설계의 열쇠라고 할 수 있다.
본 연구에서는 독일의 초공동 어뢰인 Barracuda에서 채택하고 있는 캐비테이터의 기본 형상인 원뿔에서 발생하는 정상 상태의 초공동 유동에 대한 수치계산 결과를 분석하고 항력공식을 도출하여 원뿔 캐비테이터의 항력을 추정하는 방법을 제안하였다.
캐비티 유동(cavity flows)을 수치해석하기 위하여 혼합유체(mixture fluid) 모델을 적용하면 지배방정식인 연속방정식과 RANS(Reynolds averaged Navier-Stokes) 방정식은 식 (1), (2)와 같이 나타낼 수 있다.
여기서
여기서
한편, 상과 상 사이의 경계는 체적분율(0.5)로 결정되며, 먼저 2차상(secondary phase)에 대한 체적분율 방정식(volume fraction equation)으로부터 2차상의 체적분율을 구하고, 식 (4)로부터 1차상(primary phase)의 체적분율을 구한다. 2차상
캐비티 유동을 해석을 위해서는 체적분율 방정식 뿐 아니라 증기에 대한 수송방정식도 풀어야 하며 증기 수송방정식은 식 (8)과 같다. 식 (8)의 우변은 기화 압력에 대한 상대 압력에 따라서 기화(evaporation) 또는 응축(condensation)하는 상변화를 모델링한 소스(source) 항이다.
여기서
[Table 1] Cavitation model constants
Cavitation model constants
Table 1에서 ℜ
초공동 유동해석을 위해 선체 주위 유동계산에 자주 사용하고있는 유한체적법(FVM) 기반의 상용코드인 FLUENT(ver. 15)를 사용하였다. 지배방정식의 확산항은 2차 중심차분, 대류항은 QUICK(Quardratic Upwind Interpolation for Converctive Kinematics)방법으로 이산화 하고, 속도-압력 연성은 SIMPLEC(Semi Implicit Method for Pressure Linked Equation Consistent) 방법을 사용하였다. 난류모형은 realizable
본 장에서는 원뿔 캐비테이터의 초공동 발생 및 항력에 대한 기본적인 특성을 파악하기 위해 광범위한 원뿔 각의 원뿔을 대상으로 축대칭 정상 초공동 유동에 대한 수치계산을 수행하고 초공동 유동의 계산결과를 살펴보았다.
Fig. 1과 같이 원뿔 베이스의 중심을 원점으로 잡고, 균일류
유동 계산을 위한 수치격자는 상용코드인 Gridgen을 사용하여 축대칭 조건을 적용한 2차원 정렬 격자를 생성하였다(Fig. 2). 유입경계는 원뿔 직경(
원뿔의 베이스 직경은 1.0 인치(″)로 동일하게 고정하고, 원뿔각은 15° ~ 180° 범위에서 총 8개를 선정하였으며, 원뿔 각이 180°인 경우는 원판이므로 편의상 원판 두께를 직경의 10%로 정했다. 또한, 원뿔 크기에 따른 캐비티의 상사성을 알아보기 위해 원뿔 각 45°의 경우에는 직경이 0.25″와 2.00″인 원뿔에 대해 추가적인 계산을 수행하였다. 원뿔 직경을 사용하여 정의한 Reynolds수는 2.53×105 이며, 캐비테이션 수(
[Table 2] Calculation conditions for supercavitating cones
Calculation conditions for supercavitating cones
원뿔 각 45°인 원뿔에 대해 캐비테이션 수 0.2인 경우의 압력분포와 유선 및 캐비티 형상을 Fig. 3에 나타내었다.
원뿔 전방의 감속 영역과 경사면을 따라 형성되는 가속 영역 그리고 캐비티 내부 압력계수(
계산결과를 검증하기 위해 캐비티의 길이(Fig. 4)와 직경(Fig. 5)을 실험 데이터 (Self & Ripken, 1955)와 비교하여 함께 나타내었으며, 캐비티 길이와 직경은 원뿔 직경으로 무차원화 하였다. 실험은 미국 미네소타 대학의 10인치 자유분사류수조(free-jet water tunnel)에서 직경이 서로 다른 5개의 45° 원뿔 상사모형을 대상으로 수행되었으며 캐비티 길이 및 최대직경에 미치는 척도효과를 조사한 것으로서 원뿔의 크기와 상관없이 기하학적으로 상사인 원뿔에서 발생하는 캐비티의 길이와 최대직경은 원뿔 직경에 비례함을 보여준다.
본 연구에서 계산된 캐비티 길이 및 직경은 대부분의 캐비테이션 수에서 실험결과 보다 작게 계산되었지만 캐비테이션 수가 0.05 이하인 경우에는 모형실험 결과와 잘 일치하고 있다. 또한 계산결과 역시 기하학적으로 상사인 원뿔에서 길이와 직경이 원뿔 직경에 비례하는 캐비티가 발생하는 것을 보여준다.
Fig. 6에는 초공동이 발생하는 원뿔의 항력계수를 실험 데이터(Reichardt, 1946; Kermeen, 1954; O’Neill, 1954)와 비교하여 나타내었으며, 원전에서 직접 확인하지 못하고 다른 참고문헌(Knapp, et al., 1970)으로부터 인용한 것임을 밝혀둔다. 항력계수는 식 (12)와 같이 항력을 정체압력과 베이스 면적의 곱으로 나누어 무차원화 하였고, 실험에서는 원뿔 높이(altitude)와 베이스 직경(diameter of base)의 비를 기준으로 형상을 정의하였기 때문에 비교의 편의상 해당하는 원뿔 각을 괄호 안에 표기하였다.
원뿔 각 90°와 180°(원판)외에는 계산과 실험의 원뿔 각이 조금씩 다르기 때문에 직접 비교하기는 어려우나 본 계산결과는 실험에 비해 항력계수를 약간 크게 추정하는 경향을 보이고 있다. 하지만 원뿔 각에 대한 항력계수의 차이와 캐비테이션 수에 의한 항력계수의 변화를 근사하게 재현하는 것으로 판단된다.
초공동 캐비테이터의 항력을 추정하는 기존의 방법은 포텐셜 유동이론에 기초한 점근 근사(asymptotic approximation) 또는 포텐셜 유동의 수치해나 실험 데이터를 사용해서 보정한 준 경험식을 사용하고 있다. 점근적 근사에 의한 항력 공식은 식 (13)과 같이 표현되며 캐비테이션 수가 작은, 즉 매우 긴 캐비티가 발생한 축대칭 및 2차원 캐비테이터에 적용할 수 있는 것으로 알려져 있다 (Semenenko, 2001; Franc & Michel, 2005).
이 식은 무한공동 상태(
무한공동 상태의 항력계수
[Table 3] Drag coefficients of disk and cones for zero cavitation number
Drag coefficients of disk and cones for zero cavitation number
한편, 초공동 원판과 원뿔의 포텐셜 유동에 대한 엄밀한 경계치 문제의 수치계산 결과를 기반으로 도출한 항력공식(Guzevsky, 1979; 1983)은 식 (15)와 같다.
항력공식 (13)과 (15)는 모두 캐비테이션 수의 1차 함수로서 절편은
각 항력공식에 사용된
Fig. 8은 식 (13) 및 (14) 그리고 식 (15)을 사용해서 추정한 원뿔의 항력계수를 실험 데이터와 비교하여 보여준다. 예상한 바와 같이 기존의 항력 공식을 사용하여 추정한 항력계수는 대체로 실험 데이터와 근사한 결과를 주며, 특히 Guzevsky의 항력공식 (15)는 원뿔 각이 작은 경우에도 실험 데이터와 잘 일치하는 것으로 보인다. 점근 근사에 의한 항력공식 (13)은 원뿔 각이 작을수록 실험 데이터와 잘 맞지 않는데, 이는
기존 항력 공식의 큰 취약점은 무한공동 상태의 항력계수
초공동 원뿔에 작용하는 항력은 경사면(slant surface)에 작용하는 압력항력 성분과 베이스면(base surface)에 작용하는 압력항력 성분 그리고 경사면에 작용하는 마찰항력의 합으로서 식 (16)과 같이 표현할 수 있다.
원뿔의 경사면에 작용하는 압력은 원뿔 각에 따라서 달라지며 발생한 공동(cavity)의 영향을 받아 변화한다. 원뿔 각이 서로 다른 8개 원뿔의 초공동 유동 계산으로부터 구한 원뿔 경사면 및 바닥면의 압력분포를 각각 Fig. 9와 Fig. 10에 나타내었다.
원뿔 경사면의 압력(Fig. 9)은 원뿔 각이 클수록 높아져 원뿔각이 180°일 때, 즉 유동에 수직인 원판의 경우 가장 높은 압력 분포를 나타낸다. 경사면의 압력분포와 길이(slant height)의 변화를 함께 고려할 때, 원뿔 각이 작을수록 순(favorable) 압력구배의 절대치가 작아짐을 알 수 있다. 또 캐비테이션 수가 작을수록 압력이 더 높아지며 베이스에 가까울수록 압력차가 커진다.
베이스 압력(Fig. 10)은 증기압으로서 일정하고 압력계수는 음의 캐비테이션 수와 동일한 값을 갖는다. 단, 원뿔 각이 작고 캐비테이션 수가 큰 일부 경우, 베이스 압력이 일정하지 않은 것은 초공동이 충분히 형성되지 못한 때문이다.
Fig. 11에는 경사면(slant surfaces)에 작용하는 압력항력 성분을 정체압력과 베이스 면적의 곱으로 무차원화 하여 캐비테이션 수에 대해 나타내었다. 캐비테이션 수가 커짐에 따라서 경사면의 압력항력계수가 선형적으로 감소하고 있다.
Fig. 11에서 볼 수 있듯이 경사면에 작용하는 압력항력 성분은 원뿔 각
식 (17)의 절편(
[Table 4] Formula of A and B coefficient (equation (17))
Formula of A and B coefficient (equation (17))
한편, 원뿔 베이스의 압력계수(Fig. 10)가 음의 캐비테이션 수이므로 베이스의 압력항력계수는 식 (18)와 같이 캐비테이션 수와 동일한 값을 갖는다.
따라서 원뿔의 압력항력계수는 식 (17)과 (18)로부터 다음과 같이 나타낼 수 있다.
Fig. 13에는 근사식 (19)와 본 연구의 초공동 유동계산으로부터 직접 구한 원뿔의 압력항력계수를 비교하여 표시하였다. 원뿔 각이 작고 캐비테이션 수가 큰 일부 경우를 제외하면 식 (19)는 본 연구의 압력항력계수를 정도 높게 근사하며 원뿔 각에 따른 항력계수의 차이와 캐비테이션 수에 의한 항력계수의 변화를 잘 재현하는 것을 알 수 있다.
원뿔 경사면에 작용하는 마찰력은 경사면과 길이가 같고 면적이 동일한 평판의 마찰력으로 근사하고자 Schulz-Grunow(1940)의 마찰 공식을 이용하였으며, 원전에서 직접 확인하지 못하고 다른 참고문헌 (Hoerner, 1965)으로부터 인용한 것임을 밝혀둔다. Schulz-Grunow의 마찰 공식은 2차원 평판에서의 마찰 공식이므로 평판의 길이와 동일한 경사면의 길이를 가진 원뿔로 확장하고 경사면의 마찰항력의
기하학적 형상이 고려된 근사식 (20)은 원뿔 각에 따른 마찰항력계수를 어느 정도 근사하게 추정하는 것으로 나타났다. 원뿔 각이 작을수록 또 캐비테이션 수가 작을수록 마찰항력의 비중이 커지며, 원뿔 각이 15°일 때, 캐비테이션 수가 0.1 이하에서는 압력항력의 30% 이상인 것으로 나타났다.
결과적으로 식 (19)와 식 (20)을 더해서 초공동 원뿔의 전체 항력계수를 식 (21)과 같이 추정할 수 있다.
식 (21)과 초공동 유동계산으로부터 직접 구한 전체 항력계수를 비교하여 Fig. 15에 나타내었다. 식 (21)은 본 연구의 초공동 유동계산을 통해 구한 항력 데이터를 기반으로 도출한 것이므로 당연히 원뿔 각에 따른 항력계수의 차이와 캐비테이션 수에 대한 항력계수의 변화에 대한 유동계산 결과를 잘 재현하고 있다.
본 연구에서 도출한 새로운 항력공식 (21)을 Guzevsky의 항력공식 (15) 및 실험 데이터와 비교하여 Fig. 16에 나타내었다.
원뿔 각이 53° 이상인 원뿔에 대해서는 본 연구에서 도출한 항력공식과 Guzevsky의 항력공식이 거의 동일한 결과를 주고 있지만, 원뿔 각이 28°인 경우에는 전자의 항력공식이 후자에 비해 항력계수를 더 크게 추정하고 있다. 한편, 두 항력공식 모두 실험 데이터와 차이를 보이며 특히 원뿔 각 127°의 경우에는 실험에 비해 최대 8% 까지 항력계수를 크게 추정하는 것으로 나타났다. 초공동 수중운동체의 초기 설계 단계에서 캐비테이터의 기본형상과 제원을 선정함에 있어서 요구되는 항력 추정의 정확도와, 또 실험 데이터가 제한적인 현재의 상황을 고려할 때, 이 정도의 차이는 크게 문제되지 않으며 항력공식 (21)은 충분히 유용한 결과를 제공할 것으로 판단된다.
결과적으로 본 연구를 통해 도출한 항력공식을 사용해서 원판 및 원뿔 캐비테이터의 전체 항력은 물론, 압력항력 및 마찰항력 성분을 간편하게 추정할 수 있게 되었다.
캐비테이터의 초공동 발생 및 항력에 대한 특성을 파악하기 위해 15° ~ 180° 범위에서 원뿔 각이 서로 다른 8개의 원뿔을 대상으로 축대칭 정상 초공동 유동에 대한 수치해석을 수행하였다.
수치해석 결과를 검증하기 위해 45° 원뿔에서 발생하는 초공동의 길이와 최대직경을 실험 데이터와 비교하였고 수치해석 결과는 실험과 전반적으로 잘 일치한다. 원뿔 각에 의한 항력계수의 차이와 캐비테이션 수에 따른 항력계수의 변화도 실험 데이터와 정성적으로 잘 일치하고 있으나, 수치해석 결과가 대체로 더 크게 나타난다.
원판과 원뿔 캐비테이터에 대한 기존의 항력공식을 비교하고 각 항력공식의 장단점을 조사하였다. 본 연구의 수치해석 결과를 기반으로 하여 새로운 항력공식을 도출하고 이를 사용해서 추정한 항력계수를 실험 데이터 및 기존의 항력공식의 추정 결과와 비교함으로써 새로운 항력공식의 타당성을 검증하였다.
결과적으로 본 연구 결과는 원판 및 원뿔 캐비테이터의 전체 항력 및 그 성분을 추정함으로써 캐비테이터의 항력 특성을 파악하는데 사용할 수 있고, 본 연구에서 도출한 항력공식은 초공동 수중운동체의 초기 설계 단계에서 캐비테이터의 기본형상과 제원을 선정함에 있어서 유용한 도구가 될 것으로 판단된다.
