CFD中文网

或者有老师知道如何将ICEM的模型输出为.stp格式吗？

分别将小球和圆柱建立两个body，再采用结构化网格划分球体，再设置交界面节点对应生产圆柱体网格

@evensun 对的

@liuwenchao 感谢大佬的回复，但是出来的id和boundary conditions那里显示的id不一样欸

f4e03b50-4bf0-415a-aff4-7a48433a5f66-image.png 这是boundary conditions那里的id

06a8dd9d-b32e-42ef-a664-954169e4b990-image.png 这是surface list显示的id
:xinlei:

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

......

0_1512380644852_捕获.JPG

\begin{equation}
m_0=w_1+w_2
\end{equation}
\begin{equation}
m_1=z_1w_1\xi_1+z_2w_2\xi_2
\end{equation}
\begin{equation}
m_2=z_1^4w_1\xi_1^2+z_2^4w_2\xi_2^2
\end{equation}
\begin{equation}
m_3=z_1^9w_1\xi_1^3+z_2^9w_2\xi_2^3
\end{equation}
\begin{equation}
m_4=z_1^{16}w_1\xi_1^4+z_2^{16}w_2\xi_2^4
\end{equation}
\begin{equation}
m_5=z_1^{25}w_1\xi_1^5+z_2^{25}w_2\xi_2^5
\end{equation}
左边的$m$都已知，右边的6个未知数$z_1$, $z_2$, $w_1$, $w_2$, $\xi_1$, $\xi_2$需要求。

除了Newton-Raphson迭代方法？

按照如下网址的方法编译安装成功了
https://github.com/wyldckat/FluidStructureInteraction/issues/4

各位大佬们，本人现在研究方向是slamming，能推荐一些书籍看看么，并且最近才接触openfoam，感觉要学的好多！

@浪迹天大 师兄可以加个微信吗？wallong95:laughing:

@luofq-sysu 好的，非常谢谢你！

哈哈，好的

IP访问

注册用户可以自行创建、修订词条，试试看？谢谢

我又回来了，这次发现虽然写了用坐标去计算不同点的阻力系数，但是实际计算中发现还是均一的，以为是定义的变量被重新赋值了，重写了所有的中间变量发现也没办法解决，目前认为是坐标计算的不对或者openfoam porosityModel不支持非均一的阻力系数设置
如果想用场自定义阻力需要如何操作呢？直接改求解器加入源项，然后在0文件夹写初始场？

forAll(cells, i) { gamma_[zoneI][i] = Zero; scalar x = mesh_.C()[i].x(); scalar y = mesh_.C()[i].y(); scalar z = mesh_.C()[i].z(); scalar dir = y; if (xxx) { gamma_[zoneI][i].xx() = xxx; gamma_[zoneI][i].yy() = xxx; gamma_[zoneI][i].zz() = xxx; } D_[zoneI][i] = Zero; D_[zoneI][i].xx() += pow(gamma_[zoneI][i].xx(), 1/3) * dXYZ_.value().x(); D_[zoneI][i].yy() += pow(gamma_[zoneI][i].yy(), 1/3) * dXYZ_.value().y(); D_[zoneI][i].zz() += pow(gamma_[zoneI][i].zz(), 1/3) * dXYZ_.value().z(); D_[zoneI][i] = coordSys_.R().transformTensor(D_[zoneI][i]); // leading 0.5 is from 1/2*rho F_[zoneI][i] = Zero; F_[zoneI][i].xx() += 0.5 * fXYZ_.value().x(); F_[zoneI][i].yy() += 0.5 * fXYZ_.value().y(); F_[zoneI][i].zz() += 0.5 * fXYZ_.value().z(); F_[zoneI][i] = coordSys_.R().transformTensor(F_[zoneI][i]); }

如题，fluent中UDF 能实现气体密度能随位置变化吗，不打开能量方程计算，会不会导致压力计算异常？

笔记本15.6寸的是4k屏，光看2k分辨率或者4k分辨率的时候好像还能适应，就是一切换的时候，感觉2k分辨率简直都是颗粒贼大的点....

您好！问题解决了不？我遇到了相似的问题，系统耦合中只有一条数据传输，是位移，fluent的耦合界面没有输出力。不知您后来怎么解决的呢？

UDF如何读取网格热流（W/m2）...为啥fluent只能通过UDF写入边界热流不能读取呀

@液固两相流 目前我能想到施加的速度函数就是$\nabla\cdot\bfU=0$.

在做液滴撞击的模拟，因为要用到动态接触角所以需要获取三相点位置以及速度。目前是想用两步迭代的位移差来计算速度，但是现在不知道怎么去算这个三相点的位置。代码是这么写的

Thread* thread = Lookup_Thread(domain, 9); Thread** pt = THREAD_SUB_THREADS(thread); cell_t cell; face_t f; real x[ND_ND], area[ND_ND]; double max_x = 0, f_Hoff_inverse, x_hoff, temp, Ca, volume; int n; sum = 0.0; begin_c_loop_all(cell, pt[1]) { if (C_VOF(cell, pt[1]) != 0) { C_CENTROID(x, cell, pt[1]); if (x[0] > max_x) max_x = x[0]; } sum += C_VOF(cell, pt[1]); } end_c_loop_all(cell, pt[1])

我想问一下怎么把SUM变成面积啊