해양환경에 대한 연구를 수행할 때, 일반적으로 유체 밀도는 단일하다고 가정하지만, 실제 해양 내부는 여러 자연현상의 영향으로 수온이나 염분이 일정치 않아 밀도차가 존재하게 된다. 이때, 밀도가 서로 다른 유체층 사이에서 파가 발생하게 되는데 이를 통상 내면파(Interfacial wave)라고 한다. 즉, 해양내부 경계층에서 생성되는 파를 의미한다. 지구상 여러 곳에서 해양 내면파가 관측되고 있는데, 지브롤터 해협에서는 수십킬로미터의 파장을 갖는 내면파가 관측되기도 하였고(Alpers and Violette, 1993), 남중국해 동쪽 연안에서 SAR(Synthetic aperture radar)를 이용한 내면파 관측이 보고되었다(Liu et al., 1998). 특히, 우리나라에서도 15m 파고를 갖는 내면파가 동해안에서 관측되었고 (Kim et al., 2001), 제주도 북동쪽에서도 내면파가 관측 보고 되었다(Kim and Won, 1999). 해양 내면파는 표층에서는 보이지 않지만, 파의 규모가 매우 커서 여러 해상 사고를 유발할 가능성이 충분하다. 예를 들어, 미국 핵잠수함 트레셔(USS Thresher)와 스콜피온(USS Scorpion)의 갑작스런 침몰 사건의 경우 아직까지 여러 가설들이 존재하는데, 그 중 내면파에 의한 영향이라는 가설이 유력하게 추측 되었다(Pinet, 1992). 이외에도, Osborne et al.(1978)은 내면파에 의해 심해시추선의 계류선 장력이 변하는 것을 조사하였다.

그러나 실제 해양 내면파의 경우 표면파와 달리 직접 눈으로 관측하기 어려우며, 내면파에 의한 부유체의 직접적인 피해나 구조물에 작용하는 영향을 명확히 알기 어렵고, 해저지형이나 해류 등의 영향을 받아 균일한 내면파의 특성을 파악하기 어렵다. 이에 따라, 해양 내면파의 물리적 특성과 영향을 파악하기 위해, 밀도가 상이한 유체층에 대한 여러 연구들이 수행되었다. Yeung and Nguyen(1999)는 복층 유체에서 삼 차원 물체에 대해 방사 및 산란문제를 풀기 위해 그린함수를 정의하여 경계 적분식을 개발하였고, 이후 유사한 연구가 꾸준히 진행 되었다(Ten and Kashiwagi, 2004; Kashiwagi et al., 2006; Kim and Koo, 2010).

유체가 밀도에 따라 복층으로 존재할 때, Barotropic 모드와 Baroclinic 모드라는 서로 다른 특징을 가지는 두 개의 해양파 모드가 존재한다(Yeung and Nguyen, 1999). 따라서 복층 유체에 대한 해양파 분산 관계식을 통해, 한 주파수에 각각의 모드에 해당하는 두 개의 파장을 얻을 수 있다. 이때, 일반적으로 Barotropic 모드에서의 파장이 Baroclinic 모드에 비해 길다. 두 해양파 모드의 주요 특징은 Fig. 1에 나타낸 것처럼, Barotropic 모드는 표면파고가 내면파고보다 더 크고 두 파의 위상이 같은 반면, Baroclinic 모드는 내면파고가 더 크고 두 파의 위상차는 180도 이다.

본 연구에서는 경계요소법을 기반으로 자체 개발한 주파수영역의 수치조파수조를 이용하여, 복층 유체가 존재하는 상태에서 부유체의 상하동요에 따른 방사문제를 해석하였다. 방사문제 계산을 위해서는 적절한 방사경계 조건의 개발과 적용이 필수적인데, 기존 연구에서는 Barotropic 모드에서의 파수를 이용하여 수치조파수조 끝단의 방사 경계조건을 설정하였다(Kim and Koo, 2010). 그러나 부유체가 복층 유체 내에서 비교적 높은 주파수로 운동을 하게 되면 생성파는 일반적으로 Barotropic 모드가 우세하여 Baroclinic 모드는 거의 나타나지 않지만(Koo and Min, 2015), 낮은 주파수로 운동하게 되면 앞서 언급한 두 해양파 모드가 모두 존재하게 된다. 따라서 어느 한 모드의 파수로 방사조건을 설정하는 기존의 방사조건은 만족할만한 결과를 도출하는데 한계가 있다.

본 연구에서는 이러한 문제점을 보완하기 위해, 주파수 영역에서 복층 유체가 존재하는 수치파동수조에 적용할 수 있는 인공감쇠영역 수치 모델링을 새롭게 개발하여 적용하였다.

수치 모델링에서 인공감쇠영역의 설치와 적용은 보통 시간영역 수치조파수조 계산에서 사용되지만, 이번 연구에서는 주파수 영역에 적용할 수 있는 수치 해석 모델링을 개발하여, 유체 영역내의 반사파 감쇠효과를 알아보고, 기존 계산결과(Ten and Kashiwagi, 2004)와 비교 검증하였다. 주파수영역에서 방사문제 해석을 통해, 복층 유체의 수면에서 운동하는 부유체의 부가질량 및 조파 감쇠계수를 각 유체층의 밀도비 변화에 따라 그 차이를 조사하였다.

수치 계산 모델링을 위해, 각 유체는 이상유체로 가정하고 지배방정식으로는 라플라스 방정식을 사용하였다. 또한, 계산 영역 내 해저면은 편평하다고 가정하였다.

본 연구에서는 밀도가 다른 복층 유체에서 부유체 상하동요에 의한 방사문제를 주파수 영역에서 해석하였다. 계산 결과의 검증을 위해 Ten and Kashiwagi(2004)와 Kim and Koo(2010)의 결과와 비교하였다. 이때 부유체는 앞선 수렴도 검사에서 진행된 모델과 동일하다.

Fig. 5와 6은 표면파고와 내면파고의 정상상태(Steady state)를 정 위치별로 나타낸 그래프로써, 가로축과 세로축은 각각 Barotropic mode에서의 파장과 부유체의 운동 진폭 크기(η₀)로 무차원화 하였으며, 기존 방사조건을 적용한 결과(Damping X: Kim and Koo, 2010)와 인공감쇠영역을 적용한 본 계산 결과(Damping O)를 비교하였다. 이때, 모든 계산 조건은 밀도 비(γ = 0.75), 수심 (h₁/h = 0.5)을 동일하게 두고 계산을 진행하였다.

기존 방사조건은 Barotropic 모드에서의 파수를 이용하여 진행파를 방사경계면에서 수치적으로 통과시키는 방법으로 방사조건 경계식 (∂ϕ^(m)/∂x = ikϕ^(m))을 사용하였다.

부유체가 낮은 주파수로 운동하는 경우(Fig. 5), 상대적으로 긴 파장을 갖는 Barotropic 모드의 생성파와 상대적으로 짧은 파장을 갖는 Baroclinic 모드의 생성파가 혼합되어 진행하게 된다. 그러나 인공 감쇠영역에서 진행파를 모두 소멸시켜 계산 영역 양쪽 벽에서 반사되지 않기 때문에, 재 반사파로 인한 결과의 간섭 없이 각 수면 위치별로 규칙적인 파고 변화 형태를 보이는 반면, 기존 방사조건을 적용한 경우는 파고 측정 위치에 따라 불규칙한 파고 형태를 나타낸다. 이는 다시 말해 생성된 파가 계산 유체영역 외부로 완전히 빠져 나가지 못하고 영역내로 재유입 되어 생긴 반사파들의 영향으로 볼 수 있다. 반면에 부유체가 고주파수 운동 시(Fig. 6), 두 방사조건에 따른 계산 결과가 거의 일치하는 것을 확인 할 수 있다. 이는 부유체 운동 주파수가 증가할수록 생성파의 Barotropic 모드가 우세해져, 기존 방사경계 조건식도 유효하다고 판단된다. 인공감쇠영역을 적용한 결과는 계산영역 끝에서 파고가 0으로 수렴함을 알 수 있다. 이는 진행파가 감쇠영역 내에서 효과적으로 제거됨을 보여준다. 따라서 본 연구에서 개발된 인공감쇠영역은 부유체 운동 주파수 변화와 무관하게 반사파의 감쇠가 복층 유체 내에서 효과적으로 이루어짐을 알 수 있다.

Fig. 7은 부유체 운동에 의한 부가질량 및 방사감쇠계수를 기존 계산 결과(Ten and Kashiwagi, 2004)와 비교한 결과이다.

부유체 운동 주파수가 큰 영역에서는 모든 비교 결과들이 대체적으로 잘 일치하는 것을 확인할 수 있다. 이는 Fig. 6에서 확인한 바와 같이, 고주파수 운동에서 Barotropic 모드가 우세한 파가 생성되기 때문에 이를 적용한 기존 방사 경계조건도 효과적임을 의미한다.

그러나 운동 주파수가 작을수록, 기존 방사조건을 적용한 결과는 상당한 차이가 발생함을 확인할 수 있다. 이는 부유체의 저주파수 운동 시 발생하는 파가 Barotropic 모드와 Baroclinic 모드가 혼합되어 있기 때문에, 기존 방사 경계조건식으로는 진행파를 유체 영역 경계면에서 완전하게 통과시키지 못하기 때문이다. 하지만, 인공감쇠영역을 적용한 계산 결과는 부유체 운동 주파수와 무관하게 Ten and Kashiwagi(2004)의 결과와 잘 일치함을 확인할 수 있다.

그러나 주파수가 0에 가까워지게 되면, 수치조파수조의 계산영역이 매우 커지게 되어, 본 연구에서 적용한 인공 감쇠영역의 역할과 무관하게 계산 값의 신뢰성이 떨어질 수밖에 없는 한계점이 존재한다.

Fig. 8은 인공감쇠영역을 적용한 상태에서, 두 유체의 밀도 차에 따른 부유체의 상하동요 부가질량과 방사감쇠계수를 비교하였다. 부가질량과 감쇠계수는 밀도비(γ)가 작아질수록, 즉 유체층간 밀도차가 클수록 상대적으로 작은 값을 나타낸다. 이는, 부유체가 상부 유체에 위치하고 있어 상부유체의 밀도가 하부 유체에 비해 상대적으로 작아지면서, 부유체 운동이 상하부 유체에 미치는 영향이 변화됨을 의미한다. 부유체 운동 주파수가 증가함에 따라 부가질량은 점점 증가하는 반면, 방사감쇠 계수는 감소하여 밀도 비에 의한 차이가 점차 작아지는 것을 볼 수 있다. 부가질량의 경우 저주파수 영역에서 주파수에 따른 변동이 생긴 후 계속 증가하는 모습을 보이는데, 이 변동이 일어나는 주파수 위치가 밀도비가 증가할수록 점점 작아져서 단일 층에서의 특징과 점차 비슷해지는 것을 확인 할 수 있다. 방사감쇠 계수의 경우도 저주파수 영역에서 약간의 변동을 보인 뒤 계속해서 감소하는데, 밀도비가 증가할수록 변동이 생기는 주파수 위치가 작아지고 점점 단일 층의 결과로 수렴함을 확인할 수 있다.

Fig. 9는 전체 수심(h)은 고정시킨 채, 상부 유체의 수심(h₁)의 변화에 따른 부가질량과 방사감쇠 계수를 나타내었다. 저주파수 부근에서, 부가질량은 상하로 변동한 뒤 일정 값에 수렴하고, 방사감쇠 계수는 특정 주파수 근처에서 상승했다가 서서히 감소한다. 이러한 특징은 상부유체의 깊이가 얕을수록 더 명확하게 나타나고, 이는 부유체와 두 유체의 경계면과의 거리가 상대적으로 가깝기 때문에 복층 유체의 특성을 더 잘 보여주는 것으로 보인다. 반면, 상부유체의 깊이가 일정 깊이 이상으로 증가하면 결과 값의 차이가 크지 않은 모습을 보인다. 또한, 주파수가 증가함에 따라, 부가질량이나 감쇠계수 모두 수심 변화에 따른 차이가 작아지는 것을 확인 할 수 있다.

Fig. 8과 9에서 확인된 저주파수 부근에서의 부가질량과 감쇠계수의 변화는 두 해양파 모드가 혼합되어 진행하는 방사파의 특징으로 인해, 단일 층에서의 결과와 다른 특징을 보여주는 것으로 판단된다. 반면, 부유체 운동 주파수가 증가함에 따라 Barotropic 모드가 우세해져, 계수 값들의 변화는 상대적으로 완만하게 수렴하는 것을 확인할 수 있다.

본 연구에서는 경계 요소법 기반의 주파수 영역 이차원 수치조파수조 기법을 이용하여, 복층 유체가 존재하는 계산영역에서 부유체의 상하동요시 발생하는 방사문제를 해석하였다. 부유체 운동으로 생성되는 파를 효과적으로 감쇠하기 위해, 주파수영역 계산에 적용할 수 있는 인공감쇠영역 경계 조건식을 개발 및 적용하여 그 효과를 검증하였다. 부유체는 이차원 Lewis 형상을 이용하였고, 상하동요에 따른 부가질량과 방사감쇠계수를 각각의 방사 경계조건에 따라 비교하였다. Barotropic 모드의 파수를 이용한 기존 방사 경계조건은, 부유체가 고주파수 운동시에는 효과적이지만, 저주파수 운동시에는 그 유효성이 제한적이다. 반면에, 인공감쇠영역을 적용시킨 결과는 부유체 운동주파수와 무관하게 계산 영역내 반사파를 효과적으로 제어할 수 있음을 확인하였다. 또한, 두 유체의 밀도비와 전체 수심에 대한 상부유체의 수심비에 따른 부유체의 부가질량과 방사감쇠계수의 변화를 확인하였다.

각 유체의 밀도차가 클수록 부가질량과 감쇠계수는 대체로 감소함을 알 수 있었고, 저주파수 부근에서 변동 진폭이 증가하였다. 그러나 이러한 차이는 주파수가 증가할수록, 상부유체의 수심이 깊어질수록 변동폭이 감소하였다. 이 같은 이유는 부유체가 낮은 주파수 운동시에는 두 해양파 모드가 혼합된 방사파가 생성되며, 높은 주파수 운동시에는 Barotropic 모드에서의 생성파가 지배적이기 때문으로 판단할 수 있다.

본 연구에서 개발 검증된 인공감쇠영역 경계조건을 적용한 수치모델링 방법은, 밀도가 다른 복층 유체에서 다양한 조건하에 생성되는 진행파의 모드별 특성을 주파수 영역 해석을 통해 정밀하게 파악하고 분석하는데 유용한 방법이 될 것이다.