diff --git a/test/test_mpi.jl b/test/test_mpi.jl
index 001d9bff86e..f36b669d3c3 100644
--- a/test/test_mpi.jl
+++ b/test/test_mpi.jl
@@ -23,6 +23,7 @@ CI_ON_WINDOWS = (get(ENV, "GITHUB_ACTIONS", false) == "true") && Sys.iswindows()
     # P4estMesh and T8codeMesh tests
     include("test_mpi_p4est_2d.jl")
     include("test_mpi_t8code_2d.jl")
+    include("test_mpi_t8code_fv.jl")
     if !CI_ON_WINDOWS # see comment on `CI_ON_WINDOWS` above
         include("test_mpi_p4est_3d.jl")
         include("test_mpi_t8code_3d.jl")
diff --git a/test/test_mpi_t8code_fv.jl b/test/test_mpi_t8code_fv.jl
new file mode 100644
index 00000000000..65d55f6179a
--- /dev/null
+++ b/test/test_mpi_t8code_fv.jl
@@ -0,0 +1,114 @@
+module TestExamplesMPIT8codeMesh2D
+
+using Test
+using Trixi
+
+include("test_trixi.jl")
+
+const EXAMPLES_DIR = pkgdir(Trixi, "examples", "t8code_2d_fv")
+
+@testset "T8codeMesh MPI FV" begin
+#! format: noindent
+
+# Run basic tests
+@testset "Examples 2D" begin
+    # Linear scalar advection
+    @trixi_testset "elixir_advection_basic.jl" begin
+        @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_basic.jl"),
+                            l2=[0.1419061449384701],
+                            linf=[0.2086802087402776])
+
+        # @testset "error-based step size control" begin
+        #     Trixi.mpi_isroot() && println("-"^100)
+        #     Trixi.mpi_isroot() &&
+        #         println("elixir_advection_basic.jl with error-based step size control")
+
+        #     sol = solve(ode, RDPK3SpFSAL35(); abstol = 1.0e-4, reltol = 1.0e-4,
+        #                 ode_default_options()..., callback = callbacks)
+        #     summary_callback()
+        #     errors = analysis_callback(sol)
+        #     if Trixi.mpi_isroot()
+        #         @test errors.l2≈[3.3022040342579066e-5] rtol=1.0e-4
+        #         @test errors.linf≈[0.00011787417954578494] rtol=1.0e-4
+        #     end
+        # end
+
+        # Ensure that we do not have excessive memory allocations
+        # (e.g., from type instabilities)
+        let
+            t = sol.t[end]
+            u_ode = sol.u[end]
+            du_ode = similar(u_ode)
+            @test (@allocated Trixi.rhs!(du_ode, u_ode, semi, t)) < 1000
+        end
+    end
+
+    @trixi_testset "elixir_advection_gauss.jl" begin
+        @test_trixi_include(joinpath(EXAMPLES_DIR,
+                                     "elixir_advection_gauss.jl"),
+                            l2=[0.5598148317954682],
+                            linf=[0.6301130236005371])
+        # Ensure that we do not have excessive memory allocations
+        # (e.g., from type instabilities)
+        let
+            t = sol.t[end]
+            u_ode = sol.u[end]
+            du_ode = similar(u_ode)
+            @test (@allocated Trixi.rhs!(du_ode, u_ode, semi, t)) < 1000
+        end
+    end
+
+    @trixi_testset "elixir_euler_blast_wave.jl" begin
+        @test_trixi_include(joinpath(EXAMPLES_DIR,
+                                     "elixir_euler_blast_wave.jl"),
+                            l2=[
+                                0.49698968976388164,
+                                0.16934401479236502,
+                                0.16934401479236502,
+                                0.6743947137176176,
+                            ],
+                            linf=[
+                                1.1342505243873413,
+                                0.43853745700004154,
+                                0.4385374570000415,
+                                3.009703218658938,
+                            ])
+        # Ensure that we do not have excessive memory allocations
+        # (e.g., from type instabilities)
+        let
+            t = sol.t[end]
+            u_ode = sol.u[end]
+            du_ode = similar(u_ode)
+            @test (@allocated Trixi.rhs!(du_ode, u_ode, semi, t)) < 1000
+        end
+    end
+
+    @trixi_testset "elixir_euler_kelvin_helmholtz_instability.jl" begin
+        @test_trixi_include(joinpath(EXAMPLES_DIR,
+                                     "elixir_euler_kelvin_helmholtz_instability.jl"),
+                            l2=[
+                                0.2542045564471016,
+                                0.22153069577606582,
+                                0.11870840559952726,
+                                0.03626114330454897,
+                            ],
+                            linf=[
+                                0.5467901048636064,
+                                0.4156157765819209,
+                                0.26176688262532194,
+                                0.0920608815870434,
+                            ],
+                            tspan=(0.0, 1.0))
+        # Ensure that we do not have excessive memory allocations
+        # (e.g., from type instabilities)
+        let
+            t = sol.t[end]
+            u_ode = sol.u[end]
+            du_ode = similar(u_ode)
+            @test (@allocated Trixi.rhs!(du_ode, u_ode, semi, t)) < 1000
+        end
+    end
+end
+end # T8codeMesh MPI
+
+end # module
diff --git a/test/test_t8code_fv.jl b/test/test_t8code_fv.jl
index 4497b74ff14..3e9b754bd60 100644
--- a/test/test_t8code_fv.jl
+++ b/test/test_t8code_fv.jl
@@ -112,7 +112,7 @@ end
     end
 end
 
-# TODO: or better, to test
+# TODO: Add some free stream test like the following
 # @trixi_testset "elixir_euler_free_stream.jl" begin
 #     # This test is identical to the one in `test_p4est_2d.jl`.
 #     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_free_stream.jl"),