CFD中文网

在1954年，IBM创立了Fortran语言。

在1955年，多相流中的粒子法就开始进行了。在当时，其被取名为Particle In Cell (PIC)。

PIC由Los Alamos国家实验室的F.H. Harlow提出。在PIC中，拉格朗日粒子通过其坐标（x, y, z）进行定位并且具有质量。有些情况下，拉格朗日粒子还具有多个物理特征。例如在PIC中，拉格朗日粒子可能还具有大小的区别等。

@李东岳 东子老师，重装就可以了:140:

各位老哥好，有个小问题，我目前现在计算域中取一个矩形平面区域，这个矩形平面用于注入粒子，有没有相关的InjectionModel呢？

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

具体是啥阀呢，阀内流体是水还是油，我可以试试，工程问题可以用cfd计算一下阀内的噪声，但都是参考，可以通过cfd进行可视化研究，可以看到里面的流场和声场，具体咨询合作可以加微信18647841059

亲测，可以对齐~如图，word2016，自带的公式编辑器！
0_1520851049441_d3fa154b-70c5-438f-9674-805ee760bf75-image.png

@李东岳 请问在哪里下载？东岳流体？

0_1519692297906_捕获.JPG
是否是这样？应力张量我可以理解。我的疑问是速度梯度在笛卡尔坐标系下变换之后直接就变成在极坐标系中了速度梯度了？

@东岳 老师，我重新画网格后，decomposePar经常会报错，但是不影响后面计算，没有查到原因0_1537492983872_0403f9aa-ee80-4134-a228-ac6289ba58de-image.png

@东岳 内网格，尝试着接着划吧

@李东岳 http://dx.doi.org/10.1016/j.compchemeng.2015.09.008
这篇文章对三种重构方法做了对比。

@郑学习 好的，感谢回复

尖的地方可以试试将多余的节点进行合并成一个节点映射到尖尖的点上，然后对锥块做一个Y剖，但是这样的话在尖尖位置的网格要么太密要么质量差

测试出来了QAQ：
𝐵𝑖𝑛 𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛 𝑠𝑜𝑢𝑟𝑐𝑒=(𝐵𝑖𝑛 𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛 𝑐ℎ𝑎𝑛𝑔𝑒×𝐶𝑢𝑟𝑟𝑒𝑛𝑡 𝑝ℎ𝑎𝑠𝑒 𝑣𝑜𝑙𝑢𝑚𝑒 𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛)/(𝐶𝑒𝑙𝑙 𝑣𝑜𝑙𝑢𝑚𝑒 ×𝑇𝑖𝑚𝑒 𝑠𝑡𝑒𝑝 𝑠𝑖𝑧𝑒)

Eg: Bin 1 bin fraction decrease 2% in 0.1s, cell volume is 0.001 m3, inclusion phase volume fraction is 1%

322c153e-e638-4775-9e9a-4e0939981818-image.png
2 is the number that we give fluent, 0.02 is actual bin fraction change.

可能是版本不同的问题，在OF2206中，添加缓存没有成功。:mianmo:

用另一种方法，算是曲线完成任务了。

在H文件中，添加 volScalarField NN;，在C文件初始化的列表中，加入

NN ( IOobject ( "NN", this->runTime_.timeName(), this->mesh_, IOobject::NO_READ, IOobject::AUTO_WRITE ), this->mesh_ ),

然后在恰当的位置，完成赋值 NN = N2 即可。

不过有一个缺陷，需要再初始文件夹中设定 NN 的初始文件。

(先解决问题，麻烦就麻烦点吧)

链接文本

参考页面左边栏。

openfoam-v1912，求解器overInterDyMFoam，官方算例floatingObject

最近在进行重叠网格计算时，浮体采用snappyHexMesh生成复杂结构物（浮体），并拉伸为二维网格。

背景网格同样使用snappyHexMesh生成固定结构物，并拉伸为二维网格。

在进行计算时，第一步压力就发散，更改压力求解方式，网格尺寸，重叠区域插值方式等均未解决。

在背景网格使用toposet可以生成简单结构物，设置同上，可以计算。

官网算例aeroFiol中背景网格仅使用snappyHexMesh进行网格加密，未见到背景网格使用snappyHexMesh生成结构物。

请问可以大神们有没有解决思路？背景网格是否支持使用snappyHexMesh生成结构物？

如果想在背景网格生成固定的复杂结构物，该怎么操作呢？

@队长别开枪 @李东岳 @wwzhao

@郑学习 你这个是啥版本，我编译的路径不大一样64460b45-3c82-4d68-b1a3-f907b95b23c5-image.png 这个路径下面有很多拉格朗日的底层代码，我这边是在这个下面进行代码植入，然后返回上一个层级进行编译的。你这个添加的地方似乎也行，0454ef55-a345-4790-b1bf-f40fb6724612-image.png 我认为应该是这个路径下面进行编译，祝成功！

你注意一下openfoam版本，不同的版本functions写法不一样