• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.900-907
• /
• 2000

The three-dimensional turbulent flow in curved pipes susceptible to flow-accelerated corrosion has been analyzed numerically to predict the pressure and shear stress distributions on the inner surface of the pipes. The analysis employs the body-fitted non-orthogonal curvilinear coordinate system and a standard $ {\kappa}-{\varepsilon}$ turbulence model with wall function method. The finite volume method is used to discretize the governing equations. The convection term is approximated by a high-resolution and bounded discretization scheme. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of a modified version of momentum interpolation scheme. The SIMPLE algorithm is employed for the pressure and velocity coupling. The numerical calculations have been performed for two curved pipes with different bend angles and curvature radii, and discussions have been made on the distributions of the primary and secondary flow velocities, pressure and shear stress on the inner surface of the pipe to examine applicability of the present analysis method. As the result it is seen that the method is effective to predict the susceptible systems or their local areas where the fluid velocity or local turbulence is so high that the structural integrity can be threatened by wall thinning degradation due to flow-accelerated corrosion.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.6
• /
• pp.753-765
• /
• 2003

A finite element code for the numerical solution of the Navier-Stokes equation is parallelized by vertex-oriented domain decomposition. To accelerate the convergence of iterative solvers like conjugate gradient method, parallel block ILU, iterative block ILU, and distributed ILU methods are tested as parallel preconditioners. The effectiveness of the algorithms has been investigated when P1P1 finite element discretization is used for the parallel solution of the Navier-Stokes equation. Two-dimensional and three-dimensional Laplace equations are calculated to estimate the speedup of the preconditioners. Calculation domain is partitioned by one- and multi-dimensional partitioning methods in structured grid and by METIS library in unstructured grid. For the domain-decomposed parallel computation of the Navier-Stokes equation, we have solved three-dimensional lid-driven cavity and natural convection problems in a cube as benchmark problems using a parallelized fractional 4-step finite element method. The speedup for each parallel preconditioning method is to be compared using upto 64 processors.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

• Journal of Energy Engineering
• /
• v.21 no.3
• /
• pp.228-236
• /
• 2012

Large eddy simulation (LES) is increasingly used as a tool for studying the dynamics of turbulence in combustion chamber flows due to the promise of wider generality and more accurate results compared to Reynolds averaged Navier-Stokes(RANS) models. This study presents the appropriate subgrid-scale(SGS) model in LES for predicting the turbulent flow field in the internal combustion engine. The study of the effects of model and numerical parameters such as discretization scheme, initial condition, time step and SGS model was performed. The results of LES using the SGS model were found to be in the good agreement with experimental data.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2011.04a
• /
• pp.443-444
• /
• 2011

This paper discusses a model order reduction for large order rotor dynamics systems results from the finite element discretization. Typical rotor systems consist of a rotor, built-on parts, and a support system, and require prudent consideration in their dynamic analysis models because they include unsymmetric stiffness, localized nonproportional damping and frequency dependent gyroscopic effects. When the finite element model has a very large number of degrees of freedom because of complex geometry, repeated dynamic analyses to investigate the critical speeds, stability, and unbalanced response are computationally very expensive to finish within a practical design cycle. In this paper, the Krylov-based model order reduction via moment matching significantly speeds up the dynamic analyses necessary to check eigenvalues and critical speeds of a Nelson-Vaugh rotor system. With this approach the dynamic simulation is efficiently repeated via a reduced system by changing a running rotational speed because it can be preserved as a parameter in the process of model reduction. The Campbell diagram by the reduced system shows very good agreement with that of the original system. A 3-D finite element model of the Nelson-Vaugh rotor system is taken as a numerical example to demonstrate the advantages of this model reduction for rotor dynamic simulation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.5
• /
• pp.787-797
• /
• 1986

An implicit approach is employed to obtain a general boundary integral formulation of axisymmetric elastic problems in terms of a pair of singular integral equations. The corresponding kernel functions from the solutions of Navier's equation are derived by applying a three dimensional integral and a direct axisymmetrical approach. A numerical discretization schem including the evaluation of Cauchy principal values of the singular integral is described. Finally the typical axisymmetric elastic models are analyzed, i.e. the hollow sphere, the constant thickness and the V-notched round bar.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.8
• /
• pp.8-17
• /
• 2005

It is well known that preconditioning methods are efficient for convergence acceleration in the compressible low Mach number flows. In this study, the original Euler equations and three differently nondimensionalized preconditioning methods are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as numerical flux discretization and time integration, respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning allows Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.8
• /
• pp.790-797
• /
• 2010

Mode transition of the axi-symmetric screech tone in the low supersonic Mach number range from 1.0 to 1.20 is numerically analyzed. The axi-symmetric Navier-Stokes equations and the k-e turbulence model are solved in the cylindrical coordinate system. The dispersion-relation-preserving(DRP) scheme is applied for space discretization and the optimized four levels marching method are used for time integration. At low supersonic Mach numbers with an axi-symmetric A1 mode in the simulation, it is shown that acoustic propagation due to the nonlinear effects is seen in the lateral direction and the screech tone frequency is the same as the vortex passing frequency due to the generation of intense large-scale vortical motions.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.2 no.5
• /
• pp.58-65
• /
• 1994

In order to reduce particulate emissions from diesel vehicles, mathematical model is established and analyzed on ceramic wall-flow monolith filter. A wall-flow monolith filter placed in the exhaust stream of a diesel engine can effectively limit the emission of diesel particulates through the monolith. The accumulated particulates can then be periodically combusted inside the monolith by directing hot gas to the monolith while normal engine exhaust is routed around the monolith system. The resulting low flow rates through the monolith require consideration of gas dynamics through the channels as well as particulate combustion to analyze this regeneration process. A mathematical model of the regeneration is formulated as a system of nonlinear partial differential equations describing the conservation of mass, momentum and energy. Numerical solutions are obtained by using a finite difference techniques for the spatial discretization. So we can use filter simulation program for the purpose of filter design and actual filter regeneration

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.1
• /
• pp.18-30
• /
• 1999

A heat transfer analysis for the two-dimensional (2-D) steady state using finite difference method (FDM) is performed to predict the thermal behavior of the primary first-wall (FW) system of fusion reactor under various geometric and thermo-hydraulic conditions, such as the beryllium (Be) armor thickness, pitch of cooling tube, and coolant velocity. The FW consists of authentic steel (type 316 stainless steel solution annealed) for cooling tubes, Cu for cooling tubes embedding material, and Be for a protective armor, based on the International Thermonuclear Experiment Reactor (ITER) report. The present 2-D analysis, the control volume discretized with hybrid grid (rectangular grid and polar grid) and Gauss-Seidel iteration method are adapted to solve the governing equations. In the present study, geometric and thermo-hydraulic parameters are optimized with consideration of several limitations. Consequently, it is suggested that the adequate pitch of cooling tube is 22-32mm, the beryllium armor thickness is 10-12mm, and that the coolant velocity is 4.5m/s-6m/s for $100^{\circ}C$ of inlet coolant temperature. The cooling tube should locate near beryllium armor. But, it would be better for locating the center of Cu wall, considering problems of material and manufacturing. Also, 2-D analysis neglecting the axial temperature distribution of cooling tube is appropriate, regarding the discretization error in axial direction.