• Nuclear Engineering and Technology
• /
• v.54 no.12
• /
• pp.4522-4533
• /
• 2022

For severe accident assessment of reactor pressure vessel (RPV), it is important to develop an accurate model that can predict transient thermo-mechanical behavior of the RPV lower head under the given condition. The present study revisits the lower head failure with two- and three-dimensional finite element models. In particular, we aim to give clear insight regarding the effect of the three-dimensionality present in the distribution of the thickness and thermal load of the lower head. For a rigorous validation of the result, both the OLHF-1 and the OLHF-2 tests are considered in this study. The result suggests that the three-dimensional effect is not negligible as far as the failure location is concerned. The non-uniformity of the thickness distribution is found to affect the failure location and time. The thermal load, which may not be axisymmetric in general, has the most significant effect on the failure assessment. We also observe that the creep property can affect the global deformation of the lower head, depending on the applied mechanical load.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.2
• /
• pp.204-213
• /
• 2000

This study deals with jet impingement, which is extensively used in the process industries to achieve intense heating, cooling or drying rates and also widely employed as a test flow for turbulent models due to its complex flow configuration, on a flat plate by numerical methods. In this calculation, the finite volume method was employed to solve the Navier-stokes equation based on the non-orthogonal coordinate with non-staggered variable arrangement. To get a better understanding for the fluid flow and heat transfer characteristics of the turbulent jet impingements, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model was adapted and compared with the experimental data and the result of standard $k-{\varepsilon}$ turbulent model. Numerical calculations were carried out with various flow rates, nozzle to plate distances. In the case of the axisymmetric jet impingement on a flat plate, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model showed better agreement with the experimental data than the standard $k-{\varepsilon}$ turbulent model in the prediction of the mean velocity profiles, the turbulent velocity profiles. the turbulent shear stress and the heat transfer rate. The highest heat transfer rate can be obtained when the impingement occurs within the potential core..

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1149-1164
• /
• 1997

In this study, the applicability of the RNG k-.epsilon. model to the analysis of the complex flows is studied. The governing equations based on a non-orthogonal coordinate formulation with Cartesian velocity components are used and discretized by the finite volume method with non-staggered variable arrangements. The predicted results using the RNG k-.epsilon. model of three complex flows, i.e., the flow over a backward-facing step and a blunt flat plate, the flow around a 2D model car are compared to these from the standard k-.epsilon. model and experimental data. That of the unsteady axisymmetric turbulent flow within a cylinder of reciprocating model engine including port/valve assembly and the spray characteristics within a chamber of direct injection model engine are compared to these from the standard k-.epsilon. model and experimental data. The results of reattachment length, separated eddy size, average surface pressure distribution using the RNG k-.epsilon. model show more reasonable trends comparing with the experimental data than those using the modified k-.epsilon. model. Although the predicted rms velocity using the modified k-.epsilon. model is lower considerably than the experimental data in incylinder flow with poppet valve, predicted axial and radial velocity distributions at the valve exit and in-cylinder region show good agreements with the experimental data. The spray tip penetration predicted using the RNG k-.epsilon. model is more close to the experimental data than that using the modified k-.epsilon. model. The application of the RNG k-.epsilon. model seems to have some potential for the simulations of the unsteady turbulent flow within a port/valve-cylinder assembly and the spray characteristics over the modified k-.epsilon. model.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.12 no.6
• /
• pp.142-151
• /
• 2013

This study develops an automation system for metallic mold design that is applicable in forging non-axial symmetric parts. The metallic mold design program is used to design the metallic mold using two-dimensional axial symmetric metallic molds and to predict the stress concentration using finite element analyses. Then, the program redesigns the metallic mold using variables such as the optimal split diameter, maximum allowable inner pressure, fit tolerance, and stress distribution, which are calculated using the metallic mold design program. When the involute spur gear is forged, stress concentration occurs on the tooth root bounded at the symmetric surface. The SCM4 material is suitable for metallic molds because the stress is less than the yield strength of the insert and it acts on the tooth root regardless of the inner pressure. The metallic mold for forging non-axial symmetric parts can be designed through adjusting the magnitude of the contact pressure. The program developed in this study can be applied to metallic mold designs in involute spur gears of forging, which is an ordinary non-axial symmetric part.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2862-2867
• /
• 2007

The steady non-reacted compressible flow field in a symmetric micro-thruster, which is used for the accurate attitude control of a satellite, is analyzed varying the nozzle pressure ratio (NPR) to investigate the plume characteristics. The nozzle throat diameter is 0.06 inch and the area ratio is 56. The recirculation region is found just behind the normal shock at the several NPRs due to the locally adverse pressure gradient along the nozzle centerline when the environmental pressure is atmospheric. This phenomenon, the cause of flow loss, is similar to the flow behind a blunt body. As NPR increases the location of Mach disk, characteristics of the normal shock, moves downstream and its strength increases. The Mach number distribution appears in a wave-type patter after the normal shock because oblique shocks are reflected on the shock boundaries especially when NPRs are very high.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.19 no.9
• /
• pp.93-99
• /
• 2020

Recently, as the weight reduction of vehicles has been actively progressed, parts developed using aluminum 60XX series from existing steel materials are increasing. In this paper, the bushing used for the front frame rail, which is one of the parts for fixing engines and other parts in automobiles, was changed to an aluminum material of the Al60XX series, and it was intended to be produced by applying of cold forging method. The bushing is a part that secures the engine frame, and in order to produce it by cold forging, the molding limit is predicted through process design, and a multi-stage process is designed through finite element analysis. In addition, in order to verify the feasibility of the designed forging process, the limits of the multi-step process were verified based on the Cockcroft Latham theory, and the crack and overlap of the actual forging work were predicted and improved.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.11a
• /
• pp.1476-1481
• /
• 2007

Numerical investigation on radiation characteristics of discrete frequency noise from asymmetry aero-intakes was carried out. The near-field predictions were obtained by solving the linearized Euler equations with computational aeroacoustic techniques consisting of high order finite difference scheme, non-reflecting boundary conditions, oversetgrid techniques. For the prediction of far-field directivity pattern, the Kirchhoff integral method was applied. By comparing the directivities of noise radiating from the scarf and the scoop aero-intakes with that from an axisymmetric aero-intake, it is shown that noise reduction at downward peak radiation angle can be achieved. The scattering of the radiating acoustic wave by background mean flow shifts the peak lobe radiation angle toward ground and increases the amplitude of the acoustic pressure compared with the cases without mean flow effect.

• Korea-Australia Rheology Journal
• /
• v.18 no.2
• /
• pp.109-118
• /
• 2006

A finite-element method was employed to analyze axisymmetric unsteady motion of a deformable bubble near the wall. In the present study a deformable bubble in a Newtonian medium near the wall was considered. In solving the governing equations a structured mesh generator was used to describe the collapse of highly deformed bubbles with the Arbitrary Lagrangian Eulerian (ALE) method being employed in order to capture the transient bubble boundary effectively. In order to check the accuracy of the present FE analysis we compared the results of our FE solutions with the result of the collapse of spherical bubbles in a large body of fluid in which solutions can be obtained using a 1D FE analysis. It has been found that 1D and 2D bubble deformations are in good agreement for spherically symmetric problems confirming the validity of the numerical code. Non-spherically symmetric problems were also solved for the collapse of bubble located near a plane solid wall. We have shown that a microjet develops at the bubble boundary away from the wall as already observed experimentally. We have discussed the effect of Reynolds number and distance of the bubble center from the wall on the transient collapse pattern of bubble.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.12
• /
• pp.3287-3296
• /
• 1994

For the numerical computation of three-dimensional compressible flow field within a blade row in a centrifugal compressor, a quasi 3-dimensional solver which combines a reversible B-B flow and an irreversible H-S flow using finite element methods was developed. In a reversible B-B flow, the governing coordinates are modified in order to be applied to any type of turbomachinery, and two kinds of stream functions are introduced in order to make the Kutta condition exactly satisfied. In an irreversible H-S flow, the changes of entropy in the irreversible governing equations are determined not by empirical source but by the theoretical treatment of dissipation forces. The dissipation forces are obtained from the distribution of shear stresses in the flow passage which are given from the wall shear stresses using the exponential functions. A more accurate quasi-3-dimensional solver is established where the effect of body forces is involved in the non-axisymmetric H-S flow. Some numerical results obtained from authors' previous studies for axial flow machines assure that the present method is able to predict well as long as the flow is subsonic and not under strong viscous effect.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1368-1379
• /
• 2000

A numerical study on swirling pulverized coal combustion in an axisymmetric enclosure is carried out by applying the 2-phase weighted sum of gray gases model (WSGGM) approach with the discrete ordinate method (DOM) to model the radiative heat transfer equation. In the radiative transfer equation, the same polynomial equation and coefficients for weighting factors as those for gas are adopted for the coal/char particles as a function of partial pressure and particle temperature. The Eulerian balance equations for mass, momentum, energy, and species mass fractions are adopted with the standard and RNG k-${\varepsilon}$ turbulence model, whereas the Lagrangian approach is used for the particulate phase. The eddy-dissipation model is employed for the reaction rate for gaseous mixture, and the single-step and two-step first-order reaction model for the devolatilization process for coal. Special attention is given to establish the thermal boundary conditions on radiative transfer equation By comparing the numerical results with experimental ones, the radiation model used here is confirmed and found to provide an alternative for simulating the radiative transfer.