• 한국전산유체공학회지
• /
• 제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.

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

Time marching algorithms applied to compressible Navier-Stokes equation have a convergence problem at low Mach number. It is mainly due to the eigenvalue stiffness and pressure singularity as Mach number approaches to zero. Among the several methods to overcome the shortcomings of time marching scheme, time derivative preconditioning method have been used successfully. In this numerical analysis, we adopted a preconditioner of K.H. Chen and developed a two-dimensional, axisymmetric Navier-Stokes program. The steady state driven cavity flow and backward facing step flow problems were computed to confirm the accuracy and the robustness of preconditioned algorithm for low Mach number flows. And the transonic and supersonic flows insice the JPL axisymmetric nozzle internal flow is exampled to investigate the effects of preconditioning at high Mach number flow regime. Test results showed excellent agreement with the experimental data.

• 대한기계학회논문집
• /
• 제18권7호
• /
• pp.1840-1850
• /
• 1994

Enhancement of numerical algorithms for low speed compressible flow will be considered. Contemporary time-marching algorithm has been widely accepted and applied as the method of choice for transonic, supersonic and hypersonic flows. In the low Mach number regime, time-marching algorithms do not fare as well. When the velocity is small, eigenvalues of the system of compressible equations differ widely so that the system becomes very stiff and the convergence becomes very slow. This characteristic can lead to difficulties in computations of many practical engineering problems. In the present approach, the time-derivative preconditioning method will be used to control the eigenvalue stiffness and to extend computational capabilities over a wide range of flow conditions (from very low Mach number to supersonic flow). Computational capabilities of the above algorithm will be demonstrated through computation of a variety of practical engineering problems.

• 대한기계학회논문집
• /
• 제19권9호
• /
• pp.2283-2296
• /
• 1995

Experiments of normal shock wave/turbulent boundary layer interaction were conducted in a supersonic diffuser. The flow Mach number just upstream of the normal shock wave was in the range of 1.10 to 1.70 and Reynolds number based upon the turbulent boundary layer thickness was varied in the range of 2.2*10$^{[-994]}$ -4.4*10$^{[-994]}$ . The wall pressures in streamwise and spanwise directions were measured for two test cases, in which the turbulent boundary layer thickness incoming into the supersonic diffuser was changed. The results show that the interactions of normal shock wave with turbulent boundary layer in the supersonic diffuser can be divided into three patterns, i.e., transonic interaction, weak interaction and strong interaction, depending on Mach number. The weak interactions generate the post-shock expansion which its strength is strong as the Mach number increases and the strong interactions form the pseudo-shock waves. From the spanwise measurements of wall pressure, it is known that if the flow Mach number is low, the interacting flow fields essentially appear two-dimensional, but they have an apparent 3-dimensionality for the higher Mach numbers.

• 대한기계학회논문집B
• /
• 제31권9호
• /
• pp.814-821
• /
• 2007

A hybrid method is presented for efficient computation of turbulent flow noise at low Mach numbers. In this method, the turbulent flow field is computed by incompressible large eddy simulation (LES), while the acoustic field is computed with the linearized perturbed compressible equations (LPCE) derived in this study. Since LPCE is computed on the rather coarse acoustic grid with the flow variables and source term obtained by the incompressible LES, the computational efficiency of calculation is greatly enhanced. Furthermore, LPCE suppress the instability of perturbed vortical mode and therefore secure consistent and stable acoustic solutions. The proposed LES/LPCE hybrid method is applied to three low Mach number turbulent flow noise problems: i) circular cylinder, ii) isolated flat plate, and iii) interaction between cylinder wake and airfoil. The computed results are closely compared with the experimental measurements.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.251-255
• /
• 2008

Flow field around helicopter involves incompressible flow near the blade root and compressible flow at the blade tip. A problem occurs for low Mach number flow due to the stiffness of the governing equations. Time-derivative preconditioning techniques have been incorporated to reduce the stiffness that occurs at low speed region. The preconditioned form of the compressible Navier-Stokes and Euler equations is used. Computations are performed for the Caradonna-Tung's hovering and non-lifting forward flight case. Computational results are in good agreement with the experimental data.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년 추계학술대회논문집
• /
• pp.251-255
• /
• 2008

Flow field around helicopter involves incompressible flow near the blade root and compressible flow at the blade tip. A problem occurs for low Mach number flow due to the stiffness of the governing equations. Time-derivative preconditioning techniques have been incorporated to reduce the stiffness that occurs at low speed region. The preconditioned form of the compressible Navier-Stokes and Euler equations is used. Computations are performed for the Caradonna-Tung's hovering and non-lifting forward flight case. Computational results are in good agreement with the experimental data.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 추계학술대회
• /
• pp.804-809
• /
• 2003

In this study, 2D computations of the Aeolian tones for some obstacles (circular cylinder, square cylinder and NACA0012 airfoil) are simulated. First of all, we calculate the flow noise generated by a uniform flow around a two-dimensional circular cylinder at Re=150 are simulated by applying the finite difference lattice Boltzmann method (FDLBM). The third-order-accurate up-wind scheme (UTOPIA) is used for the spatial derivatives, and the second-order-accurate Runge-Kutta scheme is applied for the time marching. The results show that we successively capture very small acoustic pressure fluctuation with the same frequency of the Karman vortex street compared with the pressure fluctuation around a circular cylinder. The propagation velocity of the acoustic waves shows that the points of peak pressure are biased upstream due to the Doppler effect in the uniform flow. For the downstream, on the other hand, it is faster. To investigate the effect of the lattice dependence, furthermore, simulations of the Aeolian tones at the low Reynolds number radiated by a square cylinder and a NACA0012 airfoil with a blunt trailing edge at high incidence are also investigated.

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

The turbulence noise generated from blunt trailing-edge is numerically predicted by using the hydrodynamic/acoustic splitting method at the Reynolds number based on thickness of flat plate, $Re_h=1000$, and the freestream Mach number $M_o=0.2$. The turbulent flow field is simulated by incompressible large-eddy simulation and the acoustic field is predicted efficiently with the linearized perturbed compressible equations (LPCE) recently proposed by the authors. The turbulent flow characteristics are validated with the results of the previous experimental study and direct numerical simulation. The acoustic properties predicted from LPCE are compared with the solutions of analytical formulations.

• 한국항공우주학회지
• /
• 제37권1호
• /
• pp.17-27
• /
• 2009

압축성 영역에서 개발된 2상 유동 RoeM과 AUSMPW+를 전 마하수 영역 해석을 위해 예조건화 하였다. 압축성 영역에서 개발된 풍상차분 수치기법들은 수치점성에서 음속을 포함하는 항들을 갖고 있는데, 이는 비압축성 영역에서 관련 수치점성 크기를 지나치게 크거나 작게 만들어 해의 정확도를 떨어뜨리거나 수치 불안정 현상을 일으킨다. 따라서 예조건화된 시스템 고유값을 수치점성에 고려하여 전 마하수 영역에서 해의 정확성을 확보하였고, 점근 분석을 통해 개발된 수치기법이 비압축 극한 영역에서 보존형 지배방정식과 일관됨을 확인하였다. 수치 예제들을 통해 예조건화된 2상 유동 RoeM과 AUSMPW+가전 마하수 영역에서 합리적인 해를 도출함을 알 수 있었다.