CaseParameters.cpp

Go to the documentation of this file.

// ==============================================================================================
// Microflow 3D, http://www.microflow.pwr.edu.pl/
// Created by Roman Szafran on 28.04.19.
// Copyright (c) 2019 Wroclaw University of Science and Technology.
// Distributed under the Apache License, Version 2.0. You may obtain a copy of the License at
// http://www.apache.org/licenses/LICENSE-2.0 or see accompanying file license.txt.
// Redistributions of source code must retain the above copyright and license notice.
// ==============================================================================================

#include "CaseParameters.h"

MF::Solver_CPU::CaseParameters::CaseParameters(const std::shared_ptr<MF::Database::ConfigData>& ConfigData_Ptr) {
    try {
        //Geometry parameters
        GeometryName = ConfigData_Ptr->getCaseStringParam("GeometryParams.GeometryName");
        VoxelSize = ConfigData_Ptr->getCaseFloatParam("GeometryParams.VoxelSize");
        PhysicalVoxelSize = ConfigData_Ptr->getCaseFloatParam("GeometryParams.PhysicalVoxelSize");
        AddInnerWall = ConfigData_Ptr->getCaseFloatParam("GeometryParams.AddInnerWall");

        //Parameters of model
        GeometryDefinition = (GeometryDefinitionClass) ConfigData_Ptr->getGeometryDefinition();
        FluidFlowModel_MT = (FlowTypeClass) ConfigData_Ptr->getFlowType();
        CollisionModel_KT = (ModelTypeClass) ConfigData_Ptr->getModelType();
        LatticeType_La = (LatticeTypeClass) ConfigData_Ptr->getLatticeType();

        //Physical parameters of simulation
        FluidPhysicalViscosity_Nu_Phys = ConfigData_Ptr->getCaseFloatParam("PhysicalParams.FluidPhysicalViscosity_Nu_Phys");
        FluidPhysicalReferenceDensity_Rho0_Phys = ConfigData_Ptr->getCaseFloatParam("PhysicalParams.FluidPhysicalReferenceDensity_Rho0_Phys");
        Tau = ConfigData_Ptr->getCaseFloatParam("PhysicalParams.Tau");
        ForceX_fiz = ConfigData_Ptr->getCaseFloatVectorParam(0, "PhysicalParams.XYZForce_fiz");
        ForceY_fiz = ConfigData_Ptr->getCaseFloatVectorParam(1, "PhysicalParams.XYZForce_fiz");
        ForceZ_fiz = ConfigData_Ptr->getCaseFloatVectorParam(2, "PhysicalParams.XYZForce_fiz");

        //Case characteristic parameters
        CharactLengthPhysical_L_CH_Phys = ConfigData_Ptr->getCaseFloatParam("CharactersticParams.CharactLengthPhysical_L_CH_Phys");
        CharactVelocityPhysical_U_CH_Phys = ConfigData_Ptr->getCaseFloatParam("CharactersticParams.CharactVelocityPhysical_U_CH_Phys");

        //Termination conditions
        Sc_VelocityResidueError_ErrV = ConfigData_Ptr->getCaseFloatParam("TerminationCondition.Sc_VelocityResidueError_ErrV");
        MassFlowError_ErrM = ConfigData_Ptr->getCaseFloatParam("TerminationCondition.MassFlowError_ErrM");

        //Initial and reference conditions
        ReferenceDensityLB_Rho0_LB = ConfigData_Ptr->getCaseFloatParam("InitialCondition.ReferenceDensityLB_Rho0_LB");

        XInitialVelocityUx0_Phys = ConfigData_Ptr->getCaseFloatVectorParam(0, "InitialCondition.XYZInitialVelocityPhysical_U0");
        YInitialVelocityUy0_Phys = ConfigData_Ptr->getCaseFloatVectorParam(1, "InitialCondition.XYZInitialVelocityPhysical_U0");
        ZInitialVelocityUz0_Phys = ConfigData_Ptr->getCaseFloatVectorParam(2, "InitialCondition.XYZInitialVelocityPhysical_U0");

        //Parameters of automatic detection and addition of boundary nodes
        DefaultWallNode_BN = ConfigData_Ptr->getCaseIntParam("DefaultBoundaryCondition.DefaultWallNode_BN");
        DefaultNode_NN = ConfigData_Ptr->getCaseIntParam("DefaultBoundaryCondition.DefaultNode_NN");

        // Global parameters calculated during initialization
        XNodesNumber_Nx1 = (0);
        YNodesNumber_Ny1 = (0);
        ZNodesNumber_Nz1 = (0);
        ReynoldsNumber_Re = (0.0);
        CharacteristicVelocityLBu_U_Ch_LB = (0.0);
        DTPhysical_Dt_Phys = (0.0);
        CharacteristicLength_L_CH_LB = (0.0);
        FluidLBuViscosity_Nu_LB = (0.0);
        ForceX_LB = (0.0);
        ForceY_LB = (0.0);
        ForceZ_LB = (0.0);
        XInitialVelocityLBu_Ux0_LB = (0.0);
        YInitialVelocityLBu_Uy0_LB = (0.0);
        ZInitialVelocityLBu_Uz0_LB = (0.0);

        //MFThread parameters
        ThreadMaxSize = ConfigData_Ptr->getCaseIntParam("ThreadParams.MFThreadMaxSize");

        //Calculation and display parameters
        CheckPointNumberRead_CHP_R = ConfigData_Ptr->getCaseIntParam("RWParams.CheckPointNumberRead_CHP_R");
        VTKWriteStep_VTK_W = ConfigData_Ptr->getCaseIntParam("RWParams.VTKWriteStep_VTK_W");
        VTKFileMaxNumber_VTK_Max = ConfigData_Ptr->getCaseIntParam("RWParams.VTKFileMaxNumber_VTK_Max");
        ConsoleWriteStep_K_W = ConfigData_Ptr->getCaseIntParam("RWParams.ConsoleWriteStep_K_W");

        //Mass balance surface choice params
        AutomaticSurfaceChoice = ConfigData_Ptr->getCaseBooleanParam("BalanceParams.AutomaticSurfaceChoice");

        // CPU parameters
        CPU_ThreadsNr = ConfigData_Ptr->getCaseIntParam("CPUParams.CPU_ThreadsNr");

    }
    catch (const config4cpp::ConfigurationException &exception) {
        std::cout << "[CaseParameters::Initialize] Error during initialization of case parameter from a file: ";
        std::cout << exception.c_str() << std::endl;
        exit(EXIT_FAILURE);
    }
}