CFD中文网

@zousiyu 这个我并没有遇到过，但是我感觉你的界面不是很光滑，是网格比较粗的原因吗？

@tens 单个液滴破碎问题

用VOF做液滴破碎问题，好几个算例都出现了如下情况，不知道问题出在哪里，请大家帮忙看看。。。

遇到一个二维并行问题很奇怪，用simple方法划分域
采用以下设置可以计算，

simleCoeffs
{
       n        (16 1 1);
       delta    0.001;
}

但是采用这种设置就不能计算

simleCoeffs
{
       n        (8 2 1);
       delta    0.001;
}

而且采用第一种方式时，计算核数划分的较多的时候也不能计算

@影川风 您好，我现在也在做这个耦合的内容，您这边这个deltaX最后有什么好的处理办法吗

@李东岳 李老师，这样是可以的，

("((interpolate(grad(psi))|(mag(interpolate(grad(psi)))+(1e-06|dimChange)))&S)" none)

@李东岳 李老师，这是我自己的求解器，但是这个通量我也不知道为什么会出来
我在dynamicMeshDict加了您说的这个，结果报错，可能是代码的问题。

correctFluxes
(
    (phi none)
    (rhoPhiH none)
    (nHatf none)
    (rhoPhi none)
    (alphaPhi none)
    (ghf none)
    (alphaPhi10 none)
    (phi_0 none)
    (((interpolate(grad(psi))|(mag(interpolate(grad(psi)))+(1e-06|dimChange)))&S) none)
);

[0] --> FOAM FATAL IO ERROR: 
[0] wrong token type - expected word, found on line 56 the punctuation token '('
[0] 
[0] file: /home/oufool/OpenFOAM/oufool-5.0/run/bubblecol/bubblecol9/constant/dynamicMeshDict.correctFluxes at line 56.
[0] 
[0]     From function Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::word&)
[0]     in file primitives/strings/word/wordIO.C at line 74.
[0] 
FOAM parallel run exiting

计算过程中显示：

--> FOAM Warning : 
    From function Foam::autoPtr<Foam::mapPolyMesh> Foam::dynamicRefineFvMesh::unrefine(const labelList&)
    in file dynamicRefineFvMesh/dynamicRefineFvMesh.C at line 546
    Cannot find surfaceScalarField ((interpolate(grad(psi))|(mag(interpolate(grad(psi)))+(1e-06|dimChange)))&S) in user-provided flux mapping table 
7
(
phi none
rhoPhi none
rhoPhiH none
nHatf none
alphaPhi10 none
ghf none
alphaPhi none
)
    The flux mapping table is used to recreate the flux on newly created faces.
    Either add the entry if it is a flux or use (((interpolate(grad(psi))|(mag(interpolate(grad(psi)))+(1e-06|dimChange)))&S) none) to suppress this warning.

我在correctFluxes中加入

 (((interpolate(grad(psi))|(mag(interpolate(grad(psi)))+(1e-06|dimChange)))&S) none) 

但是结果显示格式不正确，有大佬遇到过相类似的情况吗？目前针对这块还不是很熟悉，请大家多多指教

@李东岳 好的，谢谢老师

@tens 可能是我程序问题，我去再检查检查

@李东岳 李老师，我感觉是我的程序的问题，因为有个中间的时间步出现了这种情况

大家有遇到过这种网格加密错位的情况吗？

@fubianhanshu of5,这个需要装一个Swak4foam工具箱，链接如下：http://openfoamwiki.net/index.php/Contrib/swak4Foam

初始剪切流场，funkySetField的实现方法，已测试可使用。所有参数为初始计算域条件，如之前图中所示，Uw = 1.2，H=4，图中计算域长宽比为1：2.。

 expression  
 (
Ui    
{
field u;  
expression "vector(0.3*(pos.y*-4), 0, 0)";  
condition "pos().x >= 0 &&  pos().x <=16" ;
keepPatches true;  
}

@cccrrryyy 好的老师，实在不好意思。我是想到方法了，但是还没测试。我测试好之后会传上来的。

问题已解决，谢谢大家

想请教下大家，这个初始速度场怎么直接给定，我现在的实现方式是先给定上下速度边界条件，再算到流场稳定，有没有一种方法能够直接给定到这样一个计算域中的剪切流场。

@chengan-wang 我也不是太懂您这个方面的，但是从编程角度来看，为什么您不直接定义标量ymax的值呢，而是先定义一个D的值呢

@conghongchuan 你不是用2000万算的吗

@李东岳 感谢李老师回复，我先试下