CFD中文网

@李东岳 谢谢李老师

好东西好东西！！

根据横向力做傅里叶变化即可求出频率

我觉得STARCCM+的网格功能挺6的，笛卡尔/多面体等，控制的参数有限，基本上零基础跟着教程做几遍就差不多了，网格生成速度也很快，比较复杂的几何还是很容易解决的；
其他的用point wise多一点，复杂的几何自己划分块有点麻烦，如果不能分层很容易混乱，简单几何还是非常便利的；ICEM复杂几何切起来还是比较舒心的；^_^

@bestucan 太感谢了:xiexie:

下图的变量越界了。红的过了上界。蓝色过了下界。
解得纯对流方程，我这几年...一直在和有界和守恒作斗争 :big_mouth: :big_mouth:

......

0_1512380644852_捕获.JPG

\begin{equation}
m_0=w_1+w_2
\end{equation}
\begin{equation}
m_1=z_1w_1\xi_1+z_2w_2\xi_2
\end{equation}
\begin{equation}
m_2=z_1^4w_1\xi_1^2+z_2^4w_2\xi_2^2
\end{equation}
\begin{equation}
m_3=z_1^9w_1\xi_1^3+z_2^9w_2\xi_2^3
\end{equation}
\begin{equation}
m_4=z_1^{16}w_1\xi_1^4+z_2^{16}w_2\xi_2^4
\end{equation}
\begin{equation}
m_5=z_1^{25}w_1\xi_1^5+z_2^{25}w_2\xi_2^5
\end{equation}
左边的$m$都已知，右边的6个未知数$z_1$, $z_2$, $w_1$, $w_2$, $\xi_1$, $\xi_2$需要求。

除了Newton-Raphson迭代方法？

@李东岳 您好，抱歉，很久没来了。当时在学习FVM，最近在尝试利用openfoam解决NS方程。

首先是CPU 但是CPU的优化性和效果是非常复杂的问题 一般来说 目前的主流 最好最大主频不要低于2.8 同时并行核心数不要低于32 从我的经验48核心比较好
其次是内存 内存对一些内存要求多的算例的提升效果比较明显 能够占到30%左右
硬盘一般跟软件面板操作相关 有的时候也是一些时间
显卡的话 如果你是Fluent可以适用CUDA的方法 提升效果比较明显 两张4W的卡 大概coupled的下 提升30%左右
最后就是操作系统 你的多核心下 有些操作系统能够提升速度几百倍（显然有些系统对于40以上的核心效果很差） 一般Linux下比同等的windows server要快5%

@youmengtian
你用的什么显示器 :cheeky:

@bestucan 学到了，感谢

@李东岳 666:sunglasses:

http://www.sciencedirect.com/science/article/pii/S0021999117301493

2011年Vikas提出来一种稳定的高阶矩对流格式之后，本文作者重新提出来一个新的数值格式。

新技术刚追完，马上又出新的。。

是否为容差 matchTolerance 太大导致的衔接不上~

@BlookCFD 解决了，受到下面代码的启发，即重复利用代码，更新边界名称即可

p { solver GAMG; tolerance 1e-6; relTol 0; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; cacheAgglomeration on; agglomerator faceAreaPair; nCellsInCoarsestLevel 100; mergeLevels 1; } pFinal { $p; relTol 0; } 最后的要实现结果如下图： forceCoeffs_object { // rhoInf - reference density // CofR - Centre of rotation // dragDir - Direction of drag coefficient // liftDir - Direction of lift coefficient // pitchAxis - Pitching moment axis // magUinf - free stream velocity magnitude // lRef - reference length // Aref - reference area type forceCoeffs; functionObjectLibs ("libforces.so"); //patches ("body1" "body2" "body3"); patches (C0); pName p; Uname U; rho rhoInf; rhoInf 1.0; //// Dump to file log true; CofR (0.0 0 0); liftDir (0 1 0); dragDir (1 0 0); pitchAxis (0 0 1); magUInf 1.0; lRef 1.0; // reference lenght for moments!!! Aref 2.0; // reference area 1 for 2d writeControl timeStep; writeInterval 1; } //监测新的边界条件，只需要重复利用代码，更新边界名称即可 C1 { $forceCoeffs_object patches (C1); }

@anubis 集合重构不求解方程，而是通过通量变化计算体积分数，最后再进行界面重构。得到的体积分数必然是间断（指界面处发生跃变）的一个场。

搞定，输出叶片上的压力数据（不要选择internal mesh），用slice切一下，最后用plot data即可

@liujunCFD

icemcfdbat那个文件在的，但是不知道为啥就是打不开了

小骗叽:136: :136: :136:

没币发帖