43 std::string TerminalOutflow::className
46 TerminalOutflow::create,
47 "Terminal outflow boundary condition");
72 int omega,
int offset,
int n)
88 RT = (vessel[0]->UpdateBndCondExpansion(n))->GetCoeffs()[0];
91 "RT must be comprised between -1 and 1");
92 int nq = vessel[0]->GetTotPoints();
95 A_l = inarray[0][offset+nq-1];
96 u_l = inarray[1][offset+nq-1];
100 c_0 = sqrt(beta[omega][nq-1]/(2*
m_rho))*sqrt(sqrt(A_0[omega][nq-1]));
104 c_l = sqrt(beta[omega][nq-1]/(2*
m_rho))*sqrt(sqrt(A_l));
105 u_r = (1-RT)*((u_l-u_0) + 4*(c_l-c_0)) - u_l;
108 (vessel[0]->UpdateBndCondExpansion(n))->UpdatePhys()[0] = A_r;
109 (vessel[1]->UpdateBndCondExpansion(n))->UpdatePhys()[0] = u_r;
#define ASSERTL0(condition, msg)
Array< OneD, MultiRegions::ExpListSharedPtr > m_vessels
boost::shared_ptr< SessionReader > SessionReaderSharedPtr
BoundaryFactory & GetBoundaryFactory()
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)
virtual ~TerminalOutflow()
boost::shared_ptr< PulseWavePressureArea > PulseWavePressureAreaSharedPtr
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, tDescription pDesc="")
Register a class with the factory.