如何在流场中加入体积力场

crispy

体积力场是随时间变化的量，是(x, y)的函数

学流体的小明

编程。

// 先创建一个场
volVectorField source_Vector
(
    IOobject
    (
        "source_Vector",
        mesh.time().timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedVector(dimensionSet(1,-2,-2,0,0,0,0), Zero)
);

// 根据你的需要进行赋值
forAll(mesh,cellid)
{
       source_Vector[cellid]=........;
}

// 添加方案一，添加到UEqn.H当中的动量方程中
 fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U) + fvm::div(rhoPhi, U)
      + MRF.DDt(rho, U)
      + turbulence->divDevRhoReff(rho, U)
     ==
        fvOptions(rho, U)
      + source_Vector // Add the source field
    );

// 添加方案二，添加到泊松方程里面，需要先插值为surfaceScalarField
surfaceScalarField source_phig=-fvc::interpolate(source_Vector );
// 然后添加到pEqn的phig
surfaceScalarField phig
    (
        (
            mixture.surfaceTensionForce()
          + source_phig // Add the source field
          - ghf*fvc::snGrad(rho)
        )*rAUf*mesh.magSf()
    );

我这个是改的interFoam，你类比一下就行，量纲应该不会一样。

crispy

@学流体的小明谢谢大佬，我昨天写好了。但是目前的问题是体积力场添加的位置不对，我现在改成了最简单的模型来检查错误，希望在一个长方形区域中添加固定数值的体积力场，还是不能实现。

crispy

@学流体的小明

volVectorField bForce
(
    IOobject
    (
        "bForce",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedVector("bForce", dimensionSet(1, -2, -2, 0, 0, 0, 0), Zero) // dimless
);

for (label cellI = 0; cellI < mesh.C().size(); cellI++) // only process 0?
{
    const scalar x0 = 0;
    const scalar y0 = 0.0245;
    const scalar x = mesh.Cf()[cellI][0] - x0;
    const scalar y = mesh.Cf()[cellI][1] - y0;

    if (x > 0 && x < 0.012 && y > 0 && y < 0.0005) {
        bForce[cellI][0] = 50;
    }

crispy

@学流体的小明
红线包围的是我限制的区域，黑线表示的是网格，这样看应该有满足限制条件的网格，但计算得到的结果是bForce全是0，没有任何赋值的情况。

学流体的小明

@crispy
mesh.Cf()好像是网格面心，你改用mesh.C()试试。

// only process 0?

这个无所谓的，所有processor都会执行这串代码，但是它们各自的网格中心的坐标是不一样的，不会每一个processor都会有一个(0,0,0)的坐标这样。

crispy

@学流体的小明
1、谢谢大佬！就是这个原因，mesh.Cf()改成mesh.C()就没问题了！
2、当时注释这个是因为我发现有问题后用Info输出坐标和mesh的size想看看问题出在哪，但是只输出了processer0的网格信息，所以有了这个疑问。后来我查到了是Info只会输出master processor的信息，我改成Pout就可以输出所有processor的信息。

CFD中文网

如何在流场中加入体积力场