I have a PDE system with:
- 3 PDEs
- A number of additional algebraic equations
- I have boundary and initial value conditions with a number of added algebraic equations to describe the boundary conditions
When I try to discretize the system, I get an error message on boundary conditions:
AssertionError: Boundary condition ṁ_a2t(t) ~ Y_a2t(t)*sqrt((max(Δp_a2t(t), 0)*ρ_a(t)) / (p_at2__ς*ρ_a2t__ς))*f_v(u_a2t(t))*ṁ_a2t__c is not on a boundary of the domain, or is not a valid boundary condition
Here is the relevant model of the intended boundary condition:
Does the discretizer have a problem with inserting the expression for pressure difference, etc.?