phiGradC = gradPhi(iUpwind,:, iComponent); in src/assemble/NSAssembly/Momentum/cfdAssembleMomentumDCSchemeTerm.m seems incorrect. I changed it to phiGradC = gradPhi(iUpwind, :); because we are only indexing the gradient to get the unwind gradients
phiGradC = gradPhi(iUpwind,:, iComponent); in src/assemble/NSAssembly/Momentum/cfdAssembleMomentumDCSchemeTerm.m
seems incorrect.
I changed it to
phiGradC = gradPhi(iUpwind, :);
because we are only indexing the gradient to get the unwind gradients