diff --git a/sumpy/expansion/__init__.py b/sumpy/expansion/__init__.py
index b5415b1e60250a732bcc8af141960fd57b33b9c0..08754a10240061ac9858034a561fb301c7db6925 100644
--- a/sumpy/expansion/__init__.py
+++ b/sumpy/expansion/__init__.py
@@ -495,6 +495,11 @@ class DefaultExpansionFactory(ExpansionFactoryBase):
                 and base_kernel.dim == 2):
             from sumpy.expansion.local import Y2DLocalExpansion
             return Y2DLocalExpansion
+        elif (isinstance(base_kernel.get_base_kernel(), HelmholtzKernel)
+                and base_kernel.dim == 3):
+            from sumpy.expansion.local import \
+                    HelmholtzConformingVolumeTaylorLocalExpansion
+            return HelmholtzConformingVolumeTaylorLocalExpansion
         elif isinstance(base_kernel.get_base_kernel(), LaplaceKernel):
             from sumpy.expansion.local import \
                     LaplaceConformingVolumeTaylorLocalExpansion
@@ -519,6 +524,11 @@ class DefaultExpansionFactory(ExpansionFactoryBase):
             from sumpy.expansion.multipole import (
                     LaplaceConformingVolumeTaylorMultipoleExpansion)
             return LaplaceConformingVolumeTaylorMultipoleExpansion
+        elif (isinstance(base_kernel.get_base_kernel(), HelmholtzKernel)
+                and base_kernel.dim == 3):
+            from sumpy.expansion.multipole import (
+                    HelmholtzConformingVolumeTaylorMultipoleExpansion)
+            return HelmholtzConformingVolumeTaylorMultipoleExpansion
         else:
             from sumpy.expansion.multipole import VolumeTaylorMultipoleExpansion
             return VolumeTaylorMultipoleExpansion