資源描述：
《comsol地球科學-溶質運移-theory》由會員上傳分享，免費在線閱讀，更多相關內容在教育資源-天天文庫。
1、COMSOLMultiphysics“地球科學模塊” 中級培訓——溶質運移溶質運移“SaturatedPorousMedia”模式“VariableSaturatedPorousMedia”模式“SaturatedPorousMedia”模式偏微分方程“VariableSaturatedPorousMedia”模式偏微分方程應用模式屬性（ApplicationModeProperties）求解域控制方程設定（SubdomainSettings）邊界條件設定（BoundarySettings）點、邊設定（PointandEdgeSettings）內容對流吸附和阻滯水里擴散反應保守和非保守公式c
2、=soluteconcentrationcP=adsorbedconcentrationqs=fluidvolumefractionu=fluidvelocityrb=bulkdensityDL=dispersiontensorRL=aqueousreactionsRP=solidreactionsSc=solutesourceFluidandSolidrxns水動力彌散對流吸附“SaturatedPorousMedia”模式Advection—水平對流溶質隨著運動著的流體移動的過程稱為對流。追蹤溶質粒子的位置是通過平均線性流速—thevectorofaveragelinearfluidve
3、locitiesθs代表流體體積含水率，對于自由流體其值為1，對于多孔介質，為孔隙率，值在0.1-0.5之間。對于溶質輸送應用模式，用戶提供的是通過多孔介質單位面積的流體速度，比如特定的或者達西定律計算出的流網的流速。而不是孔隙內流體的流速uaSorptionandRetardation—吸附、解吸以及阻滯cP—adsorbedconcentration應用模式提供三種預定義的計算方式滯留因子RF—表征污染物在遷移過程中相對于水流被滯后的程度HydrodynamicDispersion—水動力彌散機械彌散與分子擴散共同作用的結果，形成了水動力彌散。反應—Reactions環(huán)境中的化學反應會影
4、響溶質運移，比如生物降解、放射性衰退、變質、發(fā)熱、吸熱反應等。用戶可以指定偏微分方程等式右端的反應項。通常的反應包括：保守和非保守形式在地球科學溶質運移應用模式有保守和非保守形式的控制方程可以選用保守形式的控制方程——適用可壓縮流體，流體的密度在空間是變化的非保守形式的控制方程——適用不可壓縮流體，不可壓縮流體對流項速度的散度是0“VariableSaturatedPorousMedia”模式c=aqueousconcentrationcP=adsorbedconcentrationcG=gasconcentrationqs=fluidvolumefractionrb=bulkdensity
5、av=gasvolumeu=fluidvelocityFluid,Solid,andGasrxnsDL=liquiddispersiontensorDG=gasdispersiontensorRL=aqueousreactionsRP=solidreactionsRG=gasreactionsSc=solutesource對流、吸附和阻滯同飽和溶質運移應用模式水動力彌散在前述基礎上補充了流體-氣體彌散張量2Dsystem反應增加了氣相保守非保守形式應用模式屬性（ApplicationModeProperties）求解域控制方程設定邊界條件設定（BoundarySettings）點、邊設定（P
6、ointandEdgeSettings）Example：VariablySaturatedFlowandTransport算例描述：如右圖所示，帶有化學品的圓環(huán)內的水向下方的土柱中移動，當流動穿過非飽和的土柱的時候，化學物黏附在土體顆粒上，減緩了溶質相對于水的輸送。另外，化學品濃度在水體和土體中發(fā)生生物降解而衰退。通過模擬追蹤溶質的運移過程。是一個2D軸對稱模型?？刂品匠蘁luidFlowvanGenuchtenformulas地下水流動和溶質運移通過流體流速耦合在一起SoluteTransport為了求解更方便，展開上式左項建模步驟：Physics>ConstantsPhysics>the
7、ScalarVariables定義幾何模型Draw>SpecifyObjects>RectanglePARAMETEREXPRESSIONwidth1.25height1.3r0z-1.3Draw>SpecifyObjects>RectanglePARAMETEREXPRESSIONwidth1.25height0.4r0z-0.4Draw>SpecifyObjects>Liner->00.25,z->00.
