量纲不匹配

bcp

求助各位大佬，我想要在interFoam求解器的动量方程加阻尼消波项，修改UEqn.H，编译可以通过，运行算咧的时候报错量纲不匹配，改了好多次一直都是这个问题，是我定义的有问题吗

dimensionedScalar nu1
(
    "nu1",
    dimensionSet(0,2,-1,0,0,0,0),
    1.48e-05
);
dimensionedScalar nu2
(
    "nu2",
    dimensionSet(0,2,-1,0,0,0,0),
    1.0e-6
);
// Need to store rho for ddt(rho, U)
volScalarField rho
(
    IOobject
    (
        "rho",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT
    ),
    alpha1*rho1 + alpha2*rho2
);

volScalarField nu
(
    IOobject
    (
        "nu",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT
    ),
    alpha1*nu1+alpha2*nu2
);


volScalarField muf
(
    IOobject
    (
        "muf",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT
    ),
    rho*nu
);
![3.png](/assets/uploads/files/1734663338888-3.png)

 dimensionedScalar r0
    (
        "r0",
        dimensionSet(0,0,-1,0,0,0,0),
        0
    );

    dimensionedScalar r1
    (
        "r1",
        dimensionSet(0,0,-1,0,0,0,0),
        5
    );

    dimensionedScalar x1 = max(mesh.Cf().component(vector::X))-min(mesh.Cf().component(vector::X));

    dimensionedScalar x0
    (
       "x0",
       2.4
    );

    volScalarField rField
    (
        "rField",
        max(r1*(mesh.C().component(vector::X)-x0)/(x1-x0),r0)
    );

    fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U) + fvm::div(rhoPhi, U)
        - fvm::laplacian(muf,U)
        - rho*U*rField
    );

学流体的小明

UEqn当中各项的量纲是 $[\mu][U]/L^2=M/(L^2T^2)$
rField的量纲应该是$T^{-1}$才对，但看起来你写出来的并不是这样
你可以在定义场的时候显式地指定量纲

volScalarField B                                            
(
    IOobject
    (
        "B",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedScalar(dimless, 0.0)
);

bcp

@学流体的小明感谢大佬回复，我是这样定义的，我在rField之前定义了r1和r0，这两者的量纲都是（0 0 -1 0 0 0 0），然后rField是根据r1和r0通过公式求得的，这种情况下rFiedl的量纲和r1和r0相同吗

学流体的小明

此回复已被删除！

学流体的小明

@bcp 我搞错了，之前没看对你计算rField的公式。我以为是你rField的计算公式不对。

你把报错代码发一下。
debug，分别输出rho、U、rField的量纲看一看

bcp

@学流体的小明感谢大佬，问题已经解决啦

CFD中文网

量纲不匹配