• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.528-531
• /
• 2004

For a simple one-dimensional ignition problem a mathematical model is described to investigate the difficulties in numerical simulations. Some computation results are obtained and comparison is made with analytical solution. Discussions are made on topics such as 1) coordinate transformation, 2) gas-phase and solid-phase analysis; (divergence form of the governing system, a finite-volume discretization, implicit time integration, upwind split flux, spatial accuracy improvement are described. Mass, reagent mass, and energy conservations are solved.), and 3) method to determine quantities on the burning surface (matching). Results obtained for small values of the non-dimensional pressure show a steady-combustion and good agreement with the analytical solution. Numerical instability appeared for larger values of the pressure, discussion on the cause of the problem is made. This effort is a part of a study of flame spread phenomena on solid propellant surface.

• Structural Engineering and Mechanics
• /
• v.3 no.2
• /
• pp.135-144
• /
• 1995

Since structural systems may fail in any one of several failure modes, computation of system reliability is always difficult. A method using numerical quadrature for computing structural system reliability with either one or more than one failure mode is presented in this paper. Statistically correlated safety margin equations are transformed into a group of uncorrelated variables and the joint density function of these uncorrelated variables can be generated by using the Maximum Entropy Method. Structural system reliability is then obtained by integrating the joint density function with the transformed safety domain enclosed within a set of linear equations. The Gaussian numerical integration method is introduced in order to improve computational accuracy. This method can be used to evaluate structural system reliability for Gaussian or non-Gaussian variables with either linear or nonlinear safety boundaries. It is also valid for implicit safety margins such as computer programs. Both the theory and the examples show that this method is simple in concept and easy to implement.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2012.05a
• /
• pp.187-195
• /
• 2012

This paper deals with numerical computations of cryogenic flows around turbopump inducer. Firstly, we introduce numerical methods to compute compressible/incompressible cryogenic two-phase flow. As a validation problem, computation results of 2 dimensional/axi-symmetric cryogenic flow will be presented. In this process, various cavitation model will be compared. Finally, numerical simulation of 3 dimensional turbopump inducer will be presented.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.52-61
• /
• 2014

The present paper deals with the efficient computation of higher-order CFD methods for compressible flow using graphics processing units (GPU). The higher-order CFD methods, such as discontinuous Galerkin (DG) methods and correction procedure via reconstruction (CPR) methods, can realize arbitrary higher-order accuracy with compact stencil on unstructured mesh. However, they require much more computational costs compared to the widely used finite volume methods (FVM). Graphics processing unit, consisting of hundreds or thousands small cores, is apt to massive parallel computations of compressible flow based on the higher-order CFD methods and can reduce computational time greatly. Higher-order multi-dimensional limiting process (MLP) is applied for the robust control of numerical oscillations around shock discontinuity and implemented efficiently on GPU. The program is written and optimized in CUDA library offered from NVIDIA. The whole algorithms are implemented to guarantee accurate and efficient computations for parallel programming on shared-memory model of GPU. The extensive numerical experiments validates that the GPU successfully accelerates computing compressible flow using higher-order method.

• Journal of the Society of Naval Architects of Korea
• /
• v.44 no.1 s.151
• /
• pp.8-15
• /
• 2007

In this paper, numerical calculations are performed to analyze the unsteady flow of NACA airfoil sections. In order to ease the flow computation for the fluid region changing in time, improve the quality of solution and simplify the grid generation for the oscillating foil flow, the computational method adopts a moving and deforming mesh with the multi-block grid topology. The multi-block, structured-unstructured hybrid grid is generated using the commercial meshing software Gridgen V15. The MDM (Moving & Deforming Mesh) and the UDF (User Define function) function of FLUENT 6 are adopted for computing turbulent flows of the foil in pitching motion. Computed unsteady lift and drag forces are compared with experimental data. in general, the characteristics of unsteady lift and drag of the experiments are reproduced well in the numerical analysis.

• Journal of Hydrogen and New Energy
• /
• v.31 no.3
• /
• pp.293-301
• /
• 2020

This paper describes a numerical study on the improvement of performance of the 2 vane pump impellers. The design of these impellers was optimized using a commercial computation fluid dynamics code and design of experiments. Geometric design variables were defined by the impeller blade angle distribution. The objective functions were defined as the total head, total efficiency and solid material size of the impellers. The importance of the geometric design variables was analyzed using 2^k factorial designs. The interaction between the total head, total efficiency and solid material size, according to the impeller blade angle distribution, is discussed by analyzing the 2^k factorial design results.

• Journal of Korea Water Resources Association
• /
• v.34 no.1
• /
• pp.91-103
• /
• 2001

In this studs, hydraulic routing model has been developed to predict the water level and discharge in each river section with considering the full interaction between surface runoff and river flow. It improved the computation of flood runoff by reflecting the shape of hydrograph that was determined by the geological and flood characteristics, and the excessive computation of the peak discharge was eliminated by considering the effect of infiltration. The Inflow from surface runoff to river flow was applied to the equation of continuity by implementing effectively the flow in a number of river section, and resulted in a numerical stability at the rapid variation of rainfall. Measurements were conducted during heavy rain in the watershed area of Yang-Yang Namdae-Chun. The present model was tested to the field, and the computed results were compared to the observed data. Its applicability was confirmed with its verification.

• The Transactions of the Korea Information Processing Society
• /
• v.5 no.1
• /
• pp.97-102
• /
• 1998

In the most cases of eigen value problems in the social sciences, the object matrix to analyze is real-valued symmetric matrix. And many cases of eigen value problems in this field needs 2-4 eigen pairs according to the magnitude of their absolute values. The methods to obtain eigen pairs by numerical computation using computer, we would face the problem of round off error because matrix computation needs a number of calculations. In this paper, an algorithm which make us to get some needed eigcn pairs according to the magnitude of their absolute values is designed. And in this algorithm, the power method is used to obtain some eigen pairs. This algorithm is expected to be effective by the reduction of the number of calculations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.609-616
• /
• 2005

While the SIMPLE algorithm is most widely used for the simulations of flow phenomena that take place in the industrial equipments or the manufacturing processes, it is less adopted for the simulations of the free surface flow. Though the SIMPLE algorithm is free from the limitation of time step, the free surface behavior imposes the restriction on the time step. As a result, the explicit schemes are faster than the implicit scheme in terms of computation time when the same time step is applied to, since the implicit scheme includes the numerical method to solve the simultaneous equations in its procedure. If the computation time of SIMPLE algorithm can be reduced when it is applied to the unsteady free surface flow problems, the calculation can be carried out in the more stable way and, in the design process, the process variables can be controlled based on the more accurate data base. In this study, a modified SIMPLE algorithm is presented fur the free surface flow. The broken water column problem is adopted for the validation of the modified algorithm (MoSIMPLE) and for comparison to the conventional SIMPLE algorithm.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1509-1517
• /
• 1993

The purpose of this study is to develop a program for the automatic derivation of the symbolic equations necessary for the kinematic and dynamic analyses of the spatial mechanism. For this purpose, a symbolic manipulation package called MCSYMA is used. Every symbolic equation is formulated using relative joint coordinate to obtain the numerically efficient system equations. These equations are produced in FORTRAN statements and linked to a FORTRAN program for numerical analysis. Several examples are taken for comparison with the commercial package called DADS which is using Cartesian coordinate approach. Also, this symbolic formulation approach is compared with a conventional numerical approach for an example. The results show that this symbolic approach with relative joint coordinate system is most efficient in computational time among three and is recommended for the derivation of macro elements frequently used.