大涡模拟中如何产生湍流？

[stratification]

最近研究自由表面与湍流的相互作用，具体情况跟一个实验类似：一个水箱底部装有一个振栅，通过振栅振动产生湍流，以此来研究湍流与自由表面的运动。
许多用DNS的文章中都是在N-S方程中添加源项来产生湍流，不知道这种方法LES是否可用?此外是否有别的数值方法（适用于LES）产生或模拟这种zero-mean velocity的湍流

[李东岳]

OpenFOAM里面的DNSFoam是通过快速傅里叶变换产生的随机湍流场。

许多用DNS的文章中都是在N-S方程中添加源项来产生湍流

有链接么？

LES在网格足够细致之后趋向于DNS，或许也可以？不确定。

[stratification]

@李东岳 DNS的链接http://xinguo-fluid.com/uploads/3/1/0/6/3106919/guo_shen_-2010-_interaction_of_a_deformable_free_surface_with_statistically_steady_homogeneous_turbulence.pdf

另外，我也看了下dnsFoam.C的代码，确实在N-S方程中增加了外力项force

想问一下这一段代码是如何定义force的？

我认为这个force应该是用来维持湍流的并不是产生湍流，否则没有这一项的话湍流会被粘性耗散掉。如果初始的速度场没有用boxTurb产生随机场的话，而是均匀的速度为零的速度场的话，根本算不出湍流。不知我的理解是否正确？

[李东岳]

我认为这个force应该是用来维持湍流的并不是产生湍流，否则没有这一项的话湍流会被粘性耗散掉。如果初始的速度场没有用boxTurb产生随机场的话，而是均匀的速度为零的速度场的话，根本算不出湍流。

这个我还不确定，你可以不用boxturb，然后用dnsFoam算一下他那个tutorial看看有没有湍流出来。

想问一下这一段代码是如何定义force的？

目前只知道基于快速傅里叶变换，更深的代码还没看 :cheeky: