@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