• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1213-1222
• /
• 1999

A finite element method is developed for calculating the temperature and enthalpy distribution and accordingly the solid, liquid and mushy zone in a three-dimensional body subjected to any heat boundary conditions. The method concurrently consider both temperature and enthalpy for consideration of the latent heat effect, differently from other methods of using a special energy balance equation for solving a mushy zone. The developed brick element has eight nodes with one degree of freedom at each node. The numerical method and procedure are verified using the results of one and two dimensional analytic solutions and by other researchers. It is shown that the present method presents a consistent and stable results in either abrupt or ranged phase change problems. Moreover, the numerical results by the present method are hardly effected by the calculation time steps which otherwise are difficult to determine in most phase change problems. Finally, as a three-dimensional application, a T-shaped body of a phase change is presented and the temperature and enthalpy variation along the time are solved.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.4
• /
• pp.433-440
• /
• 2004

Since the numerical analysis was adopted in the mold design, lots of computational methods have been proposed for the simulations of casting processes for the various shaped molds. Today, it is possible to simulate the filling and solidification processes of most casts using the VOF technique. Though the three-dimensional numerical model based on the Cartesian coordinate system can be applied to any shape of cast, it becomes very inefficient when the three-dimensional model is applied to the cast of axi-symmetrical shape since the control volume includes at least 11 of the physical model. In addition, the more meshes should be distributed along the circumferential boundaries of curved shape in the Cartesian coordinate system fur the better results, while such curved circumferential boundary does not need to be considered in the two-dimensional cylindrical coordinate system. This motivates the present study i.e. developing a two-dimensional numerical model for the axi-symmetrically shaped casts. The SIMPLER algorithm, the VOF method, and the equivalent specific heat method have been adopted in the combined algorithm for the flow calculation, the free surface tracking, and the phase change heat transfer, respectively. The numerical model has been applied to the casting process of a pulley, and it was proven that the mesh and time effective calculation was accomplished comparing to the calculation using three-dimensional model.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.04a
• /
• pp.195-198
• /
• 2000

The process of squeeze casting for metal matrix composites (MMCs) has been simulated numerically by using finite difference method. The governing equations to describe fluid flow through porous medium and heat transfer are applied to two dimensional model which is similar to a real system. A computational code has been developed to solve this problem. The influence on infiltration kinetics and solidification time of several parameters is investigated. Cooling curves and temperature distribution with time and position is also shown. The result can be used to design the squeeze casting for MMCs.

• Journal of Korea Foundry Society
• /
• v.11 no.1
• /
• pp.79-84
• /
• 1991

The process of solidification of metal is accompanied by liquid-solid change and known as Stefan's heat conduction problem on the moving boundary. In this study the temperature distribution in ingot and metallic mould during casting was analyzed by the two dimensional heat conduction theory. The transient temperature distribution was numerically calculated using a finite element method on the nodal point of mesh screen representing ingot and mould cross section. The theory was applied on the casting of aluminum(purity ; 99%) in flat ingot mould of GC25. The analysis will make it possible to calculate an optimum mould shape of which temperature gradient becomes minimum.

• Journal of Welding and Joining
• /
• v.16 no.5
• /
• pp.56-66
• /
• 1998

During solidification or welding of alloys, the solute redistribution brings out microsegregation. The microsegregation causes the formation of non-equilibrium second phases, shrinkage and porosity degrading mechanical/chemical properties Therefore, it has been required to predict microsegregation quantitatively. To predict the degree of microsegregation, more exact and appropriate computer simulation technique has been actively used during last two decades. To predict the degree of microsegregation in weld metal, an advanced two dimensional model was suggested. In the new model, both primary and secondary arm regions were defined for the analysis region. The growth in the primary arm regina was assumed to be a planar for effective calculation. Especially, for the growth of a secondary arm, a simple and effective mathematical function was established to show the growing pattern, the solute diffusion in the solid phase was calculated by finite difference method (FDM). The solid-liquid interface movement was considered to be in local equilibrium state. The experiments for welding of 310S stainless steel were carried out in order to examined the reasonability and feasibility of this model. The concentration profiles of the solute predicted by this model were compared with those obtained from experimental works.

• Nuclear Engineering and Technology
• /
• v.54 no.10
• /
• pp.3909-3918
• /
• 2022

In the event of a severe accident, the reactor core may melt due to insufficient cooling. the high-temperature core melt will have a strong interaction (FCI) with the coolant, which may lead to steam explosion. Steam explosion would pose a serious threat to the safety of the reactors. Therefore, the study of steam explosion is of great significance to the assessment of severe accidents in nuclear reactors. This research focuses on the development of a two-dimensional steam explosion integrated analysis code called SEINA. Based on the semi-implicit Euler scheme, the three-phase field was considered in this code. Besides, the influence of evaporation drag of melt and the influence of solidified shell during the process of melt droplet fragmentation were also considered. The code was simulated and validated by FARO L-14 and KROTOS KS-2 experiments. The calculation results of SEINA code are in good agreement with the experimental results, and the results show that if the effects of evaporation drag and melt solidification shell are considered, the FCI process can be described more accurately. Therefore, it is proved that SEINA has the potential to be a powerful and effective tool for the analysis of steam explosions in nuclear reactors.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.14 no.5
• /
• pp.897-903
• /
• 2011

The objective of this study is to investigate the thermal behavior on material surface and the variation in the surface morphology of aluminum 6061 alloy by the Nd:YAG pulsed-laser irradiation. First, we predicted the surface temperature variation during pulsed-laser irradiation by using the two dimensional finite element analysis. When the pulsed-laser of 133 mJ energy and 5 ns pulse duration is irradiated on the surface of aluminum alloy, the material surface is thought to be melting because the surface temperature rises steadily up to about $660^{\circ}C$ exceeding the melting point. Also, the experimental results show that the solidification microstructure has been developed clearly after surface melting. Second, the diameter of melted zone was analysed by finite element analysis and measured by OM(Optical Microscopy). It increased logarithmically with increase in the number of laser irradiation. In addition, AFM(Atomic Force Microscopy) measurement showed an increase in the average surface roughness during pulsed-laser irradiation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.1
• /
• pp.142-155
• /
• 1992

Numerical analysis is performed for the unsteady axisymmetric two dimensional phase change problem of freezing of water around a vertical tube. Heat conduction in the tube wall and solid phase, natural convection in liquid phase and volume expansion caused by density difference between solid and liquid phases are included in the numerical analysis. Existing correlation is used for estimating density-temperature relation of water, and the effect of volume expansion is reflected as fluid velocity at the interface and the free surface. As pure water has maximum density at 4.deg. C, it is found that there exists an initial temperature at which the flow direction reverses near the interface and by this effect the slope of interface becomes reversed depending on the initial temperature of water. By considering natural convection and solid-liquid density difference in the calculation, their effects on phase change process are studied and the effects of various parameters are also studied quantitatively.