Hi Jeff. Thank you for the clarification (I realized this shortly before your reply and have implemented the necessary changes). After defining RO_xs0 for all relevant species, however, the declared variable didn’t return the appropriate value when called - this might be something on my end, but for the time being, I defined the densities as parameters within a UDF file. The other error I was running into was likely due to a Boolean condition that I had commented out (i.e., there was a negative field variable or a close equivalent). Upon implementing these changes, I got a partially working model!
However, one of the rates (specifically for the defined “Decompf” reaction) was on the order of 10^26 (this might be attributed to a Boolean condition I invoke to ensure a I don’t run into a divide by zero error, so I may need to increase the threshold value if that’s the case). However, I’m now running into an error that likely has to do with variable declarations within my reaction rates. Specifically, I’m getting the error “float divide by zero in __leq_gmres_mod_MOD_leq_gmres0” around 30 seconds after starting the simulation. This might be attributed to a variable declaration within the reaction rates (i.e., a negative field variable). Any insight?
Here’s an updated “.zip” file with all recorded updates:
clc_main_2024-08-13T142353.947723.zip (1.4 MB)