I can confirm the reported behavior. VTK output #2 is configured as Particle data and RX1 rate is checked - this sets the keyword vtk_part_rrate(2,1) = .True.
But looking at the output file BACKGROUND_IC_1_0000.vtp
strings BACKGROUND_IC_1_0000.vtp |grep Name
<DataArray type="Float32" Name="coordinates" NumberOfComponents="3" format="appended" offset=" 0 " />
<DataArray type="Float32" Name="Diameter" format="appended" offset=" 3172" />
<DataArray type="Float32" Name="Velocity" NumberOfComponents="3" format="appended" offset=" 4232" />
<DataArray type="Float32" Name="Temperature" format="appended" offset=" 7404" />
<DataArray type="Float32" Name="X_s(1)" format="appended" offset=" 8464" />
the reaction rate is not there.
As an experiment I increased nrr - this shouldn’t matter, and it didn’t.
I also tried saving the reaction rate as cell data, as part of VTK output region #1 - this worked. This sets vtk_des_rrate(1,1) = .True. and RRates_RX appears in the VTU file:
strings BACKGROUND_IC_0000.vtu |grep Name
<DataArray type="Float32" Name="coordinates" NumberOfComponents="3" format="appended" offset=" 0" />
<DataArray type="Int32" Name="connectivity" format="appended" offset=" 976" />
<DataArray type="Int32" Name="offsets" format="appended" offset=" 2004" />
<DataArray type="Int32" Name="types" format="appended" offset=" 2136" />
<DataArray type="Float32" Name="EP_G" format="appended" offset=" 2268" />
<DataArray type="Float32" Name="P_G" format="appended" offset=" 2400" />
<DataArray type="Float32" Name="Gas_Velocity" NumberOfComponents="3" format="appended" offset=" 2532" />
<DataArray type="Float32" Name="Gas_temperature" format="appended" offset=" 2920" />
<DataArray type="Float32" Name="H2O_Gas_mass_fractions_1" format="appended" offset=" 3052" />
<DataArray type="Float32" Name="CH4_Gas_mass_fractions_2" format="appended" offset=" 3184" />
<DataArray type="Float32" Name="RRates_RX1" format="appended" offset=" 3316" />
So it looks like VTK rrates works for cell data but not particle data.
– Charles
