diff --git a/examples/fluid/static_sphere_shape.jl b/examples/fluid/static_sphere_shape.jl
deleted file mode 100644
index 4c414407c..000000000
--- a/examples/fluid/static_sphere_shape.jl
+++ /dev/null
@@ -1,92 +0,0 @@
-# In this example we try to approach the static shape of a water droplet on a horizontal plane.
-# The shape of a static droplet can be calculated from the Young-Laplace equation.
-using TrixiParticles
-using OrdinaryDiffEq
-
-# ==========================================================================================
-# ==== Resolution
-fluid_particle_spacing = 0.0025
-
-boundary_layers = 4
-spacing_ratio = 1
-
-# ==========================================================================================
-# ==== Experiment Setup
-gravity = 9.81
-tspan = (0.0, 0.3)
-
-# Boundary geometry and initial fluid particle positions
-initial_fluid_size = (0.0, 0.0)
-tank_size = (1.0, 0.1)
-
-fluid_density = 1000.0
-sound_speed = 120.0
-state_equation = StateEquationCole(; sound_speed, reference_density=fluid_density,
-                                   exponent=1)
-
-tank = RectangularTank(fluid_particle_spacing, initial_fluid_size, tank_size, fluid_density,
-                       n_layers=boundary_layers, spacing_ratio=spacing_ratio,
-                       faces=(true, true, true, false),
-                       acceleration=(0.0, -gravity), state_equation=state_equation)
-
-sphere_radius = 0.05
-
-sphere1_center = (0.5, sphere_radius)
-sphere1 = SphereShape(fluid_particle_spacing, sphere_radius, sphere1_center,
-                      fluid_density, sphere_type=VoxelSphere())
-
-# ==========================================================================================
-# ==== Fluid
-fluid_smoothing_length = 1.0 * fluid_particle_spacing
-fluid_smoothing_kernel = SchoenbergCubicSplineKernel{2}()
-
-fluid_density_calculator = ContinuityDensity()
-
-# for perfect wetting
-nu = 0.00025
-# for no wetting
-# nu = 0.001
-
-alpha = 8 * nu / (fluid_smoothing_length * sound_speed)
-viscosity = ArtificialViscosityMonaghan(alpha=alpha, beta=0.0)
-
-sphere_surface_tension = WeaklyCompressibleSPHSystem(sphere1, fluid_density_calculator,
-                                                     state_equation, fluid_smoothing_kernel,
-                                                     fluid_smoothing_length,
-                                                     viscosity=viscosity,
-                                                     density_diffusion=nothing,
-                                                     acceleration=(0.0, -gravity),
-                                                     surface_tension=SurfaceTensionAkinci(surface_tension_coefficient=0.01),
-                                                     correction=AkinciFreeSurfaceCorrection(fluid_density))
-
-# ==========================================================================================
-# ==== Boundary
-boundary_density_calculator = AdamiPressureExtrapolation()
-boundary_model = BoundaryModelDummyParticles(tank.boundary.density, tank.boundary.mass,
-                                             state_equation=state_equation,
-                                             boundary_density_calculator,
-                                             fluid_smoothing_kernel,
-                                             fluid_smoothing_length,
-                                             viscosity=ViscosityAdami(nu=2.0 * nu))
-
-# adhesion_coefficient = 1.0 and surface_tension_coefficient=0.01 for perfect wetting
-# adhesion_coefficient = 0.001 and surface_tension_coefficient=2.0 for no wetting
-boundary_system = BoundarySPHSystem(tank.boundary, boundary_model, adhesion_coefficient=1.0)
-
-# ==========================================================================================
-# ==== Simulation
-semi = Semidiscretization(boundary_system, sphere_surface_tension)
-ode = semidiscretize(semi, tspan)
-
-info_callback = InfoCallback(interval=50)
-saving_callback = SolutionSavingCallback(dt=0.01, output_directory="out", prefix="",
-                                         write_meta_data=true)
-
-callbacks = CallbackSet(info_callback, saving_callback)
-
-# Use a Runge-Kutta method with automatic (error based) time step size control.
-sol = solve(ode, RDPK3SpFSAL35(),
-            abstol=1e-6, # Default abstol is 1e-6
-            reltol=1e-4, # Default reltol is 1e-3
-            dt=1e-5,
-            save_everystep=false, callback=callbacks);