Nektar::UInflow Class Reference

A global linear system. More...

#include <UInflow.h>

Inheritance diagram for Nektar::UInflow:

Public Member Functions
	UInflow (Array< OneD, MultiRegions::ExpListSharedPtr > pVessel, const LibUtilities::SessionReaderSharedPtr pSession, PulseWavePressureAreaSharedPtr pressureArea)
virtual	~UInflow ()
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, 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)
	Creates an instance of this class. More...

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

Protected Member Functions
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, const NekDouble time, int omega, int offset, int n)
Protected Member Functions inherited from Nektar::PulseWaveBoundary
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, const NekDouble time, int omega, int offset, int n)=0

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

A global linear system.

Definition at line 50 of file UInflow.h.

Constructor & Destructor Documentation

UInflow()

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

Definition at line 47 of file UInflow.cpp.

    : PulseWaveBoundary(pVessel, pSession, pressureArea)
{
    // Constructor
}

~UInflow()

Nektar::UInflow::~UInflow ( )

virtual

Definition at line 58 of file UInflow.cpp.

{
    // Destructor
}

Member Function Documentation

create()

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

inlinestatic

Creates an instance of this class.

Definition at line 54 of file UInflow.h.

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

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

v_DoBoundary()

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

protectedvirtual

Definition at line 63 of file UInflow.cpp.

References Nektar::PulseWaveBoundary::m_vessels.

{
    NekDouble 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 U
    vessel[1]->EvaluateBoundaryConditions(time);

    // Read the BC values from the input file
    u   = (vessel[1]->UpdateBndCondExpansion(n))->GetCoeffs()[0];

    // Initial conditions in inlet vessel
    A_r = inarray[0][offset];
    u_r = inarray[1][offset];

    /* Fix the boundary conditions in the virtual vessel to ensure
       upwind state matches the boundary condition at the next time step */
    A_l = A_r;
    u_l = 2 * 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;
}

Member Data Documentation

className

std::string Nektar::UInflow::className

static

Initial value:

= GetBoundaryFactory().RegisterCreatorFunction(
    "U-inflow", UInflow::create, "Velocity inflow boundary condition")

Name of class.

Definition at line 64 of file UInflow.h.