CFD中文网

@bestucan hhhh谢谢您，我发现我哪里错了，我画网格的时候是xoz面上的，正确的应该是在xoy画个2D的网格，然后openfoam会在z轴增加一个面。

phi是用来调用判断速度方向的： surfaceScalarField phiBD(upwind<scalar>(psi.mesh(), phi).flux(psi));

在做液滴撞击的模拟，因为要用到动态接触角所以需要获取三相点位置以及速度。目前是想用两步迭代的位移差来计算速度，但是现在不知道怎么去算这个三相点的位置。代码是这么写的

Thread* thread = Lookup_Thread(domain, 9); Thread** pt = THREAD_SUB_THREADS(thread); cell_t cell; face_t f; real x[ND_ND], area[ND_ND]; double max_x = 0, f_Hoff_inverse, x_hoff, temp, Ca, volume; int n; sum = 0.0; begin_c_loop_all(cell, pt[1]) { if (C_VOF(cell, pt[1]) != 0) { C_CENTROID(x, cell, pt[1]); if (x[0] > max_x) max_x = x[0]; } sum += C_VOF(cell, pt[1]); } end_c_loop_all(cell, pt[1])

我想问一下怎么把SUM变成面积啊

@zp气固两相流 只有考虑面的时候才有一个own和一个nei，考虑cell的时候一个own，其他都是nei，所以考虑cell的时候，一个own，很多nei很正常啊

@李东岳 :high:

我是按照《普朗特流体力学基础》第208页进行翻译的，感兴趣你也可以看看这一节，不过这里没有将太细 :quwan:

想请教一下各位大佬关于Fluent，可不可以用来模拟电场下的流体形变过程？或者硬质的基板，可以定义为软、可以变形的薄膜吗？

无锡超算，用openfoam计算 :

国产申威单个核组（仅1主核，无从核加速），与intel商用节点的单核计算效率，有谁比较过吗？

http://www.easycae.cn/
easycae云计算平台已经有了基本功能，已经在上面算过一些fluent算例，很不错，上面提供超算资源

thermophysicalProperties中energy可以选择sensibleEnthalpy, sensibleInternalEnergy 和 absoluteEnthalpy。前两者与后者的区别，后者使用绝对能量表示包含热源的情况。而sensibleEnthalpy与sensibleInternalEnergy的区别，前者是利用焓值求解温度，后者是利用内能求解温度。热⼒学中分析闭式系统⽤内能，分析开式系统⽤焓值。那么在实际使用中，该如何选择这两种模式？例如高压离心风机。tutorials中的案例似乎找不出规律。
eb324645-1dd5-499a-bfc0-6d77667f7fa5-image.png

还有就是如果网格节点无法保证周期性一一对应，可以将对应的周期面设置为interface,然后到fluent里通过匹配两个interface面，并选择对话框下方的Periodic boundary conditions设置
另外一款网格划分软件Gambit的话就可以直接设置多对周期性边界条件

@李东岳 噢噢好，谢谢大佬指点！:sad:

之前有点问题，最终如下
\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}

我是做四方柱的涡激运动，用的是DBFI和vof，有自由液面，采用的是静水波，惯性矩和质心等参数是利用三维建模软件犀牛计算的，横荡和纵荡的自由衰减都算对了，但是艏揺的自由衰减就是不对，有大佬能帮忙检查一下看一下是哪里有问题吗，可有偿![替代文字]微信图片_20220714094457.png微信图片_20220714094605.png

在代码里加入

volScalarField rho ( IOobject ( "rho", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh );

然后直接用就行

请问您是怎么解决的

@李东岳 大神，这个我真没注意。好像不多。我是土木专业的，不是搞固体力学的。您可以再找找。

是否为容差 matchTolerance 太大导致的衔接不上~

一般什么原因比较多

@东岳 老师人太好了！再次感谢！