• 대한조선학회논문집
• /
• 제31권3호
• /
• pp.80-86
• /
• 1994

자유수면을 항주하는 선박에 의하여 발생되는 비선형 조파현상을 해석하기 위한 수치해석법을 개발하였다. 유동은 비점성, 비압축성으로 가정하고 선체 및 자유수면 형상과 일치하는 좌표계의 생성을 위하여 타원형 편미분방정식을 수치해석하여 물체적합 좌표계를 생성하였으며 변환된 정규격자 물체적합 좌표계에 대한 Euler방정식을 유한차분법(Finite Difference Method)을 이용하여 계산하였다. 수치해석을 위하여 시간에 대한 미분항은 전진차분, 공간에 대한 미분항은 중심차분법으로 이산화하였고 대류항에는 수치계산의 안정을 위해 인위적인 소산(dissipation)항을 첨가하였다. 자유수면의 형상은 매 시간 단계마다 자유수면 경계조건들을 만족하도록 다시 계산되었고 격자점들은 자유수면형상의 변화에 적합하게 다시 생성되도록 하였으며 압력에 대한 Poisson방정식은 반복연산법에 의하여 풀고 그 결과를 이용하여 속도를 외삽하였다. 개발된 수치해석법의 검증을 위해 수식선형인 Wigley 모형에 대한 계산을 Fn=0.250-0.408에 대하여 수행하고, 그 결과를 실험 결과와 비교하여 잘 일치함을 보였다.

• 대한조선학회논문집
• /
• 제33권2호
• /
• pp.127-138
• /
• 1996

고속선의 내항성능 해석 방법으로 스트립 방법, 통합이론(Unified theory), 3차원 판넬 방법등이 널리 사용되고 있다. 스트립 이론은 2차원적 해석 방법으로 전진속도가 빠른 경우나 추파중 저주파수 영역에서 유체력 계수와 운동응답이 정확하지 않으며, 통합이론은 내부영역에서는 2차원적 해석을 사용하고 외부영역은 세장체 이론을 사용하여 해석하는 방법으로 수학적으로 복잡한 단점이 있다. 3차원 판넬 방법을 이용한 해석은 계산 시간이 오래 걸리지만 고속인 경우에 모든 주파수 영역에서 정확한 해를 주는 것으로 알려져 있고, 전산기의 급속한 발달과 더불어 가장 권장되는 방법이다. 본 논문에서는 고속선의 해석에 전통적으로 사용되고 있는 스트립 법에 의한 해석과 3차원 판넬 방법을 이용한 해석법을 비교한다. 3차원 판넬 방법에 의한 해석은 쏘오스 분포법과 전진하면서 동요하는 그린 함수를 사용하고, 그린 함수의 수치, 계산은 Hoff의 방법을 이용하였고, 그린 함수는 종축에 대한 대칭 관계를 이용하여 계산 시간을 줄였다. 계산에 사용된 선종은 카타마란(Catamaran) 형태의 고속선이며 상기 두 방법에 의해 구해진 유체력 계수, 파강제력과 주파수 응답함수 등을 비교하였고, 또한 불규칙파중 운동응답의 계산 결과를 비교 검토해 보았다.

• 한국수자원학회논문집
• /
• 제55권4호
• /
• pp.257-266
• /
• 2022

국지성 호우 및 설계빈도 이상 강우의 증가로 침수피해가 매년 증가하고 있으며 이에 따라 홍수 조절 및 방어를 위한 수공구조물의 중요성이 증가하고 있다. 수공구조물은 목적과 성능에 따른 설계가 이루어지고 있고 홍수량이 중요한 산정 요소이나 국내에서는 관측자료의 신뢰성 부족 및 데이터의 부족으로 인하여 수공구조물 설계를 위한 수문해석 입력자료로 사용되는 설계강우량은 정확한 확률강우량의 산정과 시간분포가 중요한 요소로 작용한다. 실무에서는 Huff의 4분위 방법의 누가우량백분율을 이용하여 설계강우량의 시간분포 회귀식을 산정하고 있으며 분위별 곡선에 대한 회귀식은 전반적으로 정확도가 높게 나타나는 6차 다항회귀식을 일률적으로 사용하고 있다. 본 연구에서는 실무에서 일반적으로 설계강우량의 시간분포를 위해 사용하고 있는 Huff의 4분위 방법의 누가우량백분율을 이용하여 통계 모델링에서 간결함의 원리에 따라 변수선택법을 이용하여 시간분포 회귀식을 유도하였으며, 유의성 검정을 통한 시간분포 회귀식의 검증을 실시하였다. 변수선택법과 유의성 검정을 통한 시간분포 회귀식 산정 결과 전진선택법과 후방제거법의 장점을 모두 가지고 있는 단계선택법을 이용하여 시간분포 회귀식을 유도하는 것이 가장 적합한 것으로 분석되었다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 추계 학술대회논문집
• /
• pp.7-12
• /
• 1998

Optimized high-order compact(OHOC) schemes were proposed, which have high spatial order of truncation and resolution to simulate the aeroacoustic problems due to unsteady compressible flows. Generally, numerical schemes are categorized explicit or implicit by time-marching method. In this research, OHOC differences which were developed with explicit time-marching method is used to have implicit formulation and the implicit OHOC differences result in block hepta-diagonal matrix. This paper presents the comparisons between the explicit and implicit OHOC schemes with a second order accuracy of time in the 1-d linear wave convection problem, and between the explicit OHOC scheme of 4th-order accuracy in time and the implicit OHOC scheme of 1st-order accuracy in tine for the 1-d nonlinear wave convection problem. With these comparisons, the characteristics of implicit OHOC scheme are shown in the point of CFL number.

• 한국군사과학기술학회지
• /
• 제9권2호
• /
• pp.152-160
• /
• 2006

Time marching methods are well-suited for high speed compressible flow computations. However, it is well known that the time marching methods suffer a slow down in convergence due to disparity in Eigenvalues. A local preconditioning method is one of numerical methods to enhance convergence characteristics of low mach number flows by modifying Eigenvalues of the governing equations. In this paper, the local preconditioning method of Weiss is applied to a 2 dimensional Navier-Stokes code and the efficiency of the preconditioning method is shown through a number of computational examples.

• 한국해양공학회지
• /
• 제33권3호
• /
• pp.229-235
• /
• 2019

In this study, a two-dimensional oscillating foil with forward speed in a propagating wave flow field was considered. The time-mean power to maintain the heaving and pitching motions of the foil was analyzed using the perturbation theory in an ideal fluid. The power, which was a non-linear quantity of the second-order, was expressed in terms of the quadratic transfer functions related to the mutual product of the heaving and pitching motions and incoming vertical flow. The effects of the pivot point and phase difference among the disturbances were studied. The negative power, which indicates energy extraction from the fluid, is shown as an example calculation.

• 대한조선학회논문집
• /
• 제28권2호
• /
• pp.77-94
• /
• 1991

일정속도로 항주하는 배에 치하여 생성되는 파를 유한차분법을 이용하여 수치적으로 시뮬레이션 하였다. 지배방정식으로는 3차원 운동량방정식(Euler Eq.)과 연속방정식을 사용하였으며, 직교 좌표 계상에서의 수치계산은 FTCS(forward time/central space)로 차분화하고 시간전진 방법(time marching scheme)으로 수행하였다. 좌표계로 staggered variable mesh system을 채용하여 기존의 계산조직을 개선하므로써, 유체 유동의 변화가 심한 선체주위에 계산 정도를 높이기 위하여 cell을 집중시켰다. 개선된 전산프로그램을 이용하여 Wigley, Series 60($C_{b}$=0.6), 그리고 MS416B 산적화물선 선형등을 대상으로 수치계산을 수행하였으며, 그 결과들을 실험값과 비교하였다.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2004년도 추계학술대회논문집
• /
• pp.938-941
• /
• 2004

Sound generation and radiation from the duct-rotor system are calculated numerically. The wake geometries of a two-bladed rotor are calculated by using a time-marching fiee-wake method without a non-physical model of the far wake. Acoustic free field due to a rotating rotor is obtained by Lowson's equation. Using Kirchhoff source, rotating sources are modeled as stationary ones and can be inserted in the thin body boundary element method. The Kirchhoff source is validated through calculation of acoustic pressure due to a rotating point force. The thin body boundary element method (thin body BEM) is validated through calculation of acoustic radiation of ducted dipole. Using Kirchhoff source and thin body BEM, acoustic radiation of a ducted rotating source is calculated. Acoustic shielding is observed by inserting a duct and shows different phenomena at each major frequency. Acoustic radiation of a real duct-rotor system is also calculated using this method and the ducted acoustic field is significantly different from rotor only.

• 한국추진공학회지
• /
• 제16권5호
• /
• pp.37-46
• /
• 2012

로켓엔진 연소실에서 막냉각을 고려한 열전달 해석을 위한 모델을 수립하였다. 연소실 내에서의 대류, 복사, 막냉각 효과는 로켓엔진을 위한 경험식을 사용하였으며, 벽 외부는 일반적인 대류, 복사 식을 적용하였다. 또한 벽 내부는 유한체적의 격자를 구성하여 시간전진법에 의해 비정상 열해석을 수행하였다. 연소압 50 기압의 액체산소/케로신 엔진을 예로 해석하였으며, 이 경우 4% 유량의 막냉각까지는 벽온도를 낮추는데 크게 기여하였으나, 그 이상의 공급은 막냉각 유량에 비해 효과가 미비함을 확인하였다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2002년도 제18회 학술발표대회 논문초록집
• /
• pp.59-59
• /
• 2002

램제트 추진기관은 압축과정을 별도의 부품 없이 형상에 의해서 감속하여 연소 압력비를 얻는다. 따라서 구동 마하수와 형상에 의해 흡입과정의 압축 효율이 결정된다. 설계점은 충분한 유량을 확보 할 수 있는 유량과 충격파 각을 조절하여 전압력 손실을 줄이도록 고려되어야 한다. 또한 연소가 일어나면 연소실 압력이 배압으로 작용하고 비행시에 받음각은 변하므로 이에 따른 성능 분석도 고려 되어야 할 사항이다. 본 연구는 국내에서 실험한 형상에 대해 수치계산을 수행하여 코드의 검증과 아울러 램제트 유동장의 수치적 시뮬레이션도 설계단계에서 하나의 도구로 이용할 수 있음을 보여준다. 실험에서는 배압 조건을 얻기 위해 유동 블록키지를 유로 내에 두어 상응하는 배압을 얻었지만 본 계산에서는 압력 경계조건을 직접 사용하였다. 유동이 비정상 특성을 가지므로 시간 정확도를 이차로 가지도록 이중시간 전진법을 사용하였다. 사용한 압력비는 충격파가 카울 끝에 닿는 임계상태에 가까운 12, 13, 14에 대해 계산을 수행하였고 부스터모드로 흡입구 끝이 막혀 있다가 램제트 모드로 바뀌어 연소실 압력이 위의 압력비라고 가정할 때의 비정상 천이 과정을 계산해 보았다. 본 계산은 흡입구 부분만을 떼어놓고 적절한 가정 하에 수행되었지만 연소실 내부도 비정상 특성을 가지므로 흡입구와 연소실을 동시에 같이 계산해야한다. 추후에 전체적인 계산을 진행하기 위한 전 단계로 흡입구 계산만을 수행하여 실험과 잘 일치하는 계산 결과를 얻었고 전체 계산을 위한 연구는 진행 중에 있다.