superboundarycondition3d.hh

来自「open lattice boltzmann project www.open」· HH 代码 · 共 1,230 行 · 第 1/4 页

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addPressureBoundary2P( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge0NN(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge0NN( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge0NP(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge0NP( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge0PN(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge0PN( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge0PP(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge0PP( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge1NN(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge1NN( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge1NP(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge1NP( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge1PN(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge1PN( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge1PP(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge1PP( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge2NN(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge2NN( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge2NP(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge2NP( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge2PN(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge2PN( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addExternalVelocityEdge2PP(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addExternalVelocityEdge2PP( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addInternalVelocityEdge0NN(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addInternalVelocityEdge0NN( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }
    }
}

template<typename T, template<typename U> class Lattice>
void sOnLatticeBoundaryCondition3D<T,Lattice>::
        addInternalVelocityEdge0NP(T x0, T x1, T y0, T y1, T z0, T z1, T omega){

    int overlap = _sLattice.get_overlap();
    int nC = _sLattice.get_load().size();
    int locX0, locX1, locY0, locY1, locZ0, locZ1;
    for (int iC=0; iC<nC; iC++) {
        int iCglob = _sLattice.get_load().glob(iC);
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters( x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1, overlap)) {
            _blockBCs[iC]->addInternalVelocityEdge0NP( 
                locX0, locX1, locY0, locY1, locZ0, locZ1, omega);
        }
        if (_sLattice.get_cGeometry().get_cuboid(iCglob).checkInters(x0, x1,
            y0, y1, z0, z1, locX0, locX1, locY0, locY1, locZ0, locZ1)) {
            addPoints2CommBC(locX0, locX1, locY0, locY1, locZ0, locZ1, iCglob);
        }