@wwzhao 不是,我程序中有一句写的不对
const vectorField& Cf = patch.Cf()/delta1VD;
vectorField meanVelocity(Cf.size());
我改成:
vectorField Cf = patch.Cf()/delta1VD;
vectorField meanVelocity(Cf.size());
或者
tmp<vectorField> tCf = patch.Cf()/delta1VD;
const vectorField& Cf = tCf();
vectorField meanVelocity(Cf.size());
就可以了