43 std::string QInflow::className
47 "Inflow boundary condition");
83 vessel[0]->EvaluateBoundaryConditions(time);
90 Q = (vessel[0]->UpdateBndCondExpansion(n))->GetCoeffs()[0];
92 A_r = inarray[0][offset];
93 u_r = inarray[1][offset];
103 (vessel[0]->UpdateBndCondExpansion(n))->UpdatePhys()[0] = A_l;
104 (vessel[1]->UpdateBndCondExpansion(n))->UpdatePhys()[0] = u_l;
132 W2 = u_r - 4*sqrt(beta/(2*rho))*sqrt(sqrt(A_r));
136 while ((proceed) && (iter < MAX_ITER))
140 fa = Q - W2*A_calc - A_calc*4*sqrt(beta/(2*rho))*sqrt(sqrt(A_calc));
141 dfa = -W2 - 5*sqrt(beta/(2*rho))*sqrt(sqrt(A_calc));
142 delta_A_calc = fa/dfa;
143 A_calc = A_calc - delta_A_calc;
145 if (sqrt(delta_A_calc*delta_A_calc) < Tol)
150 uu = W2+4*sqrt(beta/(2*rho))*sqrt(sqrt(A_calc));
void Q_RiemannSolver(NekDouble Q, NekDouble A_r, NekDouble u_r, NekDouble A_0, NekDouble beta, NekDouble &Au, NekDouble &uu)
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)
boost::shared_ptr< PulseWavePressureArea > PulseWavePressureAreaSharedPtr
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, tDescription pDesc="")
Register a class with the factory.