CFD中文网

@程迪 666啊老铁，听君一席话，胜读百篇论文:146: :146: :146: :146: :146:

@cfd-china :sunglasses: ，其实我也懒得搞了。。我就问问:cheeky:

自己写的代码终于可以出结果了，但是用cgns存坐标数据的时候发现只能存一半，很奇怪，代码举例如下
三维的面网格，节点数2x11，因此是1x10个网格, 坐标储存在x,y,z数组中

real, dimesion(2,11) :: x,y,z CALL CG_BASE_WRITE_F(INDEX_FILE, BASENAME, 2, 3, INDEX_BASE, IER) ISIZE(1,1) = 2 ISIZE(2,1) = 11 ISIZE(1,2) = 1 ISIZE(2,2) = 10 ISIZE(:,3) = 0 CALL CG_ZONE_WRITE_F(INDEX_FILE, INDEX_BASE, ZONENAME, ISIZE, Structured, INDEX_ZONE, IER) CALL CG_COORD_WRITE_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, RealSingle, & 'CoordinateX', x, INDEX_COORD, IER) YZ坐标同上

读取时代码为

real :: node(2,11) node = 0.0 call CG_COORD_READ_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, 'CoordinateX', RealSingle, (/1,1/),(/2,11/),node,IER ) write(*,*) node

node的读取结果为

x(1,1) x(2,1) ... x(1,6) 0.0 0.0 ... 0.0

请问有没有大佬能指点一下

$\tau$这个雷诺应力，对角线的是normal stress，非对角线的是shear stress。不过你说的更详细，感谢分享！

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

同意之前的回复，还是要看实验结果。可能的误差包括三方面：

本身物理模型是否准确。比如Species Transport, 湍流模型以及反应方程，都有一定程度的近似。 离散误差。由于离散化的计算格式（N阶），会有人工粘性的存在，对于计算会有误差。另外网格的大小和质量，虽然可以做网格无关性，不过还是存在误差的。 边界条件的误差。这一点很少有人强调 ，不过由于实际条件的限制，不同壁面条件（粗糙度），以及入口湍流度等，和实际也是有一定出入的。

由于这三方面的原因，所以还是离不开实验的验证。

@李东岳 好的，谢谢老师

这个放在他们搞涡那面，估计能讨论一上午 :duang:

我个人也不太熟悉。在我的文章里面，也不太关心涡的情况。

最近要做一些跨临界的问题，但是在跨临界中，如下图，由于伪沸腾的影响，流体的热力学性质会在伪沸腾区呈现剧烈的变化，形成类似于激波的区域。但是在亚临界区域，压力对密度的变化极为敏感，所以传统的全守恒格式会引起压力震荡。很多研究者针对这个问题做了很多修正。比如：求解压力传播方程并加入人工扩散、计算左右通量求黎曼解并对最后的结果熵修正。但是目前的算法总体来说很繁琐，且由于通过压力反推能量会引起能量不守恒。主要有两个问题，非常感谢您的解答：

求解压力传播方程的那个算法很类似simple算法的pressure based方法，为什么现在几乎所有的计算空气动力学论文里面几乎没有使用过simple算法，更倾向于TVD RK这种显示迭代？ 显式格式很适合GPU加速计算，想开发一套GPU的代码用来计算。但是GPU对双精度计算的支持很差。当在激波附近使用高阶离散格式的时候如果使用单精度浮点数计算，截断误差是否会对计算结果产生比较大的影响？
3.png

@yhdthu 好的 谢谢大佬

@望山跑 多谢，已经解决了。。

他这个聊天记录我没太懂。什么是拖布结构？

你好，请问现在有basilisk的交流群吗，

@coolhhh 发现问题了，因为我的是两相流，需要设定field 为U.water

@手艺白胶
立方式的加倍？是怎么加倍的？:confused:

谢谢啦

@望山跑 好吧，多谢啦

这个我不大会，看编程算的同学总用~不过最近开新坑，当然可以搂草打兔子一起学了

工科末流985硕士毕业，想在武汉找CFD相关的工作，家人们有没有推荐的公司。

@cfd-china 哈哈，没有多久不见吧:cheeky:

操作压力（operating pressure）是不是就是参考压力？ 0_1466035493090_2016-06-16_080420.png 我在模拟的时候，都把操作压力设置成了零，这样得到的压力都是绝对压力？不用怎么换算了，这样没事吧:confused: 《ANSYS Fluent理论指南》已完成约90%，同步校对中什么时候能够完稿呢？