$E=mc^2$
youmengtian
帖子
-
-
@Foamer24 你可以贴一个fvSchemes完整的设置看一下,很有可能是离散格式的问题。
-
@Foamer24 冒昧回复仅供参考,我觉得还是在于fvSchemes中离散格式的选用。对于skew较大等网格质量不好的问题其实就两个思路可供选择,第一重画网格,第二如果的确无法提高网格质量只好退而求其次选择数值耗散更强但是有界性相对较好的的离散格式。其实tutorial中给出的很多示例很难拓展到工程背景下的计算中。对于计算效率的问题往往跟选择的线性矩阵求解方法有关,目前尚未有普适的定论。您给出的结论也并不能涵盖所有的计算情况,您给出的截图也是OpenFOAM更快一些,但是往往需要选择合理的求解方法。根据Fluent的官方手册,它广泛使用的是AMG的求解方法,特别详细的内容我也并不清楚。
以上,希望没有传递特别离谱的错误信息,如有不当请及时联系,共同学习。
-
@东岳 老哥具体说一说,让我们也学习一下
-
@东岳 正好这两天下载安装了一下,icem还在但是基本没变化。感觉ansys这几个版本还是在做整合,大的改变很难搞了。
-
有的呀,还有的
-
@random_ran 您说的太客气了,虽然咱们这儿叫CFD中文网但是也没说不让使用英文。你之前的阅读逻辑是没有什么错误的,但是我不知道你怎么读到polyModifyPoint.H这个头文件去,因为我觉得貌似跟这个头文件没有关系。但是我感觉你对错误的基本定位应该是准确的,应该就出现在polyModifyFace.H第190行开辟新的智能指针出错了
polyModifyFace.H 188 virtual autoPtr<topoAction> clone() const { return autoPtr<topoAction>(new polyModifyFace(*this)); }
-
@yuhx 你想加option的时候加,可以参考/src/OpenFOAM/include中相关内容
-
@Yongbo 看预计多大的计算量。CFD这东西感觉也没个上限吧,不过个人工作站还是有极限,我觉得貌似个人财力是上限,在这个限制之下买最后的。
-
@lllwonderliquid 方法可能不止一种,我给你提供一个思路吧,可能需要写一些程序,通过fvOptions实现。具体可以参考limitTemperature。
./system/fvOptions limitT { type limitTemperature; active yes; selectionMode all; min 200; max 500; phase gas; //optional }
代码实现在 src/fvOptions/corrections/limitTemperature/ 目录。
-
另,建议改小matchToTolerance。相对于你的网格尺寸,这个值明显偏大。
-
@random_ran polyModifyFace是一个类,而非函数。另外在polyModifyFace.H同级文件夹下有polyModifyPoint.H 。至于抛出的异常之前未见过,不过看样子应该是在分配内存时出现了异常,非常有可能是你说的问题。但是实际上可能和你演示的例子还是有差别,固定数组能够使用的栈本身大小有限很容易出现这种异常。但是实际中像createPatch之类需要对网格数据进行处理的往往是通过动态数组开辟内存空间的,其上线与计算机内存有直接关系。createPatch在运行中需要额外生成相关patch的数据信息,需要占用相应大小的内存,如果网格过大很有可能超出计算机内存上限抛出错误。
另,过于createPatch具体实现我也没有进行深入的分析,仅提供一个思路,正确与否请谨慎。
-
@lllwonderliquid 可能你对插值格式的理解有误??vanLeer01这里的0-1指的是相对值。你想给定k的范围可能是绝对值?那就需要使用bound进行限制。
-
@random_ran 从check的结果来看,原始生成的网格没有什么问题。您说的周期性边界失败是指?可是尝试使用createPatch这个utility创建cyclic边界实现周期性边界。
-
可以贴个详细的case来大家都跑一下,这样比较容易定位问题。这样问的话还是比较模糊。
-
@东岳 vim or emacs
-
这个我也挺好奇的。看OpenFOAM中的一些说明文档大约是用emacs org模式写的,看样子应该有很多开发者在用emacs。另外vim也是很好的一个选择,ide的话就很多了。我现在用eclipse结合git插件看代码貌似也不错,但是eclipse貌似有点儿太累赘,已经不能算是编辑器的范畴。我觉得,这个因人而异吧,用自己最习惯的就好。
-
gdb跑一下,这一个错误提示根本看不出来个啥
-
@李东岳 最近也在做打算。数据越来越多,做一个nas貌似挺有必要的
-
对,最近也是发现这个改动。不太明白出于什么目的,之前函数重载用着挺舒服的。不知道为什么要改成这样
-
@李东岳 有钱还是做一个NAS私有云,配上无线局域网,备份同步都很好
-
@李东岳 感觉并没好多少,至少对于耦合求解来说依然很难进行下去
-
@程迪 对的,这个模板类是个残疾。。。前几天给esi的开发组发了邮件问了一下,他们说比较“健全”的耦合求解矩阵在18-19年的版本中发布
-
@程迪 最近在cfd-online上看到你很多的帖子回复。我最近也在关注LduMatrix这个模板类,有时间可以交流一下。
-
@李东岳 可以看,但是稳定性还是不比付费vpn
-
@李东岳 我用xx-net
-
@流卜 的确是接口差别很大,目前准备把extend上的部分内容移到v1706,感觉崩溃了快
-
foam-extend可以参考OpenFOAM-1.7版本。各种接口和模型各版本基本没有区别。
-
@程迪 程博士说的对呀,痛苦呀
-
@李东岳 妥妥的,不过感觉超宽屏并没啥,不如上多屏
-
@李东岳 统一李老大的说法。这个是一个现实问题,如果没有商业驱动,这种专业级的开源软件真的很难推进。
-
@影川风 的确差别很大
-
@nuaawqf 今天测试了dbnsTurbFoam自带的naca0012算例,最大CFL数在100依然能保持稳定,最大自适应步长7.14e-5。更大的CFL目前暂时没有测试。感觉使用隐式时间推进之后求解器的鲁棒性提高了不少。。
-
最好贴一下你的错误信息,不然的话这点儿信息貌似也很不好判断。
-
@nuaawqf 是这样,很不稳定。特别是设计粘性后更是很复杂。我写的东西暂时也很初级。如果你已经接触dbns库了,那就有的聊了。我的思路是在numericflux类中通过成员函数实现对流通量雅各比矩阵,基本思路与沈淳老师的实现方式一致。之后可以给你部分代码。目前主要想怎么实现预处理形式的隐式时间推进格式',可以做深入探讨。我了解南航一直在cfd数值计算方面还是挺有自主开发能力的,向'你们学习
-
@sibo 其实大约都是基于虚拟化技术来的。我是给docker装的vncserver,虚拟一个主机,在的形式windows通过vnc访问。这样server和client
-
@李东岳 我感觉可能和bash的性能相似。但是bash配置完就只能在自己的电脑上使用,docker生成的硬盘镜像和虚拟机类似可以随意上传或者传播。并且通过dockerfile可以很方便创建一个openfoam的环境
-
@青瓷碎了夏天 面积分可以分两部分,第一步先使用面两边的体变量插值出面上的物理量,比如一个面k,它左右两边为i,j两个控制体,则P_k=interpolate(P_i,P_j)。这个值相当于面上的平均值,与面矢量内积就能够得到面积分了。
-
这个问题使用fluent算貌似同样会出现这个问题,其实这时候已经出现了部分空化,处理不当产生负压。
-
明显粘性系数不是通过这个成员函数得到的。这个函数是为动量方程提供粘性部分。粘性系数定义是mueff,通过湍流粘性系数与物性粘性系数加和得到的。
-
如题,发现win10下docker挺完备的,感觉上比虚拟机方便。有没有感兴趣的同学
-
@100yearsalone 恩,最近没过来逛逛咱们这个论坛。由于OpenFOAM官方版本一直以pressure-based为主导,所以会产生很多迷惑。
-
@100yearsalone 这一套理论是density-based部分的,openfoam基本都是pressure-based。不一样
-
去报个名,把最近的工作跟各位大牛交流一下。话说咱们的这个workshop论文是神马级别的?还是根本不算会议?
-
@nuaawqf 不知道你是用的OpenFOAM的什么版本,如果恰好使用的是foam-extend的话,实现vector5目前并不是非常复杂的问题,并且extend中包含了完整的块矩阵以及相关的数值解法。可以现在这里说一下求解思路:
1、基本思路可以参考东岳大神主持翻译的fluent理论手册的求解器理论部分,目前手头没有这个资料,应该在20章左右。其中有关于隐式时间推进法的相关内容;
2、其中比较复杂的部分就是构建对流通量Jacobian矩阵$A_c$,主要参考Blazek《Computational Fluid Dynamics Principles and Applications》附录中的表达式;
3、系统稀疏矩阵的构建可以参考fvm::UDiv的格式实现,并且基本思路一致。目前就这些实质内容,这两天测试CFL数能够在2000内保持稳定,当然只针对可压缩无粘流动,主要包括rhoCentralFoam中给出的几个经典算例,并且在32核的工作站上能保持相对良好的并行加速比。不知道你这边目前的侧重点在哪儿,可以进一步交流。
-
@daniel_lee 不认识,只是看他的文章。
-
@李东岳 是的,貌似openfoam上做超声速高超声速的很少。定位问题吧可能
-
看来大家做这方面工作的很少。
-
刚好最近在看thermophysicalmodels这一部分。感觉的确复杂,大家可以参考userguide中的定义部分。特别是1.x版本更接近C++的书写习惯,一步步查下去就能找到了。特别需要搞清楚的就是每一层到底定义了些什么
-
东岳大神说了开发的问题。我来说说我的看法?清楚物理模型的物理意义貌似是最基本的也是最容易出错的部分。openfoam方面可以参考The Finite Volume Method in Computational Fluid Dynamics。还有就是最笨的办法:读代码
请问代码块内部如何插入公式的?
关于rhoCentralFoam计算高超声速问题时温度出现负值的原因
讨论 Fluent计算效率优于OpenFOAM
论坛有人搞矩方法么?
ansys 2019版本icem模块没了?
ansys 2019版本icem模块没了?
OpenFOAM周期边界生成不成功,cyclic的问题
argList::addOption
Openfoam计算要选什么配置电脑比较好?
openfoam的fvScheme中是否可以给某个变量指定范围?
OpenFOAM周期边界生成不成功,cyclic的问题
OpenFOAM周期边界生成不成功,cyclic的问题
openfoam的fvScheme中是否可以给某个变量指定范围?
OpenFOAM周期边界生成不成功,cyclic的问题
关于rhoCentralFoam计算高超声速问题时温度出现负值的原因
成立党支部
openfoam developer
自定义Solver出现“浮点数例外,核心已转储”
vmware虚拟机打不开了,数据怎么弄出来?
internalField() 和 Ref()
vmware虚拟机打不开了,数据怎么弄出来?
最快速的上google的方式是什么?
最快速的上google的方式是什么?
最快速的上google的方式是什么?
最快速的上google的方式是什么?
最快速的上google的方式是什么?
OpenFOAM与extend命令区别有前辈总结吗
OpenFOAM与extend命令区别有前辈总结吗
一个不得不承认的事实
搞CFD的都能玩点啥?
一个不得不承认的事实
单核与并行时输出的区别
openfoam implicit density-based solver
buoyantPimpleFoam的边界条件使用问题
openfoam implicit density-based solver
有没有打算基于docker在windows下开发openfoam的
有没有打算基于docker在windows下开发openfoam的
OpenFOAM 中function object 对一个面上的P×V的值积分如何设置啊
compressiblemultiphaseinterfoam计算的bug
多相流求解器的Ueqn中 turbulence->divDevRhoReff 在哪段代码求了平均粘性系数?
有没有打算基于docker在windows下开发openfoam的
可压缩求解器中的音速的定义和离散格式
可压缩求解器中的音速的定义和离散格式
13th OpenFOAM Workshop 2018
openfoam implicit density-based solver
openfoam implicit density-based solver
openfoam implicit density-based solver
openfoam implicit density-based solver
怎么才能在代码里快速找到函数的定义?
OpenFOAM programming学习问题