저자는 기존의 연구에서 대용량-비선형성을 가지는 유체의 최적화를 수행하기 위해 몇 가지 강력한 방법들을 제시한 바 있다. 즉, 최적화 과정에서 수렴성을 높이기 위해 step by step기법을 사용하였고, 또한 수렴속도를 높이기 위하여 최적화이터레이션 과정에서 얻어지는 민감도정보를 이용하여 시스템 평형방정식의 해석을 위한 좋은 초기치를 제공하는 방법과, 평형방정식을 구속조건으로 사용하는 동시기법(simultaneous technique)에서 착안하여 해석과 최적화 수렴 판정치를 조작하는 방법을 제시한 바 있다. 그러나 그들 기법은 기본적으로 유사뉴턴법에 기본을 두고 있다. 현재까지 최적화에서 SQP기법을 사용할 때는 정확한 헤시안 매트릭스의 유도가 매우 까다롭고 힘들기 때문에 유사뉴턴법을 사용하고 있는 실정이다. 그러나 3차원 문제와 같이 더욱 큰 용량의 문제를 위해서는 진정한 의미에서의 뉴턴법, 트루 뉴턴법(true Newton method)을 사용할 필요가 있다. 본 연구에서는 트루 뉴턴법을 사용하기 위해 헤시안 매트릭스의 정확치를 얻는 과정을 유도하고 이를 기본으로 트루 뉴턴법을 이용한 최적화 루틴을 만들었다. 그리고 이를 3차원 문제에 적용하여 그 효과를 검증하였다.
트랜섬 선미 형상변화에 따른 후류 유동특성을 알아보기 위해 Re=3.5×103 및 Re=7.0×103에서 수면하부 유동을 계측하였다. 선미형상변화에 따른 유동장 계측을 위해 2-프레임 그레이레벨 상호상관 PIV기법을 이용하여 분석하였다. 선미형상은 실선조사 결과를 바탕으로 선저와 이루는 각을 각각 45˚(모델 A), 90˚(모델 B) 및 135˚(모델 C)를 실험에 적용하였다. 모델의 선미 침수면의 깊이는 자유 수면으로부터 40mm를 기초하였다. 레이놀즈수가 증가함에 따라, 와의 규모가 증가하고 하류로 길게 형성되는 경향을 나타냈다.
평행한 벽 사이에 원형배관(circular pipeline)을 놓고 그 주위의 유동특성에 대한 수치연구를 수행하였다. 비압축성 유체를 가지고, Navier-Stokes 방정식을 풀었고 3차 풍상(upwind) 차분의 수치해법을 이용하였다. 한쪽 벽과의 거리가 매우 작아질 때, 볼텍스 떨어짐이 상당히 억압되는데 이것은 벽 경계와의 상호 박리 작용 때문으로 간수된다. 본 연구는 레이놀드 수의 변함과 물체가 벽에 접근함에 따른 볼텍스 떨어짐의 특성을 규명하는데 있다. 원형배관 후류와 평행벽내 유기된 박리의 상호작용을 집중적으로 다루며 서로 다른 조건에서 박리와 와역한(vorticity dynamics)의 특성을 해석하였다.
This paper presents the numerical analysis results of the viscous flow for a full ship model. The mass and momentum conservation equations are used for governing equations, and the flow field is discretized by the Finite-Volume Method for the numerical calculation. An algebraic grid and elliptic grid generation techniques are adopted for generation of the body-fitted coordinates system, which is suitable to ship's hull forms. Time-marching procedure is used to solve the three-dimensional unsteady problem, where the convection terms are approximated by the QUICK scheme and the 2nd-order central differencing scheme is used for other spatial derivatives. A Sub-Grid Scale turbulence model is used to approximate the turbulence, and the wall function is used at the body surface. Pressure and velocity fields are calculated by the simultaneous iteration method. Numerical calculations were accomplished for the Crude Oil Tanker(DWT 95,000tons, Cb=0.805) model. Calculation results are compared to the experimental results and show good agreements.