I have been seeing an error when running simulations for a system containing two components (pH and salt) which are both treated as non-binding (isotherm is NoneBinding) where the inlet profile is given. The simulation is running these components through a PFR (lumped rate without pores) and column (GRM). The error I receive (see below) is not consistently reproducible - for example, ~2/10 simulations (with identical inputs) crash. It happens at the very start of the simulation when assembling the initial conditions. I am pretty confident all the code (using MATLAB for this) is correct because my simulations perform as expected when it doesn’t crash. Worth noting that after the simulation crashes, additional simulations will not run in that MATLAB window - maybe because the instance of CADET it is tied to gets stuck? So then I open a new window and try again, hoping for the best that it runs, which is inconvenient. I have tried changing the tolerances and other time integrator settings, but this didn’t help address the problem.
Error in simulation: Error in IDASolve: IDA_CONV_FAIL at t = 0.000000
Any tips on how to address this would be greatly appreciated.