TWI464603B - 熱分析模型建立方法及其非暫時性電腦可讀取紀錄媒體 - Google Patents
熱分析模型建立方法及其非暫時性電腦可讀取紀錄媒體 Download PDFInfo
- Publication number
- TWI464603B TWI464603B TW100120735A TW100120735A TWI464603B TW I464603 B TWI464603 B TW I464603B TW 100120735 A TW100120735 A TW 100120735A TW 100120735 A TW100120735 A TW 100120735A TW I464603 B TWI464603 B TW I464603B
- Authority
- TW
- Taiwan
- Prior art keywords
- thermal
- reference point
- power density
- density function
- temperature difference
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Description
本揭示內容是有關於一種模型建立方法,且特別是有關於一種熱分析模型建立方法及其非暫時性電腦可讀取紀錄媒體。
半導體晶片隨著製程的演進,包含愈來愈多電子元件於其中。因此,精確的熱分析模型將是半導體晶片設計中不可或缺的一環,以確保晶片運作時不會由於其空間中的熱能過多無法順利釋出而燒燬或損壞。為了進行分析與模擬,將需要對晶片內溫度的分佈進行計算,而晶片內溫度的分佈是由半導體晶片中各元件及各連接件的溫度來決定。然而,現行的分析及模擬技術常需要以非常大的運算量及複雜度,來對整個晶片中所有的元件及連接件進行計算,既耗時亦耗資源,對講求效率的晶片設計程序中,將是不利的。
因此,如何設計一個新的熱分析模型建立方法及其非暫時性電腦可讀取紀錄媒體,以對半導體晶片快速而有效地建立熱分析模型,乃為此一業界亟待解決的問題。
因此,本揭示內容之一態樣是在提供一種熱分析模型建立方法,包含下列步驟:根據電子系統所包含之複數電子模組建立電熱(electrothermal)π
網路模型,以定義電
子系統之熱源、熱傳遞路徑以及共基恆溫散熱面。定義電子系統中之複數個觀察點,其中各觀察點位於相對參考點之外圍空間之等溫表面上,其中參考點為熱源或觀察點其中之一。根據參考點發散至各觀察點之熱功率密度函數、參考點至各觀察點之熱導係數及參考點至觀察點之距離向量,計算各觀察點與參考點間之熱傳導溫度差或熱對流溫度差。根據共基恆溫散熱面之定義溫度,由各觀察點與參考點間之熱傳導溫度差或熱對流溫度差建立電子系統之溫度分佈。
依據本揭示內容一實施例,其中當參考點為熱源,係藉由熱功率密度函數、熱導係數及距離向量計算熱傳導溫度差。計算熱傳導溫度差更包含藉由熱功率密度函數、熱導係數及距離向量對觀察點所位於參考點對應之外圍空間進行體積分以及位於參考點對應之等溫表面進行面積分計算。其中熱功率密度函數包含體積熱功率密度函數qv以及表面熱功率密度函數qs,熱導係數包含體積熱導係數kv以及表面熱導係數ks,觀察點之位置為r,熱源之位置為r`,外圍空間之該等溫表面之一單位法向量為,該表面為a,該熱傳導溫度差係表示為:
依據本揭示內容另一實施例,其中當參考點為觀察點其中之一,係藉由熱功率密度函數、熱導係數及距離向量計算熱對流溫度差。計算熱對流溫度差更包含藉由熱功率密度函數、熱導係數及距離向量對觀察點所位於對應參考點之
溫表面進行面積分計算。其中熱功率密度函數包含表面熱功率密度函數qs,熱導係數包含表面熱導係數ks,觀察點之位置為r,熱源之位置為r`,外圍空間之等溫表面之單位法向量為,表面為a,熱傳導溫度差係表示為:
依據本揭示內容又一實施例,其中參考點至各觀察點係為熱傳遞路徑之一部份,且包含電子模組至少其中之一,熱導係數係與電子模組至少其中之一之材質相關。
依據本揭示內容再一實施例,其中根據熱分析模型建立方法更包含下列步驟:根據電子模組以及電熱π
網路模型建立包含複數熱阻之等效電熱電路模型以及根據熱功率密度函數、熱導係數及距離向量,計算等效電熱電路模型之熱阻。
依據本揭示內容更具有之一實施例,其中熱源為電子系統之最高溫處,共基恆溫散熱面為電子系統之最低溫處。
一種非暫時性電腦可讀取紀錄媒體,用以儲存一應用程式,該應用程式使一電腦裝置執行一種熱分析模型建立方法,包含下列步驟:根據電子系統所包含之複數電子模組建立電熱(electrothermal)π
網路模型,以定義電子系統之熱源、熱傳遞路徑以及共基恆溫散熱面。定義電子系統中之複數個觀察點,其中各觀察點位於相對參考點之外圍空間之等溫表面上,其中參考點為熱源或觀察點其中之一。根據參考點發散至各觀察點之熱功率密度函數、參考點至各觀察點之熱導係數及參考點至觀察點之距離向量,計算各觀察點與參考點間之熱傳導溫度差或熱對流溫度
差。根據共基恆溫散熱面之定義溫度,由各觀察點與參考點間之熱傳導溫度差或熱對流溫度差建立電子系統之溫度分佈。
依據本揭示內容一實施例,其中當參考點為熱源,係藉由熱功率密度函數、熱導係數及距離向量計算熱傳導溫度差。計算熱傳導溫度差更包含藉由熱功率密度函數、熱導係數及距離向量對觀察點所位於參考點對應之外圍空間進行體積分以及位於參考點對應之等溫表面進行面積分計算。其中熱功率密度函數包含體積熱功率密度函數qv以及表面熱功率密度函數qs,熱導係數包含體積熱導係數kv以及表面熱導係數ks,觀察點之位置為r,熱源之位置為r`,外圍空間之該等溫表面之一單位法向量為,該表面為a,該熱傳導溫度差係表示為:
依據本揭示內容另一實施例,其中當參考點為觀察點其中之一,係藉由熱功率密度函數、熱導係數及距離向量計算熱對流溫度差。計算熱對流溫度差更包含藉由熱功率密度函數、熱導係數及距離向量對觀察點所位於對應參考點之溫表面進行面積分計算。其中熱功率密度函數包含表面熱功率密度函數qs,熱導係數包含表面熱導係數ks,觀察點之位置為r,熱源之位置為r`,外圍空間之等溫表面之單位法向量為,表面為a,熱傳導溫度差係表示為:
依據本揭示內容又一實施例,其中參考點至各觀察點
係為熱傳遞路徑之一部份,且包含電子模組至少其中之一,熱導係數係與電子模組至少其中之一之材質相關。
依據本揭示內容再一實施例,其中根據熱分析模型建立方法更包含下列步驟:根據電子模組以及電熱π
網路模型建立包含複數熱阻之等效電熱電路模型以及根據熱功率密度函數、熱導係數及距離向量,計算等效電熱電路模型之熱阻。
依據本揭示內容更具有之一實施例,其中熱源為電子系統之最高溫處,共基恆溫散熱面為電子系統之最低溫處。
應用本揭示內容之優點係在於藉由電熱π
網路模型之建立確定熱源、熱傳遞路徑以及共基恆溫散熱面後,藉由Green Theorem所發展出來的溫度差積分函數計算電子系統中的觀察點溫度差,建立電子系統中的溫度分佈,而輕易地達到上述之目的。
請參照第1圖。第1圖為本揭示內容一實施例中,熱分析模型建立方法之流程圖。請同時參照第2A圖。第2A圖為本揭示內容一實施例之電子系統2之立體圖。第1圖中之熱分析模型建立方法可用以建立第2A圖所繪示之電子系統2之熱分析模型。於一實施例中,熱分析模型建立方法可藉由一非暫時性電腦可讀取紀錄媒體中,如唯讀記憶體(read only memory;ROM)、快閃記憶體(Flash Memory)、軟碟、硬碟、光碟、隨身碟、磁帶、可由網路存取之資料庫或熟悉此技術者可輕易思及具有相同功能之
儲存媒體,所儲存一應用程式來執行。
於本實施例中,第2A圖所繪示之電子系統2包含垂直共振腔面射型雷射20(Vertical-Cavity Surface-Emitting Laser;VCSEL)。於不同實施例中,垂直共振腔面射型雷射20所包含之通道數目可依情形調整。於本實施例中,垂直共振腔面射型雷射20具有四個通道。隨著電腦的運算速度高達數GHz以上,傳統利用銅線傳遞電訊號的頻寬將不敷使用,因此利用光作為高速傳遞媒介之光連結技術將大幅取代電子產品中之各種訊號的連結。垂直共振腔面射型雷射20藉由覆晶接合(flip-chip bonding)技術形成於矽基光學平台(Silicon optical bench;SiOB)22上,以做為傳輸速度高達Gbps的光學訊號傳輸裝置,並透過光纖24來進行傳輸。光纖24於不同實施例中,可嵌於矽基光學平台22上的V形槽或是U形槽(未繪示)中。請參照第2B圖,為第2A圖的電子系統2中,虛框部份沿D方向之側視圖。如第2B圖所示,垂直共振腔面射型雷射20底部具有一介電層26,介電層26上可形成散熱孔26’(形成於垂直共振腔面射型雷射20正下方者由於角度關係未繪示),並藉由連接墊28連接至外部,以達到散熱之功效。
請同時參照第3圖。熱分析模型建立方法之步驟101,是根據電子系統2所包含之電子模組建立如第3圖所示之電熱π
網路模型3。電子系統2之電子模組即包含前述的垂直共振腔面射型雷射20、矽基光學平台22、光纖24、介電層26、散熱孔26’以及連接墊28等等。如第3圖所示,電熱π
網路模型3可定義出電子系統2之熱源30、熱傳遞
路徑32以及共基恆溫散熱面34。其中,以虛線框繪示者表示絕熱面。各個阻抗Z1、Z2、Z3與Z4表示熱源30內的熱阻、熱對流、熱電感以及其他邊界情形造成的熱阻效應。邊界A及邊界B可為熱流的共基端、外部邊界或是與其他系統之連接面。一般來說,熱源30為電子系統2之最高溫處,共基恆溫散熱面34為電子系統2之最低溫處。於本實施例中,熱源30即為垂直共振腔面射型雷射20,熱傳遞路徑32可包含光纖24、介電層26、散熱孔26’以及連接墊28等,而共基恆溫散熱面34則為矽基光學平台22。
藉由電熱π
網路模型3的建立,電子系統2對應的熱模型可以迅速地建立,而降低熱分析之複雜度。因此,熱流行為可以藉由一些數學方法計算溫度分佈而被推導出。需注意的是,於其他的電子系統中,所據以建立的電熱π
網路模型可能與第3圖所繪示之有所差異,不為第3圖所繪示的所限。
請參照第4圖。第1圖之步驟102,是根據第2A圖之電子系統2中的各電子模組以及第3圖的電熱π
網路模型3建立如第4圖所示包含複數熱阻之等效電熱電路模型4(equivalent electrothermal circuit model;EETCM)。
第3圖中的阻抗Z1及阻抗Z2為垂直共振腔面射型雷射20的熱阻,意即第4圖中的熱阻RVCSEL以及R’VCSEL。阻抗Z3是由於忽略空氣熱對流而為一無限大阻值之熱阻之熱阻。阻抗Z4為矽基光學平台22之熱電感,即第4圖之CSiOB。第3圖之熱傳遞路徑32可包含串聯之金層及BCB層(RGOLD及RBCB),並與空氣的熱阻並聯(RAir)。第3
圖之共基恆溫散熱面34則包含共基阻抗RSiOB。第4圖中所繪示的開關則表示各通道的垂直共振腔面射型雷射20可經由切換而啟動或關閉。由於熱源30之熱行為複雜度高,因此需要以新的方式來計算溫度分佈,以瞭解其行為。
格林定律(Green’s Theorem)為用以處理特定體積之邊界問題。當此體積內的場源及圍繞此特定體積之表面的淨場通量均已知,則此特定體積內的場以及表面之通量之物理特性可以被決定。由第二格林定律可知:
其中以及可為兩個任意的向量場。、V
'、以及a
'分別為場源之位置(參考點)、欲處理的特定體積、圍繞此特定體積之表面之單位法向量以及圍繞此特定體積之表面。為處理熱分佈問題,上述之兩個向量場可分別定義為:
其中T及分別為溫度分佈以及參考點。請參照第5圖。第5圖為本揭示內容一實施例中,特定體積5之示意圖。需注意的是,於本實施例中,此特定體積5為球體,並於球體之中央包含有熱源。如第5圖所示,O點為原點,標示為之處為熱源(參考點)之位置,標示為之處則為觀察點之位置。因此,觀察點所觀察到的溫度,是受到熱源經由熱傳導(conduction)形式傳達至觀察點之熱的影響。
將式(2)及(3)代入式(1),並利用下列等式:
溫度之分佈函數可以推導為下式:
為簡化邊界條件,此表面,亦即圍繞特定體積5之表面50可設為一等溫表面。因此,式(5)中的第三項可成為一常數分佈,式(5)可變為一溫度差之函數:
因此,藉由體積分及面積分的處理,此特定體積5內,參考點與觀察點間的溫度差可由上式推出。
式(6)中的熱源之體積分及面積分可以下列式子表示:
其中qv及qs分別為欲處理之特定體積5內的體積功率密度(W
/m 3
)以及等溫表面之表面功率密度(W
/m 2
)。Kv及ks則為體積熱導係數以及表面熱導係數(均為W
/m
.K
)。就物理意義上來說,於式(7)及式(8)中的功率密度為一特定觀察點上所觀察到的熱源的體積功率密度以及
表面功率密度。將式(7)及式(8)代入式(6)後,成為:
如熱源之功率密度為常數且熱之傳導為等向(isotropic),式(9)可進一步簡化為:
式(10)可再藉由假設3kv=ks=k而進一步簡化為一維的問題,而上述之觀察點與參考點間的距離向量差即為熱傳導路徑△L
。因此式(10)可變為:
或是
式(11)暗示了自等溫表面50散發出的功率密度減去特定體積5所產生的功率密度,將等於溫度梯度含數以及一維熱傳導係數的乘積。而在一個沒有熱源的特定體積中,式(12)將等同於熱傳導之傅立葉定律。
同時,如第4圖所繪示之等效電熱電路模型4中的各熱阻R,可藉由上述之計算推導出:
請參照第6圖。第6圖為本揭示內容一實施例中,特
定體積6之示意圖。需注意的是,於本實施例中,此特定體積6為球體,球體內並無熱源。如第6圖所示,O點為原點,標示為之處為參考點之位置,標示為T
∞之處則為觀察點之位置。因此,觀察點所觀察到的溫度,是受到參考點經由熱對流(convection)形式傳達至觀察點之熱的影響。
如式(6)中並未有造成熱傳導之熱源的存在,將變為:
熱對流機制對表面60的影響可以延伸至其周圍的空間,其中等溫表面之溫度T
∞相當於一個常數,且不會受到內部熱系統之影響。因此,觀察點之溫度分佈亦可由溫度差所表示出:
或是
其中h相當於k
/△L
,並被定義為熱對流常數(W
/m 2
.K
)。
而熱的第三個傳導機制,即熱幅射,由於僅在非常高溫的情形下才有較大的影響力,因此在此忽略不計。
請再次參考第1圖。因此,第1圖之步驟103中,定義電子系統2中之複數個觀察點,其中各觀察點位於相對參考點之外圍空間之等溫表面上,其中參考點為熱源或觀
察點其中之一。接著於步驟104,根據參考點發散至各觀察點之熱功率密度函數、參考點至各觀察點之熱導係數及參考點至觀察點之距離向量,計算各觀察點與參考點間之熱傳導溫度差或熱對流溫度差。其中熱傳導溫度差可藉由式(9)計算,而熱對流溫度差可藉由式(10)計算。其中,觀察點與參考點間即為熱傳遞路徑之一部份,其介質即為觀察點與參考點間的電子模組如第2圖中的光纖24、介電層26、散熱孔26’以及連接墊28等等,且其熱導係數將與此電子模組的材質相關。
於步驟105,根據第3圖中的共基恆溫散熱面34之定義溫度,亦即電子系統2之最低溫度,可由各觀察點與參考點間之熱傳導溫度差或熱對流溫度差建立電子系統2之溫度分佈。於一實施例中,共基恆溫散熱面34即為電子系統2中的散熱模組,可控制在此表面上的溫度為一個電子系統2之最低溫度的定值。由於熱源(於本實施例中為垂直共振腔面射型雷射20)之溫度較難直接決定,因此可藉由共基恆溫散熱面34的溫度以及熱傳導溫度差或熱對流溫度差反推整個電子系統2的溫度。
於步驟106,根據熱功率密度函數、熱導係數及距離向量,可進一步根據式(13)計算第4圖中,等效電熱電路模型4之熱阻。因此,如將電子系統2的空間切割為較大的數個單位體積,整體電子系統2的溫度分佈將可藉由上述步驟快速地計算出。如欲進行更精密的溫度分佈計算,可將電子系統2以更小的單位體積進行切割。愈小的單位體積意謂需要愈複雜的計算,第4圖之等效電熱電路模型
4可能將不足以表達其精確度,因此可另設計其他演算法輔助計算,然而其核心仍為上述之電熱π
網路模型的建立以及熱對流、熱傳導溫度差之計算,任何熟習此技藝者,在不脫離本揭示內容之精神和範圍內,當可作各種之更動與潤飾。需注意的是,單位體積之形狀亦不限於上述實施例所揭露之球體,亦可以立方體、梯形體等形狀進行電子系統2空間的切割。
於一實施例中,四個垂直共振腔面射型雷射形成於625微米厚的矽基光學平台上,並運作於一般的空氣中而未有其他形式的熱對流。矽基光學平台底部為等溫表面,其溫度為攝氏溫度75度。其他表面均為絕熱面,因此均無熱流於其上。由於空氣之熱阻過大,且矽基光學平台之熱阻過小,因此在第4圖中的等效電熱電路模型4中可忽略不計。因此,關於其材質的運算均可忽略而降低整體運算的時間。
垂直共振腔面射型雷射的體積功率密度可藉由垂直共振腔面射型雷射每單位體積的輸入功率以及其運作效能之乘積而得,亦即:qv
=(I 2 R h
/V h
)η
(17)
其中I
、R h
、V h
及η
分別為熱源之輸入電流、電阻、體積以及電熱能轉換效能常數。請參照第7圖,第7圖為一實施例中,單一垂直共振腔面射型雷射啟動時,電子系統2之等溫表面之示意圖。如第7圖所示,當僅有一個垂直共振腔面射型雷射啟動時,等溫表面之結構顯示了電子系統中的最高溫之處為A點。其中,熱源區域是由「+」記
號繪示。於本實施例中,輸入電流約為8微安培,電熱能轉換效能常數為85.7%。此最高溫處的溫度約為攝氏溫度78.4度。
在決定最高溫點後,溫度分佈即可由式(9)及(14)決定。而各元件之熱阻亦可由式(13)決定。請參照第8A圖。第8A圖為本揭示內容一實施例中,忽略空氣及矽基光學平台之熱阻的情形下,僅啟動單一垂直共振腔面射型雷射時的等效電熱電路模型8。由於垂直共振腔面射型雷射的熱傳導性弱於散熱孔,且運作中的垂直共振腔面射型雷射與其鄰接的垂直共振腔面射型雷射間的厚度,相較於散熱孔有兩個數量級的差距,因此主要的熱傳導路徑即為散熱孔。
第8B圖為本揭示內容一實施例中的電子系統模擬結果示意圖。模擬結果顯示電子系統在單一垂直共振腔面射型雷射以8微安培輸入電流以及2伏特偏壓之情形下運作,於P1點、P2點及P3點之溫度分別為攝氏溫度78.9度、76.8度以及75度。最高溫的P1點接近如第7圖所示,由等效電熱電路模型所計算之結果。最高溫之點將位於P1點,而主要的熱流方向將由P1點傳至P3點。因此,根據電熱π
網路模型建立的等效電熱電路模型,可提供快速而正確的熱分析及模型之建立。
雖然本揭示內容已以實施方式揭露如上,然其並非用以限定本揭示內容,任何熟習此技藝者,在不脫離本揭示內容之精神和範圍內,當可作各種之更動與潤飾,因此本揭示內容之保護範圍當視後附之申請專利範圍所界定者為
準。
101-106‧‧‧步驟
2‧‧‧電子系統
20‧‧‧垂直共振腔面射型雷射
22‧‧‧矽基光學平台
24‧‧‧光纖
26‧‧‧介電層
26’‧‧‧散熱孔
28‧‧‧連接墊
3‧‧‧電熱π
網路模型
30‧‧‧熱源
32‧‧‧熱傳遞路徑
34‧‧‧共基恆溫散熱面
4、8‧‧‧等效電熱電路模型
5、6‧‧‧特定體積
50、60‧‧‧表面
為讓本揭示內容之上述和其他目的、特徵、優點與實施例能更明顯易懂,所附圖式之說明如下:第1圖為本揭示內容一實施例中,熱分析模型建立方法之流程圖;第2A圖為本揭示內容一實施例之電子系統之立體圖;第2B圖為第2A圖的電子系統中,虛框部份沿D方向之側視圖;第3圖為本揭示內容一實施例中,根據第2A圖之電子系統所包含之電子模組建立之電熱π
網路模型;第4圖為本揭示內容一實施例中,根據第2圖之電子系統中的各電子模組以及第3圖之電熱π
網路模型建立之等效電熱電路模型;第5圖為本揭示內容一實施例中,特定體積之示意圖;第6圖為本揭示內容一實施例中,特定體積之示意圖;第7圖為一實施例中,單一垂直共振腔面射型雷射啟動時,電子系統之等溫表面之示意圖;第8A圖為本揭示內容一實施例中,忽略空氣及矽基光學平台之熱阻的情形下,僅啟動單一垂直共振腔面射型雷射時的等效電熱電路模型;以及第8B圖為本揭示內容一實施例中的模擬結果示意圖。
101-106...步驟
Claims (14)
- 一種熱分析模型建立方法,包含下列步驟:根據一電子系統所包含之複數電子模組建立一電熱(electrothermal)π 網路模型,以定義該電子系統之一熱源、一熱傳遞路徑以及一共基恆溫散熱面;定義該電子系統中之複數個觀察點,其中各該等觀察點位於相對一參考點之一外圍空間之一等溫表面上,其中該參考點為該熱源或該等觀察點其中之一;根據該參考點發散至各該等觀察點之一熱功率密度函數、該參考點至各該等觀察點之一熱導係數及該參考點至該觀察點之一距離向量,計算各該等觀察點與該參考點間之一熱傳導溫度差或一熱對流溫度差;以及根據該共基恆溫散熱面之一定義溫度,由各該等觀察點與該參考點間之該熱傳導溫度差或該熱對流溫度差建立該電子系統之一溫度分佈,其中當該參考點為該熱源,係藉由該熱功率密度函數、該熱導係數及該距離向量計算該熱傳導溫度差,其中計算該熱傳導溫度差更包含藉由該熱功率密度函數、該熱導係數及該距離向量對該觀察點所位於該參考點對應之該外圍空間進行一體積分以及位於該參考點對應之該等溫表面進行一面積分計算,其中該熱功率密度函數包含一體積熱功率密度函數qv以及一表面熱功率密度函數qs,該熱導係數包含一體積熱導係數kv以及一表面熱導係數ks,該觀察點之位置為r,該熱源之位置為r`,該外圍空間之該等溫表面之一單位法 向量為,該表面為a,該熱傳導溫度差係表示為:
- 如請求項1所述之熱分析模型建立方法,其中當該參考點為該等觀察點其中之一,係藉由該熱功率密度函數、該熱導係數及該距離向量計算該熱對流溫度差。
- 如請求項2所述之熱分析模型建立方法,其中計算該熱對流溫度差更包含藉由該熱功率密度函數、該熱導係數及該距離向量對該觀察點所位於對應該參考點之該等溫表面進行一面積分計算。
- 如請求項3所述之熱分析模型建立方法,其中該熱功率密度函數包含一表面熱功率密度函數qs,該熱導係數包含一表面熱導係數ks,該觀察點之位置為r,該熱源之位置為r`,該外圍空間之該等溫表面之一單位法向量為,該表面為a,該熱傳導溫度差係表示為:
- 如請求項1所述之熱分析模型建立方法,其中該參考點至各該等觀察點係為該熱傳遞路徑之一部份,且包含該等電子模組至少其中之一,該熱導係數係與該等電子模組至少其中之一之材質相關。
- 如請求項1所述之熱分析模型建立方法,更包含下列步驟:根據該等電子模組以及該電熱π 網路模型建立一等效電熱電路模型,該等效電熱電路模型包含複數熱阻;以及根據該熱功率密度函數、該熱導係數及該距離向量,計算該等效電熱電路模型之該等熱阻。
- 如請求項1所述之熱分析模型建立方法,其中該熱源為該電子系統之最高溫處,該共基恆溫散熱面為該電子系統之最低溫處。
- 一種非暫時性電腦可讀取紀錄媒體,用以儲存一應用程式,該應用程式使一電腦裝置執行一種熱分析模型建立方法,包含下列步驟:根據一電子系統所包含之複數電子模組建立一電熱(electrothermal)π 網路模型,以定義該電子系統之一熱源、一熱傳遞路徑以及一共基恆溫散熱面;定義該電子系統中之複數個觀察點,其中各該等觀察點位於相對一參考點之一外圍空間之一等溫表面上,其中該參考點為該熱源或該等觀察點其中之一;根據該參考點發散至各該等觀察點一熱功率密度函數、該參考點至各該等觀察點之一熱導係數及該參考點至該觀察點之一距離向量,計算各該等觀察點與該參考點間之一熱傳導溫度差或一熱對流溫度差;以及 根據該共基恆溫散熱面之一定義溫度,由各該等觀察點與該參考點間之該熱傳導溫度差或該熱對流溫度差建立該電子系統之一溫度分佈,其中當該參考點為該熱源,係藉由該熱功率密度函數、該熱導係數及該距離向量計算該熱傳導溫度差,其中計算該熱傳導溫度差更包含藉由該熱功率密度函數、該熱導係數及該距離向量對該觀察點所位於該參考點對應之該外圍空間進行一體積分以及位於該參考點對應之該等溫表面進行一面積分計算,其中該熱功率密度函數包含一體積熱功率密度函數qv以及一表面熱功率密度函數qs,該熱導係數包含一體積熱導係數kv以及一表面熱導係數ks,該觀察點之位置為r,該熱源之位置為r`,該外圍空間之該等溫表面之一單位法向量為,該表面為a,該熱傳導溫度差係表示為:
- 如請求項8所述之非暫時性電腦可讀取紀錄媒體,其中當該參考點為該等觀察點其中之一,係藉由該熱功率密度函數、該熱導係數及該距離向量計算該熱對流溫度差。
- 如請求項9所述之非暫時性電腦可讀取紀錄媒體,其中計算該熱對流溫度差更包含藉由該熱功率密度函數、該熱導係數及該距離向量對該觀察點所位於對應該參 考點之該等溫表面進行一面積分計算。
- 如請求項10所述之非暫時性電腦可讀取紀錄媒體,其中該熱功率密度函數包含一表面熱功率密度函數qs,該熱導係數包含一表面熱導係數ks,該觀察點之位置為r,該熱源之位置為r`,該外圍空間之該等溫表面之一單位法向量為,該表面為a,該熱傳導溫度差係表示為:
- 如請求項8所述之非暫時性電腦可讀取紀錄媒體,其中該參考點至各該等觀察點係為該熱傳遞路徑之一部份,且包含該等電子模組至少其中之一,該熱導係數係與該等電子模組至少其中之一之材質相關。
- 如請求項8所述之非暫時性電腦可讀取紀錄媒體,該熱分析模型建立方法更包含下列步驟:根據該等電子模組以及該電熱π 網路模型建立一等效電熱電路模型,該等效電熱電路模型包含複數熱阻;以及根據該熱功率密度函數、該熱導係數及該距離向量,計算該等效電熱電路模型之該等熱阻。
- 如請求項8所述之非暫時性電腦可讀取紀錄媒體,其中該熱源為該電子系統之最高溫處,該共基恆溫散熱面為該電子系統之最低溫處。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW100120735A TWI464603B (zh) | 2011-06-14 | 2011-06-14 | 熱分析模型建立方法及其非暫時性電腦可讀取紀錄媒體 |
US13/204,945 US8543361B2 (en) | 2011-06-14 | 2011-08-08 | Method and non-transitory computer readable medium thereof for thermal analysis modeling |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW100120735A TWI464603B (zh) | 2011-06-14 | 2011-06-14 | 熱分析模型建立方法及其非暫時性電腦可讀取紀錄媒體 |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201250488A TW201250488A (en) | 2012-12-16 |
TWI464603B true TWI464603B (zh) | 2014-12-11 |
Family
ID=47354369
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW100120735A TWI464603B (zh) | 2011-06-14 | 2011-06-14 | 熱分析模型建立方法及其非暫時性電腦可讀取紀錄媒體 |
Country Status (2)
Country | Link |
---|---|
US (1) | US8543361B2 (zh) |
TW (1) | TWI464603B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9740804B2 (en) * | 2014-11-03 | 2017-08-22 | Mentor Graphics Corporation | Chip-scale electrothermal analysis |
EP3203250B1 (en) * | 2016-02-03 | 2023-05-24 | Mitsubishi Electric R&D Centre Europe B.V. | Method and device for estimating a level of damage or a lifetime expectation of a power semiconductor module |
US9990454B2 (en) | 2016-06-03 | 2018-06-05 | International Business Machines Corporation | Early analysis and mitigation of self-heating in design flows |
US11691343B2 (en) | 2016-06-29 | 2023-07-04 | Velo3D, Inc. | Three-dimensional printing and three-dimensional printers |
US9910947B1 (en) * | 2016-06-30 | 2018-03-06 | Cadence Design Systems, Inc. | Methods, systems, and articles of manufacture for implementing an electronic design with solid-fluid analysis driven techniques |
CN106449453B (zh) * | 2016-09-29 | 2019-02-22 | 通富微电子股份有限公司 | 一种半导体封装体的自然对流换热系数及热阻的侦测方法 |
US20200118012A1 (en) * | 2017-04-18 | 2020-04-16 | Hewlett-Packard Development Company, L.P. | Monitoring the Thermal Health of an Electronic Device |
CN109187630A (zh) * | 2018-09-11 | 2019-01-11 | 海南中航特玻科技有限公司 | 一种复合材料导热系数的预测方法 |
TWI685417B (zh) * | 2019-03-14 | 2020-02-21 | 國立虎尾科技大學 | 以熱阻抗迴路設計散熱片配置於石膏模仁之方法 |
CN110362892B (zh) * | 2019-06-28 | 2022-09-27 | 南京理工大学 | 相变型自然冷却散热器肋片最佳尺寸参数确定方法 |
CN113032984B (zh) * | 2021-03-10 | 2022-04-29 | 浙江大学 | 一种基于傅里叶级数的功率模块热建模方法 |
CN115950916B (zh) * | 2023-03-14 | 2023-05-26 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种物体表面热流密度检测方法、装置以及设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW200306009A (en) * | 2002-01-28 | 2003-11-01 | Semiconductor Energy Lab | Semiconductor device and method of manufacturing the same |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7401304B2 (en) * | 2004-01-28 | 2008-07-15 | Gradient Design Automation Inc. | Method and apparatus for thermal modeling and analysis of semiconductor chip designs |
-
2011
- 2011-06-14 TW TW100120735A patent/TWI464603B/zh not_active IP Right Cessation
- 2011-08-08 US US13/204,945 patent/US8543361B2/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW200306009A (en) * | 2002-01-28 | 2003-11-01 | Semiconductor Energy Lab | Semiconductor device and method of manufacturing the same |
Non-Patent Citations (2)
Title |
---|
"Thermal Modeling, Analysis, and Management in VLSI Circuits: Principles and Methods " I "High-Efficiency Green Function-Based Thermal Simulation Algorithms " IEEE 2007 * |
"Thermal Modeling, Analysis, and Management in VLSI Circuits: Principles and Methods " IEEE2006 * |
Also Published As
Publication number | Publication date |
---|---|
US20120323539A1 (en) | 2012-12-20 |
TW201250488A (en) | 2012-12-16 |
US8543361B2 (en) | 2013-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI464603B (zh) | 熱分析模型建立方法及其非暫時性電腦可讀取紀錄媒體 | |
Choi et al. | A new CPU cooler design based on an active cooling heatsink combined with heat pipes | |
Wu et al. | Design and optimization of 3D printed air-cooled heat sinks based on genetic algorithms | |
Culham et al. | Thermal characterization of electronic packages using a three-dimensional Fourier series solution | |
Deng et al. | Experimental study on cooling of high-power laser diode arrays using hybrid microchannel and slot jet array heat sink | |
Choobineh et al. | Experimental and numerical investigation of interdie thermal resistance in three-dimensional integrated circuits | |
US20240028801A1 (en) | Systems and methods for designing a module semiconductor product | |
Tamdogan et al. | Natural convection immersion cooling with enhanced optical performance of light-emitting diode systems | |
Kaplan et al. | Fast thermal modeling of liquid, thermoelectric, and hybrid cooling | |
Solovitz et al. | Micro-channel thermal management of high power devices | |
Grimes et al. | Air flow and heat transfer in fan cooled electronic systems | |
Barabadi et al. | Multiscale transient thermal analysis of microelectronics | |
Zahn | Evaluating thermal characterization accuracy using CFD codes-a package level benchmark study of IcePak/sup TM/and Flotherm/sup R | |
Mohan et al. | Thermal analysis of CPU with CCC and copper base plate heat sinks using CFD | |
KR20210091107A (ko) | 모터제어기의 수냉식 냉각장치의 냉각성능 계산 방법 | |
Shidore et al. | A comparative study of the performance of compact model topologies and their implementation in CFD for a plastic ball grid array package | |
Liu et al. | Analytical modeling for prediction of chip package-level thermal performance | |
Khan et al. | Transient heat transfer in a heat‐generating fin with radiation and convection with temperature‐dependent heat transfer coefficient | |
Yang et al. | Thermal-stress analysis of 3D-IC based on artificial neural network | |
Gess et al. | Single-and two-phase particle image Velocimetry characterization of fluid flow within a liquid immersion cooled electronics module | |
Cai et al. | An internal penalty discontinuous Galerkin method for simulating a thermoelectric cooler | |
Santos et al. | Thermal modeling methodology for efficient system-level thermal analysis | |
Wits et al. | Understanding Thermal Interface Conditions of Chip-on-Board LED Light Sources | |
Regmi et al. | Heat Dissipation Optimization for the Electronic Device in Enclosure Condition | |
Nakayama | Study on heat conduction in a simulated multicore processor chip—Part II: Case studies |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | Annulment or lapse of patent due to non-payment of fees |