• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.6
• /
• pp.1921-1930
• /
• 1996

To design an optimum engine intake system, a flow model for the intake manifold was developed by the finite difference method. The flow in the intake manifold was one-dimensional, and the finite difference equations were derived from governing equations of flow, continuity, momentum and energy. The thermodynamic properties of the cylinder were found by the first law of thermodynamics, and the boundary conditions were formulated using steady flow model. By comparing the calculated results with experimental data, the appropriate boundary conditions and convergence limits for the flow model were established. From this model, the optimum manifold lengths at different engine operating conditions were investigated. The optimum manifold length became shorter when the engine speeds were increased. The effect of intake valve timings on inlet air mass was also studied by this model. Advancing intake valve opening decreased inlet air mass slightly, and the optimum intake valve closing was found. The difference in inlet air mass between cylinders was very small in this engine.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1624-1629
• /
• 2008

The governing equations in generalized curvilinear coordinates for a 3D pulsatile flow are the Incompressible Navier-Stokes (INS) equations with the artificial dissipative terms and continuity equation discretized using a second-order accurate, finite volume method on the nonstaggered computational grid. This method adopts a dual or pseudo time-stepping Artificial Compressibility (AC) method integrated in pseudo-time. The computational technique implements the implicit approximate factorization method of the Beam and Warming method (1978), which is the extension of the Alternate Direction Implicit (ADI) method. The algorithm yields practically identical velocity profiles and secondary flows that are in excellent overall agreement with an experimental measurement (Rindt & Steenhoven, 1991).

• Water for future
• /
• v.24 no.1
• /
• pp.73-81
• /
• 1991

The Saint-Venant equations of the continuity and momentum principles of one-dimensional, unsteady, open-channel flow are expressed in terms of the phase velocities of constant depth, velocity, and discharge, which results in unique relationships between these phase velocities and channel velocity. A case study shows that these unique relationships developed in this study can be used as an indicator of precision of numerical solutions of the Saint-Venant equations. Further physical interpretation of these relationships and utilization to the numerical analyses of the Saint-Venant equations are to be investigated.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.5
• /
• pp.36-42
• /
• 2003

A CFD simulation technique has been developed to handle the unsteady body motion with large amplitude by use of overlapping multi-block grid system. The three-dimensional, viscous and incompressible flow around body is investigated by solving the Navier-Stokes equations, and the motion of body is represented by moving effect of the grid system. Composite grid system is employed in order to deal with both the body motion with large amplitude and the condition of numerical wave maker in convenience at the same time. The governing equations, Navier-Stokes (N-S) and continuity equations, are discretized by a finite volume method, in the framework of an O-H type boundary-fitted grid system (inner grid system including test model) and a rectangular grid system (outer grid system including simulation equipments for generation of wave environments). If this study, several flow configurations, such as an oscillating cylinder with large KC number, are studied in order to predict and evaluate the hydrodynamic forces. Furthermore, the motion simulation of a Series 60 model advancing in a uniform flow under the condition of enforced roll motion of angle 20$^{\circ}$ is performed in the developed numerical wave tank.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.8
• /
• pp.1950-1963
• /
• 1995

A modification of the approximate-factorization method is made to accelerate the convergency rate and to take sufficiently large Courant number without loss of accuracy. And a stable implicit finite-difference scheme for solving the incompressible Navier-Stokes equations employed above modified method is developed. In the present implicit scheme, the volume fluxes with contravariant velocity components and the pressure formulation in curvilinear coordinates is adopted. In order to satisfy the continuity condition completely and to remove spurious errors for the pressure, the Navier-Stokes equations are solved by a modified SMAC scheme using a staggered gird. The upstream-difference scheme such as the QUICK scheme is also employed to the right hand side. The implicit scheme is unconditionally stable and satisfies a diagonally dominant condition for scalar diagonal linear systems of implicit operator on the left hand side. Numerical results for some test calculations of the two-dimensional flow in a square cavity and over a backward-facing step are obtained using both usual approximate-factorization method and the modified one, and compared with each other. It is shown that the present scheme allows a sufficiently large Courant number of O(10$^{2}$) and reduces the computing time.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.4
• /
• pp.839-846
• /
• 1994

This paper presents a simultaneous solution algorithm for one-dimensional unsteady flow routing through a dendritic channel system. This simulations solution algorithm is based on the double-sweep method and utilizes separate recursion equations for continuity, momentum and energy equations for each of the individual components of a dendritic channel system. Through separate recursion equations for each of the components. the new algorithm converts a dendritic channel network problem into a single-channel problem. The new algorithm is utilized in conjunction with a linearized unsteady flow model using full dynamic flow equations. The required computer storage for the coefficient matrix of the whole system is reduced significantly from the $2N{\times}2N$ matrix to a $2N{\times}4$ matrix, where N is the number of cross sections used in the computation of flow variables in a dendritic channel system. The algorithm presented in this paper provides an efficient and accurate modeling of unsteady flow events through a dendritic channel system.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2005.06a
• /
• pp.965-968
• /
• 2005

Roll drafting operation causes variations in the linear density of bundles because the bundle flow cannot be controlled completely by roll pairs. Defects occurring in this operation bring about many problems successively in the next processes. In this paper, we attempt to analyze the draft dynamics and the linear density irregularity based on the governing equation of a bundle motion that has been suggested in our previous studies. For analyzing the dynamic characteristics of the roll drafting operation, it is indispensable to investigate a transient state in time domain before the bundle flux reaches a steady state. However, since governing equations of bundle flow consisting of continuity and motion equations turn out to be nonlinear, and coupled between variables, the solutions for a transient state cannot be obtained by an analytical method. Therefore, we use the Finite Difference Method(FDM), particularly, the FTBS(Forward-Time Backward-Space) difference method. Then, the total equations system yields to an algebraic equations system and is solved under given initial and boundary conditions in an iterative fashion. From the simulation results, we confirm that state variables show different behavior in the transient state; e.g., the velocity distribution in the flow field changes more quickly the linear density distribution. During a transient flow in a drafting zone, the output irregularity is influenced differently by the disturbances, e.g., the variation in input bundle thickness, the drafting speed, and the draft ratio.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.8 no.3
• /
• pp.277-300
• /
• 2016

The motion of a floating body and the free surface flow are the most important design considerations for ships and offshore platforms. In the present research, a numerical method is developed to simulate the motion of a floating body and the free surface using a fixed rectilinear grid system. The governing equations are the continuity equation and Naviere-Stokes equations. The boundary of a moving body is defined by the interaction points of the body surface and the centerline of a grid. To simulate the free surface the Modified Marker-Density method is implemented. Ships advancing in regular waves, the interaction of waves by a fixed circular cylinder array and the response amplitude operators of an offshore platform are simulated and the results are compared with published research data to check the applicability. The numerical method developed in this research gives results good enough for application to the initial design stage.

• Journal of Mechanical Science and Technology
• /
• v.19 no.12
• /
• pp.2263-2269
• /
• 2005

The ship wave phenomena in the restricted waterway were investigated by a numerical analysis. The Euler and continuity equations were employed for the present study. The boundary fitted and moving grid system was adopted to enhance the computational efficiency. The convective terms in the governing equations and the kinematic free surface boundary condition were solved by the Constrained Interpolated Profile (CIP) algorithm in order to solve accurately wave heights in far field as well as near field. The advantage of the CIP method was verified by the comparison of the computed results by the CIP and the Maker and Cell (MAC) method. The free surface flow simulation around Wigley hull was performed and compared with the experiment for the sake of the validation of the numerical method. The present numerical scheme was applied to the free surface simulation for various canal sections in order to understand the effect of the sectional shape of waterways on the ship waves. The wave heights on the side wall and the shape of the wave patterns with their characteristics of flow are discussed.

• Transactions of the Society of Information Storage Systems
• /
• v.5 no.1
• /
• pp.19-35
• /
• 2009

The present work is an experimental and analytical study on a flexible disk rotating close to a rigid rotating disk in open air. In the analytical study, the air flow in the gap between the flexible disk and the rigid disk is modeled using Navier-Stokes and continuity equations while the flexible disk is modeled using the linear plate theory. The flow equations are discretized using the cell centered finite volume method (FVM) and solved numerically with semi-implicit pressure-linked equations (SIMPLE algorithm). The spatial terms in the disk equation are discretized using the finite difference method (FDM) and the time integration is performed using fourth-order Runge-Kutta method. An experimental test-rig is designed to investigate the dynamics of the flexible disk when rotating close to a co-rotating, a counter-rotating and a fixed rigid disk, which works as a stabilizer. The effects of rotational speed, initial gap height and inlet-hole radius on the flexible disk displacement and its vibration amplitude are investigated experimentally for the different types of stabilizer. Finally, the analytical and experimental results are compared.