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

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

#include <advectionDiffusionReactionCouplingPostProcessor2D.h>

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

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

Public Member Functions
	ConcentrationAdvectionDiffusionCouplingPostProcessor2D (int x0_, int x1_, int y0_, int y1_, const std::vector< T > &stochiometricCoeff_, const std::vector< T > latticeReactionCoeff_, const std::vector< T > &react_order_, std::vector< BlockStructureD< 2 > * > 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_) override
	Execute post-processing step on a sublattice.
void	computeSources (T sources[], const T concentrations[])
Public Member Functions inherited from olb::PostProcessor2D< T, DESCRIPTOR >
	PostProcessor2D ()
virtual	~PostProcessor2D ()
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::PostProcessor2D< T, DESCRIPTOR >
int	_priority

Detailed Description

template<typename T, typename DESCRIPTOR>
class olb::ConcentrationAdvectionDiffusionCouplingPostProcessor2D< T, DESCRIPTOR >

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

Definition at line 45 of file advectionDiffusionReactionCouplingPostProcessor2D.h.

Constructor & Destructor Documentation

ConcentrationAdvectionDiffusionCouplingPostProcessor2D()

template<typename T , typename DESCRIPTOR >

olb::ConcentrationAdvectionDiffusionCouplingPostProcessor2D< T, DESCRIPTOR >::ConcentrationAdvectionDiffusionCouplingPostProcessor2D ( int x0_, int x1_, int y0_, int y1_, const std::vector< T > & stochiometricCoeff_, const std::vector< T > latticeReactionCoeff_, const std::vector< T > & react_order_, std::vector< BlockStructureD< 2 > * > partners_ )

inline

Definition at line 49 of file advectionDiffusionReactionCouplingPostProcessor2D.h.

   : x0(x0_), x1(x1_), y0(y0_), y1(y1_),
     stochiometricCoeff(stochiometricCoeff_),
     latticeReactionCoeff(latticeReactionCoeff_),
     react_order(react_order_), partners(partners_) {
    this->getName() = "ConcentrationAdvectionDiffusionCouplingPostProcessor2D";
    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::PostProcessor2D< T, DESCRIPTOR >::getName().

Here is the call graph for this function:

Member Function Documentation

computeSources()

template<typename T , typename DESCRIPTOR >

void olb::ConcentrationAdvectionDiffusionCouplingPostProcessor2D< T, DESCRIPTOR >::computeSources ( T sources[], const T concentrations[] )

inline

Definition at line 113 of file advectionDiffusionReactionCouplingPostProcessor2D.h.

                                                             {
    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 *= util::pow(concentrations[iter_component],
          react_order[iter_reaction*component_number+iter_component]);
      }
      lambda[iter_reaction] = reaction_rate*latticeReactionCoeff[iter_reaction];
    }
 
    for (int iter_component = 0; iter_component<component_number; ++iter_component) {
      sources[iter_component] = 0;
      for (int iter_reaction = 0; iter_reaction<reaction_number; ++iter_reaction) {
        sources[iter_component]
          += stochiometricCoeff[iter_reaction*component_number+iter_component]*lambda[iter_reaction];
      }
    }
  }

References olb::util::pow().

Here is the call graph for this function:

Here is the caller graph for this function:

extent() [1/2]

template<typename T , typename DESCRIPTOR >

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

inlineoverridevirtual

Extent of application area (0 for purely local operations)

Implements olb::PostProcessor2D< T, DESCRIPTOR >.

Definition at line 68 of file advectionDiffusionReactionCouplingPostProcessor2D.h.

                              {
    return 0;
  }

extent() [2/2]

template<typename T , typename DESCRIPTOR >

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

inlineoverridevirtual

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

Implements olb::PostProcessor2D< T, DESCRIPTOR >.

Definition at line 72 of file advectionDiffusionReactionCouplingPostProcessor2D.h.

                                                {
    return 0;
  }

process()

template<typename T , typename DESCRIPTOR >

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

inlineoverridevirtual

Execute post-processing step.

Implements olb::PostProcessor2D< T, DESCRIPTOR >.

Definition at line 76 of file advectionDiffusionReactionCouplingPostProcessor2D.h.

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

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

Here is the call graph for this function:

processSubDomain()

template<typename T , typename DESCRIPTOR >

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

inlineoverridevirtual

Execute post-processing step on a sublattice.

Implements olb::PostProcessor2D< T, DESCRIPTOR >.

Definition at line 80 of file advectionDiffusionReactionCouplingPostProcessor2D.h.

                                                                     {
 
    int newX0, newX1, newY0, newY1;
    if ( util::intersect (
           x0, x1, y0, y1,
           x0_, x1_, y0_, y1_,
           newX0, newX1, newY0, newY1 ) ) {
 
      for (int iX=newX0; iX<=newX1; ++iX) {
        for (int iY=newY0; iY<=newY1; ++iY) {
 
          //std::cout << "iX: " << iX << "     iY: " << iY << std::endl;
 
          T conc[component_number];
          conc[0] = blockLattice.get(iX,iY).computeRho();
          for (int iter_component = 1; iter_component<component_number; ++iter_component) {
            conc[iter_component] = tpartners[iter_component-1]->get(iX,iY).computeRho();
          }
 
          T sources[component_number];
          computeSources(sources, conc);
 
          blockLattice.get(iX,iY).template setField<descriptors::SOURCE>(sources[0]);
          for (int iter_component = 1; iter_component<component_number; ++iter_component) {
            tpartners[iter_component-1]->get(iX,iY).template setField<descriptors::SOURCE>(
              sources[iter_component]);
          }
        }
      }
    }
  }

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

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/advectionDiffusionReactionCouplingPostProcessor2D.h