I recently noticed the flux boundary condition in the new version is different from the old versions. For the flux boundary n.j=j_0 is j in the domain always considered as the diffusive j=-D(grad c) in the newer version? If I enable convection in the domain, does this the flux automatically switch to j=-D(grad c)+uc or I will have to manually adjust for the convective flux (uc) in terms of expression of j_0?
↧