diff --git a/applications/modules/multiphaseEuler/phaseSystem/phaseSystem/phaseSystemSolve.C b/applications/modules/multiphaseEuler/phaseSystem/phaseSystem/phaseSystemSolve.C
index d21965f54c..6718204b62 100644
--- a/applications/modules/multiphaseEuler/phaseSystem/phaseSystem/phaseSystemSolve.C
+++ b/applications/modules/multiphaseEuler/phaseSystem/phaseSystem/phaseSystemSolve.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2013-2023 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2013-2024 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -154,7 +154,7 @@ void Foam::phaseSystem::solve(const PtrList<volScalarField>& rAs)
                 (
                     IOobject
                     (
-                        "Sp",
+                        IOobject::groupName("Sp", phase.name()),
                         mesh_.time().name(),
                         mesh_
                     ),
@@ -168,7 +168,7 @@ void Foam::phaseSystem::solve(const PtrList<volScalarField>& rAs)
                 phasei,
                 new volScalarField::Internal
                 (
-                    "Su",
+                    IOobject::groupName("Su", phase.name()),
                     min(alpha.v(), scalar(1))
                    *fvc::div(fvc::absolute(phi_, phase.U()))->v()
                 )
diff --git a/applications/modules/multiphaseEuler/phaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C b/applications/modules/multiphaseEuler/phaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
index 9c6308af97..1e8f0c775f 100644
--- a/applications/modules/multiphaseEuler/phaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
+++ b/applications/modules/multiphaseEuler/phaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2015-2023 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2015-2024 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -708,7 +708,7 @@ void Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::invADVs
                     addField
                     (
                         i,
-                        "HVm",
+                        IOobject::groupName("HVm", phase.name()),
                         VmPhase*HDUDts[i],
                         HVms
                     );
@@ -726,7 +726,7 @@ void Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::invADVs
                     addField
                     (
                         i,
-                        "HVm",
+                        IOobject::groupName("HVm", phase.name()),
                         -VmPhase*HDUDts[j],
                         HVms
                     );
@@ -887,7 +887,7 @@ Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::invADVfs
                 addField
                 (
                     i,
-                    "HVmf",
+                    IOobject::groupName("HVmf", phase.name()),
                     VmPhasef
                    *byDt
                     (
@@ -911,7 +911,7 @@ Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::invADVfs
                     addField
                     (
                         i,
-                        "HVmf",
+                        IOobject::groupName("HVmf", phase.name()),
                        -VmPhasef
                        *byDt
                         (
@@ -1174,7 +1174,7 @@ void Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::dragCorrs
                 addField
                 (
                     i,
-                    "dragCorr",
+                    IOobject::groupName("dragCorr", phase.name()),
                     K1*(j == -1 ? -Uphis[i] : (Uphis[j] - Uphis[i])),
                     dragCorrs
                 );
@@ -1182,7 +1182,7 @@ void Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::dragCorrs
                 addField
                 (
                     i,
-                    "dragCorrf",
+                    IOobject::groupName("dragCorrf", phase.name()),
                     fvc::interpolate(K1)
                    *(j == -1 ? -phase.phi() : (otherPhase.phi() - phase.phi())),
                     dragCorrfs
diff --git a/applications/modules/multiphaseEuler/populationBalance/populationBalanceModel/populationBalanceModel.C b/applications/modules/multiphaseEuler/populationBalance/populationBalanceModel/populationBalanceModel.C
index d3766298ea..9435b9d415 100644
--- a/applications/modules/multiphaseEuler/populationBalance/populationBalanceModel/populationBalanceModel.C
+++ b/applications/modules/multiphaseEuler/populationBalance/populationBalanceModel/populationBalanceModel.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2017-2023 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2017-2024 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -843,7 +843,12 @@ Foam::diameterModels::populationBalanceModel::populationBalanceModel
             i,
             new volScalarField
             (
-                IOobject("Su", fluid_.time().name(), mesh_),
+                IOobject
+                (
+                    IOobject::groupName("Su" + Foam::name(i), this->name()),
+                    fluid_.time().name(),
+                    mesh_
+                ),
                 mesh_,
                 dimensionedScalar(inv(dimTime), 0)
             )
@@ -854,7 +859,12 @@ Foam::diameterModels::populationBalanceModel::populationBalanceModel
             i,
             new volScalarField
             (
-                IOobject("Sp", fluid_.time().name(), mesh_),
+                IOobject
+                (
+                    IOobject::groupName("Sp" + Foam::name(i), this->name()),
+                    fluid_.time().name(),
+                    mesh_
+                ),
                 mesh_,
                 dimensionedScalar(inv(dimTime), 0)
             )
@@ -871,7 +881,7 @@ Foam::diameterModels::populationBalanceModel::populationBalanceModel
             (
                 IOobject
                 (
-                     "coalescenceRate",
+                     IOobject::groupName("coalescenceRate", this->name()),
                      mesh_.time().name(),
                      mesh_
                 ),
@@ -897,7 +907,7 @@ Foam::diameterModels::populationBalanceModel::populationBalanceModel
             (
                 IOobject
                 (
-                    "breakupRate",
+                    IOobject::groupName("breakupRate", this->name()),
                     fluid_.time().name(),
                     mesh_
                 ),
@@ -915,7 +925,7 @@ Foam::diameterModels::populationBalanceModel::populationBalanceModel
             (
                 IOobject
                 (
-                    "binaryBreakupRate",
+                    IOobject::groupName("binaryBreakupRate", this->name()),
                     fluid_.time().name(),
                     mesh_
                 ),
@@ -951,7 +961,7 @@ Foam::diameterModels::populationBalanceModel::populationBalanceModel
             (
                 IOobject
                 (
-                    "driftRate",
+                    IOobject::groupName("driftRate", this->name()),
                     fluid_.time().name(),
                     mesh_
                 ),
@@ -969,7 +979,7 @@ Foam::diameterModels::populationBalanceModel::populationBalanceModel
             (
                 IOobject
                 (
-                    "nucleationRate",
+                    IOobject::groupName("nucleationRate", this->name()),
                     fluid_.time().name(),
                     mesh_
                 ),
diff --git a/src/finiteVolume/fvMesh/wallDist/wallDist/wallDist.C b/src/finiteVolume/fvMesh/wallDist/wallDist/wallDist.C
index a75c16dc1b..be2fecc9a4 100644
--- a/src/finiteVolume/fvMesh/wallDist/wallDist/wallDist.C
+++ b/src/finiteVolume/fvMesh/wallDist/wallDist/wallDist.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2015-2023 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2015-2024 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -74,8 +74,7 @@ Foam::wallDist::wallDist(const fvMesh& mesh, const word& patchTypeName)
     (
         patchDistMethod::New
         (
-            static_cast<const fvSchemes&>(mesh)
-           .subDict(patchTypeName_ & "Dist"),
+            mesh.schemes().subDict(patchTypeName_ & "Dist"),
             mesh,
             patchIndices_
         )
@@ -94,7 +93,9 @@ Foam::wallDist::wallDist(const fvMesh& mesh, const word& patchTypeName)
     ),
     nRequired_
     (
-        static_cast<const fvSchemes&>(mesh).subDict(patchTypeName_ & "Dist")
+        mesh
+       .schemes()
+       .subDict(patchTypeName_ & "Dist")
        .lookupOrDefault<Switch>("nRequired", false)
     )
 {
@@ -124,8 +125,7 @@ Foam::wallDist::wallDist
     (
         patchDistMethod::New
         (
-            static_cast<const fvSchemes&>(mesh)
-           .subDict(patchTypeName_ & "Dist"),
+            mesh.schemes().subDict(patchTypeName_ & "Dist"),
             mesh,
             patchIndices_
         )
@@ -144,7 +144,9 @@ Foam::wallDist::wallDist
     ),
     nRequired_
     (
-        static_cast<const fvSchemes&>(mesh).subDict(patchTypeName_ & "Dist")
+        mesh
+       .schemes()
+       .subDict(patchTypeName_ & "Dist")
        .lookupOrDefault<Switch>("nRequired", false)
     )
 {