Multi-Objective Optimization of a Microturbine Compact Recuperator

Author(s):  
Diego Micheli ◽  
Valentino Pediroda ◽  
Stefano Pieri

An automatic approach for the multi-objective shape optimization of microgas turbine heat exchangers is presented. According to the concept of multidisciplinary optimization, the methodology integrates a CAD parametric model of the heat transfer surfaces, a three-dimensional meshing tool, and a CFD solver, all managed by a design optimization platform. The repetitive pattern of the surface geometry has been exploited to reduce the computational domain size, and the constant flux boundary conditions have been imposed to better suit the real operative conditions. A new approach that couples cold and warm fluids in a periodic unitary cell is introduced. The effectiveness of the numerical procedure was verified comparing the numerical results with available literature data. The optimization objectives are maximizing the heat transfer rate and minimizing both friction factor and heat transfer surface. The paper presents the results of the optimization of a 50kWMGT recuperator. The design procedure can be effectively extended and applied to any industrial heat exchanger application.

Author(s):  
Diego Micheli ◽  
Valentino Pediroda ◽  
Stefano Pieri

An automatic approach for the multi-objective shape optimization of micro gas turbine heat exchangers is presented. According to the concept of Multi Disciplinary Optimization (MDO), the methodology integrates a CAD parametric model of the heat transfer surfaces, a three dimensional meshing tool and a CFD solver, all managed by a design optimization platform. The repetitive pattern of the surface geometry has been exploited to reduce the computational domain size, and constant flux boundary conditions have been imposed to better suit the real operative conditions. A new approach that couples cold and warm fluids in a periodic unitary cell is introduced. The effectiveness of the numerical procedure was verified comparing the numerical results with available literature data. The optimization objectives are maximizing the heat transfer rate and minimizing both friction factor and heat transfer surface. The paper presents the results of the optimization of a 50 kW MGT recuperator. The design procedure can be effectively extended and applied to any industrial heat exchanger application.


Author(s):  
Takeharu Misawa ◽  
Hiroyuki Yoshida ◽  
Hidesada Tamai ◽  
Kazuyuki Takase

The three-dimensional two-fluid model analysis code ACE-3D is developed in Japan Atomic Energy Agency for the thermal design procedure on two-phase flow thermal-hydraulics of light water-cooled reactors. In order to perform thermal hydraulic analysis of SCWR, ACE-3D is enhanced to supercritical pressure region. As a result, it is confirmed that transient change in subcritical and supercritical pressure region can be simulated smoothly using ACE-3D, that ACE-3D can predict the results of the past heat transfer experiment in the supercritical pressure condition, and that introduction of thermal conductivity effect of the wall restrains fluctuation of wall.


2018 ◽  
Vol 941 ◽  
pp. 2278-2283
Author(s):  
Nima Bohlooli Arkhazloo ◽  
Farzad Bazdidi-Tehrani ◽  
Morin Jean-Benoit ◽  
Mohammad Jahazi

Simulation and analysis of thermal interactions during heat treatment is of great importance for accurate prediction of temperature evolution of work pieces and consequently controlling the final microstructure and mechanical properties of products. In the present study, a three-dimensional CFD model was employed to predict the heating process of large size forged ingots inside an industrial gas-fired heat treatment furnace. One-ninth section of a loaded furnace, including details such as fixing bars and high-momentum cup burners, was employed as the computational domain. The simulations were conducted using the ANSYS-FLUENT commercial CFD package. The k-ε, P-1 and Probability Density Function (PDF) in the non-premix combustion, as low computational cost numerical approaches were employed to simulate the turbulent fluid flow, thermal radiation, combustion and conjugate heat transfer inside the furnace. Temperature measurement at different locations of the forged ingot surfaces were used to validate the transient numerical simulations. Good agreement was obtained between the predictions of the CFD model and the experimental measurements, demonstrating the reliability of the proposed approach and application of the model for process optimization purposes. Detailed analysis of conjugate heat transfer together with the turbulent combustion showed that the temperature evolution of the product was significantly dependant on the furnace geometry and the severity of turbulent flow structures in the furnace.


Author(s):  
Tariq Amin Khan ◽  
Nasir Mehdi Gardezi ◽  
Wei Li ◽  
Yang Zhou ◽  
Zahid Ayub

Abstract The performance on the air side flow is often limited due to its lower heat transfer coefficient. This work is related to numerical simulation to study the significance of employing delta winglets in flat finned and wavy finned-tube heat exchangers. For this purpose, three-dimensional simulation data and a multi-objective genetic algorithm are employed. The angle of attack (α) of delta winglets and Reynolds number varied from 15° to 75° and 500 to 1300, respectively. Employing delta winglets has increased the heat transfer per unit temperature and per unit volume (Z) and the fan power per unit core volume (E) for both flat finned and wavy finned-tube heat exchangers. To achieve a maximum heat transfer enhancement and a minimum friction factor, the optimal values of these parameters (Re and α) are calculated using the Pareto optimal strategy. For this purpose, CFD data, a surrogate model (neural network) and a multi-objective optimization genetic algorithm are combined. Results show that the performance of wavy finned-tube heat exchangers is higher than flat-finned tube heat exchangers which signify the importance of delta winglets in the wavy finned-tube heat exchangers.


2006 ◽  
Vol 128 (8) ◽  
pp. 819-828 ◽  
Author(s):  
Jixiang Yin ◽  
Guojun Li ◽  
Zhenping Feng

This paper reported three-dimensional numerical simulations of the steady laminar flow and heat transfer in corrugated-undulated channels with sinusoidal waves, aiming to investigate the effects of intersection angles (θ) between corrugated and undulated plate and Reynolds number (Re) on the flow and heat transfer. The simulations are conducted by using multi-channel computational domain for three different geometries. The code is validated against experimental results and then data for Nusselt number (Nu) and friction factor (f) are presented in a Re range of 100-1500, and intersection angle range of 30-150deg. The simulation confirms the changes of Nuu (averaged over undulated plate) and Nuc (averaged over corrugated plate) with θ representing different characteristics. As θ increases, Nu (Nuu or Nuc) is about 2–16 times higher for the corrugated-undulated configurations CP-UH1 and CP-UP1 and the concomitant f is about 4–100 higher, when compared to a straight channel having square cross section. The minimum of local Nu ( Nuu or Nuc ) is situated at the four contact points where the top plate touches the bottom one, and the high Nu is located upstream of the crest of the conjugate duct. Performance evaluation for the CP-UH1 channel shows that the goodness factors (G) are larger than 1 with the straight channel having a square cross section as a reference, and the 30deg geometry channel has optimal flow area goodness.


Author(s):  
Ramesh Narayanaswamy ◽  
Tilak T. Chandratilleke ◽  
Andrew J. L. Foong

A numerical study was conducted to investigate the fluid flow and heat transfer characteristics of a square microchannel with four longitudinal internal fins. Three-dimensional numerical simulations were performed on the microchannel with variable fin height ratio in the presence of a hydrodynamically developed, thermally developing laminar flow. Constant heat flux boundary conditions were assumed on the external walls of the square microchannel. Results of local Nusselt number distribution along the channel length were obtained as a function of the fin height ratio. The analysis was carried out for different fin heights and flow parameters. Interesting observations that provide more physical insight on this passive enhancement technique, and the existence of an optimum fin height is brought out in the present study.


1999 ◽  
Vol 121 (4) ◽  
pp. 1002-1010 ◽  
Author(s):  
E. Yu ◽  
Y. K. Joshi

A three-dimensional investigation of combined conduction, natural convection, and radiation in a side-vented compact enclosure is carried out. The focus of the study is on the enhancement of overall heat transfer through the opening, and the roles of the various modes in achieving it. A discrete heat source, flush-mounted centrally on a vertical substrate, is placed in the enclosure with a single rectangular opening on the opposite vertical wall. Steady-state computations are carried out for Rayleigh numbers, Ra, at 2.6 × 106 and 2.0 × 107. The results show that radiation plays a significant role in the overall heat transfer, and the radiative transport is even more pronounced for lower Ra. It is found that natural convection is weakened by radiation, however, contrary to the existing studies on top vented enclosures, the overall heat transfer is enhanced when radiation is included in the computations. Flow recirculation by radiative heating of enclosure walls is predicted, and is also observed experimentally. Heat spreading in the substrate is found to effect both convection and radiation. The numerical solutions on an extended computational domain are found in good agreement with the experimental data, when the conjugate effects are accounted for.


Author(s):  
Vijay K. Garg ◽  
Reza S. Abhari

The predictions from a three-dimensional Navier-Stokes code have been compared to the Nusselt number data obtained on a film-cooled, rotating turbine blade. The blade chosen is the ACE rotor with five rows containing 93 film cooling holes covering the entire span. This is the only film-cooled rotating blade over which experimental heat transfer data is available for the present comparison. Over 2.25 million grid points are used to compute the flow over the blade. Usually in a film cooling computation on a stationary blade, the computational domain is just one spanwise pitch of the film-cooling holes, with periodic boundary conditions in the span direction. However, for a rotating blade, the computational domain consists of the entire blade span from hub to tip, as well as the tip clearance region. As far as the authors are aware of, the present work offers the first comparison of the prediction of surface heat transfer using a three dimensional CFD code with film injection and the measured heat flux on a fully film-cooled rotating transonic turbine blade. In a detailed comparison with the measured data on the suction surface, a reasonably good comparison is obtained, particularly near the hub section. On the pressure surface, however, the comparison between the data and the prediction is poor. A potential reason for the discrepancy on the pressure surface could be the presence of unsteady effects due to stator-rotor interaction in the experiments which are not modeled in the present numerical computations.


Author(s):  
Seyed Saied Bahrainian

Most three dimensional tetrahedral grid generators can refine an initial grid in matter of seconds. But making an initial tetrahedral grid for complex geometry can be a tedious and time consuming task. This paper describes a novel procedure for generation of starting tetrahedral cells using hexahedral block topology. Hexahedral blocks are arranged around an aerodynamic body to fill-up a computational flow domain. Each of the hexahedral blocks is then decomposed into six tetrahedral elements to obtain an initial tetrahedral grid around the same aerodynamic body. This resulted in an algorithm that enables users to produce starting tetrahedral grids for variety of aerodynamic configurations. To construct an initial starting tetrahedral grid suitable for computational flow simulations, representing a solid surface geometry (fuselage or a wing section) attached to a plane-of-symmetry, a topology containing at least 5 hexahedral blocks is required. This results in an initial starting grid consisting of 30 tetrahedral cells with 74 faces and 16 vertices, which is the same number of vertices as for the hexahedral blocks. Since the number of vertices and their coordinate locations are kept the same, a connectivity matrix can be produced to describe the forming faces of the tetrahedral grid. This procedure was performed for a single block, 5-block, and 9-block topologies to produce starting tetrahedral cells for numerous domain size and shapes.


1996 ◽  
Vol 118 (1) ◽  
pp. 11-19 ◽  
Author(s):  
D. Bohn ◽  
R. Emunds ◽  
V. Gorzelitz ◽  
U. Kru¨ger

Increasing the thermal efficiency by higher turbine inlet temperatures is one of the most important aims in the area of gas turbine development. Because of the high temperatures, the turbine vanes and blades have to be cooled, and also knowledge of the mechanically and thermally stressed parts in the hottest zones of the rotor is of great interest. The prediction of the temperature distribution in a gas turbine rotor containing closed, gas-filled cavities, for example, in between two disks, has to account for the heat transfer conditions encountered in these cavities. In an entirely closed annulus, forced convection is not present, but a strong natural convection flow exists, induced by a nonuniform density distribution in the centrifugal force field. In Bohn et al. (1994), experimental and numerical investigations on rotating cavities with pure centripetal heat flux had been carried out. The present paper deals with investigations on a pure axially directed heat flux. An experimental setup was designed to realize a wide range of Ra numbers (2·108< Ra < 5·1010) usually encountered in cavities of gas turbine rotors. Parallel to the experiments, numerical calculations have been conducted. The numerical results are compared with the experimental data. The numerical scheme is also used to account for the influence of Re on heat transfer without changing Ra. This influence could not be pointed out by experiments, because a variation of the Re–Ra characteristic of the employed annuli was not possible. It was found that the numerical and experimental data are in quite good agreement, with exception of high Ra, where the numerical scheme predicts higher heat transfer than the experiments show. One reason may be that in the experiments the inner and outer cylindrical walls were not really adiabatic, an assumption used in the numerical procedure. Moreover, the assumption of a two-dimensional flow pattern may become invalid for high Ra. The influence of three-dimensional effects was studied with the three-dimensional version of the numerical code. In contrast to the radial directed heat transfer, it was found that Nu is much smaller and depends strongly on Re, whereas the radial heat transfer is only weakly influenced by Re.


Sign in / Sign up

Export Citation Format

Share Document