A multiscale thermo-fluid computational model for a two-phase cooling system

Abstract In this paper, we describe a mathematical model and a numerical simulation method for the condenser component of a novel two-phase thermosiphon cooling system for power electronics applications. The condenser consists of a set of roll-bonded vertically mounted fins among which air flows by either natural or forced convection. In order to deepen the understanding of the mechanisms that determine the performance of the condenser and to facilitate the further optimization of its industrial design, a multiscale approach is developed to reduce as much as possible the complexity of the simulation code while maintaining reasonable predictive accuracy. To this end, heat diffusion in the fins and its convective transport in air are modeled as 2D processes while the flow of the two-phase coolant within the fins is modeled as a 1D network of pipes. For the numerical solution of the resulting equations, a Dual Mixed-Finite Volume scheme with Exponential Fitting stabilization is used for 2D heat diffusion and convection while a Primal Mixed Finite Element discretization method with upwind stabilization is used for the 1D coolant flow. The mathematical model and the numerical method are validated through extensive simulations of realistic device structures which prove to be in excellent agreement with available experimental data.

[1]  Frank P. Incropera,et al.  Fundamentals of Heat and Mass Transfer , 1981 .

[2]  Christopher E. Brennen,et al.  Fundamentals of Multiphase Flow , 2005 .

[3]  Davide Del Col,et al.  Condensation inside and outside smooth and enhanced tubes: a review of recent research , 2003 .

[4]  J. Gillis,et al.  Matrix Iterative Analysis , 1961 .

[5]  Jean E. Roberts,et al.  Global estimates for mixed methods for second order elliptic equations , 1985 .

[6]  H. Gummel,et al.  Large-signal analysis of a silicon Read diode oscillator , 1969 .

[7]  Michel Fortin,et al.  A mixed finite element for the stokes problem using quadrilateral elements , 1995, Adv. Comput. Math..

[8]  M. Shah.  A general correlation for heat transfer during film condensation inside pipes , 1979 .

[9]  Ludmil T. Zikatanov,et al.  A monotone finite element scheme for convection-diffusion equations , 1999, Math. Comput..

[10]  Rebecca J. Shipley,et al.  Multiscale Modelling of Fluid and Drug Transport in Vascular Tumours , 2010, Bulletin of mathematical biology.

[11]  O. García-Valladares,et al.  Review of In-Tube Condensation Heat Transfer Correlations for Smooth and Microfin Tubes , 2003 .

[12]  Michael Herty.  Coupling Conditions for Networked Systems of Euler Equations , 2008, SIAM J. Sci. Comput..

[13]  Francesco Agostini,et al.  Simulation aided design of a two-phase thermosyphon for power electronics cooling , 2011, IECON 2011 - 37th Annual Conference of the IEEE Industrial Electronics Society.

[14]  Zhangxin Chen,et al.  On the implementation of mixed methods as nonconforming methods for second-order elliptic problems , 1995 .

[15]  J. Thome,et al.  Convective Boiling and Condensation , 1972 .

[16]  M. McLinden,et al.  NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport Properties-REFPROP, Version 9.1 | NIST , 2013 .

[17]  H. Baehr,et al.  An International Standard Formulation for the Thermodynamic Properties of 1,1,1,2‐Tetrafluoroethane (HFC‐134a) for Temperatures from 170 K to 455 K and Pressures up to 70 MPa , 1994 .

[18]  Riccardo Sacco,et al.  Stabilized Mixed Finite Volume Methods for Convection-Diffusion Problems , 1997 .

[19]  Cezar O.R. Negrão,et al.  A numerical simulation model for plate-type, roll-bond evaporators , 2008 .

[20]  P. Whalley.  Two-phase flow and heat transfer , 1996 .

[21]  Stefano Micheletti,et al.  Stability and error analysis of mixed finite-volume methods for advection dominated problems , 2006, Comput. Math. Appl..

[22]  Jean E. Roberts,et al.  Mixed and hybrid methods , 1991 .

[23]  P. Jenny,et al.  A Coupled Discrete/Continuum Model for Describing Cancer-Therapeutic Transport in the Lung , 2012, PLoS ONE.

[24]  Michael Herty,et al.  Gas Pipeline Models Revisited: Model Hierarchies, Nonisothermal Models, and Simulations of Networks , 2011, Multiscale Model. Simul..

[25]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[26]  Axel Klar,et al.  Gas flow in pipeline networks , 2006, Networks Heterog. Media.

[27]  Salomon Levy,et al.  Two-phase flow in complex systems , 1999 .

[28]  C. Dangelo,et al.  Multiscale modelling of metabolism and transport phenomena in living tissues , 2007 .

[29]  Lutz Tobiska,et al.  Numerical Methods for Singularly Perturbed Differential Equations , 1996 .

[30]  Martin Vohralík.  Equivalence between mixed finite element and multi-point finite volume methods , 2004 .

[31]  W. Rohsenow,et al.  Thermally Optimum Spacing of Vertical, Natural Convection Cooled, Parallel Plates , 1984 .

[32]  T. Gradinger,et al.  Roll-bond Condenser In A Two-phaseThermosyphon Loop For Power ElectronicsCooling , 2012 .

[33]  Alessio Fumagalli,et al.  A reduced model for Darcy’s problem in networks of fractures , 2014 .

[34]  R. R. P. Van Nooyen,et al.  A Petrov-Galerkin mixed finite element method with exponential fitting , 1995 .

[35]  Roland Span,et al.  Short Fundamental Equations of State for 20 Industrial Fluids , 2006 .

[36]  A. M. Zhuravlev,et al.  Thermophysical properties of air and air components , 1971 .