diff --git a/examples/helmholtzFlow/freefield/README.md b/examples/helmholtzFlow/freefield/README.md
index 2121bb85479f8580f1e86f44532e74ee480320bf..df6beb8c3ffaa05212110a98e592fbbf3245df2f 100644
--- a/examples/helmholtzFlow/freefield/README.md
+++ b/examples/helmholtzFlow/freefield/README.md
@@ -31,7 +31,7 @@ with `[PARAM]`:
 By default a circular domain with a Padé condition is used, for a centered source at the frequency $k=6\pi$, with a mean flow at $M=0.8$ oriented at $\theta=\pi/4$.
 
 ## Reference
-> To be updated
+> P. Marchner, H. Bériot, S. Le Bras, X. Antoine and C. Geuzaine. A domain decomposition solver for large scale time-harmonic flow acoustics problems. (2023), Preprint [HAL-04254633](https://hal.science/hal-04254633).
 
 ## Results reproducibility
 We assess the performance of different ABCs in the case of a radiating point source set in a circular domain.