openfoam 10 中怎么按某场的特征动态加密网格？

JqBu

想在of10里做些工作，但是发现of10的动态网格加密和以前版本有很大不同，求助大家一下：
比如我想按alpha.water的大小加密网格，以前可以写：

dynamicFvMesh      dynamicRefineFvMesh;

但是现在of规则变了，不能这样写了，src里没有dynamicfvMesh了，只有这些是和动网格相关的，但是要不是网格是移动的要不就是要topo的，

请问大家of10里怎么给alpha加密网格呀？

JqBu

@JqBu 找到个of官方算例动态加密网格的，是下面这样，不过在我的算例里还是跑不了，贴出来记录一下把

topoChanger
{
    type            refiner;

    libs            ("libfvMeshTopoChangers.so");

    // How often to refine
    refineInterval  1;

    // Field to be refinement on
    field           alpha.water;

    // Refine field in between lower..upper
    lowerRefineLevel 0.001;
    upperRefineLevel 0.999;

    // Have slower than 2:1 refinement
    nBufferLayers   1;

    // Refine cells only up to maxRefinement levels
    maxRefinement   2;

    // Stop refinement if maxCells reached
    maxCells        200000;

    // Flux field and corresponding velocity field. Fluxes on changed
    // faces get recalculated by interpolating the velocity. Use 'none'
    // on surfaceScalarFields that do not need to be reinterpolated.
    correctFluxes
    (
        (phi none)
        (nHatf none)
        (rhoPhi none)
        (alphaPhi.water none)
        (ghf none)
    );

    // Write the refinement level as a volScalarField
    dumpLevel       true;
}

JqBu

@JqBu 报错如下，贴出来给大家看看，有没有大佬指教一下呢```
Courant Number mean: 0.0440313 max: 0.39994
Interface Courant Number mean: 0.00208161 max: 0.39994
deltaT = 1.54321e-06
Selected 1971 cells for refinement out of 125250.
Refined from 125250 to 139047 cells.
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::fvMeshTopoChangers::refiner::refineUfs(Foam::HashSet<int, Foam::Hash<int> > const&, Foam::polyTopoChangeMap const&) at ??:?
#4 Foam::fvMeshTopoChangers::refiner::refine(Foam::List<int> const&) at ??:?
#5 Foam::fvMeshTopoChangers::refiner::update() at ??:?
#6 Foam::fvMesh::update() at ??:?
#7 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/interFoam"
#8 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/interFoam"

尚善若水

@JqBu 请问这个问题解决了吗？我最近做了点二维AMR，运行也是如此报错，正在debug

李东岳

OpenFOAM没有2D AMR，你是把empty变类型了么，类似一个quasi-2D?

尚善若水

@李东岳李老师，我根据别人支持2D与3D的 AMR代码更新到了最新版本，3D跑着没问题，2D还是报错了

李东岳

我当时好奇你怎么做的2D AMR，就问了你一句。默认的AMR，你要是要跑2D，你把empty面全部改成patch，然后把前后面设置成slip边界，也可以。

不过你那个发散了我不清楚什么情况。

jxzs_zy

请问下现在版本3D四面体可以AMR吗？

尚善若水

@jxzs_zy 在开发，好像还没正式发布。Added dynamicPolyRefinementFvMesh

CFD中文网

openfoam 10 中怎么按某场的特征动态加密网格？