CFD中文网

我想建立一个水槽，水流从水槽一边冲过去，然后又从另一边冲过来，两边的入口和出口边界要随时间变化，有没有大神指导一下，怎么实现？谢谢！

@禾斗石开 fluent

页数按照pdf的页数，不是书上的页数

上册：

https://www.jianguoyun.com/p/DVcobNoQ9s3ZBhjejV8

73页【直线运动的速度】：导数的直观描述，CFD方程这面全是导数，大部分为偏导数，后面会介绍 75页、76页【公式1-4】 138页【公式3-3】：泰勒公式，CFD很多方程都是基于泰勒公式推导出来的

下册

https://www.jianguoyun.com/p/DQnanmcQ9s3ZBhjijV8

23页：流量计算公式 71-75页：偏导数，CFD这面都是偏导，重点看一下，理解一下 77页提到了个拉普拉斯方程，随便扫一眼 78页全微分，跟CFD的物质导数关联 79页3-3 81页有个例子 112页介绍了梯度，但是感觉描述的不是很简单，可以看看 141-144页 重积分，为之后理解有限体积法做准备 224-225页 曲面积分的定义，为之后理解有限体积法做准备。其中225页有个图，比较形象 230-231页，看高数中对通量的定义 243-245页，高数中明确的对速度散度的定义

我这边需要仿真一个构件在烘箱中逐渐被加热的过程，结构非常简单，就是一个构件周围用壁面包裹，其截面如下图所示。
IGBT fixture-温度2.jpg
四周壁面设定为固定温度358K，初始温度为293K。由于工作中计算资源有限，且目前只需要定性得到构件各部位升温的快慢情况，因此仿真时只打开了能量方程，上图就是仿真至380秒左右时的状态。

然后在仿真中，初始的时间步长是0.001s，因为太小了导致仿真进度比较慢，而且能量方程的残差一直保持在10e-7级别，我这边通过不断地放大时间步长，至380s时时间步长为0.1s，至400s左右时时间步长改为1s。在这个过程中，我发现每次增大时间步长，监测点的温升曲线斜率都会变小，也就是说温升变慢了，比如下图就是时间步长由0.1s改为1s时，监测点温升情况的变化：
IGBT fixture-温升.jpg

在我更改时间步长的过程中，方程的残差均保持在10e-7级别，那么我想知道，温升速率为什么会因为时间步长的增加而变小呢？是不是因为网格不够精细造成的，因为为了快速出结果，网格确实比较粗糙，还是说还有其他可能的因素？谢谢大家了

请问我想使用dpm模型模拟液滴的凝固，应该选什么类型得到颗粒呢？我看几种类型的粒子没有考虑凝固过程，求教！！！

各位大佬，我要使用fluent进行山地管道泄露扩散研究 。现在不知道怎么能将山地地形DEM高程数据，导入到icem中，请教各位有没有走过这个流程的，（在arcGis中应该转化成什么格式）这应该是怎么个搞法。

@李东岳 wifi用5G就不会干扰了啊

张量标识预习测验

下面的方程取自各个不同的CFD教材。方程并没有采用代码的形式而是直接截图，这样更加直观的可以感受到文献里面方程这种多样性的书写方式。各位老师和同学需要能够拆分为分量形式。

0_1529193840362_1.JPG

0_1529193880752_2.JPG

0_1529194030583_3.JPG

0_1529194061311_4.JPG

0_1529194127887_5.JPG

0_1529195669196_6.JPG

0_1529195734307_7.JPG

0_1529195783721_8.JPG

0_1529196217099_9.JPG

0_1529196245498_10.JPG

0_1529196323392_11.JPG

0_1529196403565_12.JPG

0_1529196500098_13.JPG

0_1529197325295_14.JPG

这个是理论大坑啊。

如果压力基的求解器，求解的是速度和压力场，那么得用连续性方程来求压力。
如果是密度基的，连续性可以求密度，然后用状态方程求压力。

@bestucan 这个混乱是可以解决的，在引用的库文件前面加入路径就行了

-L$(FOAM_USER_LIBBIN) \ -llagrangianIntermediate \

而且我在其他论坛上也看到，lagrangian的代码修改时必须把名字改成一致的，不然会出现很多warning，

虽然我知道我弄不成（硬盘没有没分区的剩余空间了），但我还是想试试，然后根目录炸了:jingya:
我说我的根目录怎么一点空间都不剩了，还报警。生成了一个 swapfile 占满了。操作有风险，尝试需谨慎。

这是占根目录的空间当作swap用，绕过了重新分区，挺好，就是空间没的有点……突然 _(┐「ε:)_

各位大佬，最近想用OF计算固定的多体装置，由于在大域内生成的结构过于粗糙 需在小域里生成 然后用mergeMesh合并到大域里去，但是小域的边界没办法删除
请问各位大佬，有没有什么设置建议？删除边界或指定特殊的边界？

@东岳 这个边界条件已经整合到exteranlHeatFluxTemperature里面了，算的也是温度梯度，也就是$q=-\kappa \frac{\partial T}{\partial n}$，这种边界按照FDS技术手册里面应该是DNS的计算方法。
LES计算时给定的方法是$q=h \Delta T$,$h$是对流换热系数。openfoam里面的边界条件要么是固定值要么是梯度，这种差值不知道在of里面能不能实现

周老师，这里我测试发现好像这种结构网格的加密方式对于三维模型可以，对于二维网格，在经过局部加密后再导入到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中

为什么我这个推导的反了？应该是大于0 。

0_1529644803725_微信图片_20180622131807.jpg

@西湖冷月 这就是天秀

今天解决了这个问题
先在$WM_THIRD_PARTY_DIR编译paraFoam并执行wmRefresh
到$WM_PROJECT_DIR/modules/visualization路径下执行 ./Allwclean && ./Allwmake && wmRefresh
转到$WM_PROJECT_DIR路径下执行 ./Allwmake && wmRefresh

命令：

cd $WM_THIRD_PARTY_DIR
./makeParaView
wmRefresh

cd $WM_PROJECT_DIR/modules/visualization
./Allwclean && ./Allwmake && wmRefresh

cd $WM_PROJECT_DIR
./Allwmake && wmRefresh

有大佬用cfx-pre算过二维旋转对称模型吗？我用二维旋转对称算的结果跟实验值相比差了蛮多，但用三维模型算的结果就对的上，但三维模型计算时间太长了，想用二维算。

@fireztw

https://openfoamwiki.net/index.php/Parameter_Definitions_-_dynamicMotionSolverFvMesh

https://forum.cfdac.com/t/topic/660/2

这两个网站可以参考一下~~

@XIN 请问您解决了吗，我也遇到这个问题了