• 한국전산유체공학회지
• /
• 제9권3호
• /
• pp.57-65
• /
• 2004

The incompressible Navier-Stokes equations in two dimensions are stabilized by a modified residual method, and then discretized by hierarchical elements. The stabilization is necessary to escape from the Ladyzhenskaya-Babuska-Brezzi(LBB) constraint and hence to achieve an equal order formulation. To expedite a standard iterative method such as the conjugate gradient squared(CGS) method, a preconditioning technique called the Hierarchical Iterative Procedure(HIP) has been applied. In this paper, we increased the order of interpolation within an element up to cubic. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far The numerical results by the present HIP for the lid driven cavity flow and others showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• 대한기계학회논문집B
• /
• 제27권4호
• /
• pp.458-465
• /
• 2003

The main objective of the research is to develop a research code solving transient incompressible Navier-Stokes equation. In this research code, Adams-Bashforth method was applied to the convective terms of the navier stokes equation and the splitted equations were discretized spatially by finite element methods to solve the complex geometry problems easily. To reduce the divergence on the boundaries of pressure poisson equation due to the unsuitable pressure boundary conditions, multi step approximation pressure boundary conditions derived from the boundary linear momentum equations were used. Simulations of Lid Driven Flow and Flow over Cylinder were conducted to prove the accuracy by means of the comparison with results of the previous workers.

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

• International Journal of Automotive Technology
• /
• 제7권2호
• /
• pp.139-143
• /
• 2006

As an automobile tends to be high grade, the needs for more luxurious interior and comfortable HVAC system are emerged. The defrosting ability is another major factor of the performances of HVAC system. The present work is to simulate the flow and the temperature field of cabin interior during the defrost mode. The three-dimensional incompressible Navier-Stokes equations and energy equation were solved on the multi blocked grid system by the iterative time marching method and AF scheme, respectively. The present computations were validated by the comparison of the temperature field of a driven cavity and velocity field of 1/5 model scale of an automobile. Generally good agreements were obtained. By the present computation, the complicated features of flow and temperature within the automotive cabin interior could be well understood.

• 대한기계학회논문집B
• /
• 제21권3호
• /
• pp.341-349
• /
• 1997

The characteristics of defect correction method are discussed in a sample heat conduction problem showing the numerical solution of the error correction equation can predict the error of the numerical solution of the original governing equation. A way of using defect correction method combined with the existing algorithm for the incompressible fluid flow, is proposed and subsequently tested for the driven square cavity problem. The error correction equations for the continuity equation and the momentum equations are considered to estimate the errors of the numerical solutions of the original governing equations. With this new approach, better velocity and pressure fields can be obtained by correcting the original numerical solutions using the estimated errors. These calculated errors also can be used to estimate the orders of magnitude of the errors of the original numerical solutions.

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

The present work was undertaken to numerically analyze the defrosting phenomena of windshield glass. In order to analysis the phase change from frost to water on windshield glass by discharging hot air from a defroster nozzle, the flow and the temperature field of the cabin interior, the heat transfer through the windshield glass, and the phase change of frost should be solve simultaneously. In the present work, the flow field was obtained by solving 3-D incompressible Navier-Stokes equations, and the temperature field was computed from the incompressible energy equation. The phase change process was solved by the enthalpy method. For the code validation, the temperature and the phase change of the driven cavity were calculated. The calculation showed a good agreement with other numerical results. Then, the present code was applied to the defrosting problem of a real automobile, and a good agreement with the experimental data was also obtained.

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

This paper is a continuation of the recent development on the hermite-based divergence free basis function and deals with a non-isothermal fluid flow thru the buoyancy driven flow in a square cavity with temperature difference across the two sides. The basis functions for the velocities consist of the hermite function and its curl. However, the basis for the temperature are the hermite function and its gradienst. Hence, the number of degrees of freedom at a node becomes 6, which are the stream function, two velocities, the temperature and its x- and y-derivatives. Numerical results for the streamlines, the temperatures, the x-velocities and the y-velocities show good agreements with those of De vahl Davis[7].

• 대한기계학회논문집B
• /
• 제27권4호
• /
• pp.534-541
• /
• 2003

To investigate the combustion characteristics in a low-pollutant municipal waste incinerator, the generalized multi-block simulation code that can be applied to turbulent reacting flow with gaseous hydrocarbon fuel in a 3D complex geometry has been developed with nongray radiation effects. To deal with the complex geometry, structured multi-block method and the scheme which treats interfaces implicitly are adopted. The developed code is validated through various engineering problems such as curved duct flow, driven cavity flow, gray multi-block radiation, nongray radiation. and combustion in a incinerator.

• 한국전산유체공학회지
• /
• 제18권3호
• /
• pp.79-86
• /
• 2013

Preconditioned compressible Navier-Stokes equations are solved for almost incompressible flows. Unstructured meshes are utilized for spatial discretization of complex flow domain. Effectiveness of the current preconditioning algorithm, with respect to various Reynolds numbers and Mach numbers, is demonstrated by the solution of canonical problems for incompressible flows, e.g. driven cavity flows.

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

Objective of the present study is to examine the convergence characteristics of the various multi-stage Runge-Kutta methods in solving the incompressible Navier-Stokes equations of a time-marching from casted by the artificial compressibility method. Convergence characteristics are examined over 2-stage, 4-stage and hybrid type (using 4-, 3-, 2-stages sequentially) Runge-Kutta methods for a laminar lid-driven cavity flow, and also for a turbulent bump channel flow using Chien's low-Reynolds number turbulence model. Efforts are made to establish a stable and fast convergent multi-stage Runge-Kutta method with minimal artificial dissipations.