@浪迹天大
需要一个动网格求解器,然后添加动网格字典,就是动网格~
CFD中文网 发布的最佳帖子
-
RE: rhoCentralFoam求解器的问题
@vivian
错误显示迭代步数超出了。我采用了debug模式在rhoCentralFoam.C的215行加入断点,想输出e,结果屏幕输出了大量信息,我应该如何输出才能得到想要的场信息呢?
e是个场,你要输出场的话,就会输出所有的场的值,不清楚你要看什么,也可以单独输出场的最大值最小值看看越界否。如果一定要查看整个场的信息,输出到log再看,这样log文件会变得很大。
rhoCentralFoam &> log
如果我想梳理出程序运行的主线,搞清楚从场的离散点到通过离散方法,获得最终求解的代数方程矩阵,我认为最好的办法是从建立场开始,关注场中存储的离散点的信息变化。那么请问如何输出存储点压力,面通量信息的矩阵呢?
从离散的矩阵系数如
fvScalarMatrix
组建我们常规的稀疏矩阵不是一俩行代码可以搞定(估计不到100行?)因为openfoam里面存储的是稀疏矩阵,并用指针+指针。所以,可能得用点时间。更简单的方法,是回到问题本身,很少是因为矩阵系统解不出来的原因,而是设置有问题导致发散。
-
RE: 关于两相流计算某一相流过某一个面上的流量
对于一个面,要计算体积流量,另一种方法是通过phi乘以面积来获得。
对于面速度矢量,压力和每相流量是怎么给出的?
面速度=
fvc::interpolate(U)
是插值过来的。不需要后处理,再求解器里面已经处理好了。所以显然的phi=fvc::interpolate(U) & S
,Q=fvc::interpolate(U) & S * mag(S)
,phia=alpha1f * fvc::interpolate(U) & S
,Qa=phia * mag(S)
多谢反馈。手边目前没有linux系统。你用的什么浏览器?没有这个标志?
-
RE: controlDict文件中关于自动调节库朗数
maxCo取决于时间格式、流动状态、网格质量、求解模型。
显性格式maxCo可以最大为1;
隐形格式maxCo可以适当扩大如2;
规则几何和流型库朗数可以适当扩大;
多相流库朗数可以最大为1,levelset可以为0.5或者更低;
动网格要确保点移动求解器每个时间步不能过大引起负体积;等等所以是个很复杂的话题,并没有统一标准。可以逐步加大库朗数,知道发散为止。然后再降低。:sunglasses:
-
RE: adjustPhi的作用是检查边界条件?
泊松方程如果给定Neumann边界条件,需要给定一个附加条件。在NS方程中,这个附加条件就是就是要保证速度边界条件是守恒的。
例如:我们有1D压力方程:
\begin{equation}
\frac{\partial }{\partial x} \frac{\partial p}{\partial x}=k
\end{equation}
其解为:
\begin{equation}
p=\frac{1}{2}kx^2+C_1 x+C_2
\end{equation}
如果给定Neumann条件,$p'(0)=\alpha$,$p'(1)=\beta$
带入后我们有:
\begin{equation}
C_1=\alpha, C_1=\beta-k
\end{equation}
这个就是附加条件,在此情况下,$C_2$无解,因此在压力均为Neumann边界条件的时候,需要设定参考压力,并且满足公式3. 公式3在NS方程中即为adjustphi()
的目的。从代码看这个函数只修正了边界上的流量,没有碰内部场,这样有什么意义吗?就是为了强制边界上总的通量是守恒的
修正边界条件确保压力泊松方程有解。
它不检查存在给定压力边界(p.needReference()为false)的算例,是默认这种算例能最终达到边界通量的守恒?
如果存在Dirichlet边界条件,则
p.needReference()
为false
,因此不需要满足附加条件。 -
RE: 大窝模拟中的大,多大算大?
我同意你的观点,但是les模拟和网格关联性很大。也可以在大尺度下做les模拟,比如,同一个工况,不同的网格尺度,执行les模拟,在这种情况下,过滤的涡是不同的:cheeky:
-
RE: 如何验证fluent计算结果的合理性
好像目前还是最后一点用的比较多。模拟都要和实验对比,才能说明你算得准。
模拟和模拟对比的文章我没怎么遇到。模拟终究逃脱不了实验啊!
CFD中文网 发布的最新帖子
-
公众号影响力分析
2018年第一季度CFD/CAE类自媒体影响力分析
本统计只收录2018年第一季度的消息,统计日期为2018年3月30日。
相关数据有证可查(来源于各个公众号),本帖只做事实统计。
- 原创号:内容大部分为原创;很少转载;
- 营销号:内容大部分为转载;偶尔原创;
%(#ff0000)[平均每篇浏览量]
本指数统计公众号发布的头条平均浏览量。计算方法为:总浏览量/篇数
公众号名称 篇浏览量 原创比例 转载比例 硬广比例 1 CFD界(原创号) 4243 %(#ff0000)[100%] 0% 0% 2 CFD之道(原创号) 1612 %(#ff0000)[100%] 0% 0% 3 ANSYS学习与应用(营销号) 1433 14% 78% 8% 4 shipCFD(原创号) 1152 %(#ff0000)[100%] 0% 0% 5 工程事(营销号) 884 36% 64% 0% 6 南流坊(原创号) 664 55% 18% 27% 7 风流知音(原创号) 570 %(#ff0000)[100%] 0% 0% 8 AutoAero(原创号) 424 %(#ff0000)[100%] 0% 0% %(#ff0000)[总浏览量前三名]
本指数统计公众号发布的头条的总浏览量。计算方法为:总浏览量
公众号名称 总浏览量 头条数量 公众号发布的消息数量(非转载) 备注 1 ANSYS学习与应用(营销号) 118947 83 2 CFD之道(原创号) 45144 28 51 3 CFD界(原创号) 38912 9 43 第四名为工程事(浏览量为20343),第五名为风流知音(浏览量为10266),后略。
CFD界头条备案
标题 截止至20180330日浏览量 评论数 备注 CFD名人堂:Dimitri Gidaspow 1446 7 原创 LBM-格子玻尔兹曼:Are you OKAY 3194 18 原创 碾压你智商:液体弹射蹦蹦床 1962 5 原创 CFD预测的雾霾北京醇 3319 6 原创 一起在CFD中吃鸡:CFD虚拟现实 9878 7 原创 开年创奇迹:蚂蚁才是流体力学大师 4030 9 原创 地表最强CFD课题组汇总!考研党福利 4542 21 原创 CFD深度学习,无法阻挡 7523 8 原创 不为人知的OpenFOAM 2298 6 原创 884.4782609118947 CFD之道头条概览
ANSYS学习与应用头条概览
南流坊头条
工程事头条
-
RE: of中的多相流求解器
$\alpha_1=1-\alpha_2$是用来保证alpha有界的,因为:
\begin{equation}
0<\alpha<1
\end{equation}
守恒体现在连续性方程中:
\begin{equation}
\frac{\partial \rho\alpha}{\partial t}+\nabla\cdot(\rho\mathbf{U}\alpha)=0
\end{equation}of中多相流求解器中的alpha方程是套了循环
具体指的是什么?MULES么?
-
RE: 请教几个关于壁面函数的问题,谢谢
@lsprxd 在 请教几个关于壁面函数的问题,谢谢 中说:
近 壁 区 高 Re数 采 用 的紊 动 各 向同 性 假 设 不 再 适 用,因而k-e紊 流 模式在 此 不 能 应 用,故 应 选 取 距 固 壁某 一 距 离 的 边 界作 为 计 算 边界。
文献出自哪里?
-
RE: 请教fluent后处理
还有在这种计算条件下应该如何考虑湍流模型呢?考虑到壁面附近主要为低速高粘的液体,可以认为是层流?而主要空间是较为高速且低粘的空气,会形成涡流影响破碎后的液滴运动轨迹,又还是应该选择湍流模型比如 SST, SAS 这种过渡的模型?求解释。
很难说,很难说具体用什么湍流模型,CFD中很多模型的选用都是很难提前知道,因此好多都是尝试各种模型进行对比。ke和kOmega对比,SST和SA对比,PIV和CFD对比,你应该庆幸幸好没有公认准确的CFD模型,CFD工作者才能有点事干。所以,你可以看看文献用的是什么模型,如果你模拟的数据要和文献对比,那么采用文献的模型,如果文献纯时间,碰巧你研究两相流的湍流模型,对比一下没准能发个Nature :cheeky:
-
RE: 如何获得cell信息?
OpenFOAM和我们的想法不太一样,回想有限体积法一般面上的值和周围的网格有关,因此OpenFOAM主要定义面的owner和neighber就足够进行任意迎风的插值格式了。因此你不会看到cell文件,你会看到owner和neighbour,如果想知道cell的信息,你需要知道包围网格的面,然后通过这些面确定cell的点,