new target for Non-Overbeck-Boussinesq convection

src/GksGpu/Analyzer/HeatFluxAnalyzer.cu

+#include "HeatFluxAnalyzer.h"
+
+#include <cuda.h>
+#include <cuda_runtime.h>
+#include <helper_cuda.h>
+
+#include <cmath>
+#include <sstream>
+
+#include <thrust/device_vector.h>
+#include <thrust/reduce.h>
+#include <thrust/device_ptr.h>
+
+#include <iomanip>
+
+#include "Core/Logger/Logger.h"
+
+#include "DataBase/DataBase.h"
+
+#include "GksGpu/BoundaryConditions/BoundaryCondition.h"
+
+#include "FlowStateData/AccessDeviceData.cuh"
+#include "FlowStateData/FlowStateDataConversion.cuh"
+
+#include "FluxComputation/SutherlandsLaw.cuh"
+
+#include "CudaUtility/CudaRunKernel.hpp"
+
+__global__                 void heatFluxKernel  ( DataBaseStruct  dataBase, GksGpu::BoundaryConditionStruct  boundaryCondition, Parameters  parameters, real* heatFlux, uint startIndex, uint numberOfEntities );
+__host__ __device__ inline void heatFluxFunction( DataBaseStruct& dataBase, GksGpu::BoundaryConditionStruct& boundaryCondition, Parameters& parameters, real* heatFlux, uint startIndex, uint index );
+
+////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
+
+bool HeatFluxAnalyzer::run(uint iter, Parameters parameters)
+{
+    if( iter % this->analyzeIter != 0 ) return false;
+
+    uint numberOfCells = this->boundaryCondition->numberOfCellsPerLevel[ dataBase->numberOfLevels - 1 ];
+
+    thrust::device_vector<real> heatFlux( numberOfCells );
+
+    CudaUtility::CudaGrid grid( numberOfCells, 32 );
+
+    for( uint level = 0; level < dataBase->numberOfLevels - 1; level++ ) parameters.dx *= c1o2; 
+
+    runKernel( heatFluxKernel,
+               heatFluxFunction,
+               dataBase->getDeviceType(), grid, 
+               dataBase->toStruct(),
+               boundaryCondition->toStruct(),
+               parameters,
+               heatFlux.data().get(),
+               boundaryCondition->startOfCellsPerLevel[ dataBase->numberOfLevels - 1 ] );
+
+    getLastCudaError("HeatFluxAnalyzer::run(uint iter)");
+
+    real q = thrust::reduce( heatFlux.begin(), heatFlux.end(), c0o1, thrust::plus<real>() ) * parameters.dx * parameters.dx;
+
+    real qIdeal = c1o4 * (parameters.K + c5o1) * ( parameters.mu / parameters.Pr ) * ( c1o1 / lambdaHot - c1o1 / lambdaCold );
+
+    this->heatFluxTimeSeries.push_back( q / qIdeal );
+
+    if( iter % this->outputIter == 0 ) *logging::out << logging::Logger::INFO_HIGH << "q = " << q / qIdeal << "\n";
+}
+
+//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
+
+__global__ void heatFluxKernel(DataBaseStruct  dataBase, GksGpu::BoundaryConditionStruct  boundaryCondition, Parameters  parameters, real* heatFlux, uint startIndex, uint numberOfEntities)
+{
+    uint index = blockIdx.x * blockDim.x + threadIdx.x;
+
+    if( index >= numberOfEntities ) return;
+
+    heatFluxFunction( dataBase, boundaryCondition, parameters, heatFlux, startIndex, index );
+}
+
+__host__ __device__ void heatFluxFunction(DataBaseStruct& dataBase, GksGpu::BoundaryConditionStruct& boundaryCondition, Parameters& parameters, real* heatFlux, uint startIndex, uint index)
+{
+    uint ghostCellIndex  = boundaryCondition.ghostCells [ startIndex + index ];
+    uint domainCellIndex = boundaryCondition.domainCells[ startIndex + index ];
+
+    if( isCellProperties( dataBase.cellProperties[ domainCellIndex ], CELL_PROPERTIES_GHOST ) )
+    {
+        heatFlux[ startIndex + index ] = c0o1;
+        return;
+    }
+
+    //////////////////////////////////////////////////////////////////////////
+
+    ConservedVariables ghostCons;
+
+    readCellData(ghostCellIndex, dataBase, ghostCons);
+
+    ConservedVariables domainCons;
+
+    readCellData(domainCellIndex, dataBase, domainCons);
+
+    PrimitiveVariables ghostPrim  = toPrimitiveVariables(ghostCons,  parameters.K);
+    PrimitiveVariables domainPrim = toPrimitiveVariables(domainCons, parameters.K);
+
+    //////////////////////////////////////////////////////////////////////////
+
+    real lambda = c1o2 * (ghostPrim.lambda + domainPrim.lambda);
+
+    real r   = parameters.lambdaRef / lambda;
+
+    real mu;
+    if ( parameters.viscosityModel == ViscosityModel::constant ){
+        mu = parameters.mu;
+    }
+    else if ( parameters.viscosityModel == ViscosityModel::sutherlandsLaw ){
+        mu = sutherlandsLaw( parameters, r );
+    }
+
+    heatFlux[ startIndex + index ] = c1o4 * (parameters.K + c5o1) * ( mu / parameters.Pr ) / parameters.dx * ( c1o1 / domainPrim.lambda - c1o1 / ghostPrim.lambda );
+}
+
+//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
+
+HeatFluxAnalyzer::HeatFluxAnalyzer( SPtr<DataBase> dataBase, SPtr<GksGpu::BoundaryCondition> boundaryCondition, uint analyzeIter, uint outputIter, real lambdaHot, real lambdaCold, real L )
+    : dataBase(dataBase), 
+      boundaryCondition(boundaryCondition), 
+      analyzeIter(analyzeIter), 
+      outputIter(outputIter), 
+      lambdaHot(lambdaHot), 
+      lambdaCold(lambdaCold), 
+      L(L)
+{
+}
+
+void HeatFluxAnalyzer::writeToFile(std::string filename)
+{
+    *logging::out << logging::Logger::INFO_INTERMEDIATE << "HeatFluxAnalyzer::writeToFile( " << filename << " )" << "\n";
+
+    std::ofstream file;
+
+    file.open(filename + ".dat" );
+
+    for( auto& EKin : this->heatFluxTimeSeries )
+        file << std::setprecision(15) << EKin << std::endl;
+
+    file.close();
+
+    *logging::out << logging::Logger::INFO_INTERMEDIATE << "done!\n";
+}
+
+

src/GksGpu/Analyzer/HeatFluxAnalyzer.h

+#ifndef  HeatFluxAnalyzer_H
+#define  HeatFluxAnalyzer_H
+
+#include <vector>
+#include <string>
+
+#include "VirtualFluidsDefinitions.h"
+
+#include "Core/PointerDefinitions.h"
+#include "Core/DataTypes.h"
+
+#include "GksGpu/BoundaryConditions/BoundaryCondition.h"
+
+#include "FlowStateData/FlowStateData.cuh"
+
+#include "Parameters/Parameters.h"
+
+struct DataBase;
+
+class VF_PUBLIC HeatFluxAnalyzer
+{
+private:
+
+    SPtr<DataBase> dataBase;
+    SPtr<GksGpu::BoundaryCondition> boundaryCondition;
+
+    uint outputIter;
+
+    uint analyzeIter;
+
+    std::vector<real> heatFluxTimeSeries;
+
+    real lambdaHot;
+    real lambdaCold;
+
+    real L;
+
+public:
+
+    HeatFluxAnalyzer( SPtr<DataBase> dataBase, SPtr<GksGpu::BoundaryCondition> boundaryCondition, uint analyzeIter, uint outputIter, real lambdaHot, real lambdaCold, real L );
+
+    bool run( uint iter, Parameters parameters );
+
+    void writeToFile( std::string filename );
+
+};
+
+#endif

targets/apps/GKS/RayleighBenardMultiGPU/3rdPartyLinking.cmake

+include (${CMAKE_SOURCE_DIR}/${cmakeMacroPath}/MPI/Link.cmake)
+linkMPI(${targetName})
+include (${CMAKE_SOURCE_DIR}/${cmakeMacroPath}/Cuda/Link.cmake)
+linkCuda(${targetName})
+#include (${CMAKE_SOURCE_DIR}/${cmakeMacroPath}/Metis/Link.cmake)
+#linkMetis(${targetName})
+
+#if(HULC.BUILD_JSONCPP)
+#  include (${CMAKE_SOUR#CE_DIR}/${cmakeMacroPath}/JsonCpp/Link.cmake)
+#  linkJsonCpp(${targetName})
+#endif()

targets/apps/GKS/RayleighBenardMultiGPU/CMakeLists.txt

+setTargetNameToFolderName(${CMAKE_CURRENT_LIST_DIR})
+
+set(linkDirectories "")
+set(libsToLink Core GridGenerator GksMeshAdapter GksVtkAdapter GksGpu)
+set(includeDirectories "${CMAKE_SOURCE_DIR}/src"
+                       "${CMAKE_SOURCE_DIR}/src/Core"
+                       "${CMAKE_SOURCE_DIR}/src/GridGenerator"
+                       "${CMAKE_SOURCE_DIR}/src/GksMeshAdapter"
+                       "${CMAKE_SOURCE_DIR}/src/GksVtkAdapter"
+                       "${CMAKE_SOURCE_DIR}/src/GksGpu")
+
+#glob files and save in MY_SRCS
+include(CMakePackage.cmake)
+
+buildExe(${targetName} "${MY_SRCS}" "${linkDirectories}" "${libsToLink}" "${includeDirectories}")
+groupTarget(${targetName} ${gksAppFolder})
+
+# Specify the linking to 3rdParty libs
+include(3rdPartyLinking.cmake)

targets/apps/GKS/RayleighBenardMultiGPU/CMakePackage.cmake

+#FILE ENDINGS
+resetFileEndingsToCollect()
+addCAndCPPFileTypes()
+addFileEndingToCollect("*.cu")
+addFileEndingToCollect("*.cuh")
+
+#GLOB SOURCE FILES IN MY_SRCS
+unset(MY_SRCS)
+includeRecursiveAllFilesFrom(${targetName} ${CMAKE_CURRENT_LIST_DIR})
\ No newline at end of file