#include <frameChangeF3D.h>

Inheritance diagram for olb::RectangleTrigonometricPoiseuille3D< T >:

Collaboration diagram for olb::RectangleTrigonometricPoiseuille3D< T >:

Public Member Functions
	RectangleTrigonometricPoiseuille3D (std::vector< T > &x0_, std::vector< T > &x1_, std::vector< T > &x2_, std::vector< T > &maxVelocity, int numOfPolynoms_=1)

	RectangleTrigonometricPoiseuille3D (SuperGeometry< T, 3 > &superGeometry_, int material_, std::vector< T > &maxVelocity_, T offsetX, T offsetY, T offsetZ, int numOfPolynoms_=1)
	constructor from material numbers

T	getProfilePeak ()

T	getRatioAvgToPeak ()

bool	operator() (T output[], const T x[]) override

Public Member Functions inherited from olb::AnalyticalF< D, T, S >
AnalyticalF< D, T, S > &	operator- (AnalyticalF< D, T, S > &rhs)

AnalyticalF< D, T, S > &	operator+ (AnalyticalF< D, T, S > &rhs)

AnalyticalF< D, T, S > &	operator* (AnalyticalF< D, T, S > &rhs)

AnalyticalF< D, T, S > &	operator/ (AnalyticalF< D, T, S > &rhs)

Public Member Functions inherited from olb::GenericF< T, S >
virtual	~GenericF ()=default

int	getSourceDim () const
	read only access to member variable _m

int	getTargetDim () const
	read only access to member variable _n

std::string &	getName ()
	read and write access to name

std::string const &	getName () const
	read only access to name

virtual bool	operator() (T output[], const S input[])=0
	has to be implemented for 'every' derived class

bool	operator() (T output[])
	wrapper that call the pure virtual operator() (T output[], const S input[]) from above

bool	operator() (T output[], S input0)

bool	operator() (T output[], S input0, S input1)

bool	operator() (T output[], S input0, S input1, S input2)

bool	operator() (T output[], S input0, S input1, S input2, S input3)

Protected Member Functions
void	calcPeakAndAvg ()

Protected Member Functions inherited from olb::AnalyticalF< D, T, S >
	AnalyticalF (int n)

Protected Member Functions inherited from olb::GenericF< T, S >
	GenericF (int targetDim, int sourceDim)

Protected Attributes
OstreamManager	clout

std::vector< T >	x0

std::vector< T >	x1

std::vector< T >	x2

std::vector< T >	maxVelocity

int	numOfPolynoms

T	profilePeak

T	profileRatioAvgToPeak

Additional Inherited Members
Public Types inherited from olb::AnalyticalF< D, T, S >
using	identity_functor_type = AnalyticalIdentity<D,T,S>

Public Types inherited from olb::GenericF< T, S >
using	targetType = T

using	sourceType = S

Public Attributes inherited from olb::GenericF< T, S >
std::shared_ptr< GenericF< T, S > >	_ptrCalcC
	memory management, frees resouces (calcClass)

Static Public Attributes inherited from olb::AnalyticalF< D, T, S >
static constexpr unsigned	dim = D

Detailed Description

template<typename T>
class olb::RectangleTrigonometricPoiseuille3D< T >

Definition at line 264 of file frameChangeF3D.h.

Constructor & Destructor Documentation

◆ RectangleTrigonometricPoiseuille3D() [1/2]

template<typename T >

olb::RectangleTrigonometricPoiseuille3D< T >::RectangleTrigonometricPoiseuille3D	(	std::vector< T > &	x0_,
		std::vector< T > &	x1_,
		std::vector< T > &	x2_,
		std::vector< T > &	maxVelocity,
		int	numOfPolynoms_ = 1 )

Definition at line 518 of file frameChangeF3D.hh.

  : AnalyticalF3D<T,T>(3), clout(std::cout,"RectangleTrigonometricPoiseille3D"), x0(x0_), x1(x1_),
    x2(x2_), maxVelocity(maxVelocity_), numOfPolynoms(numOfPolynoms_)
{
  calcPeakAndAvg();
}

References olb::RectangleTrigonometricPoiseuille3D< T >::calcPeakAndAvg().

Here is the call graph for this function:

◆ RectangleTrigonometricPoiseuille3D() [2/2]

template<typename T >

olb::RectangleTrigonometricPoiseuille3D< T >::RectangleTrigonometricPoiseuille3D	(	SuperGeometry< T, 3 > &	superGeometry_,
		int	material_,
		std::vector< T > &	maxVelocity_,
		T	offsetX,
		T	offsetY,
		T	offsetZ,
		int	numOfPolynoms_ = 1 )

constructor from material numbers

offsetX,Y,Z is a positive number denoting the distance from border voxels of material_ to the zerovelocity boundary

Definition at line 528 of file frameChangeF3D.hh.

  : AnalyticalF3D<T,T>(3), clout(std::cout, "RectangleTrigonometricPoiseuille3D"), maxVelocity(maxVelocity_),
    numOfPolynoms(numOfPolynoms_)
{
  std::vector<T> min = superGeometry_.getStatistics().getMinPhysR(material_);
  std::vector<T> max = superGeometry_.getStatistics().getMaxPhysR(material_);
  // set three corners of the plane, move by offset to land on the v=0
  // boundary and adapt certain values depending on the orthogonal axis to get
  // different points
  x0 = min;
  x1 = min;
  x2 = min;
  if ( util::nearZero(max[0]-min[0]) ) { // orthogonal to x-axis
    x0[1] -= offsetY;
    x0[2] -= offsetZ;
    x1[1] -= offsetY;
    x1[2] -= offsetZ;
    x2[1] -= offsetY;
    x2[2] -= offsetZ;
 
    x1[1] = max[1] + offsetY;
    x2[2] = max[2] + offsetZ;
  }
  else if ( util::nearZero(max[1]-min[1]) ) {   // orthogonal to y-axis
    x0[0] -= offsetX;
    x0[2] -= offsetZ;
    x1[0] -= offsetX;
    x1[2] -= offsetZ;
    x2[0] -= offsetX;
    x2[2] -= offsetZ;
 
    x1[0] = max[0] + offsetX;
    x2[2] = max[2] + offsetZ;
  }
  else if ( util::nearZero(max[2]-min[2]) ) {   // orthogonal to z-axis
    x0[0] -= offsetX;
    x0[1] -= offsetY;
    x1[0] -= offsetX;
    x1[1] -= offsetY;
    x2[0] -= offsetX;
    x2[1] -= offsetY;
 
    x1[0] = max[0] + offsetX;
    x2[1] = max[1] + offsetY;
  }
  else {
    clout << "Error: constructor from material number works only for axis-orthogonal planes" << std::endl;
  }
  calcPeakAndAvg();
}

References olb::RectangleTrigonometricPoiseuille3D< T >::calcPeakAndAvg(), olb::RectangleTrigonometricPoiseuille3D< T >::clout, olb::SuperGeometry< T, D >::getStatistics(), olb::util::nearZero(), olb::RectangleTrigonometricPoiseuille3D< T >::x0, olb::RectangleTrigonometricPoiseuille3D< T >::x1, and olb::RectangleTrigonometricPoiseuille3D< T >::x2.

Here is the call graph for this function:

Member Function Documentation

◆ calcPeakAndAvg()

template<typename T >

void olb::RectangleTrigonometricPoiseuille3D< T >::calcPeakAndAvg ( )

protected

Definition at line 582 of file frameChangeF3D.hh.

{
  //2006_Bruus_Lecture Theoretical microfluidics
  T sumMax = 0;
  T sumAvg = 0;
  for (int i=1; i<=numOfPolynoms; ++i) { //TODO: higher polynoms still to be tested
    T n = i*2-1;   //Creating only odd indices
    T A = util::sin(n*M_PI*0.5);
    T max = (1./util::pow(n,3)) * (1.-(1 / util::cosh(n*M_PI / 2 ))) * A;
    //TODO: only for y,z=1 //TODO: Maximum calculation needs to be adapted for different edge ratios
    T avg = 2. * A * A * (n*M_PI-2*util::tanh(n*M_PI*0.5)) / (util::pow(n,5)*M_PI*M_PI);
    sumMax += max;
    sumAvg += avg;
  }
  profilePeak = sumMax;
  profileRatioAvgToPeak = sumAvg / sumMax;
}

References olb::util::cosh(), M_PI, olb::util::pow(), olb::util::sin(), and olb::util::tanh().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getProfilePeak()

template<typename T >

T olb::RectangleTrigonometricPoiseuille3D< T >::getProfilePeak ( )

Definition at line 602 of file frameChangeF3D.hh.

{
  return profilePeak;
}

◆ getRatioAvgToPeak()

template<typename T >

T olb::RectangleTrigonometricPoiseuille3D< T >::getRatioAvgToPeak ( )

Definition at line 608 of file frameChangeF3D.hh.

{
  return profileRatioAvgToPeak;
}

◆ operator()()

template<typename T >

bool olb::RectangleTrigonometricPoiseuille3D< T >::operator()	(	T	output[],
		const T	x[] )

override

Definition at line 614 of file frameChangeF3D.hh.

{
  OstreamManager clout( std::cout,"RectPois" ); //TODO: to be removed
 
  // create plane normals and supports, (E1,E2) and (E3,E4) are parallel
  std::vector<std::vector<T> > n(4,std::vector<T>(3,0)); // normal vectors
  std::vector<std::vector<T> > s(4,std::vector<T>(3,0)); // support vectors
  for (int iD = 0; iD < 3; iD++) {
    n[0][iD] =  x1[iD] - x0[iD];
    n[1][iD] = -x1[iD] + x0[iD];
    n[2][iD] =  x2[iD] - x0[iD];
    n[3][iD] = -x2[iD] + x0[iD];
    s[0][iD] = x0[iD];
    s[1][iD] = x1[iD];
    s[2][iD] = x0[iD];
    s[3][iD] = x2[iD];
  }
  for (int iP = 0; iP < 4; iP++) {
    n[iP] = util::normalize(n[iP]);
  }
 
  // determine plane coefficients in the coordinate equation E_i: Ax+By+Cz+D=0
  // given form: (x-s)*n=0 => A=n1, B=n2, C=n3, D=-(s1n1+s2n2+s3n3)
  std::vector<T> A(4,0);
  std::vector<T> B(4,0);
  std::vector<T> C(4,0);
  std::vector<T> D(4,0);
 
  for (int iP = 0; iP < 4; iP++) {
    A[iP] = n[iP][0];
    B[iP] = n[iP][1];
    C[iP] = n[iP][2];
    D[iP] = -(s[iP][0]*n[iP][0] + s[iP][1]*n[iP][1] + s[iP][2]*n[iP][2]);
  }
 
  // determine distance to the walls by formula
  std::vector<T> d(4,0);
  T aabbcc(0);
  for (int iP = 0; iP < 4; iP++) {
    aabbcc = A[iP]*A[iP] + B[iP]*B[iP] + C[iP]*C[iP];
    d[iP] = util::fabs(A[iP]*x[0]+B[iP]*x[1]+C[iP]*x[2]+D[iP])*util::pow(aabbcc,-0.5);
  }
 
  // length and width of the rectangle
  std::vector<T> l(2,0);
  l[0] = d[0] + d[1] ;
  l[1] = d[2] + d[3] ;
 
  std::vector<T> dl(2,0);
  dl[0] = d[0] / l[0];
  dl[1] = d[2] / l[1];
 
  //Polynomial creation for odd indices (2i-i) from:
  //2006_Bruus_Lecture Theoretical microfluidics
  T sumPoly = 0;
  for (int i=1; i<=numOfPolynoms; ++i) { //TODO: higher polynoms still to be tested
    T n = i*2-1;
    T poly = (1./util::pow(n,3)) * (1.-(util::cosh(n*M_PI*(dl[0]-0.5)) / util::cosh(n*M_PI/ (2) ))) * util::sin(n*M_PI*dl[1]);
    sumPoly += poly;
  }
 
  output[0] = maxVelocity[0] * ( sumPoly / profilePeak);
  output[1] = maxVelocity[1] * ( sumPoly / profilePeak);
  output[2] = maxVelocity[2] * ( sumPoly / profilePeak);
  return true;
}