CFD中文网

:wolaile: :wolaile: :wolaile: :wolaile: :wolaile:

:xiabanle: :xiabanle: :xiabanle: :xiabanle:

:tishizi:

@东岳 谢谢 岳老师

loop的是面网格上的所有网格。你的需求应该能实现，最简单的办法是，每个面网格中心去读一遍文档，把散点距离自己最近的那个点的速度赋给所在网格。但是会有其他问题，散点疏密是不是和面网格疏密比较匹配，也许一个面网格配多个散点或者反之。

@桎梏 那青轴可以说是爆炸 :xinxin:

鼓励鼓励！:high:

@李东岳 车灯热分析。

国内形势变化太快了，一年一年的水涨船高，就像我们看现在的教授，当年竟然博士毕业就留校了。以后的年青一代看我们，估计是“这家伙当年竟然博士发了5篇文章就留校了，现在我们博士发10篇都留不下来”..

@mohui 您也用interfoam算过水渠流动问题吗

@浪迹天大 在mathjax里面要设置跳过某些关键词，比如skipTags: ["<pre>","<span>","<code>"]

谢谢各位，解决了:mihu:

请参考无痛苦ns方程笔记的这一节：R、sigma、prime2Mean、以及雷诺应力相关量

@东岳 好的谢谢！
我想问下圆周沿线要怎么画？:chouchou: ，另外流向的平均就是流向沿线每个点的u+相加起来然后除以点数吗？
那y方向上的位置要选多少个？谢谢！

我这边需要仿真一个构件在烘箱中逐渐被加热的过程，结构非常简单，就是一个构件周围用壁面包裹，其截面如下图所示。
IGBT fixture-温度2.jpg
四周壁面设定为固定温度358K，初始温度为293K。由于工作中计算资源有限，且目前只需要定性得到构件各部位升温的快慢情况，因此仿真时只打开了能量方程，上图就是仿真至380秒左右时的状态。

然后在仿真中，初始的时间步长是0.001s，因为太小了导致仿真进度比较慢，而且能量方程的残差一直保持在10e-7级别，我这边通过不断地放大时间步长，至380s时时间步长为0.1s，至400s左右时时间步长改为1s。在这个过程中，我发现每次增大时间步长，监测点的温升曲线斜率都会变小，也就是说温升变慢了，比如下图就是时间步长由0.1s改为1s时，监测点温升情况的变化：
IGBT fixture-温升.jpg

在我更改时间步长的过程中，方程的残差均保持在10e-7级别，那么我想知道，温升速率为什么会因为时间步长的增加而变小呢？是不是因为网格不够精细造成的，因为为了快速出结果，网格确实比较粗糙，还是说还有其他可能的因素？谢谢大家了

@cfd中文网 东岳老师您好，我对第三点不是很明白。不太理解这个源项是什么意思，是数值耗散项吗？（如果是数值耗散项的话，这是什么原因产生的呢）可以有具体一点的方程可以参考吗？

http://www.mafengwo.cn/gonglve/ziyouxing/99502.html

就是如何从LES计算epsilon吧？目前据我所知只能通过epsilon的定义来求。你需要写个后处理程序

@东岳 是的，我看的文献大多数是关注速度，温度和涡量场变化，几乎没有涉及压力场变化，于是我编程计算的时候也没有考虑，结果老师问我为什么，我都不知道如何回答了:xinlei:

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

采用interFoam修改的求解器计算算例，遇到网格太小带来的壁面处压力问题（图中为第一个输出时间步（约0.0001 s）的压力分布，出现了不正常最大和最小压力）。
算例的壁面（BOTTOM）为锯齿状壁面，图中所示部分为倾斜向上的一小部分，壁面附近的网格精度约为1微米，压力场初始为0，压力和速度壁面处边界条件如下：
0001.png
p_rgh

BOTTOM { type fixedFluxPressure; value uniform 0; }

U

BOTTOM { type fixedValue; value uniform (0 0 0); }

大家有没有遇到过这样的问题，有什么解决思路吗？

需要一些大流量流量计300Nm3/h的，网上能找到的都是一些小品牌的，有没有推荐的，热式、涡街、压差都可以