[Looby219]

The only way i can get the PEM module to function is by turning off reaction/joule heating and also liquid phase.

Activating reaction/joule heating causes the O2 residual to rise to 1e+2 and not converge. It also causes I_anode to be smaller than I_cathode even after many iterations.

Activating liquid phase causes Capillary Pressure residual to rise to 1e+10 and stay there.

I quadruple checked all boundary conditions. I'm only simulating a 3mm x 3mm cell so the following inlet conditions are used:

Air inlet: 2e-8 kg/s, 0.3 O2, 0.1 H20 mass fracts
H2 inlet: 4e-9 kg/s, 0.7 H2, 0.3 H20 mass fracts

Using V_cell = 0.5V, Temp is 353K, 1 ATM pressure, Mesh is certainly not an issue.