源項的案例
CFX中質量源項的加載
在電磁流計算中，往往需要將電磁場的計算結果電磁力的分布加載在流場計算域內；或者需要單獨加載某區(qū)域內的質量源。首先需要將電磁力處理成CFX用戶程序可識別的CSV格式文件，導入到CFX中，然后設置子域動量源，數(shù)據(jù)會按照插分方法分布到各個位置上。
以下是CFX可識別文件的格式。
[Name]
djz
[Spatial Fields]
X, Y, Z
[Data]
X [ m ], Y [ m ], Z [ m ], Force X [ N/m^3 ], Force Y [ N/m^3 ], Force Z [ N/m^3 ]
,,,,,
,,,,,
這些數(shù)據(jù)CFX軟件通過采用“點云圖”的方法來實現(xiàn)插分。整個過程包括對被插分點最近的三個原始數(shù)據(jù)點的快速定位，以及根據(jù)它們離被插分點距離遠近的一個反向加權的平均過程。在求解過程當中，根據(jù)離散和數(shù)值積分過程的具體要求，求解器需要各個不同位置上的值，比如在積分點上，節(jié)點上和各個面的中心點上。在所有的情況下，這個需要的位置都會被確定，原始的數(shù)據(jù)就會被插分到該位置。
展開 施工安全危險源426項,史上最全!
以下是施工安全危險源清單，共426項，請收好！
FEKO中激勵源與求解項設置細節(jié)（持續(xù)完善）
feko中可以設置多個求解進程configuration
feko中的求解項分為兩種，一種是standralconfiguration，主要是求解方向圖，RCS，傳輸系數(shù)等電磁參數(shù)；還有一種是S-parameter configuration，主要求解端口網(wǎng)絡/陣列天線等端口之間的S參數(shù)。天線方向圖farfield與多端口的S參數(shù)的求解不能同時放在一個configuration中，需要添加一個configuration，但是不同的congfiguration可以放在一個CFX文件中（只要模型一致）。
3. configuration之間的差異有幾種選擇
frequency per config 每個求解項的頻率不同;
source per config每個求解項的激勵設置不同（幅度，相位）;
load per config（每個求解項的負載不同）;
power per config（每個求解項的功率不同）
4. 對稱面的設置
理想電面和理想磁面的使用可以有效的降低未知量數(shù)，結構對稱的面的使用則可以簡化建模。理想電面和理想磁面的使用依據(jù)，主要看電場，電流或磁場和磁流的方向。電場，電流的方向與理想電面垂直，與理想磁面平行；磁場和磁流的方向垂直于理想磁面，與理想電面平行。
展開 四十四、Fluent 收斂標準-質量和能量守恒
需要先將內部邊界類型更改為wall(用于定位面積向量)，不用進行任何操作，然后再更改為interior，計算Fluxes</p><p> </p><p><br></p><p><span style="color: rgb(255, 129, 36);">如果不涉及質量源項，那么模型進出口的mass flow rate應該相等，凈值Net Results應該等于0；</span></p><p><span style="color: rgb(255, 129, 36);">如果不涉及能量源項，模型進出口及壁面處的Total Heat Transfer Rate的凈值Net Results應該等于0</span></p><p> </p><p><strong>2.3 涉及源項問題</strong></p><p><br></p><p>如果在Cell Zone Conditions---Fluid中設置了質量源項（kg/(m3*s)），則Net results顯示的為添加質量源項計算域的積分值。即每個網(wǎng)格的質量源項與計算域網(wǎng)格體積相乘之和，單位恰好為kg/s。</p><p> <img src="https://mmbiz.qpic.cn/mmbiz_png/8tJMdLVYZy9aicLdtZXwZOLDDLnP3ic2M2L1nBwOvOXJw9DAYOutbTKT7eWiaRNDvxlOs43CrnQpcYeHSun0iauEKw/640?wx_fmt=png" width="100%"> </p><p><br></p><p>同理如果在Cell Zone Conditions---Fluid中設置了能量源項（W/m3），則Net results顯示的為添加能量源項計算域的積分值。
展開 
湍流-化學作用的噴霧燃燒模擬 | 基于OpenFOAM的FGM模型實現(xiàn)與分析
在3D求解器中，從FGM表差值得到的最重要的變量是進度變量源項，它決定了反應進程和點火延遲時間。
圖 6 三個算例中進度變量源項變化與間的關系
盡管在βz–βc和βz–δc情況之間，整體點火特性（點火延遲時間）差異很小，但差異在反應進度變量源項的時間演變中更為明顯（圖 7（a））。圖 7（b）中顯示了提取進度變量源項峰值的單元格中的混合分數(shù)及其時間演變。圖 8給出了三個階段中進度變量源項與進度變量的關系。
圖 7 （a）反應進度變量源項最大值隨時間演變過程；（b）提取進度變量源項峰值的單元格中的混合分數(shù)及時間演變
圖 8 三個階段中進度變量源項與進度變量的關系
3.4 TCI對火焰結構的影響
本研究發(fā)現(xiàn)，在穩(wěn)態(tài)條件下，混合分數(shù)和進度變量的變化也會改變火焰的結構。圖 11顯示了從βz–βc，δz–βc和βz–δc情況獲得的結果在1.5 ms ASI時計算出的OH質量分數(shù)分布。
圖 11 三個階段中進度變量源項與進度變量的關系
如圖 11所示，在βz–βc和δz–βc案例之間，OH質量分數(shù)分布存在顯著差異。與δ-PDF相比，Z-的β-PDF的OH質量分數(shù)的峰值要低得多。另一方面，考慮到混合比例的變化，β-PDF預測了空間中更多的OH質量分數(shù)分布。而PDF的火焰相對較薄。在使用mRIF，F(xiàn)GM，TFM和TPDF模型的其他ECN噴霧（正庚烷或正十二烷）燃燒模型中也報告了相同的發(fā)現(xiàn)，這歸因于這些研究中的湍流波動。至于進度變量方差的影響，圖 11中的OH質量分數(shù)分布表明，在βz–βc情況和βz–δc情況下，火焰結構和OH質量分數(shù)最大值相似。也就是說，進度變量的方差對火焰結構的影響很小。然而，從圖 11(c)中可以明顯看出，它使火焰浮起長度變短。
展開 基于部分浸潤效應的歐拉壁膜流動形態(tài)演變模擬仿真
液滴與壁面相互作用決策圖
歐拉壁膜的質量、動量、能量守恒
質量守恒
等式左邊：非穩(wěn)態(tài)項和對流項；
等式右邊：單位面積下的質量源項，如液滴收集、壁膜分離、壁膜脫落、相變等行為下，需更新壁膜質量源項。
動量守恒
等式左邊：非穩(wěn)態(tài)項和對流項
等式右邊：
l 第一項：擴散項：其中，壓力=氣流壓力+垂直于壁膜的重力分量+表面張力
l 第二項：重力源項：平行于壁膜的重力分量
l 第三&四項：凈粘性切應力源項：氣流-壁膜之間&壁面與壁膜之間
l 第五項：動量源項：液滴的收集與分離
l 第六項：表面力源項：壁膜的表面張力、壁膜與壁面間的接觸角
能量守恒
等式左邊：非穩(wěn)態(tài)項和對流項
等式右邊：
l 第一項：氣流-壁膜間與壁膜-壁面間的凈熱通量
l 第二項：液體從氣流中撞擊到壁面的能量源項
l 第三項：相變能量源項：是蒸發(fā)/冷凝的質量變化率，是相變潛熱值
控制方程的應用與求解
l 質量守恒和動量守恒：歐拉壁膜模型的核心基礎
l 能量守恒：涉及壁膜的傳熱需要啟用
l 由于歐拉壁膜模型適用的液膜很薄，故采用潤滑近似(平行流動)，因此這些控制方程是在平行于壁面的局部當?shù)刈鴺酥星蠼獾摹?壁膜速度的獲取
壁膜速度可以假設為兩個組成：外部流動產(chǎn)生的切應力驅動壁膜速度、由于重力產(chǎn)生的重力驅動壁膜速度。
切應力驅動-速度線性變化
由氣液交界面處外部氣體流動側，切應力與壁膜側切應力平衡，得出外部流動產(chǎn)生的切應力驅動壁膜速度：
▲ 圖3.
展開 fluent中的湍流阻尼
添加以下源項到ω方程：
其中，
A_i是i相的界面面積密度；
Δn是單元到界面的法向高度；
β是封閉k-ω模型破壞系數(shù)項,等于0.075；
B是阻尼因子；
μ_i是相i的粘度；
ρ_i是相i的密度。
相i的界面面積密度計算為：
其中，
α_i是相i的體積分數(shù)；
|Δα_i|是體積分數(shù)梯度的大小。
網(wǎng)格大小Δn是使用網(wǎng)格信息在內部計算的。您可以在粘性模型對話框中指定阻尼因子B。阻尼因子的默認值為10。
湍流阻尼是可用的對VOF和混合模型。注意，當使用非混合相流體模型時，它也適用于歐拉多相流模型。
如果啟用了歐拉多相模型，則可以指定湍流多相模型。如果每相都使用湍流模型,那么ω方程的源項添加到每個相。如果啟用了VOF或混合模型,或歐拉多相模型與混合湍流模型,這時所有相求和作為源項添加到混合水平的ω方程。
展開 煉焦過程中VOCs無組織排放的計算方法探討
4.1 世界范圍內的煉焦過程VOCs 無組織排放量的計算方法仍存在較大差異，計算結果存在數(shù)量級差距；歐洲環(huán)境署排放系數(shù)法存在較多的源項遺漏；美國國家環(huán)保局排放系數(shù)法、公式法分類較細，但未將推焦環(huán)節(jié)納入估算；上海地區(qū)排放系數(shù)法對推焦環(huán)節(jié)VOCs 無組織排放占整體比例的估計值達到了51%，與美國國家環(huán)保局有一定差異。
4.2 值得注意的是，對于年產(chǎn)250 萬t 產(chǎn)能級別的焦爐組，即使按照美國國家環(huán)保局危險大氣污染物國家排放標準（NESHAP）的最優(yōu)控制措施運行，其VOCs年無組織排放量估計值仍然達到了40 t 級，十分可觀。對于控制措施一般情形下的估計值則達到了150 t 級。因此，煉焦過程的VOCs 無組織排放不容小覷。
4.3 煉焦過程的VOCs 無組織排放源項、過程都很復雜，有毒有害物質多、排放量大，需要引起重視；但目前國內外的無組織排放量計算方法仍存在較大差異，計算結果存在數(shù)量級差距。后續(xù)研究可以以具體排放源項和過程為切入點，對排放量計算方法進一步細化研究。
精彩推薦，點擊進入 ！
◆28歲當縣長，36歲升至省長，一路驚心動魄九死一生
◆權色仕途：走近女領導
◆從鄉(xiāng)鎮(zhèn)到省委的官場筆記，讀懂受益無窮！
?本文適用本平臺“免責聲明”請回復“免責聲明”查詢
▼點擊“閱讀原文”進入“優(yōu)秀小說選讀頻道"
展開 CFD理論|網(wǎng)格體積的求解
背景
在CFD計算中，源項的求解需要用到網(wǎng)格體積，比如在湍動能的輸運方程中：
將湍動能在網(wǎng)格進行積分，可以得到：
篇幅所限，這里只討論最后一項源項，在使用有限體積法求解輸運方程時，首先需要對控制體積上的方程進行積分，由于積分是滿足加減法運算，因此可以對每一項單獨積分再相加，因此可以得到源項的對控制體積的積分：
因此，源項的體積積分可以等于網(wǎng)格中心的值乘上網(wǎng)格體積：
燃油噴射器潛在空蝕的數(shù)值預測
對于氣相，該等式具有以下形式：
方程（1）
其中α是汽相的體積分數(shù)，是蒸汽密度，是氣相的速度矢量，S是凈質量源項。
方程（2）
其中，μ是粘度混合物，ρ是混合物的密度，是蒸汽密度，是液體密度，p是壓力。
2 空化模型
方程（1）中的質量源項S可以寫為：
方程（3）
方程（4）
蒸汽體積分數(shù)與氣泡數(shù)密度和氣泡半徑有關，方程（5）：
根據(jù)Rayleigh-Plesset方程忽略二階項和表面張力，氣泡半徑的變化率為方程（6）：
由方程（3）至（6），質量源項可以改寫為方程（7）：
蒸發(fā)過程涉及蒸汽形成和冷凝過程，涉及蒸汽坍塌、質量源項分別定義：
時，蒸發(fā)源項為方程（8）：
時，縮合源項為方程（9）：
上述空化模型未考慮溶解氣體在液相中的影響。
3 流體屬性
液相和氣相都被建模為可壓縮流體。
展開 如何在 COMSOL 中進行粒子計數(shù)
累加基于：設為單元時，累加變量將只是計算給定時間點駐留在單元內的所有粒子的源項之和。設為單元和時間時，粒子會基于其在單元駐留時間的長短對經(jīng)過的單元產(chǎn)生貢獻。
源項：這是在粒子上定義并希望投射到底層網(wǎng)格中的表達式。當進行粒子計數(shù)時，只需將源項設為 “1”，但也支持粒子上的任意表達式，比如電荷或動能；還可以基于在粒子駐留域內定義的變量。
單位：當為累積變量選擇單位時，源項所需的單位將相應變化。
為了計數(shù)粒子總數(shù)，您可以在累加器的駐留域中增加積分組件耦合。邊界累加器會自動在選定邊界上增加組件耦合。在我們的示例中，粒子總數(shù)由 <integration_operator_name>((pt.count))得到。這可以通過全局計算節(jié)點計算。每個網(wǎng)格單元內的粒子數(shù)也可以于其他物理場耦合，因為它也是一個自由度。我們可以在累加變量和底層網(wǎng)格繪圖中繪制粒子的位置，以此顯示軟件如何進行粒子計數(shù)。
繪圖顯示了粒子在底層網(wǎng)格（灰線）上的位置（黑點）。每個單元內顏色表示累加變量的值。
從上圖中可以清楚看到，累加器的確會統(tǒng)計每個網(wǎng)格單元內的粒子數(shù)。對不包含粒子的網(wǎng)格單元，累加變量為零，網(wǎng)格單元顯示為藍色。大部分網(wǎng)格單元中包含一個粒子，這類網(wǎng)格顯示為綠色。其中一個網(wǎng)格內恰好包含兩個粒子，顯示為紅色。
您還可以使用累加器來統(tǒng)計經(jīng)過一條內部邊界的粒子數(shù)。為此，您只需在粒子經(jīng)過的邊界上增加一個壁條件，并設為穿過。向壁節(jié)點增加累加器子特征，并設為：
當粒子經(jīng)過邊界時，累加器會增加對應邊界網(wǎng)格單元中的自由度；得到通過內部邊界的粒子數(shù)的空間分布，如下方動畫所示。
我們可以輕松以時間函數(shù)的形式繪制經(jīng)過邊界的粒子總數(shù)；只需增加一維繪圖組和全局繪圖特征。累加器會創(chuàng)建多個預定義變量，并加成到所有網(wǎng)格單元的累加變量中。
展開 
k-ε模型中的K和ε物理意義
在求解標準k-e雙方程湍流模型時（采用渦粘假設，求湍流粘性系數(shù)，然后和N－S方程耦
合求解粘性流場），發(fā)現(xiàn)湍動能產(chǎn)生項（雷諾應力和一個速度張量相乘組成的項）出現(xiàn)負
值，請問是不是一種錯誤現(xiàn)象？如果是錯誤現(xiàn)象一般怎樣避免。另外處理湍動能產(chǎn)生項采
用什么樣的差分格式最好。而且因為源項的影響，使得程序總是不穩(wěn)定，造成k,e值出現(xiàn)負
值，請問有什么辦法克服這種現(xiàn)象。
你可以試試這里計算的時候加一個判斷，出現(xiàn)負值的時候強制為一個很小的正值。
這可能是因為你采用的數(shù)值格式的問題，一般計算程序對k方程都要做一定處理，
以保證k的正定。比如，強制規(guī)定源項與0的關系，以使數(shù)值計算穩(wěn)定。
就ke模型而言。
它是problem dependent.對簡單的無彎曲無旋轉無...的湍流問題，它能算而且能給出好的結果，但對復雜的流動問題，它就不能使用了。
出現(xiàn)負的ke不僅僅是計算格式的問題，
更重要的是模型問題，
沒有誰能證明ke模型在任何流動問題中都能保證ke是正的。
有這么一些辦法避免ke出現(xiàn)負值
1。對K＝ln(k)和E=ln(e)求解，問題：壁面ke=0難處理，
2。先用層流計算500步，然后再用ke算，
3。各種強制限制辦法
4。源項局部線性化。
5。算到一定程度，如果k值趨勢對了，就干脆不求ke方程了。
可以參考：《湍流的計算模型》 陳義良 1991 中國科技大學出版社
展開 PHOENICS程序應用-理論基礎部分
壓力梯度項從源項中分離出來。例如對ue的控制容積:
?(pp-pe)Dy
這里假設在ue的控制容積的東、西界面上壓力是各自均勻的，分別為pE、pp。于是關于ue的離散方程具有以下形式:
aeue=?anbunb+b+(pp-pe)Ae
類似地，對vn的控制容積作積分可得:
anvn=?anbvnb+b+(pp-pN)An 3.1.3 計算方法 3.1.3.1 SIMPLE算法的計算步驟 采用SIMPLE算法實施關于u、v、p代數(shù)方程的分離式求解時，計算步驟如下: (1) 假定一個速度分布，記為u0，v0，以次計算動量離散方程的系數(shù)及常數(shù)項;
(2) 假定一個壓力場p*;
(3) 依次求解兩個動量方程，得u*、v* ;
(4) 求解壓力修正值方程，得p’ ;
(5) 據(jù)p’改進速度值 ;
(6) 利用改進后的速度場求解那些通過源項物性等與速度場耦合的F變量。如果F并不影響流場，則應在速度場收斂后再求解 ;
(7) 利用改進后的速度場重新計算動量離散方程的系數(shù)，并用改進后的壓力場作為下一層次迭代計算的初值。重復上述步驟，直到獲得收斂的解。 PHOENICS程序計算方法采用的是SIMPLEST算法，與SIMPLE相比，它主要有以下兩個特點: (1) 對流項采用迎風格式，因為這是一個絕對穩(wěn)定的格式，且擴散項與對流項的影響系數(shù)可以分離開來，不象指數(shù)(或乘方)格式那樣綜合在一起，至于由迎風差分所引起的假擴散問題，則采用逐步加密網(wǎng)格、以獲得與網(wǎng)格稀密程度無關的解這種做法加以克服。
展開 VirtualFlow | 熱管相變換熱仿真，支持不同尺度的氣液兩相相變計算
蒸發(fā)和冷凝過程中的相變通過UDF在體積分數(shù)方程、能量方程和組分輸運方程中分別添加質量源項、能量源項和相等的質量源項實現(xiàn)。
這種算法能夠精確地模擬吸液芯的毛細現(xiàn)象、蒸發(fā)管的沸騰、冷凝器的冷凝等復雜現(xiàn)象，為熱管的設計與優(yōu)化提供了堅實的技術支持。
（二）準確、可靠的計算結果
在實際案例中，VirtualFlow軟件展現(xiàn)了優(yōu)秀的計算精度和可靠性。以某物理研究所的環(huán)路熱管項目為例，在50W功率下2D軸對稱條件下，蒸發(fā)器內的流場最終達到穩(wěn)態(tài)，其液相體積分數(shù)、相變速率、液體/氣體總體積、質量流量等參數(shù)的計算結果與實驗趨勢高度一致。
蒸發(fā)器壁溫計算結果與測量結果的偏差基本控制在1.5℃以內，冷凝器部件仿真結果同樣表現(xiàn)出色。整機仿真結果顯示，隨著熱流密度的增高，冷凝器中的液體體積先減后增，這一結果與實際物理現(xiàn)象相符，充分驗證了軟件的準確性和可靠性。
蒸發(fā)器內各統(tǒng)計量隨時間的變化：a)相變速率；b)液體/氣體總體積；c)質量流量
（三）解決方案優(yōu)勢
VirtualFlow軟件在熱管領域的應用具有諸多獨特優(yōu)勢。
它具備氣液兩相模型，能夠模擬微納米尺度如空隙尺度的多孔介質、微納結構等吸液芯的毛細潤濕和蒸發(fā)過程，預測毛細能力及蒸發(fā)換熱性能。
支持在微通道納米尺度中計算兩相相變，可用于表面凝結和核態(tài)沸騰的相變過程計算，以及計算在相變過程中的換熱情況。
軟件支持熱限制模型與RPI壁面沸騰模型，并開發(fā)有先進的的壁面冷凝模型，可根據(jù)此對池沸騰、大空間冷凝相變、壁面相變等進行數(shù)值模擬。
在處理熱虹吸問題時，通過模擬蒸發(fā)相變，觸發(fā)熱虹吸效應，進而研究熱邊界及固體結構對虹吸過程流量、流速的影響。
軟件能夠根據(jù)計算的兩相流動狀態(tài)自動切換所采用的兩相流模型，適用的多相流典型形態(tài)包括界面流、離散相以及混合流，提升計算準確性。
展開 使用用戶自定義標量（UDS）計算流場中的電場強度 ￥20
UDS案例 msh 文件 cas 文件
概述
UDS：用戶自定義標量
操作
Define-> user define -> Scalar
設置種類
對流項
時間項
擴散系數(shù)
源項
案例
二維槽道內電場求解
案例描述
二維槽道內，長L,入口水流v，在槽道兩端加電壓V，入口為Vin，出口為Vout，求解整個槽道內的電壓分布和電場強度
幾何模型和物理模型
UDS 設置
Define-> user define -> Scalar
Number of UDS ：1
沒有對流項和時間項
沒有源項
Define->material
定義擴散系數(shù) 為1
入口設置UDS 為 Vin
出口設置UDS 為Vout
求解
展開 