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

Free Surface Processor 7 Finishes up left over cell conversions and prepares the state for the next simulation step. More...

#include <freeSurfacePostProcessor2D.h>

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

Public Member Functions
int	getPriority () const
template<typename CELL >
void	apply (CELL &cell) any_platform
	Free Surface Processor 7: Finalize Conversion.

Static Public Attributes
static constexpr OperatorScope	scope = OperatorScope::PerCell

Detailed Description

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

Free Surface Processor 7 Finishes up left over cell conversions and prepares the state for the next simulation step.

Definition at line 146 of file freeSurfacePostProcessor2D.h.

Member Function Documentation

apply()

template<typename T , typename DESCRIPTOR >

template<typename CELL >

void olb::FreeSurfaceFinalizeConversionPostProcessor2D< T, DESCRIPTOR >::apply

(

CELL &

cell

)

Free Surface Processor 7: Finalize Conversion.

Definition at line 481 of file freeSurfacePostProcessor2D.hh.

                                                                                  {
 
  using namespace olb::FreeSurface;
 
  // Convert flagged cells to appropriate cell types
  if (hasCellFlags(cell, FreeSurface::Flags::ToFluid)) {
    setCellType(cell, FreeSurface::Type::Fluid);
    cell.template setField<FreeSurface::EPSILON>(T(1));
  }
  else if (hasCellFlags(cell, FreeSurface::Flags::ToGas)) {
    setCellType(cell, FreeSurface::Type::Gas);
    cell.template setField<FreeSurface::EPSILON>(T(0));
  }
  else if (hasCellFlags(cell, FreeSurface::Flags::NewInterface)) {
    setCellType(cell, FreeSurface::Type::Interface);
  }
  else {
    // if (hasCellFlags(cell, FreeSurface::Flags::None) == true),
    // no further action is needed, as no conversion is triggered.
  }
 
  // Mass excess is distributed to old and new interface cells only,
  // thus skip the non-interface cells.
  if (!isCellType(cell, FreeSurface::Type::Interface)) {
    return;
  }
 
  // Collection of mass excess in a pulling step for the interface cells
  T collected_excess = T(0);
  for (int iPop = 1; iPop < DESCRIPTOR::q; ++iPop) {
    auto nbrCell = cell.neighbor(descriptors::c<DESCRIPTOR>(iPop));
 
    if (hasCellFlags(nbrCell, FreeSurface::Flags::ToFluid | FreeSurface::Flags::ToGas)) {
      auto tempMassExchange = nbrCell.template getField<FreeSurface::TEMP_MASS_EXCHANGE>();
      collected_excess += tempMassExchange[descriptors::opposite<DESCRIPTOR>(iPop)];
    }
  }
 
  T mass_tmp = cell.template getField<FreeSurface::MASS>();
  mass_tmp += collected_excess;
  T rho = T(0);
  T u_tmp[DESCRIPTOR::d] = {T(0), T(0)};
  cell.computeRhoU(rho, u_tmp);
  Vector<T,DESCRIPTOR::d> u_vel{u_tmp[0], u_tmp[1]};
  cell.template setField<FreeSurface::EPSILON>(mass_tmp / rho);
  cell.template setField<FreeSurface::MASS>(mass_tmp);
  cell.template setField<FreeSurface::PREVIOUS_VELOCITY>(u_vel);
}

References olb::descriptors::c(), olb::FreeSurface::Fluid, olb::FreeSurface::Gas, olb::FreeSurface::Interface, olb::FreeSurface::NewInterface, olb::descriptors::opposite(), olb::FreeSurface::ToFluid, and olb::FreeSurface::ToGas.

Here is the call graph for this function:

getPriority()

template<typename T , typename DESCRIPTOR >

int olb::FreeSurfaceFinalizeConversionPostProcessor2D< T, DESCRIPTOR >::getPriority ( ) const

inline

Definition at line 150 of file freeSurfacePostProcessor2D.h.

                          {
    return 6;
  }

Member Data Documentation

scope

template<typename T , typename DESCRIPTOR >

OperatorScope olb::FreeSurfaceFinalizeConversionPostProcessor2D< T, DESCRIPTOR >::scope = OperatorScope::PerCell

staticconstexpr

Definition at line 148 of file freeSurfacePostProcessor2D.h.

---

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

• src/dynamics/freeSurfacePostProcessor2D.h
• src/dynamics/freeSurfacePostProcessor2D.hh