• Journal of Mechanical Science and Technology
• /
• 제15권5호
• /
• pp.671-680
• /
• 2001

Boundary condition is one of the major factors to influence the numerical stability and solution accuracy in numerical analysis. One of the most important physical boundary conditions in the flowfield analysis is the wall boundary condition imposed on the body surface. To solve a two-dimensional Euler equation, totally four numerical wall boundary conditions should be prescribed. Two of them are supplied by the flow tangency condition. The other two conditions, therefore, should be prepared additionally in a suitable way. In this paper, four different sets of wall boundary conditions are proposed and then applied to solve high-speed flowfields around a quarter circle geometry. A two-dimensional compressible Euler solver is prepared based on the finite volume method. This solver hires three different upwind schemes; Steger-Warmings flux vector splitting, Roes flux difference splitting, and Lious advection upstream splitting method. It is found that the way to specify the additional numerical wall boundary conditions strongly affects the overall stability and accuracy of the upwind schemes in high-speed flow calculation. The optimal wall boundary conditions should be also chosen very carefully depending on the numerical schemes used to solve the problem.

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

높은 효율의 풍력터빈 블레이드을 위해 와류 셀이 장착된 에어포일의 정지상태 및 동실 속 상태에서의 유동제어 특성을 수치적으로 연구하였다. 수치기법은 Roe의 flux-difference-splitting을 사용한 격자점 중심 유한체적법과 이중시간 전진 기법을 사용하는 내재적 시간적분법을 사용하였다. 계산결과 와류 셀을 장착한 경우 셀 내부의 부압으로 인해 양항비증가를 얻을 수 있음을 확인하였다. 동실속의 경우 셀 내부의 와류에 의해 hysterisis 현상을 상당히 감소시킬 수 있음을 확인하였다.

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

An Euler solver is developed to predict accurate aerodynamic data such as lift coefficient, drag coefficient, and moment coefficient. The conservation law form of the compressible Euler equations are used in the generalized curvilinear coordinates system. The Euler solver uses a finite volume method and the second order Roe's flux difference splitting scheme with min-mod flux limiter to calculate the fluxes accurately. An implicit scheme which includes the boundary conditions is implemented to accelerate the convergence rate. The multi-block grid is integrated into the flow solver for complex geometry. The flowfields are analyzed around NACA 0012 airfoil in the cases of $M_{\infty}=0.75,\;\alpha=2.0\;and\;M_{\infty}=0.80,\;\alpha=1.25$. The numerical results are compared with other numerical results from the literature. The final goal of this research is to prepare a robust and an efficient Navier-Stokes solver eventually.

• 한국전산유체공학회지
• /
• 제15권2호
• /
• pp.55-61
• /
• 2010

A dependence of weighting parameter in preconditioning method for solving low Mach number flow with incompressible flow nature is investigated. The present preconditioning method employs a finite-difference method applied Roe‘s flux difference splitting approximation with the MUSCL-TVD scheme and 4th-order Runge-Kutta method in curvilinear coordinates. From the computational results of benchmark flows through a 2-D backward-facing step duct it is confirmed that there exists a suitable value of the weighting parameter for accurate and stable computation. A useful method to determine the weighting parameter is introduced. With this method, high accuracy and stable computational results were obtained for the flow with low Mach number in the range of Mach number less than 0.3.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.23-26
• /
• 2002

An evaluation of one algebraic and two one-equation eddy viscosity-transport turbulence closure models as implemented to the CFDS(Characteristic Flux Difference Splitting) scheme is presented for the efficient computation of the turbulent flow. Comparisons of Baldwin-Lomax model as algebraic turbulence model and Baldwin-Barth and Spalart-Allmaras model as one-equation turbulence model are presented for three test cases for 3-dimensional flow. The numerical result of the CFDS schem is examined through comparison with the experimental data.

• Journal of Mechanical Science and Technology
• /
• 제18권12호
• /
• pp.2190-2203
• /
• 2004

A combined procedure for two-dimensional Delaunay mesh generation algorithm and an adaptive remeshing technique with higher-order compressible flow solver is presented. A pseudo-code procedure is described for the adaptive remeshing technique. The flux-difference splitting scheme with a modified multidimensional dissipation for high-speed compressible flow analysis on unstructured meshes is proposed. The scheme eliminates nonphysical flow solutions such as the spurious bump of the carbuncle phenomenon observed from the bow shock of the flow over a blunt body and the oscillation in the odd-even grid perturbation in a straight duct for the Quirk's odd-even decoupling test. The proposed scheme is further extended to achieve higher-order spatial and temporal solution accuracy. The performance of the combined procedure is evaluated on unstructured triangular meshes by solving several steady-state and transient high-speed compressible flow problems.

• 한국전산유체공학회지
• /
• 제3권2호
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.

• International Journal of Aeronautical and Space Sciences
• /
• 제8권2호
• /
• pp.44-53
• /
• 2007

A three-dimensional viscous flow solver has been developed for the prediction of the aerodynamic performance of hovering helicopter rotor blades using unstructured hybrid meshes. The flow solver utilized a vertex-centered finite-volume scheme that is based on the Roe's flux-difference splitting with an implicit Jacobi/Gauss-Seidel time integration. The eddy viscosity are estimated by the Spalart- Allmaras one-equation turbulence model. Calculations were performed at three operating conditions with varying tip Mach number and collective pitch setting for the Caradonna-Tung rotor in hover. Additional computations are made for the UH-60A rotor in hover. Reasonable agreements were obtained between the present results and the experiment in both blade loading and overall rotor performance. It was demonstrated that the present vertex-centered flow solver is an efficient and accurate tool for the assessment of rotor performance in hover.

• 대한기계학회논문집
• /
• 제19권11호
• /
• pp.3031-3040
• /
• 1995

Two-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with the unstructured triangular meshes. The pressure-velocity coupling is handled by the artificial compressibility algorithm due to its computational efficiency associated with the hyperbolic nature of the resulting equations. The convective fluxes are obtained by the Roe's flux difference splitting scheme using edge-based connectivities and higher-order differences are achieved by a reconstruction procedure. The time integration is based on an explicit four-stage Runge-Kutta scheme. Numerical procedures with local time stepping and implicit residual smoothing have been implemented to accelerate the convergence for the steady-state solutions. Comparisons with experimental data and other numerical results have proven accuracy and efficiency of the present unstructured approach.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 1999년도 유체기계 연구개발 발표회 논문집
• /
• pp.179-186
• /
• 1999

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67), and to evaluate the performances of k-$\epsilon$ and Baldwin-Lomax turbulence models. A finite volume method is used for spatial discretization. And, the equations are solved implicitly in time with the use of approximate factorization. Upwind difference scheme is used for inviscid terms, but viscous terms are centrally differenced. The flux-difference-splitting of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. And, the results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, we conclude that k-$\epsilon$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost same.