• 대한기계학회논문집
• /
• 제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.

• International Journal of Aeronautical and Space Sciences
• /
• 제1권1호
• /
• pp.1-12
• /
• 2000

The wake geometries of a two-bladed rotor in axial flights using a time-marching free-wake method without a non-physical model of the far wake are calculated. The computed free-wake geometries of AH-1G model rotor in climb flight are compared with the experimental visualization results. The time-marching free-wake method can predict the behavior of the tip vortex and the wake roil-up phenomena with remarkable agreements. Tip vortices shed from the two-bladed rotor can interact with each other significantly. The interaction consists of a turn of the tip vortex from one blade rolling around the tip vortex from the other. Wake expansion of wake geometries in radial direction after the contraction is a result of adjacent tip vortices begging to pair together and spiral about each other. Detailed numerical results show regular pairing phenomenon in the climb flights, the hover at high angle of attack and slow descent flight too. On the contrary, unstable motions of wake are observed numerically in the hover at low angle of attack and fast descent flight. It is because of the inherent wake instability and blade-vortex-interaction rather then the effect of recirculation due to the experimental equipment.

• 한국군사과학기술학회지
• /
• 제9권2호
• /
• pp.152-160
• /
• 2006

Time marching methods are well-suited for high speed compressible flow computations. However, it is well known that the time marching methods suffer a slow down in convergence due to disparity in Eigenvalues. A local preconditioning method is one of numerical methods to enhance convergence characteristics of low mach number flows by modifying Eigenvalues of the governing equations. In this paper, the local preconditioning method of Weiss is applied to a 2 dimensional Navier-Stokes code and the efficiency of the preconditioning method is shown through a number of computational examples.

• International Journal of Aeronautical and Space Sciences
• /
• 제5권2호
• /
• pp.62-70
• /
• 2004

CFD method has been coupled with a time-marching free-wake model by usingfield velocity approach suggested by J. D. Baeder (Ref. 1). The coupled method hasbeen applied to rectangular and BERP-like blades and the calculated perfonnance dataare compared with the experimental results.For hovering analysis, the present method could yield sufficiently good resultswith reasonable computation time and is particularly suitable for the flow fieldanalysis with the complex shaped blade.

• Structural Engineering and Mechanics
• /
• 제73권1호
• /
• pp.65-73
• /
• 2020

The present work describes a new time marching procedure for structural dynamics analyses. In this novel technique, time integration parameters are automatically evaluated according to the properties of the model. Such parameters are locally defined, allowing the user to input a numerical dissipation property for each element, which defines the amount of numerical dissipation to be introduced. Since the integration parameters are locally defined as a function of the structural element itself, the time marching technique adapts according to the model, providing enhanced accuracy. The new methodology is based on displacement-velocity relations and no computation of accelerations is required. Furthermore, the method is second order accurate, it has guaranteed stability, it is truly self-starting and it allows highly controllable algorithm dissipation in the higher modes. Numerical results are presented and compared to those provided by the Newmark and the Bathe methods, illustrating the good performance of the new time marching procedure.

• 한국전산유체공학회지
• /
• 제20권1호
• /
• pp.32-38
• /
• 2015

A 3-dimensional compressible turbulent boundary layer solver has been developed. A time marching method is used to integrate the turbulent boundary layer equations. While the direct integration of the boundary layer equations is performed for unseparated flow regions, the inverse integration is performed for separated flow regions. The program is verified for flows that have analytical solutions or other numerical results. The solver will be merged with an Euler solver for viscous-inviscid interaction.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2011년도 추계학술대회 논문집
• /
• pp.171-177
• /
• 2011

In this paper, a square simply supported panel flutter have been considered at high supersonic flow by using coupled fluid-structure (FSI) analysis that based on time domain method. The Reynolds-Average Navier Stokes (RANS) equation with Spalart-Allmaras turbulent model were applied for unsteady flow problems of panel flutter. A fully implicit time marching schemed based on the Newmark direct integration method is used for calculating the coupled aeroelastic governing equations of it. In addition, the SOL 145 solver of MSC.NASTRAN was used to investigate flutter velocity based on PK-method of Piston theory. Our numerical results indicated that there is a good agreement result between Piston Theory in MSC.NASTRAN and coupled fluid-structure analysis.

• 한국항공우주학회지
• /
• 제37권4호
• /
• pp.329-335
• /
• 2009

비정상 패널법과 시간 전진 자유 후류법을 연계하여 제자리 및 수직 비행 하는 헬리콥터 로터의 공력 및 익단와류를 검증하였다. 검증된 수치 방법을 이용하여 저속 하강 비행하는 로터에서 발생하는 BVI(Blade Vortex Interaction)의 비정상 공력을 해석 하고 이때 나타나는 익단 와류에 의한 공력의 변화를 고찰하였다.

• 한국전산유체공학회지
• /
• 제1권1호
• /
• pp.26-34
• /
• 1996

The performance of the SIMPLE, SIMPLE-C and PISO algorithms for the treatment of the pressure-velocity coupling in fluid flow problems were examined by comparing the computational effort required to obtain the same level of the convergence. Example problems are circular duct and 90-degree bent square-duct. For circular duct case, laminar and turbulent flow were computed. For 90-degree bent square-duct case, laminar flow was simulated by the time-marching method as well as the iterative method. The convergence speed of the other two algorithms are not always superior to SIMPLE algorithm. SIMPLE algorithm is faster than SIMPLE-C algorithm in the simple laminar flow calculations. The application of the PISO algorithm in three dimensional general coordinates is not so effective as in two-dimensional ones. Since computational time of PISO algorithm is increased at each time step(or iterative step) in three dimension, the total convergence speed is not decreased. But PISO algorithm is stable for large time step by using time marching method,.

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

Two representative finite volume methods, i.e., the time marching method and the pressure correction method, were applied to 8 centrifugal compressor impeller flows, with low to very high level of pressure ratio, among which 7 impellers' experimental performance is given in the open literature. The present study is focused on the prediction differences from both methods, developed by the authors, in the pressure correction method's point of view. In all cases, the time marching method gives a satifactory solution, but the pressure correction method does not. Up to about $18\%$ less level of total-to-total pressure ratio is predicted by the pressure correction method as the level of the impeller pressure ratio increases up to about 10. The drop of total pressure ratio is caused by the underestimation of static pressure rise which seems to be attributed to inappropriate linearization and discretization of the pressure/density coupling terms in the pressure correction equation.