OpenFOAM有方法能够使一部分网格不参与计算吗？

cccrrryyy

所以说你这个做的应该不是这个体积力的东西，是另外一种了。有什么参考文献可以分享一下么？我很好奇别的处理是什么样子的

Samuel-Tu

@cccrrryyy 我现在做的就是体积力的方法，我在想有没有其他办法。。

我是河滩

@马乔在 OpenFOAM有方法能够使一部分网格不参与计算吗？中说：

set the value of p and U as 0.

可以给出一些有关矩阵修改的参考文献吗？目前也需要对一些给定单元赋值（通过周围正常求解网格单元的值进行插值），OpenFoam中的fixedinternalValue边界条件也有类似功能，但不是很理解。

random_ran

(Sorry, I don't have chinese input right now)

I once had the exactly same issue!

Check my thesis Chapter 6.3 :

https://scholar.uwindsor.ca/cgi/viewcontent.cgi?article=8585&context=etd

I wrote a paper in this: https://ascelibrary.org/doi/pdf/10.1061/9780784415153.ch06

Figure 2 is the original plate. (Sorry due to copy-right reason, I cannot put it here.)

In my thesis, I had to open the plate with some holes in ICEM. So, I maual dig out the block out of the computation domain.

It's really a painful process.

But it works!

This is the original plate:

This is after the "dig":

李东岳

        fvVectorMatrix UEqn
        (   
            fvm::ddt(U)
          + fvm::div(phi, U)
          - fvm::laplacian(nu, U)
        );

        // 固定值
        labelList cells = mesh.cellZones()["null"];
        vectorField value(cells.size(), Zero);
        UEqn.setValues(cells, value);

        if (piso.momentumPredictor())
        {
            solve(UEqn == -fvc::grad(p));
        }

可以通过设定一个null区域，然后对区域的网格点强行赋值，但就像我之前说的，总觉得边界处需要处理。也就像 @cccrrryyy 说的，变成了浸没边界法。这个代码只是大体上是个思路，需要更深的研究。可以看看fvOptions里面一些源项的处理，都是操纵矩阵的，比如fixedTempresure，meanVelocityForce等，

Samuel-Tu

@东岳嗯嗯，我试一下看看能不能行。我记得setVAlue是切断了相邻面上的联系。我处理下剪切力试试

Samuel-Tu

@random_ran 主要的，我的边界随时要改变。这样的话就会一直改网格，有点像自适应网格了。。

OItoCFD

@cccrrryyy 他现在这个方法属于cut cell方法的一种吧我也今天刚要尝试在固体切割网格的地方要特殊处理然后当固体是运动时候为了保持质量守恒要对连续性方程修正一下减去一个固体相对速度详情见这篇论文：A fixed-grid model for simulation of a moving body in free surface flows

李东岳

最近把理论这部分弄了一下，

捕获.PNG

CFD中文网

OpenFOAM有方法能够使一部分网格不参与计算吗？