43 "Power law pressure area relationship for the arterial system");
67 (2 *
m_rho * c0 * c0 / b) *
68 (pow((
A / A0), b / 2) - 1)
69 - A0 * gamma * dAUdx / (
A *
sqrt(
A));
82 c = c0 * pow((
A / A0), b / 4);
117 A = A0 * pow(((b / (8 * c0)) * (W1 - W2) + 1), 4 / b);
139 I = (4 / b) * (c - c0);
146 const std::string &type)
149 if (type ==
"Bifurcation")
154 for (
int i = 0; i < 3; ++i)
156 GetC(c[i],
beta[i], Au[i], A0[i], alpha[i]);
162 J.SetValue(0, 3, c[0] / Au[0]);
170 J.SetValue(1, 4, -c[1] / Au[1]);
178 J.SetValue(2, 5, -c[2] / Au[2]);
180 J.SetValue(3, 0, Au[0]);
181 J.SetValue(3, 1, -Au[1]);
182 J.SetValue(3, 2, -Au[2]);
183 J.SetValue(3, 3, uu[0]);
184 J.SetValue(3, 4, -uu[1]);
185 J.SetValue(3, 5, -uu[2]);
187 J.SetValue(4, 0, 2 * uu[0]);
188 J.SetValue(4, 1, -2 * uu[1]);
190 J.SetValue(4, 3, 2 * c[0] * c[0] / Au[0]);
191 J.SetValue(4, 4, -2 * c[1] * c[1] / Au[1]);
194 J.SetValue(5, 0, 2 * uu[0]);
196 J.SetValue(5, 2, -2 * uu[2]);
197 J.SetValue(5, 3, 2 * c[0] * c[0] / Au[0]);
199 J.SetValue(5, 5, -2 * c[2] * c[2] / Au[2]);
204 else if (type ==
"Merge")
209 for (
int i = 0; i < 3; ++i)
211 GetC(c[i],
beta[i], Au[i], A0[i], alpha[i]);
217 J.SetValue(0, 3, -c[0] / Au[0]);
225 J.SetValue(1, 4, c[1] / Au[1]);
233 J.SetValue(2, 5, c[2] / Au[2]);
235 J.SetValue(3, 0, Au[0]);
236 J.SetValue(3, 1, -Au[1]);
237 J.SetValue(3, 2, -Au[2]);
238 J.SetValue(3, 3, uu[0]);
239 J.SetValue(3, 4, -uu[1]);
240 J.SetValue(3, 5, -uu[2]);
242 J.SetValue(4, 0, 2 * uu[0]);
243 J.SetValue(4, 1, -2 * uu[1]);
245 J.SetValue(4, 3, 2 * c[0] * c[0] / Au[0]);
246 J.SetValue(4, 4, -2 * c[1] * c[1] / Au[1]);
249 J.SetValue(5, 0, 2 * uu[0]);
251 J.SetValue(5, 2, -2 * uu[2]);
252 J.SetValue(5, 3, 2 * c[0] * c[0] / Au[0]);
254 J.SetValue(5, 5, -2 * c[2] * c[2] / Au[2]);
259 else if (type ==
"Interface")
264 for (
int i = 0; i < 2; ++i)
266 GetC(c[i],
beta[i], Au[i], A0[i], alpha[i]);
271 J.SetValue(0, 2, c[0] / Au[0]);
277 J.SetValue(1, 3, -c[1] / Au[1]);
279 J.SetValue(2, 0, Au[0]);
280 J.SetValue(2, 1, -Au[1]);
281 J.SetValue(2, 2, uu[0]);
282 J.SetValue(2, 3, -uu[1]);
284 J.SetValue(3, 0, 2 * uu[0]);
285 J.SetValue(3, 1, -2 * uu[1]);
286 J.SetValue(3, 2, 2 * c[0] * c[0] / Au[0]);
287 J.SetValue(3, 3, -2 * c[1] * c[1] / Au[1]);
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
void v_GetCharIntegral(NekDouble &I, const NekDouble &beta, const NekDouble &A, const NekDouble &A0, const NekDouble &alpha=0.5) override
void GetC0(NekDouble &c0, const NekDouble &beta, const NekDouble &A0)
void v_GetUFromChars(NekDouble &u, const NekDouble &W1, const NekDouble &W2) override
PowerPressureArea(Array< OneD, MultiRegions::ExpListSharedPtr > pVessel, const LibUtilities::SessionReaderSharedPtr pSession)
void v_GetW2(NekDouble &W2, const NekDouble &u, const NekDouble &beta, const NekDouble &A, const NekDouble &A0, const NekDouble &alpha=0.5) override
static std::string className
void v_GetPressure(NekDouble &P, const NekDouble &beta, const NekDouble &A, const NekDouble &A0, const NekDouble &dAUdx, const NekDouble &gamma=0, const NekDouble &alpha=0.5) override
void v_GetJacobianInverse(NekMatrix< NekDouble > &invJ, const Array< OneD, NekDouble > &Au, const Array< OneD, NekDouble > &uu, const Array< OneD, NekDouble > &beta, const Array< OneD, NekDouble > &A0, const Array< OneD, NekDouble > &alpha, const std::string &type) override
void GetB(NekDouble &b, const NekDouble &c0)
void v_GetAFromChars(NekDouble &A, const NekDouble &W1, const NekDouble &W2, const NekDouble &beta, const NekDouble &A0, const NekDouble &alpha=0.5) override
void v_GetC(NekDouble &c, const NekDouble &beta, const NekDouble &A, const NekDouble &A0, const NekDouble &alpha=0.5) override
void v_GetW1(NekDouble &W1, const NekDouble &u, const NekDouble &beta, const NekDouble &A, const NekDouble &A0, const NekDouble &alpha=0.5) override
static PulseWavePressureAreaSharedPtr create(Array< OneD, MultiRegions::ExpListSharedPtr > &pVessel, const LibUtilities::SessionReaderSharedPtr &pSession)
void GetCharIntegral(NekDouble &I, const NekDouble &beta, const NekDouble &A, const NekDouble &A0, const NekDouble &alpha=0.5)
void GetC(NekDouble &c, const NekDouble &beta, const NekDouble &A, const NekDouble &A0, const NekDouble &alpha=0.5)
LibUtilities::SessionReaderSharedPtr m_session
std::shared_ptr< SessionReader > SessionReaderSharedPtr
@ beta
Gauss Radau pinned at x=-1,.
@ P
Monomial polynomials .
PressureAreaFactory & GetPressureAreaFactory()
scalarT< T > sqrt(scalarT< T > in)