CFD中文网

这几天在看有限体积法的书，大佬的论文、《数值传热学》和 《The FVM inCFD OpenFOAM..Matlab》都看了，逛论坛发现都在推荐《Computational Methods for Fluid Dynamics》这本；去SpringerLink一搜正好今年出了新版，网上资源好像还不多，在这里分享一下： 点击下载

你可以读取openFOAM处理好的网格文件啊，它的网格文件还是很有规律的，至少我现在已经实现读取和处理openFOAM的网格数据了，自己在编写求解器。

https://codepen.io/pkra/pen/EPeKjo

@东岳 好的，李老师，我先按照《computational methods for fluid dynamics》这本书上说的那种把这一项也处理成梯度的形式试着去推导一下，看能否推导出来

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

列表

@Samuel-Tu 会的

@Samuel-Tu 可行蓝色可以设置为internal面，顺便统计一下风量

@wwzhao

简单的感受了一下两种赋值方式，效率确实很大提高。

#include <stdio.h> #include <string> #include <chrono> #include <iostream> using namespace std ; class Timer { private: // Type aliases to make accessing nested type easier using clock_t = std::chrono::high_resolution_clock; using second_t = std::chrono::duration<double, std::ratio<1> >; std::chrono::time_point<clock_t> m_beg; public: Timer() : m_beg(clock_t::now()) { } void reset() { m_beg = clock_t::now(); } double elapsed() const { return std::chrono::duration_cast<second_t>(clock_t::now() - m_beg).count(); } }; class SomethingBefore11 { private: int m_array[5]; public: SomethingBefore11() // zero the member array { // If we want the array to have values, we'll have to use assignment here m_array[0] = 1; m_array[1] = 2; m_array[2] = 3; m_array[3] = 4; m_array[4] = 5; } }; class SomethingAfter11 { private: int m_array[5]; public: SomethingAfter11(): m_array { 1, 2, 3, 4, 5 } //zero the member array { } }; int main(){ Timer tBefore11; SomethingBefore11 m_array_before11; std::cout << "Time elapsed: " << tBefore11.elapsed() << ‘n’; Timer tAfter11; SomethingAfter11 m_array_after11; std::cout << "Time elapsed: " << tAfter11.elapsed() << ‘n’; return 0; } [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.506e-06 Time elapsed: 1.47e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 8.664e-06 Time elapsed: 1.9e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.646e-06 Time elapsed: 1.89e-07

@bestucan 谢谢大神的解答，受教了

@东岳 :xiexie:

@东岳 感谢李老师！

@李东岳 谢谢李老师

按网格类型分类的软件列表：
Mesh Generation & Grid Generation on the Web

@李东岳 好吧，谢谢大佬！

0_1516841984913_捕获.PNG
动图1：河床沙石冲刷运动

河流的流动会改变河床的形状。底部的沙子以及大石头通常会分为两个区：

顶层的bed-load zone（有研究河海的么，这个怎么翻译？），其中的颗粒会随之移动； 底层：沙子以及石头会被困住，但是会缓慢的爬行；

在碎石较多的河床中，较大的沙石倾向于堆积在顶层，这种现象被称之为armoring。个人感觉之所以称之为armoring是因为较大的沙石趋向于在表面，类似一种武装。

实验表明，在这个区域，大沙石有一个净得向上的运动速度，小颗粒具有一个小的向下的移动速度。 目前对于这种现象的原因尚不清楚。不过现存了一些潜在的理论。例如支持程度比较高的巴西果效应（Brazil Nut Effect）。

0_1516841984913_捕获.PNG
动图2：有人晃荡大米里面的密度差不多的小石头

替代文字
动图3：进行的计算模拟

巴西果效应是指如果把两种颗粒的混合物置于容器中，然后施加外加的振荡，体积比较大的颗粒会上升到表层，而较小的颗粒会沉降到底部。关于这种古老的效应的动力学机制至今仍众说纷纭。

DEM模拟：文中作者采用LIGGGHTS对这个现象进行了模拟，下图可以看出，随着时间的推移，大颗粒漂浮出来。

0_1516841984913_捕获.PNG
图4：Nature Communication中的LIGGGHTS模拟

看起来很有意思 :cheeky:

本文来自Nature Communication

嗯 AP 很多AP都很强力 强力AP！

@东岳 谢谢大佬~

最近看到了一个新的并行分解算法，KaHIP，依然存在疑问。不清楚这种不同的挂载方式的数学本质是什么。

http://algo2.iti.kit.edu/kahip/

替代文字

我试试，不过感觉你的算法没毛病，M3守恒。
你考虑过用Wheeler求节点么？比PD要强健。