선택한 단어 수는 입니다.
최소 단어 이상 선택하여야 합니다.
최소 단어 이상 선택하여야 합니다.
선택한 단어 수는 30입니다.
최대 10 단어까지만 선택 가능합니다.
최대 10 단어까지만 선택 가능합니다.
문제 정의
- 따라서, 본 연구는 다양한 해양 환경에 따라 변하는 부유체의 부가질량 계수를 보다 간편하게 추정할 수 있도록, 단순 형상 부유체에 대한 파 주파수 및 해저면 영향이 고려된 부가질량 계수를 추정할 수 있는 약산식을 개발하고자 하였다. 무한 수심에서 단순 형상 부유체의 운동 주파수에 따른 부가질량 계수는 이미 많은 연구 결과를 통해 비교적 용이하게 그 값을 찾아 볼 수 있다(Lewis, 1989).
- 본 연구는 최근 증가하고 있는 천해에서의 선박의 내항성능, 조종성 해석 및 부유체의 파랑중 운동 성능평가를 위해, 선체 해석의 기본 형상이 되는 반원형 실린더 부유체의 수직운동에 따른 부가질량 계수를 추정하는 약산식을 개발하는데 있다. 기존 연구 결과를 통해 알려져 있는 심해조건의 고주파수 영역에서 단순 형상 부유체의 부가질량 계수에, 해저면 영향으로 인한 부가질량의 변화량을 계수화한 보정계수를 개발하여 대입함으로써 각 운동 주파수와 해저면 영향에 따른 부가질량 계수를 간편하게 계산하였다.
- 여기서 k는 방사파수(Radiated wave number)를, S는 부유체의 표면을 나타내며, V는 부유체를 강제 진동시키는 속도를 나타내며, n은 단위 법선벡터를 나타낸다. 본 연구에서는 부유체의 수직운동에 기인한 부가질량과 감쇠계수를 구하기 위해 강제 수직 운동을 대입하였다.
제안 방법
- 1은 단순 반원형 부유체가 포함된 2차원 수치 파동 수조의 계산 영역을 나타내는데, 부가질량과 감쇠계수를 계산하기 위한 방사문제를 풀기위해 부유체의 강제 수직 운동을 나타내었고, 해저면의 영향을 고려할 수 있도록 평탄 해저면을 계산 영역에 포함시켰다. 2차원 수치 파동 수조를 이용하여 고주파수 영역에서 해저면 영향이 거의 없는 심해조건(h/T=30, 여기서 h는 수심, T는 부유체 흘수)의 부가질량 계수와 특정 기준 수심에서의 계수값들을 계산하여 기준값으로 삼아, 해저면 영향에 따른 부유체의 부 가질량 변화량을 비교함으로써 부가질량 변화에 작용하는 해저면의 영향을 계수화 하였다. 계산에 적용한 반원형 부유체의 형상은 폭(B)과 흘수(T)의 비가 2이고, 해저면의 효과를 고려하기 위해서 수심(h)과 흘수(T)의 비(h/T)를 변화시키면서 각 조건에 따른 부가질량 기준값을 계산하였다.
- 2차원 수치 파동 수조를 이용하여 고주파수 영역에서 해저면 영향이 거의 없는 심해조건(h/T=30, 여기서 h는 수심, T는 부유체 흘수)의 부가질량 계수와 특정 기준 수심에서의 계수값들을 계산하여 기준값으로 삼아, 해저면 영향에 따른 부유체의 부 가질량 변화량을 비교함으로써 부가질량 변화에 작용하는 해저면의 영향을 계수화 하였다. 계산에 적용한 반원형 부유체의 형상은 폭(B)과 흘수(T)의 비가 2이고, 해저면의 효과를 고려하기 위해서 수심(h)과 흘수(T)의 비(h/T)를 변화시키면서 각 조건에 따른 부가질량 기준값을 계산하였다.
- 본 연구는 최근 증가하고 있는 천해에서의 선박의 내항성능, 조종성 해석 및 부유체의 파랑중 운동 성능평가를 위해, 선체 해석의 기본 형상이 되는 반원형 실린더 부유체의 수직운동에 따른 부가질량 계수를 추정하는 약산식을 개발하는데 있다. 기존 연구 결과를 통해 알려져 있는 심해조건의 고주파수 영역에서 단순 형상 부유체의 부가질량 계수에, 해저면 영향으로 인한 부가질량의 변화량을 계수화한 보정계수를 개발하여 대입함으로써 각 운동 주파수와 해저면 영향에 따른 부가질량 계수를 간편하게 계산하였다.
- 따라서 본 연구에서는 선박 및 해양 구조물의 운동 성능 해석에 필요한 유체동역학 계수 중 고주파수 영역에서 일정한 값으로 수렴하는 부가질량 계수를 기준값으로 하여 각 운동 주파수와 수심 변화에 따른 부가질량 계수의 변화를 추정하였다.
- 반원형 부유체의 해저면의 영향에 따른 부가질량 변화량을 계수화하기 위해 각 수심비에서 부가질량 기준값을 계산하였다. 앞장에서 설명한 주파수 영역의 2차원 수치 파동 수조 기법을 이용하여 단순 반원형 부유체의 수직 운동에 의한 부가질량 계수와 방사 감쇠 계수를 각 대표 수심 조건에서 계산하고 이를 비교하였다.
- 본 연구에서 개발한 2차원 수치 파동 수조 기법의 검증을 위해 무한 수심에서 반원형 부유체의 수직운동에 따른 부가질량과 감쇠 계수를 기존 문헌값들과 비교하였다. Fig.
- 수심대 흘수비(h/T)에 따른 2차원 반원형 부유체의 수직 운동에 대한 부가질량 계수를 예측하기 위해 심해 조건(h/T=30)에서 부터 수심비가 낮은 지점까지 부가질량 계수를 계산하고 이를 기본값으로 이용하여, 각 운동 주파수에서 해저면 영향을 계수화한 보정계수를 구하고 약산식을 개발하였다. 상대적인 고주파수 영역(#)에서 반원형 부유체의 수직 운동에 대한 부가질량이 일정한 값을 유지하는 점을 이용하여, 2차원 수치 파동 수조 기법을 사용하여 Table 1에서 제시한 흘수대 수심비(h/T) 조건에서의 부가질량 계수를 계산하고 이를 기본값으로 하여 보간법을 이용하여 해저면 영향을 계수화 하였다. 식 (11)은 본 연구에서 개발한 천해에서의 수직운동의 부가질량 약산식을 나타낸다.
- 수심 변화에 따른 부가질량의 변화를 해저면의 상대적 영향력으로 계수화 시킨 보정계수를 산출하기 위해, 해저면 경계면을 포함한 2차원 수치 파동 수조 기법을 이용하여 고주파수 조건의 부가질량 계수를 기준값으로 계산하였다. 개발된 약산식을 사용하여 반원형 부유체의 각 운동 주파수에서 수심비 변화에 따른 부가질량 계수를 추정하고 그 정확성을 수치 계산 결과와 비교하여 추정치의 유효성을 분석하였다.
- 수심대 흘수비(h/T)에 따른 2차원 반원형 부유체의 수직 운동에 대한 부가질량 계수를 예측하기 위해 심해 조건(h/T=30)에서 부터 수심비가 낮은 지점까지 부가질량 계수를 계산하고 이를 기본값으로 이용하여, 각 운동 주파수에서 해저면 영향을 계수화한 보정계수를 구하고 약산식을 개발하였다. 상대적인 고주파수 영역(#)에서 반원형 부유체의 수직 운동에 대한 부가질량이 일정한 값을 유지하는 점을 이용하여, 2차원 수치 파동 수조 기법을 사용하여 Table 1에서 제시한 흘수대 수심비(h/T) 조건에서의 부가질량 계수를 계산하고 이를 기본값으로 하여 보간법을 이용하여 해저면 영향을 계수화 하였다.
- 반원형 부유체의 해저면의 영향에 따른 부가질량 변화량을 계수화하기 위해 각 수심비에서 부가질량 기준값을 계산하였다. 앞장에서 설명한 주파수 영역의 2차원 수치 파동 수조 기법을 이용하여 단순 반원형 부유체의 수직 운동에 의한 부가질량 계수와 방사 감쇠 계수를 각 대표 수심 조건에서 계산하고 이를 비교하였다. Fig.
- 경계면에 연속인 경계 적분 방정식(식 (4))을 풀기 위해 경계면을 이산화하고 각 요소의 대표값을 요소의 중심점으로 치환 하는 일정요소법(Constant panel method)를 사용하였다. 연속 경계 적분 방정식을 각 요소로 구별된 이산화된 행렬 방정식으로 변화시키고, 여기에 각 요소의 경계 조건을 대입하여 유체 방사 문제를 풀어 아래와 같은 유체 동역학적 계수를 얻었다.
- 따라서, 고주파수 영역에서 해저면 변화에 따른 부가질량 계수의 산정은 Table 1의 비교에서처럼 매우 성공적임을 알 수 있다. 하지만, 고주파수 영역에서 구한 보정계수가 중저 주파수 영역에서의 부가질량 산정에도 효과적인지 검증하기 위해, 모든 주파수영역에서 약산식을 이용하여 부가질량을 계산하고 이를 Fig. 6에서 비교하였다. 다양한 수심대 흘수비(h/T)의 변화에 따라, 반원형 부유체의 부가질량 계수는 주파수가 클수록 약산식 결과와 수치 계산 결과(경계 요소법 기반의 수치해석법)가 거의 일치함을 알 수 있다.
데이터처리
- 수심 변화에 따른 부가질량의 변화를 해저면의 상대적 영향력으로 계수화 시킨 보정계수를 산출하기 위해, 해저면 경계면을 포함한 2차원 수치 파동 수조 기법을 이용하여 고주파수 조건의 부가질량 계수를 기준값으로 계산하였다. 개발된 약산식을 사용하여 반원형 부유체의 각 운동 주파수에서 수심비 변화에 따른 부가질량 계수를 추정하고 그 정확성을 수치 계산 결과와 비교하여 추정치의 유효성을 분석하였다.
- 본 연구를 통해 개발된 반원형 부유체의 수직운동에 대한 부가질량 계수 약산식의 정확성을 검증하기 위해 경계 요소법 (BEM, boundary element method)을 사용하여 계산한 부가질량 계수와 Table 1에서 비교하였다. Table 1은 반원형 부유체의 해저 각 수심비(h/T)에 따른 부가질량의 보정 계수, 약산식을 통한 부가질량 추정 값, 수치 파동수조를 이용한 수치계산 결과 및 추정치와 계산 결과의 오차(%)를 나타내었다.
- 해저면 영향을 기준값 (Table 1의 각 수심별 수치해석 값)을 이용하여 계수화한 보정계수를 수심대 흘수비(h/T)의 함수로 개발하여 식 (12)에 나타내었다. 여기서 보정계수의 산출은 최소자승법을 이용하여 회귀분석 하였다.
이론/모형
- 또한 계산영역을 비압축성, 비회 전성, 비점성의 포텐셜 유동(Potential flow)으로 간주하고 경계 면을 랭킨소스(Rankine source)법을 이용하여 이산화한 경계 적분 방정식을 풀어서 방사 문제의 해를 계산 하였다. 개발된 수치 파동 수조 기법은 해저면이 포함된 모든 경계면 내의 계산 영역 유체를 속도 포텐셜로 나타내었고(식 (1)), 유체가 연속 방정식을 만족한다는 조건을 대입하여 계산 영역 내 지배방정식을 라플라스 방정식(Laplace equation) (식 (2))으로 하였다.
- 경계면에 연속인 경계 적분 방정식(식 (4))을 풀기 위해 경계면을 이산화하고 각 요소의 대표값을 요소의 중심점으로 치환 하는 일정요소법(Constant panel method)를 사용하였다. 연속 경계 적분 방정식을 각 요소로 구별된 이산화된 행렬 방정식으로 변화시키고, 여기에 각 요소의 경계 조건을 대입하여 유체 방사 문제를 풀어 아래와 같은 유체 동역학적 계수를 얻었다.
- 계산 유체 영역 각 경계면의 경계 조건들을 이용하여 지배방정식을 풀기 위해 2차원에서의 그린 함수(Green function) (식 (3))를 이용하여 라플라스 방정식을 경계 적분 방정식(식 (4))으로 변화시켰다.
- 약산식 산출을 위한 기준값을 구하기 위해 본 연구에서 개발한 2차원 수치 파동 수조 기법은 경계 요소법을 근간으로 경계 적분 방정식을 사용하여 부유체 주변의 자유 표면뿐만 아니라, 해저면의 영향까지 고려할 수 있다. 또한 계산영역을 비압축성, 비회 전성, 비점성의 포텐셜 유동(Potential flow)으로 간주하고 경계 면을 랭킨소스(Rankine source)법을 이용하여 이산화한 경계 적분 방정식을 풀어서 방사 문제의 해를 계산 하였다. 개발된 수치 파동 수조 기법은 해저면이 포함된 모든 경계면 내의 계산 영역 유체를 속도 포텐셜로 나타내었고(식 (1)), 유체가 연속 방정식을 만족한다는 조건을 대입하여 계산 영역 내 지배방정식을 라플라스 방정식(Laplace equation) (식 (2))으로 하였다.
- 이를 위해 2차원 수치 파동 수조(Numerical wave tank)기법을 이용하여 주파수 영역에서 방사 문제(Radiation problem)를 풀어 2차원 반원형 부유체의 수직운동에 대한 동유체력 계수를 계산하였다. 수치해석을 통해 고주파수 운동 영역의 부가질량 계수값(기준값)을 이용하여부유체의 흘수대 수심비와 주파수에 따른 부가질량 보정 계수를 보간법을 이용하여 산출하였다. 이 보정계수는 고주파수 영역에서 각 수심비에 따른 해저면 영향을 계수화한 값으로, 여기에 무한수심에서의 부가질량 계수를 대입하여 각 주파수와 수심 조건에 따른 반원형 부유체의 부가질량을 구할 수 있다.
- 약산식 산출을 위한 기준값을 구하기 위해 본 연구에서 개발한 2차원 수치 파동 수조 기법은 경계 요소법을 근간으로 경계 적분 방정식을 사용하여 부유체 주변의 자유 표면뿐만 아니라, 해저면의 영향까지 고려할 수 있다. 또한 계산영역을 비압축성, 비회 전성, 비점성의 포텐셜 유동(Potential flow)으로 간주하고 경계 면을 랭킨소스(Rankine source)법을 이용하여 이산화한 경계 적분 방정식을 풀어서 방사 문제의 해를 계산 하였다.
- 파 주파수와 해저면 영향에 따른 단순 반원형 부유체의 부가질량 계수를 추정하기 위해서는 각 조건에 따른 부가질량의 변화량을 계수화 하기 위한 기준값들을 구해야 한다. 이를 위해 2차원 수치 파동 수조(Numerical wave tank)기법을 이용하여 주파수 영역에서 방사 문제(Radiation problem)를 풀어 2차원 반원형 부유체의 수직운동에 대한 동유체력 계수를 계산하였다. 수치해석을 통해 고주파수 운동 영역의 부가질량 계수값(기준값)을 이용하여부유체의 흘수대 수심비와 주파수에 따른 부가질량 보정 계수를 보간법을 이용하여 산출하였다.
성능/효과
- 본 연구에서 개발한 2차원 수치 파동 수조 기법의 검증을 위해 무한 수심에서 반원형 부유체의 수직운동에 따른 부가질량과 감쇠 계수를 기존 문헌값들과 비교하였다. Fig. 2와 Fig. 3에서 경계 요소법을 사용하여 계산한 본 연구결과와 Vugts(1968)의 실험값, 포텐셜 이론을 바탕으로 한 Ursell(1948)의 계산 결과, 그리고 Lewis form 근사치(Vugts, 1968)를 이용한 이론값과 비교하여 본 연구에 적용한 수치파동 수조 기법의 정확성을 입증하였다.
- 0이하인 저주파수 영역에서는 추정치의 오차가 증가함을 알 수 있 었다. 개발된 약산식을 통해 반원형 부유체의 수직 운동에 따른 부가질량 계수를 고주파수 영역에서 매우 정확하게 계산할 수 있으며, 중저 주파수영역에서도 상대적으로 정확한 값(오차범위 10%이내)을 추정할 수 있음을 알 수 있다. 그러므로 본 약산식은 수심의 영향이 큰 지역에서 선체의 과도 운동 및 특성을 알기위해 스트립 이론등을 사용한 비선형 시간 영역 해석을 빠르게 수행하고자 할 때, 매 순간 변화하는 2차원 선체 단면의 부가질량을 매우 간편하게 계산할 수 있는 유용한 방법이 될 수 있다.
- 개발된 약산식을 통해 반원형 부유체의 수직 운동에 따른 부가질량 계수를 고주파수 영역에서 매우 정확하게 계산할 수 있으며, 중저 주파수영역에서도 상대적으로 정확한 값(오차범위 10%이내)을 추정할 수 있음을 알 수 있다. 그러므로 본 약산식은 수심의 영향이 큰 지역에서 선체의 과도 운동 및 특성을 알기위해 스트립 이론등을 사용한 비선형 시간 영역 해석을 빠르게 수행하고자 할 때, 매 순간 변화하는 2차원 선체 단면의 부가질량을 매우 간편하게 계산할 수 있는 유용한 방법이 될 수 있다. 이는 복잡한 해양 동역학 프로그램을 사용하지 않고도 천해에서 유체 동역학적 계수를 상당히 정밀하고 손쉽게 계산할 수 있다는 가능성을 제시한다.
- 그러므로 본 연구를 통해 개발된 반원형 부유체의 수직운동에 따른 부가질량 계수 약산식은 무차원 운동 주파수가 1.0이상일때 사용하는 것이 유효하며(오차범위 10% 이내), 운동 주파수가 1.0 이하인 저주파수 영역에서는 오차가 점점 증가함을 알 수 있다. 이는 Fig.
- 6에서 비교하였다. 다양한 수심대 흘수비(h/T)의 변화에 따라, 반원형 부유체의 부가질량 계수는 주파수가 클수록 약산식 결과와 수치 계산 결과(경계 요소법 기반의 수치해석법)가 거의 일치함을 알 수 있다. 이는 이미 Table 1에서 확인한 바와 같이 상대적으로 고주파수 영역에서 보정계수가 더 유효하기 때문이다.
- 4의 주파수에 따른 부가질량 계수의 변화에서 볼 수 있듯이 고주파수 영역에서 부가질량이 대체로 일정한 값을 가진다는 것을 이용하여 해저면 영향도를 나타내는 보정계수를 산출하였기 때문에 부유체의 수직 운동 주파수가 클수록 약산식의 정확도가 향상됨을 알 수 있다. 따라서 본 약산식은 반원형 부유체의 수직 운동에 따른 부가질량 계수를 고주파수 영역에서 매우 정확하게 계산할 수 있으며, 중저 주파수영역에서도 상대적으로 정확한 값(오차범위 10% 이내)을 제시한다.
- 또한, 계산하려는 부유체의 수직 운동 주파수가 1.0보다 작은 경우에는 추정된 부가질량 계수와 수치 계산 결과는, 본 연구에 적용한 저주파수(ω# = 0.7)에서 최대 20%정도로 나타났다.
- 이 비교를 통해 고주파수영역(#)에서 부가질량 계수는 대체로 일정한 값으로 수렴함을 알 수 있다. 또한, 고주파수 영역에서 심해 조건의 부가질량 계수는 해저면의 영향이 어느 정도 증가할 때까지 (h/T=3) 계수 변화가 크지 않지만, 본 계산의 최소 수심의 조건 (h/T=1.5)에서는 부가질량이 급격하게 증가함을 알 수 있다. 이는 부유체에 수직 운동에 작용하는 해저면의 영향이 수심이 어느 한계치를 넘어서면 급격 하게 증가함을 의미한다.
- 이 보정계수는 고주파수 영역에서 각 수심비에 따른 해저면 영향을 계수화한 값으로, 여기에 무한수심에서의 부가질량 계수를 대입하여 각 주파수와 수심 조건에 따른 반원형 부유체의 부가질량을 구할 수 있다. 무한수심에서 단순 형상 부유체의 부가질량은 이미 많은 연구 결과를 통해 알려져 있기 때문에, 본 연구를 통해 개발된 약산식을 이용하면 각 수심 조건 및 주파수에 따른 부가질량 계수를 매우 쉽게 구할 수 있다.
- 반면, 부유체의 수직 운동 주파수가 작아짐에 따라 약산식을 통한 부가질량계수는 수치 계산 결과에 비해 점점 더 큰 값을 가지게 되고 그 오차가 서서히 증가하는 것을 확인할 수 있다. 이는 해저면 영향에 대한 보정계수 산정이 고주파수에서 계산되었기 때문에 저주파수 영역에서는 그 정확도가 떨어지기 때문이다.
- 본 연구에서 개발된 부가질량 계수 약산식은 무차원 운동 주파수가 1.0이상일 때 유효함을 확인하였고, 운동 주파수가 1.0이하인 저주파수 영역에서는 추정치의 오차가 증가함을 알 수 있 었다. 개발된 약산식을 통해 반원형 부유체의 수직 운동에 따른 부가질량 계수를 고주파수 영역에서 매우 정확하게 계산할 수 있으며, 중저 주파수영역에서도 상대적으로 정확한 값(오차범위 10%이내)을 추정할 수 있음을 알 수 있다.
- 수심대 흘수비 (h/T)가 30인 경우는 부유체 흘수에 비해 수심이 30배 더 큰 경우로 심해조건을 의미하며, 비교한 결과들은 모두 고주파수 영역 (ω# =2)에서 계산한 값이다. 약산식으로 계산한 부가질량과 수치파동 수조 기법을 이용하여 해저면 경계를 포함한 경계 적분 방정식을 푼 수치해석 결과와의 오차는 최대 0.5%정도로 나타났으며 이는 약산식에 포함된 보정계수가 해저면의 영향을 정확히 반영하고 있다고 판단할 수 있다.
- 그러므로 본 약산식은 수심의 영향이 큰 지역에서 선체의 과도 운동 및 특성을 알기위해 스트립 이론등을 사용한 비선형 시간 영역 해석을 빠르게 수행하고자 할 때, 매 순간 변화하는 2차원 선체 단면의 부가질량을 매우 간편하게 계산할 수 있는 유용한 방법이 될 수 있다. 이는 복잡한 해양 동역학 프로그램을 사용하지 않고도 천해에서 유체 동역학적 계수를 상당히 정밀하고 손쉽게 계산할 수 있다는 가능성을 제시한다.
후속연구
- 본 연구에서 제시한 반원형 부유체의 수직운동에 따른 부가질량 계수의 약산식 산출방법과 유효성 확인은 차후 다양한 형상의 부유체 단면에 대한 수평 및 수직 운동에 따른 유체 동역학 계수를 간편하면서도 정확하게 추정할 수 있는 약산식 개발의 기초연구로서 그 의미가 있다고 할 수 있다.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.