@东岳 发现知乎那篇文章里提到,早期的物质点法计算速度也是很慢的,需要上集群才行。
.J.
帖子
-
-
昨天看了清华姚班的一篇知乎文章,“99行代码的《冰雪奇缘》”,介绍他在MIT读博期间做的CG相关的成果。为了计算加速,手写了编译器,实在厉害。
原文:https://zhuanlan.zhihu.com/p/97700605?utm_source=qq&utm_medium=social&utm_oi=44817232953344 -
@东岳 李博是不是在植入这个算法,我等你
-
@东岳 感觉这个问题还挺麻烦的,出现奇异矩阵后算一会就发散了
-
最近用MFiX试了下DQMOM模型,发现在求解积分点和权重输运的时候,涉及到的矩阵求逆操作会出现奇异矩阵的情况,然后程序就崩溃了。Rong Fan 2004, Application of the direct quadrature method of moments to polydisperse gas–solid fluidized beds这篇文章也确实提到了会有奇异矩阵的情况,并且提到一些解决方法。
我想知道的是,这个问题是算法本身就有的,还是代码实现的不够好?有人试过fluent里的DQMOM有没有这个情况吗?或者OpenQBMM呢?如果是软件问题可能考虑换软件了。
-
@东岳 什么时候可以像李博一样优秀
-
@东岳 非常感谢李博耐心回答。也就是说QMOM计算量比DQMOM小对吧?因为现在对MFiX平台熟悉点,所以想的是在这个平台做点工作,而恰好MFiX有DQMOM的植入。这两个算法的速度差很多吗?如果DQMOM计算量大的话,应该会有别的方面的优势吧?
我们这边是基于EMMS做一些介尺度曳力方面的工作,但是感觉有点做到头了,所以想换点思路,博士搞原创还是挺难。。
-
最近在看了下Rong Fan (2004) Application of the direct quadrature method of moments to polydisperse gas–solid fluidized beds这篇介绍DQMOM植入MFiX中的论文,这个算法和QMOM,还有李博你说的QBMM有什么应用场景上的区别吗?只是大概知道DQMOM和QMOM在求解数密度矩的策略上有不同。因为现在有个想法是利用PBM加入合适的聚并破碎模型,或者粒径增长项去描述颗粒聚团的输运,从而修正曳力,类似的工作目前看到过程所的Shanwei Hu有研究。他是用的QMOM,我在想用DQMOM是不是也行?
-
@玉玺 我觉得用mfix比较方便
-
@玉玺 可以用双流体也可以用ddpm,但是fluent这两个模型都不是很完善,现在一般都不用这个做流化床了。
-
@东岳 好~
-
@东岳 之前看你微信发过一篇交通堵塞的论文,说是可以类比激波。论文题目是啥?谢谢~
-
@fluidized_bed 就按照质量分数去定义四个组分的量就可以。
-
@fluidized_bed 现成算例没有,都是自己摸索。
-
是对的,然后定义一个mixture对应到biomass那一相就可以了。另外,你算流化床的话建议用mfix。
-
@light 我做流化床,但是这个卸料过程我没接触过。但是总体来说,双流体想算准一些过程比DEM费劲很多,因为除了曳力的影响,还有本构模型的影响,以及介尺度效应的影响。
-
@light 我看了原文,原文指的应该是出口边界条件不用设置了,但是你在fluent里画几何还是需要的。
-
@东岳 还没呢,明年申请
-
我也在盼着~
-
@东岳 明年三月才申请,现在方向还没确定好。不过倾向于找欧洲那边的,感觉目前中美形式不适合去美国~
-
@东岳 哈哈哈~自从李博大规模退QQ群以后,江湖再无群秀~
-
@东岳 太秀了~~
-
@东岳 这不明年3月份申请csc吗,现在想抓紧整点文章出来,到时候有点竞争力。感觉这时间一天天太快,每天龟速进展。。
-
@东岳 那还真的是挺纠结的~人生啊
-
@东岳 你过去应该是做博后这类的吧,高校允许留职然后读博后这样吗?如果可以感觉也挺爽啊,出去赚个外快相当于
-
工作得辞了?
-
@东岳 最近也在看这个,就发现和玻尔兹曼方程很像,因为之前看过颗粒动理学和经典分子动力学相关理论,都是基于玻尔兹曼方程推导的。看起来相空间的思想把这几者都贯穿起来了。
-
@bestucan 学习了,我们做颗粒动理学,只知道老师这么叫我们也跟着这么叫,原来还有这样一段故事。
-
@东岳 多谢李博指点
-
@东岳 世界第一~~强悍
-
李博好像博后的外导和他有往来?@东岳
也做气固,而且看了论文感觉深度及高,发的也多是JFM这类论文,但是却很低调的样子。有点想套个磁,明年CSC想确提前定下目标~~ -
@书生kao 用mfix的多孔介质模型试过类似的效果,不知道fluent多孔介质加双流体行不行
-
@嵇何劫康 你好,我看了下你说的几篇文献,应该都是用的DEM方法研究的吧,你们这个领域液固曳力是已经非常准确了吗?研究脉动对颗粒运动影响的过程中是否需要排除曳力不准确的情况,去分离考虑呢?例如用浸润边界法(IBM),会不会更好一些?如想法不成熟,望见谅!
-
@东岳 这个用的浸润边界法
-
试了下IBM的算例,后处理动画显示有点问题,先更新流场动画,结束后才更新颗粒的动画,请问怎么解决?
-
可以加mfix群:147181147
-
@东岳 艾特我居然没有消息通知,还是无意间看到这个帖子的
-
mfix可以做颗粒传热,但是颗粒相边界只能用顶温,不能用热流,气固传热也是可以实现的,所以大体上感觉可以做
-
@陈宏洋 至少得有目录页才能制作吧
-
@东岳 李博有情怀,用真金白银为作者打call。
-
@东岳 你曾经分享的版本是有书签的,我后来在springer官网下载发现没有书签,然后就想自己做一下。
-
很久没水论坛了。。最近摸索了一下给PDF书籍自动添加目录标签的方法,感觉挺实用。
https://www.cnblogs.com/Jay-CFD/p/9099812.html -
@李东岳 又试了下,用命令安装成功了,源码包直接编译有问题。
不过又有个新问题paraFoam
这个命令运行的时候说无法连接到x serverCreated temporary 'cavity.OpenFOAM' paraview: cannot connect to X server
x server是啥东西?
-
@李东岳 但是blockMesh命令找不到
ywj123450@DESKTOP-7HCQCPV:/mnt/f/BaiduYunDownload$ blockMesh blockMesh: command not found
我看编译信息最后有一段
touch: setting times of ‘/mnt/f/OpenFOAM-5.x/platforms/linux64GccDPInt32OptSYSTEMOPENMPI/src/Pstream/mpi/using:openmpi-system’: No such file or directory
是不是代表缺少文件啊?
-
@李东岳 就是win10下的bash
-
所有依赖库都安装完毕
在~/.bashrc末尾添加了
错误信息为:
ywj123450@DESKTOP-7HCQCPV:/mnt/f/OpenFOAM-5.x$ source ~/.bashrc ywj123450@DESKTOP-7HCQCPV:/mnt/f/OpenFOAM-5.x$ ./Allwmake Allwmake /mnt/f/OpenFOAM-5.x make: Nothing to be done for `all'. ======================================== Start ThirdParty Allwmake ======================================== ======================================== Build MPI libraries if required ======================================== Build Scotch decomposition library scotch_6.0.3 /mnt/f/ThirdParty-5.x/platforms/linux64GccDPInt32/scotch_6.0.3 scotch header in /mnt/f/ThirdParty-5.x/platforms/linux64GccDPInt32/scotch_6.0.3/include scotch libs in /mnt/f/ThirdParty-5.x/platforms/linux64GccDPInt32/lib ======================================== Build PTScotch decomposition library scotch_6.0.3 (uses MPI) /mnt/f/ThirdParty-5.x/platforms/linux64GccDPInt32/scotch_6.0.3 ptscotch header in /mnt/f/ThirdParty-5.x/platforms/linux64GccDPInt32/scotch_6.0.3/include/openmpi-system ptscotch libs in /mnt/f/ThirdParty-5.x/platforms/linux64GccDPInt32/lib/openmpi-system ======================================== Build Metis decomposition optional component Metis was not found ======================================== Build CGAL optional component was not found ======================================== Done ThirdParty Allwmake ======================================== Allwmake src no git description found Allwmake src/Pstream wmake dummy wclean mpi wmake mpi ....省略中间部分(这部分应该没问题)... UPstream.C: In static member function ‘static bool Foam::UPstream::finishedRequest(Foam::label)’: /usr/lib/openmpi/include/mpi.h:632:43: warning: use of old-style cast [-Wold-style-cast] #define MPI_STATUS_IGNORE ((MPI_Status *) 0) ^ UPstream.C:787:9: note: in expansion of macro ‘MPI_STATUS_IGNORE’ MPI_STATUS_IGNORE ^ g++ -std=c++11 -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -O3 -DNoRepository -ftemplate-depth-100 -DOMPI_SKIP_MPICXX -I/usr/lib/openmpi/include -I/usr/lib/openmpi/include/openmpi -pthread -IlnInclude -I. -I/mnt/f/OpenFOAM-5.x/src/OpenFOAM/lnInclude -I/mnt/f/OpenFOAM-5.x/src/OSspecific/POSIX/lnInclude -fPIC -shared -Xlinker --add-needed -Xlinker --no-as-needed /mnt/f/OpenFOAM-5.x/platforms/linux64GccDPInt32OptSYSTEMOPENMPI/src/Pstream/mpi/UOPwrite.o /mnt/f/OpenFOAM-5.x/platforms/linux64GccDPInt32OptSYSTEMOPENMPI/src/Pstream/mpi/UIPread.o /mnt/f/OpenFOAM-5.x/platforms/linux64GccDPInt32OptSYSTEMOPENMPI/src/Pstream/mpi/UPstream.o /mnt/f/OpenFOAM-5.x/platforms/linux64GccDPInt32OptSYSTEMOPENMPI/src/Pstream/mpi/PstreamGlobals.o -L/mnt/f/OpenFOAM-5.x/platforms/linux64GccDPInt32OptSYSTEMOPENMPI/lib \ -pthread -L/usr//lib -L/usr/lib/openmpi/lib -lmpi -ldl -lhwloc -o /mnt/f/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/openmpi-system/libPstream.so touch: setting times of ‘/mnt/f/OpenFOAM-5.x/platforms/linux64GccDPInt32OptSYSTEMOPENMPI/src/Pstream/mpi/using:openmpi-system’: No such file or directory
求助大家!感谢!
-
Weller是一个狂热的开源软件抵制者
这句话是不是有笔误
-
@李东岳 赞
计算机图形学中的流体模拟为什么计算那么快,可否借鉴其优势。
计算机图形学中的流体模拟为什么计算那么快,可否借鉴其优势。
DQMOM模型奇异矩阵的问题
DQMOM模型奇异矩阵的问题
DQMOM模型奇异矩阵的问题
多相流与矩方法的数学模型
多相流与矩方法的数学模型
多相流与矩方法的数学模型
生物质流化床气化
生物质流化床气化
李博之前发过一篇关于交通堵塞类比激波的论文叫什么?
李博之前发过一篇关于交通堵塞类比激波的论文叫什么?
生物质流化床气化
生物质流化床气化
生物质流化床气化
ICEM敞开边界的网格如何绘制
ICEM敞开边界的网格如何绘制
CFD神书开始预定咯!!
CFD神书开始预定咯!!
FLUENT欧拉-欧拉模型模拟双组份鼓泡流化床中的一些问题
有对Rodney O. Fox比较了解的吗?
胖猫要去美国了
胖猫要去美国了
胖猫要去美国了
胖猫要去美国了
胖猫要去美国了
胖猫要去美国了
论坛有人搞矩方法么?
论坛有人搞矩方法么?
《无痛苦N-S方程笔记》勘误
有对Rodney O. Fox比较了解的吗?
有对Rodney O. Fox比较了解的吗?
有对Rodney O. Fox比较了解的吗?
多层流化床的数值模拟如何实现?
拉格朗日中的湍流分散力模型
CFDEM结果用paraview后处理,流场与颗粒显示不同步
CFDEM结果用paraview后处理,流场与颗粒显示不同步
MFIX传热问题
MFIX传热问题
MFIX传热问题
批量自动化生成PDF目录标签 了解一下?
批量自动化生成PDF目录标签 了解一下?
批量自动化生成PDF目录标签 了解一下?
批量自动化生成PDF目录标签 了解一下?
在WSL里用源码包安装openfoam失败
在WSL里用源码包安装openfoam失败
在WSL里用源码包安装openfoam失败
在WSL里用源码包安装openfoam失败
你要的OpenFOAM编程杂七杂八
陶哲轩这篇论文大概讲了什么内容,有没有大神能讲解一下?