已解决，辐射换热的文件里没有将limestone去掉，谢谢大家

phi一般是流场解出来之后，用速度在网格面心插值得到的surfaceScalarField。也就是说，phi是求解出来的量，不能随意改动吧？如果要改温度方程对流项，应该再叠加一个phir，类比interFoam求解器相分数方程的interface compression term。

@yhdthu

感谢，没想到知乎上也有大神研究这个！

@东岳 那个我添加成了，可能是我考虑多了，加和不加区别不大

H. Weller, A New Approach to VOF-based Interface Capturing Methods for Incompressible and Compressible Flow, OpenCFD Ltd., Report TR/HGW/04, (2008).

@东岳 是想获得颗粒轨迹

@dzw05 好的 谢谢 我再试试 看看是不是我设置的问题

paraFoam需要编译一下，个人不是很喜欢paraFoam，启动慢，显示网格也慢。不过有空我更新一下编译方式

你好！可以参考foam-extend版本，相对而言，extend版本的固体力学求解器选择更多。我看了下，extend版本的考虑温度影响的弹性求解器是用的温度差更新应力的。

bounded不建议用于瞬态计算哦

本人目前在用rhoPorousSimpleFoam求解器进行多孔介质模型的计算，对angleDuctImplicit算例进行修改。计算条件：固相提供恒定热流密度边界条件。目的：获得固相与液相的温度场。
现在遇到两个问题：
1 无法设置恒定的热流密度边界条件。算例中给的是恒定温度，如果要设置成热流密度需要定义换热区域（固相&液相）。
2 rhoPorousSimpleFoam求解器中的thermophysicalProperties只能设置流体的参数，无法设置固体参数。

请问各位大神这个问题有什么好建议么

请问老师，VanLeer格式在处理相方程的对流项时，可以保证质量守恒和有界吗？那MULES的作用是什么呢，对于使用对流项使用WENO格式，您有什么看法？求指教

@东岳 东岳老师，问题已解决。非定常计算一开始要求时间步长要比较小，先得到一个较稳定的流场，在别的求解器里可能会先用稳态求解器算一个初场。而我上面的算例一开始时间步长较大，到后面越来越不稳定最终发散。开始用变时间步长让对应的库郎数从1以下慢慢增加即可。

我看了下，确实不是连续的。在ManualInjection.C文件中：

template<class CloudType> Foam::scalar Foam::ManualInjection<CloudType>::timeEnd() const { // Injection is instantaneous - but allow for a finite interval to // avoid numerical issues when interval is zero return ROOTVSMALL; }

我看了一下，可以用cellZoneInject：

Description Injection positions specified by a particle number density within a cell set. User specifies: - Number density of particles in cell set (effective) - Total mass to inject - Initial parcel velocity Properties: - Parcel diameters obtained by PDF model - All parcels introduced at SOI

你试试这个？

PS. 你需要用topoSet做一个cellSet出来toturial里面有

@zhangxc0223 可以用cyclic边界

@东岳 成功了，谢谢东岳老师了

@东岳
在interfaceProperties.C中编程如下
fvc::interpolate(fvc::grad(sigma) - fvc::grad(alpha)*(fvc::grad(alpha)&fvc::grad (sigma))*mag(fvc::grad (sigma))) & mesh.Sf()
将此项添加到速度和压力方程。其中当计算压力项时发散了。

@东岳 好的，明白了。谢谢李老师！

//- Breakage kernel SBre_[kth].ref() += breakagePtr_->g()[i]()*weis_[i] *( DSDPtr_->beta(abs_[i], kth)() - pow(abs_[i], kth) ); //- DSD TO DO implemented //- Coalescence kernel if (coalescence_) { for (int j

上周五自己做的评论，DSD函数还没植入，周一来的时候完全忘了。吭哧吭哧算了好几天，感觉没有破碎的效果呢？

检查代码发现函数还没植入，返回的0..

对于不可压缩流动，密度/压力的大小是无关紧要的。重要的是压力差。因此你不需要保证求解的密度是你要的密度。

\begin{equation}
\frac{\partial \rho \mathbf{U} }{\partial t}+\nabla \cdot \left(\rho \mathbf{U} \mathbf{U} \right)-\nabla \cdot \left(\mu \nabla \mathbf{U} \right)=-\nabla p
\label{NS2}
\end{equation}

比如这个方程，如果密度不变，你直接除掉密度就可以。因此在不可压缩求解器中你不需要设置密度的大小