Abstract:In this study, we develop vector flow imaging techniques for multi-layered models with a high wavespeed contrast using photoacoustic and ultrasonic imaging. We use refraction-corrected delay-and-sum image reconstruction (RC-DAS), which enforces Snell's law to accurately calculate time delays within each layer. We compare RC-DAS against conventional delay-and-sum for vector flow imaging in benchtop phantoms made of transparent polymethyl methacrylate (PMMA) in a water bath. We study the flow beneath a PMMA layer using two phantoms, where the PMMA layer has different shapes and thicknesses. We image a slow-moving suspension of carbon microspheres (~4 mm/s) using interleaved photoacoustic and multi-angle plane wave ultrasound acquisitions measured with a 7.6 MHz linear ultrasound array. Photoacoustic waves are generated by a 1064 nm wavelength nanosecond-pulsed laser at 50 Hz, and multi-angle plane wave ultrasound data are acquired at 100 Hz for eleven steering angles between $\pm$10$^\circ$. RC-DAS improves the flow speed accuracy, reducing the mean absolute error by 0.41-0.63 mm/s compared to the expected flow profile. The error in direction estimates improves when we use RC-DAS, with the interdecile range reducing by up to 17$^\circ$. This work emphasises the importance of refraction correction for accurate flow measurements in layered media with photoacoustics and ultrasonic imaging. While both imaging modalities can quantify flow in these multi-layered models, the modality best suited for a specific application will depend on the imaging target and flow dynamics. These techniques show promise for biomedical applications such as intraosseous and transcranial blood flow quantification, and in nondestructive testing to monitor fluid motion.