[strakey]

Has anyone noticed that the species transport in rhoreactingFoam works differently with

mvConvection->fvmDiv(phi, Yi)

than if you use the more conventional

fvm::div(phi, Yi)

and specify the div schemes (limitedLinear 1) one-by-one? I see more diffuse gradients (by about 35%) with the mvConvection approach in a simple shear-layer problem when I have chemical Reactions turned on, even when the reactions are all taking place far downstream of the zone I am looking at. With reactions turned off, the gradients look fine. With fvm::div I don't see that which is consistent with my Fluent results. I also see this behavior in detonationFoam which is a rhoCentralFoam derivative with transport and reactions. Using OF6

Thanks,
Pete