Nektar::QInflow Class Reference

#include <QInflow.h>

Inheritance diagram for Nektar::QInflow:

Public Member Functions
	QInflow (Array< OneD, MultiRegions::ExpListSharedPtr > pVessel, const LibUtilities::SessionReaderSharedPtr pSession, PulseWavePressureAreaSharedPtr pressureArea)
	~QInflow () override
Public Member Functions inherited from Nektar::PulseWaveBoundary
	PulseWaveBoundary (Array< OneD, MultiRegions::ExpListSharedPtr > &pVessel, const LibUtilities::SessionReaderSharedPtr &pSession, PulseWavePressureAreaSharedPtr &pressureArea)
virtual	~PulseWaveBoundary ()
void	DoBoundary (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &A_0, Array< OneD, Array< OneD, NekDouble > > &beta, Array< OneD, Array< OneD, NekDouble > > &alpha, const NekDouble time, int omega, int offset, int n)

Static Public Member Functions
static PulseWaveBoundarySharedPtr	create (Array< OneD, MultiRegions::ExpListSharedPtr > &pVessel, const LibUtilities::SessionReaderSharedPtr &pSession, PulseWavePressureAreaSharedPtr &pressureArea)

Static Public Attributes
static std::string	className
	Name of class. More...

Protected Member Functions
void	v_DoBoundary (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &A_0, Array< OneD, Array< OneD, NekDouble > > &beta, Array< OneD, Array< OneD, NekDouble > > &alpha, const NekDouble time, int omega, int offset, int n) override
virtual void	v_DoBoundary (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &A_0, Array< OneD, Array< OneD, NekDouble > > &beta, Array< OneD, Array< OneD, NekDouble > > &alpha, const NekDouble time, int omega, int offset, int n)=0

Private Member Functions
void	Q_RiemannSolver (NekDouble Q, NekDouble A_r, NekDouble u_r, NekDouble A_0, NekDouble beta, NekDouble alpha, NekDouble &Au, NekDouble &uu)

Additional Inherited Members
Protected Attributes inherited from Nektar::PulseWaveBoundary
Array< OneD, MultiRegions::ExpListSharedPtr >	m_vessels
LibUtilities::SessionReaderSharedPtr	m_session
PulseWavePressureAreaSharedPtr	m_pressureArea
NekDouble	m_pext
NekDouble	m_pout
NekDouble	m_rho

Detailed Description

Definition at line 52 of file QInflow.h.

Constructor & Destructor Documentation

QInflow()

Nektar::QInflow::QInflow	(	Array< OneD, MultiRegions::ExpListSharedPtr >	pVessel,
		const LibUtilities::SessionReaderSharedPtr	pSession,
		PulseWavePressureAreaSharedPtr	pressureArea
	)

Definition at line 47 of file QInflow.cpp.

    : PulseWaveBoundary(pVessel, pSession, pressureArea)
{
}

~QInflow()

Nektar::QInflow::~QInflow ( )

override

Definition at line 54 of file QInflow.cpp.

{
}

Member Function Documentation

create()

static PulseWaveBoundarySharedPtr Nektar::QInflow::create ( Array< OneD, MultiRegions::ExpListSharedPtr > &  pVessel, const LibUtilities::SessionReaderSharedPtr &  pSession, PulseWavePressureAreaSharedPtr &  pressureArea  )

inlinestatic

Definition at line 56 of file QInflow.h.

    {
        return MemoryManager<QInflow>::AllocateSharedPtr(pVessel, pSession,
                                                         pressureArea);
    }

References Nektar::MemoryManager< DataType >::AllocateSharedPtr().

Q_RiemannSolver()

void Nektar::QInflow::Q_RiemannSolver ( NekDouble  Q, NekDouble  A_r, NekDouble  u_r, NekDouble  A_0, NekDouble  beta, NekDouble  alpha, NekDouble &  Au, NekDouble &  uu  )

private

Q-inflow Riemann solver for pulse wave propagation. This Riemann solver is called by DoOdeProjection in case of a Q-inflow boundary condition. It is based on the conservation of mass and total pressure and on the characteristic information. For further details see "Pulse wave propagation in the human vascular system", section 3.4.1 Returns the upwinded quantities \((A_u,u_u)\) and stores them into the boundary values

Definition at line 111 of file QInflow.cpp.

{
    NekDouble W2           = 0.0;
    NekDouble A_calc       = 0.0;
    NekDouble FA           = 0.0;
    NekDouble dFdA         = 0.0;
    NekDouble delta_A_calc = 0.0;
    NekDouble I            = 0.0;
    NekDouble c            = 0.0;
 
    int proceed  = 1;
    int iter     = 0;
    int MAX_ITER = 200;
 
    // Tolerances for the algorithm
    NekDouble Tol = 1.0e-10;
 
    // Riemann invariant \f$W_2(Ar,ur)\f$
    m_pressureArea->GetW2(W2, u_r, beta, A_r, A_0, alpha);
 
    // Newton Iteration (Area only)
    A_calc = A_r;
    while ((proceed) && (iter < MAX_ITER))
    {
        iter += 1;
 
        m_pressureArea->GetCharIntegral(I, beta, A_calc, A_0, alpha);
        m_pressureArea->GetC(c, beta, A_calc, A_0, alpha);
 
        FA           = Q - A_calc * (W2 + I);
        dFdA         = -W2 - I - c;
        delta_A_calc = FA / dFdA;
        A_calc -= delta_A_calc;
 
        if (sqrt(delta_A_calc * delta_A_calc) < Tol)
        {
            proceed = 0;
        }
    }
 
    m_pressureArea->GetCharIntegral(I, beta, A_calc, A_0, alpha);
 
    // Obtain u_u and A_u
    uu = W2 + I;
    Au = A_calc;
}

References Nektar::LibUtilities::beta, Nektar::PulseWaveBoundary::m_pressureArea, and tinysimd::sqrt().

Referenced by v_DoBoundary().

v_DoBoundary()

void Nektar::QInflow::v_DoBoundary ( const Array< OneD, const Array< OneD, NekDouble > > &  inarray, Array< OneD, Array< OneD, NekDouble > > &  A_0, Array< OneD, Array< OneD, NekDouble > > &  beta, Array< OneD, Array< OneD, NekDouble > > &  alpha, const NekDouble  time, int  omega, int  offset, int  n  )

overrideprotectedvirtual

Implements Nektar::PulseWaveBoundary.

Definition at line 58 of file QInflow.cpp.

{
    NekDouble Q;
    NekDouble A_u;
    NekDouble u_u;
    NekDouble A_r;
    NekDouble u_r;
    NekDouble A_l;
    NekDouble u_l;
 
    Array<OneD, MultiRegions::ExpListSharedPtr> vessel(2);
 
    // Pointers to the domains
    vessel[0] = m_vessels[2 * omega];
    vessel[1] = m_vessels[2 * omega + 1];
 
    // Evaluates the time-dependent Q
    vessel[0]->EvaluateBoundaryConditions(time);
 
    // Q is contained as A in the input file
    Q = (vessel[0]->UpdateBndCondExpansion(n))->GetCoeffs()[0];
 
    // Initial conditions in the inlet vessel
    A_r = inarray[0][offset];
    u_r = inarray[1][offset];
 
    // Call the Q-inflow Riemann solver
    Q_RiemannSolver(Q, A_r, u_r, A_0[omega][0], beta[omega][0], alpha[omega][0],
                    A_u, u_u);
 
    /* Fix the boundary conditions in the virtual region to ensure
       upwind state matches the boundary condition at the next time step */
    A_l = A_r;
    u_l = 2 * u_u - u_r;
 
    // Store the updated values in the boundary condition
    (vessel[0]->UpdateBndCondExpansion(n))->UpdatePhys()[0] = A_l;
    (vessel[1]->UpdateBndCondExpansion(n))->UpdatePhys()[0] = u_l;
}

References Nektar::LibUtilities::beta, Nektar::PulseWaveBoundary::m_vessels, and Q_RiemannSolver().

Member Data Documentation

className

std::string Nektar::QInflow::className

static

Initial value:

= GetBoundaryFactory().RegisterCreatorFunction(
    "Q-inflow", QInflow::create, "Inflow boundary condition")

Name of class.

Definition at line 66 of file QInflow.h.