CFD中文网

openfoam3太老了，之前拉格朗日粒子总有这个bug，新版本好了一些。

您好，已经2024年了，我在使用DPMFoam时候也遇到了同样问题，请问您已经解决了么。我目前用的of300，版本太老了，而且昨天使用时候都没有这个报错，今天就改了一下粒子数和初速度就有这个报错。是不是of新版本的拉格朗日的求解器稍微好用一点

感谢分享 楼主好人 :xiexie: :xiexie:

之前在德国的时候，从来没发生过这个现象。

回国之后，用的工作站，配置更强劲，但是经常会掉帧，重装了4 5次也不知道问题所在。后来只能凑合用。用了4个月。

这几天胖猫把他的电脑拿回来，发现连接我的显示器同样掉帧明显。我这才感觉到应该不是工作站的问题，可能是显示器的问题。

但是之前在德国的时候也没问题啊？

有考虑了一下，德国的时候用的是DP线连接的，现在是用的HDMI线。于是京东买了一根DP线。

可算不掉帧了。。。

@李东岳 我先将网格细化做一下测试，排除或者确认一下网格的原因。

我采用的稳态MRF模型进行模拟，如果桨距角改变，同一来流速度下给旋转域设置的旋转速度需要改变吗

@东岳 如果能够通过代码让边界层网格按湍流模型和y+自行调整就好了

@李东岳 最后还是用了后处理的方法进行计算，又要做实验又要做模拟有点赶，用的是老板之前研究的方法几何平均粒径、几何标准差与浓度的关系，算了浓度再来得到颗粒数关系。不过布朗运动要加进去我感觉得用群体平衡方法。

请问解决了吗？也遇到了同样的问题。

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

......

0_1512380644852_捕获.JPG

用STARCCM做船舶阻力分析，考虑升沉和纵倾，之前流域底部设置为速度入口，结果会收敛。考虑浅水效应即流域底部的影响时，我把底部的边界条件修改为壁面 结果不收敛 且提示如图警告-the mesh may have non-positive volume cells，请问前辈是壁面条件不能修改为壁面而是其他，还是修改为壁面后需要再修改其他的地方，谢谢

之前有点问题，最终如下
\par $C_p$采用4次多项式分段拟合
\begin{equation}\label{equ:NS_Cp}
C_p(T)=a_1+a_2T+a_3T^2+a_4T^3+a_5T^4
\end{equation}
\par静焓
\begin{equation}\label{equ:NS_H}
H(T)=\int_{T_{0}}^{T_{x}}C_p(T)dT=a_1T+\frac{a_2}{2}T^2+\frac{a_3}{3}T^3+\frac{a_4}{4}T^4+\frac{a_5}{5}T^5+a_6
\end{equation}
\par熵
\begin{equation}\label{equ:NS_S}
S(T)=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}=a_1\ln{T}+a_2T+\frac{a_3}{2}T^2+\frac{a_4}{3}T^3+\frac{a_5}{4}T^4+a_7
\end{equation}
\par$\bullet$求解静温(已知总温和马赫数)
\begin{equation}\label{equNSUs}
U_s^2=2\left(H(T_{tot})-H(T_{sta})\right)
\end{equation}
\begin{equation}\label{equNSmach}
Mach^2=\frac{U^2}{\gamma(T)R_gT}=\frac{U^2}{\frac{C_p(T)}{C_p(T)-R_g}R_gT}
\end{equation}
\par由(\ref{equ:NSmach})和(\ref{equ:NSUs})得
\begin{equation}\label{equNSTsta}
T_{sta}=\frac{2\left(H(T_{tot})-H(T_{sta})\right)}{Mach^2\frac{C_p(T_{sta})R_g}{C_p(T_{sta})-R_g}}
\end{equation}
$\bullet$求解静压(已知总温、总压和静温)
\par由p等熵过程
\begin{equation}
ds = C_p(T)\frac{dT}{T} -R_g\frac{d p}{p}=0
\end{equation}
\par两边同时积分有
\begin{equation}
\int_{T_{tot}}^{T_{sta}}C_p(T)\frac{dT}{T} =\int_{p_{tot}}^{p_{sta}} R_g\frac{d p}{p}
\end{equation}
\par记
\begin{equation}
S(T_{x})=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}
\end{equation}
\par则
\begin{equation}
S(T_{sta}) - S(T_{tot}) = R_g\ln\frac{p_{sta}}{p_{tot}}
\end{equation}
\par那么
\begin{equation}\label{equNS_psta}
p_{sta}=p_{tot}e^{\left(\frac{S(T_{sta})-S(T_{tot})}{R_g}\right)}
\end{equation}

https://www.2025.unsw.edu.au/apply/scientia-phd-scholarships/explore-particle-scale-nature-drug-transport-blood-vessels-cells

The effective delivery of drug particles or agglomerates from blood vessels to tumour cells is a critical step for their effectiveness in treating the target tumour. This project aims to understand the particle-scale fundamentals in the transport processes by means of unique state-of-the-art numerical techniques. Of particular interests are non-spherical particles transports inside blood vessels, through vessel walls and inside tumour cells, agglomerates breakage, and the effects of key variables. The outcomes will help improve the understanding of underlying micro-scale physical processes that control nanoparticle drug delivery in basic terms, and then provide an effective tool for designing effective drug delivery to targeted regions of clinical importance.

The applicant is expected to have proven research skills including multiphysics modelling and drug delivery related experiments, research experience in drug delivery or biological process, and demonstrated track record of independent research for example published first-author research papers in high-impact journals.

0_1516841984913_捕获.PNG
动图1：河床沙石冲刷运动

河流的流动会改变河床的形状。底部的沙子以及大石头通常会分为两个区：

顶层的bed-load zone（有研究河海的么，这个怎么翻译？），其中的颗粒会随之移动； 底层：沙子以及石头会被困住，但是会缓慢的爬行；

在碎石较多的河床中，较大的沙石倾向于堆积在顶层，这种现象被称之为armoring。个人感觉之所以称之为armoring是因为较大的沙石趋向于在表面，类似一种武装。

实验表明，在这个区域，大沙石有一个净得向上的运动速度，小颗粒具有一个小的向下的移动速度。 目前对于这种现象的原因尚不清楚。不过现存了一些潜在的理论。例如支持程度比较高的巴西果效应（Brazil Nut Effect）。

0_1516841984913_捕获.PNG
动图2：有人晃荡大米里面的密度差不多的小石头

替代文字
动图3：进行的计算模拟

巴西果效应是指如果把两种颗粒的混合物置于容器中，然后施加外加的振荡，体积比较大的颗粒会上升到表层，而较小的颗粒会沉降到底部。关于这种古老的效应的动力学机制至今仍众说纷纭。

DEM模拟：文中作者采用LIGGGHTS对这个现象进行了模拟，下图可以看出，随着时间的推移，大颗粒漂浮出来。

0_1516841984913_捕获.PNG
图4：Nature Communication中的LIGGGHTS模拟

看起来很有意思 :cheeky:

本文来自Nature Communication

@吴优 现在的问题是如何将我的我的结构求解模块和fluent联立起来。是如何获取wall的表面力用于插值到固体域网格。

啊哦，第一次发帖好像还把格式弄错了非常抱歉，，
图在这里
1.png

另外请教一下如何发文件，想发例子的压缩包的貌似也没成功

81a56a41-f2c6-4e5c-8504-2cb7df0a1f28-image.png ，我看帮助文档中DPM对droplet的描述是换热、蒸发和沸腾，这都是液滴材料从液相传质到气相，请问能不能实现从气相到液相的传质（就是吸收）？

@李东岳 这点我清楚，实际上这么多参数里也就总摩擦系数对不上，其他的基本上全对上了。LES的话得在文献里是说在入口用vortex method施加扰动，而且模拟的是有肋的，所以基本上不担心涡的问题。流向周期性的模型也可以在RANS模拟之后用TUI命令施加扰动，之前试过了涡能长期存在，不过那个case没算完结果对不对就不清楚了。