CFD中文网

最近要做一些跨临界的问题，但是在跨临界中，如下图，由于伪沸腾的影响，流体的热力学性质会在伪沸腾区呈现剧烈的变化，形成类似于激波的区域。但是在亚临界区域，压力对密度的变化极为敏感，所以传统的全守恒格式会引起压力震荡。很多研究者针对这个问题做了很多修正。比如：求解压力传播方程并加入人工扩散、计算左右通量求黎曼解并对最后的结果熵修正。但是目前的算法总体来说很繁琐，且由于通过压力反推能量会引起能量不守恒。主要有两个问题，非常感谢您的解答：

求解压力传播方程的那个算法很类似simple算法的pressure based方法，为什么现在几乎所有的计算空气动力学论文里面几乎没有使用过simple算法，更倾向于TVD RK这种显示迭代？ 显式格式很适合GPU加速计算，想开发一套GPU的代码用来计算。但是GPU对双精度计算的支持很差。当在激波附近使用高阶离散格式的时候如果使用单精度浮点数计算，截断误差是否会对计算结果产生比较大的影响？
3.png

对以空气为介质，中间有孔板的管路进行模拟，使用周期性边界条件。
取了完整模型的1/18，也就是20°，用icem画的结构网格，边界条件和其他设置和完整模型都是一样的，但是云图结果和完整模型差的有点多，会有一坨高温区域。起初以为是网格不够，毕竟才21w网格，后面网格加到53w，结果没什么太大变化，想问一下，周期性边界是不是不适合我这种可压缩高速流动，或者我网格的画法有问题？
下面是几何模型和网格，底部弄了个y切，网格质量0.63。
d8c1d2d1-1867-484b-aeb9-db3965f704bd-image.png
6d636e0c-67f4-4578-ab38-28d3489b82f5-image.png
63ea7a28-ed09-48b1-993e-20487c0e1f7b-image.png 2ef7c281-fa1b-44e1-ac40-069826a8d2e1-image.png
53548316-414a-4941-80f9-d70a3df3e59b-image.png

@adcpk 超算每个节点20核，是因为只有20核还是只用了20核？如果是只用了20核，其他核在跑别的作业吗？如果只有20核，建议换成15+15这样来测试。

一般情况下超算平台的CPU性能肯定没有你工作站的好，工作站肯定比超算新吧，快一些的话其实也可以理解。

用的windAroundBuildings算例，840万网格的二维模型，网格质量检查过了，计算域是14000*4000米。epsilonInlet 0.0003; // Cmu^0.75 * k^1.5 / L ; L =10

dimensions [0 2 -3 0 0 0 0];

internalField uniform $epsilonInlet;

boundaryField
{
INLET
{
type fixedValue;
value uniform $epsilonInlet;
}

OUTLET { type zeroGradient; } WALL { type epsilonWallFunction; value uniform $epsilonInlet; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
kInlet 1.5; // approx k = 1.5*(I*U)^2 ; I = 0.1

dimensions [0 2 -2 0 0 0 0];

internalField uniform $kInlet;

boundaryField
{
INLET
{
type fixedValue;
value uniform $kInlet;
}
OUTLET
{
type zeroGradient;
}

WALL { type kqRWallFunction; value uniform $kInlet; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
nut
internalField uniform 0;

boundaryField
{
INLET
{
type calculated;
value uniform 0;
}

OUTLET { type calculated; value uniform 0; } WALL { type nutUSpaldingWallFunction; value uniform 0; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
p
internalField uniform 0;

boundaryField
{
INLET
{
type zeroGradient;
}
OUTLET
{
type fixedValue;
value uniform 0;
p0 uniform 0;
}
WALL
{
type zeroGradient;
}
frontAndBackPlanes
{
type empty;
}

#includeEtc "caseDicts/setConstraintTypes"

}
Uinlet (10 0 0);

dimensions [0 1 -1 0 0 0 0];

internalField uniform (10 0 0);

boundaryField
{
INLET
{
type fixedValue;
value uniform $Uinlet;
}
OUTLET
{
type zeroGradient;
}

WALL { type noSlip; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
时间步长试过1 0.1 0.01 其他的都没动

对于给定的NDF，划分为$i$个$N_{pp}$，对每个$i$上定义$k$阶矩$m_k^i$，给定$m_k^i$，可以计算第$i$区间的节点$d^i_0,d^i_1$以及权重$w^i_0,w^i_1$：
\begin{equation}
\begin{split}
w^i_0&=w^i_1=0.5
\\
d^i_0&=m_1^i-\frac{1}{\sqrt{3}}\sqrt{\frac{m_3^i}{m_1^i}-{m_1^i}^2}
\\
d^i_1&=m_1^i+\frac{1}{\sqrt{3}}\sqrt{\frac{m_3^i}{m_1^i}-{m_1^i}^2}
\end{split}
\end{equation}
对于仅考虑破碎的PBE：
\begin{equation}\label{pbe}
\frac{\p n(d)}{\p t}=\int_d^{d_{max}}g(d')\beta(d|d')n(d')\rd d'-g(d)n(d)
\end{equation}
对方程\eqref{pbe}在$i$上取$k$阶矩：
\begin{equation}\label{m}
\frac{\p m_k^i}{\p t}=\int_{d_{i-1/2}}^{d_{i+1/2}}\int_d^{d_{max}}g(d')d^k\beta(d|d')n(d')\rd d'\rd d-\sum^2_{j=0} g(d_j^i)w_j^i(d_j^i)^k
\end{equation}
\begin{equation}
\begin{split}
\int_{d_{i-1/2}}^{d_{i+1/2}}\int_d^{d_{max}}g(d')\beta(d|d')n(d')\rd d'\rd d&=
\int_{d_{i-1/2}}^{d_{max}}g(d')n(d')\left(\int_{d_{i-1/2}}^{d'}\beta(d|d')\rd d\right)\rd d'
\\&=
\sum_{m=i}^{N}\sum_{j=0}^2g(d_j^m)w_j^m\left(\int_{d_{i-1/2}}^{d_j^m}d^k\beta(d|d_j^m)\rd d\right)
\end{split}
\end{equation}
Therefore
\begin{equation}
\frac{\p m_k^i}{\p t}=\sum_{m=i}^{N}\sum_{j=0}^2g(d_j^m)w_j^m\left(\int_{d_{i-1/2}}^{d_j^m}d^k\beta(d|d_j^m)\rd d\right)-\sum^2_{j=0} g(d_j^i)w_j^i(d_j^i)^k
\end{equation}

@dxl 噢，那可以在相应的位置加一些测试语句试试。

386b6db0-5254-4f1c-917a-606c2fa41539-image.png
377452dc-2a13-489e-b437-19f5411e4443-image.png
1、图中是fluent2019R1版本，湍流模型选择为transition SST k-w模型，保持默认参数，Re_thetac,Re_thetat,F_length等参数取值均为none；
2、但是在fluent理论手册中转捩模型与SST模型耦合明确要用到Re_thetac这个参数，如果保持默认none，那么这个参数是有个默认值吗？还是有其他的处理方法。
恳请各位前辈指点，十分感谢。

topoSet有一些比如寻找点最近的面，寻找距离stl最近的面之类的，你试试

都柏林大学发来贺电！

thread_loop_c 获取subthread层的值似乎有点问题，见https://muchong.com/html/201211/5150612_2.html，最好是用mp_thread_loop_c，配合**pt和C_VOF(c, pt[0])来进行取值。

@李东岳 谢谢东哥，就是我把resolveFeatrueAngle调大之后画网格没有捕捉这条特征线，画几何的时候因为在两个不同平面上，这条特征线是存在的

这种边界条件不是特别稳定。很容易发散。所以总是报错。

@HITSC30 在 用Fluent做大涡模拟LES，用CDRFG或者NSRFG实现大气边界湍流入口 中说：

@shwjz 在 用Fluent做大涡模拟LES，用CDRFG或者NSRFG实现大气边界湍流入口 中说：

@HITSC30 在 用Fluent做大涡模拟LES，用CDRFG或者NSRFG实现大气边界湍流入口 中说：

@shwjz 在 用Fluent做大涡模拟LES，用CDRFG或者NSRFG实现大气边界湍流入口 中说：

@HITSC30 你好，udf的逻辑是不是这样：将生成的入流面每个节点的风速时程，读取当前时间步，提取对应时间步的风速赋予节点，相当于每个时间步都要更新一次入流面所有节点的u、v、w三向风速

这个程序我用的是之前师兄写过的程序，我自己没咋看过逻辑。最近在忙点其他事情，忙完这阵，我看一下代码，再和你交流。

谢谢，希望能交流一下，请教你是用哪种方法生成湍流入口？

你好，我用的是CDRFG。

好的，等您有时间，希望能与您邮箱交流一下（z15392569@163.com），我现在也在尝试CDRFG，但是遇到一些问题想请教

李博的桌面？

SA是Spalart-Allmaras，Reynolds-Averaged Navier-Stokes Spalart-Allmaras turbulence和Spalart-Allmaras Detached Eddy
Spalart-Allmaras采用的是标准的SA，不是高雷诺数的。DES也是标准的DES模型

@Junren-Hou 在 ICEM划分网格 中说：

可以试试check block，有些时候其实没有问题，就它软件自己的毛病，然后block里最后有一个check block，可以让他自己检查一下。里面有好多选项，都试试，什么check/fixed block

回忆了一下，我的旋转区域用ICEM画的应该没有问题。出问题的应该是我画的外流场区域，一方面追求地面边界层，致使网格尺度过小，一方面为了平衡网格数量，过渡不均匀。猜测应该是外流区域的地面边界层设置导致了相关报错的产生。

求解答

请问大家为什么我用meshing划分流体域网格时,明明两个面没有重叠，模型树上还会出现图片上这种接触区域呢？![6392L5}X9$J}C{GTK_ZZ7F.jpg

@李东岳 感谢李老师！现在对这个求解器理解了一些，对于$\mathbf{U}_1$的控制方程, $K_d \mathbf{U}_2$一项作为显式项被直接纳入pEqn.H的界面通量的修正。