CFD中文网

您好，我也遇到了相同的问题，请问您解决了吗

残差图IMG_20200601_201953.jpg

@zhangxc0223
Hi
因为并没有做deterministic的模拟，所以code也基本没看过。。
stochastic collision 在计算collision probablity的时候是有限制的，具体你可以看一下O‘roucke的论文。具体是怎么定义的，怎么得到一个non-dimensional collision probablity的，还有如何决定这两个parcel发生了碰撞。

虽然碰撞在这个instantaneously 随机的，但是在统计学上，比如你有1M次碰撞的话，总的碰撞的结果就是相似的。
这也就是为什么stochastic碰撞会比deterministic相对来说’便宜'的原因，通过统计学的方法达到相似的碰撞结果，但是不需要追踪所有的parcel。

lagnrangian maxCo 存在的意义，个人认为，在一定程度上可以理解为再一次离散了每个euler time step，所以maxCo一定程度上决定了每个lagrangian time step的大小，和Euler 的Co 定义是相似的。

碰撞是不会被错过的，因为deterministic理论上来说是追踪所有的parcel，所以下一步这个parcel会在哪，走的路径应该也是计算的，所以不会存在错过的，stochastic collision 会计算两个parcel的碰撞概率，当概率很高的时候，也是不会错过的，而且就算是这次错过了，在茫茫parcel海中，总会碰到的。。统计学上是不会错过的

以上仅为个人理解

各位老师好：
我用Fluent计算流道进出口之间的总压降，发现几个问题如下：

我设置压力入口总压为1MPa，压力出口静压为0，计算结束后发现出口的总压反而比入口的还高，在没有外界做功的情况下我认为这是不合理的，总压平白无故增加了，为什么？ 我详细看了看出口的压力分布，采用Mass-Weighted Average表面积分方法，发现出口处的总压和动压差距很大，按理来说，出口静压设置为0以后，总压和静压应该是相等的，但结果显示并不相等（见图1）。换用Area-Weighted Average方法或Facet Average 方法，问题依然存在，进一步导出出口的ASCII文件，发现出口处各单元的总压和动压也不相等（见图2）。总压≠静压+动压这个问题仅存在于计算结果中的出口处，入口完全没问题。 改变边界条件，换用速度入口和outflow出口的组合时上述问题2得以解决。但只要使用的是压力出口，不管是速度入口压力出口的组合，还是压力入口压力出口的组合，都存在这个问题。

Fluent的具体设置如图3所示，因为模型比较简单，收敛的也很快，计算在87次迭代后收敛，收敛情况及监控物理量如图4所示。
我想知道为什么会出现上述问题，出口的总压为什么会比入口高？出口处计算出来的压力和我预想的不一致的原因是什么？恳请各位老师指教。

d2425a1d-5980-4190-952f-e3b2e246058f-image.png
图 1 采用Mass-Weighted Average表面积分方法时的出口压力结果

9a171ad1-dde7-4ed6-a22e-5f3958cdd438-image.png
图 2 ASCII文件中的出口压力结果

7163498f-eed6-450d-8255-6d377f2c2630-image.png
图3 Fluent的具体设置

a703bcd8-ae07-4b9f-a66f-4ca91b040653-image.png
图4 计算收敛情况及监控物理量的变化情况

目前openfoam不支持这个。不太清楚fluent支不支持。持续保持关注。

感谢分享！！！收录在此 https://www.cfd-china.com/topic/3499/openfoam小代码

@东岳 又出现了owner和neighbour的概念，以前还以为和i,j,k一样的概念，看来又有知识点需要补咯，谢谢东岳老师

@huoyanxiadebing 在 说实话，根本不知道选择什么湍流模型 中说：

看到你的帖子，发现你做的都是外流场，对于外部流动，沿表面位置分布的雷诺数大于500000时，通常认为流动状态为湍流；沿障碍物的雷诺数大于20000时认为流动状态为湍流，如果你的模型中雷诺数不是很大的话，你用层流模型计算即可，这就是有时候你用层流模型算的比较准确的原因；对于雷诺数的计算，最重要的一点就是你要选好特征长度，对于没有结构物时，你要以水槽的长度或高度来算雷诺数，加上了结构物，就要判断是否产生了湍流现象，以结构物的某个长度及其附近的速度来计算雷诺数。

前辈您好，请问这里提到的两个临界值是否有相应文献的支撑？

如果要计算防波堤（例如沉箱防波堤和水下防波堤）附近的雷诺数来选取相应的湍流模型（适用于低雷诺数或高雷诺数的），这个特征长度应该取防波堤的高度吗？对于二维和三维情况，特征长度的选取是否有不同呢？

另外，结构物附近的速度是否有较为通用的估算方式？对于防波堤附近的速度，能否用基于浅水波理论的公式Uh=cη进行估算呢（水深h，波速c，自由面高度η）？

还有为什么在相同边界条件下，NIXTURE模型加表面张力和不加表面张力，收敛性差距很大？有人说MIXTURE可以不用开启表面张力正确吗？请问各位大佬！

这个估计很难了。你这个涉及流固耦合的问题。CFD这边做不鸟。但是你可以这么做，流场获得你需要的压强时程曲线，然后你在固体场中把压强加上去。这样是可以的。

问个小问题，就是平时后处理分析时说的某条线上的压降曲线，纵坐标指的是这条线上各点的压力值吗？还是这条线上各点的压力值与某个参考点的差值？

HTML文件，内容如下

<div component="composer" class="composer<!-- IF resizable --> resizable<!-- ENDIF resizable --><!-- IF !isTopicOrMain --> reply<!-- ENDIF !isTopicOrMain -->"> <div class="composer-container"> <nav class="navbar navbar-fixed-top mobile-navbar hidden-md hidden-lg"> <div class="btn-group"> <button class="btn btn-sm btn-primary composer-discard" data-action="discard" tabindex="-1"><i class="fa fa-times"></i></button> <button class="btn btn-sm btn-primary composer-minimize" data-action="minimize" tabindex="-1"><i class="fa fa-minus"></i></button> </div> <!-- IF isTopic --> <div class="category-name-container"> <span class="category-name"></span> <i class="fa fa-sort"></i> </div> <!-- ENDIF isTopic --> <!-- IF !isTopicOrMain --> <h4 class="title">[[topic:composer.replying_to, "{title}"]]</h4> <!-- ENDIF !isTopicOrMain --> <div class="btn-group"> <button class="btn btn-sm btn-primary composer-submit" data-action="post" tabindex="-1"><i class="fa fa-chevron-right"></i></button> </div> </nav> <div class="row title-container"> <!-- IF showHandleInput --> <div data-component="composer/handle"> <input class="handle form-control" type="text" tabindex="1" placeholder="[[topic:composer.handle_placeholder]]" value="{handle}" /> </div> <!-- ENDIF showHandleInput --> <div data-component="composer/title"> <!-- IF isTopicOrMain --> <input class="title form-control" type="text" tabindex="1" placeholder="[[topic:composer.title_placeholder]]" value="{title}"/> <!-- ELSE --> <span class="title form-control">[[topic:composer.replying_to, "{title}"]]</span> <!-- ENDIF isTopicOrMain --> <ul class="dropdown-menu quick-search-results hidden"> {{{each posts}}} <li> <a href="{config.relative_path}/post/{posts.pid}"> {buildAvatar(posts.user, "sm", true)} <span class="quick-search-title">{posts.topic.title}</span> <br/> <p class="snippet"> {posts.snippet} </p> <small class="post-info pull-right"> <span class="fa-stack" style="{function.generateCategoryBackground, posts.category}"><i style="color:{posts.category.color};" class="fa {posts.category.icon} fa-stack-1x"></i></span> {posts.category.name} &bull; <span class="timeago" title="{posts.timestampISO}"></span> </small> </a> </li> <!-- IF !@last --> <li role="separator" class="divider"></li> <!-- ENDIF --> {{{end}}} <!-- IF multiplePages --> <li class="text-center"> <a href="{url}"> [[search:see-more-results, {matchCount}]] </a> </li> <!-- ENDIF multiplePages --> </ul> </div> <!-- IF isTopic --> <div class="category-list-container hidden-sm hidden-xs"></div> <!-- ENDIF isTopic --> <div class="pull-right draft-icon hidden-xs hidden-sm"></div> <div class="btn-group pull-right action-bar hidden-sm hidden-xs"> <button class="btn btn-default composer-discard" data-action="discard" tabindex="-1"><i class="fa fa-times"></i> [[topic:composer.discard]]</button> <button class="btn btn-primary composer-submit" data-action="post" tabindex="6"><i class="fa fa-check"></i> [[topic:composer.submit]]</button> </div> </div> <div class="category-tag-row"> <div class="btn-toolbar formatting-bar"> <ul class="formatting-group"> <!-- BEGIN formatting --> <!-- IF formatting.spacer --> <li class="spacer"></li> <!-- ELSE --> <!-- IF !formatting.mobile --> <li tabindex="-1" data-format="{formatting.name}" title="{formatting.title}"><i class="{formatting.className}"></i></li> <!-- ENDIF !formatting.mobile --> <!-- ENDIF formatting.spacer --> <!-- END formatting --> <!--[if gte IE 9]><!--> <!-- IF privileges.upload:post:image --> <li class="img-upload-btn hide" data-format="picture" tabindex="-1" title="[[modules:composer.upload-picture]]"> <i class="fa fa-file-image-o"></i> </li> <!-- ENDIF privileges.upload:post:image --> <!-- IF privileges.upload:post:file --> <li class="file-upload-btn hide" data-format="upload" tabindex="-1" title="[[modules:composer.upload-file]]"> <span class="fa-stack"> <i class="fa fa-file-o fa-stack-1x"></i> <i class="fa fa-arrow-up fa-stack-1x"></i> </span> </li> <!-- ENDIF privileges.upload:post:file --> <!--<![endif]--> <!-- IF allowTopicsThumbnail --> <li tabindex="-1"> <i class="fa fa-th-large topic-thumb-btn topic-thumb-toggle-btn hide" title="[[topic:composer.thumb_title]]"></i> </li> <div class="topic-thumb-container center-block hide"> <form id="thumbForm" method="post" class="topic-thumb-form form-inline" enctype="multipart/form-data"> <img class="topic-thumb-preview"></img> <div class="form-group"> <label for="topic-thumb-url">[[topic:composer.thumb_url_label]]</label> <input type="text" id="topic-thumb-url" class="form-control" placeholder="[[topic:composer.thumb_url_placeholder]]" /> </div> <div class="form-group"> <label for="topic-thumb-file">[[topic:composer.thumb_file_label]]</label> <input type="file" id="topic-thumb-file" class="form-control" /> </div> <div class="form-group topic-thumb-ctrl"> <i class="fa fa-spinner fa-spin hide topic-thumb-spinner" title="[[topic:composer.uploading]]"></i> <i class="fa fa-times topic-thumb-btn hide topic-thumb-clear-btn" title="[[topic:composer.thumb_remove]]"></i> </div> </form> </div> <!-- ENDIF allowTopicsThumbnail --> <form id="fileForm" method="post" enctype="multipart/form-data"> <!--[if gte IE 9]><!--> <input type="file" id="files" name="files[]" multiple class="gte-ie9 hide"/> <!--<![endif]--> <!--[if lt IE 9]> <input type="file" id="files" name="files[]" class="lt-ie9 hide" value="Upload"/> <![endif]--> </form> </ul> </div> </div> <div class="row write-preview-container"> <div class="write-container"> <div class="help-text"> <span class="help hidden">[[modules:composer.compose]] <i class="fa fa-question-circle"></i></span> <span class="toggle-preview hide">[[modules:composer.show_preview]]</span> </div> <div class="pull-right draft-icon hidden-md hidden-lg"></div> <textarea class="write" tabindex="4"></textarea> </div> <div class="hidden-sm hidden-xs preview-container"> <div class="help-text"> <span class="toggle-preview">[[modules:composer.hide_preview]]</span> </div> <div class="preview well"></div> </div> </div> <!-- IF isTopicOrMain --> <div class="tag-row"> <div class="tags-container"> <div class="btn-group dropup <!-- IF !tagWhitelist.length -->hidden<!-- ENDIF !tagWhitelist.length -->" component="composer/tag/dropdown"> <button class="btn btn-default dropdown-toggle" data-toggle="dropdown" type="button"> <span class="visible-sm-inline visible-md-inline visible-lg-inline"><i class="fa fa-tags"></i></span> <span class="caret"></span> </button> <ul class="dropdown-menu"> <!-- BEGIN tagWhitelist --> <li data-tag="@value"><a href="#">@value</a></li> <!-- END tagWhitelist --> </ul> </div> <input class="tags" type="text" class="form-control" placeholder="[[tags:enter_tags_here, {minimumTagLength}, {maximumTagLength}]]" tabindex="5"/> </div> </div> <!-- ENDIF isTopicOrMain --> <!-- IF isTopic --> <ul class="category-selector visible-xs visible-sm"> </ul> <!-- ENDIF isTopic --> <div class="imagedrop"><div>[[topic:composer.drag_and_drop_images]]</div></div> <div class="resizer"><div class="trigger text-center"><i class="fa"></i></div></div> </div> </div>

86f6fd36-007a-42a5-b105-de5d4fcefffa-image.png
cited:http://openfoamwiki.net/index.php/See_the_MRF_development

图为在旋转坐标系内解绝对速度的控制方程，其中 \vec \omega \times \vec U_i 这项称为科氏力，旋转坐标系内独有的力。他的效果如下面这个视频所示https://haokan.baidu.com/v?vid=4227630901810560602&pd=bjh&fr=bjhauthor&type=video
我的理解是可以将上面的式子理解为绝对速度产生组成的动量的输运方程，其中的通量是相对通量（\vec U_R \vec U_I）,同时由于动量是一个矢量，因此在旋转坐标系中还受到科氏力的作用。等式坐标是压力的贡献，和由于运输的物理量U_I的扩散贡献。

因此，我在思考是不是诸如湍动能，湍流耗散率这些标量，也可以认为在旋转坐标系内，运输物理量为湍动能，湍流耗散率这些标量，通量为相对速度的通量，但是由于这些量都是标量，因此没有多的项需要添加。另外比如SA模型中的发生项需要计算到voricity的大小，此处的voricity是不是也应该用相对速度？

是欧拉-欧拉，确实如果引入边界条件的话不知道压力场怎么处理。谢谢分析，我再想想。加多孔介质主要就是物理上好解释，我现在觉得可以在wall的区域再附加一个指向wall外面的源项（根据具体问题可以指定方向），类似重力，以此抵消重力以及其它数值问题造成的穿透问题，相当于在这个源项下，wall区域的固体会被很缓慢地排出wall。我估计一开始如果就加上的话，可能就可以阻止固相进入wall，而且有多孔介质巨大的阻力系数在，排出的固相应该速度非常小，对wall外面的计算应该影响有限。

@东岳 这位是我后来一步一步检查，用checkMesh得到得具体的问题
30bf2ade-fc6e-412f-9486-f0083a60212b-image.png
如果感觉没有那么复杂的不规则边界，那么基本就是几何拓扑关系出了问题吧

在做单气泡尺寸的时候我都是用db = C_PHASE_DIAMETER(cell, thread_g)来获取气泡直径，最近想尝试双气泡直径，但是在参考的曳力计算式中有大气泡分数（fraction of large bubbles）这一项，不知道这项的值应该用什么程序句子来获取呢？

还请大佬能帮忙解惑，万分感谢！

ok，感谢。

请教大家一个问题，笛卡尔网格，已知网格8个顶点的坐标，以及网格中心的物理量的值，但网格排序是杂乱无章的，有什么好方法进行后处理呀？

@TINGHAIK 谢谢您，问题已解决