CFD中文网

不是很清楚你说的那种涡脱落现象。不过用户指南中我记得有一章曾经表示，大体意思为：“如果用稳态求解器模拟比较强的非稳态问题会出现残差的震荡”。

能否附稳态求解涡街的文献？谢谢。

@cfd-china link text!![alt text]0_1492654992588_QQ截图20170420102158.png 这是这个文献上说的在固壁上的条件；还有一篇文献直接在壁面设置为0了link text0_1492655090308_QQ截图20170420101901.png

thinkpad t480s i7 8550u 16G内存
平时写论文、做ppt用，也可以开虚拟机改求解器算算小规模问题，真要数值计算还得远程服务器:140:

各位老师，最近在研究双欧拉模型模拟液固两相流动，遇到以下两个疑惑，希望老师给予帮助:xiexie:
①想利用UDF对颗粒相剪切粘性系数和固相压力等参数进行修正（如26fbd023-2fb6-48f6-842a-733347fa08d1-图片.png ），公式中包含颗粒拟温度项，在编写UDF的过程中不知道如何获取该参数？？:134:
是否可以利用颗粒脉动速度C表示的颗粒拟温度7bf0c5a9-9020-484b-ba1e-2ba848de0117-图片.png 该方法中如何获取颗粒脉动速度？？？曾尝试用速度获取宏 C_U(c,t); C_V(c,t); C_W(c,t)得到颗粒速度，所得结果明显错误
②湍流模型中的湍流黏度是否要定义:141: 378f2ea3-f977-4a04-904a-05db21608467-图片.png
k-e湍流模型中湍流黏度定义为e94ab208-99c9-49db-9e65-664f5af4574b-图片.png
在输出湍流黏度时，是默认按照该公式进行计算的吗？？？只有有别于改公式时才需要写UDF定义吗:141:

是欧拉-欧拉，确实如果引入边界条件的话不知道压力场怎么处理。谢谢分析，我再想想。加多孔介质主要就是物理上好解释，我现在觉得可以在wall的区域再附加一个指向wall外面的源项（根据具体问题可以指定方向），类似重力，以此抵消重力以及其它数值问题造成的穿透问题，相当于在这个源项下，wall区域的固体会被很缓慢地排出wall。我估计一开始如果就加上的话，可能就可以阻止固相进入wall，而且有多孔介质巨大的阻力系数在，排出的固相应该速度非常小，对wall外面的计算应该影响有限。

这几天备课OKS，今年大幅度增加湍流的内容。前几天也在朋友圈发了个图。简单说一下。LES处理壁面几种方法：

LES-RANS混合方程：用RANS求解壁面区域，LES求解主流区。也就是你需要调用DES、DDES等之类

Wall-modelled LES：对壁面区域进行模化，主流区采用LES。这个进一步还区分不同的模型，比如代数类、PDE类、ODE类等。

http://www.dyfluid.cn/theory.pdf 混合方程我在这里更新了一点内容。

另外没理解你说的壁面雷诺数什么意思。

捕获.JPG

Ref: Journal of Computational Physics 305 (2016) 589–603

1.个人思路是在DPM中把雾滴设置成一个injection，粉尘颗粒设置成另一个injection，然后运算，想看到这个喷雾开启前后粉尘浓度的变化，但是冥冥之中感觉没有这么简单，希望各位老师能指点一下如何实现这个过程?
2.DEM可以模拟两种不同颗粒的碰撞过程吗？比如水滴和粉尘碰撞。

这个提醒就是有负质量了，或许是负网格造成的，check mesh 一下

感谢各位大佬！:chitang:

无意间发现这个帖子，里面提到了normal turbulent model,
https://www.cfd-online.com/Forums/openfoam/125473-low-reynolds-turbulence-models.html

@arctan 非常感谢！

@赵一铭 润滑理论的雷诺方程就是NS方程简化了，不考虑厚度方向的压力变化。因为在润滑中动静环之间的液膜很薄，可以认为在厚度方向没有压力变化，是雷诺提出的，就叫雷诺方程了。不知道和那个雷诺是不是一个雷诺。

@百善孝为先0626 可以的话还是想麻烦您把我的联系方式告知他，这边看能不能合作一下

首先感谢李老师的热情帮助~~
这个PPT中的第9页给出公式（1.5.9）。如果这个公式是正确的，则我的疑问就没有了，因为我原图中（4）式的右端项，第二项与第三项抵消了。但可疑的，这个公式（1.5.9）似乎并不正确。贴上来再麻烦李老师帮我看看：
0_1502119393399_7fd74579-5b34-4c15-ac3c-1cc3aae56236-image.png

@dyb 你可以参考自带的案例——airfoil2D，里面有升阻力系数监测的方法。

水中Re=40000，c=100mm，第一层网格高度给0.5mm以下，尾部圆弧边界层就挺好的，由y+算出来的y是1.25mm，尾端边界层就乱了
alt text

回复: 两个相切的圆如何进行关联
@kscandal
同学，我说的话，你确定，你懂？

问题背景是在simpleFoam计算中添加被动标量运输，因此在controlDict里添加了function的scalarTransport。但是在设置参数的时候发现resetOnStartUp。
在源代码中的discription是Reset scalar to zero on start-up。
现在比较疑惑的是，如果打开resetOnStartUp:
１．是在每个时间增量计算开始时把标量reset为０
２．还是一整个计算结束后，再次重新计算时把标量设为０
大家有经验吗

Von Neumann Stability Analysis 在普通的CFD教材里面讲的太简单了
但是到复杂非线性问题（比如Euler方程）中分析时间推进格式的阻尼特性 (damping characteritics) 的时候就比较复杂
而且一般期刊文章中只给应用结果而且缺乏细节
它背后的逻辑是什么，我一直没有搞清楚。。。
希望有大神能推荐一些详细一点的教材，能把思路讲的比较透彻的

@liqiwen 还可以对结构网格进行二次开发， [基于FLUENT软件的导弹意外点火情况下舱室维护仿真平台开发

](http://d.g.wanfangdata.com.cn/Thesis_D429353.aspx)这个论文的3.1.2节有讲用index来对ICEM CFD二次开发的功能。