43 std::string ROutflow::className
47 "Resistive outflow boundary condition");
72 int omega,
int offset,
int n)
87 NekDouble RT=((vessel[0]->GetBndCondExpansions())[n])->GetCoeffs()[0];
89 int nq = vessel[0]->GetTotPoints();
92 A_l = inarray[0][offset+nq-1];
93 u_l = inarray[1][offset+nq-1];
97 beta[omega][nq-1],pout,A_u,u_u);
105 (vessel[0]->UpdateBndCondExpansion(n))->UpdatePhys()[0] = A_r;
106 (vessel[1]->UpdateBndCondExpansion(n))->UpdatePhys()[0] = u_r;
130 c_l = sqrt(beta/(2*rho))*sqrt(sqrt(A_l));
137 while ((proceed) && (iter < MAX_ITER))
141 fa = R*W1*A_calc-4*R*sqrt(beta/(2*rho))*A_calc*sqrt(sqrt(A_calc))-pext-beta*(sqrt(A_calc)-sqrt(A_0))+pout;
142 dfa = R*W1-5*R*sqrt(beta/(2*rho))*sqrt(sqrt(A_calc))-beta/(2*sqrt(A_calc));
143 delta_A_calc = fa/dfa;
144 A_calc = A_calc - delta_A_calc;
146 if (sqrt(delta_A_calc*delta_A_calc) < Tol)
151 u_u=(pext+beta*(sqrt(A_calc)-sqrt(A_0))-pout)/(R*A_calc);
Array< OneD, MultiRegions::ExpListSharedPtr > m_vessels
void R_RiemannSolver(NekDouble R, NekDouble A_l, NekDouble u_l, NekDouble A_0, NekDouble beta, NekDouble pout, NekDouble &A_u, NekDouble &u_u)
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)
boost::shared_ptr< PulseWavePressureArea > PulseWavePressureAreaSharedPtr
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, tDescription pDesc="")
Register a class with the factory.