• 한국항공우주학회지
• /
• 제36권6호
• /
• pp.525-531
• /
• 2008

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.199-201
• /
• 2011

NACA0012 Airfoil was simulated with Computational Fluid Dynamics(CFD) and the aerodynamic characteristics was analyzed for various far-field boundary distances ranging from 10 airfoil chord to 50 chord Drag coefficient distribution was dependent on the far-field distance and circulation, integrated along the loop inside the flow region, was also dependent. It was turned out that some corrections based on the circulation should be added to the far-field boundary condition for accurate airfoil simulation.

• 한국음향학회:학술대회논문집
• /
• 한국음향학회 2000년도 하계학술발표대회 논문집 제19권 1호
• /
• pp.299-302
• /
• 2000

본 연구에서는 원형실린더에 의한 음향파의 산란현상을 전산공력음향학 기법을 이용하여 계산하였다. 특히 전산공력음항학에서 정확도를 위해 요구되는 좌표의 직교성을 유지하기 위해서 그에 대한 적절한 관계식을 유도하였으며 정확성의 검증을 위해서 수치적인 해를 이론적인 해와 비교, 분석하였다. 공간차분법으로는 Taylor 전개를 통하여 차 정확도를 가진 차분법을 바탕으로 주파수 공간에서 최적화 된 DRP(Dispersion Relation Preserving) 기법을 사용하였으며, 시간차분법으로는 Adams-Bashford 방법을 기준으로 최적화된 4단계 외재적(explicit) 적분방법을 사용하였다 벽면 경계조건으로는 가상점 개념을 이용한 경계조건을 사용하였으며 원방 경계조건으로서는 선형화 된 Euler 방정식의 점근해(Asymptotic Solution)을 이용한 방사경계조건(Radiation Boundary Condition)을 사용하였다.

• 대한조선학회논문집
• /
• 제29권4호
• /
• pp.114-131
• /
• 1992

본 논문은 2차원 자유표면파문제에서 시간영역해법을 이용하여 2차원 운동문제에 적용할 수 있는 수치해석을 하였다. 경계조건으로는 엄밀한 물체표면 경계조건과 비선형 자유표면경계 조건을 부과했다. 수치해를 구하는데는 코시이론을 이용하여 제2종 프레드흘름 경계적분방정식을 도출하고 이를 이산화시켜 처리하였다. 수치계산을 위해 전영역을 유한한 영역으로 제한하여야 한다. 제한된 영역에서 방사해의 부과를 위해 전영역을 수치해석영역과 외부영역으로 나누고, 외부영역의 해는 그린 제2정리를 이용하고, 선형자유표면조건을 만족하는 과도그린함수를 사용한다. 위의 그린 제2정리를 이용한 식으로 부터 초기조건, 선형 자유표면조건, 무한원방조건을 이용하여 단순화시킨 다음 포텐셜과 유동함수의 관계식으로 치환하면 비선형해와 정합할 수 있는 정합행렬을 구할 수 있다. 본 논문에서 개발한 정합방법을 이용하여 적용할 문제로서 첫째는 주상체가 상하동요, 수평동요하는 경우 계산이고 두번째로는 수면하에서 타원형실린더가 일정속도로 항진할 때 계산을 수행한 결과를 고차스펙트럴방법과 비교하였다.

• 한국항공우주학회지
• /
• 제37권3호
• /
• pp.215-223
• /
• 2009

본 연구에서는 헬리콥터 로터 블레이드의 움직임을 모사하기 위해 중첩 격자 기법을 적용하여 헬리콥터 로터의 전진 및 제자리 비행을 모사하였다. 제자리 및 무양력 전진 비행은 Caradonna & Tung의 로터 블레이드를 적용하였으며 전진 비행은 AH-1G 로터 블레이드를 적용하여 수치해석 하였다. 전진 비행 시 cyclic pitch각에 대해서 Newton-Raphson 수렴 방법으로 수치 트림을 수행하였으며 수치 트림에 의한 결과를 실험 및 다른 수치해석 결과와 비교하였을 때 실험값과 유사한 결과를 얻었다. 또한 수치 트림에 의한 결과는 로터 전진면에서 나타나는 BVI 현상을 잘 모사하였다. 지배 방정식은 3차원 비정상 오일러 방정식을 사용하였으며 원방 경계 조건으로 리만 불변치 경계조건을 적용하였다.

• 대한기계학회논문집
• /
• 제16권11호
• /
• pp.1999-2010
• /
• 1992

본 연구에서는 (1)의 방법의 문제점.즉, 두꺼운 판에 강체 균열을 열매립하는 과정에서 균열 첨단에 열변형이 일어나, 변위측정시 측정 오차가 발생하는 점을 해결 하기 위해, 매유 얇은 P.V.C.판(t=0.15mm)에 강체 균열을 열매립 하는 대신 접착제로 부착하는 방법을 채택하여 강체균열 경계상의 곡률 영향을 검토하였으며, (2)의 방법 의 해결책으로서, 강체균열에서 충분히 먼곳의 곡률은 얇은 합금공구강판(alloy tool steel plate`t=1mm)으로 제작한 프레임(Frame)으로 원방곡률을 미리 슬랩 경계조건과 일치하도록 하였고, 또한 변위도 고정(동결)할수있는 새로운 고정 하중장치를 개발하 여 실험 측정치의 정도를 높였다.

• 한국농업기계학회:학술대회논문집
• /
• 한국농업기계학회 2017년도 춘계공동학술대회
• /
• pp.59-59
• /
• 2017

무인 헬리콥터의 양력을 개선하기 위한 익형 설계 단계로서 두꺼운 익형(V1505A)과 얇고 처진 익형(V2008B)의 기본 두 익형의 특성을 예측하는데 있어 회전하는 블레이드의 현실적 조건을 반영한 3D 모델을 마련하고 성능을 예측하였다. Fluent를 이용한 400 mm 선형모델의 시뮬레이션에서는 V1505A 익형은 높은 받음각에서 안정적인 특성을 보인 반면 V2008B는 비교적 높은 동력효율 특성을 보였으나, 높은 받음각에서는 실속 이후 양력이 급락하는 특성을 나타낸다. 형성된 노드 수는 약 870,000개로 하였다. 시위길이 135 mm인 익형 V2008B의 형상은 ANSYS (Fluent v16.2)를 이용해 반경(길이) 1,502 (1,380) mm 의 로터 블레이드를 구성하였다. 충분하지 않은 유동장이 익형 표면에서의 유동의 영향에 영향을 주지 않도록 직경 20 m의 원방경계(far field)를 형성하였다. 사용된 매쉬의 형태는 정사면체 형태로 로터 표면으로부터의 첫 번째 두께 높이는 0.001 m이고 10개의 층으로 형성하였다. 정지 비행하는 헬리콥터의 상태를 가정하여 회전좌표계를 이용하여 정상상태의 유동을 해석하고 사용된 난류모델은 넓은 영역에서의 유동을 고려하여 Realizable $k-{\varepsilon}$ 모델을 사용하였다. 내측그립 받음각 $6{\sim}22^{\circ}$에 대하여 현실적인 회전속도를 연동하여 600~1000 rpm을 적용하였다. 반복수(iteration)는 2000으로 하여 잔차값(residual)이 충분히 수렴하도록 하였다. 전체적으로 실제 헬리콥터가 발휘하는 양력보다는 낮은 수치로 예측되었으며 모델 및 해석 조건에 대한 검토가 필요해 보인다. 양력 값은 받음각 $10^{\circ}$에서 자중(약 68 kgf)을 극복하였고 받음각 $12^{\circ}$에 유상하중 20 kgf을 발휘하며 888 N의 양력을 보였다, 이어 받음각 $22^{\circ}$에서 실속 현상이 발생하였다. 받음각이 증가함에 따라 항력 역시 증가하였으며 받음각 $12^{\circ}$에서 121 N이었고 실속에 이르며 항력은 갑자기 증가할 것으로 예측된다. 본 연구는 변이 익형 개발의 선행 단계로 기본 익형에 대한 공력특성을 CFD 시뮬레이션을 통하여 예측하였다. 예측 값은 현실적 실험방법을 통하여 검증이 되어야 하며 이후 변이익형에 대한 예측과 설계가 가능하다.