安全模擬與仿真概述安全模擬與仿真

模擬與仿真的概念所謂仿真就是建立系統(tǒng)的模型（數(shù)學模型、物理效應(yīng)模型或數(shù)學-物理效應(yīng)模型），并在模型上進行實驗和研究一個存在的或設(shè)計中的系統(tǒng)。模擬，即是外形仿真、操作仿真、視覺感受仿真，使用真實的汽車模型或其他等比例的飛機、飛船等模型作為參與者的操控平臺，利用VR技術(shù)（虛擬現(xiàn)實技術(shù)），通過實際操作，使參與者有身臨其境的切身體會。安全專業(yè)有哪些方面涉及到模擬與仿真？《計算機模擬與仿真在XXX方面的應(yīng)用進展》小論文11月7日提交仿真模擬技術(shù)的三大組成部分對一個工程技術(shù)系統(tǒng)進行模擬仿真，包括了建立模型、實驗求解和結(jié)果分析三個主要步驟。幾何模型數(shù)學－物理模型數(shù)值計算的軟件

計算流體力學CFD(1)引言流體力學的三種研究方法流體力學的控制方程組基本物理學原理基本物理學原理流體力學基本控制方程連續(xù)性方程質(zhì)量守恒定律動量方程牛頓第二定律能量方程能量守恒定律流動模型流動模型1）有限控制體模型對于有連續(xù)性的流體，有下面兩種模型：2）無窮小流體微團我們不是同時觀察整個流場，而是將物理學基本原理用在這些流動模型上，從而得到流體流動方程。流動模型有限控制體模型空間位置固定的有限控制體，流體流過控制體隨流體運動的有限控制體，同一批流體質(zhì)點始終位于同一控制體內(nèi)流動模型無窮小流體微團模型空間位置固定的無窮小流體微團，流體流過微團沿流線運動的無窮小流體微團，其速度等于流線上每一點的當?shù)厮俣任镔|(zhì)導數(shù)（運動流體微團的時間變化率）流動控制方程經(jīng)常用物質(zhì)導數(shù)來表達。物質(zhì)導數(shù)（運動流體微團的時間變化率）沿流線運動的無窮小流體微團，其速度等于流線上每一點的當?shù)厮俣炔捎昧黧w微團模型來理解物質(zhì)導數(shù)的概念：物質(zhì)導數(shù)（運動流體微團的時間變化率）流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖物質(zhì)導數(shù)（運動流體微團的時間變化率）流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖考慮非定常流動：物質(zhì)導數(shù)（運動流體微團的時間變化率）流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖考慮非定常流動：物質(zhì)導數(shù)（運動流體微團的時間變化率）流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖在1點做如下的泰勒級數(shù)展開：物質(zhì)導數(shù)（運動流體微團的時間變化率）流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖物質(zhì)導數(shù)（運動流體微團的時間變化率）流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖這里D

/Dt代表流體微團通過1點時，流體微團密度變化的瞬時時間變化率。我們把D

/Dt定義為密度的物質(zhì)導數(shù)。物質(zhì)導數(shù)（運動流體微團的時間變化率）流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖注意D

/Dt是給定的流體微團在空間運動時，其密度的時間變化率。我們必須跟蹤運動的流體微團，注意它通過點1時密度的變化。物質(zhì)導數(shù)（運動流體微團的時間變化率）流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖物質(zhì)導數(shù)D

/Dt與偏導數(shù)

/

t不同

，

/

t是在固定點1時觀察密度變化的時間變化率，該變化由流場瞬間的起伏所引起。物質(zhì)導數(shù)（運動流體微團的時間變化率）物質(zhì)導數(shù)（運動流體微團的時間變化率）向量算子物質(zhì)導數(shù)（運動流體微團的時間變化率）D/Dt是物質(zhì)導數(shù)，它在物理上是跟蹤一個運動的流體微團的時間變化率；流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖物質(zhì)導數(shù)（運動流體微團的時間變化率）

/

t叫做當?shù)貙?shù)，它在物理上是固定點處的時間變化率；流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖物質(zhì)導數(shù)（運動流體微團的時間變化率）

叫做遷移導數(shù)，它在物理上表示由于流體微團從流場中的一點運動到另一點，流場的空間不均勻性而引起的時間變化率。流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖物質(zhì)導數(shù)（運動流體微團的時間變化率）物質(zhì)導數(shù)可用于任何流場變量，比如Dp/Dt、DT/Dt等流體微團在流場中的運動－物質(zhì)導數(shù)的示意圖物質(zhì)導數(shù)（運動流體微團的時間變化率）人進入山洞，洞內(nèi)溫度比洞外溫度低，正經(jīng)過洞口向里進時，同時被雪球擊中。洞內(nèi)溫度比洞外溫度低所引起的溫降遷移導數(shù)物質(zhì)導數(shù)當?shù)貙?shù)遷移導數(shù)被雪球擊中所引起的溫降當?shù)貙?shù)總的溫降物質(zhì)導數(shù)物質(zhì)導數(shù)（運動流體微團的時間變化率）物質(zhì)導數(shù)全微分：對時間的全導數(shù)：物質(zhì)導數(shù)（運動流體微團的時間變化率）物質(zhì)導數(shù)物質(zhì)導數(shù)在本質(zhì)上與對時間的全導數(shù)相同。對時間的全導數(shù)：速度散度及其物理意義速度散度這一表達式也經(jīng)常出現(xiàn)在流體動力學方程中。隨流體運動的有限控制體，同一批流體質(zhì)點始終位于同一控制體內(nèi)速度散度及其物理意義考慮如圖所示隨流體運動的控制體。這個控制體在運動中，總是由相同的流體粒子組成，因此它的質(zhì)量是固定的，不隨時間變化。隨流體運動的有限控制體，同一批流體質(zhì)點始終位于同一控制體內(nèi)速度散度及其物理意義但是，當它運動到流體不同的區(qū)域，由于密度不同，它的體積和控制面會隨著時間改變。隨流體運動的有限控制體，同一批流體質(zhì)點始終位于同一控制體內(nèi)速度散度及其物理意義也就是說，隨著流場特性的變化，這個質(zhì)量固定的、運動著的控制體，體積不斷地增大或減小，形狀也在不斷地改變著。速度散度及其物理意義速度散度的物理意義：是每單位體積運動著的流體微團，體積相對變化的時間變化率。連續(xù)性方程空間位置固定的有限控制體模型空間位置固定的有限控制體模型空間位置固定的有限控制體模型連續(xù)性方程質(zhì)量守恒定律通過控制面S流出控制體的凈質(zhì)量流量＝控制體內(nèi)質(zhì)量減少的時間變化率空間位置固定的有限控制體模型空間位置固定的有限控制體模型通過控制面S流出控制體的凈質(zhì)量流量＝控制體內(nèi)質(zhì)量減少的時間變化率或空間位置固定的有限控制體模型空間位置固定的有限控制體模型連續(xù)性方程：隨流體運動的有限控制體模型隨流體運動的有限控制體模型隨流體運動的有限控制體模型連續(xù)性方程質(zhì)量守恒定律有限控制體的總質(zhì)量為：隨流體運動的有限控制體模型隨流體運動的有限控制體模型連續(xù)性方程：空間位置固定的無窮小微團模型空間位置固定的無窮小微團模型空間位置固定的無窮小微團模型連續(xù)性方程質(zhì)量守恒定律流出微團的質(zhì)量流量＝微團內(nèi)質(zhì)量的減少空間位置固定的無窮小微團模型空間位置固定的無窮小微團模型X方向的凈流出量為：流出微團的質(zhì)量流量＝微團內(nèi)質(zhì)量的減少空間位置固定的無窮小微團模型空間位置固定的無窮小微團模型Y方向的凈流出量為：流出微團的質(zhì)量流量＝微團內(nèi)質(zhì)量的減少空間位置固定的無窮小微團模型空間位置固定的無窮小微團模型Z方向的凈流出量為：流出微團的質(zhì)量流量＝微團內(nèi)質(zhì)量的減少空間位置固定的無窮小微團模型空間位置固定的無窮小微團模型微團內(nèi)質(zhì)量增加的時間變化率為：流出微團的質(zhì)量流量＝微團內(nèi)質(zhì)量的減少空間位置固定的無窮小微團模型空間位置固定的無窮小微團模型流出微團的質(zhì)量流量＝微團內(nèi)質(zhì)量的減少或空間位置固定的無窮小微團模型空間位置固定的無窮小微團模型或連續(xù)性方程：隨流體運動的無窮小微團模型隨流體運動的無窮小微團模型隨流體運動的無窮小微團模型流體微團的質(zhì)量：連續(xù)性方程質(zhì)量守恒定律隨流體運動的無窮小微團模型隨流體運動的無窮小微團模型連續(xù)性方程質(zhì)量守恒定律方程不同形式之間的轉(zhuǎn)換空間位置固定的有限控制體模型空間位置固定的無窮小微團模型積分形式與微分形式的重要注釋積分形式的方程允許出現(xiàn)間斷，微分形式的方程要求流動參數(shù)是連續(xù)的。因此，積分形式的方程比微分形式的方程更基礎(chǔ)、更重要。在流動包含真實的間斷（如激波）時，這一點尤其重要。動量方程動量方程動量方程牛頓第二定律動量方程力的兩個來源：1）體積力：直接作用在流體微團整個體積微元上的力，而且作用是超距離的，比如重力，電場力，磁場力。隨流體運動的無窮小微團模型動量方程力的兩個來源：2）表面力：直接作用在流體微團的表面。隨流體運動的無窮小微團模型動量方程表面力的兩個來源：1）壓力2）粘性力動量方程粘性力的兩個來源：1）正應(yīng)力2）切應(yīng)力動量方程切應(yīng)力：與流體剪切變形的時間變化率有關(guān)，如下圖中的

xy動量方程正應(yīng)力：與流體微團體積的時間變化率有關(guān)，如下圖中的

xx動量方程作用在單位質(zhì)量流體微團上的體積力記做，其X方向的分量為隨流體運動的無窮小微團模型動量方程作用在流體微團上的體積力的X方向分量＝隨流體運動的無窮小微團模型動量方程作用在流體微團上的X方向的壓力＝動量方程作用在流體微團上的X方向的正應(yīng)力＝動量方程作用在流體微團上的X方向的切應(yīng)力＝動量方程作用在流體微團上的X方向總的表面力＝隨流體運動的無窮小微團模型動量方程作用在流體微團上的X方向總的力：隨流體運動的無窮小微團模型動量方程作用在流體微團上的X方向總的力：動量方程運動流體微團的質(zhì)量：隨流體運動的無窮小微團模型動量方程運動流體微團的X方向的加速度：隨流體運動的無窮小微團模型動量方程由牛頓第二定理得粘性流X方向的動量方程：隨流體運動的無窮小微團模型動量方程類似地，可得Y方向和Z方向的動量方程：動量方程三個方向的動量方程：以上為非守恒形式的納維－斯托克斯方程(Navier-Stokes方程)，簡稱非守恒形式的N－S方程。動量方程非守恒形式的的N－S方程可以轉(zhuǎn)化為如下守恒形式的N－S方程動量方程牛頓流體：流體的切應(yīng)力與應(yīng)變的時間變化率(也就是速度梯度)成正比。在空氣動力學的所有實際問題中，流體都可以看成牛頓流體。動量方程對牛頓流體，有動量方程完整的N－S方程守恒形式：能量方程能量方程隨流體運動的無窮小微團的能量通量能量方程能量守恒定律能量方程隨流體運動的無窮小微團的能量通量流體微團內(nèi)能量的變化率流入微團內(nèi)的凈熱流量＝＋體積力和表面力對微團做功的功率能量方程隨流體運動的無窮小微團的能量通量作用于速度為V的流體微團上的體積力，做功的功率為：能量方程隨流體運動的無窮小微團的能量通量對比下圖作用在面adhe和面bcgf上的壓力，則壓力在X方向上做功的功率為：能量方程隨流體運動的無窮小微團的能量通量類似地，在面abcd和面efgh上，切應(yīng)力在X方向上做功的功率為：能量方程隨流體運動的無窮小微團的能量通量所有表面力（包括壓力、正應(yīng)力、切應(yīng)力）在X方向上做功的功率為：能量方程所有力（包括體積力、表面力）做功的功率總和（包括X方向、Y方向、Z方向）為：能量方程隨流體運動的無窮小微團的能量通量流體微團內(nèi)能量的變化率流入微團內(nèi)的凈熱流量＝＋體積力和表面力對微團做功的功率能量方程隨流體運動的無窮小微團的能量通量流入微團的凈熱流量來源兩個方面：1）體積加熱，如吸收或釋放的熱輻射。能量方程隨流體運動的無窮小微團的能量通量流入微團的凈熱流量來源兩個方面：2）由溫度梯度導致的跨過表面的熱輸運，即熱傳導。能量方程隨流體運動的無窮小微團的能量通量定義為單位質(zhì)量的體積加熱率；運動流體微團的質(zhì)量為，因此，微團的體積加熱為能量方程隨流體運動的無窮小微團的能量通量考慮面adhe和面bcgf，熱傳導在X方向?qū)α黧w微團的加熱為：能量方程隨流體運動的無窮小微團的能量通量熱傳導在X、Y、Z三個方向?qū)α黧w微團的加熱為：能量方程隨流體運動的無窮小微團的能量通量因此，流入微團內(nèi)的凈熱流量為：能量方程根據(jù)傅立葉熱傳導定律，熱傳導產(chǎn)生的熱流與當?shù)氐臏囟忍荻瘸烧?，設(shè)k為熱導率，則能量方程隨流體運動的無窮小微團的能量通量因此，流入微團內(nèi)的凈熱流量可寫為：能量方程隨流體運動的無窮小微團的能量通量流體微團內(nèi)能量的變化率流入微團內(nèi)的凈熱流量＝＋體積力和表面力對微團做功的功率能量方程隨流體運動的無窮小微團的能量通量跟隨流體運動的微團的能量有兩個來源：1）由分子隨機運動而產(chǎn)生的內(nèi)能，定義單位質(zhì)量內(nèi)能為e能量方程隨流體運動的無窮小微團的能量通量跟隨流體運動的微團的能量有兩個來源：2）流體微團平動時具有的動能，單位質(zhì)量的動能為能量方程隨流體運動的無窮小微團的能量通量運動流體微團的質(zhì)量為，因此，流體微團內(nèi)能量的變化率為能量方程隨流體運動的無窮小微團的能量通量流體微團內(nèi)能量的變化率流入微團內(nèi)的凈熱流量＝＋體積力和表面力對微團做功的功率根據(jù)能量守恒定律，有能量方程流體微團內(nèi)能量的變化率流入微團內(nèi)的凈熱流量＝＋體積力和表面力對微團做功的功率于是能量方程（非守恒形式）為：能量方程只用內(nèi)能e表示的能量方程（非守恒形式）為：只用內(nèi)能e表示的能量方程中不包含體積力項。能量方程只用內(nèi)能e表示的能量方程（非守恒形式）可寫為：根據(jù)，，能量方程對牛頓流體，有能量方程只用內(nèi)能e表示的能量方程（非守恒形式）可寫為：能量方程只用內(nèi)能e表示的能量方程（守恒形式）為：能量方程用總能表示的能量方程（守恒形式）為：流體力學控制方程的總結(jié)與注釋粘性流動的納維－斯托克斯(Navier-Stokes)方程粘性流動的納維－斯托克斯(Navier-Stokes)方程非定常三維可壓縮粘性流動的控制方程總結(jié)如下：1.連續(xù)性方程非守恒形式：守恒形式：粘性流動的納維－斯托克斯(Navier-Stokes)方程非定常三維可壓縮粘性流動的控制方程總結(jié)如下：2.動量方程非守恒形式：X方向：Y方向：Z方向：粘性流動的納維－斯托克斯(Navier-Stokes)方程非定常三維可壓縮粘性流動的控制方程總結(jié)如下：2.動量方程守恒形式：X方向：Y方向：Z方向：粘性流動的納維－斯托克斯(Navier-Stokes)方程非定常三維可壓縮粘性流動的控制方程總結(jié)如下：3.能量方程非守恒形式：粘性流動的納維－斯托克斯(Navier-Stokes)方程非定常三維可壓縮粘性流動的控制方程總結(jié)如下：3.能量方程守恒形式：無粘流歐拉(Euler)方程非定常三維可壓縮無粘流動的控制方程總結(jié)如下：1.連續(xù)性方程非守恒形式：守恒形式：無粘流歐拉(Euler)方程非定常三維可壓縮無粘流動的控制方程總結(jié)如下：2.動量方程非守恒形式：X方向：Y方向：Z方向：無粘流歐拉(Euler)方程非定常三維可壓縮無粘流動的控制方程總結(jié)如下：2.動量方程守恒形式：X方向：Y方向：Z方向：無粘流歐拉(Euler)方程非定常三維可壓縮無粘流動的控制方程總結(jié)如下：3.能量方程非守恒形式：無粘流歐拉(Euler)方程守恒形式：關(guān)于控制方程的注釋關(guān)于控制方程的注釋連續(xù)性方程、動量方程、能量方程共有5個，但有六個未知的流場變量：關(guān)于控制方程的注釋在空氣動力學中，通常假設(shè)氣體是完全氣體（分子間作用力可忽略），狀態(tài)方程是：狀態(tài)方程提供了第6個方程，但引進了第七個未知量：溫度T關(guān)于控制方程的注釋用以封閉整個方程組的第七個方程必須是狀態(tài)參量之間的熱力學關(guān)系。比如：對常比熱容完全氣體，這個關(guān)系可以是：其中的是定容比熱。這個方程有時候也被稱為量熱狀態(tài)方程。物理邊界條件物理邊界條件無論流動是波音747飛機周圍的流動、亞聲速風洞內(nèi)的流動，還是流過一個風車流動，控制方程都是相同的。然而，盡管流動的控制方程是相同的，可這些情形中流動卻是完全不同的。為什么會這樣的呢？差異是哪里產(chǎn)生的呢？物理邊界條件答案是邊界條件。不同的邊界條件，有時還包括初始條件，使得同一個控制方程得到不同的特解。物理邊界條件對于粘性流動，物面上的物理邊界條件有物面速度無滑移邊界條件和物面溫度邊界條件。物面速度無滑移邊界條件指：緊挨物面的氣流與物面之間的相對速度為零。即：在物面（對于粘性流動）物理邊界條件大部分粘性流動的物面溫度邊界條件要么給定一個常數(shù)作為壁面溫度，即在物面要么假設(shè)壁面為絕熱壁，即在物面物理邊界條件對于無粘流動，物面上唯一的物理邊界條件是法向速度為零邊界條件。也就是說物面上的流動與物面相切。在物面（對于無粘流動）物理邊界條件無論是粘性流還是無粘流，根據(jù)問題的不同，流場中不是物面的地方有多種不同類型的邊界條件。比如對于流過固定形狀管道的流動，應(yīng)該在管道的入口和出口有適合的入流和出流邊界條件。比如對于已知來流中的飛行物，則給定自由來流條件作為物體四周無窮遠處的邊界條件。適合CFD使用的控制方程適合CFD使用的控制方程守恒變量：非守恒變量：適合CFD使用的控制方程非守恒變量可以由守恒變量求出：適合CFD使用的控制方程守恒形式的控制方程：流動控制方程中的因變量是守恒變量。非守恒形式的控制方程：流動控制方程中的因變量是非守恒變量。適合CFD使用的控制方程守恒形式的控制方程相比非守恒形式控制方程的第一個優(yōu)點：守恒形式的控制方程為算法設(shè)計和編程計算提供了方便。守恒形式的連續(xù)性方程、動量方程和能量方程可以用同一個通用方程來表達，這有助于計算程序的簡化和程序結(jié)構(gòu)的組織。適合CFD使用的控制方程守恒形式的控制方程組都可以表達成如下形式：U,F,G,H,J都是列向量。適合CFD使用的控制方程守恒形式的控制方程組都可以表達成如下形式：對于無粘或粘性流動：適合CFD使用的控制方程守恒形式的控制方程組都可以表達成如下形式：對于無粘流動：適合CFD使用的控制方程守恒形式的控制方程組都可以表達成如下形式：對于粘性流動：適合CFD使用的控制方程守恒形式的控制方程組都可以表達成如下形式：對于粘性流動：