CFD中文网

只要原来的网格点都在就没问题

@supersoldier 好久没来了，2d的问题进展的如何啊？我最近在尝试做AMR的并行运算，level一高网格数过大nodes之间的数值传达总会出error，我想可能是跟内存大小和nodes之间的通讯设置有关系。你有试过高level的并行运算吗？我的case背景网格用的uniform 1mm,到level6（15.625um）就算不动了，你有什么建议吗？

nevermind,用gMax代替max就可以了

hey there!
我想请教一下关于max（）函数的机理。
以下是我的代码

//----------------------------sensor-------------------------------------//
        volScalarField Uy = mag(U.component(vector::Y));
        volVectorField gradUy = fvc::grad(Uy);
        volScalarField gradUyy = mag(gradUy.component(vector::Y));
        gradUyyN = gradUyy/max(gradUyy);

//        using tensor
        volTensorField gradU0 = fvc::grad(U);
        volScalarField gradU1 = mag(gradU0.component(tensor::YY)); 
        gradUtN = gradU1/max(gradU1);
//
        volScalarField gradrho = mag(fvc::grad(rho));
        gradrhoN = gradrho/max(gradrho);
//----------------------------sensor-------------------------------------//

因为要作为检测器使用，我需要将gradrho和gradU做无次元化处理，使其的值范围在0-1之间。
这里我用了两种方法算了gradU的值，结果显示它们是一样的。

问题：
gradrhoN具有0-1的值范围。
然而 gradU的范围却在 下限：一个很小的值10的负20次方左右 ------上限：0.5左右

下限不为0的倒是可以理解，但是上限为什么不是1？
max/max!=1???????????

谢谢

恩恩，我也希望能早点弄出来，共勉。我现在在用rhoCentralDyFoam做二维的计算，frontAndBack可以用empty边界条件。另外，在4.0以前的版本中，我用的这个solver有很多bug，基本没一个能跑AMR的。至于别的DyFoam我不是太了解。感觉关键就是meshcutting engine的问题，hex4要是能搞出来就算解决了一个大问题了。btw，AMR的sensor的设定也是一个巨大的坑，那个dynamicFvReineDict里的refinement相关的parameter感觉也有问题，比如加密范围的指定等。

你用什么做的网格？如果是icem的话改一下scale factor就可以。

我用vanleer算Shock相关的问题效果感觉还不错的，你算得是关于什么的case？

前面提到的边界条件我也都试过了，没有能防止z被不分割的，如果你成功了能分享一下吗？

关于openfoam 2D AMR的成功的论文目前就三篇，有源代码的除了这篇，还有https://www.academia.edu/16217705/Two-dimensional_adaptive_meshing_in_OpenFOAM这一篇，不过里面的hex4的源代码没有公开，CFDonline上有个人写了hex4的代码公开在github上，不过有很多毛病也用不了。很好奇什么边界条件可以不使z方向被切割。