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}

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

我目前不清楚怎么用paraview处理这个
但我感觉paraview应该可以弄，可能你得外挂paraview的代码脚本，我不太熟悉这个

我也忘了..

这个site是什么意思

目前我也不清楚 知道了可以告诉我一下 :mihu:

@youmengtian 能请教一下你是怎么测试tut中的算例的吗？我计算算例的时候先用scalarTransportFoam计算出一个时间段的结果，再使用PODSolver就报错了。请问你是怎么计算的呀

@东岳 谢谢东岳老师

@东岳 谢谢东岳老师！

分离区大小的问题一般是和粘性有关。看一下用的是什么湍流壁面函数，个人经验，湍流壁面函数对分离区的预测通常是不准确的。我计算的时候一般y+小于1，不采用壁面函数（当然在openfoam里是所谓的低雷诺数模型）。

不太清楚这个，只能建议看看源文件了..

@cfdngu 这是哪篇文献，想拜读一下:chitang: