• Journal of the Korean Society for Precision Engineering
• /
• v.18 no.12
• /
• pp.60-65
• /
• 2001

To simulate the real molding conditions, the effects of phase change and compressibility of the resin were considered in the present investigation. A modified Cross model with either an Arrhenius-type or WLF-type functional form was used for modeling viscosity of the resin. A double-domain Tait equation of state was employed to describe the compressibility of the resin during molding. The energy balance equation including latent-heat dissipation fur semi-crystalline materials was solved in order to predict the solidified layer and temperature profile. Injection molding experiments were carried out using polypropylene(PP) in the present study. Based on the comparison between experiments and simulations, it was found out the predicted pressure distributions and melt front propagations were accurate. Thus it was concluded that the program developed in this study was proved to be useful in simulations of injection molding process.

• KIEAE Journal
• /
• v.16 no.1
• /
• pp.81-87
• /
• 2016

Purpose: Heat transfers phenomena are described by the second order partial differential equation and its boundary conditions. In a three-dimensional structure of a building, the heat transfer phenomena generally include more than one material, and thus, become complicate. The analytic solutions are useful to understand heat transfer phenomena, but they can hardly be applied in engineering or design problems. Engineers and designers have generally been forced to use numerical methods providing reliable results. Finite volume methods with the unstructured grid system is only the suitable means of the analysis for the complex and arbitrary domains. Method: To obtain an numerical solution, a discretization method, which approximates the differential equations, and the interpolation methods for temperature and heat flux between two or more materials are required. The discretization methods are applied to small domains in space and time, and these numerical solutions form the descretized equations provide approximated solutions in both space and time. The accuracy of numerical solutions is dependent on the quality of discretizations and size of cells used. The higher accuracy, the higher numerical resources are required. The balance between the accuracy and difficulty of the numerical methods is critical for the success of the numerical analysis. A simple and easy interpolation methods among multiple materials are developed. The linear equations are solved with the BiCGSTAB being a effective matrix solver. Result: This study provides an overview of discretization methods, boundary interface, and matrix solver for the 3-dimensional numerical heat transfer including two materials.

• KIEAE Journal
• /
• v.16 no.6
• /
• pp.159-166
• /
• 2016

Purpose: The architectural structures in the engineering field include more than one material, and the heat transfer through these multiple materials becomes complicated. More or less, the analytic solutions obtained by the hand calculation can provide the limited information of heat transfer phenomena. However, the engineers have generally been forced to obtain reliable results than those of the hand calculation. The numerical calculation such as a finite volume methods with the unstructured grid system is only the suitable means of the analysis for the complex and arbitrary domains that consists of multiple materials. In this study, a new numerical code is developed to provide temperature distributions in the multiple material domains, and the results of this code are compared with the validation cases in ISO10211. Method: Finite volume methods with an unstructured grid is employed. In terms of numerical methods, the heat transfer conduction coefficient is not defined on the surface of the cell between different material cells. The heat transfer coefficient is properly defined to accurately mimic the heat transfer through the surface. The boundary conditions of heat flux considering radiation or heat convection are also developed. Result: The comparison between numerical results and ISO 10211 cases. We are confirmed that the numerical method provides the proper temperature distributions, and the heat transfer equation and its boundary conditions are developed properly.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.7
• /
• pp.486-494
• /
• 2009

One-dimensional modeling of $CH_4$/air premixed flame was conducted to validate the heat loss model and investigate NOx formation characteristics in the postflame region. The predicted temperature and NO concentration were compared to experimental data and previous heat loss model results using a constant gradient of temperature (100 K/cm). The following conclusions were drawn. In the heat loss model using steady-state heat transfer equation, the numerical results using the effective heat loss coefficient ($h_{eff}$) of $1.0\;W/m^2K$ were in very good agreement with the experiments in terms of temperature and NO concentration. On the other hand, the calculated values using the constant gradient of temperature (100 K/cm) were lower than that in the experiments. Although the effects of heat loss suppress NO production near the flame region, a significant difference in NO concentration was not found compared to that under adiabatic conditions. In the postflame region, however, there were considerable differences in NO emission index as well as the contribution of NO formation mechanisms. In particular, in the range of ${\phi}\;{\geq}\;0.8$, the prompt NO mechanism plays an important role in the NO reduction under the adiabatic condition. On the other hand, the mechanism contributes to the NO production under the heat loss conditions.

• Fire Science and Engineering
• /
• v.11 no.1
• /
• pp.21-35
• /
• 1997

The natural convection and combined heat transfer induced by fire in a rectangular enclosure is numerically studied. The model for this numerical analysis is partially opened right wall. The solution procedure includes the standard k-$\varepsilon$ model for turbulent flow and the discrete ordinates method (DOM) is used for the calculation of radiative heat transfer equation. In numerical study, SIMPLE algorithm is applied for fluid flow analysis, and the investigations of combustion gas induced by fire is performed by FAST model of HAZARD I program. In this study, numerical simulation on the combined naturnal convection and radiation is carried out in a partial enclosure filled with absorbed-emitted gray media, but is not considered scattering problem. The streamlines, isothermal lines, average radiation intensity and kinetic energy are compared the results of pure convection with those of the combined convection-radiation, the combined heat transfer. Comparing the results of pure convection with those of the combined convection-radiation, the combined heat transfer analysis shows the stronger circulation than those of the pure convection. Three different locations of heat source are considered to observe the effect of heat source location on the heat transfer phenomena. As the results, the circulation and the heat transfer in the left region from heating block are much more influenced than those in the right region. It is also founded that the radiation effect cannot be neglected in analyzing the building in fire. And as the results of combustion gas analysis from FAST model, it is found that O2 concentration is decreased according to time. While CO and CO2 concentration are rapidly increased in the beginning(about 100sec), but slowly decreased from that time on.

• Journal of Advanced Marine Engineering and Technology
• /
• v.30 no.5
• /
• pp.572-579
• /
• 2006

Accurate heat release analysis based on the cylinder pressure trace is important for evaluating combustion process of diesel engines. However, traditional single-zone heat release models (SZM) have significant limitations due mainly to their simplified assumptions of uniform charge and homogeneity while neglecting local temperature distribution inside cylinder during combustion process. In this study, a heat release analysis based on single-zone model has been evaluated by comparison with computational simulation result using Fire-code, which is based on multidimensional model (MDM). The limitations of the single-zone assumption have been estimated, To overcome these limitations, an improved model that includes the effects of spatial non-uniformity has been applied. From this improved single-zone heat release model (Improved-SZM), two effective values of specific heat ratios, denoted by ${\gamma}_V$ and ${\gamma}_H$ in this study, have been introduced. These values are formulated as the function of charge temperature changing rate and overall equivalence ratio. Also, it is applied that each equation of ${\gamma}_V$ and ${\gamma}_H$ has respectively different slopes according to several meaningful periods during combustion progress. The heat release analysis results based on improved single-zone model gives a good agreement with FIRE-code results over the whole range of operating conditions of target engine, Hyundai HiMSEN H21/32.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.14 no.11
• /
• pp.5317-5322
• /
• 2013

This study presents a prediction method for heat transfer performance of a gas cooler of $CO_2$ heat pump using ${\epsilon}$-NTU method, and compared the results with the experimental data from the open literature. The heat transfer rate, refrigerant side outlet temperature and water side outlet temperature were calculated by using EES(Engineering Equation Solver)program in multi-tube-in-tube type $CO_2$ heat pump gas cooler. Analysis was performed in two methods : The first method performed without dividing into the test section by applying an analysis of the mean properties(mean analysis). The second method, tube length divided into 50 sections, was applied to the local properties(local analysis). From the present study, a good agreement at the local analysis was obtained between the analytical and experimental results by 0.3~1.1%, 1.31~1.88% and 3.12~5.18% for heat transfer rate, water and refrigerant side outlet temperatures, respectively.

• Fire Science and Engineering
• /
• v.30 no.1
• /
• pp.49-56
• /
• 2016

This paper numerically examines, straight and twisted electrical heat trace installations for their anti-freezing effects on water inside a pipe. The unsteady incompressible Navier-Stokes equations coupled with an energy equation were solved to compare the two installation methods. The heat conduction of the pipe with a heat source interacts with the natural convection of the water, and the conjugate heat transfer was considered using a commercial code (ANSYS-FLUENT) based on a SIMPLE-type algorithm. Numerical experiments, were done to investigate the isotherms and the vector fields in the water region to extract the evolutions of the minimum and maximum temperatures of the water inside the pipe. There was no substantial difference in the anti-freezing effects between the straight and twisted. Therefore, the straight installation is recommended after considering the damage and short circuit behavior of the electrical heat trace.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2696-2701
• /
• 2007

For the analysis of a two-phase flow, the interaction between two phases such as the interfacial momentum or heat transfer is proportional to the interfacial area. So the interfacial area concentration (IAC) is one of the most important parameters governing the behavior of each phase. This study focuses on the development of a computational fluid dynamics (CFD) code for investigating a boiling flow with a one-group IAC transport equation. It was based on the two-fluid model and governing equations were calculated by SMAC algorithm. For checking the robustness of the developed code, the experiment of a subcooled boiling in a vertical annulus channel was analyzed to validate the capability of the IAC transport equation. As the results, the developed code was confirmed to have the capability in predicting multi-dimensional phenomena of vapor generation and propagation in a subcooled boiling.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.7-13
• /
• 2011

A comparative analysis of thermal models in the lattice Boltzmann method(LBM) for the simulation of laminar natural convection in a square cavity is presented. A HYBRID method, in which the thermal equation is solved by the Navier-Stokes equation method while the mass and momentum conservation are resolved by the lattice Boltzmann method, is introduced and its merits are explained. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with a deferred correction method to ensure stability of the solutions. The HYBRID method and the double-population method are applied to the simulation of natural convection in a square cavity and the predicted results are compared with the benchmark solutions given in the literatures. The predicted results are also compared with those by the conventional Navier-Stokes equation method. In general, the present HYBRID method is as accurate as the Navier-Stokes equation method and the double-population method. The HYBRID method shows better convergence and stability than the double-population method. These observations indicate that this HYBRID method is an efficient and economic method for the simulation of incompressible fluid flow and heat transfer problem with the LBM.