@浪迹天大
需要一个动网格求解器,然后添加动网格字典,就是动网格~
CFD中文网
帖子
-
-
@vivian
错误显示迭代步数超出了。我采用了debug模式在rhoCentralFoam.C的215行加入断点,想输出e,结果屏幕输出了大量信息,我应该如何输出才能得到想要的场信息呢?
e是个场,你要输出场的话,就会输出所有的场的值,不清楚你要看什么,也可以单独输出场的最大值最小值看看越界否。如果一定要查看整个场的信息,输出到log再看,这样log文件会变得很大。
rhoCentralFoam &> log
如果我想梳理出程序运行的主线,搞清楚从场的离散点到通过离散方法,获得最终求解的代数方程矩阵,我认为最好的办法是从建立场开始,关注场中存储的离散点的信息变化。那么请问如何输出存储点压力,面通量信息的矩阵呢?
从离散的矩阵系数如
fvScalarMatrix
组建我们常规的稀疏矩阵不是一俩行代码可以搞定(估计不到100行?)因为openfoam里面存储的是稀疏矩阵,并用指针+指针。所以,可能得用点时间。更简单的方法,是回到问题本身,很少是因为矩阵系统解不出来的原因,而是设置有问题导致发散。
-
So what is your question?
-
@MX
可以通过壁面函数nutkroughwallfunction来实现 -
对于一个面,要计算体积流量,另一种方法是通过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系统。你用的什么浏览器?没有这个标志?
-
maxCo取决于时间格式、流动状态、网格质量、求解模型。
显性格式maxCo可以最大为1;
隐形格式maxCo可以适当扩大如2;
规则几何和流型库朗数可以适当扩大;
多相流库朗数可以最大为1,levelset可以为0.5或者更低;
动网格要确保点移动求解器每个时间步不能过大引起负体积;等等所以是个很复杂的话题,并没有统一标准。可以逐步加大库朗数,知道发散为止。然后再降低。:sunglasses:
-
泊松方程如果给定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
,因此不需要满足附加条件。 -
我同意你的观点,但是les模拟和网格关联性很大。也可以在大尺度下做les模拟,比如,同一个工况,不同的网格尺度,执行les模拟,在这种情况下,过滤的涡是不同的:cheeky:
-
这个是不是更简单点
-
好像目前还是最后一点用的比较多。模拟都要和实验对比,才能说明你算得准。
模拟和模拟对比的文章我没怎么遇到。模拟终究逃脱不了实验啊!
-
Hi @simthere ,好久不见:sunglasses:
并没有太细看流沙的原文,针对这几个关键词,简单说一下我的理解:
- 相对压力是绝对压力-大气压力;
- 表压=相对压力;
上述俩个压力在工程类的书里面都应该有介绍。不同的是参考压力:
- 参考压力是CFD下的产物。和表压、大气压力无关。参考压力的出现是不可压缩流体中的特例。某些情况下,求解NS方程的时候,如果不给定压力的参考压力,压力方程无解。因此需要给定一个参考压力。求解出来的压力和参考压力有关;
-
延长了出口和进口,计算的结果有很大差异.是否说明没有延长进口和出口的结果不可行了?
我更倾向于延长出口和进口。因为这很好解释。
残差震荡正常吗?残差一直下降正常吗?
很正常,尤其是一个固有的瞬态现象用稳态模拟的时候。看这个第二问:链接文本
从这个图片来看,网格无关性已经做好了。即88万之后,压差达到稳定且不随网格变化而变化。
另一方面,从88万网格的流出系数来看,值为0.800. 实验值为0.72. 误差在10%。已经可以接受。
其他网格的数据,误差均为20%之内,也都可以接受。
可以换一种思路:本文基于45/68/88/112/140万网格进行了无关性分析,结果如下。在网格88万之后,压差达到稳定且不随网格变化而变化。流出系数和实验值误差在10%之内,提高网格数量并不能减少这个误差。因此在下面的模拟中采用了88万的网格。
-
-
网格数量小更接近于实验,有可能是巧合造成的结果。这种网格无关性分析放在文章中,可能并不会很好的解释。
不过确定的是,网格数量越多,迭代收敛越困难,这在教科书上一般都有提及,多重网格求解器因此而生。
另外,楼主12万的网格最接近真实解,误差0.162%。 7万的误差3.59%。 尝试过5万的网格么?
还有,你的监测点,是进口和流场中的某个点?如果不是进口,延长进口,通过充分发展的湍流可以获取更稳定的解:sunglasses:
-
真是一个大话题。非结构网格,即时最标准的几何,比如圆柱体。质量都不会很高。这是由于四面体本身特性决定的。看国外一些高手使用非常复杂的技巧处理非结构网格,依然质量一般般。提高网格质量,六面体网格莫属。不过生成很复杂。
-
-
金榜题名应该是最荣耀的时候。可是。事情远非这么简单。
如果你认为进了清华北大,随便一个专业都是可以的,那么你大错特错
-
一个好的专业的本科,完全可以秒杀清华、北大的某些专业的博士;
-
以我gf的经历来看,毕业于"北京某IT技校",真名就不说了,大家都知道。本科纯物理学,硕士改行计算机。2016年3月毕业,2015年10月份实验室大部分同学就签出工作。gf工资最低,年薪23万+京户。其他大部分为30万左右无京户;
-
以我同学的经历来看。广大同学对不住了我要开扒了。2014年6月毕业的同学,同样也是硕士,大部分也在混月薪的阶段。待遇要比gf同学低的多。并且经常各种出差;
-
以我同学本科的经历来看。大体参照硕士同学的情况。因为专业相同。需要提及的是,本科里面有一个奇才。我俩关系甚好。他本科是生物材料专业,但是对计算机却产生了浓厚的兴趣,本科专业课全部低空飘过。毕业后找了一个月薪4000的工作,当时公司(某业内知名网络安全公司)还只有10多个人。我这同学在这公司一干就是现在。目前公司已经200多人。我同学也做上了技术副总裁。并且迎娶了白富美(富不富不知道,但是确实又白又美,请联想网络上的车模脑补);
其实,如果关注就业的同学很明显的都知道。目前最火的莫过于计算机和金融。说实话,金融不了解。网传有些投行年薪三五十万。但是由于我的所有同学和投行均无接触。因此并不十分了解金融行业的就业。
单说计算机,不得不承认。IT的本科可以秒杀某些专业的博士。你需要深刻的理解这一点。
- 举个栗子,硕士舍友的同学东北某大学毕业,工程类本科。毕业去某鸟学了计算机。是的,没错,就是那个某鸟。然后来北京找了一个游戏公司的工作,后来1-2年就做到了主程序员。工资1万+。
有些同学不高兴了说:别的行业越老越吃香
当然了,我给自己找的理由就是这个。不过话说回来,别人起点就比你高了2,3倍,后来居上也不是容易的事情
有的同学又说了:IT民工命短
然而,化学行业不命短?之所以大家说IT民工命短的原因在于,比较累。确实,gf所在的IT行业加班现象比较多。但是在哪不都是坐一天?坐在电脑之前coding是码农最幸福的事情。
最后,大家一定听说过这句话:全民CS,是的,就是这个意思(CS就是计算机科学)。
作为一个非计算机专业的博士,我拥护我的专业。然而不得不承认的是,行业的重要,远高于你手里的一纸学位证书。
-
-
不匿名,中航工业624所搞航空发动机的,工作9年,本科生,明年评高工,普通职工,不是领导干部。
税后年收入十万上下,在成都市郊区,房子单价五千左右。住140㎡的房子,开十万左右的车。自己觉得还行,就是每隔一段时间工作压力略大。9年年收入10万年。
看来还是得CS -
比如爱荷华州立的Anson Marston Distinguished Professor
二维自适应网格
rhoCentralFoam求解器的问题
边界条件的问题
OpenFOAM 粗糙度
关于两相流计算某一相流过某一个面上的流量
controlDict文件中关于自动调节库朗数
adjustPhi的作用是检查边界条件?
大窝模拟中的大,多大算大?
全导和偏导,怎么分清楚?
如何验证fluent计算结果的合理性
讨论Fluent中的各种压力
越多的网格得到越接近实验的解?
用满足网格无关性的网格数计算得到的结果是否最符合实验结果?
用满足网格无关性的网格数计算得到的结果是否最符合实验结果?
带尖角的几何体怎怎么处理比较好
使用ICEM绘制非结构网格时,如何提高网格质量?
Aspect ratio 达到300多,如何调小呢
这个论坛是用什么做的?
行业的重要,远高于你手里的一纸学位证书
南方科技大学力学与航空航天系诚聘海内外优秀人才
David Croll Fellow Professor