CFD中文网

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

拉格朗日一般要从粒子碰撞的角度来考虑破碎聚丙啥的

可参考该篇博文：http://liusuanyatong.blog.163.com/blog/static/2155951572016527422386/

太厉害了。。。。。。又做实验又写代码， 全才啊

打电话警告怕啥 先用再说

很有帮助。

0_1481532599819_7447183924@chatroom_1481532220845_75.png
AeroFoam: an accurate inviscid compressible solver for aerodynamic applications

多谢两位老师

请教各位，该怎么对边界附近的单元赋值，比如下面这个UDF.我想将一个区域的 C_UDMI(c, tt, 0)全部设置为0，然后边界附近的 C_UDMI(c, tt, 0)设置为1，该怎么操作？谢谢！

void set_zone(Thread* t) { cell_t c; real vol; Thread *tt = THREAD_SUPER_THREAD(t); begin_c_loop(c, t) { if (THREAD_ID(t) == 1047) { C_UDMI(c, tt, 0) = 1; } end_c_loop(c, lct) } void set_cells(Thread *t) { Thread *tt = THREAD_SUPER_THREAD(t); face_t f; Thread *t_nbr; cell_t c0, c1; Domain *mix_domain = Get_Domain(1); for (int i = 0; i < num_boundaries; i++) { t_nbr = Lookup_Thread(mix_domain, 1); begin_f_loop(f, t_nbr){ c0 = F_C0(f, t_nbr); c1 = F_C1(f, t_nbr); if (c0 != -1 && THREAD_ID(THREAD_T0(t_nbr)) == 1047) { C_UDMI(c0, tt, 1) = 0.0; } if (c1 != -1 && THREAD_ID(THREAD_T1(t_nbr)) == 1047) { C_UDMI(c1, tt, 1) = 0.0; } } end_f_loop(f, t_nbr) } } DEFINE_ADJUST(aaa, mix_domain) { Thread* mix_thread; int phase_domain_index = 0; thread_loop_c(mix_thread, mix_domain) { Thread* t = THREAD_SUB_THREAD(mix_thread, phase_domain_index); if (t != NULL) { set_capillary_pressure_sources(t); set_boundary_cells(t); } } }

https://f.guge.gq/

国内用户试试？

有个问题想问问大伙，如题

CFD界的竞赛：https://mp.weixin.qq.com/s/0PjOnv09LerabVSZyp-g9Q

@cfd-china Patient-specific computational fluid dynamics structured mesh
以上是关键字，搜一下应该就有

好像只要有网格尺寸分割切换的地方，中间过渡不够，都会出现这种问题，是不是还得尝试调整网格。

均匀各向同性湍流，其雷诺应力矩阵非对角项等于0 竖向非均匀各向异性湍流，比如槽道流（水平向均匀），只有$ \left \langle u_{1} u_{2} \right \rangle $有值，其他非对角项均为0。Pope的书中也给出了解释。 现在问题是，均匀各向异性湍流，其雷诺应力矩阵非对角项是否为0？如果参照Pope书中类似推导，在空间3个方向都是均匀的，理论上$ \left \langle u_{1} u_{2} \right \rangle $，$ \left \langle u_{1} u_{3} \right \rangle $，$ \left \langle u_{2} u_{3} \right \rangle $都为0吗？

《Pope, S.B., Turbulent Flows. Vol. 12. 2000: Cambridge University Press.》89页
e0bada21-7ac0-4ae4-acc3-3319a012dfc7-image.png

情况描述：一个叶片模型；操作压力设置为0，来流使用Pressure-inlet入口，总压9027.6Pa，总温299K；出口使用Pressure-outlet，使用静压出口。上下壁面使用平移周期边界，在壁面上的静压结果下壁面前缘发现高于9027.6Pa，来流Ma数大概0.3左右。
总压云图：272e8ab6-c0ca-40a6-a8ae-85289baeec81-b6de291b644f1eae2029ad5d035c97c.png
壁面静压输出结果：71799acb-1afc-4de1-9b0b-92d1214bef24-eb82740d448428364bb1ed13c86ea1d.jpg
虽然差距很小，但是却比来流总压还要高；
出问题的位置：
aaeda908-0e6a-4dd9-b51c-6c1a826e3282-f3816906934847d300a0cdfbdf0f738.png
还有很好奇的一点，理想气体的总压-静压公式应该是和Ma数相关的，在壁面上Ma=0，按理说总压应该等于静压，但是实际上并不是；
计算公式：839a2402-39a1-47b7-822c-dceda0b70b88-1f796c8655497c12e0a5455dd059f99.png
壁面总静压差△P：
d4869df2-580e-4051-ac51-1341ca071771-image.png
所以现在很疑惑。

z 和 x 用得最顺。

其他的还在熟悉中。

也不能说是复合，2016年Slash 51岁生日，Axl和Slash同台表演，时隔20年的分道扬镳。
http://www.gunsnroses.com/tour
0_1485104660529_20057521-mmmain.jpg

@chon01 这也只是我的猜测，因为我按了一下计算器，lg和ln函数算出来的结果差距确实比较明显，希望能帮到你吧。

对于这个问题，作为一个重度OpenFOAM v4.1, 轻度OpenFOAM v1606+,和目前轻度但在不远未来的重度OpenFOAM extend-4.1的使用者，我很想去了解这其中的一些渊源。我想从一个最普通用户的角度，来谈谈这个问题。

OpenFOAM v4.x 这个系列的最早的开创者是 Henry Weller. 来源 他在 1989年将FOAM,也就是OpenFOAM的前身创造成OpenFOAM，他被认为是OpenFOAM的创始者之一。目前这个系列的代码由 OpenFOAM Foundation 有限责任公司所推动。网站的页角列出了三名主要的董事会成员:

Henry Weller Chris Greenshields Cristel de Rouvray

OpenFOAM+ 系列的基地则对应这里。 它是由ESI公司的子公司OpenCFD有限责任公司所推动。

OpenFOAM extend 对应的是一个“延展”计划，它的基地在这里 目前的管理员主要有 来源：

Bernhard Gschaider Henrik Rusche Hrvoje Jasak Håkan Nilsson Martin Beaudoin Robert Keser

以上是仅列举了一些主要三大分支的一些基本信息。

============================================================================================

下面我想谈谈基于这些基本信息的个人感受：

个人来讲，我使用OpenFOAM的最主要原因是Licence的问题。Fluent的费用是无法承担的。 可以这么说，如果研究组里能用钱能解决Licence的问题，我恐怕很难会像今天这么依赖OpenFOAM。

Henry Weller对于我来说，一直是一个神秘一般的存在。作为OpenFOAM这样一款优秀的C++库的作者之一，我在互联网上很难找到有关他的个人信息。如果有一天，当面见到他，我恐怕都无法认出他。对Hrvoje Jasak 则是我在这个帖子中才听说的一个人，可能是由于我过分依赖OpenFOAM vx.x 系列。但是，我简单搜索就能找到关于Hrvoje Jasak的个人信息。我不想也不愿意做过多的猜测Henry Weller为什么如此“低调”。

另外，我想谈谈自己对开源精神的理解。

我觉得开源最核心的思想是这样：当一个东西被开源出来，它就不属于与任何一个人，而存在于人们共同的某种“想象”。人们觉得在自己力所能及的范围内，贡献出自己的东西，能够让这个“想象”更好，是推动这个“想象”不断前进的动力。至于推动的动机何在，不同的人有不同的理解。商业运作并没有被排斥在外。维基百科的理念，我个人觉得就是对开源精神的最好实现之一。但是人是复杂的，社会是复杂的。一个人的想法都会随着年龄段的不同而呈现不同想法，更何况千千万万人组成的社会呢？Linux各种各样的发行版，何尝不是这种对开源精神的不同解读而产生的结果。

============================================================================================