CFD中文网

对简单的长方体进行网格划分，再检查网格质量时发现在长方体的周围有很多有问题的face产生，也不知原因，很疑惑，求大神指点，万分感谢！

0T5OC4QVS4OWJO998JSUAU7.png

@江山不如凉竹 我就补充下求解器中保存组分的化学反应源项吧。
不知道你用的什么chemistryModel。假设你用的是StandardChemistryModel，那就将StandardChemistryModel.C文件中下面代码里的NO_WRITE改成AUTO_WRITE，编译完成后算例再运行一步保存就可以了。

// Create the fields for the chemistry sources forAll(RR_, fieldi) { RR_.set ( fieldi, new volScalarField::Internal ( IOobject ( "RR." + Y_[fieldi].name(), this->mesh().time().timeName(), this->mesh(), IOobject::NO_READ, //IOobject::NO_WRITE IOobject::AUTO_WRITE ), this->mesh(), dimensionedScalar(dimMass/dimVolume/dimTime, Zero) ) ); }

注：要拷贝到自己的工程路径下修改编译，不建议直接在原始代码里改。

想设置一个壁面接触角随时间增加减小，用公式描述
QMF9EZM9NI7B{$AS17ASF{M.jpg
但是模拟过程中不到0.01s就已经完全铺开了
8d0dfcc4-1ada-4bec-8a31-d4000f59d818-0~FKZ}KBX0WEF3X7@2%@NTD.png
接触角为定值的时候没有问题，请问接触角随时间增加减小的设置哪里出了问题，该如何设置

3965ae40-ebda-4f61-8653-2c5c5e4bbc2e-0e3073cf46e92704523c1c59afe2762.png

缺这个包：binutils

或者路径没配置对，可以用echo $PATH看看有没有/usr/bin这一项

Unknown rhoReactionThermo type Valid rhoReactionThermo types :

@李东岳 谢谢！

@香柏树 你这个问题解决了吗

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

各位老师，近期在研究立管涡激振动问题，想采用切片法，主要思路是针对具有大长细比特征的立管，我想沿展长方向设置多个切片，通过求解每个切片上二维圆柱的升力和阻力来获得立管不同高度处的外力，然后再将这些外力作为输入项，自己编一段有限元程序来求解立管结构动力学问题。
现在存在的主要问题和疑惑是：我如果想设置100个二维切片，应该采用什么样的方式实现循环求解这100个二维圆柱涡激振动问题，以及如何将这100个切片上求解得到的外力与有限元程序耦合实现流固耦合求解呢？
想请问各位老师有什么思路吗？

ICEM画二维网格导出时报错：
WARNING: Mesh has uncovered edges. ANSYS Fluent needs a complete boundary (lines in 2D) or it will give a variety of errors and not read in the mesh! If this was 2D Hexa, perhaps your edges are not associated with perimeter curves

模型描述：二维管道流体的测试网格，使用ICEM导入模型后：
1、划分块
2、在preMesh 里 使用Edge Paras进行边界层的绘制，
3、preMesh
4、在preMesh里Convert to Unstruct Mesh
导出网格选二维，其他未改变，出现上述报错
新手从B站学了一点划分网格操作就尝试了这一个，不知道哪里出错了，还希望前辈们解答

@李东岳 好的~谢谢东岳老师:xinxin:

论坛里面有没有高手遇到overset 的后处理问题：
问题描述：fluent overset 方法计算翼型表面转拟，但是在显示翼型面上切应力时，数据分布呈现锯齿状分布。
下图1是overset cell type.
9de13be5-dde7-4e34-ba97-95264d78b5b8-image.png
图2是翼型表面的切应力
5b284d9b-e3bc-45a4-b92e-d749155a96c3-image.png

谢谢各位，解决了:mihu:

在1954年，IBM创立了Fortran语言。

在1955年，多相流中的粒子法就开始进行了。在当时，其被取名为Particle In Cell (PIC)。

PIC由Los Alamos国家实验室的F.H. Harlow提出。在PIC中，拉格朗日粒子通过其坐标（x, y, z）进行定位并且具有质量。有些情况下，拉格朗日粒子还具有多个物理特征。例如在PIC中，拉格朗日粒子可能还具有大小的区别等。

请问大佬们知道怎么用starccm实现传送带（带凸起）的运动吗。
目前能想到的就只有给凸起物一个个建模设置好运动路径，但应该有更简单的方法吧。提前感谢各位大佬。
0ae2d1b3-d3ba-4382-8cac-a470e7232877-image.png

各位前辈，大佬们，最近在学习中遇到了几个问题希望各位前辈指教一下：
目前在大多数文献中我看到气液两相流中的曳力模型系数都是以下形式：
$k=\frac{3}{4}\alpha_d\rho_c\frac{C_D}{d_d}\left|U_c-U_d\right|$
但是最近学习文献的时候在核口的一些CFD代码文章比如SIMMER-III，AFDM，还有一些其他程序里出现了另外一种形式：
$k=A+B*\left|U_c-U_d\right|$
B项一般是跟$C_d$相关，在SIMMER-III中其形式为：
$B=\frac{1}{8}a_i\rho_cC_D,a_i是交界面面积$
而A项给出的解释一般是comes from Stokes law for viscous flow，形式一般是：
$A=\frac{3\alpha_c\mu_c}{d_d}$
比较好奇这种形式的曳力系数是怎么得到的，尤其是A项，在文献里也没看到相关的引用。
是否是由于这些代码一般是用来计算液态金属-气体两相流，因此比较特殊呢？
希望各位大佬不吝赐教，非常感谢！:xiexie:

@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

loop的是面网格上的所有网格。你的需求应该能实现，最简单的办法是，每个面网格中心去读一遍文档，把散点距离自己最近的那个点的速度赋给所在网格。但是会有其他问题，散点疏密是不是和面网格疏密比较匹配，也许一个面网格配多个散点或者反之。