• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.666-670
• /
• 2008

The low-speed preconditioning technique is applied to solve the compressible Reynolds averaged Navier-Stokes equations for low-speed flows. The space discretization is based on Roe's flux-difference splitting with third-order-accurate MUSCL extrapolation. Time integration is performed employing a diagonal approximate factorization algorithm. The dual-time stepping has been incorporated to solve the unsteady flows. Full multigrid method is implemented to accelerate the convergence rate. To verify the algorithms several cases have been tested. Demonstrated the improvement on convergence and quality of the solution.

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

A dual-time stepping algorithm combined with a parallelized multigrid DADI method is presented to predict the dynamic damping coefficients. The Basic Finner model is chosen to validate the prediction capability of the present unsteady Euler method. The linearity of the pitch- and roll-damping coefficients is shown in the low angular rates and the interesting large drop and stiff increment in transonic region for roll-damping coefficients are explained in detail. Through the analysis for the pressure distributions at Mach number 1.0 to 1.2, the sudden drop results from the normal shock and the stiff increment of roll-damping reflects the transition of the normal shock to the oblique shock. The results also show that the Euler equations can give the damping coefficients with a comparable accuracy.

• 한국항공우주학회지
• /
• 제38권12호
• /
• pp.1153-1161
• /
• 2010

본 연구에서는 주기적 비정상 유동 해석을 위해 푸리에 변환을 이용하는 조화 균형법의 효율적인 해법을 제안한다. 내재적으로 유속항을 처리하고 외재적으로 조화 원천항을 처리하였다. 외재적 조화 균형법 보다 더 빠르게 수렴 시킬 수 있으며 내재적 조화 균형법을 적용할 때 추가되는 자코비안 행렬을 처리할 필요가 없다. 또한 완전 내재적 기법에 상응하는 수준의 수렴안정성을 확인할 수 있었다. 2차원 비정상 유동 문제로 피칭하는 NACA0012 익형에 적용하였으며 이중 시간 적분법 및 외재적 Runge-Kutta기법의 해와 매우 일치하는 결과를 얻었다.

• 한국전산유체공학회지
• /
• 제17권1호
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.389-396
• /
• 2010

In this paper, a new computational code was developed using Chorin's artificial compressibility method to solve the two-dimensional incompressible Navier-Stokes equations. In spatial derivatives, Roe's flux difference splitting was used for the inviscid flux, while central differencing was used for the viscous flux. Furthermore, AF-ADI with dual time stepping method was implemented for accurate unsteady computations. Two-equation turbulence models, Menter's $k-{\omega}$ SST model and Coakley's $q-{\omega}$ model, hae been adopted to solve high-Reynolds number flows. A number of numerical simulations were carried out for steady laminar and turbulent flow problems as well as unsteady flow problem. The code was verified and validated by comparing the results with other computational results and experimental results. The results of numerical simulations showed that the present developed code with the artificial compressibility method can be applied to slve steady and unsteady incompressible flows.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2003년도 제21회 추계학술대회 논문집
• /
• pp.157-163
• /
• 2003

A systematic numerical experiment has been conducted to study droplet gasification in high pressure environments with pressure oscillations. The general frame of previous rigorous model[1] is retained but tailored for flash equilibrium calculation of vapor-liquid interfacial thermodynamics. Time-dependent conservation equations of mass, momentum, energy, and species concentrations are formulated in axisymmetric coordinate system for both the droplet interior and ambient gases. In addition, a unified property evaluation scheme based on the fundamental equation of state and empirical methods are used to find fluid thermophysical properties over the entire thermodynamic domain of interest. The governing equations with appropriate physical boundary conditions are numerically time integrated using an implicit finite-difference method with a dual time-stepping technique. A series of calculation have been carried out to investigate the gasification of an isolated n-pentane droplet in a nitrogen gas environment over a wide range of ambient pressures and frequencies. Results show that the mean pressures and frequencies of the ambient gas have strong influences on the characteristics of the droplet gasification. The amplitude of the response increases with increasing pressure, and the magnitude of the vaporization response increases with the frequency.

• 한국항공우주학회지
• /
• 제40권11호
• /
• pp.948-956
• /
• 2012

본 논문에서는 SDM 형상의 세로와 방향의 안정성 미계수를 예측하였다. 피치와 요 방향에 대한 강제조화 진동운동을 이용하여 정적 및 동적 미계수를 한 번에 계산하였다. 계산은 비정상 해석을 위한 이중시간 적분법을 적용한 3차원 Euler 해석자를 사용하여 수행하였다. 본 연구에서는 마하수뿐만 아니라 다양한 운동 변수에 따른 미계수를 예측하였다. 예측된 결과는 이전에 발표된 수치적, 실험적 연구 결과들과 비교하여 검증하였다.

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

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• 한국항공우주학회지
• /
• 제34권1호
• /
• pp.24-31
• /
• 2006

본 연구에서는 완전 내재적 기법을 이용한 2계 자유도 모델에 대한 플러터 해석을 수행하였다. 유동해석을 위하여 2차원 Navier-Stokes 지배방정식을 ε-SST 난류모델과 DP-SGS 병렬화 기법을 이용해 구성하였다. 구조해석을 위하여 피치 와 플런지의 2계 자유도를 갖는 모델을 구성하였으며 시간영역에서의 해석을 위하여 유동해석과 마찬가지로 이중 시간 전진 기법을 이용하였다. 가상 시간 전진에서 유체-구조 결합을 통해서 강성결합을 구현하였으며 이를 천음속 플러터 예측에 적용하였다. 플러터 해석의 전형적인 응답인 감쇠응답, 중립응답, 발산응답 및 limit cycle oscillation (LCO) 등을 계산하였으며, 더불어 플러터 속도 경계곡선을 작성하였다. 천음속 플러터 해석은 난류모델 뿐만 아니라 유체-구조 결합 방법에 따라 다른 특성을 보임을 확인하였다.

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

본 논문에서는 전산유체역학을 이용하여 3차원 비행체 형상에 대한 동안정 미계수를 예측하였다. 피치에 대한 미계수는 피치방향의 조화진동운동을 통하여 계산하였고 롤 감쇠계수는 비관성 좌표계에 대한 정상해석을 통하여 계산하였다. 계산은 Basic Finner와 SDM 형상에 대해 수행했으며 다른 연구자의 실험적/수치적 결과와 비교하였다. 유동 계산을 위해 비관성 좌표계와 관성 좌표계에서 모두 사용할 수 있는 3차원 Euler 해석자를 개발하였다. 시간 정확성을 유지한 비정상 해석을 위해 이중시간적분법을 적용하였다. 동안정 미계수계산 결과는 다른 수치 및 실험적 연구 결과들과 잘 일치하는 것을 알 수 있었다.