Skip to content

GitLab

Closed
Created by Kilian Weishaupt@kweisOwner0 of 1 task completed0/1 task

[md][stokesdarcy] Upwinding of component enthalpy is probably wrong

const bool insideDiffFluxIsUpstream = std::signbit(diffusiveFlux[domainICompIdx]);

should probably be

const bool insideDiffFluxIsUpstream = diffusiveFlux[domainICompIdx] > 0.0;
  • Furthermore, we should check if all fluidsystems implement componentEnthalpy.
Edited by Kilian Weishaupt