1 / 30

Mahmood Silieti Eduardo Divo Alain Kassab

An Inverse BEM/GA Approach to Determining Heat Transfer Coefficient Distributions Within Film Cooling Holes/Slots. Mahmood Silieti Eduardo Divo Alain Kassab. Mechanical, Materials, and Aerospace Engineering Department University of Central Florida, Orlando, FL, USA. Overview :. Motivation

keiki
Télécharger la présentation

Mahmood Silieti Eduardo Divo Alain Kassab

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. An Inverse BEM/GA Approach to Determining Heat Transfer Coefficient Distributions Within Film Cooling Holes/Slots Mahmood Silieti Eduardo Divo Alain Kassab Mechanical, Materials, and Aerospace Engineering Department University of Central Florida, Orlando, FL, USA IPE 2003 Tuscaloosa, Alabama

  2. Overview: • Motivation • Procedure • Problem Setup • Conjugate Heat Transfer Solution • Direct BEM Conduction Solution/Verification • Inverse Problem and Objective Function • Optimization Technique: Genetic Algorithms • Numerical Results • Conclusions and Extensions IPE 2003 Tuscaloosa, Alabama

  3. Motivation: Find end Wall Film Cooling Effectiveness and heat transfer coefficients (HTC) Can measure film effectiveness using optical thermography: which also provides To define endwall HTC. IPE 2003 Tuscaloosa, Alabama

  4. Closed loop Transonic Test Rig at UCF funded by SWPC IPE 2003 Tuscaloosa, Alabama

  5. To find endwall HTC: will solve 3D inverse conduction problem using endwall temperature measurements, however, HTC in film hole is unknown? There are 10 types of film cooling holes in this experiment, several of these are shaped and all are inclined. Objective of this feasibility study is to find a means of determining HTC correlation in film hole to be used in later 3D inverse problem analysis for endwall HTC Each type of film cooling hole is subject of a single hole calibration Experiment that will yield to be used in correlation IPE 2003 Tuscaloosa, Alabama

  6. h (or q) = ? Measured T & q Measured T & q Procedure: • Conjugate Heat Transfer (CHT) simulation of 2-D film cooling slot of end wall • Results from CHT simulation used to model experimentally measured surface heat flux and temperature. • T measured using temperature sensitive paint (TPS) • q measured using an optical thermographic technique under development at UCF • Inverse Problem: Input: T and q at exposed endwall surfaces Output: h (or q) at slot surface using the boundary element method (BEM) and a genetic algorithm (GA) IPE 2003 Tuscaloosa, Alabama

  7. Cooling Slot Setup for CHT Simulation: 2-D Film Cooling Slot IPE 2003 Tuscaloosa, Alabama

  8. Mesh has been created using Gambit (FLUENT grid generator) Fluid Grid Nodes=41,112 Solid Grid Nodes=2,104 IPE 2003 Tuscaloosa, Alabama

  9. CHT Simulation conditions chosen to match experiment to be carried out in wind-tunnel Main Air Flow: Turbulent Boundary Layer profile (1/7)th. Temperature= 350 K Coolant Air Flow: Uniform Pressure =105800 Pa Temperature= 300 K Fluid is Air: compressible, other properties are function of temperature Solid is Steel: properties are linear function of temperature IPE 2003 Tuscaloosa, Alabama

  10. CHT Solver: • Commercial Code “Fluent” Finite Volume • Full Navier-Stokes Equation for compressible turbulent flow “RNG “ CHT Results: • Results are converged at least for all residuals ( mass, momentum, energy, & ) IPE 2003 Tuscaloosa, Alabama

  11. IPE 2003 Tuscaloosa, Alabama

  12. IPE 2003 Tuscaloosa, Alabama

  13. IPE 2003 Tuscaloosa, Alabama

  14. IPE 2003 Tuscaloosa, Alabama

  15. IPE 2003 Tuscaloosa, Alabama

  16. measured T Direct BEM Conduction Solution: • Numerical consistency check of BEM (in-house) and CHT (commercial) code. • BEM surface mesh and CHT surface mesh are different radial basis function (RBF) interpolation used to pass information from one grid to the other. • Input CHT wall temperatures at solid surfaces to BEM and check BEM computed heat fluxes. IPE 2003 Tuscaloosa, Alabama

  17. BEM Methodology • Surface Mesh only: (we use quadratic discontinuous elements) • Governing Equation: (Boundary Integral Equation for Laplace Eqn.) q(x) = -kT(x) / n C() = 1 if    C() = 1/2 if    Where: G(x,) = (-1/ 2k) ln r(x,) in 2D H(x,) = -kG(x,) / n IPE 2003 Tuscaloosa, Alabama

  18. BEM Methodology • Discretized BIE is collocated at the boundary points, leading to • Introducing Boundary Conditions: • Contour plot of direct BEM temperature distribution: IPE 2003 Tuscaloosa, Alabama

  19. 55 30 61 1 1 21 70 135 91 Direct BEM Results: BEM fluxes consistent with FLUENT fluxes IPE 2003 Tuscaloosa, Alabama

  20. 30 21 61 70 Direct BEM Results: Heat Fluxes @ Cooling Slot IPE 2003 Tuscaloosa, Alabama

  21. h (or q) = ? Inverse Problem: • Cauchy conditions (T and q) imposed at the surfaces exposed to hot and cold gases. • Both temperature and flux are unknown on the surfaces of the cooling slot. IPE 2003 Tuscaloosa, Alabama

  22. Anchor point BEM node Inverse Problem: • Identification of heat fluxes in the cooling slot to match over-specified boundary data at the exposed surfaces. • Parametric representation of heat flux in cooling slot using radial basis functions (RBF) • Objective function is to minimize IPE 2003 Tuscaloosa, Alabama

  23. 1 1 0 1 0 1 1 0 1 0 1 0 1 1 0 0 1 0 1 1 1 1 0 0 1 1 1 1 0 1 1 1 1 0 1 0 1 0 0 0 • Optimization Technique:Genetic Algorithms • Non-gradient-based global search technique based on Darwinian evolution and operated by rules of natural selection: “Survival of the fittest” genes • Represent the design variables by a string of binary bits. • Generate a population of individuals genetically characterized by one chromosome or binary string. • Evaluate the fitness of each individual to identify its likelihood of propagating its genetic material. • Select and reproduce pairs of individuals to generate new generation subject to a probability of mutation. Chromosome IPE 2003 Tuscaloosa, Alabama

  24. Optimization Technique:Genetic Algorithms Advantages: - Very robust - Almost guaranteed global optimal - Inherent regularization Disadvantage: - Very slow Solution: - Parallelize process in a Computer Cluster by assigning different individuals to different nodes in the cluster. (Very efficient parallelization as very little communication is necessary) IPE 2003 Tuscaloosa, Alabama

  25. Optimization Technique: Parallel Genetic Algorithms Parameters: - population size = 50 - probability of jump mutation = 4% - probability of creep mutation = 20% - number of bits per parameter = 8 (255 steps) - number of children = 1 - ellitistic generation = 1 - parameter bound = searches for q between qmin and qmax block#1 (-15,000 to 15,000) block#2 (-2,000 to 2,000) IPE 2003 Tuscaloosa, Alabama

  26. (T, q) =? (T, q) =? • Inverse BEM Results: Evolution of objective function for Heat fluxes IPE 2003 Tuscaloosa, Alabama

  27. 61 70 30 21 Inverse BEM Results: Temperature Distribution @ Cooling Slot IPE 2003 Tuscaloosa, Alabama

  28. 41 61 1 70 135 91 55 30 1 21 Inverse BEM Results: Heat Fluxes IPE 2003 Tuscaloosa, Alabama

  29. 61 70 30 21 Inverse BEM Results: Heat Fluxes @ Cooling Slots 1 5 0 0 Q G A Q C F D + Q A P 1 0 0 0 + + 5 0 0 + + + 0 6 1 6 2 6 3 6 4 6 5 6 6 6 7 6 8 6 9 7 0 E l e m e n t IPE 2003 Tuscaloosa, Alabama

  30. Conclusions and Extensions: • Methodology shows promise in predicting the temperature and the heat fluxes within the slot. • Add more anchor points to capture the changes in heat fluxes. Add a regularization term to reduce unwanted oscillations associated with more anchor points. • Need to study the effect of input error in temperature and heat flux on resolution. • Apply the methodology to multiple slots. • Apply the methodology to 3-d single and multiple film-cooling holes. IPE 2003 Tuscaloosa, Alabama

More Related