• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2010년도 제35회 추계학술대회논문집
• /
• pp.363-366
• /
• 2010

여러 형태의 풍동장치 중에서 Ludwieg Tube(LT)는 높은 레이놀즈수 유동을 실현할 수 있는 가장 적합한 장치로 알려져 있다. 일반적으로 LT는 격막의 위치에 따라 두 종류로 구분된다. 본 연구에서는 상류막 방식과 하류막 방식 LT의 작동특성을 조사하기 위하여, 수치해석적인 방법을 적용하였다. 수치해석에서는 2차원 축대칭, 비정상, 압축성 Navier-Stokes 방정식에 유한체적법을 적용하였다. 그 결과, 본 수치해석은 시동과정의 유동 메커니즘을 파동선도를 이용하여 상세하게 설명하였으며, 시동 시간과 작동시간의 특성을 조사하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1999년도 춘계 학술대회논문집
• /
• pp.98-105
• /
• 1999

The Pressure-based Unstructured-grid Navier-Stokes Solver PUNS-2/3D for incompressible steady and unsteady viscous flows has been developed. It is based on nonstaggered cell-centered finite volume method. Second-order upwind scheme with least-square reconstruction is used for convective fluxes. The SIMPLE method is implemented to couple the pressure and velocity fields. And the time derivatives in the momentum equations are discretised using a second-order Euler backward-differencing scheme. The discretised linear equations are solved by the preconditioned Biconjugate Gradient Stabilized method(Bi-CGSTAB). The developed solver is applied to validation problems using hybrid meshes.

• 한국항공우주학회지
• /
• 제30권4호
• /
• pp.44-52
• /
• 2002

비정상, 비압축성 Navier-Stokes 코드를 이용하여, 저 레이놀즈수 유동에서 flapping 운동을 하는 익형의 공력특성을 수치해석적인 방법으로 연구하였다. 비정상 유동장의 효율적인 계산을 위하여, 개발된 코드는 MPI 프로그래밍 기법을 이용하여 병렬처리 되었으며, 난류 유동장의 계산을 위해 2방정식 난류모델의 하나인 k-$\omega$ SST 모델을 적용하였다. 익형의 3가지 운동모드 즉, pitching, plunging, flapping과 주파수 및 진폭의 변화 그리고 두께와 캠버의 변화에 의한 공력특성을 살펴보았고, 이를 위해 NACA4자 계열의 익형을 이용하였다. 해석 결과는 실험치와 비교하여 보았을 때 잘 일치하였으며, 각 운동모드에서의 공기역학적 특성을 파악할 수 있었다.

• 지질공학
• /
• 제3권2호
• /
• pp.149-165
• /
• 1993

지하매질의 유동특성을 알기 위한 직접적인 방법은 아직도 현장수리시험에 의존하고 있다. 열극암반에서 유동모델의 개념정립을 위하여는 열극체계분포와 특히, 수리시험해석결과가 중요시되고 있다. 국내에서 수행되고 있는 수리시험은 조사공구간별에 따른 수리전도도를 방사상의 정상류로 가정하여 결과를 해석하고 있다. 또한 시험대상 매질을 균질한 대수층 또는 단순한 기하학적 형태의 매질로 간주하며 구간별 투수량계수는 시험구간을 가로지르는 각 열극의 투수량계수의 합과 같다고 가정한다. 국내 기존의 수리시험 및 해석방법으로는 수리학적 경계면의 영향(boundary effects)이나 유동로의 기하학적 형태(flow geometry)에 대한 자료를 얻기가 힘들며, 일정한 조사공 주변을 벗어나면 투수성열극의 연결성에 대한 정보를 구할 수 없다. 열극특성을 고려한 지하수유동 해석을 위하여 단일공 정압주입시험을 실시하였으며, 비정상류해석방법을 통하여 유동차원(flow dimension)에 대한 분석을 시도하였다. 시추시에는 단일패커시험을 일정구간별로 시행하였으며, 조사공의 시추후에는 이중패커에 의한 구간별 시험을 실시하였다. 비정상류해석으로 구한 수리전도도값은 정상류해석의 결과와 큰 차이(10배 이내)는 없었으나, 정압하에서 도출된 유동량변화곡선에서 유동차원분석이 가능하였다. 상부구간(<10m 깊이)의 단일 및 이중패커시험결과는 모두 정상류의 유동차원이 나타났으며, 이는 영향반경의 경계면이 open system임을 알 수 있다. 15m 깊이에서 도출된 유동량변화곡선은 1차원 유동상태에서 3차원(구상유동)으로 변화하였다. 하부구간(25m 깊이)의 시험결과는 closed system 특성이 관찰되었으며, 이는 조사구간에서 연결된 열극이 수리적으로 격리되어 있음을 알 수 있다. 현장수리시험으로부터 보다 많은 자료를 도출하기 위하여는 무엇보다도 수리시험장비의 보완이 필요하다. 특히 조사구간을 완벽하게 분리할 수 있는 패커장비와 미세한 유동량변화를 계측할 수 있는 유량계의 확보가 필수적이며, 조사구간의 압력변화를 자동기록할 수 있는 계측기기가 필요하였다.

• 대한조선학회논문집
• /
• 제33권4호
• /
• pp.48-59
• /
• 1996

비정상 점성유동 수치해석단계는 연속방정식을 만족시키는 공간해석단계와 시간전진단계로 구분할 수 있다. 본 연구에서는 공간해석단계의 압력 Poisson 방정식의 해를 구하는데 스팩트럴법을 이용하였다. 압력 Poisson 방정식의 최고차미분항을 Fourier 급수로 전개하면 압력 및 압력의 1차미분항의 Fourier 급수의 적분으로 표현되므로 Gibb's 현상을 제거할 수 있어, 비주기성인 경우에도 스팩트럴법을 적용할 수 있다. 수치해법의 검증을 위하여 2차원 원주상체 및 날개주위 비정상 점성유동을 수치해석하였고, 그 결과를 비교하여 보았다.

• 한국항공우주학회지
• /
• 제33권12호
• /
• pp.1-8
• /
• 2005

타원형 실린더의 두께와 받음각 및 레이놀즈수가 실린더에 작용하는 항력과 양력에 어떤 영향을 미치는가를 고찰하기 위해서 수치적 연구를 수행하였다. 타원형 실린더 주위를 흐르는 비정상 점성 유동을 해석하기 위하여 SIMPLER 기법을 이용한 2차원 비압축성 Navier-Stokes 유동 해석 프로그램을 개발하였으며, 두께-시위길이 비가 0.2, 0.4, 0.6인 타원형 실린더 형상에 대해서 레이놀즈수가 400, 600인 조건, 그리고 받음각이 10도, 20도, 30도인 조건하에서 유동을 해석하였다. 본 연구를 통해서 실린더 두께 비와 받음각 및 레이놀즈수가 항력과 양력 계수의 시간 평균값 및 진폭의 크기, 그리고 진동 주기에 크게 영향을 미침을 확인할 수 있었다.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2011년도 정기 학술발표대회
• /
• pp.75-75
• /
• 2011

본 논문에서는 진동중인 교량 거더에 작용하는 풍하중을 산정하고 그에 따른 플러터 발생풍속을 예측하기 위하여 분산형 전산환경을 활용한 수치해석 연구를 수행하였다. 분산형 전산환경은 웹 포탈을 기반으로 수치해석 환경을 제공하는 수치풍동 시스템으로서, 전산유체역학(CFD : Computational Fluid Dynamics)에 대한 전문지식이 부족한 사용자들도 격자생성, 수치해석자를 이용한 계산, 가시화 등의 전 과정을 편리하게 수행할 수 있는 차세대 토목분야 연구 환경이다. 본 시스템은 그리드스피어(GfidSphere)를 기반으로 구성되었으며, 기본적으로 사용자 관리, 세션 관리, 그룹 관리, 레이아웃 관리 등을 제공하여 사용자가 포탈을 통해서 다양한 서비스를 쉽게 사용할 수 있는 환경을 구축하도록 도와준다. 수치해석을 위한 유체 지배방정식은 2차원 비정상 비압축성 RANS(Reynolds-Averaged Navier-Stokes) 방정식이며, pseudo compressibility 방법을 적용하였다. 비정상 유동장을 해석하기 위하여 이중시간 전진법(dual time stepping)을 사용하였으며, 수렴가속화를 위해 Multi-grid 기법을 적용하였다. 또한 난류 유동장 해석을 위해서 $k-{\omega}$ SST 난류 모델을 사용하였으며, 난류 천이 과정에서의 유동을 모사하기 위하여 Total stress limitation 방법을 적용하였다. 교량 거더의 연직과 회전방향의 2자유도 움직임을 모사하기 위하여 동적격자 기법을 도입하였다. 교량 거더 주변의 비정상 유동해석 결과를 통해, 거더 표면에서 떨어져나가는 크고 작은 와류의 영향으로 양력 및 모멘트 계수 그래프가 중첩된 진폭과 주기를 갖고 주기적으로 나타나는 것을 확인할 수 있었다. 또한 계산된 비정상 공기력을 적용한 2자유도 플러터 방정식을 통하여 플러터 발생풍속을 산정하였다. 최종적으로 본 연구에서 계산된 결과의 타당성을 검증하기 위하여 수치적으로 구한 플러터 발생풍속과 기존의 실험 및 수치해석 결과를 비교하였으며, 결과는 잘 일치하였다.

• 한국항공우주학회지
• /
• 제40권5호
• /
• pp.413-422
• /
• 2012

덕티드팬을 추진 장치로 사용하는 소형 무인항공기는 도심 및 협소한 공간에서 정찰 및 감시에 사용 가능하며, 프로펠러에 비해 높은 추진 효율과 추력 특성을 나타낸다. 덕티드팬 무인항공기의 운용 거리와 비행 시간을 증가시키기 위해서는 정지 비행 및 전진 비행시의 추력 특성연구가 중요하며 비행 안정성 확보를 위해서는 비정상 3차원 유동 특성 연구가 필수적이다. 본 연구에서는 동익과 정익으로 구성된 덕티드팬의 설계 결과 검증과 안정적인 비행 특성을 확인하기 위해 덕티드팬의 추력 특성과 비정상 3차원 유동장을 계측하였다. 덕티드팬의 정지 및 전진 비행시의 추력 특성은 소형 아음속 풍동의 6분력 밸런스 시스템을 이용하여 측정되었고, 비정상 3차원 유동장은 $45^{\circ}$ 경사열선의 프로브 고정법에 의해 분석되었다. 덕티드팬의 덕트와 정익이 추력특성에 다소 큰 영향을 미치며, 정익에 의해 덕티드팬의 안정적인 비행이 가능함을 확인하였다.

• 한국결정성장학회지
• /
• 제4권1호
• /
• pp.46-56
• /
• 1994

균일한 자장이 Czochralski도가니에 가하여졌을 때의 유동장, 온도장 및 boron의 농도장이 수치적으로 연구되었다.이러한 Czochralski 도가니에 가하여졌을 때의 유동장, 온도장 및 boron의 농도 장이 수치적으로 연구되어 . 이러한 Caochralski 유동장을 지배하는 인자로는 온도차이에 의한 부력, 자유표면에서의 반경방향으로의 온도 구배로 인한 thermocapillarity, 도가니와 결정의 회전으로 인한 원심력, 외부에서 걸어준 자장에 의한 자성력, boron의 확산계수 및 segregation coefficient 등이 있다. 여기에서 boron의 농도가 매우 낮으므로 농도차이에 의한 부력의 효가 무시되므로 boron의 농도장은 유동장 및 온도장에 영향을 미치지 못한다. 희전방향으로의 대칭성으로 부터 먼저 정상 상태에 대한 자오평면(mericional plane)에서의 속도성분과 회전방향의 속도 성분이 구하여졌으며 온도장도 해석되었다. 이러한 정상상태에서의 유동장 및 온도장으로부터 boron의 농도에 관한 비정상상태에서의 농도구배가 Marangoni convection을 야기하고 있다. 또 비정상 상태의 boron의 농도장의 전달현상은 주로 대류작용에 의존하고 있다.존하고 있다.

• 대한기계학회논문집B
• /
• 제25권11호
• /
• pp.1518-1526
• /
• 2001

In this paper, the mechanism of unsteady potential interaction and wake interaction in one stage axial turbine is numerically investigated at design point in two-dimensional viewpoint. The numerical technique used is the upwind scheme of Van-Leer's Flux Vector Splitting (FVS) and Cubic spline interpolation is applied on zonal interface between stator and rotor. The inviscid analysis is used to embody the influence of potential interaction only and viscous analysis is used to embody the influences of both potential interaction and wake interaction at the same time. The potential-flow disturbance from the stator into a rotor passage and the periodic blockage effect of rotor produce the unsteady pressure on the blade surface in inviscid analysis. After the wake is cut by rotor, two counterrotating votical patterns flanking the wake centerline in the passage are generated. So, these phenomena magnify the unsteady pressure in viscous analysis than that in inviscid analysis. The resulting unsteady forces on the rotor, generated by the combined interaction of the two effects by potential and wake interaction, are discussed.