cuboidWall
This commit is contained in:
hamidrezanorouzi
parent
02f7f9af20
commit
a9097ad286
|
|
@ -5,6 +5,7 @@ Wall/Wall.C
|
||||||
planeWall/planeWall.C
|
planeWall/planeWall.C
|
||||||
stlWall/stlWall.C
|
stlWall/stlWall.C
|
||||||
cylinderWall/cylinderWall.C
|
cylinderWall/cylinderWall.C
|
||||||
|
cuboidWall/cuboidWall.C
|
||||||
)
|
)
|
||||||
set(link_lib phasicFlow Geometry Kokkos::kokkos)
|
set(link_lib phasicFlow Geometry Kokkos::kokkos)
|
||||||
|
|
||||||
|
|
|
||||||
|
|
@ -0,0 +1,136 @@
|
||||||
|
/*------------------------------- 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.
|
||||||
|
|
||||||
|
-----------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
|
||||||
|
#include "cuboidWall.H"
|
||||||
|
#include "planeWall.H"
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
bool pFlow::cuboidWall::readcuboidWall
|
||||||
|
(
|
||||||
|
const dictionary& dict
|
||||||
|
)
|
||||||
|
{
|
||||||
|
|
||||||
|
auto center = dict.getVal<realx3>("center");
|
||||||
|
auto edgeLength= dict.getVal<realx3>("edgeLength");
|
||||||
|
|
||||||
|
realx3 p1,p2,p3,p4;
|
||||||
|
|
||||||
|
// left plane
|
||||||
|
p1 = center + edgeLength*realx3(-0.5,-0.5,-0.5);
|
||||||
|
p2 = center + edgeLength*realx3(-0.5, 0.5,-0.5);
|
||||||
|
p3 = center + edgeLength*realx3(-0.5, 0.5, 0.5);
|
||||||
|
p4 = center + edgeLength*realx3(-0.5,-0.5, 0.5);
|
||||||
|
|
||||||
|
planeWall left(p1,p2,p3,p4);
|
||||||
|
|
||||||
|
for(const auto& t:left.triangles())
|
||||||
|
{
|
||||||
|
triangles_.push_back(t);
|
||||||
|
}
|
||||||
|
|
||||||
|
// right plane
|
||||||
|
p1 = center + edgeLength*realx3( 0.5,-0.5,-0.5);
|
||||||
|
p2 = center + edgeLength*realx3( 0.5,-0.5, 0.5);
|
||||||
|
p3 = center + edgeLength*realx3( 0.5, 0.5, 0.5);
|
||||||
|
p4 = center + edgeLength*realx3( 0.5, 0.5,-0.5);
|
||||||
|
|
||||||
|
planeWall right(p1,p2,p3,p4);
|
||||||
|
|
||||||
|
for(const auto& t:right.triangles())
|
||||||
|
{
|
||||||
|
triangles_.push_back(t);
|
||||||
|
}
|
||||||
|
|
||||||
|
// bottom plane
|
||||||
|
p1 = center + edgeLength*realx3(-0.5,-0.5,-0.5);
|
||||||
|
p2 = center + edgeLength*realx3(-0.5,-0.5, 0.5);
|
||||||
|
p3 = center + edgeLength*realx3( 0.5,-0.5, 0.5);
|
||||||
|
p4 = center + edgeLength*realx3( 0.5,-0.5,-0.5);
|
||||||
|
|
||||||
|
planeWall bottom(p1,p2,p3,p4);
|
||||||
|
|
||||||
|
for(const auto& t:bottom.triangles())
|
||||||
|
{
|
||||||
|
triangles_.push_back(t);
|
||||||
|
}
|
||||||
|
|
||||||
|
// top plane
|
||||||
|
p1 = center + edgeLength*realx3(-0.5, 0.5,-0.5);
|
||||||
|
p2 = center + edgeLength*realx3( 0.5, 0.5,-0.5);
|
||||||
|
p3 = center + edgeLength*realx3( 0.5, 0.5, 0.5);
|
||||||
|
p4 = center + edgeLength*realx3(-0.5, 0.5, 0.5);
|
||||||
|
|
||||||
|
planeWall top(p1,p2,p3,p4);
|
||||||
|
|
||||||
|
for(const auto& t:top.triangles())
|
||||||
|
{
|
||||||
|
triangles_.push_back(t);
|
||||||
|
}
|
||||||
|
|
||||||
|
// back plane
|
||||||
|
p1 = center + edgeLength*realx3(-0.5,-0.5,-0.5);
|
||||||
|
p2 = center + edgeLength*realx3( 0.5,-0.5,-0.5);
|
||||||
|
p3 = center + edgeLength*realx3( 0.5, 0.5,-0.5);
|
||||||
|
p4 = center + edgeLength*realx3(-0.5, 0.5,-0.5);
|
||||||
|
|
||||||
|
planeWall back(p1,p2,p3,p4);
|
||||||
|
|
||||||
|
for(const auto& t:back.triangles())
|
||||||
|
{
|
||||||
|
triangles_.push_back(t);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// fron plane
|
||||||
|
p1 = center + edgeLength*realx3(-0.5,-0.5, 0.5);
|
||||||
|
p2 = center + edgeLength*realx3(-0.5, 0.5, 0.5);
|
||||||
|
p3 = center + edgeLength*realx3( 0.5, 0.5, 0.5);
|
||||||
|
p4 = center + edgeLength*realx3( 0.5,-0.5, 0.5);
|
||||||
|
|
||||||
|
planeWall front(p1,p2,p3,p4);
|
||||||
|
|
||||||
|
for(const auto& t:front.triangles())
|
||||||
|
{
|
||||||
|
triangles_.push_back(t);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
return true;
|
||||||
|
|
||||||
|
}
|
||||||
|
|
||||||
|
pFlow::cuboidWall::cuboidWall()
|
||||||
|
{}
|
||||||
|
|
||||||
|
pFlow::cuboidWall::cuboidWall
|
||||||
|
(
|
||||||
|
const dictionary& dict
|
||||||
|
)
|
||||||
|
:
|
||||||
|
Wall(dict)
|
||||||
|
{
|
||||||
|
if(!readcuboidWall(dict))
|
||||||
|
{
|
||||||
|
fatalExit;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
@ -0,0 +1,60 @@
|
||||||
|
/*------------------------------- 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 __cuboidWall_H__
|
||||||
|
#define __cuboidWall_H__
|
||||||
|
|
||||||
|
|
||||||
|
#include "Wall.H"
|
||||||
|
#include "types.H"
|
||||||
|
|
||||||
|
namespace pFlow
|
||||||
|
{
|
||||||
|
|
||||||
|
class cuboidWall
|
||||||
|
:
|
||||||
|
public Wall
|
||||||
|
{
|
||||||
|
protected:
|
||||||
|
|
||||||
|
bool readcuboidWall(const dictionary& dict);
|
||||||
|
|
||||||
|
public:
|
||||||
|
|
||||||
|
TypeName("cuboidWall");
|
||||||
|
|
||||||
|
cuboidWall();
|
||||||
|
|
||||||
|
cuboidWall(const dictionary& dict);
|
||||||
|
|
||||||
|
add_vCtor
|
||||||
|
(
|
||||||
|
Wall,
|
||||||
|
cuboidWall,
|
||||||
|
dictionary
|
||||||
|
);
|
||||||
|
|
||||||
|
};
|
||||||
|
|
||||||
|
} // pFlow
|
||||||
|
|
||||||
|
|
||||||
|
#endif
|
||||||
|
|
@ -58,6 +58,7 @@ bool pFlow::planeWall::readPlaneWall
|
||||||
|
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
||||||
pFlow::planeWall::planeWall()
|
pFlow::planeWall::planeWall()
|
||||||
{}
|
{}
|
||||||
|
|
||||||
|
|
@ -73,3 +74,31 @@ pFlow::planeWall::planeWall
|
||||||
fatalExit;
|
fatalExit;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
|
pFlow::planeWall::planeWall(
|
||||||
|
const realx3& p1,
|
||||||
|
const realx3& p2,
|
||||||
|
const realx3& p3,
|
||||||
|
const realx3& p4)
|
||||||
|
{
|
||||||
|
|
||||||
|
if( Wall::checkTrianlge(p1,p2,p3) )
|
||||||
|
{
|
||||||
|
triangles_.push_back(realx3x3(p1,p2,p3));
|
||||||
|
}else
|
||||||
|
{
|
||||||
|
fatalErrorInFunction <<
|
||||||
|
"points p1, p2 and p3 do not form a plane wall "<<endl;
|
||||||
|
fatalExit;
|
||||||
|
}
|
||||||
|
|
||||||
|
if( Wall::checkTrianlge(p3,p4,p1) )
|
||||||
|
{
|
||||||
|
triangles_.push_back(realx3x3(p3,p4,p1));
|
||||||
|
}else
|
||||||
|
{
|
||||||
|
fatalErrorInFunction <<
|
||||||
|
"points p3, p4 and p1 do not form a plane wallendl";
|
||||||
|
fatalExit;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
@ -45,6 +45,12 @@ public:
|
||||||
|
|
||||||
planeWall(const dictionary& dict);
|
planeWall(const dictionary& dict);
|
||||||
|
|
||||||
|
planeWall(
|
||||||
|
const realx3& p1,
|
||||||
|
const realx3& p2,
|
||||||
|
const realx3& p3,
|
||||||
|
const realx3& p4);
|
||||||
|
|
||||||
add_vCtor
|
add_vCtor
|
||||||
(
|
(
|
||||||
Wall,
|
Wall,
|
||||||
|
|
|
||||||
Loading…
Reference in New Issue