您意思是只需要和学校里的论文一样,进行数值准确性验证,就可以吗?
是的。不用自己做实验
您意思是只需要和学校里的论文一样,进行数值准确性验证,就可以吗?
是的。不用自己做实验
不重叠就得插值。不过感觉插值好像也不太合适。如果都LES了,插值有失精度。如果要插值,简单的就线性插值呗。高端点的来个2次3次样条啥的。具体差不插值,我也是自己琢磨的,或许也可以吧。
你可以首先反推出网格点的位置信息。然后把相关变量处理成场数据。比如数据是
x y u
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
2 3 4
那你的网格单元体心的位置坐标就是(2 3),然后U的值是(4 0 0)
数据应该是按块分布的,这五列应该都是变量x,数据可以在tecplot里重新输出为按列分布的,
没懂
https://turbulence.pha.jhu.edu/docs/channel/profiles.txt
这是Hopkins的数据,很直观:
y+ U+ uv+ uu+ vv+ ww+ tau_uv/tau_w tau_nu/tau_w tau/tau_w P+ pp+ up+
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
0.0000000e+00 0.0000000e+00 0.0000000e+00 4.6029311e-36 0.0000000e+00 1.2527210e-35 -0.0000000e+00 1.0000000e+00 1.0000000e+00 0.0000000e+00 6.4599422e+00 -2.8535219e-20
0.1769350 0.2091558 0.2413900 0.2736352 0.3058870
数据这5列都是什么?
NASA做的2D LES模拟?
坐标,时均速度、压力、雷诺应力场
简单来几行raw数据看看
CFD的科研完全可以不做实验去做研究。我做CFD这10多年,从来没做过实验。需要实验数据的,都是找sci里面别人的实验数据。
网上搜了一下,基本没有人讨论这个。记录一下,方便其他人能搜索到。
https://www.service-public.fr/particuliers/vosdroits/F2718
法国目前给未成年人颁发的居留卡是这个DECM。虽然DECM赋予未成年人从非申根地区,返回申根地区的权利。但实际操作过程中,很有可能无法入境除了法国之外的申根地区。个人经历:2024年从洛杉矶直飞西班牙,持有DECM,不给发登机牌。说不认识这个DECM。
个人第二个经历:从马尔代夫去巴塞罗那,购票后跟航空公司合适,航空公司说不认识法语的DECM,让我们跟海关确认。海关说这个需要跟航司确认。所以基本无解。有很大风险。
群友经历:
希望以后某个人google搜索的时候,能搜到。
@WeakForm Applied Mathematical Modelling 100 (2021) 170–191
你说的对的。CFD的文章都要验证,1)要么解析解,2)要么实验,3)要么其他算法的解。
如果不验证,说明不了算法的正确性
楼上完全正确
airbubbles=air1+air2
奇怪了,是显示的网格么?
这外形跟我这个很像 http://dyfluid.com/openfoam.html#id10
@范准 前几天收到个邮件,我说这个名字咋这么熟悉,想不起来在哪,原来在这里
估计得要真正30年CFD经验的人来看了 :-)
这个算例,DES的RANS区域,基本就在壁面网格那一层,不超过2层。
时均结果都会有些区别,应该是正常的吧。
我觉得实验、数值模拟,不管什么模型,时间平均后应该都是一样的
瞬态会不同
@biaowai 那个LES的结果都时间平均了
您那边有对比过RANS与LES计算的湍流粘度么?
没有。但是不应该会这么大。
看了我那个2D Smagorinsky跟RANS对不上是Smagorinsky的问题。你那个2D WALE能跟2D rans对上?
@剑金锋 在 pitzDaily算例LES与RANS的结果不一致 中说:
先算一段时间获得平均速度,获得u_avg, v_avg, w_avg,然后用u' =u-u_avg,v' = v- v_avg, w' = w -w_avg获得脉动速度。
看起来你求解的思路是对的。你确定LES的epsilon代码里面没有bug?