The problem in outputting reaction rates

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