• Journal of the Korean Society of Propulsion Engineers
• /
• v.8 no.4
• /
• pp.109-121
• /
• 2004

A distributed nonlinear mathematical model for investigation of regenerative heat exchangers of both a continuous and periodic operation is described in the paper. The non-iterative numerical integration scheme for conjugate unsteady heat exchange problem of one dimensional flows and two dimensional matrix wall conductivity is developed. Case study of a regenerative heat exchanger with a rotary ceramic matrix is presented. The range of optimum rotation rates of the regenerator providing the greatest calorific efficiency is determined.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.10
• /
• pp.2104-2113
• /
• 2002

A finite element model is developed for the process of squeeze casting of metal matrix composites (MMCs) in cylindrical molds. The fluid flow and the heat transit. are fundamental phenomena in squeeze casting. To describe heat transfer in the solidification of molten aluminum, the energy equation is written in terms of temperature and enthalpy are applied in an axisymmetric model which is similar to the experimental system. A one dimensional flow model simulates the transient metal flow. A direct iteration technique was used to solve the resulting nonlinear algebraic equations, using a computer program to calculate the enthalpy, temperature and fluid velocity. The cooling curves and temperature distribution during infiltration and solidification were calculated fer pure aluminum. Experimentally, the temperature was measured and recorded using thermocouple wire. The measured time-temperature data were compared with the calculated cooling curves. The resulting agreement shows that the finite element model can accurately estimate the solidification time and predict the cooling process.

• Journal of Welding and Joining
• /
• v.21 no.5
• /
• pp.525-533
• /
• 2003

The dimensional accuracy and global roughness between successive layers of VLM-s, which is a new rapid prototyping process using hotwire cutter and EPS foam, depend significantly on the operating parameters of hotwire cutter. In the present study, the effect of cutting angle on the kerf width and edge shape in hotwire cutting of EPS foam for the case of single-sloped cutting with one cutting angle was investigated. Through single-sloped cutting tests, the modified relationship between kerf width and effective heat input, considering the effect of the cutting angle, and the relationship between the melted area and the cutting angle were obtained. In order to investigate the effect of cutting angles on the thermal field in EPS foam, transient heat transfer analyses using single-sloped volumetric heat flux model and locally-conformed mesh were performed. Through the comparison between experimental and numerical results, it was shown that the proposed analysis model is needed to estimate the three-dimensional temperature distribution of the EPS foam for the case of single-sloped hotwire cutting.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.1
• /
• pp.49-56
• /
• 1997

Heat-transfer enhancement is seeked through modifications of fin surface. Real life plate-fin heat exchangers have complex three-dimensional geometries. Fins can have arrays of dimples and are attached to rows of penetrating tubes. To isolate the effect of surface modification, we model the real flow by a two-dimensional channel flow with a dimple on one side. The flow is analysed by solving the incompressible Navier-Stokes equation by a finite volume method on a generalized boundary-fitted coordinate. Results show a trapped vortex inside the dimple for all cases computed. Local maximum of Nusselt number occurs near the downstream end of the dimple, due to such a vortex. Location of the vortex does not change with respect to the wall temperature change, but moved downstream when Reynolds number increases. This, together with the results that in all cases vortex core is somewhat downstream of the dimple center, suggests that the mean flow above continuously feeds the kinetic energy to the recirculating flow. Heat transfer enhancement and pressure losses are studied through analysing the relevant dimensionless parameters like, Nusselt number and friction factor. In all cases computed, dimpled channel flow experiences less pressure loss than two-dimensional Poiseuille flow.

• Korean Chemical Engineering Research
• /
• v.46 no.3
• /
• pp.604-609
• /
• 2008

This study is about the modeling and simulation of heat transfer in the box for packaging and shipping of vaccines. Comparison of the simulation results with experimental data revealed that a one-dimensional model (a spherical model of using a radius equivalent to the rectangular geometry of box) showed good agreement with experimental data during cooling process but did not successfully simulate heating process. It is considered that a rigorous boundary condition is not properly applied for outer surface of the box. However, we could successfully develop a basic algorithm for simulating heat transfer through multi-slabs combined with different materials including phase change material.

• Journal of Mechanical Science and Technology
• /
• v.20 no.9
• /
• pp.1514-1524
• /
• 2006

This work presents the three-dimensional analysis of flow and heat transfer performed for a wire-wrapped fuel assembly of liquid metal reactor using Reynolds-averaged Wavier-Stokes analysis in conjunction with 557 model as a turbulence closure. The whole fuel assembly has been analyzed for one period of the wire-spacer using periodic boundary conditions at inlet and outlet of the calculation domain. Three different assemblies, two 7-pin wire-spacer fuel assemblies and one bare rod bundle, apart from the pressure drop calculations for a 19-pin case, have been analyzed. Individual as well as a comparative analysis of the flow field and heat transfer have been discussed. Also, discussed is the position of hot spots observed in the wire-spacer fuel assembly. The flow field in the subchannels of a bare rod bundle and a wire-spacer fuel assembly is found to be different. A directional temperature gradient is found to exist in the subchannels of a wire-spacer fuel assembly Local Nusselt number in the subchannels of wire-spacer fuel assemblies is found to vary according to the wire-wrap position while in case of bare rod bundle, it's found to be constant.

• Nuclear Engineering and Technology
• /
• v.53 no.3
• /
• pp.833-841
• /
• 2021

During the postulated severe accident of nuclear reactor, eutectic reaction leads to low-temperature melting of fuel cladding and early failure of core structure. In order to model eutectic melting with the moving particle semi-implicit (MPS) method, the eutectic reaction model is developed to simulate the eutectic reaction phenomenon. The coupling of mass diffusion and phase diagram is applied to calculate the eutectic reaction with the uniform temperature. A heat transfer formula is proposed based on the phase diagram to handle the heat release or absorption during the process of eutectic reaction, and it can combine with mass diffusion and phase diagram to describe the eutectic reaction with temperature variation. The heat transfer formula is verified by the one-dimensional melting simulations and the predicted interface position agrees well with the theoretical solution. In order to verify the eutectic reaction models, the eutectic reaction of uranium and iron in two semi-infinite domains is simulated, and the profile of solid thickness decrease over time follows the parabolic law. The modified MPS method is applied to calculate Transient Reactor Test Facility (TREAT) experiment, the penetration rate in the simulations are agreeable with the experiment results. In addition, a hypothetical case based on the TREAT experiment is also conducted to validate the eutectic reaction with temperature variation, the results present continuity with the simulations of TREAT experiment. Thus the improved method is proved to be capable of simulating the eutectic reaction in the severe accident.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.5
• /
• pp.359-366
• /
• 2008

A two-dimensional thermal model of proton exchange membrane fuel cell with large active area is developed to investigate the performance of fuel cell with large active area over various thermal management conditions. The core sub-models of the two-dimensional thermal model are one-dimensional agglomerate structure electrochemical reaction model, one-dimensional water transport model, and a two-dimensional heat transfer model. Prior to carrying out the simulation, this study is contributed to set up the operating temperature of the fuel cell with large active area which is a maximum temperature inside the fuel cell considering durability of membrane electrolyte. The simulation results show that the operating temperature of the fuel cell and temperature distribution inside the fuel cell can affect significantly the total net power at extreme conditions. Results also show that the parasitic losses of balance of plant component should be precisely controlled to produce the maximum system power with minimum parasitic loss of thermal management system.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.20-26
• /
• 2008

A numerical study is carried out to analyze the steady three-dimensional turbulent flow and convective heat transfer in a staggered pin-fin array with diamond shaped elements at various geometrical configurations. Steady Reynolds-averaged Navier-Stokes equations and energy equation are solved using a finite volume based solver. Shear stress transport (SST) model is used as turbulence closure. The computational domain is composed of one pitch of pin-fin displacement with periodic boundary conditions on the surfaces normal to the streamwise direction and the cross-streamwise direction. The numerical results for Nusselt number and friction factor are validated with experimental results. The effects of pin angle, pin height and pitch on Nusselt number, friction factor and efficiency index are investigated.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.997-1005
• /
• 2005

The extinction characteristics of low strain rate normal gravity (1-g) nonpremixed methane-air flames were studied numerically and experimentally. A time-dependent axisymmetric two-dimensional (2D) model considering buoyancy effects and radiative heat transfer was developed to capture the structure and extinction limits of 1-g flames. One-dimensional (1D) computations were also conducted to provide information on 0-g flames. A 3-step global reaction mechanism was used in both the 1D and 2D computations to predict the measured extinction limit and flame temperature. A specific maximum heat release rate was introduced to quantify the local flame strength and to elucidate the extinction mechanism. Overall fractional contribution by each term in the energy equation to the heat release was evaluated to investigate the multi-dimensional structure and radiative extinction of 1-g flames. Images of flames were taken for comparison with the model calculation undergoing extinction. The two-dimensional numerical model was validated by comparing flame temperature profiles and extinction limits with experiments and ID computation results. The 2D computations yielded insight into the extinction mode and flame structure of 1-g flames. Two combustion regimes depending on the extinction mode were identified. Lateral heat loss effects and multi-dimensional flame structure were also found. At low strain rates of 1-g flame ('Regime A'), the flame is extinguished from the weak outer flame edge, which is attributed to multi-dimensional flame structure and flow field. At high strain rates, ('Regime B'), the flame extinction initiates near the flame centerline due to an increased diluent concentration in reaction zone, which is the same as the extinction mode of 1D flame. These two extinction modes could be clearly explained with the specific maximum heat release rate.