量纲不匹配

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

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

学流体的小明

@bcp OKK，可以把怎样解决的发出来吗？大家有相同问题的话也可以参考你的方案。

bcp

@学流体的小明后来我看了一下运行日志发现量纲不匹配并不是动量方程的量纲不匹配，是我在定义阻尼区起点的时候没有定义量纲所以会一直报错，后来把阻尼区起点定义量纲之后就可以啦，整个代码是这样的。```

 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",
   dimensionSet(0,1,0,0,0,0,0),
   12
);


volScalarField R
(
       "R",
       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)
    - (fvc::grad(U)&fvc::grad(muf))
    + rho*U*R
);

CFD中文网

量纲不匹配