175 lines
3.1 KiB
C++
175 lines
3.1 KiB
C++
/*------------------------------- phasicFlow ---------------------------------
|
|
O C enter of
|
|
O O E ngineering and
|
|
O O M ultiscale modeling of
|
|
OOOOOOO F luid flow
|
|
------------------------------------------------------------------------------
|
|
Copyright (C): www.cemf.ir
|
|
email: hamid.r.norouzi AT gmail.com
|
|
------------------------------------------------------------------------------
|
|
Licence:
|
|
This file is part of phasicFlow code. It is a free software for simulating
|
|
granular and multiphase flows. You can redistribute it and/or modify it under
|
|
the terms of GNU General Public License v3 or any other later versions.
|
|
|
|
phasicFlow is distributed to help others in their research in the field of
|
|
granular and multiphase flows, but WITHOUT ANY WARRANTY; without even the
|
|
implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
|
|
|
|
-----------------------------------------------------------------------------*/
|
|
|
|
#ifndef __rectangleMesh_hpp__
|
|
#define __rectangleMesh_hpp__
|
|
|
|
#include "types.hpp"
|
|
#include "error.hpp"
|
|
#include "box.hpp"
|
|
#include "IOobject.hpp"
|
|
|
|
class repository;
|
|
|
|
namespace pFlow
|
|
{
|
|
|
|
class rectangleMesh
|
|
:
|
|
public IOobject
|
|
{
|
|
private:
|
|
|
|
box meshBox_;
|
|
|
|
int32x3 numCells_;
|
|
|
|
realx3 dx_;
|
|
|
|
public:
|
|
|
|
TypeInfo("rectangleMesh");
|
|
|
|
rectangleMesh(
|
|
const box& mshBox,
|
|
int32 nx,
|
|
int32 ny,
|
|
int32 nz,
|
|
repository* rep);
|
|
|
|
rectangleMesh(const dictionary & dict, repository* rep);
|
|
|
|
|
|
rectangleMesh(const rectangleMesh&) = default;
|
|
|
|
|
|
rectangleMesh& operator = (const rectangleMesh&) = default;
|
|
|
|
rectangleMesh(rectangleMesh&&) = default;
|
|
|
|
|
|
rectangleMesh& operator = (rectangleMesh&&) = default;
|
|
|
|
~rectangleMesh() override = default;
|
|
|
|
int64 size()const
|
|
{
|
|
return numCells_.x_*numCells_.y_*numCells_.z_;
|
|
}
|
|
|
|
int32 nx()const
|
|
{
|
|
return numCells_.x();
|
|
}
|
|
|
|
int32 ny()const
|
|
{
|
|
return numCells_.y();
|
|
}
|
|
|
|
int32 nz()const
|
|
{
|
|
return numCells_.z();
|
|
}
|
|
|
|
real cellVol()const
|
|
{
|
|
return dx_.x_*dx_.y_*dx_.z_;
|
|
}
|
|
|
|
realx3 minPoint()const
|
|
{
|
|
return meshBox_.minPoint();
|
|
}
|
|
|
|
realx3 maxPoint()const
|
|
{
|
|
return meshBox_.maxPoint();
|
|
}
|
|
|
|
inline
|
|
bool isInsideIndex(const realx3 p, int32x3 & ind )const
|
|
{
|
|
if(meshBox_.isInside(p))
|
|
{
|
|
ind = (p - meshBox_.minPoint())/dx_ ;
|
|
return true;
|
|
}
|
|
else
|
|
{
|
|
return false;
|
|
}
|
|
}
|
|
|
|
bool write(iOstream& is, const IOPattern& iop)const override
|
|
{
|
|
notImplementedFunction;
|
|
return true;
|
|
}
|
|
|
|
bool read(iIstream& is, const IOPattern& iop) override
|
|
{
|
|
notImplementedFunction;
|
|
return true;
|
|
}
|
|
|
|
bool writeToVtk(iOstream& os)const
|
|
{
|
|
os<<"DATASET RECTILINEAR_GRID"<<endl;
|
|
os<<"DIMENSIONS "<<nx()+1<<" "<< ny()+1 << " "<< nz()+1 <<endl;
|
|
|
|
auto [x, y , z] = this->minPoint();
|
|
auto [dx, dy, dz] = dx_;
|
|
|
|
os<<"X_COORDINATES "<< nx()+1 <<" float\n";
|
|
for(int32 i=0; i<nx()+1; i++)
|
|
{
|
|
os<< x<<"\n";
|
|
x+= dx;
|
|
}
|
|
|
|
os<<"Y_COORDINATES "<< ny()+1 <<" float\n";
|
|
for(int32 j=0; j<ny()+1; j++)
|
|
{
|
|
os<< y <<"\n";
|
|
y+= dy;
|
|
}
|
|
|
|
os<<"Z_COORDINATES "<< nz()+1 <<" float\n";
|
|
for(int32 j=0; j<nz()+1; j++)
|
|
{
|
|
os<< z <<"\n";
|
|
z+= dz;
|
|
}
|
|
|
|
os<<"CELL_DATA "<< nx()*ny()*nz()<<endl;
|
|
|
|
return true;
|
|
}
|
|
|
|
};
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
#endif // __rectangleMesh_hpp__
|