CFD中文网

他这个聊天记录我没太懂。什么是拖布结构？

烦请各位大佬指点一下，workbench里面mesh画的带角度的翼型的网格，为什么边界层厚度不一样啊？怎么处理一下呢？上面的感觉太窄了

e9d97488-fa0f-453d-9f2b-1884d5136695-image.png

标准的限制方法是保证插值后网格顶点（vertice）的值不会出现极大值或极小值（比周围网格的最大值大或比周围网格的最小值小）因为梯度插值最容易出现越界值的就是在顶点处，从而用一个标量来scale梯度gradient；另外一种方法是保证网格的面心midponts of cell faces处不出现极大值或极小值即可，就认为满足了插值后不存在越界的问题。我想的是关于第二种方法，面心处插值后不出现最大值或最小值不能代表顶点处不出现最值，对吗

@Yan1 这个好像是采用致动线实现风机气动性能模拟的那个程序，以及后面实现风机气动水动的全耦合程序

@吴优 现在的问题是如何将我的我的结构求解模块和fluent联立起来。是如何获取wall的表面力用于插值到固体域网格。

我也想知道L怎么定义的，OpenFOAM2312，这么解释

“The order of elements of the integral scales for the vector-based condition is “(Lxu Lxv Lxw Lyu Lyv Lyw Lzu Lzv Lzw)” and for the scalar-based condition is “(Lxi Lyi Lzi)”, where \c x here indicates the streamwise components and \c y/z lateral components. The streamwise components (i.e. \c x) are Lagrangian time scales in units of seconds and the remaining components are length scales in units of meters.”

Lxi对应拉格朗日时间尺度，Lyi/Lzi对应长度尺度

@LYHua 太感谢了，解决了，给力

@李东岳 铁人都不睡觉的:136:

请问燃料燃烧时这种‘理论计算产生的平均总温’应该怎么计算啊？翻了基本燃烧学的书都没这方面的定义。温升燃烧效率.png

@李东岳 :ok3: :ok3: :ok3:

对于给定的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}

@Yan 您好，请问您是用的什么后处理软件，在导出速度的时候可以把笛卡尔坐标系改成柱坐标系

@evensun 嗯，涡粘系数是二阶张量这个按我理解其实就是非线性湍流模型了

https://www.bilibili.com/video/BV1or421x7Ut/?spm_id_from=333.999.0.0&vd_source=5dfca494f85875fa7f031a9f632b1a89
第9分钟介绍这种High Aspect Ratio对于瞬态、稳态的具体影响。

各位大佬好~向大家请教一个问题。

式（1）是大涡模拟法夫尔滤波操作以后的连续性方程
\begin{equation}
\frac{\partial{\bar{\rho}}}{t}+\frac{\partial{\bar{\rho}\tilde{u}_i}}{x_i}=0
\end{equation}

如果对$\tilde{u}_i$取时均值，即$<\tilde{u}_i>$, 这样得到的是雷诺时均值的速度吗？如果想要取法夫尔时均值，是不是需要做$<\bar{\rho} \tilde{u}_i>/<\bar{\rho}>$的处理？P.S(这里用<φ>表示对φ取时间平均)

第二个问题与第一个问题相似：基于法夫尔时均的可解尺度的湍动能的求解是否可以按照式（2）不可压缩流体的求解方式求解呢？我感觉式（2）的解更像是雷诺时均的解，因为可解尺度的湍流里面也包含了密度脉动。
\begin{equation}
k_{resolved}=\frac{1}{2}(<\tilde{u}_i\tilde{u}_i>-<\tilde{u}_i>^2)
\end{equation}

还是说需要用式（3）继续用法夫尔时均进行计算呢？
\begin{equation}
k_{resolved}=\frac{1}{2}\big(<\bar{\rho}\tilde{u}_i\tilde{u}_i>-<\bar{\rho}>(<\bar{\rho} \tilde{u}_i>/<\bar{\rho}>)^2\big)/<\bar{\rho}>
\end{equation}

@liujuncfd
本科待遇年薪15w左右，硕士年薪20w左右，博士年薪25w。
https://www.zhipin.com/job_detail/fe4f05de98c982591nN62Ni7GVZT.html?ka=hot-job-1

算不准大概率是网格问题，不仅仅是数量，还有质量，另外，系数值的获取是如何计算，自己手动计算还是直接通过软件监测？若通过软件监测，则需注意参考值的设定。

多面体明显质量好过关 过基础计算关 相对正六面体 同样参数控制 数量变小了

您好，我也在研究topset和refinemesh的问题，您知道在blockmesh生成一个梯形（具有斜面）的网格时，使用refinemesh会出现报错是什么原因吗？

求助各位老师，在扩展出口计算域计算时，发现了如下问题，具体细节如下图，不同方案中红色截面处的总压大小不同，是什么原因造成的。（尤其是方案2和方案3总压的不同是因为什么原因）
问.png