CFD中文网

是的，我有的时候也会出现这个情况。点pre-Mesh，体网格丢失，不知道还有没有其他人遇到这个问题？

@东岳 这个图是刚刚改用稳态算，迭代步数不多，迭代次数多了偏差还是挺大的

@星星星星晴 可能我得太复杂了，一共200多页，用他那个就出错

我知道了，公式里面用不了，文本里面可以用

好吧，我查了一下，附加一个yhmath这个包 :136:

https://tex.stackexchange.com/questions/100574/really-wide-hat-symbol/101136#101136

Data-Create Zone-Rectangular提取出想截取的位置
d2469623-64c0-45c3-952d-a6339a3b560b-image.png
但是提取出来的Zone没有数据，Data-Interpolate根据原数据对创建的Zone进行插值就可以
abd5cf18-caaf-4a49-a528-3b9afd696f35-image.png
这样就得到一个名为Rectangular zone的Zone
758ce629-1262-40bf-83b2-1477d15a474d-image.png

请问李老师，这是垂直管道的流型图吗？如果要在论文引用这张图片的话，引用怎么写？就写CFD中文网吗？

嘿，我还没关注这个历史。还真是。太折磨人了。

没有吧，部分情况下 2D 3D udf 有差异

666eed00-e17c-4be5-8a8c-0ec0e66e1275-image.png 是用哪个UDF函数能输出转速，能给出个例子吗？

@李东岳 @lsprxd @赵一铭 对于标准壁面函数我是这样理解的：在求解流场的时候还是用原来的k方程和epsilon方程来求解，只不过是将在壁面第一层网格（有可能是网格面上的值，有可能是网格中心的值）上的信息用壁面函数求解的结果来修正，而且壁面函数其实只求解壁面第一层网格的信息，第二层网格之后就是输运方程求解的结果了。不知我这样理解是否正确？

老铁这个整的太高端了，我只知道两个点直接的距离，不知道3个点的距离是什么意义...

@东岳 您好，我刚接触fluent不久，请问您说的sweeps数是什么意思，想请问是提高网格质量和调节松弛率吗？

@李东岳
针对 ncompressible::turbulenceModel 模板特化后的 ShihQuadraticKE 模型 ，仔细研究了一下其怎么实现非线性项的实现，发现其是在 NonlinearEddyViscosity.H 中定义 nonlinearStress_ 项，然后用 ShihQuadraticKE 作为子类继承父类 NonlinearEddyViscosity 以便在 ShihQuadraticKE.C 中直接引用omega_。

但由于我修改的是SSTkomegabase.H( 继承于TurbulenceModel类 )与 linearViscousStrss.H ( 继承于BasicturbulenceModel原则类 ) 没有特例化，不能简单类似上面的ShihQuadraticKE 模型。需要从新仿造框架，是一项大工程！

我的策略是：在现有的文件上直接硬植入，简单很多，毕竟涉及到湍流模型的类错综复杂，针对找不到omega，按照我的分析应该是对的，运行求解器时linearViscousStrss.C先运行此时并未运行SSTkomegabase.C 所以没有omega场注册，找不到他，所以我在creatfield.H中IOobject一下omega场，目前可以算了，结果还未检查，不知道我的理解有没有问题。

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

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

替代文字

直接移动到cad原点就可以了吧，cad的原点很好找的呀

加上来了

不知道大家感受到没有。换了一个开源的搜索引擎，叫meilisearch。简直屌爆了。。。。。

:146: :146: :146:

百度 google也搜不到啥。我就写在这里给大家提个醒吧。
之前尝试从华夏银行APP汇入法国账户1万欧。
应该是账户写错了。
折腾几天系统自动退回来了。
扣了90欧。
华夏银行这面没有扣钱。
就是在欧洲不知道是哪个环节扣了90欧。
大家知道就好。以后要跨境汇款心里有个数。

微信图片_20240214134432.jpg

这个问题已经解决，麻烦大家了

@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

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