받음각을 갖는 초공동 수중 운동체에서 발생하는 초월공동과 유체력 특성에 대한 수치적 연구

2.1 지배방정식

유동장의 지배방정식은 연속방정식이며, 식 (1)과 같이 나타낼 수 있다. 본 연구에서는 3차원 초공동의 수치해석을 위해 비압축성, 정상유동으로 가정하였다. 상은 물(water)과 증기(water-vapor)로 총 2개의 상으로 이루어져 있다.

여기서 $$ \overrightarrow{v}$$은 혼합유체의 속도이며, 각각의 상이 동일한 속도를 갖는 균일속도(mass averaged mean velocity)이다. ρ_m은 혼합 밀도(mixture density)이며, 2상에 대한 질량평균 혼합밀도는 식 (2)와 같이 나타낼 수 있다.

여기서 k는 혼합유체 내에 존재하는 유체의 상이며, α_k는 k상의 체적분율(volume fraction)이다.

Navier-Stokes 방정식은 식 (3)와 같이 나타낼 수 있다.

여기서 μ_m은 혼상 점성계수를 나타내며, μ_m은 식 (4)와 같이 나타낼 수 있다.

여기에서 μ_t는 난류점성계수(turbulence viscosity)이다. μ_t는 k-ω SST 난류모델을 사용하여 계산하였다. 물체 벽면 경계조건 처리는 표준 벽함수(standard wall function)를 사용하였다.

2.2 캐비테이션 모델

본 연구에서는 1차상(primary phase)을 물, 2차상(secondary phase)을 증기로 설정하였다. 물과 증기로 이루어진 다상유동(multiphase flow) 해석을 위하여 상의 경계면(phase interface)을 추적하는 방법인 Volume of fluid 방법을 적용하였다. 상과 상 사이의 경계는 식 (5)과 같이 증기에 대한 체적분율을 2차상에 대한 체적분율 방정식(volume fraction equation)을 통해 구하고, 식 (2)를 통해 1차상인 물의 체적분율을 구한다.

여기서 $$ \dot{m^{-}}$$는 상변화에 의해 액상이 기화하는 질량변화량을 의미하며, 반대로 $$ \dot{m^{+}}$$는 상변화에 의해 기상이 응축되는 질량변화량을 의미한다.

캐비테이션의 유동해석을 위해서는 체적분율 수송방적식을 풀어야 한다. 캐비테이션 수송방정식은 시스템내의 상변화에 의한 열역학 및 동역학적 항으로 구성되며, 본 연구에서는 Schnerr and Sauer (2001)의 캐비테이션 모델을 사용하였다. Schnerr and Sauer의 캐비테이션 모델은 식 (8)에서 확인할 수 있다.

여기서 R은 증기 기포의 반지름이며 식 (9)와 같이 나타낼 수 있다.

이 때 η₀는 증기의 핵밀도이며, 본 연구에서는 1×10¹³으로 설정하였다.

2.3 수치해석 모델

초공동 수중운동체의 수치해석은 상용수치해석 프로그램인 Fluent 19.2에서 수행하였다.

수중운동체 전체형상은 Fig. 1에 나타내었다. 본 연구에 활용된 초공동 수중운동체의 형상은 캐비테이터, 목(neck), 몸체, 제어핀으로 구성된다. 캐비테이터의 직경을 D_c, 전체 길이를 L_A, 몸체의 직경을 D_B, 제어핀 폭을 F_S으로 정의하였다. 제어핀의 단면은 삼각형 형상이며, 몸체에 가까이 갈수록 두께와 길이가 길어지는 Swept-wing 으로 선정하였다.

Fig. 1

External shape of an supercavity underwater vehicle

Fig. 2는 초공동 수중운동체의 유동해석을 위한 유동영역과 좌표계를 나타낸다. 캐비테이터 중심을 좌표 원점으로 잡고, 유동방향을 x축, 상방향을 z축, 운동체의 우현을 y축으로 정하였다.

Fig. 2

Definition of coordinates and domain size

유동역역 생성 시 벽효과(wall effect)와 후류효과를 방지하기 위해 격자영역을 충분히 넓게 확보하였다. 수중운동체 전체길이 L_A를 기준으로 하여 유입경계는 캐비테이터 앞면으로부터 1L_A, 유출경계는 몸체 끝에서부터 3L_A으로 설정하였다. 그 외 y축, z축의 최외각면의 거리는 제어핀 폭 F_S를 기준으로 하여 8F_S로 설정하였다.

Fig. 3은 수치해석을 위한 격자를 보여준다. 수중운동체 주위의 격자수를 다수분포시키기 위하여 27개의 유동영역으로 나누어 격자를 분포하였다. 수중운동체를 포함한 격자영역은 non-structural 격자를 분포하였으며, 격자수를 감소시킴과 동시에 정확성을 동시에 확보할 수 있도록 다면체 격자(polyhedra mesh)를 적용하였다. 또한, 그 외 나머지 영역은 structural 격자를 분포하였다. y+ 의 평균은 100 ~ 150 이며, 사용된 격자의 개수는 약 350만개이다.

Fig. 3

Polyhedral mesh in analysis domain

유동영역에 대한 경계조건으로, 유입 경계에는 균일속도 조건을 사용하였다. 유출경계에는 절대압력 조건을 사용하였으며, 나머지 모든 면은 대칭조건을 사용하였다. 수치해석에 활용한 물리조건, 경계조건, 이산화 방법은 Table 1,2,3에 요약하였다. Fig. 4는 수치해석의 경계조건을 보여준다.

Table 1

Physics

Table 2

Boundary condition

Table 3

Solution

Fig. 4

Boundary Condition

2.4 수치해석 모델 검증

본 연구의 수치해석 모델을 검증하기 위해 원판형(disk type) 캐비테이터의 단독특성 해석결과를 실험 결과를 통해 얻은 경험식과 비교하였다. 캐비테이터의 단독특성은 초공동 수중운동체의 캐비테이터를 제외한 나머지 구성부는 제외한 상태에서 속도를 0 m/s에서 100 m/s까지 10 m/s씩 올려가며 해석을 수행하였다.

원판형 캐비테이터에서 발생한 공동의 길이(L_s)와 직경(D_S) 특성을 Garabedian (1956)의 점근 근사식과 포텐셜 유동의 수치해 및 실험자료를 통해 획득한 Logvinovich (1973), Guzevsky (1983)의 경험식과 비교하여 Fig. 5와 Fig. 6에 나타내었다. Garabedian (1956), Logvinovich (1973), Guzevsky (1983)의 원판형 캐비테이터에서 발생하는 공동의 길이 및 직경 추정식은 각각 Table 4와 같다.

Fig. 5

Comparison of cavity length

Fig. 6

Comparison of cavity length

Fig. 7

Comparison of drag coefficient

Table 4

Formulas for estimating supercavity shape pf disk

각각의 식에 사용된 항력계수 공식은 Table 5에 정리하였다.

Table 5

Formulas for estimating supercavity shape pf disk

수치해석 결과 원판형 캐비테이터에 의해 발생하는 공동의 형상특성은 추정공식과 매우 유사한 것을 확인할 수 있다. 항력계수는 캐비테이션 수가 큰 구간에서는 다소 차이가 있지만, 속도가 빨라져 캐비테이션 수가 낮은 구간에서는 해석결과와 실험결과가 매우 유사한 것을 확인할 수 있다. Table 6는 수치해석을 통해 얻은 공동의 길이(L_s), 직경(D_S), 항력계수(C_D)를 Garabedian의 추정식을 통해 획득한 값과 비교하여 표로 나타낸 결과이다. 길이, 직경, 항력계수에 대한 가장 큰 오차는 각각 8%, 9%, 6.9%로 수치해석 결과와 추정공식 및 실험결과가 매우 잘 일치하는 것을 보여준다.

Table 6

Comparison between simulation and Garabedians formulas

3.1 캐비테이터 받음각에 따른 공동 특성

캐비테이터 받음각(α)에 따른 공동의 특성과 초공동 수중운동체에 작용하는 외력특성을 수치해석을 통해 분석하였다. 캐비테이터 받음각은 0°에서 -20°까지 -5°씩 증가하며 해석을 수행하였다. Fig. 9는 캐비테이터 받음각에 따라 발생하는 공동의 형상r과 윤곽도를 나타낸다. 물과 증기의 체적분율 0.5인 면을 공동의 경계면으로 정의하였다.

Fig. 9

Predicted Cavity Shape & Vapo phase according to α

캐비테이터 받음각(α)이 0°일 때, 캐비테이터에서 발달한 공동이 몸체 전체를 감싸고 있다. 제어핀의 일부도 공동에 잠겨 물과 접하는 부분(이하 액침)의 길이가 감소한 것을 확인할 수 있다. α가 0°일 때, 발생한 공동의 형상은 수중운동체 몸체의 중심으로부터 동일한 반경을 이루고 있으며, 공동 바깥으로 노출된 제어핀 액침의 길이는 제어핀 길이의 약 65%인 것을 확인할 수 있다. 공동의 후미로 갈수록 부력의 영향을 받아 +z방향으로 약간 기울어져 진행하는 것을 확인할 수 있다.

α가 커짐에 따라, 발생하는 공동의 형상이 비대칭적으로 변형된다. +z 방향의 공동 체적은 줄어들고, -z 방향의 공동 체적은 점점 증가하는 것을 확인할 수 있다. 그 결과, 초공동 수중운동체의 목(Neck)부분이 점점 공동 밖으로 노출되고, 상부 제어핀과 하부 제어핀의 액침길이가 차이가 있음을 볼 수 있다. 또한, 비대칭적인 공동의 형상으로 인하여 공동이 몸체 전체를 덮지 못한채 마감되는 것을 확인할 수 있다.

α가 -20°일 때 상부 제어핀의 액침은 약 75% 노출되어 있으며, 하부 제어핀의 액침은 50%만 노출이 된다. 이 상태에서는, 액침의 길이가 길어진 상부 제어핀의 제어력은 상승하게 되고, 액침의 길이가 짧아진 하부 제어핀의 제어력은 저하될 것으로 보인다.

Fig. 9의 단면 A는 초공동 수중운동체의 어깨부를 나타내고, 단면 B는 길이방향 95%지점을 나타낸다. Fig. 10은 단면 A,B에서 공동의 윤곽도를 보여준다. 0°에서 원형을 이루고 있던 공동의 형상이 α가 증가함에 따라 공동의 형상이 점점 비대칭적으로 변형된다. α가 -20°일 때, +z 방향의 공동 반경이 줄어들어 수중운동체의 어깨를 공동으로 감싸지 못하게 되며, 몸체의 끝부분 또한 공동 바깥으로 노출되는 것을 확인할 수 있다.

Fig. 10

Vapor phase at Section A and Section B

Table 8은 캐비테이터 받음각(α)의 절대값에 따른 초공동 수중운동체의 항력, 양력 그리고 모멘트를 보여준다. Fig. 11은 α의 절대값에 따른 항력을 압력항력과 마찰항력 성분으로 나누어 곡선으로 나타내었으며, Fig. 12와 Fig. 13은 각각 α에 따른 양력과 y축 모멘트를 곡선으로 나타내었다.

Table 8

Drag, Lifg, and Moment according to α

Fig. 11

Drag according to the α of the cavitator

Fig. 12

Lift according to the α of the cavitator

Fig. 13

Moment according to the α of the cavitator

α 절대값이 커짐 따라 항력은 감소한다. 이는, α가 커짐에 따라 공동 바깥으로 노출되는 면적이 늘어나 마찰항력이 증가하지만, 캐비테이터의 압력면적이 감소하여 줄어드는 압력항력 비율이 더 크기 때문이다. 항력이 감소하는 반면, 양력과 모멘트는 증가하는 것을 볼 수 있다. α가 증가함에 따라 캐비테이터에 력과 모멘트를 발생시킨다. 캐비테이터에 따른 항력, 양력, y축 작용하는 항력이 분산되어 +z 방향으로 작용하게 되고, 이는 양모멘트의 경향성을 통해 유체력 계수를 도출할 수 있고, 추후 제어기 설계에 활용할 수 있다.

Table 9는 캐비테이터 받음각(α)의 절대값이 증가함에 따른 초공동 수중운동체의 전체항력과 각 구성부 별 항력성분을 보여준다. Fig. 14는 α 절대값에 따른 구성부 별 항력을 그래프로 나타내었다.

Table 9

Drag according to α

Fig. 14

Drag by component according to the α of the cavitator

캐비테이터 받음각(α)이 0°일 때, 캐비테이터에 85% 이상, 제어핀에 14% 이상의 항력이 작용하며, 몸체에는 0.01% 이하의 항력이 작용하는 것을 확인할 수 있다.

α 절대값이 증가함에 따라 초공동 수중운동체에 작용하는 전체 항력은 감소한다. 이는 항력의 대부분을 차지하고 있는 캐비테이터의 유동방향 투영면적이 감소하여, 캐비테이터의 항력이 감소하였기 때문이다. 반면, 몸체에 작용하는 항력은 증가하는 것을 볼 수 있다. 이는, α가 증가함에 따라 공동의 형상이 비대칭적으로 발달하여, +z 방향의 캐비테이터 반경이 감소함으로 인해, 몸체의 일부가 공동 밖으로 노출되었기 때문이다.

α 절대값이 증가함에 따라, 제어핀에 작용하는 항력이 일정하게 유지되다가 15°일 때, 약 12%의 항력 상승이 이루어지는 것을 확인할 수 있다. α가 0°에서 10°까지는 공동이 비대칭적으로 발달하면서 상부 제어핀의 노출면적이 증가하고, -z방향의 노출면적이 감소하면서 제어핀에 작용하는 항력은 유사한 값을 가진다. 하지만, α가 15°일 때부터는 상부제어핀의 노출면적이 증가한 반면, 하부제어핀의 노출면적은 크게 감소하지 않았기 때문에 제어핀에 작용하는 항력이 증가한 것으로 보인다.

α 절대값이 20°일 때, 몸체에 작용하는 항력은 전체 항력의 0.02%로 대세에 영향을 미치지는 않는다. 하지만, 캐비테이터의 형상설계 및 초기 제어자세를 설정할 때, 캐비테이터의 받음각에 의해 몸체가 공동 밖으로 노출되어 항력에 영향을 미칠 수 있는 것과 상하부의 제어핀의 제어력이 상이할 수 있음을 고려하여야 한다.

3.2 몸체 받음각에 따른 공동 특성

초공동 수중운동체의 몸체 받음각(AOA)에 따른 공동의 형상 특성과 초공동 수중운동체에 작용하는 외력특성을 수치해석적으로 분석하였다. AOA는 0°에서 5°까지 1°씩 증가하며 해석을 수행하였다. Fig. 15는 AOA에 따라 발생하는 공동의 형상과 공동의 상을 윤곽도로 나타낸다.

Fig. 15

Predicted Cavity Shape & Vapor phase according to AOA

AOA이 커짐에 따라, 발생하는 공동의 형상이 비대칭적으로 변형된다. 캐비테이터에서 형성되는 공동의 형상은 AOA이 증가하면서 +z 방향으로 체적이 증가하고, -z 방향으로 체적이 점점 감소하는 것을 확인할 수 있다. 따라서, AOA이 증가하면서 몸체의 하부가 공동 밖으로 노출되는 면적이 늘어남과 동시에, 상부 제어핀의 액침은 감소하는 반면, 하부 제어핀의 액침길이는 커지고 있음을 확인할 수있다.

AOA가 1°일 때, 초공동 수중운동체의 하부 어깨가 공동 밖으로 노출되고, 길이 방향 85% 위치의 몸체 하부에서 공동의 형상이 마감되는 것을 확인할 수 있다. 초공동 수중운동체의 하부형상이 공동 밖으로 노출되는 면적이 늘어나고, +z방향 공동의 체적이 늘어남에 따라 하부 제어핀의 액침은 약 95%인 반면, 상부 제어핀은 30%만 노출되어 있는 것을 확인할 수 있다.

AOA가 2°일 때, 초공동 수중운동체의 목(neck)이 공동 밖으로 노출되고, 어깨에서 공동이 생성되는 것을 확인할 수 있다. 생성된 공동은 길이방향 80% 위치의 몸체 하부에서 마감하였다. AOA이 1°와 비교하였을 때, 공동 밖으로 노출된 면적이 늘어난 것으로 확인할 수 있다. 이 때, 하부 제어핀의 액침은 100%인 반면, 상부 제어핀의 액침은 약 15%만 노출된다.

AOA가 5°일 때, 초공동 수중운동체의 Neck이 50% 이상이 공동 밖으로 노출되어 있으며, 길이방향 60% 위치의 몸체 하부에서 공동이 마감된다. 하부 제어핀의 액침은 100% 노출되어 있으며, 상부 제어핀의 액침은 공동에 다 잠겨져 제어력을 완전히 상실하였을 것으로 보인다.

AOA가 증가함에 공동의 직경은 비대칭적으로 발달하게 된다. AOA가 증가함 따라 +z방향의 공동의 반경은 커지고, -z방향의 공동의 반경은 작아진다. 이는, 캐비테이터 받음각에 의해 공동의 형상이 비대칭적으로 형성될 뿐만 아니라, 캐비테이터에서 생성된 공동이 몸체 하부의 형상과 간섭됨에 따라 완전하게 발달하지 못하고 마감되기 때문이다.

AOA에 따른 공동의 형상을 길이의 관점으로 보면, 초공동 수중운동체 몸체를 중심으로 +z방향의 공동의 길이는 줄어들고, -z방향의 공동의 길이는 늘어난다. 이는, AOA이 0°일 때는 캐비테이터에서 생성된 공동이 상하부의 제어핀까지 대칭을 이루며 도달하여, 상하부의 제어핀의 액침면적이 동일하기 때문에, 제어핀에 의해 생성되는 공동의 크기가 동일하다. 반면, AOA이 증가함에 따라, 상부 제어핀의 액침은 감소하고, 하부 제어핀의 액침은 증가한다. 따라서, 상부 제어핀에 의해 생성되는 공동의 크기는 감소하는 반면, 하부 제어핀에 의해 생성되는 공동의 크기는 증가하여 -z방향의 공동의 길이가 늘어난다.

Fig. 15에서 단면A는 초공동 수중운동체의 어깨부를 나타내고, 단면 B는 길이방향 95% 지점을 나타낸다. Fig. 16은 단면 A와 B에서 공동의 상을 윤곽도로 보여주고 있다. 0°에서 몸체를 중심으로 원형을 이루고 있던 공동의 형상이 AOA가 증가함에 따라 공동의 형상이 점점 변형되어, 몸체 아래쪽의 체적은 감소하고 몸체 위쪽의 체적이 증가하는 것을 볼 수 있다.

Fig. 16

Vapor phase at Section A and Section B

Table 10은 초공동 수중운동체의 몸체 받음각(AOA)에 따른 항력, 양력 및 y축 모멘트 값을 보여준다. Fig. 17은 AOA에 따른 전체 항력과 압력항력 및 마찰항력 성분을 나타낸 것이다. AOA가 증가함에 따라 초공동 수중운동체의 저항이 크게 증가하는 것을 볼 수 있다. AOA가 0°에서 5°로 증가할 때, 전체 항력은 약 110% 증가하는 것을 확인할 수 있는데, 압력항력은 82% 마찰항력은 2000% 증가하였다. 압력항력이 증가한 이유는, AOA가 0°일 때는 유동방향과 초공동 수중운동체의 진행방향이 동일하여, 캐비테이터에만 압력항력이 작용하였지만, AOA가 증가함에 따라 몸체의 압력면적이 증가하여, 몸체에 작용하는 압력항력이 증가하였기 때문이다.

Table 10

Drag, Lifg, and Moment according to AOA

Fig. 17

Drag according to the angle of attack

AOA가 증가함에 따라 마찰향력이 비약적으로 증가하는 것을 볼 수 있다. 이는, 0°일 때는, 캐비테이터에서 발생하는 초월공동이 몸체 전체를 덮으면서 일부 노출된 제어핀에서만 마찰항력이 발생하였지만, AOA가 증가함에 따라 공동 밖으로 노출된 몸체의 표면이 늘어남에 따라 전체 마찰항력이 증가하였다.

Fig. 18은 AOA에 따른 초공동수중운동체의 전체양력과 압력과 마찰면에 의한 양력성분을 그래프로 나타낸 것이다. AOA가 증가함에 따라, 초공동 수중운동체에 작용하는 양력이 증가하는 경향을 확인할 수 있으며, 전체 양력의 98% 이상이 압력차에 의해 양력이 발생하는 것을 확인할 수 있다. 이는, 항력성분으로만 작용하던 압력항력이 초공동 수중운동체의 받음각에 의하여 양력성분으로 나뉘어졌을 뿐만 아니라, 공동 바깥으로 노출된 몸체가 물에 의해 부력을 받기 때문인 것으로 보인다.

Fig. 18

Lift according to the angle of attack

Table 11는 초공동 수중운동체의 몸체 받음각(AOA)에 따른 각 구성부에 작용하는 항력을 표로 보여주고 있으며, Fig. 19는 AOA에 따른 각 구성부의 항력을 그래프로 나타낸 것이다.

Table 11

Drag according to AOA

Fig. 19

Drag by component according to AOA

AOA가 0°일 때, 항력의 약 86%를 캐비테이터가 차지하고, 제어핀에 작용하는 항력이 약 14%를 차지한다. 캐비테이터에서 발생한 공동이 몸체 전체를 덮고 있기 때문에, 몸체에 작용하는 항력은 거의 없고, 캐비테이터 전면부와 공동 바깥으로 일부 노출된 제어핀의 액침이 전체 저항의 99% 이상을 차지한다.

AOA가 증가함에 캐비테이터에 작용하는 항력은 큰 변화가 없는 반면, 몸체에 작용하는 항력은 만배 이상 증가하였고, 제어핀에 작용하는 항력은 약 130% 증가한 것을 볼 수 있다. 몸체에 작용하는 항력이 증가하는 것은 받음각에 의해 몸체가 유동과 직접 닿음으로 인해 몸체의 압력면적이 증가하여 몸체에 작용하는 압력항력이 발생하였으며, 몸체가 공동 바깥으로 노출되는 면적이 증가하여 마찰항력이 증가하였기 때문이다. AOA이 5°일 때는, 캐비테이터에 작용하는 항력보다 몸체에 작용하는 항력이 더 큰 영향을 미치는 것을 알 수 있다.

AOA가 증가함에 따라 제어핀에 작용하는 항력이 증가하는 것은, +z방향 공동의 체적이 증가하고, -z방향 공동의 체적이 감소하였기 때문이다. 이로 인해, 상부 제어핀은 공동에 의해 덮히는 반면, 하부 제어핀은 공동 바깥으로 다 노출되었다. 또한, 좌우 제어핀 역시 공동 바깥으로 노출된 상태에서 압력면적이 증가하였기 때문에 제어핀에 작용하는 항력이 증가한 것으로 보인다.

Table 12는 AOA에 따른 각 구성부에 작용하는 y축 모멘트를 표로 보여주고 있으며, Fig. 20은 AOA에 따른 각 구성부의 y축 모멘트를 그래프로 나타낸 것이다. y축 모멘트는 AOA가 증가함에 따라 일정한 경향성을 가지지 않는다. 일반적으로 수중운동체가 받음각을 가지면, 받음각이 커짐에 따라 모멘트가 증가하는 경향성을 가진다. 하지만, 본 연구의 해석결과에서는 0°에서 2°일 때 까지는 y축 모멘트가 증가하다가, 3°일 때 급격히 감소한 뒤, 다시 AOA이 증가함에 따라 y축 모멘트가 증가하는 것을 확인할 수 있다. 또한, 구성부에 따른 영향을 살펴보면, 초공동 수중운동체의 AOA에서는 캐비테이터가 모멘트에는 큰 영향을 주지 않는 것을 볼 수 있다. 오히려, 초공동 수중운동체에 발생하는 모멘트는 몸체와 제어핀에 의해 큰 영향을 받는 것을 알 수 있다.

Table 12

Moment according to AOA

Fig. 20

Moment by component according to AOA

AOA가 0°에서 2°일 때, 받음각이 커짐에 따라 y축 모멘트는 증가한다. 이는 몸체와 제어핀이 받음각을 받음에 따라 발생하는 모멘트가 커지기 때문이다.

AOA가 3°일 때, y축 모멘트가 급격히 감소하는 것을 볼 수 있는데, 이는 초고동 수중운동체의 어깨와 몸체 아래부분이 공동 바깥으로 노출되는 면적이 커짐에 따라 항력이 커지게 되고, 몸체 아래부분에 작용하는 항력이 –y축 방향의 모멘트를 발생시키기 때문이다. 또한, AOA가 3°일 때 상부 제어핀의 액침이 공동에 90% 이상 잠기게 되어, 상부제어핀은 항력을 거의 받지 못하게 된다. 따라서 전체 제어핀에 작용하는 모멘트는 받음각이 증가하였음에도 불구하고 줄어드는 것을 확인할 수 있다.

AOA가 4°일 때부터는 초공동 수중운동체의 y축 모멘트가 다시 증가하는 것을 확인할 수 있다. 몸체는 하부에 작용하는 항력에 의해 –y축 방향 모멘트가 증가하고 있지만, 좌우 측면의 제어핀의 받음각이 커짐에 따라 전체 제어핀에 작용하는 y축 모멘트가 증가한 것으로 보인다. 몸체에 의해 감소하는 모멘트의 크기보다 제어핀에 의해 증가하는 모멘트의 크기가 더 크기 때문에, 전체적으로 y축 모멘트가 증가하는 것을 확인할 수 있다.

수중운동체의 조종성능을 분석하기 위해서는 AOA와 제어핀의 받음각에 따른 항력, 양력, 그리고 모멘트의 경향성을 통해 유체력계수를 도출한다. 일반적인 수중운동체는 AOA와 제어핀의 받음각에 따라 각 항목들이 일정한 경향성을 가진다. 하지만, 초공동 수중운동체는 캐비테이터의 받음각, AOA에 따라 공동의 형상이 변하게 되어 공동의 형상이 항력과 양력, 그리고 모멘트에 영향을 미치는 것을 알 수 있다. 제어기 설계를 위한 유체력계수를 도출하기 위해서는 초기제어 자세를 정의하고 정의된 범위 내에서 유체력계수를 도출할 필요가 있음을 알 수 있다.