CFD中文网

欧拉欧拉那面由于已经在统计学方面已经处理过了，没有噪音。欧拉拉格朗日模拟这面目前有什么文章在处理么？找了找，还没看到。

As an alternative, Lagrangian methods ‘‘discretize” the density function into ‘‘parcels” that are simulated using Monte-Carlo methods. While quite accurate, as in any statistical approach, Lagrangian methods require a relatively large number of parcels to control statistical noise, and thus are computationally expensive.

@李东岳 大神，这个我真没注意。好像不多。我是土木专业的，不是搞固体力学的。您可以再找找。

在1954年，IBM创立了Fortran语言。

在1955年，多相流中的粒子法就开始进行了。在当时，其被取名为Particle In Cell (PIC)。

PIC由Los Alamos国家实验室的F.H. Harlow提出。在PIC中，拉格朗日粒子通过其坐标（x, y, z）进行定位并且具有质量。有些情况下，拉格朗日粒子还具有多个物理特征。例如在PIC中，拉格朗日粒子可能还具有大小的区别等。

特征长度弦长

@wwzhao

简单的感受了一下两种赋值方式，效率确实很大提高。

#include <stdio.h> #include <string> #include <chrono> #include <iostream> using namespace std ; class Timer { private: // Type aliases to make accessing nested type easier using clock_t = std::chrono::high_resolution_clock; using second_t = std::chrono::duration<double, std::ratio<1> >; std::chrono::time_point<clock_t> m_beg; public: Timer() : m_beg(clock_t::now()) { } void reset() { m_beg = clock_t::now(); } double elapsed() const { return std::chrono::duration_cast<second_t>(clock_t::now() - m_beg).count(); } }; class SomethingBefore11 { private: int m_array[5]; public: SomethingBefore11() // zero the member array { // If we want the array to have values, we'll have to use assignment here m_array[0] = 1; m_array[1] = 2; m_array[2] = 3; m_array[3] = 4; m_array[4] = 5; } }; class SomethingAfter11 { private: int m_array[5]; public: SomethingAfter11(): m_array { 1, 2, 3, 4, 5 } //zero the member array { } }; int main(){ Timer tBefore11; SomethingBefore11 m_array_before11; std::cout << "Time elapsed: " << tBefore11.elapsed() << ‘n’; Timer tAfter11; SomethingAfter11 m_array_after11; std::cout << "Time elapsed: " << tAfter11.elapsed() << ‘n’; return 0; } [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.506e-06 Time elapsed: 1.47e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 8.664e-06 Time elapsed: 1.9e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.646e-06 Time elapsed: 1.89e-07

装ANSYS18不就可以了

@bestucan
感谢您的回复！“不同方向的流动计算差异就大”这个解释很直接明了，我突然就有点明白这个网格线出现的原因了:xinxin:

@shrine 没关系，无关紧要这些文件

OpenFOAM v2006中，adiabaticPerfectFluid状态方程如下：

$\rho (p)=\rho_{ref}(\frac{p+B}{p_{ref}+B})^{\frac{1}{n_{T}}}$

根据熵的计算公式：

$dS=c_{p}\frac{dT}{T}-(\frac{\partial v}{\partial T})_{p}dp$

求解器中将这一公式分为两部分进行积分计算，其中第一部分$c_{p}\frac{dT}{T}$ 的积分计算在thermophysicalModels/specie/thermo/hConst/hConstThermoI.H中实现，即$c_{p}*log(\frac{T}{T_{std}})$
但第二部分的积分，由于公式是绝热状态的，没有温度作为参数，想请教一下官方求解器中为什么会进行以下处理：

通过源代码公式的求导反推，官方求解器在计算第二部分积分时，在状态方程的分母上直接添加了温度T，如下式：
CodeCogsEqn (1).gif
公式源代码详见adiabaticPerfectFluidI.H文件的S函数：

template<class Specie> inline Foam::scalar Foam::adiabaticperfectFluid<specie>::S ( scalar p, scalar T ) const { scalar n=1-1.0/gamma_; return -pow(p0_+B_,1.0/gamma_)*(pow((p+B_),n)-pow((Pstd+B_),n)) /(rho0_*T*n); }

除上述绝热状态方程外，线性变化的状态方程linear也采用了同样的处理方式，想知道为什么可以这样处理？
另外一个问题是，如果采用上述处理之后，对于气体的$c_{p}-c_{v}$，根据公式:
CodeCogsEqn (2).gif 该如何处理？

@李东岳 感谢老师 发现是我Co计算公式写错辽~已经修改并且成功运行了 谢谢~

@下里巴人 看到这个帖子想起我之前做那个流化床模拟，加了跟管子就死活不收敛了，估计和你这个一样的问题。
话说有尺度差异的时候出来的这个网格好呆……明明圆形的楞给画成多边形了。
我想起Gmsh他们组最近有个博士答辩的东西就是做curve mesh的，不知道能不能用在计算上了。。

我的理解是

Advection Scheme 是动量方程对流项的离散格式 Turbulence Numerics 是湍流输运方程对流项的离散格式

See "15.5.5. Advection Scheme for Turbulence Equations"

Irrespective of the Advection Scheme setting, the advection scheme for the turbulence model equations depends on the Turbulence Numerics option: ...

Note that Turbulence Numerics option also specifies the transient scheme for the turbulence model equations (see Turbulence Numerics in the CFX-Pre User's Guide). This can be overridden in CFX-Pre in the Solver Control details view on the Equation Class Settings tab.

你的图就是节点对齐。

有文献里是中心六面体网格，边界层为四面体问题网格。

你可以试试hex-core，在做体网格的时候，可以成圣hex-core

有的时候icem就是会报错，只能没做一些操作之后多保存，然后再进行。。

@西湖冷月 借帖求教大佬，狭窄裂缝中的湍流用哪种方法模拟好呢？

http://blog.sina.com.cn/s/blog_5fdfa7e601010rkx.html 苏老师的博客

@李东岳 非常感谢。
又读了几篇类似的文献，真是用时均处理的仿真模型找不到，高雷诺数下的这种流动都是在用实验测试了，但是我需要计算一个比较准确的结果，合作单位的同志希望能找到统一的规律，形成经验公式类似的结论，然后用作产品开发:crying: 任重道远啊，有个师兄推荐用DES试试，目前在测试，希望别掉进坑了。

@李东岳 ansys18版本 linux
链接：https://pan.baidu.com/s/1QpXN_GhW5tCv1rNJrlBUYA
提取码：1234

没办法，挪走吧，找个高性能路由器，远程连。服务器又不能自己拆😂

@浪迹天大 非常感谢！今天仔细看了下，终于搞定了。