[cig-commits] commit 1865 by buerg to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 1865

Change formula for temperature source term.

U   trunk/aspire/source/material_model/cantera.cc
U   trunk/aspire/source/simulator/assembly.cc


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=1865&peg=1865

Diff:
Modified: trunk/aspire/source/material_model/cantera.cc
===================================================================
--- trunk/aspire/source/material_model/cantera.cc	2013-08-30 16:47:08 UTC (rev 1864)
+++ trunk/aspire/source/material_model/cantera.cc	2013-08-30 16:47:55 UTC (rev 1865)
@@ -60,10 +60,10 @@
           for (unsigned int i = 0; i < n_species; ++i)
             {
               out.compositional_sources[q][i] = molar_mass[i] * net[i];
-              out.thermal_sources[q] -= in.temperature[q] * enthalpy_RT[i] * out.compositional_sources[q][i];
+              out.thermal_sources[q] += enthalpy_RT[i] * out.compositional_sources[q][i];
             }
           
-          out.thermal_sources[q] *= ::Cantera::GasConstant / out.specific_heat[q];
+          out.thermal_sources[q] *= -1.0;
           out.thermal_expansion_coefficients[q] = ideal_gas_mix->thermalExpansionCoeff ();
           out.viscosities[q] = transport->viscosity ();
         }

Modified: trunk/aspire/source/simulator/assembly.cc
===================================================================
--- trunk/aspire/source/simulator/assembly.cc	2013-08-30 16:47:08 UTC (rev 1864)
+++ trunk/aspire/source/simulator/assembly.cc	2013-08-30 16:47:55 UTC (rev 1865)
@@ -1027,7 +1027,7 @@
         const double lambda = scratch.material_model_outputs.thermal_conductivities[q];
         const double old_value = scratch.old_temperature_values[q];
         const double rho = scratch.material_model_outputs.densities[q];
-        const double source_term = scratch.material_model_outputs.thermal_sources[q];
+        const double source_term = scratch.material_model_outputs.thermal_sources[q] / c_P;
         const Tensor<1, dim> current_u = scratch.current_velocity_values[q];
 
         for (unsigned int i = 0; i < dofs_per_cell; ++i)