Hi,
I couldn’t guess why you were not able to measure any voltage in what’s meant to be a similar setup to what’s given in the ion_burner_flame.py example. This example is fairly similar to the one shown in the work of Han et al (2015), where there are plots showing H3O+ and E- concentrations with comparison to some experimental measurements.
In older Cantera versions, the “stage=2” solver option is the one in which the electric field calculation was enabled, and diffusion of ionized species was influenced by this field, while the “stage=1” solution neglected diffusion of ionized species. My understanding of the reason for taking this approach in the first stage is that if the electric field is neglected, diffusion of electrons would be grossly overestimated. In Cantera 3.2.0, these “stages” are now replaced by the electric_field_enabled property, which can be toggled in a manner similar to other solver options like the soret_enabled flag.
I’m not sure why the “frozen ion” method would ever be better at describing any physical phenomena. My interpretation was that it was just a feature to help with numerical convergence in some cases, similar to how we allow solving with the temperature profile fixed first followed by allowing the temperature to evolve.
There is currently no way to set boundary condition for the electric field (i.e., the electric field at the burner surface), but I don’t think that would be terribly difficult to implement.
Regards,
Ray