Follow-up from "[PL/CT] Molar flux and molar flow node-wise output"
The following discussion from !4006 (merged) should be addressed:
-
@chleh started a discussion: (+1 comment)
I'm fine for this MR. Two questions for future developments:
The current MR enables molar flow rate only for the Newton scheme (or did I overlook something?). Should we also have it for the Picard scheme?
Looking into the local assembler the molar flow rate will contain the following contributions:
local_rhs.noalias() -= w * rho * (N.transpose() * phi * R * N * (cdot + alpha * c) + dNdx.transpose() * D * dNdx * c); local_rhs.noalias() -= w * rho * N.transpose() * q.transpose() * dNdx * c;
One contribution is related to the time derivative
cdot
. Should these individual contributions be split in the output?