• Transactions of the Korean Society of Automotive Engineers
• /
• v.22 no.5
• /
• pp.13-19
• /
• 2014

Urea-SCR system is the selective catalytic reduction to reduce nitrogen oxides ($NO_x$) emitted from diesel vehicles. The objective of this study is numerical analysis of 3-dimensional unsteady melting problems of frozen urea by using an electric heater. It can be applied to determine capacity of power with respect to time and the location of the urea suction pipe in urea storage tank. The study includes the change of liquid volume fraction, temperature profiles and a influence of natural convection by using the commercial software STAR-CCM+(v7.06). The accuracy of the numerical analysis is estimated by comparisons with experimental data. After validation, a numerical analysis for freezing urea is conducted with four different heating power. From the results, it was found that relation of velocity of phase interface and amount of melting urea by increasing heating power in a container. There is also a difference in trend between velocity of phase interface and amounts of melting urea because of effect of natural convection.

• Korea-Australia Rheology Journal
• /
• v.13 no.3
• /
• pp.131-139
• /
• 2001

A coupled approach is proposed for the prediction of sheet profile in sheet casting process, which combines one-dimensional analytical method on planar elongational flow region and three-dimensional numerical method on the other region. The strategy is constructed from the observations that the flow domain of sheet casting process can be separated into two parts based old the flow kinematics. The flow field in the central region of sheet, over which the planar elongational flow dominates, is possibly replaced by one-dimensional analytical solution. Then only a partial flow domain near the edge region of sheet, where the flow kinematics cannot be described by the planar elongational flow itself, requires three-dimensional numerical simulation. Good agreement is observed between the coupled approach developed in this study and the full three-dimensional numerical simulation previously developed and reported by the authors. This coupled approach may have provided flexibility with low costs to accommodate a wide range of die sizes in sheet casting process.

• Nuclear Engineering and Technology
• /
• v.16 no.1
• /
• pp.21-28
• /
• 1984

The reduction of the truncation error including numerical diffusion, has been one of the most important tasks in the development of numerical methods. The stream line method is used to cancel cross numerical diffusion and some of the non-diffusion type truncation error. The two-step stream line method which is the combination of the stream line method and finite difference methods is developed in this work for the solution of the govern ing equations of incompressible buoyant turbulent flow. This method is compared with the finite difference method. The predictions of both classes of numerical methods are compared with experimental findings. Truncation error analysis also has been performed in order to the compare truncation error of the stream line method with that of finite difference methods.

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.13-22
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second-order scheme.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.11 no.2
• /
• pp.1-6
• /
• 1974

The three-dimensional correction factor of the added mass of finite-length elliptic cylinders in vertical vibration in a free surface was calculated. This problem has already been dealt by T. Kumai[5] to contribute to analytical prediction of the three-dimensional correction factor for the added mass in vertical vibration of ships. In Kumai's work, the body boundary condition involved in the appropriate boundary value problem was approximately treated in the course of obtaining the solution. In this work, obtaining the solution derived from mathematically exact treatment of the body boundary condition, the author recalculated the three-dimensional correction factor for length-beam ratio $4{\sim}8$, beam-draught ratio $2.00{\sim}4.50$ and number of nodes from 2 to 7. And the numerical results were compared with both Kumai's results and the author's experimental data for two and three-noded vibrations of the cylinder of beam-draught ratio 2.40 The comparison of the numerical results shows that the author's are always higher than the Kumai's as expected. And the comparison of the numerical results with experimental data shows that the Kumai's numerical results have less deviation in case of two-noded vibration, and that, in case of three-noded vibration, the author's numerical results are in fairly good correspondence.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.4
• /
• pp.58-65
• /
• 1994

It is very important to understand characteristics of solution due to the variation of computational grid sizes, especially when turbulence model not incorporating wall-function is used. The present paper performs numerical investigation on the grid dependency of numerical solution for three dimensional turbulent flow field around a ship. In the present study a finite volume method with a modified sub-grid scale turbulence model and a numerically constructed non-orthogonal curvilinear coordinate system capable of conforming complex ship geometries are used. Numerical studies are then performed for a mathematical Wigley hull and the Series 60, $C_B=0.8$ hull forms. The results for various grid sizes are compared with each other and with measured data to show grid dependencies of numerical solutions.

• Geomechanics and Engineering
• /
• v.37 no.1
• /
• pp.21-29
• /
• 2024

The effective management of damage in tunnels is crucial for ensuring their safety, longevity, and operational efficiency. In this paper, we propose an educational management model tailored specifically for addressing damage in tunnels, utilizing numerical solution techniques. By leveraging advanced computational methods, we aim to develop a comprehensive understanding of the factors contributing to tunnel damage and to establish proactive measures for mitigation and repair. The proposed model integrates principles of tunnel engineering, structural mechanics, and numerical analysis to facilitate a systematic approach to damage assessment, diagnosis, and management. Through the application of numerical solution techniques, such as finite element analysis, we demonstrate the efficacy of the proposed model in simulating various damage scenarios and predicting their impact on tunnel performance. Additionally, the educational component of the model provides valuable insights and training opportunities for tunnel management personnel, empowering them to make informed decisions and implement effective strategies for ensuring the structural integrity and safety of tunnel infrastructure. Overall, the proposed educational management model represents a significant advancement in tunnel management practices, offering a proactive and knowledge-driven approach to addressing damage and enhancing the resilience of tunnel systems.

• Water Engineering Research
• /
• v.2 no.2
• /
• pp.75-87
• /
• 2001

A new moving boundary technique for leap-frog finite difference numerical mode is proposed for the resonable simulation of coastal inundation over discontinuous topography. The new scheme improves the moving boundary technique developed by Imamura(1996). The present scheme is tested using the analytical solution of Thacker(1981) for the case of free oscillation with moving boundary in a parabolic bowl. Finally, a numerical simulation is conducted for the flooding over a tidal barrier constructed on a simple concave geometry. A general feature of inundation over a discontinuous topography is well described by the numerical model.

• Structural Engineering and Mechanics
• /
• v.18 no.1
• /
• pp.59-78
• /
• 2004

This paper is devoted to illustrate some numerical procedures to solve the boundary equilibrium problems of three-dimensional solids that are subjected to thermal strains. The constitutive equations take into account the bounded tensile strength of the material and they are presented in the framework of non-linear elasticity and small strains. The associated equilibrium problem is solved numerically by means of the finite element method and the numerical techniques, i.e. the Newton-Raphson method and the secant method, are revised in order to assure the solution convergence of the discretized problem. Some numerical examples are illustrated.

• Journal of applied mathematics & informatics
• /
• v.24 no.1_2
• /
• pp.17-30
• /
• 2007

A new scheme for solving the Vlasov equation using a compactly supported wavelets basis is proposed. We use a numerical method which minimizes the numerical diffusion and conserves a reasonable time computing cost. So we introduce a representation in a compactly supported wavelet of the derivative operator. This method makes easy and simple the computation of the coefficients of the matrix representing the operator. This allows us to solve the two equations which result from the splitting technique of the main Vlasov equation. Some numerical results are exposed using different numbers of wavelets.