BoundaryConditions

||Boundary Condition||OLD||NEW|| ||$\phi|_{x=L} = \phi_L$||`FixedValue(L, phiL)`||`phi.contrain(phiL, L)`|| ||$\hat{n}\cdot\nabla\phi|_{x=L} = \phi'_L$||~~no way~~||`phi.getFaceGrad().contrain(gradPhiL * mesh.faceNormals, L)`|| ||$\hat{n}\cdot\vec{J}_\phi|_{x=L} = 0$||default|| default (in absence of constraints)|| ||$\hat{n}\cdot\vec{J}_\phi|_{x=L} = J_L$|| `FixedFlux(L, JL)` ||$\hat{n}\cdot\vec{v}|_{x=L} = 0$ && `JL.getDivergence()`|| ||$\hat{n}\cdot\nabla^2\vec{J}_\phi|_{x=L} = J_L$||`NthOrderBoundaryCondition()`||`NthOrderBoundaryCondition()`|| ||(in/out)flow||`Implicit/ExplicitSourceTerm()` (a mess)|| $\hat{n}\cdot\vec{v}|_{x=L}$ determined by constraints ($\hat{n}\cdot\vec{v}|_{x=L} = 0$ in absence of constraints) ||