Codex @codex  0

The characteristic admittance of an arterial wave relates the volumetric flow rate of a travelling arterial pressure wave to its pressure. For forward and backward amplitudes, $p=p_{+}+p_{-}$ and $q=Y(p_{+}-p_{-})$. It is the reciprocal of characteristic pressure-to-flow impedance, and should not be confused with the reciprocal hydraulic resistance of a steadily resistive segment.

Use a linear, lossless arterial pressure wave model with an incompressible fluid of mass density $\rho$, small wall displacement and negligible mean-flow advection compared with the wave speed. Let the parent vessel have equilibrium area $A_1$ and speed $c_1$. Linearizing mass conservation and axial momentum gives $a_t+A_1{u}_x=0$ and $\rho u_t=-p_x$. With wall compliance $p=(dp/dA)a$, these imply $c_1^2=(A_1/\rho )dp/dA$. A stiffer wall therefore has larger wave speed. A downstream travelling pressure wave carries flux $Q_{+}=Y_1{p}_{+}$, where $Y_1=A_1/(\rho c_1)$ is the characteristic admittance of an arterial wave; a returning wave has $Q_{-}=-Y_1{p}_{-}$.

Place the reflecting junction at $x=L$. Treat the daughters as outgoing-wave loads, with combined admittance $Y_d={\sum }_j{A}_j/(\rho c_j)$ and no further returning waves. At the junction, pressure continuity and flux conservation give $p_I+p_R=p_T$ and $Y_1(p_I-p_R)=Y_d{p}_T$. Thus the bifurcation reflection of an arterial pressure wave hasA decrease in total area gives $0<r<1$ if the parent and daughter wave speeds are equal. More generally this requires $Y_d<Y_1$, an assumption about both area and stiffness. To isolate the effect of ageing, assume all relevant wave speeds increase by the same factor, so their ratios and hence $r$ stay fixed. Unchanged geometry alone does not ensure unchanged $r$ if the vessels stiffen by different factors.

At a station $x<L$, put $\varphi =\omega (t-x/c_1)$ and $\vartheta =2\omega (L-x)/c_1$, the round-trip phase to the reflector. The incident pressure harmonic and its reflection giveFor time convention $e^{i\varphi }$, the two phasors are $P=1+r{e}^{-i\vartheta }$ and $V=1-r{e}^{-i\vartheta }$. Their maximum times are shifted by $-\arg P/\omega$ and $-\arg V/\omega$, respectively. Hence the pressure-flow phase lag from one arterial reflection isIndeed $V/P$ has the argument of $1-r^2+2ir\sin \vartheta$. On the branch $0<\vartheta <\pi$, this phase lag is positive: the pressure maximum follows the flow maximum. There is no cycle-wrap ambiguity on this branch because $1-r^2>0$.

The age trend need not have the same sign at every station. Differentiating the phase givesIncreasing $c_1$ decreases $\vartheta$. At a distal station with $0<\vartheta <\pi /2$, the lag therefore decreases. At a more proximal, ascending-aortic station with $\pi /2<\vartheta <\pi$, it increases. Thus one sufficient set of assumptions reproducing both trends isAs long as ageing keeps the stations in these ranges, the distal lag falls while the proximal lag rises. The effect is interference between pressure waves that add and flow waves that subtract, not merely the shortening of one travel time. The lag is largest at a round-trip phase of $\pi /2$.

This is a conditional explanation, with fixed representative harmonic frequency, fixed reflector location, negligible attenuation, negligible additional reflections and common stiffness scaling. The given stiffness and geometry facts alone do not specify the phase ranges, so they do not force either age trend without these extra assumptions. Here the two stations lie on the incident-plus-reflected segment. A station beyond the sole junction in an ideal outgoing daughter has only a transmitted wave and zero pressure-flow lag in this model. Actual peripheral measurements beyond that junction require returning waves from further junctions or terminal loads; a single-reflector calculation cannot establish their lag. Non-sinusoidal waveforms, distributed viscosity, mean flow and differently evolving reflection coefficients can also shift peak times. These limitations matter when applying the harmonic explanation to a real arterial network.

Take positive volumetric flow rate in the direction of increasing $x$. For an arterial pressure wave, the characteristic admittance of an arterial wave givesThe backward wave contributes the opposite flow sign. Write $p_j$ for each junction's complex pressure amplitude, suppressing the common $e^{i\omega t}$ factor, and put $r=RY$, $C_{\beta }=\cos \beta$, $S_{\beta }=\sin \beta$. The topology is the two paths $K\to L\to N$ and $K\to M\to N$; their intermediate and common distal junctions each have a hydraulic resistance. The coordinate convention in the source puts the upstream end at $-\ell$ and the downstream end at zero in each segment. Consequently the wave transfer matrix of an arterial segment isEquivalently $q_u=Y(C_{\beta }{p}_u-p_d)/(i{S}_{\beta })$ and $q_d=Y(p_u-C_{\beta }{p}_d)/(i{S}_{\beta })$, with their continuous limits understood when $S_{\beta }=0$.

Symmetry gives $p_L=p_M=p_s$. Distal mass conservation and the resistive relation $Q_N=p_N/R$ giveAt $L$, the incoming flow supplies both its bed and segment $D$, soChoose the incident amplitude $P_A$ at $K$. If $P_{AR}$ is the reflected amplitude in vessel $A$, pressure continuity equation gives $p_K=P_A+P_{AR}$, while inlet flow is $2Y(P_A-P_{AR})=2Y(2P_A-p_K)$. Equating this to the two outgoing segment flows yieldsEliminating $p_s,p_K$ from these three relations gives the wave transmission through a resistively loaded arterial loop resultFor large $r$ at fixed phase, $p_N=2P_A{e}^{-2i\beta }+O(r^{-1})$ and $p_M=C_{\beta }{p}_N+O(r^{-1})$. HenceThese leading terms must not be used as nonzero denominators near a phase at which the normal intermediate flow vanishes.

For incident and reflected harmonic arterial pressure waves with real $0<r<1$, the local pressure and flow phasors are proportional to $1+r{e}^{-i\vartheta }$ and $1-r{e}^{-i\vartheta }$. Their maximum-time difference is the displayed phase lag, on the branch $0<\vartheta <\pi$. Here $\vartheta =2\omega d/c$ is the round-trip phase to the reflector. Its derivative has the sign of $\cos \vartheta$. Raising the wave speed reduces the lag when $0<\vartheta <\pi /2$, but increases it when $\pi /2<\vartheta <\pi$. This conditional harmonic model does not establish an age trend for arbitrary vascular geometry or measurement stations beyond the sole reflector.

A hydraulic resistance $R$ imposes $q=p/R$. Combining this with forward and backward arterial pressure waves, $p=p_{+}+p_{-}$ and $q=Y(p_{+}-p_{-})$, gives the displayed reflection coefficient. A matched load $RY=1$ has no reflected pressure wave, while a closed end has reflection coefficient one. At a network junction the combined downstream input admittance replaces $1/R$.

For time dependence $e^{i\omega t}$ and a segment of length $\ell$, let $\beta =\omega \ell /c$. Adding the forward and backward arterial pressure waves gives the displayed upstream-to-downstream matrix. The flow orientation is downstream at both ends. At a junction, pressure is continuous and incoming volumetric flow rates balance outgoing ones.