• 한국전산유체공학회지
• /
• 제3권2호
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.

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

An implicit turbulent flow solver is developed for 2-D unstructured hybrid meshes. Spatial discretization is accomplished by a cell-centered finite volume formulation using an upwind flux differencing. Time is advanced by an implicit backward Euler time stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one equation model, which is coupled with wall function. The numerical method is applied for flows on a flat plate, the NACA 0012 airfoil, and the Douglas 3 element airfoil. The results are compared with experimental data.

• Kyungpook Mathematical Journal
• /
• 제52권3호
• /
• pp.327-346
• /
• 2012

In this paper, we present a simple explicit type numerical method for discretizations in time for solving one dimensional Burgers' equations. The proposed method does not need an iteration process that may be required in most implicit methods and have good convergence and efficiency in computational sense compared to other known numerical methods. For evidences, several numerical demonstrations are also provided.

• International Journal of Precision Engineering and Manufacturing
• /
• 제3권4호
• /
• pp.54-63
• /
• 2002

In recent years, mechanical systems such as high speed vehicles and railway trains moving on elastic beam structures have become a very important issue to consider. In this paper, a general approach, which can predict the dynamic behavior of a constrained mechanical system moving on a flexible beam structure, is proposed. Various supporting conditions for the foundation support are considered for the elastic beam structure. The elastic structure is assumed to be a non-uniform and linear Bernoulli-Euler beam with a proportional damping effect. Combined differential-algebraic equation of motion is derived using the multi-body dynamics theory and the finite element method. The proposed equations of motion can be solved numerically using the generalized coordinate partitioning method and predictor-corrector algorithm, which is an implicit multi-step integration method.

• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.66-72
• /
• 1997

Optimum nozzle design exploiting the method of characteristic(M.O.C) has been in application as an efficient design methodology targeting a less weighted and short expansion nozzle. This paper treats the optimum nozzle design and the analysis of the inviscid compressible flow inside. Based on traditional Rao's method, the optimum nozzle design is coded with minor modifications for the identification of the control surface across which the mass flux should be conserved. Internal flow field is simulated numerically by M.O.C and implicit/explicit Taylor-Galerkin finite element method(F.E.M) with the aid of adaptive remeshing to capture the shock wave, hence improve the accuracy. Designed and calculated flow fields due to the separate analyses show that the mass flux predicted by optimum nozzle design with M.O.C is not conserved across the control surface and the sonic line should be located upstream of the nozzle throat. Rao's optimum nozzle design methodology exaggerates the momentum thrust and tends to overemphasize the engine performance loss.

• 한국전산유체공학회지
• /
• 제4권2호
• /
• pp.67-73
• /
• 1999

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler, thin layer Navier-Stokes and full Navier-Stokes ones. are solved using implicit LU-ADI decomposition scheme. The gradient projection method with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

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

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler , thin layer Navier- Stokes and full Navier-Stokes ones, are solved using implicit LU-ADI decomposition scheme. The feasible direction algorithm with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

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

The delta-formulation of the Navier-Stokes equations has been popularly used in the aerodynamics area. Implicit algorithm can be easily implemented in that by using Taylor series expansion. This formulation is extended for an unsteady analysis by using a dual-time integration. In the meanwhile, the incompressible flows with heat transfers which occur in the area of thermo-hydraulics have been solved by a segregated algorithm such as the SIMPLE method, where each equation is discretised by using an under-relaxed deferred correction method and solved sequentially. In this study, the dual-time delta formulation is implemented in the segregated Navier-Stokes solver which is based on the collocated cell-centerd scheme with un unstructured mesh FVM. The pressure correction equation is derived by the SIMPLE method. From this study, it was found that the Euler dual-time method in the delta formulation can be combined with the SIMPLE method.

• Journal of Mechanical Science and Technology
• /
• 제20권11호
• /
• pp.1961-1971
• /
• 2006

A soft recovery system (SRS) is a device that stops a high speed projectile without damaging the projectile. The SRS is necessary to verify the shock resistant requirements of microelectronics and electro-optic sensors in smart munitions, where the projectiles experience over 20,000 g acceleration inside the barrel. In this study, a computer code for the performance evaluation of a SRS based on ballistic compression decelerator concept has been developed. It consists of a time accurate compressible one-dimensional Euler code with use of deforming grid and a projectile motion analysis code. The Euler code employs Roe's approximate Riemann solver with a total variation diminishing (TVD) method. A fully implicit dual time stepping method is used to advance the solution in time. In addition, the geometric conservation law (GCL) is applied to predict the solutions accurately on the deforming mesh. The equation of motion for the projectile is solved with the four-stage Runge-Kutta time integration method. A small scale SRS to catch a 20 mm bullet fired at 500 m/s within 1,600 g-limit has been designed with the proposed method.

• 한국정밀공학회지
• /
• 제17권11호
• /
• pp.165-175
• /
• 2000

In recent years, it becomes a very important issue to consider the mechanical systems such as high-speed vehicles and railway trains moving on elastic beam structures. In this paper, a general approach, which can predict the dynamic behavior of constrained mechanical system and elastic beam structure, is proposed. Also, various supporting conditions of a foundation support are considered for the elastic beam structures. The elastic structure is assumed to be a nonuniform and linear Bernoulli-Euler beam with proportional damping effect. Combined Differential-Algebraic Equations of motion are derived using multibody dynamics theory and Finite Element Method. The proposed equations of motion can be solved numerically using generalizd coordinate partitioning method and Predictor-Corrector algorithm, which is an implicit multi-step integration method.