From aee3215101013494a31ec7372ef65a53c2991b63 Mon Sep 17 00:00:00 2001 From: "lisa.pizzo" Date: Mon, 26 Jan 2026 20:38:22 +0100 Subject: [PATCH] Delete Project/CalculateLaplace_mult_rot.m --- Project/CalculateLaplace_mult_rot.m | 40 ----------------------------- 1 file changed, 40 deletions(-) delete mode 100644 Project/CalculateLaplace_mult_rot.m diff --git a/Project/CalculateLaplace_mult_rot.m b/Project/CalculateLaplace_mult_rot.m deleted file mode 100644 index a205c5b..0000000 --- a/Project/CalculateLaplace_mult_rot.m +++ /dev/null @@ -1,40 +0,0 @@ -function [K, F] = CalculateLaplace_mult_rot(model, lambda_wall, lambda_fluid, lambda_air) -mesh = model.Mesh; -nodes = mesh.Nodes; -elements = mesh.Elements; - -Nnodes = size(nodes,2); -Nelems = size(elements,2); - -K = sparse(Nnodes, Nnodes); -F = zeros(Nnodes,1); - -regions = zeros(Nelems,1); -regions(findElements(mesh,'region','Face',1)) = 1; -regions(findElements(mesh,'region','Face',2)) = 2; -regions(findElements(mesh,'region','Face',3)) = 3; - -for e = 1:Nelems - vert = elements(:,e); - - x = nodes(1,vert); % r-coordinates - y = nodes(2,vert); % z-coordinates - - Ae = polyarea(x,y); - - b = [y(2)-y(3); y(3)-y(1); y(1)-y(2)]; - c = [x(3)-x(2); x(1)-x(3); x(2)-x(1)]; - - rbar = mean(x); % <-- axisymmetric weight - switch regions(e) - case 1 - lambda = lambda_wall; - case 2 - lambda = lambda_fluid; - case 3 - lambda = lambda_air; - end - Ke = rbar * (lambda/(4*Ae)) * (b*b.' + c*c.'); - K(vert,vert) = K(vert,vert) + Ke; -end -end