new temperature limiter

src/GksGpu/FluxComputation/FluxComputation.cu

@@ -234,20 +234,20 @@ __host__ __device__ inline void fluxFunction(DataBaseStruct dataBase, Parameters
 
             ////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
             // this one works for some time
-            real S = parameters.dx * parameters.dx * ( fabsf(dTdx1) + fabsf(dTdx2) + fabsf(dTdx3) );
+            //real S = parameters.dx * parameters.dx * ( fabsf(dTdx1) + fabsf(dTdx2) + fabsf(dTdx3) );
             //k += real(0.00002) / real(0.015625) * S;
 
             //real T = getT(facePrim);
             //if( T > 20 )
-                k += parameters.temperatureLimiter * S;
+                //k += parameters.temperatureLimiter * S;
 
             ////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
             
             //real S = parameters.dx * ( fabsf(dTdx1) + fabsf(dTdx2) + fabsf(dTdx3) );
             //k += real(0.00001) * real(0.0025) * S * S;
             
-            //real S = parameters.dx * parameters.dx * ( dTdx1 * dTdx1 + dTdx2 * dTdx2 + dTdx3 * dTdx3 );
-            //k += real(0.00001) * real(0.001) * S;
+            real S = parameters.dx * parameters.dx * ( dTdx1 * dTdx1 + dTdx2 * dTdx2 + dTdx3 * dTdx3 );
+            k += fminf(real(0.001), parameters.temperatureLimiter * S);
 
             ////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
             // this one works for some time

targets/apps/GKS/Flame7cm/Flame7cm.cpp

@@ -139,13 +139,13 @@ void thermalCavity( std::string path, std::string simulationName, uint _gpuIndex
 
     parameters.useHeatReleaseRateLimiter = true;
     parameters.useReactionLimiter        = true;
-    parameters.useTemperatureLimiter     = false;
+    parameters.useTemperatureLimiter     = true;
     parameters.usePassiveScalarLimiter   = true;
     parameters.useSmagorinsky            = true;
 
     parameters.heatReleaseRateLimiter = 5000000.0;
     parameters.reactionLimiter        = 1.005;
-    parameters.temperatureLimiter     = 1.0e-3;
+    parameters.temperatureLimiter     = 1.0e-8;//1.0e-3;
 
     parameters.useSpongeLayer = true;
     parameters.spongeLayerIdx = 0;