incompatible fields for operation

[danao]

在interFoam方程加入一个加速度aB.

volVectorField aB
(
    IOobject
    (
        "aB",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
mesh,
A*2*pi/T*Foam::sin(2*pi/T*runTime.value())*vector(1,0,0)
);

fvm::ddt(rho, U)  + fvm::div(rhoPhi, U)+... fvm::Sp(rho, aB) 
  == fvOptions(rho, U)

一开始计算就报错：

incompatible fields for operation
 [U]+[aB]

请问该如何处理？

[danao]

源项改为显式可以计算，但结果不太对。