@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());就可以了