#include <navierStokesAdvectionDiffusionCouplingPostProcessor3D.h>

Inheritance diagram for olb::VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D< T, DESCRIPTOR, POROSITY, ADLattice, FIELD_A, FIELD_B >:

Collaboration diagram for olb::VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D< T, DESCRIPTOR, POROSITY, ADLattice, FIELD_A, FIELD_B >:

Public Member Functions
	VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D (int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, int iC_, std::vector< BlockStructureD< 3 > * > partners_, std::vector< std::reference_wrapper< AdvectionDiffusionForce3D< T, DESCRIPTOR, ADLattice > > > forces_)

int	extent () const override
	Extent of application area (0 for purely local operations)

int	extent (int whichDirection) const override
	Extent of application area along a direction (0 or 1)

void	process (BlockLattice< T, DESCRIPTOR > &blockLattice) override
	Execute post-processing step.

void	processSubDomain (BlockLattice< T, DESCRIPTOR > &blockLattice, int x0_, int x1_, int y0_, int y1_, int z0_, int z1_) override
	Execute post-processing step on a sublattice.

Public Member Functions inherited from olb::PostProcessor3D< T, DESCRIPTOR >
	PostProcessor3D ()

virtual	~PostProcessor3D ()

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

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

int	getPriority () const
	read only access to priority

Protected Attributes
std::vector< std::reference_wrapper< AdvectionDiffusionForce3D< T, DESCRIPTOR, ADLattice > > >	_forces

Protected Attributes inherited from olb::PostProcessor3D< T, DESCRIPTOR >
int	_priority

Detailed Description

template<typename T, typename DESCRIPTOR, typename POROSITY, typename ADLattice, typename FIELD_A, typename FIELD_B>
class olb::VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D< T, DESCRIPTOR, POROSITY, ADLattice, FIELD_A, FIELD_B >

Definition at line 359 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.h.

Constructor & Destructor Documentation

◆ VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D()

template<typename T , typename DESCRIPTOR , typename POROSITY , typename ADLattice , typename FIELD_A , typename FIELD_B >

olb::VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D< T, DESCRIPTOR, POROSITY, ADLattice, FIELD_A, FIELD_B >::VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D	(	int	x0_,
		int	x1_,
		int	y0_,
		int	y1_,
		int	z0_,
		int	z1_,
		int	iC_,
		std::vector< BlockStructureD< 3 > * >	partners_,
		std::vector< std::reference_wrapper< AdvectionDiffusionForce3D< T, DESCRIPTOR, ADLattice > > >	forces_ )

Definition at line 998 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.hh.

  :  _forces(forces_),
     x0(x0_), x1(x1_), y0(y0_), y1(y1_), z0(z0_), z1(z1_), iC(iC_),
     _partnerLattice(static_cast<BlockLattice<T,ADLattice> *>(partners_[0])),
     _cell(_partnerLattice->get(x0,y0,z0)),
     _cellXp(_partnerLattice->get(x0+1,y0,z0)),
     _cellXn(_partnerLattice->get(x0-1,y0,z0)),
     _cellYp(_partnerLattice->get(x0,y0+1,z0)),
     _cellYn(_partnerLattice->get(x0,y0-1,z0)),
     _cellZp(_partnerLattice->get(x0,y0,z0+1)),
     _cellZn(_partnerLattice->get(x0,y0,z0-1))
{
  this->getName() = "VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D";
}

References olb::PostProcessor3D< T, DESCRIPTOR >::getName().

Here is the call graph for this function:

Member Function Documentation

◆ extent() [1/2]

template<typename T , typename DESCRIPTOR , typename POROSITY , typename ADLattice , typename FIELD_A , typename FIELD_B >

int olb::VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D< T, DESCRIPTOR, POROSITY, ADLattice, FIELD_A, FIELD_B >::extent ( ) const

inlineoverridevirtual

Extent of application area (0 for purely local operations)

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 366 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.h.

  {
    return 1;
  }

◆ extent() [2/2]

template<typename T , typename DESCRIPTOR , typename POROSITY , typename ADLattice , typename FIELD_A , typename FIELD_B >

int olb::VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D< T, DESCRIPTOR, POROSITY, ADLattice, FIELD_A, FIELD_B >::extent ( int direction ) const

inlineoverridevirtual

Extent of application area along a direction (0 or 1)

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 370 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.h.

  {
    return 1;
  }

◆ process()

template<typename T , typename DESCRIPTOR , typename POROSITY , typename ADLattice , typename FIELD_A , typename FIELD_B >

void olb::VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D< T, DESCRIPTOR, POROSITY, ADLattice, FIELD_A, FIELD_B >::process ( BlockLattice< T, DESCRIPTOR > & blockLattice )

overridevirtual

Execute post-processing step.

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 1153 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.hh.

{
  processSubDomain(blockLattice, x0, x1, y0, y1, z0, z1);
}

◆ processSubDomain()

template<typename T , typename DESCRIPTOR , typename POROSITY , typename ADLattice , typename FIELD_A , typename FIELD_B >

void olb::VolumeAveragedNavierStokesAdvectionDiffusionParticleCouplingPostProcessor3D< T, DESCRIPTOR, POROSITY, ADLattice, FIELD_A, FIELD_B >::processSubDomain	(	BlockLattice< T, DESCRIPTOR > &	blockLattice,
		int	x0_,
		int	x1_,
		int	y0_,
		int	y1_,
		int	z0_,
		int	z1_ )

overridevirtual

Execute post-processing step on a sublattice.

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 1017 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.hh.

{
  auto vel     = par ? _cell.template getField<FIELD_B>() : _cell.template getField<FIELD_A>();
  auto vel_new = par ? _cell.template getFieldPointer<FIELD_A>() : _cell.template getFieldPointer<FIELD_B>();
 
  auto velXp = par ? _cellXp.template getFieldPointer<FIELD_B>() : _cellXp.template getFieldPointer<FIELD_A>();
  auto velXn = par ? _cellXn.template getFieldPointer<FIELD_B>() : _cellXn.template getFieldPointer<FIELD_A>();
  auto velYp = par ? _cellYp.template getFieldPointer<FIELD_B>() : _cellYp.template getFieldPointer<FIELD_A>();
  auto velYn = par ? _cellYn.template getFieldPointer<FIELD_B>() : _cellYn.template getFieldPointer<FIELD_A>();
  auto velZp = par ? _cellZp.template getFieldPointer<FIELD_B>() : _cellZp.template getFieldPointer<FIELD_A>();
  auto velZn = par ? _cellZn.template getFieldPointer<FIELD_B>() : _cellZn.template getFieldPointer<FIELD_A>();
 
  auto forceNS = blockLattice.get(x0,y0,z0).template getFieldPointer<descriptors::FORCE>();
 
  auto _cellNSXp = blockLattice.get(x0+1,y0,z0);
  auto _cellNSXn = blockLattice.get(x0-1,y0,z0);
  auto _cellNSYp = blockLattice.get(x0,y0+1,z0);
  auto _cellNSYn = blockLattice.get(x0,y0-1,z0);
  auto _cellNSZp = blockLattice.get(x0,y0,z0+1);
  auto _cellNSZn = blockLattice.get(x0,y0,z0-1);
 
  int newX0, newX1, newY0, newY1, newZ0, newZ1;
  T porosityForce[3];
  T pressure;
 
  if ( util::intersect (
         x0, x1, y0, y1, z0, z1,
         x0_, x1_, y0_, y1_, z0_, z1_,
         newX0, newX1, newY0, newY1, newZ0, newZ1 ) ) {
 
    for (int iX=newX0; iX<=newX1; ++iX) {
      for (int iY=newY0; iY<=newY1; ++iY) {
        for (int iZ=newZ0; iZ<=newZ1; ++iZ) {
      auto nsCell = blockLattice.get(iX,iY,iZ);
      T porosityNS = 1. - _cell.computeRho();
      nsCell.template setField<POROSITY>(porosityNS);
    }
      }
    }
 
    for (int iX=newX0; iX<=newX1; ++iX) {
      for (int iY=newY0; iY<=newY1; ++iY) {
        for (int iZ=newZ0; iZ<=newZ1; ++iZ) {
          int latticeR[4] = {iC, iX, iY, iZ};
          T velGrad[3] = {0.,0.,0.};
          T forceValue[3] = {0.,0.,0.};
          T velF[3] = {0.,0.,0.};
 
          auto nsCell = blockLattice.get(iX,iY,iZ);
          pressure = nsCell.computeRho() / descriptors::invCs2<T,DESCRIPTOR>();
 
          if (_forces.begin() != _forces.end()) {
            // calculating upwind Gradient
            // vel contains velocity information on ADlattice
            // velGrad contains upwind vel on ADlattice
            if (vel[0]<0.) {
              velGrad[0] = vel[0]*(velXp[0]-vel[0]);
              velGrad[1] = vel[0]*(velXp[1]-vel[1]);
              velGrad[2] = vel[0]*(velXp[2]-vel[2]);
            }
            else {
              velGrad[0] = vel[0]*(vel[0]-velXn[0]);
              velGrad[1] = vel[0]*(vel[1]-velXn[1]);
              velGrad[2] = vel[0]*(vel[2]-velXn[2]);
            }
            if (vel[1]<0.) {
              velGrad[0] += vel[1]*(velYp[0]-vel[0]);
              velGrad[1] += vel[1]*(velYp[1]-vel[1]);
              velGrad[2] += vel[1]*(velYp[2]-vel[2]);
            }
            else {
              velGrad[0] += vel[1]*(vel[0]-velYn[0]);
              velGrad[1] += vel[1]*(vel[1]-velYn[1]);
              velGrad[2] += vel[1]*(vel[2]-velYn[2]);
            }
            if (vel[2]<0.) {
              velGrad[0] += vel[2]*(velZp[0]-vel[0]);
              velGrad[1] += vel[2]*(velZp[1]-vel[1]);
              velGrad[2] += vel[2]*(velZp[2]-vel[2]);
            }
            else {
              velGrad[0] += vel[2]*(vel[0]-velZn[0]);
              velGrad[1] += vel[2]*(vel[1]-velZn[1]);
              velGrad[2] += vel[2]*(vel[2]-velZn[2]);
            }
 
            for (AdvectionDiffusionForce3D<T, DESCRIPTOR, ADLattice>& f : _forces) {
              // writes force in forceValues, vel refers to ADlattice
              auto adCell = _partnerLattice->get(x0,y0,z0);
              f.applyForce(forceValue, &nsCell, &adCell, vel.data(), latticeR);
              if (par) {
                _cell.template setField<FIELD_B>(vel);
              } else {
                _cell.template setField<FIELD_A>(vel);
              }
            }
 
            // compute new particle velocity and opposite fluid force
            for (int i=0; i < DESCRIPTOR::d; i++) {
              vel_new[i] = vel[i] + forceValue[i] - velGrad[i];
              forceNS[i] = -forceValue[i];
            }
          }
          else {
 
        T porXp = _cellNSXp.template getField<POROSITY>();
            T porXn = _cellNSXn.template getField<POROSITY>();
            T porYp = _cellNSYp.template getField<POROSITY>();
            T porYn = _cellNSYn.template getField<POROSITY>();
            T porZp = _cellNSZp.template getField<POROSITY>();
            T porZn = _cellNSZn.template getField<POROSITY>();
 
            porosityForce[0] = 0.5 * pressure * (porXp - porXn);
        porosityForce[1] = 0.5 * pressure * (porYp - porYn);
            porosityForce[2] = 0.5 * pressure * (porZp - porZn);
 
        for (int i = 0; i < DESCRIPTOR::d; i++) {
          forceNS[i] += porosityForce[i];
        }
 
            nsCell.computeU(velF);
 
            for (int i = 0; i < DESCRIPTOR::d; i++) {   // set particle velocity to fluid velocity
              vel_new[i] = velF[i];
            }
          }
        }
      }
    }
  }
  par = !par;
}