Initialization.h

Go to the documentation of this file.

// ==============================================================================================
// Microflow 3D, http://www.microflow.pwr.edu.pl/
// Created by Roman Szafran on 03.06.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.
// ==============================================================================================

# pragma once

// OpenVDB
#include <openvdb/openvdb.h>

#include <memory>
#include <vector>

#include "MFDatabase/ThreadArray.h"
#include "GridUtilities/LatticeParametersD3Q19.h"
#include "GeometryBuilder/NodeID.h"
#include "MFDatabase/ConfigData.h"
#include "MFDatabase/Vec3.h"
#include "MFDatabase/Vec4.h"
#include "LBParametersConverter/ParametersConversion.h"
#include "MFSolver_CPU/CaseParameters.h"
#include "MFSolver_CPU/Collision.h"

namespace MF::Solver_CPU {

        class Initialization {
        public:
            Initialization () = default;
            ~Initialization() = default;

            static void ThreadDataInitialize (const std::shared_ptr<MF::Database::ThreadArray>& ThreadArray_Ptr, const std::shared_ptr<MF::Database::ConfigData>& ConfigData_Ptr);

            static void ThreadDefaultDataInitialize (const std::shared_ptr<MF::Database::ThreadArray>& ThreadArray_Ptr, const std::shared_ptr<MF::Solver_CPU::CaseParameters>& CaseParameters_Ptr);

            static void EquilibriumInitialization(const std::shared_ptr<MF::Database::ThreadArray>& ThreadArray_Ptr, const std::shared_ptr<MF::Solver_CPU::CaseParameters>& CaseParameters_Ptr){
                if (CaseParameters_Ptr->FluidFlowModel_MT == MF::Solver_CPU::FlowTypeClass::QuasiCompressible)
                    EquilibriumQcompr(ThreadArray_Ptr, CaseParameters_Ptr);
                else if (CaseParameters_Ptr->FluidFlowModel_MT == MF::Solver_CPU::FlowTypeClass::Incompressible)
                    EquilibriumIncompr(ThreadArray_Ptr, CaseParameters_Ptr);
                else {
                    std::cout << "Error in Initialization.h -> Unknown fluid flow model " << std::endl;
                    exit(EXIT_FAILURE);
                }
            };

        private:
            static void EquilibriumQcompr(const std::shared_ptr<MF::Database::ThreadArray>& ThreadArray_Ptr, const std::shared_ptr<MF::Solver_CPU::CaseParameters>& CaseParameters_Ptr);
            static void EquilibriumIncompr(const std::shared_ptr<MF::Database::ThreadArray>& ThreadArray_Ptr, const std::shared_ptr<MF::Solver_CPU::CaseParameters>& CaseParameters_Ptr);

        };
} /* namespace MF::Solver_CPU */