ParametersConversion.cpp

Go to the documentation of this file.

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

// static parameters initialisation

double MF::LBPC::ParametersConversion::Re;
double MF::LBPC::ParametersConversion::Cu;
double MF::LBPC::ParametersConversion::u_ch_phys;
double MF::LBPC::ParametersConversion::l_ch_phys;
double MF::LBPC::ParametersConversion::nu_phys;
double MF::LBPC::ParametersConversion::dx_phys;
double MF::LBPC::ParametersConversion::l_ch_LB;
double MF::LBPC::ParametersConversion::u_ch_LB;
double MF::LBPC::ParametersConversion::dt_phys;
double MF::LBPC::ParametersConversion::nu_LB;
double MF::LBPC::ParametersConversion::TAU;
double MF::LBPC::ParametersConversion::csq;
double MF::LBPC::ParametersConversion::rho0_LB;
double MF::LBPC::ParametersConversion::rho0_phys;
int MF::LBPC::ParametersConversion::FluidFlowModel_MT;

void MF::LBPC::ParametersConversion::ComputeBasicParameters(const std::shared_ptr<MF::Solver_CPU::CaseParameters>& CaseParameters_Ptr,const openvdb::Int64Grid::Ptr& VDBGrid_Ptr) {

    openvdb::CoordBBox bbox;

    Re = u_ch_phys * l_ch_phys / nu_phys;
    l_ch_LB = l_ch_phys/dx_phys;
    Cu = 1.0 / ((TAU - 0.5) * csq) * (nu_phys / dx_phys);
    u_ch_LB = u_ch_phys / Cu;
    dt_phys = l_ch_phys / u_ch_phys * u_ch_LB / l_ch_LB;
    nu_LB = u_ch_LB * l_ch_LB / Re;

    CaseParameters_Ptr->ReynoldsNumber_Re = Re;
    CaseParameters_Ptr->CharacteristicVelocityLBu_U_Ch_LB = u_ch_LB;
    CaseParameters_Ptr->CharacteristicLength_L_CH_LB = l_ch_LB;
    CaseParameters_Ptr->DTPhysical_Dt_Phys = dt_phys;
    CaseParameters_Ptr->FluidLBuViscosity_Nu_LB = nu_LB;
    CaseParameters_Ptr->ForceX_LB = MF::LBPC::ParametersConversion::FPhys_toFLB(CaseParameters_Ptr->ForceX_fiz);
    CaseParameters_Ptr->ForceY_LB = MF::LBPC::ParametersConversion::FPhys_toFLB(CaseParameters_Ptr->ForceY_fiz);
    CaseParameters_Ptr->ForceZ_LB = MF::LBPC::ParametersConversion::FPhys_toFLB(CaseParameters_Ptr->ForceZ_fiz);
    CaseParameters_Ptr->XInitialVelocityLBu_Ux0_LB = MF::LBPC::ParametersConversion::UPhys_toLB(CaseParameters_Ptr->XInitialVelocityUx0_Phys);
    CaseParameters_Ptr->YInitialVelocityLBu_Uy0_LB = MF::LBPC::ParametersConversion::UPhys_toLB(CaseParameters_Ptr->YInitialVelocityUy0_Phys);
    CaseParameters_Ptr->ZInitialVelocityLBu_Uz0_LB = MF::LBPC::ParametersConversion::UPhys_toLB(CaseParameters_Ptr->ZInitialVelocityUz0_Phys);

    VDBGrid_Ptr->treePtr()->evalActiveVoxelBoundingBox(bbox);
    CaseParameters_Ptr->XNodesNumber_Nx1 = bbox.dim().x();
    CaseParameters_Ptr->YNodesNumber_Ny1 = bbox.dim().y();
    CaseParameters_Ptr->ZNodesNumber_Nz1 = bbox.dim().z();
}