修改自https://www.xsim.info/articles/OpenFOAM/Switch-multiple-versions.html，致谢大牛。

要添加其他版本只需仿照5-7行进行修改即可。

将该脚本拷贝到/usr/local/bin并赋予执行权限：

sudo cp ./switchFoamVersion /usr/local/bin/ && sudo chmod +x /usr/local/bin/switchFoamVersion

切换OF版本时一定要在命令前加"."或"source"，不然环境变量不能正确地传递到当前shell。

3967c521-2e7a-46bf-81e8-d0bb3be3fe33-屏幕截图 2023-09-12 115521.png

switchFoamVersion

自带的ensightWriter就能整体输出结果

@nuaawqf您好，请问您知道如何安装sowfa6吗

最近在看DPMFoam的代码，但是发现动量方程好像和其数学表达式对不上

fvVectorMatrix UcEqn ( fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc) - fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc) + continuousPhaseTurbulence->divDevTau(Uc) == (1.0/rhoc)*cloudSU ); UcEqn.relax(); fvOptions.constrain(UcEqn); volScalarField rAUc(1.0/UcEqn.A()); surfaceScalarField rAUcf("Dp", fvc::interpolate(rAUc)); surfaceScalarField phicForces ( fvc::flux(rAUc*cloudVolSUSu/rhoc) + rAUcf*(g & mesh.Sf()) ); if (pimple.momentumPredictor()) { solve ( UcEqn == fvc::reconstruct ( phicForces/rAUcf - fvc::snGrad(p)*mesh.magSf() ) ); fvOptions.correct(Uc); }

此处通过fvc::reconstrct返回的

fvc::snGrad(p)*mesh.magSf()

是否需要乘上alpha

http://dyfluid.com/class.html 你看我这个里面有个动图是不是你想要的？

我对造波消波不是特别熟悉，我只是拿着现成的方程植入一下。你最后说的是有道理的。一方面这个不是特别难（我都可以整。。可能比较深的比较难吧），另一方面目前openfoam自带了这个功能。所以wave2foam就没有开发了。

不过也有可能开发人员忙别的事情去了。他们这波人我不熟悉。谁知道呢 。哇咔咔

@队长别开枪
感谢前辈，采用几何重构方法后确实解决了这个问题~:ok3:

@李东岳 intensive variable是密集变量的意思么？图中那个截面的网格进行了网格加密，加密尺寸1mm。处理的变量是气相温度、组分质量分数这些。

@LiuWenchao 好的，非常感谢！

库朗数减小试一下 也就是时间步长减小试一下

看起来可以用。但是你需要先装openfoam，然后告诉visualCFD的openfoam的安装路径，然后他才能调用

谢谢你的回复。我自己试了几次基本上弄明白。没有物性的源文件不能进行计算，软件恢自动报错

为什么我用的

forAll(prticles,i) { particles[i] }

这个particles[i]取不到particle

感谢，有一个问题，为什么不把入口直接设置为湍流入口，是因为计算的雷诺数不高嘛

更新科氏力源项：$\mathbf{f}\times\bfU=2\Omega\sin \lambda \times\bfU$，其中$\mathbf{f}$表示科氏力参数，一些文献会直接给出$\mathbf{f}$的值。如果绕z旋转，那么$\mathbf{f}$在$z$方向存在值。在这种情况下，$\mathbf{f}\times\bfU$的值为：
\begin{equation}
\mathbf{f}\times\bfU=\left[
\begin{matrix}
f_2u_3-f_3u_2\\
f_3u_1-f_1u_3\\
f_1u_2-f_2u_1\
\end{matrix}
\right]=\left[
\begin{matrix}
-f_3u_2\\
f_3u_1\\
0\
\end{matrix}
\right]
\end{equation}
如果进一步考虑速度方向是x方向的，那么变为：
\begin{equation}
\mathbf{f}\times\bfU=\left[
\begin{matrix}
0\\
f_3u_1\\
0\
\end{matrix}
\right]
\end{equation}
可以看出如果地球绕z旋转，对于x方向的速度，科氏力导致y方向的偏离。

$\Omega$表示地球旋转速率，单位为[弧度$/\mathrm{s}$]=$1/\mathrm{s}$，$\lambda$表示维度[无单位]。因此$f$的单位为[$1/\mathrm{s}$]

@李东岳
谢谢李老师的回复。。。
我的 internal Fields 没 map 上，可能有地方没设置对

谢谢老师

@coolhhh 感谢！

@xjwang

CFDRFG的附录code还有个问题，求解$K$，不是按照论文中的原理计算的。他先用$PQ$没有代入$S(f)$的值，算了一组$K$，然后所有坐标点都用了同一组$K$。而根据理论，应该是要根据每个坐标点的非均匀$PQ$，求对应该点的$K$。DSRFG如果生成非均匀，理论上也要这么计算。
739858e5-ac4f-4736-94bd-9f4ec9f07a9d-image.png

关于CDRFG方法，他们的code里是对每一个点都生成一组PQK

最近看到一篇JFM文章《Guo, H., et al., An efficient and low-divergence method for generating inhomogeneous and anisotropic turbulence with arbitrary spectra. Journal of Fluid Mechanics, 2023. 970: p. A2.》，用的就是三维能谱，不是频谱，生成非均匀湍流场也是low-divergence。

@luofq-sysu 这个靠谱 差不多 嗯

通过vorticity求解涡量； 写个程序求解$\nabla w$
2.1 在这个程序里面，$\nabla w$的边界场$\cdot\bfn$再乘以粘性系数

说简单也简单，只不过繁琐了一些，好像没有几个命令就能出来的东西

可以尝试用functionObject中的surfaces提取物体表面的压力