A global linear system. More...

#include <TerminalOutflow.h>

Inheritance diagram for Nektar::TerminalOutflow:

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

Protected Member Functions
	TerminalOutflow (Array< OneD, MultiRegions::ExpListSharedPtr > pVessel, const LibUtilities::SessionReaderSharedPtr pSession, PulseWavePressureAreaSharedPtr pressureArea)

	~TerminalOutflow () override=default

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

Friends
class	MemoryManager< TerminalOutflow >

Additional Inherited Members
Public Member Functions inherited from Nektar::PulseWaveBoundary
	PulseWaveBoundary (Array< OneD, MultiRegions::ExpListSharedPtr > &pVessel, const LibUtilities::SessionReaderSharedPtr &pSession, PulseWavePressureAreaSharedPtr &pressureArea)

virtual	~PulseWaveBoundary ()=default

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)

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

A global linear system.

Definition at line 52 of file TerminalOutflow.h.

Constructor & Destructor Documentation

◆ TerminalOutflow()

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

protected

Definition at line 46 of file TerminalOutflow.cpp.

    : PulseWaveBoundary(pVessel, pSession, pressureArea)
{
}

◆ ~TerminalOutflow()

Nektar::TerminalOutflow::~TerminalOutflow ( )

overrideprotecteddefault

Member Function Documentation

◆ create()

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

inlinestatic

Definition at line 58 of file TerminalOutflow.h.

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

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

◆ v_DoBoundary()

void Nektar::TerminalOutflow::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 54 of file TerminalOutflow.cpp.

{
    NekDouble A_r;
    NekDouble A_l;
    NekDouble u_r;
    NekDouble u_l;
    NekDouble RT;
    NekDouble W1;
 
    Array<OneD, MultiRegions::ExpListSharedPtr> vessel(2);
 
    // Pointers to the domains
    vessel[0] = m_vessels[2 * omega];
    vessel[1] = m_vessels[2 * omega + 1];
 
    /* Find the terminal resistance boundary condition and
     * calculate the reflection. We assume A_r = A_l and
     * apply the reflection in u_r after Jordi Alastruey's
     * PhD thesis. Stems from R = -W2 / W1.
     */
 
    // Note: The R_t value is contained in A in the inputfile
    RT = (vessel[0]->UpdateBndCondExpansion(n))->GetCoeffs()[0];
 
    ASSERTL0((-1 <= RT && RT <= 1), "RT must be between -1 and 1");
    int nq = vessel[0]->GetTotPoints();
 
    // Get the left values A and u, and the forward characteristic
    A_l = inarray[0][offset + nq - 1];
    u_l = inarray[1][offset + nq - 1];
 
    m_pressureArea->GetW1(W1, u_l, beta[omega][nq - 1], A_l, A_0[omega][nq - 1],
                          alpha[omega][nq - 1]);
 
    // Calculate the boundary values
    A_r = A_l;
    u_r = (1 - RT) * W1 - u_l;
 
    // Store the new values in the boundary condition
    (vessel[0]->UpdateBndCondExpansion(n))->UpdatePhys()[0] = A_r;
    (vessel[1]->UpdateBndCondExpansion(n))->UpdatePhys()[0] = u_r;
}

References ASSERTL0, Nektar::LibUtilities::beta, Nektar::PulseWaveBoundary::m_pressureArea, and Nektar::PulseWaveBoundary::m_vessels.

Friends And Related Function Documentation

◆ MemoryManager< TerminalOutflow >

friend class MemoryManager< TerminalOutflow >

friend

Definition at line 49 of file TerminalOutflow.h.

Member Data Documentation

◆ className

std::string Nektar::TerminalOutflow::className

static

Initial value:

=
    GetBoundaryFactory().RegisterCreatorFunction(
        "Terminal", TerminalOutflow::create,
        "Terminal outflow boundary condition")

Definition at line 68 of file TerminalOutflow.h.

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Friends

Additional Inherited Members