CFD中文网

@yfclark 可以去试试，最近我全天都在上课呢:zoule:

基类fvPatchField<Type>的snGrad()定义就是返回face normal gradient... 也就是$\nabla_{\perp}\phi=\nabla\phi\cdot\vec n $

已解决 谢谢啦

周老师，这里我测试发现好像这种结构网格的加密方式对于三维模型可以，对于二维网格，在经过局部加密后再导入到OF里面是会报错的，错误是：--> FOAM FATAL ERROR:
face 11528151 does not belong to cell 210100. Face owner: -1 neighbour: 210099

From Foam::cellShape Foam::extrudedQuadCellShape(Foam::label, const labelList&, const faceList&, const labelList&, const labelList&, Foam::label, Foam::faceList&) in file extrudedQuadCellShape.C at line 105.

FOAM aborting
感觉好像是在z方向拉伸出来的面没有识别到具体的cell中

国内形势变化太快了，一年一年的水涨船高，就像我们看现在的教授，当年竟然博士毕业就留校了。以后的年青一代看我们，估计是“这家伙当年竟然博士发了5篇文章就留校了，现在我们博士发10篇都留不下来”..

@队长别开枪 您好，我在OF8中编译了ccm26ToFoam，想模拟波物相互作用，转化后的网格checkMesh也是ok的，转化之后的网格我放到OF v1912里面进行计算的(还没编译ccmToFoam)，但是相同的设置用OpenFOAM自带的sHM划分的网格可以计算，ccm26ToFoam转化的网格却在0.02s就突然发散了。想问一下有什么解决办法吗？:135:

使用的软件不限，主要做算法还是实现也不限。课题组自己一个人搞这东西真的头大。要是有了解的拜托拉一下，多谢啦！

@youmengtian
你用的什么显示器 :cheeky:

@cloud1990 在 VOF求解自由水面问题，求解空气区域吗？ 中说：

Image

Left: Correct jet shape predicted by TruVOF technique used in FLOW-3D. Right: Incorrect jet shape predicted by pseudo-VOF technique used by other CFD codes.

Image
Left: FLOW-3D‘s TruVOF technique predicts jet impingement on wall and some outflow. Right: Pseudo-VOF methods don’t predict realistic jetting of fluid on side walls.

从结果上看，flow3D 所谓的TruVOF是有问题的，air velocity 和　air resistance = 0?　这个可能比较适合粗网格计算，作为一个tricky的存在．但理论上，这个TruVOF是无法解决wind wave的问题．

@东岳 @wwzhao 谢谢两位老师 哈哈，，，我按照东岳老师的去试试

求教岳神，我想做反应熔渗的仿真，涉及及多孔介质毛细现象，碳与硅液生成碳化硅的化学反应，反应后又会填充多孔区域孔隙，结果需要填充过程趋势，有好的方法做这方面的仿真吗

http://www.sciencedirect.com/science/article/pii/S0021999117301493

2011年Vikas提出来一种稳定的高阶矩对流格式之后，本文作者重新提出来一个新的数值格式。

新技术刚追完，马上又出新的。。

@东岳 您好！cyclicAMI和renumberMesh之间还存在一些bug，我将cyclicAMI改为cyclic，还是类似的操作，就可以了。:ok: 需要注意的是：Pointwise导出含周期性边界的网格时，需要先将周期边界设置为patch，然后在openfoam中用createPatchDict将其设为cyclic，就能renumberMesh了。当遇到下面的错误时，提高boundary文件的matchTolerance的值就可以了，不过网格也不能太差。

--> FOAM FATAL ERROR : face 0 and 65 areas do not match by 0.000222148% -- possible face ordering problem From function cyclicFvPatch::makeWeights(scalarField& w) const in file meshes/fvMesh/fvPatches/derivedFvPatches/cyclicFvPatch/cyclicFvPatch.C at line 62.

@东岳 其实我是想看一下加了重力会不会对CdA和Cp做成影响 如果没有多大影响的话 就可以不考虑重力了

当你的学生肯定很幸福

@李东岳 最近在潜水

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

不懂Fluent，不过EMMS类似这些模型的Cd，其实都大差不差。我只能点评到这了..

@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

BTW, you can SAVE 20% when you order online and enter Promo Code ENG19

You can order it here: https://www.crcpress.com/Multiphase-Particulate-Systems-in-Turbulent-Flows-Fluid-Liquid-and-Solid-Liquid/Podgorska/p/book/9781482235685

:baobao:

批注 2019-10-20 163534.png