From 358bc227a3d0c07b871f8e44afbfd29afbe98a44 Mon Sep 17 00:00:00 2001
From: "lisa.pizzo" <lisa.pizzo@noreply.imsc.uni-graz.at>
Date: Mon, 26 Jan 2026 20:38:07 +0100
Subject: [PATCH] Delete Project/AddMass_mult_rot.m

---
 Project/AddMass_mult_rot.m | 35 -----------------------------------
 1 file changed, 35 deletions(-)
 delete mode 100644 Project/AddMass_mult_rot.m

diff --git a/Project/AddMass_mult_rot.m b/Project/AddMass_mult_rot.m
deleted file mode 100644
index 0a0736b..0000000
--- a/Project/AddMass_mult_rot.m
+++ /dev/null
@@ -1,35 +0,0 @@
-function M = AddMass_mult_rot(model, M, c_wall, c_fluid, c_air)
-mesh = model.Mesh;
-nodes = mesh.Nodes;
-elements = mesh.Elements;
-
-Nelems = size(elements,2);
-
-regions = zeros(Nelems,1);
-regions(findElements(mesh,'region','Face',1)) = 1; % wall
-regions(findElements(mesh,'region','Face',2)) = 2; % fluid
-regions(findElements(mesh,'region','Face',3)) = 3; % air
-
-for e = 1:Nelems
-    vert = elements(:,e);
-
-    x = nodes(1,vert);   % r-coordinates
-    y = nodes(2,vert);   % z-coordinates
-
-    Ae = polyarea(x,y);  % element area
-    rbar = mean(x);      % mean radius (axisymmetric weight)
-
-    switch regions(e) % Select heat capacity
-        case 1
-            c = c_wall;
-        case 2
-            c = c_fluid;
-        case 3
-            c = c_air;
-    end
-    % Axisymmetric element mass matrix
-    Me = rbar * c * Ae / 12 * [2 1 1; 1 2 1; 1 1 2];
-    % Assemble
-    M(vert,vert) = M(vert,vert) + Me;
-end
-end