In this work we provide a novel energy-consistent formulation for the classical 1D formulation of blood flow in an arterial segment. The resulting reformulation is shown to be suitable for the coupling with a lumped (0D) model of the heart that incorporates a reduced formulation of the actin-myosin interaction. The coupling being consistent with energy balances, we provide a complete heart-circulation model compatible with thermodynamics hence stable numerically and informative physiologically. These latter two properties are verified by numerical experiments.

The goal of this work is to derive a reliable – stable and accurate – inverse problem strategy for reconstructing cardiac output – blood flow entering the ascending aorta – from pressure measurements at a distal site of the arterial tree, assumed here to be the descending aorta. We assume that a reduced one-dimensional model of the aorta can be linearized around its steady state, resulting in a wave system with absorbing boundary condition at the outlet. Using this model, we attempt to reconstruct the inlet flow from a pressure measurement at the distal outlet. First, we investigate the observability of the problem and prove that the inversion of the input-output operator for the flow and pressure in the space of time-periodic solutions is ill-posed of degree one. We then develop a variational approach where we minimize the discrepancy between measurements and a simulated state and penalize the error with respect to a periodic state. It is shown that the penalty strategy is convergent and provides an efficient solution for the minimization. Numerical results illustrate the robustness of our approach to noise and the potential of our method to reconstruct inlet flow from real pressure recordings during anesthesia.