Collision.h

Go to the documentation of this file.

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

#include <memory>
#include <utility>
#include "MFSolver_CPU/CaseParameters.h"
#include "GridUtilities/LatticeParametersD3Q19.h"
#include "MFDatabase/ThreadArray.h"
#include "MFDatabase/Vec4.h"

namespace MF::Solver_CPU {

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

            inline static std::shared_ptr<MF::Solver_CPU::CaseParameters> m_CaseParameters_Ptr = nullptr;
            inline static std::shared_ptr<MF::GU::LatticeParametersD3Q19> m_LatticeParameters_Ptr = nullptr;

            static void Initialize (const std::shared_ptr<MF::Solver_CPU::CaseParameters>& CaseParameters_Ptr, const std::shared_ptr<MF::GU::LatticeParametersD3Q19>& LatticeParameters_Ptr) {
                m_CaseParameters_Ptr = CaseParameters_Ptr;
                m_LatticeParameters_Ptr = LatticeParameters_Ptr;
            };

            static void SetCollisionPointerToFunc(const std::shared_ptr<MF::Database::ThreadArray>& ThreadArray_Ptr); 

            // FEQ BGK getters used during initialization of computations ----------------------------------
            static double GetBGKFeq_Qcompr(const MF::Database::Vec4<double>& Vector4, const unsigned char& k); 
            static double GetBGKFeq_Incompr(const MF::Database::Vec4<double>& Vector4, const unsigned char& k); 

        private:
            // BGK -------------------------------------------------------------------------------
            static double Feq_Qcompr(const MF::Database::Vec4<double>& Vector4, const unsigned char& k); 
            static double Feq_Incompr(const MF::Database::Vec4<double>& Vector4, const unsigned char& k); 
            static void Coll_BGK_Qcompr(MF::Database::Node * pNode); 
            static void Coll_BGK_Incompr(MF::Database::Node * pNode); 

            // FBGK with external force F ---------------------------------------------------------
            static void Coll_FBGK_Qcompr(MF::Database::Node * pNode); 
            static void Coll_FBGK_Incompr(MF::Database::Node * pNode); 
            static double Feq_Qcompr_F(const MF::Database::Vec4<double>& Vector4, const unsigned char& k, const MF::Database::Vec3<double> & XYZForce_LB); 
            static double Feq_Incompr_F(const MF::Database::Vec4<double>& Vector4, const unsigned char& k, const MF::Database::Vec3<double> & XYZForce_LB); 

            // MRT and MRT2 -------------------------------------------------------------------------------
            static void Coll_MRT_Qcompr(MF::Database::Node * pNode); 
            static void Coll_MRT_Incompr(MF::Database::Node * pNode); 
            static void Coll_MRT2_Qcompr(MF::Database::Node * pNode); 
            static void Coll_MRT2_Incompr(MF::Database::Node * pNode); 
            static double Meq_Qcompr(const MF::Database::Vec4<double>& Vector4, const unsigned char& k); 
            static double Meq_Incompr(const MF::Database::Vec4<double>& Vector4, const unsigned char& k); 

            // Bounce-back 61. Do nothing since propagation function includes swap f(i) procedure -----------------------------------------
            static void Coll_BounceBackForNode61( MF::Database::Node * pNode) {}; 
        };
    } /* namespace MF */