arxiv: 2604.08307 · v1 · submitted 2026-04-09 · 💻 cs.ET

Recognition: unknown

Analytical Modeling of Dispersive Closed-loop MC Channels with Pulsatile Flow

Theofilos Symeonidis , Fardad Vakilipoor , Robert Schober , Nunzio Tuccitto , Maximilian Sch\"afer

Authors on Pith no claims yet

Pith reviewed 2026-05-10 17:35 UTC · model grok-4.3

classification 💻 cs.ET

keywords molecular communicationpulsatile flowchannel impulse responsewrapped normal distributionclosed-loop channelsdispersive transportcardiovascular system

0 comments

The pith

An analytical model expresses the impulse response of pulsatile-flow molecular communication channels as a wrapped normal distribution with time-variant parameters.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper creates a mathematical description for molecular communication inside the cardiovascular system, where blood flow pulses with each heartbeat. It provides a closed-form expression for the channel impulse response in a looping channel under dispersive conditions and pulsatile flow. The expression uses a wrapped normal distribution whose mean and variance change with time to capture the periodic pulsing. A reader would care because it offers a fast way to predict molecule concentrations without running complex simulations, helping design systems for health monitoring or targeted delivery inside the body.

Core claim

We derive an analytical expression for the channel impulse response (CIR), which follows a wrapped Normal distribution with time-variant mean and variance. The obtained model reveals the cyclostationary nature of the channel and quantifies the influence of pulsation on the temporal concentration profile compared to steady-flow systems. The model is validated by three-dimensional particle-based simulations showing excellent agreement.

What carries the argument

The time-variant wrapped Normal distribution used to model the channel impulse response in the one-dimensional approximation of the closed-loop channel with pulsatile flow.

If this is right

The model shows the channel is cyclostationary, with behavior repeating over cardiac cycles.
Pulsation alters the concentration profile in measurable ways relative to constant flow.
The analytical form supports efficient computation of molecule propagation.
Three-dimensional simulations confirm the one-dimensional model's accuracy for this setting.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Designers of in-body molecular systems could time molecule releases to exploit the pulsation cycles for better signal strength.
The model opens the door to analyzing multi-hop communication in branching blood vessels under realistic flow conditions.
Extensions to non-ideal vessel shapes or variable heart rates could be tested using the same wrapped-normal framework.

Load-bearing premise

The one-dimensional approximation combined with the specific incorporation of pulsatile flow into a time-variant wrapped normal distribution accurately captures the particle transport in the closed-loop channel.

What would settle it

A direct comparison of the analytical CIR predictions against experimental measurements of molecule concentrations in a physical closed-loop setup with pulsatile flow that shows large deviations would falsify the model.

Figures

Figures reproduced from arXiv: 2604.08307 by Fardad Vakilipoor, Maximilian Sch\"afer, Nunzio Tuccitto, Robert Schober, Theofilos Symeonidis.

**Figure 1.** Figure 1: Closed-loop molecular communication system and model reduction. Signaling molecules are released by the transmitter and propagate [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗

**Figure 2.** Figure 2: Velocity waveforms used for validating the proposed ana [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

**Figure 4.** Figure 4: Normalized received signal for different mean velocities [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗

**Figure 5.** Figure 5: Physiologically motivated pulsatile flow: Normalized re [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗

read the original abstract

Molecular communication (MC) is a communication paradigm in which information is conveyed through the controlled release, propagation, and reception of molecules. Many envisioned healthcare applications of MC are expected to operate inside the human body. In this environment, the cardiovascular system ( CVS) acts as the physical channel, which forms a closed-loop network where particle transport is mainly governed by the combined effects of diffusion and flow. Despite the fact that physiological flows in many parts of the human body are inherently pulsatile due to the cardiac cycle, most existing models for dispersive closed-loop MC channels assume a constant flow velocity. In this paper, we present a time-variant one-dimensional (1D ) channel model for dispersive closed-loop MC systems with pulsatile flow. We derive an analytical expression for the channel impulse response (CIR ), which follows a wrapped Normal distribution with time-variant mean and variance. The obtained model reveals the cyclostationary nature of the channel and quantifies the influence of pulsation on the temporal concentration profile compared to steady-flow systems. Finally, the model is validated by three-dimensional ( 3D ) particle-based simulations (PBS s), showing excellent agreement and enabling an efficient analytical characterization of the channel.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a usable analytical wrapped-normal CIR for pulsatile flow in closed-loop MC by making the mean time-dependent while keeping variance linear in t, but the 1D uniform-advection step likely misses radial dispersion from the Womersley profile.

read the letter

The new piece is the explicit time-variant wrapped normal for the impulse response under periodic flow. They integrate the pulsatile velocity for the mean and keep the usual diffusion term for the variance, then wrap it for the closed loop. This captures the cyclostationary character of the channel and lets you see how the pulsation stretches or compresses the concentration peaks compared with constant flow. The validation is the strongest part. Three-dimensional particle-based simulations line up closely with the analytical curves, which suggests the expressions are accurate enough for the tested conditions and parameters. The soft spot sits in the modeling assumptions. The derivation starts from the one-dimensional advection-diffusion equation on a circle. That only holds exactly when the flow velocity is uniform across the cross-section and the diffusion coefficient does not change with the flow rate. Physiological pulsatile flow follows the Womersley profile, so radial shear produces extra unsteady dispersion that a constant effective D does not automatically include. If the paper treats D as fixed, the model could underpredict spreading in real vessels even while the mean transport stays right. The abstract claims excellent agreement, but without seeing how they chose the effective parameters it is hard to judge how general the match is. This work is for molecular communication researchers focused on cardiovascular or in-body channels who need fast analytical predictions rather than full numerical runs. It is coherent and grounded enough to go to peer review, where the main questions will be about the dispersion term and the range of Womersley numbers it covers. I would bring this to a reading group to discuss the derivation steps and the simulation setup. I would not cite it myself in the next year since it is outside my area, but the contribution looks worth the time for the right audience.

Referee Report

2 major / 1 minor

Summary. The paper claims to derive an analytical expression for the channel impulse response (CIR) in dispersive closed-loop molecular communication channels with pulsatile flow. The CIR is modeled as following a wrapped Normal distribution with a time-variant mean (the time-integral of the pulsatile velocity) and time-variant variance (growing as 2Dt). The model is presented as capturing the cyclostationary nature of the channel and is validated against 3D particle-based simulations with claimed excellent agreement.

Significance. If the result holds, this analytical model offers a significant advance by providing a closed-form, computationally efficient way to characterize time-varying MC channels in physiological pulsatile flow environments, such as the cardiovascular system. The explicit derivation of the wrapped-normal form and the demonstration of cyclostationarity, supported by 3D simulation validation, would enable better design and analysis of in-body MC systems compared to steady-flow assumptions.

major comments (2)

[Model derivation] The derivation of the time-variant mean and variance expressions assumes a one-dimensional advection-diffusion process with spatially uniform velocity and constant diffusion coefficient. This leads to the wrapped Normal distribution. However, the skeptic note highlights that physiological pulsatile flow involves radially varying velocity profiles, leading to unsteady Taylor-Aris dispersion that may not be captured by a constant D independent of flow rate. This assumption underpins the central claim of an exact analytical CIR for the closed-loop channel.
[Simulation validation] The abstract states 'excellent agreement' with 3D particle-based simulations, but the full paper should provide quantitative error metrics, the range of parameters tested (e.g., Womersley number), and how the effective D was determined to allow assessment of the model's robustness.

minor comments (1)

[Abstract] Minor typographical issues include inconsistent spacing around acronyms, such as ' ( CVS)' and ' ( 3D )', which should be standardized.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments, which help improve the clarity and robustness of our manuscript. We address each major comment below and indicate the planned revisions.

read point-by-point responses

Referee: [Model derivation] The derivation of the time-variant mean and variance expressions assumes a one-dimensional advection-diffusion process with spatially uniform velocity and constant diffusion coefficient. This leads to the wrapped Normal distribution. However, the skeptic note highlights that physiological pulsatile flow involves radially varying velocity profiles, leading to unsteady Taylor-Aris dispersion that may not be captured by a constant D independent of flow rate. This assumption underpins the central claim of an exact analytical CIR for the closed-loop channel.

Authors: Our derivation is exact under the stated 1D assumptions of uniform cross-sectional velocity and constant effective D, which enable the closed-form wrapped Normal solution with time-variant mean (time-integral of pulsatile velocity) and variance (linear growth 2Dt). This is a deliberate modeling choice common in MC literature to obtain analytical tractability for closed-loop cyclostationary channels. We acknowledge that real physiological pulsatile flows exhibit radial velocity profiles and unsteady Taylor-Aris effects that can render effective dispersion time- or flow-dependent. The manuscript already includes a brief discussion of these modeling assumptions; we will expand this section in revision to explicitly note the limitations and position the model as a 1D approximation suitable for initial analysis of longitudinal transport in pulsatile MC systems. revision: partial
Referee: [Simulation validation] The abstract states 'excellent agreement' with 3D particle-based simulations, but the full paper should provide quantitative error metrics, the range of parameters tested (e.g., Womersley number), and how the effective D was determined to allow assessment of the model's robustness.

Authors: We agree that quantitative metrics and parameter details would strengthen the validation section. The current manuscript relies on visual overlays of analytical and simulated concentration profiles to claim excellent agreement. In the revised version we will add quantitative error measures (e.g., mean absolute percentage error or integrated squared difference between analytical and simulated CIRs), explicitly list the tested parameter ranges including Womersley numbers representative of cardiovascular flows, and describe the procedure used to select or fit the effective diffusion coefficient D (either from literature values or by matching steady-flow dispersion in the 3D simulations). These additions will be placed in a new subsection of the numerical results. revision: yes

Circularity Check

0 steps flagged

Derivation from 1D advection-diffusion PDE is mathematically independent

full rationale

The paper derives the CIR by solving the one-dimensional advection-diffusion equation with time-dependent pulsatile velocity on a closed circular domain, yielding a wrapped normal whose mean is the time-integral of velocity and variance grows linearly as 2Dt. This is the direct analytic solution of the PDE under the uniform-velocity and constant-D assumptions; no parameter is fitted to the target CIR data inside the derivation, and the form is not smuggled via self-citation or ansatz. The 3D PBS validation is external and post-hoc rather than an input that forces the result. No load-bearing step reduces to a tautology or to the authors' prior fitted quantities.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The model rests on standard advection-diffusion transport adapted to periodic flow and a closed loop; no new particles or forces are introduced.

axioms (2)

domain assumption Particle transport in the channel is adequately captured by a one-dimensional advection-diffusion equation with time-periodic velocity.
Invoked to justify the 1D model and the resulting wrapped-normal form.
domain assumption The closed-loop geometry combined with periodic pulsation produces a cyclostationary concentration profile that can be represented by a wrapped normal distribution.
Central modeling choice stated in the abstract.

pith-pipeline@v0.9.0 · 5530 in / 1278 out tokens · 64344 ms · 2026-05-10T17:35:49.255923+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

27 extracted references · 18 canonical work pages

[1]

I. F. Akyildiz, M. Pierobon, and S. Balasubramaniam. 2019. Moving forward with molecular communication: from theory to human health applications.Proc. IEEE 107, 5 (2019), 858–865. doi:10.1109/JPROC.2019.2913890

work page doi:10.1109/jproc.2019.2913890 2019
[2]

R. Aris. 1956. On the dispersion of a solute in a fluid flowing through a tube. Proc. R. Soc. Lond. A.235, 1200 (1956), 67–77

1956
[3]

Gabriel U Talasso, Allan M de Souza, Luis FG Gonzalez, Eduardo Cerqueira, Antonio AF Loureiro, and Leandro A Villas

L. Brand et al . 2024. Closed Loop Molecular Communication Testbed: Setup, Interference Analysis, and Experimental Results. InProc. IEEE Int. Conf. on Commun.doi:10.1109/ICC51166.2024.10622231

work page doi:10.1109/icc51166.2024.10622231 2024
[4]

Brand et al

L. Brand et al. 2025. Closed-Loop Molecular Communication with Local and Global Degradation: Modeling and ISI Analysis. InProc. 12th ACM Int. Conf. Nanoscale Comput. Commun.(Chengdu, China)

2025
[5]

Chahibi, M

Y. Chahibi, M. Pierobon, S. Song, and I. F. Akyildiz. 2013. A Molecular Communi- cation System Model for Particulate Drug Delivery Systems.IEEE Trans. Biomed. Eng.He60, 12 (2013), 3468–3483. doi:10.1109/TBME.2013.2271503

work page doi:10.1109/tbme.2013.2271503 2013
[6]

Dincau, E

B. Dincau, E. Dressaire, and A. Sauret. 2020. Pulsatile Flow in Microfluidic Systems.Small16, 9 (2020), 1904032

2020
[7]

Felicetti, M

L. Felicetti, M. Femminella, G. Reali, and P. Liò. 2016. Applications of molecular communications to medicine: A survey.Nano Commun. Netw.7 (2016), 27–45. doi:10.1016/j.nancom.2015.08.004

work page doi:10.1016/j.nancom.2015.08.004 2016
[8]

Furubayashi, Y

T. Furubayashi, Y. Sakatani, T. Nakano, A. W. Eckford, and N. Ichihashi. 2017. Design and wet-laboratory implementation of reliable end-to-end molecular communication.Wireless Networks24 (2017), 1809 – 1819

2017
[9]

Gay and L

M. Gay and L. T. Zhang. 2009. Numerical studies of blood flow in healthy, stenosed, and stented carotid arteries.Int. Num. Methods Fluids61, 4 (2009), 453–472. doi:10.1002/fld.1966

work page doi:10.1002/fld.1966 2009
[10]

Grebenstein, J

L. Grebenstein, J. Kirchner, W. Wicke, A. Ahmadzadeh, V. Jamali, G. Fischer, R. Weigel, A. Burkovski, and R. Schober. 2019. A Molecular Communication Testbed Based on Proton Pumping Bacteria: Methods and Data.IEEE Trans. Mol., Biol., Multi-Scale Commun.5, 1 (2019), 56–62. doi:10.1109/TMBMC.2019.2957783

work page doi:10.1109/tmbmc.2019.2957783 2019
[11]

D W Holdsworth, C J D Norley, R Frayne, D A Steinman, and B K Rutt. [n. d.]. Characterization of common carotid artery blood-flow waveforms in normal human subjects.Physiological Measurement20, 3 ([n. d.]), 219. doi:10.1088/0967- 3334/20/3/301

work page doi:10.1088/0967-
[12]

Jakumeit, L

T. Jakumeit, L. Brand, J. Kirchner, R. Schober, and S. Lotter. 2025. Molecular Signal Reception in Complex Vessel Networks: The Role of the Network Topology. arXiv:2410.15943 [cs.ET]

work page arXiv 2025
[13]

Jakumeit, B

T. Jakumeit, B. Heinlein, N. Tuccitto, R. Schober, S. Lotter, and M. Schäfer. 2026. Mixture of Inverse Gaussians for Hemodynamic Transport (MIGHT) in Multiple- Input Multiple-Output Vascular Networks. arXiv:2510.11743 [q-bio.QM] Symeonidis et al

work page arXiv 2026
[14]

Jamali et al

V. Jamali et al. 2019. Channel Modeling for Diffusive Molecular Communica- tion—A Tutorial Review.Proc. IEEE107, 7 (2019), 1256–1301

2019
[15]

Abouei, A

G Kianfar, M Azadi, J. Abouei, A. Mohammadi, and K. N. Plataniotis. 2024. Wireless Body Area Nanonetworks via Vascular Molecular Communication. IEEE Trans. NanoBiosc.23, 2 (2024), 355–367. doi:10.1109/TNB.2024.3365737

work page doi:10.1109/tnb.2024.3365737 2024
[16]

Y. Lee. 2026. Corrected-Inverse-Gaussian First-Hitting-Time Modeling for Molec- ular Communication Under time-variant Drift. arXiv:2602.15335 [cs.IT]

work page arXiv 2026
[17]

Tunç Mcguiness et al

D. Tunç Mcguiness et al. 2019. Molecular-Based Nano-Communication Network: A Ring Topology Nano-Bots for In-Vivo Drug Delivery Systems.IEEE Access7 (2019), 12901–12913

2019
[18]

Mosayebi, A

R. Mosayebi, A. Ahmadzadeh, W. Wicke, V. Jamali, R. Schober, and M. Nasiri- Kenari. 2019. Early Cancer Detection in Blood Vessels Using Mobile Nanosensors. IEEE Trans. on NanoBiosc.18, 2 (2019), 103–116. doi:10.1109/TNB.2018.2885463

work page doi:10.1109/tnb.2018.2885463 2019
[19]

Schäfer, W

M. Schäfer, W. Wicke, L. Brand, R. Rabenstein, and R. Schober. 2021. Transfer Function Models for Cylindrical MC Channels With Diffusion and Laminar Flow. IEEE Trans. Mol., Biol., Multi-Scale Commun.7, 4 (2021), 271–287. doi:10.1109/ TMBMC.2021.3061030

work page arXiv 2021
[20]

T. W. Secomb. 2016. Hemodynamics.Comprehensive Physiology6, 2 (2016), 975–1003. doi:10.1002/j.2040-4603.2016.tb00700.x

work page doi:10.1002/j.2040-4603.2016.tb00700.x 2016
[21]

Vakilipoor et al

F. Vakilipoor et al. 2025. The CAM Model: An in Vivo Testbed for Molecular Communication Systems.IEEE Trans. Mol. Biol. Multi-Scale Commun.11, 4 (2025), 618–638. doi:10.1109/TMBMC.2025.3601432

work page doi:10.1109/tmbmc.2025.3601432 2025
[22]

Wang et al

Y. Wang et al. 2025. Modeling of Longitudinal Solute transport in Pulsatile Blood Flow Through A Slender Artery Using Invariant Manifold Method.Phys. Fluids 37, 9 (2025)

2025
[23]

Wicke, T

W. Wicke, T. Schwering, A. Ahmadzadeh, V. Jamali, A. Noel, and R. Schober
[24]

Modeling Duct Flow for Molecular Communication. InProc. IEEE Global Commun. Conf.206–212. doi:10.1109/GLOCOM.2018.8647632

work page doi:10.1109/glocom.2018.8647632 2018
[25]

In: ICC 2025 - IEEE International Conference on Communications

A. Wietfeld et al. 2025. Advanced Plaque Modeling for Atherosclerosis Detection Using Molecular Communication. InProc. IEEE Int. Conf. Commun. (ICC 2025). 6056–6062. doi:10.1109/ICC52391.2025.11161338

work page doi:10.1109/icc52391.2025.11161338 2025
[26]

Wille, C

L. Wille, C. Pfannenmüller, and J. Kirchner. 2025. From Steady to Pulsatile Flow in Molecular Communication: Propagation of Nanoparticles in Mid-Sized Arteries.IEEE Trans. on Mol., Biol., Multi-Scale Commun.11, 4 (2025), 531–536. doi:10.1109/TMBMC.2025.3608558

work page doi:10.1109/tmbmc.2025.3608558 2025
[27]

J. R. Womersley. 1955. Method for the Calculation of Velocity, Rate of Flow and Viscous Drag in Arteries when the Pressure Gradient is Known.J. Physiol. (Lond.)127, 3 (1955), 553–563

1955