olb::TotalEnthalpyPhaseChangeCouplingPostProcessor3D< T, DESCRIPTOR, DYNAMICS > Class Template Reference

#include <navierStokesAdvectionDiffusionCouplingPostProcessor3D.h>

Inheritance diagram for olb::TotalEnthalpyPhaseChangeCouplingPostProcessor3D< T, DESCRIPTOR, DYNAMICS >:

Collaboration diagram for olb::TotalEnthalpyPhaseChangeCouplingPostProcessor3D< T, DESCRIPTOR, DYNAMICS >:

Public Member Functions
	TotalEnthalpyPhaseChangeCouplingPostProcessor3D (int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, T gravity_, T T0_, T deltaTemp_, std::vector< T > dir_, std::vector< BlockStructureD< 3 > * > partners_)
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

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

Detailed Description

template<typename T, typename DESCRIPTOR, typename DYNAMICS>
class olb::TotalEnthalpyPhaseChangeCouplingPostProcessor3D< T, DESCRIPTOR, DYNAMICS >

Definition at line 41 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.h.

Constructor & Destructor Documentation

TotalEnthalpyPhaseChangeCouplingPostProcessor3D()

template<typename T , typename DESCRIPTOR , typename DYNAMICS >

olb::TotalEnthalpyPhaseChangeCouplingPostProcessor3D< T, DESCRIPTOR, DYNAMICS >::TotalEnthalpyPhaseChangeCouplingPostProcessor3D	(	int	x0_,
		int	x1_,
		int	y0_,
		int	y1_,
		int	z0_,
		int	z1_,
		T	gravity_,
		T	T0_,
		T	deltaTemp_,
		std::vector< T >	dir_,
		std::vector< BlockStructureD< 3 > * >	partners_ )

Definition at line 43 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.hh.

  :  x0(x0_), x1(x1_), y0(y0_), y1(y1_), z0(z0_), z1(z1_),
     gravity(gravity_), T0(T0_), deltaTemp(deltaTemp_),
     dir(dir_), partners(partners_)
{
  this->getName() = "TotalEnthalpyPhaseChangeCouplingPostProcessor3D";
  // we normalize the direction of force vector
  T normDir = T();
  for (unsigned iD = 0; iD < dir.size(); ++iD) {
    normDir += dir[iD]*dir[iD];
  }
  normDir = util::sqrt(normDir);
  for (unsigned iD = 0; iD < dir.size(); ++iD) {
    dir[iD] /= normDir;
  }
 
  for (unsigned iD = 0; iD < dir.size(); ++iD) {
    forcePrefactor[iD] = gravity * dir[iD];
  }
 
  tPartner = static_cast<BlockLattice<T,descriptors::D3Q7<descriptors::VELOCITY,descriptors::TEMPERATURE>> *>(partners[0]);
}

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

Here is the call graph for this function:

Member Function Documentation

extent() [1/2]

template<typename T , typename DESCRIPTOR , typename DYNAMICS >

int olb::TotalEnthalpyPhaseChangeCouplingPostProcessor3D< T, DESCRIPTOR, DYNAMICS >::extent ( ) const

inlineoverridevirtual

Extent of application area (0 for purely local operations)

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 46 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.h.

  {
    return 1;
  }

extent() [2/2]

template<typename T , typename DESCRIPTOR , typename DYNAMICS >

int olb::TotalEnthalpyPhaseChangeCouplingPostProcessor3D< T, DESCRIPTOR, DYNAMICS >::extent ( int direction ) const

inlineoverridevirtual

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

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 50 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.h.

  {
    return 1;
  }

process()

template<typename T , typename DESCRIPTOR , typename DYNAMICS >

void olb::TotalEnthalpyPhaseChangeCouplingPostProcessor3D< T, DESCRIPTOR, DYNAMICS >::process ( BlockLattice< T, DESCRIPTOR > & blockLattice )

overridevirtual

Execute post-processing step.

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 114 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.hh.

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

processSubDomain()

template<typename T , typename DESCRIPTOR , typename DYNAMICS >

void olb::TotalEnthalpyPhaseChangeCouplingPostProcessor3D< T, DESCRIPTOR, DYNAMICS >::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 70 of file navierStokesAdvectionDiffusionCouplingPostProcessor3D.hh.

{
 
  int newX0, newX1, newY0, newY1, newZ0, newZ1;
  if ( util::intersect (
         x0, x1, y0, y1, z0, z1,
         x0_, x1_, y0_, y1_, z0_, z1_,
         newX0, newX1, newY0, newY1, newZ0, newZ1 ) ) {
    auto* dynamics = static_cast<DYNAMICS*>(tPartner->template getDynamics<DYNAMICS>());
    auto& parameters = static_cast<ParametersOfDynamicsD<DYNAMICS>&>(
      tPartner->template getData<OperatorParameters<DYNAMICS>>());
 
    for (int iX=newX0; iX<=newX1; ++iX) {
      for (int iY=newY0; iY<=newY1; ++iY) {
        for (int iZ=newZ0; iZ<=newZ1; ++iZ) {
          auto cell = blockLattice.get(iX,iY,iZ);
          auto partnerCell = tPartner->get(iX,iY,iZ);
 
          T enthalpy = partnerCell.computeRho();
 
          cell.template setField<descriptors::POROSITY>(
            dynamics->template computeLiquidFraction<T>(parameters, enthalpy));
          auto temperature = partnerCell.template getFieldPointer<descriptors::TEMPERATURE>();
          temperature[0] = dynamics->template computeTemperature<T>(parameters, enthalpy);
 
          // computation of the bousinessq force
          auto force = blockLattice.get(iX,iY,iZ).template getFieldPointer<descriptors::FORCE>();
          T temperatureDifference = temperature[0] - T0;
          for (unsigned iD = 0; iD < L::d; ++iD) {
            force[iD] = forcePrefactor[iD] * temperatureDifference;
          }
          // Velocity coupling
          FieldD<T,DESCRIPTOR,descriptors::VELOCITY> u;
          blockLattice.get(iX,iY,iZ).computeU(u.data());
          tPartner->get(iX,iY,iZ).template setField<descriptors::VELOCITY>(u);
        }
      }
    }
  }
}

References olb::Vector< T, D >::data(), olb::BlockLattice< T, DESCRIPTOR >::get(), and olb::util::intersect().

Here is the call graph for this function:

---

The documentation for this class was generated from the following files:

• src/dynamics/navierStokesAdvectionDiffusionCouplingPostProcessor3D.h
• src/dynamics/navierStokesAdvectionDiffusionCouplingPostProcessor3D.hh