CFD中文网

这个主要看cfd之道公众号流沙的教程就可以学会了（哈哈哈哈，你是我拿到注册码之后回复的第一个人）

Examinination of bubble collisions and coalescence in bubbly flows
D. Broder and M. Sommerfeld
还有这一篇，这好像都是小的那个做的，大的具体尺寸我也不知道，不过什么时候能在做实验不好说，实验仪器快被转移了，大的那个。。

对简单的长方体进行网格划分，再检查网格质量时发现在长方体的周围有很多有问题的face产生，也不知原因，很疑惑，求大神指点，万分感谢！

0T5OC4QVS4OWJO998JSUAU7.png

@bestucan 不太懂这个，就是用google浏览器打开网页，之前都没问题，最近突然出现这种打不开的情况。:136:

报错了，我检查不出错误，希望有前辈能够指导一二

我想建立一个水槽，水流从水槽一边冲过去，然后又从另一边冲过来，两边的入口和出口边界要随时间变化，有没有大神指导一下，怎么实现？谢谢！

@leo_chou 看过一个挺不错的视频：https://www.bilibili.com/video/BV1Yt411W7td

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>

以前简单看过，好像是研究表面活性剂的方法。大概是在三维曲面上用类似FVM的方法离散控制方程并求解。

snappyHexMesh边界层这个问题有的时候确实不太好处理。目前我也没有什么太好的建议。看看有没有其他大佬关注关注 :jingya:

如果只是查看，emacs 有个 whitespace mode。
空格、tab、换行之类的特殊字符都可以显示
vim按说应该也有，没搜到

Screenshot from 2021-05-22 15-16-48.png

各位大佬好，我是CFDEM新手，最近尝试参照twoSpheresGlowinskiMPI构建自己的算例，采用dynamicRefine，在计算过程中断构建的算例基本与参照算例一致，但试跑算例出现中断，不知道是否是哪里没有注意到，求大佬指导
报错描述
当不使用自适应网格时，算例正常计算；当采用自适应网格功能interDyMFoam，网格未加密前（未达到设置加密条件），正常计算;网格开始加密计算中断。
报错代码如下：

Time = 0.0001 Selected 0 cells for refinement out of 1613716. Selected 0 split points out of a possible 0. Courant Number mean: 0.000199582 max: 0.0804916 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' Setting up run at Fri Jul 26 12:59:28 2024 Memory usage per processor = 6.75634 Mbytes Step Atoms KinEng rke Volume dragtota dragtota dragtota 1 1 2.5261073e-08 0 0.0260508 0 0 0 CFD Coupling established at step 10 11 1 3.5918126e-08 0 0.0260508 0 0 0 Loop time of 0.000962037 on 16 procs for 10 steps with 1 atoms, finish time Fri Jul 26 12:59:28 2024 Pair time (%) = 4.70437e-06 (0.489001) Neigh time (%) = 0 (0) Comm time (%) = 4.65044e-06 (0.483395) Outpt time (%) = 6.52643e-05 (6.78397) Other time (%) = 0.000887418 (92.2436) Nlocal: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Total # of neighbors = 0 Ave neighs/atom = 0 Neighbor list builds = 0 Dangerous builds = 0 LIGGGHTS finished Foam::cfdemCloudIB::reAllocArrays() nr particles = 1 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.06984e-07, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 4.99282e-07, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 6.95219e-10, No Iterations 2 DICPCG: Solving for p, Initial residual = 1, Final residual = 9.26695e-07, No Iterations 573 time step continuity errors : sum local = 3.69904e-10, global = -1.28591e-13, cumulative = -1.28591e-13 DICPCG: Solving for p, Initial residual = 0.0476412, Final residual = 9.26405e-07, No Iterations 513 time step continuity errors : sum local = 3.16922e-08, global = -3.49546e-12, cumulative = -3.62406e-12 DICPCG: Solving for p, Initial residual = 0.0055667, Final residual = 9.51869e-07, No Iterations 500 time step continuity errors : sum local = 3.3706e-08, global = -6.10801e-12, cumulative = -9.73206e-12 DICPCG: Solving for p, Initial residual = 0.00178549, Final residual = 9.94251e-07, No Iterations 473 time step continuity errors : sum local = 3.53268e-08, global = -5.71629e-12, cumulative = -1.54484e-11 No finite volume options present particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 1, Final residual = 9.63924e-07, No Iterations 501 ExecutionTime = 30.7 s ClockTime = 32 s Time = 0.0002 Selected 784 cells for refinement out of 1613716. Refined from 1613716 to 1619204 cells. Selected 0 split points out of a possible 784. Courant Number mean: 0.126921 max: 9.46381 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' Setting up run at Fri Jul 26 12:59:56 2024 Memory usage per processor = 6.75634 Mbytes Step Atoms KinEng rke Volume dragtota dragtota dragtota 11 1 3.5918126e-08 0 0.0260508 0 0 0.017347056 CFD Coupling established at step 20 21 1 3.6972707e-08 0 0.0260508 0 0 0.017347056 Loop time of 0.00149142 on 16 procs for 10 steps with 1 atoms, finish time Fri Jul 26 12:59:56 2024 Pair time (%) = 3.75399e-06 (0.251706) Neigh time (%) = 0 (0) Comm time (%) = 3.65471e-06 (0.24505) Outpt time (%) = 8.7527e-05 (5.86871) Other time (%) = 0.00139648 (93.6345) Nlocal: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Total # of neighbors = 0 Ave neighs/atom = 0 Neighbor list builds = 0 Dangerous builds = 0 LIGGGHTS finished nr particles = 1 -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 15 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. --------------------------------------------------------------------------

谁能帮忙填充一下没考虑到的内容么？

捕获.PNG

Direct numerical simulation of reactor two-phase flows enabled by high performance computing

捕获.JPG

@bestucan
感谢您的回复！“不同方向的流动计算差异就大”这个解释很直接明了，我突然就有点明白这个网格线出现的原因了:xinxin:

请教各位大佬：
openfoam有没有自带的插值可以将volField转化成pointField，使用tecplot进行后处理时插值操作会将cell data转化成point data，操作比较麻烦于是自行编写openfoam程序，很多操作是基于volField，然后结果和tecplot有比较大的区别，不知道是不是插值的影响

@李东岳 请问这可能是是网格太小导致吗？

论坛里面有没有高手遇到overset 的后处理问题：
问题描述：fluent overset 方法计算翼型表面转拟，但是在显示翼型面上切应力时，数据分布呈现锯齿状分布。
下图1是overset cell type.
9de13be5-dde7-4e34-ba97-95264d78b5b8-image.png
图2是翼型表面的切应力
5b284d9b-e3bc-45a4-b92e-d749155a96c3-image.png

国内外有哪些流体力学经典著作？ - TurbPhys的回答 - 知乎
https://www.zhihu.com/question/335644582/answer/822877093