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

• 한국유체기계학회 논문집
• /
• 제7권3호
• /
• pp.32-38
• /
• 2004

Shape of a multi-blades centrifugal fan is optimized by response surface method based on three-dimensional Navier-Stokes analysis. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard $k-{epsilon}$ turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in this centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Optimizations with and without constraints are carried out. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. The correlation of efficiency with relative size of inactive zone at the exit of impeller is discussed as well as with average momentum fluxes in the scroll.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 추계 학술대회논문집
• /
• 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.

• 한국항해항만학회지
• /
• 제27권5호
• /
• pp.473-478
• /
• 2003

평행한 벽 사이에 원형배관(circular pipeline)을 놓고 그 주위의 유동특성에 대한 수치연구를 수행하였다. 비압축성 유체를 가지고, Navier-Stokes 방정식을 풀었고 3차 풍상(upwind) 차분의 수치해법을 이용하였다. 한쪽 벽과의 거리가 매우 작아질 때, 볼텍스 떨어짐이 상당히 억압되는데 이것은 벽 경계와의 상호 박리 작용 때문으로 간수된다. 본 연구는 레이놀드 수의 변함과 물체가 벽에 접근함에 따른 볼텍스 떨어짐의 특성을 규명하는데 있다. 원형배관 후류와 평행벽내 유기된 박리의 상호작용을 집중적으로 다루며 서로 다른 조건에서 박리와 와역한(vorticity dynamics)의 특성을 해석하였다.

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

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.

• 한국전산유체공학회지
• /
• 제11권4호
• /
• pp.75-80
• /
• 2006

Navier-Stokes simulation of the flow in a micro viscous pump is carried out. The micro viscous pump consists of a rotating circular rotor placed in a two dimensional channel. All simulation is carried out by using a finite volume approach, at the Reynolds number of 0.5, to study the performance of the micro viscous pump. Length of channel of the pump is varied to simulate the effects of the pumping load. Numerical solutions show that the net flow of the pump is realized by two counter rotating vortices formed on both sides of the rotor. The volume flow rate of the pump is decreased as length of the channel is increased, while the static pressure difference across the rotor is increased. The static pressure difference across the rotor is observed to be inversely proportional to the volume flow rate as inertia effects are negligibly small. The efficiency of the pump is found to reach a maximum when two counter rotating vortices on both sides of the rotor becomes to merge forming an outer enveloping vortex.

• 한국전산유체공학회지
• /
• 제19권1호
• /
• pp.27-33
• /
• 2014

A new and efficient implicit scheme is proposed to obtain a steady-state solution in time integration and the comparison of characteristics with the approximation ways for the implicit method to solve the incompressible Navier-Stokes equations is provided. The conservative, finite-volume cell-vertex upwind scheme and artificial compressibility method using dual time stepping for time accuracy is applied in this paper. The numerical results obtained indicate that the direct application of Jacobian matrix to the Lower and upper sweeps of implicit LU-SGS leads to better performance as well as convergence regardless of CFL number and true time step than explicit scheme and approximation of Jacobian matrix. The flow simulation around box in uniform flow with unstructured meshes is demonstrated to check the validity of the current formulation.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2002년도 춘계학술대회논문집
• /
• pp.140-146
• /
• 2002

The fluid induced vibration (FIV) phenomena of a 2-D.O.F airfoil system have been investigated in low Reynolds number incompressible flow region. Unsteady flows with viscosity are computed using two-dimensional incompressible Navier-stokes code. To validate developed Navier-Stokes code, steady and unsteady flow fields around airfoil are analyzed. The present fluid/structure interaction analysis is based on the most accurate computational approach with computational fluid dynamics (CSD) and computational structural dynamics (CSD) techniques. The highly nonlinear fluid/structure interaction phenomena due to severe flow separations have been analyzed fur the low Reynolds region (R$_{N}$ =500~5000) that has a dominancy of flow viscosity. The effect of R$_{N}$ on the fluid/structure coupled vibration instability of 2-DOF airfoil system is presented and the effect of initial angle of attack on the dynamic instability are also shown.own.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• 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.

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

Shape of a multi-blades centrifugal fan is optimized by response surface method based on three-dimensional Navier-Stokes analysis. For numerical analysis, Reynolds-averaged Wavier-Stokes equations with standard $k-{\varepsilon}$ turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in this centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Optimizations with and without constraints are carried out. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. The correlation of efficiency with relative size of inactive zone at the exit of impeller is discussed as well as with average momentum fluxes in the scroll.