• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2329-2338
• /
• 1993

An accurate analysis procedure to solve the flow about a flat plate at various incidences has been developed. The Navier-Stokes equations of stream function and vorticity form are solved in a sufficiently large computational domain, in which the grid lines are mutually orthogonal. The details of the flow near the singularity at the tip of the plate is well captured by the analytic solution which is asymptotically matched to the numerically generated outer solution. The solution for each region is obtained iteratively : the solution of one (inner or outer) region uses that of the other as the boundary condition after each cycle. The resulting procedure is accurate everywhere and also computationally efficient as the singularity has been removed. It is applied to the flat plate for a wide range of Re : the results agree very well with the existing computation and experiment.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.91-97
• /
• 2001

The thin-layer Navier-Stokes equations are solved for the hypersonic flow over blunt cone configurations with applications to laminar as well as turbulent flows. The equations are expressed in the forms of flux-vector splitting and explicit algorithm. The upwind schemes of Steger-Warming and van Leer are investigated in their ability to accurately predict the heating loads along the surface of the body. A comparison with the second order extensions of these schemes is made and a hybrid scheme incorporating a combination of central differencing and flux-vector-splitting is presented. This scheme is also investigated in its ability to accurately predict heat transfer distributions.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.91-98
• /
• 2000

The three-dimensional unsteady compressible Full Navier-Stokes equation solver with sliding multi-block method has been applied to analyze three dimensional characteristics of the flow field and compression wave around the high speed train which Is entering into a tunnel. The numerical scheme of AF + ADI was used to efficiently solve Navier-Stokes equations in the curvilinear coordinate system. The vortex formation around the nose region was found and the generation of compression wave due to the blockage effects was observed ahead of the train in the form of plane wave. The three dimensional characteristics of the flow field compared to the analytic results were discussed in detail. The variation of pressure of tunnel wall surface and velocity profile of the train are identified as the train enters into a tunnel. The changes in aerodynamic forces and streamlines of each specific sections are also discussed and presented.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.35-46
• /
• 1996

원심형 임펠러 내부 유로등 큰 곡률을 수반하는 터보기계 요소의 유동해석을 위한 계산코드를 개발하였다. 이 코드에서는 곡선좌표계에 유도된 3차원 비압축성 Navier-Stokes의 운동 방정식을 SMAC 음해법으로 푼다. 이 코드를 이용하여 유로의 단면이 정사각형이고 90도로 굽은 덕트내부의 층류 입구유동을 해석하고, 굽은 관 특유의 유동현상을 수치모사하였다 또한 곡관부 입구에서 충분히 발달한 유동, 또는 발달중인 유동이 유입될 경우에 이것이 곡관부 내부의 유동에 미치는 영향을 상·하류의 계산영역이 서로 다를 몇몇 유동장에 대하여 조사하고, 본 계산에서 얻어진 결과와 실형결과와의 비교로 본 3차원 유동해석 코드의 유효성을 검토 하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.3
• /
• pp.1015-1027
• /
• 1996

In this paper, a numerical procedure for the calculation of the incompressible Navier-Stokes equations in a generalized nonorthogonal body fitted coordinate system is proposed and is validated through three test problems. Present numerical procedure derives the pressure equation by using the pressure substitution method on the regular grid system, and discretized momentum equations are based on the covariant velocity components. Cavity flow, backward facing step flow, and two dimensional channel flow with a sinusoidal wavy wall are chosen as three test problems. Numerical solutions obtained by present procedure shows a good agreement with previous numerical and/or experimental results. Convergence rate is also satisfactory.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1122-1133
• /
• 2003

Numerical characteristics of implicit upwind schemes, such as upwind ADI, line Gauss-Seidel (LGS) and point Gauss-Seidel (LU) algorithms, for Navier-Stokes equations have been investigated. Time-derivative preconditioning method was applied for efficient convergence at low Mach/Reynolds number regime as well as at large grid aspect ratios. All the algorithms were expressed in approximate factorization form and von Neumann stability analysis was performed to identify stability characteristics of the above algorithms in the presence of high grid aspect ratios. Stability analysis showed that for high aspect ratio computations, the ADI and LGS algorithms showed efficient damping effect up to moderate aspect ratio if we adopt viscous preconditioning based on min-CFL/max-VNN time-step definition. The LU algorithm, on the other hand, showed serious deterioration in stability characteristics as the grid aspect ratio increases. Computations for several practical applications also verified these results.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.117-124
• /
• 2005

Splitting algorithms of the incompressible Navier-Stokes equations using P1P1/P2P1 finite element formulation are newly proposed. P1P1 formulation allocates velocity and pressure at the same nodes, while P2P1 formulation allocates pressure only at the vertex nodes and velocity at both the vertex and mid nodes. For comparison of the elapsed time and accuracy of the two methods, they have been applied to the well-known benchmark problems. The three cases chosen are the two-dimensional steady and unsteady flows around a fixed cylinder, decaying vortex, and impinging slot jet. It is shown that the proposed P2P1 semi-splitting method performs better than the conventional P1P1 splitting method in terms of both accuracy and computation time.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.9 no.4
• /
• pp.1043-1049
• /
• 2008

In this paper, we investigate the static response. natural frequencies and buckling loads of functionally graded material (FGM) plates, using a Navier method. The eigenvalues of the FGM plates and shells are calculated by varying the volume fraction of the ceramic and metallic constituents using a sigmoid function, but their Poisson's ratios of the FGM plates and shells are assumed to be constant. The expressions of the membrane. bending and shear stiffness of FGM plates art more complicated combination of material properties than a homogeneous element. In order to validate the present solutions, the reference solutions of rectangular plates based on the classical theory are used. The various examples of composite and FGM structures are presented. The present results are in good agreement with the reference solutions.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1726-1731
• /
• 2004

This paper presents the response surface optimization method using three-dimensional Navier-Stokes analysis to optimize the shape of a forward-curved blades centrifugal fan. For numerical analysis, Reynolds-averaged Navier-Stokes equations with $k-{\varepsilon}$ turbulence model are discretized with finite volume approximations. In order to reduce huge computing time due to a large number of blades in forward-curved blades centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models. Three geometric variables, i.e., location of cut off, radius of cut off, and width of impeller, and one operating variable, i.e., flow rate, were selected as design variables. As a main result of the optimization, the efficiency was successfully improved. And, optimum design flow rate was found by using flow rate as one of design variables. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time.