Local Source Based CFD Modeling of Effusion Cooling Holes: Validation and Application to an Actual Combustor Test Case

Author(s):  
Antonio Andreini ◽  
Riccardo Da Soghe ◽  
Bruno Facchini ◽  
Lorenzo Mazzei ◽  
Salvatore Colantuoni ◽  
...  

State-of-the-art liner cooling technology for modern combustion chambers is represented by effusion cooling (or full-coverage film cooling). Effusion is a very efficient cooling strategy typically based on the use of several inclined small diameter cylindrical holes, where liner temperature is controlled by the combined protective effect of coolant film and heat removal through forced convection inside each hole. A CFD-based thermal analysis of such components implies a significant computational cost if the cooling holes are included in the simulations; therefore many efforts have been made to develop lower order approaches aiming at reducing the number of mesh elements. The simplest approach models the set of holes as a uniform coolant injection, but it does not allow an accurate assessment of the interaction between hot gas and coolant. Therefore higher order models have been developed, such as those based on localized mass sources in the region of hole discharge. The model presented in this paper replaces the effusion hole with a mass sink on the cold side of the plate, a mass source on the hot side, whereas convective cooling within the perforation is accounted for with a heat sink. The innovative aspect of the work is represented by the automatic calculation of the mass flow through each hole, obtained by a run time estimation of isentropic mass flow with probe points, while the discharge coefficients are calculated at run time through an in-house developed correlation. In the same manner, the heat sink is calculated from a Nusselt number correlation available in literature for short length holes. The methodology has been applied to experimental test cases of effusion cooling plates and compared to numerical results obtained through a CFD analysis including the cooling holes, showing a good agreement. A comparison between numerical results and experimental data was performed on an actual combustor as well, in order to prove the feasibility of the procedure.

Author(s):  
A. Andreini ◽  
R. Da Soghe ◽  
B. Facchini ◽  
L. Mazzei ◽  
S. Colantuoni ◽  
...  

State-of-the-art liner cooling technology for modern combustion chambers is represented by effusion cooling (or full-coverage film cooling). Effusion is a very efficient cooling strategy typically based on the use of several inclined small diameter cylindrical holes, where liner temperature is controlled by the combined protective effect of coolant film and heat removal through forced convection inside each hole. A CFD-based thermal analysis of such components implies a significant computational cost if the cooling holes are included in the simulations, therefore many efforts have been made to develop lower order approaches aiming at reducing the number of mesh elements. The simplest approach models the set of holes as a uniform coolant injection, but it does not allow an accurate assessment of the interaction between hot gas and coolant. Therefore higher order models have been developed, such as those based on localized mass sources in the region of hole discharge. The model presented in this paper replaces the effusion hole with a mass sink on the cold side of the plate, a mass source on the hot side, whereas convective cooling within the perforation is accounted for with a heat sink. The innovative aspect of the work is represented by the automatic calculation of the mass flow through each hole, obtained by a run time estimation of isentropic mass flow with probe points, while the discharge coefficients are calculated at run time through an in-house developed correlation. In the same manner the heat sink is calculated from a Nusselt number correlation available in literature for short length holes. The methodology has been applied to experimental test cases of effusion cooling plates and compared to numerical results obtained through a CFD analysis including the cooling holes, showing a very good agreement. A comparison between numerical results and experimental data was performed on an actual combustor as well, in order to prove the feasibility of the procedure.


Author(s):  
Sebastian Scholl ◽  
Catherine Gorle ◽  
Farzad Houshmand ◽  
Tanya Liu ◽  
Hyoungsoon Lee ◽  
...  

This study considers CFD simulations with conjugate heat transfer performed in the framework of designing a complex micro-scale cooling geometry. The numerical investigation of the three-dimensional, laminar flow (Reynolds number smaller than 480) and the solid conduction is done on a reduced model of the heat sink micro-structure to enable exploring a variety of configurations at a limited computational cost. The reduced model represents a unit-cell, and uses periodic and symmetry boundary conditions to mimic the conditions in the entire cooling manifold. A simulation of the entire heat sink micro-structure was performed to verify the unit-cell set-up, and the comparison demonstrated that the unit-cell simulations allow reducing the computational cost by two orders of magnitude while retaining accurate results. The baseline design for the unit-cell represents a configuration used in traditional electronic heat sinks, i.e. a simple channel geometry with a rectangular cross section, with a diameter of 50 μm, where the fluid flows between two cooling fins. Subsequently three types of modified geometries with feature sizes of 50 μm were considered: baffled geometries that guide the flow towards the hotspot region, geometries where the fins are connected by crossbars, and a woodpile structure without cooling fins. Three different mass-flow rates were tested. Based on the medium mass-flow rate considered, the woodpile geometry showed the highest heat transfer coefficient with an increase of 70% compared to the baseline geometry, but at the cost of increasing the pressure drop by more than 300%. The crossbar geometries were shown to be promising configurations, with increases in the heat transfer coefficient of more than 20% for a 70% increase in pressure drop. The potential for further optimization of the crossbar configurations by adding or removing individual crossbars will be investigated in a follow up study. The results presented demonstrate the increase in performance that can be obtained by investigating a variety of designs for single phase cooling devices using unit-cell conjugate heat transfer simulations.


Author(s):  
Yufang Zhang ◽  
Ke Wang

Full-coverage film cooling, based on use of closely spaced inclined small diameter cylindrical holes, is considered as a very efficient cooling approach for modern combustion chambers. However, a CFD-based thermal analysis of such components is very expensive if mesh is generated in numerous cooling holes. In order to reduce the computational cost, several models have been developed by researchers. Among them, source term model is one of the most potential models. In this model, a film cooling hole is replaced by a mass source and a mass sink in the fluid domain, and an energy sink is used to simulate the convective cooling inside hole. In this paper, the source term model was used for RANS simulation of a multi-perforated plate case, and the effect of grid setup was analyzed. It has been concluded that the height of first prism layer has an important effect on the numerical results, whereas effect of maximum grid size on coupled walls is relatively small. Then, the source term model was used for numerical simulations of a single dome combustor case. Numerical results of flow split and pressure losses in non-reacting case agreed well with experimental results. In reacting case, a good agreement between numerical and experimental results, including liner temperature and pressure loss, was also obtained.


2000 ◽  
Author(s):  
X. Wei ◽  
Y. Joshi

Abstract A novel heat sink based on a multi-layer stack of liquid cooled microchannels is investigated. For a given pumping power and heat removal capability for the heat sink, the flow rate across a stack of microchannels is lower compared to a single layer of microchannels. Numerical simulations using a computationally efficient multigrid method [1] were carried out to investigate the detailed conjugate transport within the heat sink. The effects of the microchannel aspect ratio and total number of layers on thermal performance were studied for water as coolant. A heat sink of base area 10 mm by 10 mm with a height in the range 1.8 to 4.5 mm (2–5 layers) was considered with water flow rate in the range 0.83×10−6 m3/s (50 ml/min) to 6.67×10−6 m3/s (400 ml/min). The results of the computational simulations were also compared with a simplified thermal resistance network analysis.


2016 ◽  
Vol 78 (10-2) ◽  
Author(s):  
Nik Ahmad Faiz Nik Mazlam ◽  
Normah Mohd-Ghazali ◽  
Thierry Mare ◽  
Patrice Estelle ◽  
Salma Halelfadl

The microchannel heat sink (MCHS) has been established as an effective heat removal system in electronic chip packaging. With increasing power demand, research has advanced beyond the conventional coolants of air and water towards nanofluids with their enhanced heat transfer capabilities. This research had been carried out on the optimization of the thermal and hydrodynamic performance of a rectangular microchannel heat sink (MCHS) cooled with carbon nanotube (CNT) nanofluid, a coolant that has recently been discovered with improved thermal conductivity. Unlike the common nanofluids with spherical particles, nanotubes generally come in cylindrical structure characterized with different aspect ratios. A volume concentration of 0.1% of the CNT nanofluid is used here; the nanotubes have an average diameter and length of 9.2 nm and 1.5 mm respectively. The nanofluid has a density of 1800 kg/m3 with carbon purity 90% by weight having lignin as the surfactant. The approach used for the optimization process is based on the thermal resistance model and it is analyzed by using the non-dominated sorting multi-objective genetic algorithm. Optimized outcomes include the channel aspect ratio and the channel wall ratio at the optimal values of thermal resistance and pumping power. The optimized results show that, at high operating temperature of 40°C the use of CNT nanofluid reduces the total thermal resistance by 3% compared to at 20°C and consequently improve the thermal performance of the fluid. In terms of the hydrodynamic performance, the pumping power is also being reduced significantly by 35% at 40°C compared to the lower operating temperature.  


2021 ◽  
Author(s):  
Andisheh Tavakoli ◽  
Kambiz Vafai

Abstract The present study analyzes the optimal distribution of a limited amount of high thermal conductivity material to enhance the heat removal of circular 3D integrated circuits, IC. The structure of the heat spreader is designed as a composite of high thermal conductivity (Boron Arsenide) and moderate thermal conductivity (copper) materials. The volume ratio of high-conductivity inserts to the total volume of the spreader is set at a fixed pertinent ratio. Two different boundary conditions of constant and variable temperature are considered for the heat sink. To examine the impact of adding high-conductivity inserts on the cooling performance of the heat spreader, various patterns of the single and double ring inserts are studied. A parametric study is performed to find the optimal location of the rings. Moreover, the optimal distribution of the high-conductivity material between the inner and outer rings is found. The results show that for the optimal conditions, the maximum temperature of the 3D IC is reduced up to 10%; while the size of the heat sink, and heat spreader can be diminished by as much as 200%.


Author(s):  
Дарья Владимировна Добролюбова ◽  
Элла Петровна Шурина

Рассматриваются особенности применения модифицированной вариационной постановки векторного метода конечных элементов (ВМКЭ), основанной на замене тонких сильнопроводящих объектов токонесущими поверхностями, для моделирования гармонического электрического поля в областях с криволинейными экранами при различном типе возбуждения поля. Исследуется применимость модифицированной вариационной постановки в широком диапазоне частот Purpose. The paper addresses applicability of the modified variational formulation of vector FEM for the harmonic electric field to the media with cylindrical shields. Thin highly conductive objects are treated as surfaces with the equivalent surface current density. We consider the excitation of the field by a local source (current loop) located either inside or outside the cylindrical shield. Methodology. The simulations are carried out on unstructured tetrahedral meshes. Since the modified variational formulation treats thin highly conductive objects as surfaces, only the surface of a cylinder is discretized. The results yielded by the modified variational formulation are compared with the results of the classic vector FEM. Findings. For the frequency range between 100 KHz and 100 MHz, the modified variational formulation provides correct results when the field source is located outside the cylindrical shield. The modified variational formulation reduces computational cost, since the volume of the thin shield is not discretized. When the field source is located inside the shield, the modified variational formulation gives valid results only in the proximity of the source. Originality/value. The limitations for the application of the reduced variational formulation for the modelling of harmonic electric field in the media with hollow cylindrical shields are investigated


Water ◽  
2019 ◽  
Vol 11 (10) ◽  
pp. 2033 ◽  
Author(s):  
Guodong Li ◽  
Guoding Chen ◽  
Pengfeng Li ◽  
Haixiao Jing

High-speed and accurate simulations of landslide-generated tsunamis are of great importance for the understanding of generation and propagation of water waves and for prediction of these natural disasters. A three-dimensional numerical model, based on Reynolds-averaged Navier–Stokes equations, is developed to simulate the landslide-generated tsunami. Available experiment data is used to validate the numerical model and to investigate the scale effect of numerical model according to the Froude similarity criterion. Based on grid convergence index (GCI) analysis, fourteen cases are arranged to study the sensitivity of numerical results to mesh resolution. Results show that numerical results are more sensitive to mesh resolution in near field than that in the propagation field. Nonuniform meshes can be used to balance the computational efficiency and accuracy. A mesh generation strategy is proposed and validated, achieving an accurate prediction and nearly 22 times reduction of computational cost. Further, this strategy of mesh generation is applied to simulate the Laxiwa Reservoir landslide tsunami. The results of this study provide an important guide for the establishment of a numerical model of the real-world problem of landslide tsunami.


Author(s):  
Filippo Rubechini ◽  
Michele Marconcini ◽  
Andrea Arnone ◽  
Stefano Cecchi ◽  
Federico Dacca`

A three-dimensional, multistage, Navier-Stokes solver is applied to the numerical investigation of a four stage low-pressure steam turbine. The thermodynamic behavior of the wet steam is reproduced by adopting a real-gas model, based on the use of gas property tables. Geometrical features and flow-path details consistent with the actual turbine geometry, such as cavity purge flows, shroud leakage flows and partspan snubbers, are accounted for, and their impact on the turbine performance is discussed. These details are included in the analysis using simple models, which prevent a considerable growth of the computational cost and make the overall procedure attractive as a design tool for industrial purposes. Shroud leakage flows are modeled by means of suitable endwall boundary conditions, based on coupled sources and sinks, while body forces are applied to simulate the presence of the damping wires on the blades. In this work a detailed description of these models is provided, and the results of computations are compared with experimental measurements.


Sign in / Sign up

Export Citation Format

Share Document