python interface updates for radiation models

Author: wandadars

interfaces/cython/cantera/_onedim.pxd

@@ -80,6 +80,7 @@ cdef extern from "cantera/oneD/Boundary1D.h":
 cdef extern from "cantera/oneD/Flow1D.h":
     cdef cppclass CxxFlow1D "Cantera::Flow1D" (CxxDomain1D):
         void setTransportModel(const string&) except +translate_exception
+        void setRadiationModels(const string&, const string&) except +translate_exception
         string type()
         string transportModel()
         void setPressure(double)

interfaces/cython/cantera/_onedim.pyx

@@ -790,6 +790,20 @@ cdef class FlowBase(Domain1D):
             y.push_back(t)
         self.flow.setFixedTempProfile(x, y)
 
+    def set_radiation_models(self, property_model, solver_model="OpticallyThin"):
+        """
+        Set radiation property and solver models for this flow domain.
+
+        Parameters
+        ----------
+        property_model : str
+            Name of the property calculator model. Examples include
+            "TabularPlanckMean", "RadLib.PlanckMean", "RadLib.WSGG", "RadLib.RCSLW".
+        solver_model : str, optional
+            Name of the radiation solver model. Defaults to "OpticallyThin".
+        """
+        self.flow.setRadiationModels(stringify(property_model), stringify(solver_model))
+
     def get_settings3(self):
         """
         Temporary method returning new behavior of settings getter.

interfaces/cython/cantera/onedim.py

@@ -258,6 +258,17 @@ def boundary_emissivities(self, epsilon):
             raise ValueError("Boundary emissivities must both be set at the same time.")
         self.flame.boundary_emissivities = epsilon[0], epsilon[1]
 
+    def set_radiation_models(self, property_model, solver_model="OpticallyThin"):
+        """
+        Set the radiation property and solver models.
+
+        Examples
+        --------
+        >>> f.set_radiation_models("TabularPlanckMean", "OpticallyThin")
+        >>> f.set_radiation_models("RadLib.WSGG")  # solver defaults to OpticallyThin
+        """
+        self.flame.set_radiation_models(property_model, solver_model)
+
     @property
     def grid(self):
         """ Array of grid point positions along the flame. """