Transient temperature fields with general nonlinear boundary conditions in electronic systems

• Published in: IEEE transactions on components and packaging technologies 2005-03, Vol.28 (1), p.23-33
• Main Authors: Gerstenmaier, Y.C., Wachutka, G.K.M.
• Format: Article
• Language: English
• Subjects: Algorithms; Boundary conditions; Conducting materials; Green's function methods; Green's functions; Heat transfer; Integral equations; Mathematical analysis; Mathematical models; Mesh generation; Multi-chip-modules (MCMs); Nonhomogeneous media; Nonlinear equations; nonlinear heat conduction with Green's function method; Nonlinearity; radiative and convective boundary conditions; Representations; Robustness; Studies; Temperature dependence; Temperature distribution; thermal transient modeling
• ISSN: 1521-3331; 1557-9972
• DOI: 10.1109/TCAPT.2004.843186

Green's function representations of the solution of the heat conduction equation for general boundary conditions are generalized for the nonlinear, i.e., temperature dependent case. Temperature dependent heat transfer coefficients lead to additional terms in the Green's function representation of the temperature field. For a rectangular structure with averaged homogeneous material parameters several types of Green's functions can be chosen especially simple, because of the new representation with the possibility of differing types of boundary conditions for the temperature field and the Green's function. Exact finite closed form expressions for three-dimensional-Green's functions in the time domain using elliptic theta functions are presented. The temperature field is a solution of a nonlinear integral equation which is solved numerically by iteration. The resulting algorithm is very robust, stable and accurate with reliable convergence properties and avoids matrix inversions completely. The algorithm can deal with all sizes of volume heat sources without additional grid generation. Large and small size volume heat sources are treated simultaneously in the calculations that will be presented. Heat transfer coefficients are chosen representing radiative and convective boundary conditions. An extension of the solution algorithm to composed multilayer systems of arbitrary geometry is outlined.