CFD中文网

@李东岳 不是，虽然帖子是在openfoam，但是不是用这个软件

是合作课题组开发的程序基础上计算的，不知道大佬有没有什么建议可以解决现在的问题？

@李东岳 网格是512x16x512, 对应的尺寸是52 x 2 x 52, 展向是周期边界，所以对应的尺寸和网格少一些，这些是目前测试采用的网格，求解是采用的自行编译的程序，是直接求解NS方程。。。现在遇到的问题就是轴线流向速度衰减太快，射流出口边界处的流向速度比文献数据要低一些。。。

@cccrrryyy 对比的是一篇DNS的数据，我已经更换了入口处的速度曲线，结果还是没有改善。。也尝试过选择其他的求解方式，但是发现轴线速度随流向演变，衰减的特别快，这同时也是造成展向结果对不上的原因。。。还没有搞清楚应该如何调整 跪求大佬们指点

@cccrrryyy 大佬是觉得我修改压力出口边界还是不太好吗？我现在在想可能是出口边界与法向边界的边界条件导致流场数据对不上，我已经试过了发现法向壁面如果采用0梯度的话同样会导致回流。。所以法向壁面如果不采用0梯度的话，应该对u施加什么约束比较好呢？（如果是给coflow0.1，会不会造成拖拽边壁流体向前，给无滑移的话会不会也是同样的结果？）出口也尝试了advective边界条件，存在相同的问题。。。。所以目前怀疑是出口与法向壁面边界条件组合导致的问题

最近在编程求解平面射流的DNS，因为受限于计算量，所以设置计算域流向长度为52d，法向为52d（d为射流入口直径），现在的边界条件如下：
压力边界：
流向：入口0梯度边界，出口设定为狄利克雷边界（0），为了防止回流，在出口处设置压力为-0.5*(u^2+v^2+w^2)(0.5(1-tanh(u)))
展向：周期边界
法向：狄利克雷边界（0）

速度边界：
流向：入口恒定速度曲线，出口0梯度
展向：周期边界
法向：上下壁面采用无滑移，对u来说考虑射流入口存在coflow，所以壁面u速度设定为coflow速度（实际文献中多采用0梯度边界，但是在实际模拟中发现会导致回流，这个点上也想请问各位大佬）

在修改出口处压力之后，回流问题减轻，已经可以正常流动，但是发现流场中心线速度衰减过快，不知道是什么问题？


上面是统计一小段时间的结果，分别是瞬时速度图和轴线流向速度，求大佬解惑！

@尚善若水 感谢回复！无反射边界也在一些文献里看到过，想问不可压缩流体也可以使用无反射边界吗？

@李东岳 感谢大佬的意见，我会试下大佬提到的这些边界条件！

@楼道男孩 感谢大佬的回复，因为计算量太大，目前没有考虑延长计算域，希望能够通过改变边界条件的方式来达到一个比较合理准确的流场。。。不知道大佬有没有什么建议吗关于边界条件这块

是不可压缩的流场，参数都是无量纲的

各位大佬，我最近在用DNS程序在求解平面射流，因为是DNS程序，为尽量减少计算量，所以计算域设置为（52d x 52d x 12d），分别为流向，法向，展向，d为射流入口直径，d为1，速度入口为1
速度边界条件设置为：
流向：恒定速度入口（流向速度为1，展向和法向为0），出口设置为0梯度（或二阶梯度为0）
法向：流向速度与展向速度设置为0梯度，法向速度设置为无滑移
展向：周期性边界条件
压力边界条件设置为：
流向：入口出口都为诺伊曼边界条件（0梯度）
法向：诺伊曼边界条件（0梯度）
展向：周期性边界条件
因为还在测试求解器阶段，所以实际在求解过程中没有给展向速度，可以视为是2维流场，但是求解中发现问题一般是两个，一个是流体流到出口位置发散（感觉是流体没有流不出去），一个是回流严重，最终会在出口处发散。。。
各位大佬有没有了解的，希望可以点拨一下

IMG_6563.PNG

这个case使用不同模型计算出来的结果很不一样，怀疑是模型的问题。。。

请问有没有哪位大佬使用过SA模型计算轴对称射流，在文献中，SA模型射流计算得到的spreading rate为0.248，但是我选用该模型计算，发现流场与射流结果有偏差，在60d后射流速度半宽不再呈线性变化，前部分的spreading rate也增加到了0.3附近。。。想问下哪位大佬了解过SA模型计算射流的case

使用fluent进行射流的计算，在时间步长1e-3下计算到100多秒时，残差降到1e-5以下已经很久，保持稳定，此时计算得到射流的扩张速率为0.096，但是后面继续去算这个case发现，在200多秒，400多秒后残差还是在1e-5以下，但是扩张速率却升高了到0.11.。。。。
本来已经收敛了结果，为什么流场数据还在变化？
有没有哪位大佬遇到过这种事

@veen 谢谢，我学习一下

槽道流使用fluent数值模拟，2d模型，基于标准k-ε湍流模型，考虑使用增强壁面函数，选择速度边界条件，想咨询一下各位大神，速度入口边界条件这里对应的湍流强度和长度尺度应该怎么设置？使用增强壁面函数对应处理近壁面，这里选择是否合理呢，对于槽道流，壁面函数这样选择是否合理？

@东岳 因为时间成本问题吧，现阶段还没有这方面的打算，之后应该会学。。。其实只是想问下这个求解器，如果简单方便入手的话想要了解一下。。。

@东岳 谢谢大神的回复，因为涉及到基于湍流模型的数值模拟，需要进行验证，所以在找这方面的软件，fluent的验证更改模型参数后并没有表现得如预期的结果

我查阅了相关文献，发现一部分基于湍流方程模型的计算优化并不是在ansys下fluent完成的，而是基于Wind-US Flow Solver这个软件进行完成，想请教一下有没有大神了解过这个软件~先在此感谢各路大神！

@楼道男孩 谢谢东岳大神的回复，我已经找到了问题的源头，文献中对应的是低雷诺数下ke模型，并不是标准的ke模型，是我之前知识盲区，导致一直没有找到解决的办法，现在正在尝试进行模拟验证，谢谢大家的关注！

@东岳 我又仔细对比了在两组模型参数下的计算结果，发现在使用了增强壁面函数后标准ke湍流模型更改参数之后整个流场的数据并没有很大变化，意思是使用增强壁面函数后更改模型参数对整个流场不起作用，想问问大神这个是因为什么导致的呢？