openfoam學(xué)習(xí)筆記2
OpenFOAM中的postProcess可以調(diào)用大量的函數(shù)來進(jìn)行數(shù)據(jù)處理
postProcess -list
CourantNo(輸出庫(kù)朗數(shù))
Lambda2(輸出Lambda2)
MachNo(輸出馬赫數(shù))
PecletNo(輸出Peclet數(shù))
Q(計(jì)算速度梯度第二不變量)
Qdot(計(jì)算并輸出當(dāng)前燃燒模型的放熱率)
R(計(jì)算雷諾應(yīng)力張量場(chǎng))
XiReactionRate(寫入湍流火焰速度以及反應(yīng)速率體標(biāo)量場(chǎng))
add(添加場(chǎng))
age(計(jì)算并寫出粒子從入口到位置所需的時(shí)間)
boundaryProbes(輸出場(chǎng)在某patch 上指定位置的值)
cellMax(對(duì)于一個(gè)或多個(gè)場(chǎng),輸出最大的值)
cellMin(對(duì)于一個(gè)或多個(gè)場(chǎng),輸出最小的值)
components(寫入矢量場(chǎng)的分量)eg.postProcess -fun “component(U)”
ddt(計(jì)算場(chǎng)的時(shí)間導(dǎo)數(shù))
div(計(jì)算場(chǎng)的散度)
dsmcFields(從DSMC計(jì)算結(jié)果中提取UMean、translationalT、internalT以及overallT)
enstrophy(計(jì)算場(chǎng)的渦度擬能)
faceMax(對(duì)于一個(gè)或多個(gè)場(chǎng),輸出面上最大的值)
faceMin(對(duì)于一個(gè)或多個(gè)場(chǎng),輸出面上最小的值)
faceZoneAverage(計(jì)算 faceZone 上一個(gè)或多個(gè)場(chǎng)的平均值)
fieldAverage(計(jì)算并寫入給定場(chǎng)的時(shí)間平均值)
flowRateFaceZone(計(jì)算patch 特定面區(qū)域的流量。有可能是體積通量,也有可能是質(zhì)量通量)
flowRatePatch(計(jì)算patch 的流量。有可能是體積通量,也有可能是質(zhì)量通量)
flowType(輸出速度場(chǎng)類型,-1表示有旋流;0表示剪切流;+1表示其他類型流動(dòng))
forceCoeffsCompressible(對(duì)于可壓縮求解器,通過對(duì)某patch 上的力進(jìn)行加和匯總計(jì)算升力、壓力、矩系數(shù))
forceCoeffsIncompressible(對(duì)于不可壓縮求解器,通過對(duì)某patch 上的力進(jìn)行加和匯總計(jì)算升力、壓力、矩系數(shù))
forcesCompressible(對(duì)于可壓縮求解器,計(jì)算某patch 上的壓力以及粘性力)
forcesIncompressible(對(duì)于不可壓縮求解器,計(jì)算某patch 上的壓力以及粘性力)
grad(計(jì)算場(chǎng)的梯度)
interfaceHeight(對(duì)于一系列的點(diǎn),輸出界面的高度。這個(gè)高度是在這個(gè)點(diǎn)上,界面和壁面的垂直距離)
internalProbes(寫出內(nèi)插到指定點(diǎn)云的字段值)
log(計(jì)算指定標(biāo)量場(chǎng)的自然對(duì)數(shù))
mag(對(duì)矢量進(jìn)行操作提取矢量分量以及模的大?。?/span>
magSqr(計(jì)算場(chǎng)的模的平方)
minMaxComponents(通過非標(biāo)量的形式,輸出最大最小值以及位置)
minMaxMagnitude(通過非標(biāo)量的形式,輸出最大最小值的模以及位置)
particles(跟蹤由連續(xù)相流驅(qū)動(dòng)的包云)
patchAverage(計(jì)算patch中一個(gè)或多個(gè)場(chǎng)的平均值)
patchIntegrate(計(jì)算一個(gè)patch中一個(gè)或多個(gè)場(chǎng)的曲面積分)
phaseMap(寫入相位分?jǐn)?shù)圖域alpha.map,每個(gè)相位的數(shù)值范圍遞增,例如,水的數(shù)值為0,空氣為1,油為2)
phaseScalarTransport(在多相模擬的一個(gè)階段內(nèi)求解標(biāo)量場(chǎng)的輸運(yùn)方程)
pressureDifferencePatch(計(jì)算兩個(gè)patch 的平均壓力差)
pressureDifferenceSurface(將壓力插值在三角面并計(jì)算平均壓力差)
probes(輸出場(chǎng)在指定位置距離最近的網(wǎng)格單元的值)
randomise(對(duì)場(chǎng)添加隨機(jī)分量,擾動(dòng)量可以指定)
residuals(對(duì)于指定場(chǎng),輸出每個(gè)時(shí)間步迭代的初始?xì)埐?。?duì)于矢量場(chǎng),輸出最大的分量)
scalarTransport(求解標(biāo)量場(chǎng)傳輸方程)
scale(對(duì)場(chǎng)進(jìn)行縮放)
shearStress(計(jì)算剪應(yīng)力,輸出數(shù)據(jù)為 volSymmTensorField)
singleGraph(提取某一條線上的數(shù)據(jù))
staticPressure(通過指定的密度計(jì)算動(dòng)壓)
stopAtClockTime(當(dāng)達(dá)到指定的時(shí)鐘時(shí)間(以秒為單位)時(shí),停止運(yùn)行,并在停止前可選擇寫入結(jié)果)
stopAtFile(當(dāng)文件stop在case目錄中被創(chuàng)建時(shí),停止運(yùn)行)
streamFunction(計(jì)算網(wǎng)格點(diǎn)的流函數(shù);計(jì)算面通量)
streamlines(通過VTK 格式輸出流線數(shù)據(jù))
subtract(從某個(gè)場(chǎng)中減去某個(gè)場(chǎng)的值)
surfaces(通過VTK 格式輸出切割面、等值面等)
time(寫入運(yùn)行時(shí)間、CPU 時(shí)間、鐘表時(shí)間)
totalEnthalpy(計(jì)算并寫入總焓作為 volScalarField)
totalPressureCompressible(計(jì)算可壓縮求解器的總壓)
totalPressureIncompressible(計(jì)算不可壓縮求解器的動(dòng)力總壓)
turbulenceFields(計(jì)算給定的湍流場(chǎng))
turbulenceIntensity(計(jì)算湍流強(qiáng)度場(chǎng))
volFlowRateSurface(計(jì)算三角面的體積通量,速度的計(jì)算采用對(duì)面表面進(jìn)行插值并積分。建議三角的面積足夠?。?/span>
vorticity(計(jì)算渦量場(chǎng),例如速度的旋度)
wallHeatFlux(計(jì)算壁面熱通量,輸出體矢量場(chǎng))
wallHeatTransferCoeff(計(jì)算邊界場(chǎng)的壁面換熱系數(shù))
wallShearStress(計(jì)算壁面剪切力,輸出體矢量場(chǎng))
writeCellCentres(通過體矢量場(chǎng)的方式輸出網(wǎng)格體心矢量)
writeCellVolumes(通過體標(biāo)量場(chǎng)的方式輸出網(wǎng)格單元體積)
writeObjects(輸出某些特定場(chǎng))
writeVTK(以 VTK 格式寫出指定的對(duì)象)
網(wǎng)格非均勻分布
自己體會(huì)
時(shí)間步的調(diào)整
當(dāng)使用一個(gè)非均勻網(wǎng)格的時(shí)候,blockMesh會(huì)計(jì)算網(wǎng)格大小。如果設(shè)置長(zhǎng)度為l,內(nèi)置n個(gè)網(wǎng)格單元,最末端和起始端網(wǎng)格比為R,那么最小的網(wǎng)格大小為:
r是相鄰網(wǎng)格單元的大小比:
且有:
映射場(chǎng)數(shù)據(jù)
算例文件夾下
mapFields <路徑> -consistent
復(fù)制算例
openfoam可以采用foamCloneCase命令進(jìn)行算例復(fù)制,默認(rèn)情況拷貝0文件夾,通過指定-latestTime可以拷貝最新的文件夾作為新算例的初始場(chǎng)
foamCloneCase -latestTime <舊文件夾名> <新文件夾名>
eg.foamCloneCase -latestTime cavity cavityHighRe
運(yùn)行算例
后臺(tái)運(yùn)行
nohup nice -n 19 icoFoam > log.icoFoam &
nohup當(dāng)用戶戶退出登陸的時(shí)候,這個(gè)程序依然會(huì)繼續(xù)運(yùn)行 nice 調(diào)整進(jìn)程優(yōu)先級(jí),-20 對(duì)應(yīng)最高優(yōu)先進(jìn)程,19 對(duì)應(yīng)最低優(yōu)先進(jìn)程
高雷諾數(shù)流動(dòng)
$FOAM_RUN/tutorials/incompressible/pisoFoam/RAS 的cavity 算例
Reynolds-Average Simulation(RAS)湍流模型
附帶壁面函數(shù)的標(biāo)準(zhǔn)模型來求解,0文件夾下出現(xiàn)兩個(gè)新變量
:湍流動(dòng)能場(chǎng),
?:湍流動(dòng)能耗散率場(chǎng)
當(dāng)使用附帶壁面函數(shù)的標(biāo)準(zhǔn)模型的時(shí),沒有必要引入非均勻網(wǎng)格。原因?yàn)楸诿娓浇牧餍鸵呀?jīng)被?;皇侵苯忧蠼狻?/p>
0/nut文件
0/k文件
0/epsilon文件
?
OpenFOAM還有RAS湍流模型以及LES大渦模擬