@李东岳 谢谢李老师的关注。解决方法是paraview执行D3,再cellData2PointData,不需要写ghostLevel 信息就行。
line
帖子
-
.pvtu文件ghostLevel=1 -
.pvtu文件ghostLevel=1cell-center网格,pvtu文件ghostLevel=1 文件要怎么写 paraview才能识别内边界,找了好久没找到这方法资料,包含vtk格式说明,里面也没讲清楚。也试了很多方法都没成功。或者换个其他文件格式也可以。
-
Tecplot处理叶轮内部相对速度@西湖冷月 在 Tecplot处理叶轮内部相对速度 中说:
找了一圈没发现后处理的专门板块,由于是fluent计算的结果,所以发到这儿吧。问题就是,fluent计算叶轮流场后,显示不出来相对速度,用tecplot打开也没有相对速度,但对于叶轮机械而言,相对速度又很重要,有前辈解决了该问题的,可否指教一下如何去实现tecplot中得到相对速度。初步的想法是,因为叶轮转速一定,可以在tecplot中新建一个旋转速度,得到叶轮内圆周速度分布,但如何把圆周速度和绝对速度矢量相加得到相对速度呢?
可以通过脚本文件处理计算成相对速度。速度分解因为你有x,y坐标 转速就可以通过半径转化到直角坐标系下某点得速度。
-
pointwise中如何设置周期性边界条件本人pointwise小白,一直不知道pointwise如何设置周期性边界条件,望解答,谢谢!
-
超声速湍流平板计算问题网格有没有再加密,对流通量精度是多少,一阶精度会比二阶低很多
-
变比热容计算之前有点问题,最终如下
\par $C_p$采用4次多项式分段拟合
\begin{equation}\label{equ:NS_Cp}
C_p(T)=a_1+a_2T+a_3T^2+a_4T^3+a_5T^4
\end{equation}
\par静焓
\begin{equation}\label{equ:NS_H}
H(T)=\int_{T_{0}}^{T_{x}}C_p(T)dT=a_1T+\frac{a_2}{2}T^2+\frac{a_3}{3}T^3+\frac{a_4}{4}T^4+\frac{a_5}{5}T^5+a_6
\end{equation}
\par熵
\begin{equation}\label{equ:NS_S}
S(T)=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}=a_1\ln{T}+a_2T+\frac{a_3}{2}T^2+\frac{a_4}{3}T^3+\frac{a_5}{4}T^4+a_7
\end{equation}
\par$\bullet$求解静温(已知总温和马赫数)
\begin{equation}\label{equNSUs}
U_s^2=2\left(H(T_{tot})-H(T_{sta})\right)
\end{equation}
\begin{equation}\label{equNSmach}
Mach^2=\frac{U^2}{\gamma(T)R_gT}=\frac{U^2}{\frac{C_p(T)}{C_p(T)-R_g}R_gT}
\end{equation}
\par由(\ref{equ:NSmach})和(\ref{equ:NSUs})得
\begin{equation}\label{equNSTsta}
T_{sta}=\frac{2\left(H(T_{tot})-H(T_{sta})\right)}{Mach^2\frac{C_p(T_{sta})R_g}{C_p(T_{sta})-R_g}}
\end{equation}
$\bullet$求解静压(已知总温、总压和静温)
\par由p等熵过程
\begin{equation}
ds = C_p(T)\frac{dT}{T} -R_g\frac{d p}{p}=0
\end{equation}
\par两边同时积分有
\begin{equation}
\int_{T_{tot}}^{T_{sta}}C_p(T)\frac{dT}{T} =\int_{p_{tot}}^{p_{sta}} R_g\frac{d p}{p}
\end{equation}
\par记
\begin{equation}
S(T_{x})=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}
\end{equation}
\par则
\begin{equation}
S(T_{sta}) - S(T_{tot}) = R_g\ln\frac{p_{sta}}{p_{tot}}
\end{equation}
\par那么
\begin{equation}\label{equNS_psta}
p_{sta}=p_{tot}e^{\left(\frac{S(T_{sta})-S(T_{tot})}{R_g}\right)}
\end{equation} -
变比热容计算李老师@李东岳
根据我的推导,对于可压缩理想气体 应该满足如下的关系,但是从公式看出 这是要迭代求解的 ,推导是否正确
\par$\bullet$求解静温(已知总温和马赫数)
\begin{equation}\label{equ:NSUs}
U_s^2=2\left(h_t(T_{tot})-h_t(T_{sta})\right)=2\left(C_p(T_{tot})T_{tot}-C_p(T_{sta})T_{sta}\right)
\end{equation}
\begin{equation}\label{equ:NSmach}
Mach^2=\frac{U^2}{\gamma(T)R_gT}=\frac{U^2}{\frac{C_p(T)}{C_p(T)-R_g}R_gT}
\end{equation}
\par由(\ref{equ:NSmach})和(\ref{equ:NSUs})得
\begin{equation}\label{equ:NSTsta}
T_{sta}=\frac{C_p(T_{tot})T_{tot}}{C_p(T_{sta})\left(1+\frac{1}{2}Mach^2\frac{R_g}{C_p(T_{sta})-R_g}\right)}
\end{equation}
\par$\bullet$求解静压(已知总温、总压和静温)
\par由(\ref{equ:NSdsds})得,等熵过程
\begin{equation}
ds = C_p(T)\frac{dT}{T} -R_g\frac{d p}{p}=0
\end{equation}
\par两边同时积分有
\begin{equation}
\int_{T_{tot}}^{T_{sta}}C_p(T)\frac{dT}{T} =\int_{p_{tot}}^{p_{sta}} R_g\frac{d p}{p}
\end{equation}
\par记
\begin{equation}
S(T_{x})=\int_{T_{x}}^{T_{0}}C_p(T)\frac{dT}{T}
\end{equation}
\par则
\begin{equation}
S(T_{sta}) - S(T_{tot}) = R_g\ln\frac{p_{sta}}{p_{tot}}
\end{equation}
\par那么
\begin{equation}\label{equ:NS_psta}
p_{sta}=p_{tot}e^{\left(\frac{S(T_{sta})-S(T_{tot})}{R_g}\right)}
\end{equation} -
变比热容计算理想气体比热容是温度的函数,
那么有个问题:
已知总温总压 马赫数 怎么求得静温静压?
望解答