#include <AInflow.h>

Inheritance diagram for Nektar::AInflow:

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

virtual	~AInflow ()

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

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, Array< OneD, Array< OneD, NekDouble > > &alpha, 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, Array< OneD, Array< OneD, NekDouble > > &alpha, 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

Definition at line 51 of file AInflow.h.

Constructor & Destructor Documentation

◆ AInflow()

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

Definition at line 44 of file AInflow.cpp.

     : PulseWaveBoundary(pVessel, pSession, pressureArea)
 {
 }

◆ ~AInflow()

Nektar::AInflow::~AInflow ( )

virtual

Definition at line 51 of file AInflow.cpp.

52 {

53 }

Member Function Documentation

◆ create()

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

inlinestatic

Definition at line 55 of file AInflow.h.

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

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

◆ v_DoBoundary()

void Nektar::AInflow::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
	)

protectedvirtual

Definition at line 55 of file AInflow.cpp.

 {
     NekDouble A;
     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 A
     vessel[0]->EvaluateBoundaryConditions(time);
  
     // Read the BC values from the input file
     A = (vessel[0]->UpdateBndCondExpansion(n))->GetCoeffs()[0];
  
     // Initial conditions in the 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 */
     u_l = u_r;
     // Needs to be updated to work with other tube laws
     A_l = pow(2 * sqrt(sqrt(A)) - sqrt(sqrt(A_r)), 4);
  
     // 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::PulseWaveBoundary::m_vessels, and tinysimd::sqrt().

Member Data Documentation

◆ className

std::string Nektar::AInflow::className

static

Initial value:

= GetBoundaryFactory().RegisterCreatorFunction(

"A-inflow", AInflow::create, "Area inflow boundary condition")

Nektar::AInflow::create

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

Definition: AInflow.h:55

Nektar::LibUtilities::NekFactory::RegisterCreatorFunction

tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")

Register a class with the factory.

Definition: NekFactory.hpp:200

Nektar::GetBoundaryFactory

BoundaryFactory & GetBoundaryFactory()

Definition: PulseWaveBoundary.cpp:62

Definition at line 65 of file AInflow.h.