#include <cage.h>
Public Member Functions | |
cage () | |
void | solve (int numberOfCage, int cellNumberPerPatch, int dimension, int nTruncate, real ReynoldsNumber, real *xParticleCenter, real *yParticleCenter, real *zParticleCenter, real ***xP, real ****p, int **pointNumberMeshP, int ****pCagePoint, real ***xU, real ****u, int **pointNumberMeshU, int ****uCagePoint, real ***xV, real ****v, int **pointNumberMeshV, int ****vCagePoint, real ***xW, real ****w, int **pointNumberMeshW, int ****wCagePoint, real ***xOmegaX, real ****omegaX, int **pointNumberMeshOmegaX, int ****omegaXCagePoint, real ***xOmegaY, real ****omegaY, int **pointNumberMeshOmegaY, int ****omegaYCagePoint, real ***xOmegaZ, real ****omegaZ, int **pointNumberMeshOmegaZ, int ****omegaZCagePoint, real **solvedABCoeffAllParticleNew, vector< real > &UFixedParticles, vector< real > &SFixedParticles, vector< real > &VFixedParticles, int &fixedParticles, int rotatingParticles, real **angVel, real **angVelAcc, real **torqueOld, real **torqueNew) |
Calculate coordinates of cage points and the velocity at cage points. | |
~cage () |
This class first calculate the coordinates of pressure, vorticity, and velocity at cage points. Then use their coordinates and values as the input to the class analyticalSol to calculate the velocity components at cage points. Last, use these valocity as the boundary conditions of the finite-difference scheme for the outer flows.
cage::cage | ( | ) |
constructor
cage::~cage | ( | ) |
destructor
void cage::solve | ( | int | numberOfCage, | |
int | cellNumberPerPatch, | |||
int | dimension, | |||
int | nTruncate, | |||
real | ReynoldsNumber, | |||
real * | xParticleCenter, | |||
real * | yParticleCenter, | |||
real * | zParticleCenter, | |||
real *** | xP, | |||
real **** | p, | |||
int ** | pointNumberMeshP, | |||
int **** | pCagePoint, | |||
real *** | xU, | |||
real **** | u, | |||
int ** | pointNumberMeshU, | |||
int **** | uCagePoint, | |||
real *** | xV, | |||
real **** | v, | |||
int ** | pointNumberMeshV, | |||
int **** | vCagePoint, | |||
real *** | xW, | |||
real **** | w, | |||
int ** | pointNumberMeshW, | |||
int **** | wCagePoint, | |||
real *** | xOmegaX, | |||
real **** | omegaX, | |||
int ** | pointNumberMeshOmegaX, | |||
int **** | omegaXCagePoint, | |||
real *** | xOmegaY, | |||
real **** | omegaY, | |||
int ** | pointNumberMeshOmegaY, | |||
int **** | omegaYCagePoint, | |||
real *** | xOmegaZ, | |||
real **** | omegaZ, | |||
int ** | pointNumberMeshOmegaZ, | |||
int **** | omegaZCagePoint, | |||
real ** | solvedABCoeffAllParticleNew, | |||
vector< real > & | UFixedParticles, | |||
vector< real > & | SFixedParticles, | |||
vector< real > & | VFixedParticles, | |||
int & | fixedParticles, | |||
int | rotatingParticles, | |||
real ** | angVel, | |||
real ** | angVelAcc, | |||
real ** | torqueOld, | |||
real ** | torqueNew | |||
) |
Calculate coordinates of cage points and the velocity at cage points.
Use the analytical solutions of the Stokes flows to calculate the the velocity at cage points, which acts as the inner boundary conditions for the finite-difference scheme.
numberOfCage | number of cages |
cellNumberPerPatch | cell numbers per patch surrounding the sphere |
dimension | dimension |
nTruncate | the order for truncated spectral analytical solution |
ReynoldsNumber | Reynolds Number |
xParticleCenter | x coordinate of particle center |
yParticleCenter | y coordinate of particle center |
zParticleCenter | z coordinate of particle center |
xP | coordinates of pressure |
p | pressure |
pointNumberMeshP | mesh points for pressure |
pCagePoint | cage points for pressure |
xU | coordinates of u |
u | u |
pointNumberMeshU | mesh points for u |
uCagePoint | cage points for u |
xV | coordinates of v |
v | v |
pointNumberMeshV | mesh points for v |
vCagePoint | cage points for v |
xW | coordinates of w |
w | w |
pointNumberMeshW | mesh points for w |
wCagePoint | cage points for w |
xOmegaX | coordinates of omega_x |
omegaX | omega_x |
pointNumberMeshOmegaX | mesh points for omega_x |
omegaXCagePoint | cage points for omega_x |
xOmegaY | coordinates of omega_y |
omegaY | omega_y |
pointNumberMeshOmegaY | mesh points for omega_y |
omegaYCagePoint | cage points for omega_y |
xOmegaZ | coordinates of omega_z |
omegaZ | omega_z |
pointNumberMeshOmegaZ | mesh points for omega_z |
omegaZCagePoint | cage points for omega_z |
solvedABCoeffAllParticleNew | the new values of A, B, and C |
UFixedParticles | store U for fixed partilce |
SFixedParticles | store S for fixed partilce |
VFixedParticles | store V for fixed partilce |
fixedParticles | if the partilce is fixed |
rotatingParticles | if the particle rotating |
angVel | angular velocity |
angVelAcc | angular velocity acceleration |
torqueOld | old torque |
torqueNew | new torque |