CFD中文网

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

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

替代文字

两位前辈我现在在多轴组合运动上存在一些问题，想请教一下你们是怎么设置这种组合运动的吗？如果方便能给我一个类似案例设置文件吗？我的问题是一个球体在一个圆柱流域中公转加自转的组合运动。我目前尝试了设置，但是计算出现了发散，我转速很小的情况下库朗数居然会增加到几万。@卡洛 @冷空气LQW

/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2406 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object dynamicMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dynamicFvMesh dynamicMultiMotionSolverFvMesh; // Specification of rotation around centre rotation1 { solidBodyMotionFunction rotatingMotion; rotatingMotionCoeffs { origin (0 0 0); axis (0 1 0); omega 5; } } dynamicMultiMotionSolverFvMeshCoeffs { domain { solver solidBody; cellZone domain; solidBodyCoeffs { $rotation1; } } ball { cellZone ball; solver solidBody; solidBodyCoeffs { solidBodyMotionFunction multiMotion; multiMotionCoeffs { rotation_1 { $rotation1; } rotation_2 { solidBodyMotionFunction rotatingMotion; rotatingMotionCoeffs { origin (0 0 0.018); axis (0 1 0); omega -50; } } } } } } // ************************************************************************* //

不好讲故事是比较困境了，哎..

算例是twoSpheresGlowinskiMPI

完全 - > 几乎

@jpzhao 有monitor和report的功能，具体看看相应版本的手册哈
捕获.JPG

MomentumParcelI.H

template<class ParcelType> inline Foam::scalar Foam::MomentumParcel<ParcelType>::Re ( const scalar rhoc, const vector& U, const vector& Uc, const scalar d, const scalar muc ) { return rhoc*mag(U - Uc)*d/max(muc, rootVSmall); }

啊客母~

好歹给个工况、示意图之类的瞅瞅:143:

@yfleung 设置x,y的起止范围就行

下图的变量越界了。红的过了上界。蓝色过了下界。
解得纯对流方程，我这几年...一直在和有界和守恒作斗争 :big_mouth: :big_mouth:

......

0_1512380644852_捕获.JPG

是想用cantera算一维火焰做FGM表，用在openfoam计算里ToT

@bestucan 在 关于fluentmeshing画网格 中说：

@hitsc30 我没搞过边界层 y+ 什么的。
但是，对于库郎数大于1的状况，fluent 确实比 OpenFOAM 能算，但是能算是怎么算的就不知道了。

如果是高速流动，边界处很难出现垂直界面的速度，大纵横比就没问题，建筑上吹的风，有个回流什么的就……

当然要直接解决就是减小时间步，特别小的那种，保证在垂直界面方向一个时间步也跨不过一个网格。但是估计会很慢。

snappyhexmesh确实玩的可行性比较高，我画了一个1300万的网格，运行了4h，质量还是不太行。目前对网格有点麻木了，谢谢指导！

Sep 25, 2018 , Jan 08, 2019, Required Reviews Completed

@李东岳 好的老师:xiezuoye:

@simthere
icem默认的节点参数复制是依据目标edge的长度按比例进行复制的

如果第一层网格高度都为0.8的话，可以在edge参数设置页面里勾选Copy Absolute

:cheeky: :cheeky:

问题解决

@李东岳 您好，抱歉，很久没来了。当时在学习FVM，最近在尝试利用openfoam解决NS方程。

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

@李东岳 666:sunglasses:

请问OpenFOAM有关于波浪与柔性结构相互作用的求解器吗？