olb::ConcentrationAdvectionDiffusionCouplingPostProcessor3D< T, DESCRIPTOR > Class Template Reference

Coupling of ADlattice[0] with the other AD lattices (tpartners) More...

#include <advectionDiffusionReactionCouplingPostProcessor3D.h>

Inheritance diagram for olb::ConcentrationAdvectionDiffusionCouplingPostProcessor3D< T, DESCRIPTOR >:

Collaboration diagram for olb::ConcentrationAdvectionDiffusionCouplingPostProcessor3D< T, DESCRIPTOR >:

Public Member Functions
	ConcentrationAdvectionDiffusionCouplingPostProcessor3D (int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, const std::vector< T > &stoichiometricCoeff_, const std::vector< T > latticeReactionCoeff_, const std::vector< T > &react_order_, 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>
class olb::ConcentrationAdvectionDiffusionCouplingPostProcessor3D< T, DESCRIPTOR >

Coupling of ADlattice[0] with the other AD lattices (tpartners)

Definition at line 45 of file advectionDiffusionReactionCouplingPostProcessor3D.h.

Constructor & Destructor Documentation

ConcentrationAdvectionDiffusionCouplingPostProcessor3D()

template<typename T , typename DESCRIPTOR >

olb::ConcentrationAdvectionDiffusionCouplingPostProcessor3D< T, DESCRIPTOR >::ConcentrationAdvectionDiffusionCouplingPostProcessor3D ( int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, const std::vector< T > & stoichiometricCoeff_, const std::vector< T > latticeReactionCoeff_, const std::vector< T > & react_order_, std::vector< BlockStructureD< 3 > * > partners_ )

inline

Definition at line 49 of file advectionDiffusionReactionCouplingPostProcessor3D.h.

   : x0(x0_), x1(x1_), y0(y0_), y1(y1_), z0(z0_), z1(z1_), stoichiometricCoeff(stoichiometricCoeff_), latticeReactionCoeff(latticeReactionCoeff_),
     react_order(react_order_), partners(partners_) {
    this->getName() = "ConcentrationAdvectionDiffusionCouplingPostProcessor3D";
    reaction_number = static_cast<int>(latticeReactionCoeff.size());
    component_number = static_cast<int>(partners.size())+1;
    for (int i = 0; i<component_number;i++){
      tpartners.emplace_back(
        static_cast<BlockLattice<T,DESCRIPTOR> *>(partners[i]));
    }
  }

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 >

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

inlineoverridevirtual

Extent of application area (0 for purely local operations)

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 66 of file advectionDiffusionReactionCouplingPostProcessor3D.h.

                              {
    return 0;
  }

extent() [2/2]

template<typename T , typename DESCRIPTOR >

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

inlineoverridevirtual

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

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 70 of file advectionDiffusionReactionCouplingPostProcessor3D.h.

                                                {
    return 0;
  }

process()

template<typename T , typename DESCRIPTOR >

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

inlineoverridevirtual

Execute post-processing step.

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 74 of file advectionDiffusionReactionCouplingPostProcessor3D.h.

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

References olb::ConcentrationAdvectionDiffusionCouplingPostProcessor3D< T, DESCRIPTOR >::processSubDomain().

Here is the call graph for this function:

processSubDomain()

template<typename T , typename DESCRIPTOR >

void olb::ConcentrationAdvectionDiffusionCouplingPostProcessor3D< T, DESCRIPTOR >::processSubDomain ( BlockLattice< T, DESCRIPTOR > & blockLattice, int x0_, int x1_, int y0_, int y1_, int z0_, int z1_ )

inlineoverridevirtual

Execute post-processing step on a sublattice.

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 78 of file advectionDiffusionReactionCouplingPostProcessor3D.h.

                                                                                        {
 
    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) ) {
 
      for (int iX=newX0; iX<=newX1; ++iX) {
        for (int iY=newY0; iY<=newY1; ++iY) {
          for (int iZ=newZ0; iZ<=newZ1; ++iZ) {
 
            std::vector<T> conc;
            conc.emplace_back(blockLattice.get(iX,iY, iZ).computeRho());
            for (int iter_component = 0; iter_component<component_number-1; ++ iter_component){
              conc.emplace_back(tpartners[iter_component]->get(iX,iY,iZ).computeRho());
              }
 
            T lambda[reaction_number];
            T reaction_rate;
            for (int iter_reaction = 0; iter_reaction<reaction_number; ++ iter_reaction){
              lambda[iter_reaction] = 0;
              reaction_rate = 1;
              for(int iter_component = 0; iter_component <component_number; ++ iter_component){
 
                reaction_rate = reaction_rate*(util::pow(conc[iter_component],react_order[iter_reaction*component_number+iter_component]));
              }
 
              lambda[iter_reaction] = reaction_rate*latticeReactionCoeff[iter_reaction];
           }
            T temp_source;
            for (int iter_component = 0; iter_component<component_number; ++ iter_component){
              temp_source = 0;
              for (int iter_reaction = 0; iter_reaction<reaction_number; ++ iter_reaction){
                temp_source = temp_source + stoichiometricCoeff[iter_reaction*component_number+iter_component]*lambda[iter_reaction];
              }
              if (iter_component == 0){
                blockLattice.get(iX,iY,iZ).template setField<descriptors::SOURCE>(
                temp_source);
              }
              else {
                tpartners[iter_component-1]->get(iX,iY,iZ).template setField<descriptors::SOURCE>(
                temp_source);
                }
            }
         }
         }
      }
    }
  }

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

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• src/dynamics/advectionDiffusionReactionCouplingPostProcessor3D.h