CFD中文网

@李东岳 哦哦，没问题了。:142:

只知道是按照从小到大进行排序，这个sort函数有什么作用？

void Foam::SortableList<T>::sort() { sortIndices(indices_); List<T> lst(this->size()); forAll(indices_, i) { lst[i] = this->operator[](indices_[i]); } List<T>::transfer(lst); }

谢谢，知道了。

在气液欧拉多相流的大涡模拟中，有没有什么模型是把liquid或者particle相的SGS viscosity和主相的SGS viscosity联系起来的？
就像RANS中的一种模型：
$\begin{equation}
\frac{\nu_{t,particle}}{\nu_{t,gas}}=(\frac{k_{particle}}{k_{gas}})^2=\frac{1}{1+\frac{t_p}{t_{fl}}}
\end{equation}$
其中$\nu_t$是湍流粘度，k是湍动能，$t_p$是particle relaxation time，$t_{fl}$是Lagarangian fluid time scale。
现在像找一个类似的模型用在LES上，但是翻了好久也没翻到 :135:

感谢各位大佬！:chitang:

看到个帖子不错，记录一下。关键词：householder transform

我也是相同的问题啊！！！

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>

目前想在这个求解器中添加煤燃烧的反应，将利用反应速率计算需要添加的各个方程的源项。
1 目前的问题的我在src中的thermophysicalModels找到关于固相化学反应的东西，但是找不到有关化学反应速率的求解公式，例如固相阿雷尼乌斯求解反应速率，所以感觉没有地方能下手修改成自己需要的反应速率？
2 如果要直接在求解器中编写一个函数，例如是计算反应速率，我我需要调用到的其他数据（施密特数，直径，雷诺数）那些的需要用到什么代码（原谅我是一个C++沫沫）

谢谢东岳老师，我先看下你网址里面的相关内容。

你应该看下速度分布云图，没有粘性并不代表不产生旋涡。另外上交有一个圆柱绕流的标准算例，可以在百度文库中找到。不要用Fluent去算势流的东西，势流是格林函数法或者边界元那一套，从离散开始跟fluent的fvm就不是一种东西。这种做法从概念上就有问题，因为laplace方程本人没见过fvm离散的。

Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: sixDoFRigidBodyMotion Selecting sixDoFSolver symplectic Translational constraint tensor (1 -0 -0 -0 1 -0 -0 -0 -624) Rotational constraint tensor (0 0 0 0 0 0 0 0 0) PIMPLE: No convergence criteria found PIMPLE: No corrector convergence criteria found Calculations will do 2 corrections PIMPLE: Operating solver in transient mode with 2 outer correctors Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type laminar Selecting laminar stress model Stokes No MRF models present No finite volume options present Constructing face velocity Uf Courant Number mean: 0.000625613 max: 0.10552 fieldAverage fieldAverage1: Starting averaging at time 0 forces forces1: Not including porosity effects forces forceCoeffs1: Not including porosity effects forceCoeffs forceCoeffs1: Not including porosity effects Starting time loop Courant Number mean: 0.000625613 max: 0.10552 deltaT = 0.005 Time = 0.005 PIMPLE: Iteration 1 forces forces: Not including porosity effects #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::sixDoFSolvers::symplectic::solve(bool, Foam::Vector<double> const&, Foam::Vector<double> const&, double, double) at ??:? #4 Foam::sixDoFRigidBodyMotion::update(bool, Foam::Vector<double> const&, Foam::Vector<double> const&, double, double) at ??:? #5 Foam::sixDoFRigidBodyMotionSolver::solve() at ??:? #6 Foam::motionSolver::newPoints() at ??:? #7 Foam::dynamicMotionSolverFvMesh::update() at ??:? #8 ? in "/home/z/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/pimpleFoam" #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #10 ? in "/home/z/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/pimpleFoam" 浮点数例外 (核心已转储)

是单位化的，写起来忒麻烦了没写一会我改一下。

我试了一下
\begin{equation}
P=
\left(
\begin{matrix}
1&2\\
2&3
\end{matrix}
\right)
\end{equation}
这个就可以:136:

之前也困惑问过类似的问题，李老师说复杂的情况用ICEM 之后网格转换一下，这样确实是省时又省力哦。

@Alex 太厉害了，大神，我想学习icem，您还有别的类似的icem实例教程吗，不胜感激

列表

FLUENT的lience分为高校版与企业版，目前还有高校的优惠，去年应该也是有的，以这位老师的价格能够买到ANSYS CFD（包含FLUENT CFX ICEM等模块）高校版本。

这个潘工的联系方式给我呀，500块买不了上当啊

今天空闲时间又看了一遍笔记发现有些笔误：
P9-第三章第一节第二段第五行“将手竖直的防止在水中” 应为“将手竖直的放置在水中”

@东岳 老师，这个公式（1）是他们参考另一篇文章
Wu, K., Chen, Z., Li, J., Li, X., Xu, J., & Dong, X. (2017). Wettability effect on nanoconfined water flow. Proceedings of the National Academy of Sciences, 114(13), 3358-3363. https://www.pnas.org/content/114/13/3358.short
我的理解是因为在纳米管内靠近壁面的位置，粘度变化没有统一的解析解，所以看到一些文章就这么分成“多层”粘度，算是粗略简化的方法吧 :xiezuoye: