diff --git a/sumpy/expansion/local.py b/sumpy/expansion/local.py
index bff7acee7e1003d32a023dfc04aa5bcf7d3f8781..ae408283039a3d89eafdf5f556b27dfe816754f6 100644
--- a/sumpy/expansion/local.py
+++ b/sumpy/expansion/local.py
@@ -175,7 +175,7 @@ class H2DLocalExpansion(LocalExpansionBase):
             return translated_coeffs
 
         raise RuntimeError("do not know how to translate %s to "
-                           "local 2D Hankel expansion"
+                           "local 2D Helmholtz Bessel expansion"
                            % type(src_expansion).__name__)
 
 # }}}
diff --git a/sumpy/expansion/multipole.py b/sumpy/expansion/multipole.py
index 78624ff64cedd48b2c82e87ab2cd10d5e2112b70..7deb7353e0d6e7dad9960c40fc584b7e26dd61bd 100644
--- a/sumpy/expansion/multipole.py
+++ b/sumpy/expansion/multipole.py
@@ -173,7 +173,7 @@ class H2DMultipoleExpansion(MultipoleExpansionBase):
     def translate_from(self, src_expansion, src_coeff_exprs, dvec):
         if not isinstance(src_expansion, type(self)):
             raise RuntimeError("do not know how to translate %s to "
-                               "multipole 2D Hankel expansion"
+                               "multipole 2D Helmholtz Bessel expansion"
                                % type(src_expansion).__name__)
         from sumpy.symbolic import sympy_real_norm_2
         dvec_len = sympy_real_norm_2(dvec)