ProcessLib::TES::TESFEMReactionAdaptorCaOH2 Class Reference

Detailed Description

Definition at line 109 of file TESReactionAdaptor.h.

#include <TESReactionAdaptor.h>

Inheritance diagram for ProcessLib::TES::TESFEMReactionAdaptorCaOH2:

Collaboration diagram for ProcessLib::TES::TESFEMReactionAdaptorCaOH2:

Public Member Functions
	TESFEMReactionAdaptorCaOH2 (TESLocalAssemblerData const &data)
ReactionRate	initReaction (const unsigned) override
Public Member Functions inherited from ProcessLib::TES::TESFEMReactionAdaptor
virtual bool	checkBounds (std::vector< double > const &, std::vector< double > const &)
virtual void	preZerothTryAssemble ()
virtual double	getReactionDampingFactor () const
virtual	~TESFEMReactionAdaptor ()=default

Private Types
using	Data = TESLocalAssemblerData
using	React = Adsorption::ReactionCaOH2

Static Private Member Functions
static bool	odeRhs (const double, MathLib::ODE::MappedConstVector< 1 > const y, MathLib::ODE::MappedVector< 1 > ydot, React &reaction)

Private Attributes
Data const &	_d
React &	_react
std::unique_ptr< MathLib::ODE::ODESolver< 1 > >	_ode_solver

Additional Inherited Members
Static Public Member Functions inherited from ProcessLib::TES::TESFEMReactionAdaptor
static std::unique_ptr< TESFEMReactionAdaptor >	newInstance (TESLocalAssemblerData const &data)

Member Typedef Documentation

Data

using ProcessLib::TES::TESFEMReactionAdaptorCaOH2::Data = TESLocalAssemblerData

private

Definition at line 117 of file TESReactionAdaptor.h.

React

using ProcessLib::TES::TESFEMReactionAdaptorCaOH2::React = Adsorption::ReactionCaOH2

private

Definition at line 118 of file TESReactionAdaptor.h.

Constructor & Destructor Documentation

TESFEMReactionAdaptorCaOH2()

ProcessLib::TES::TESFEMReactionAdaptorCaOH2::TESFEMReactionAdaptorCaOH2 ( TESLocalAssemblerData const & data )

explicit

Definition at line 301 of file TESReactionAdaptor.cpp.

    : _d(data),
      _react(dynamic_cast<Adsorption::ReactionCaOH2&>(*data.ap.react_sys))
{
    _ode_solver = MathLib::ODE::createODESolver<1>(_react.getOdeSolverConfig());
    // TODO invalidate config
 
    _ode_solver->setTolerance(1e-10, 1e-10);
 
    auto f = [this](const double /*t*/,
                    MathLib::ODE::MappedConstVector<1> const y,
                    MathLib::ODE::MappedVector<1>
                        ydot) -> bool
    {
        ydot[0] = _react.getReactionRate(y[0]);
        return true;
    };
 
    _ode_solver->setFunction(f, nullptr);
}

References _ode_solver, _react, MathLib::ODE::createODESolver(), Adsorption::ReactionCaOH2::getOdeSolverConfig(), and Adsorption::ReactionCaOH2::getReactionRate().

Member Function Documentation

initReaction()

ReactionRate ProcessLib::TES::TESFEMReactionAdaptorCaOH2::initReaction ( const unsigned )

overridevirtual

Implements ProcessLib::TES::TESFEMReactionAdaptor.

Definition at line 323 of file TESReactionAdaptor.cpp.

{
    // TODO if the first holds, the second also has to hold
    if (_d.ap.iteration_in_current_timestep > 1 ||
        _d.ap.number_of_try_of_iteration > 1)
    {
        return {_d.reaction_rate[int_pt], _d.solid_density[int_pt]};
    }
 
    // TODO: double check!
    // const double xv_NR  = SolidProp->non_reactive_solid_volume_fraction;
    // const double rho_NR = SolidProp->non_reactive_solid_density;
    const double xv_NR = 0.0;
    const double rho_NR = 0.0;
 
    const double t0 = 0.0;
    const double y0 =
        (_d.solid_density_prev_ts[int_pt] - xv_NR * rho_NR) / (1.0 - xv_NR);
 
    const double t_end = _d.ap.delta_t;
 
    _react.updateParam(_d.T, _d.p, _d.vapour_mass_fraction,
                       _d.solid_density_prev_ts[int_pt]);
 
    _ode_solver->setIC(t0, {y0});
    _ode_solver->preSolve();
    _ode_solver->solve(t_end);
 
    const double time_reached = _ode_solver->getTime();
    (void)time_reached;
    assert(std::abs(t_end - time_reached) <
           std::numeric_limits<double>::epsilon());
 
    auto const& y_new = _ode_solver->getSolution();
    auto const& y_dot_new = _ode_solver->getYDot(t_end, y_new);
 
    double rho_react;
 
    // cut off when limits are reached
    if (y_new[0] < ProcessLib::TES::TESFEMReactionAdaptorCaOH2::React::rho_low)
    {
        rho_react = ProcessLib::TES::TESFEMReactionAdaptorCaOH2::React::rho_low;
    }
    else if (y_new[0] >
             ProcessLib::TES::TESFEMReactionAdaptorCaOH2::React::rho_up)
    {
        rho_react = ProcessLib::TES::TESFEMReactionAdaptorCaOH2::React::rho_up;
    }
    else
    {
        rho_react = y_new[0];
    }
 
    return {y_dot_new[0] * (1.0 - xv_NR),
            (1.0 - xv_NR) * rho_react + xv_NR * rho_NR};
}

References _d, _ode_solver, _react, ProcessLib::TES::TESLocalAssemblerData::ap, ProcessLib::TES::AssemblyParams::delta_t, ProcessLib::TES::AssemblyParams::iteration_in_current_timestep, ProcessLib::TES::AssemblyParams::number_of_try_of_iteration, ProcessLib::TES::TESLocalAssemblerData::p, ProcessLib::TES::TESLocalAssemblerData::reaction_rate, Adsorption::ReactionCaOH2::rho_low, Adsorption::ReactionCaOH2::rho_up, ProcessLib::TES::TESLocalAssemblerData::solid_density, ProcessLib::TES::TESLocalAssemblerData::solid_density_prev_ts, ProcessLib::TES::TESLocalAssemblerData::T, Adsorption::ReactionCaOH2::updateParam(), and ProcessLib::TES::TESLocalAssemblerData::vapour_mass_fraction.

odeRhs()

static bool ProcessLib::TES::TESFEMReactionAdaptorCaOH2::odeRhs ( const double , MathLib::ODE::MappedConstVector< 1 > const y, MathLib::ODE::MappedVector< 1 > ydot, React & reaction )

staticprivate

Member Data Documentation

_d

Data const& ProcessLib::TES::TESFEMReactionAdaptorCaOH2::_d

private

Definition at line 119 of file TESReactionAdaptor.h.

Referenced by initReaction().

_ode_solver

std::unique_ptr<MathLib::ODE::ODESolver<1> > ProcessLib::TES::TESFEMReactionAdaptorCaOH2::_ode_solver

private

Definition at line 122 of file TESReactionAdaptor.h.

Referenced by TESFEMReactionAdaptorCaOH2(), and initReaction().

_react

React& ProcessLib::TES::TESFEMReactionAdaptorCaOH2::_react

private

Definition at line 120 of file TESReactionAdaptor.h.

Referenced by TESFEMReactionAdaptorCaOH2(), and initReaction().

---

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

• ProcessLib/TES/TESReactionAdaptor.h
• ProcessLib/TES/TESReactionAdaptor.cpp