@izumi 谢谢啦，感激不尽

双流体和多流体求解器中需要设置的各相的直径对模拟结果有什么影响

直径直接关系到各种力的计算。是计算耦合力的重要因素。并且对传质至关重要。

需要以什么样的规则设置？

直径不能任意给定，需要参考实验值严格给出。

直径大小的设定好像跟网格尺寸有关，

http://www.sciencedirect.com/science/article/pii/S1674200117300524

看这个文章，讨论了直径的限制。目前拉格朗日那面限制的更为严格。但不能反过来说直径的大小设定和网格有关，而是网格的生成应该取决于直径。

可以使用cyclicAMI边界条件

首先你需要做两套网格，可以要么用第三方网格软件来做，要么用openfoam自带的topoSet之类的来做，两个网格最重要的是，上图中的交界面要有边界。可以是wall，也可以是patch，反正后来要改。

导入到OpenFOAM之后，要把这个交界面设置为cyclicAMI，这个可以参考pimpleDyMFoam下面的propeller算例的设置过程

滑移网格是滑移和插值的过程。就不存在网格负体积了。你的运动发散是因为变形过大产生了负体积。

0_1528241930740_无标题.png

下面是在Ansys Fluent理论指南中遇到的部分描述：

The analytic scheme is very efficient. It can become inaccurate for large steps and in situations where the particles are not in hydrodynamic equilibrium with the continuous flow. The numerical schemes implicit and trapezoidal, in combination with Automated Tracking Scheme Selection, consider most of the changes in the forces acting on the particles and are chosen as default schemes. The runge-kutta scheme is recommended of nondrag force changes along a particle integration step.

还需要进一步研究。有人见过相关文献么

@肖恩曼
那可能就需要用循环一个个赋值了

请问您用的是哪个求解器对应的算例呢？风机的转速是随入口风速变化吗？

@linhan-ge 嗯嗯好的，谢谢您的回复！

of里面怎么设定一个壁面的温度沿着壁面方向线性变化啊

@东岳 @dzw05 感谢两位，明白了，那基本一样。

自己回答吧，
这个跟coblended还不一样，coblended的blendingfactor 得另外定义一个volScalarField。
但是fixedblended就加个系数就行了。

fvm返回的都是矩阵，矩阵又分为矩阵系数和矩阵源项（右边那一行）。fvm::Su返回矩阵的源项，fvm::Sp返回矩阵系数。fvc返回的都是场，比如标量场。

你的问题：fvc::Su(T, k)返回T，fvm::Su返回矩阵源项。fvm::Sp返回矩阵系数， fvc::Sp(T, k)返回Tk。fvc::Su(T, k)=T。

是的。VIM emacs ide等都需要一段时间学习。用自己习惯的就好。虽然我用VIM，但在OKS上我都用gedit，要不然别人不知道你在干什么，哈哈

@东岳 多谢东岳大神回复。之前可能我没表述清楚我的问题：我想同时fix液滴与固体表面的接触线和接触角，即壁面处的volume fraction和theta都始终保持t=0时的值，并不随时间而改变。相当于fixed value和constantAlphaContactAngle两个边界条件的结合。如果边界条件设置为constantAlphaContactAngle(zeroGradient也相当于对应theta0=90的constantAlphaContactAngle边界）, 在迭代的过程中会计算当前的curvature, 并通过与设置的theta0比较来更新壁面处的受力情况，从而重新计算边界处的volume fraction, 以使计算的theta趋近于所定义的theta0（如后面摘自interfaceProperties中的代码所示），这并不能保证接触线（即volume fraction）的固定。
我的理解是这样的，不知道对不对：接触角theta与壁面处volume fraction的梯度是相对应的。能否让壁面处的第一层和第二层网格中的volume fraction值都固定，即volume fraction方程在该区域不进行求解而是始终采用t=0时刻的初始值，这样就能实现我前面所说的目的了。不知在openfoam中能否实现。
期待您答复，谢谢！

forAll(boundary, patchi) { if (isA<alphaContactAngleFvPatchScalarField>(abf[patchi])) { alphaContactAngleFvPatchScalarField& acap = const_cast<alphaContactAngleFvPatchScalarField&> ( refCast<const alphaContactAngleFvPatchScalarField> ( abf[patchi] ) ); fvsPatchVectorField& nHatp = nHatb[patchi]; const scalarField theta ( convertToRad*acap.theta(U_.boundaryField()[patchi], nHatp) ); const vectorField nf ( boundary[patchi].nf() ); // Reset nHatp to correspond to the contact angle const scalarField a12(nHatp & nf); const scalarField b1(cos(theta)); scalarField b2(nHatp.size()); forAll(b2, facei) { b2[facei] = cos(acos(a12[facei]) - theta[facei]); } const scalarField det(1.0 - a12*a12); scalarField a((b1 - a12*b2)/det); scalarField b((b2 - a12*b1)/det); nHatp = a*nf + b*nHatp; nHatp /= (mag(nHatp) + deltaN_.value()); acap.gradient() = (nf & nHatp)*mag(gradAlphaf[patchi]); acap.evaluate(); } }

@东岳 嗯，的确是很麻烦，谢谢李老师，后期准备试试。

同问，有类似问题 http://cfd-china.com/topic/1575

湍流模型kOmegaSST是低雷诺数模型？在近壁区采用k-omega模型，在湍流区采用k-epsilon模型。
那在计算的时候应该采用什么样的边界条件？
我的case是这样的，网格较密，yplus在1附近，雷诺数为7000，存在两个壁面(wall)
我参考的是researchgate中学者的设定方法，具体方法见链接(https://www.researchgate.net/post/Should_the_RANS_model_kOmegaSST_be_fully_resolved_at_low_Re-number_Is_it_possible_to_use_it_with_a_wall_function)
方法有一定合理性，但是对于我的算例可能存在一定的问题，不知大家在应用SST模型时，有没有什么经验，如何设定边界条件？

@东岳 谢谢李老师，依照您的方法我已成功编译，但是log文件内输出的是所有边界的数据（例如icofoam-cavity这个算例输出的是靠近fixedwall和movingwall那一层的网格的速度），有没有方法对只输出靠近fixedwall那一层网格的速度。

Have no idea :chigua:

问题解决了么？是不是已经解决了？..

@ge-zhang
正常的壁面，下个月我会发布算例和代码。
你的头像很有意思啊。