TotalEnthalpyPhaseChange between a Navier-Stokes and an Advection-Diffusion lattice. More...

#include <navierStokesAdvectionDiffusionCoupling.h>

Collaboration diagram for olb::TotalEnthalpyPhaseChangeCoupling:

Classes
struct	CP_L

struct	CP_S

struct	FORCE_PREFACTOR

struct	L

struct	T0

struct	T_L

struct	T_S

Public Types
using	parameters = meta::list<FORCE_PREFACTOR,T0, CP_S, CP_L, T_S, T_L, L>

Public Member Functions
template<typename CELLS , typename PARAMETERS >
void	apply (CELLS &cells, PARAMETERS &parameters) any_platform

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

Detailed Description

TotalEnthalpyPhaseChange between a Navier-Stokes and an Advection-Diffusion lattice.

Definition at line 32 of file navierStokesAdvectionDiffusionCoupling.h.

Member Typedef Documentation

◆ parameters

using olb::TotalEnthalpyPhaseChangeCoupling::parameters = meta::list<FORCE_PREFACTOR,T0, CP_S, CP_L, T_S, T_L, L>

Definition at line 47 of file navierStokesAdvectionDiffusionCoupling.h.

Member Function Documentation

◆ apply()

template<typename CELLS , typename PARAMETERS >

void olb::TotalEnthalpyPhaseChangeCoupling::apply	(	CELLS &	cells,
		PARAMETERS &	parameters )

inline

Definition at line 50 of file navierStokesAdvectionDiffusionCoupling.h.

  {
    using V = typename CELLS::template value_t<names::NavierStokes>::value_t;
    using DESCRIPTOR = typename CELLS::template value_t<names::NavierStokes>::descriptor_t;
 
    auto& cellNSE = cells.template get<names::NavierStokes>();
    auto& cellADE = cells.template get<names::Temperature>();
 
    V enthalpy = cellADE.computeRho();
    auto& dynamics = cellADE.getDynamics();
    /*auto& dynParams = static_cast<ParametersOfDynamicsD<DYNAMICS>&>(
      tPartner->template getData<OperatorParameters<DYNAMICS>>());*/
 
    //cellNSE.template setField<descriptors::POROSITY>(
      //dynamics->template computeLiquidFraction<T>(parameters, enthalpy)
    //);
 
    auto temperature = cellADE.template getFieldPointer<descriptors::TEMPERATURE>();
    //temperature[0] = dynamics->computeTemperature<V, parameters, enthalpy>();
 
    // computation of the bousinessq force
    auto force = cellNSE.template getFieldPointer<descriptors::FORCE>();
    auto forcePrefactor = parameters.template get<FORCE_PREFACTOR>();
    V temperatureDifference = cellADE.computeRho() - parameters.template get<T0>();
    for (unsigned iD = 0; iD < DESCRIPTOR::d; ++iD) {
      force[iD] = forcePrefactor[iD] * temperatureDifference;
    }
    // Velocity coupling
    V u[DESCRIPTOR::d] { };
    cellNSE.computeU(u);
    cellADE.template setField<descriptors::VELOCITY>(u);
  }

Member Data Documentation

◆ scope

constexpr OperatorScope olb::TotalEnthalpyPhaseChangeCoupling::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 33 of file navierStokesAdvectionDiffusionCoupling.h.

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

src/dynamics/navierStokesAdvectionDiffusionCoupling.h

Classes

Public Types

Public Member Functions