xfygogo/phasicFlow
1
0
Fork 0
mirror of https://github.com/PhasicFlow/phasicFlow.git synced 2025-06-12 16:26:23 +00:00
Code Issues Packages Projects Releases Wiki Activity

correctoin for particle insertion:cylinder region

Browse Source
This commit is contained in:
hamidrezanorouzi
2022-09-28 11:31:16 +03:30
parent f36a4d7782
commit 7485ed16a8
7 changed files with 23 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
utilities/geometryPhasicFlow/CMakeLists.txt
View File

@ -4,7 +4,6 @@ geometryPhasicFlow.C
Wall/Wall.C
planeWall/planeWall.C
stlWall/stlWall.C
cylinderWall/zAxis.C
cylinderWall/cylinderWall.C
)
set(link_lib phasicFlow Geometry Kokkos::kokkos)

150
utilities/geometryPhasicFlow/cylinderWall/zAxis.C
View File

@ -1,150 +0,0 @@
/*------------------------------- 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 "zAxis.H"
pFlow::zAxis::zAxis(const realx3 &p1, const realx3 &p2) :
p1_(p1),
p2_(p2)
{
n_ = p2-p1;
auto len = pFlow::length(n_);
if(len < smallValue )
{
fatalExit;
}
n_ /= len;
makeTransMatrix();
}
pFlow::realx3 pFlow::zAxis::transferToZ(const realx3 & p)
{
real pp[4][1] = {p.x(), p.y(), p.z(), 1.0};
real pn[4][1];
MatMul(Trans_z_xz_P1_, pp, pn);
return realx3(
pn[0][0],
pn[1][0],
pn[2][0]
);
}
pFlow::realx3 pFlow::zAxis::transferBackZ(const realx3 & p)
{
real pp[4][1] = { p.x(), p.y(), p.z(), 1.0 };
real pn[4][1];
MatMul(ITrans_P1_xz_z_, pp, pn);
return realx3(
pn[0][0],
pn[1][0],
pn[2][0]
);
}
void pFlow::zAxis::makeTransMatrix()
{
// transfering point p1 to the origin
real TransP1[4][4] =
{
1.0, 0.0, 0.0, -p1_.x(),
0.0, 1.0, 0.0, -p1_.y(),
0.0, 0.0, 1.0, -p1_.z(),
0.0, 0.0, 0.0, 1.0
};
// for transformation back
real ITransP1[4][4] =
{
1.0, 0.0, 0.0, p1_.x(),
0.0, 1.0, 0.0, p1_.y(),
0.0, 0.0, 1.0, p1_.z(),
0.0, 0.0, 0.0, 1.0
};
real u = n_.x();
real v = n_.y();
real w = n_.z();
real u2v2 = sqrt(u*u + v*v);
//correcting the transformation matrix in the case of coincidence with z - axis
if ( equal(w,1.0) )
{
assignMat(TransP1 , Trans_z_xz_P1_);
assignMat(ITransP1, ITrans_P1_xz_z_);
return;
}
u2v2 = max(smallValue, u2v2);
real TransXZ[4][4] =
{
u / u2v2, v / u2v2, 0.0, 0.0,
-v / u2v2, u / u2v2, 0.0, 0.0,
0.0, 0.0, 1.0, 0.0,
0.0, 0.0, 0.0, 1.0
};
real TransZ[4][4] =
{
w, 0.0, -u2v2, 0.0,
0.0, 1.0, 0.0, 0.0,
u2v2, 0.0, w, 0.0,
0.0, 0.0, 0.0, 1.0
};
real temp[4][4];
// creat transformation matrix to transfer point from line axis to z-axis
MatMul(TransXZ, TransP1, temp);
MatMul(TransZ, temp, Trans_z_xz_P1_);
real ITransXZ[4][4] =
{
u / u2v2, -v / u2v2, 0.0, 0.0,
+v / u2v2, u / u2v2, 0.0, 0.0,
0.0, 0.0, 1.0, 0.0,
0.0, 0.0, 0.0, 1.0
};
real ITransZ[4][4] =
{
w, 0.0, +u2v2, 0.0,
0.0, 1.0, 0.0, 0.0,
-u2v2, 0.0, w, 0.0,
0.0, 0.0, 0.0, 1.0
};
// creat transformation matrix to transfer point to from z-axis to line axis
MatMul(ITransXZ, ITransZ, temp);
MatMul(ITransP1, temp, ITrans_P1_xz_z_);
}

107
utilities/geometryPhasicFlow/cylinderWall/zAxis.H
View File

@ -1,107 +0,0 @@
/*------------------------------- 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.
-----------------------------------------------------------------------------*/
/*
class name: zAxis
a class for transfering points to a cylinderical coordinates defined by points
p1 and p2
*/
#ifndef __zAxis_H__
#define __zAxis_H__
#include "types.H"
namespace pFlow
{
template <typename T, int32 nRow, int32 nInner, int32 nCol >
void MatMul(T(&A)[nRow][nInner], T(&B)[nInner][nCol], T(&C)[nRow][nCol]);
template <typename T, int32 nRow, int32 nCol >
void assignMat(T(&A)[nRow][nCol], T(&B)[nRow][nCol]);
class zAxis
{
public:
// constructors
zAxis(const realx3 &lp1, const realx3 &lp2);
real length()const
{
return pFlow::length(p2_-p1_);
}
realx3 transferToZ(const realx3 & p);
realx3 transferBackZ(const realx3 & p);
private:
void makeTransMatrix();
protected:
realx3 p1_;
realx3 p2_;
realx3 n_;
real Trans_z_xz_P1_[4][4];
real ITrans_P1_xz_z_[4][4];
};
template <typename T, int32 nRow, int32 nInner, int32 nCol >
void MatMul(T(&A)[nRow][nInner], T(&B)[nInner][nCol], T(&C)[nRow][nCol])
{
for (int32 row = 0; row < nRow; row++)
{
for (int32 col = 0; col < nCol; col++)
{
T sum = 0;
for (int inner = 0; inner < nInner; inner++)
{
sum += A[row][inner] * B[inner][col];
}
C[row][col] = sum;
}
}
}
template <typename T, int32 nRow, int32 nCol >
void assignMat(T(&A)[nRow][nCol], T(&B)[nRow][nCol])
{
for (int32 row = 0; row < nRow; row++)
{
for (int32 col = 0; col < nCol; col++)
{
B[row][col] = A[row][col];
}
}
}
}
#endif