lisa.pizzo/LisaPizzoExercises
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Delete Project/CalculateLaplace_mult_rot.m

Browse source
This commit is contained in:
Lisa Pizzo 2026-01-26 20:38:22 +01:00
commit aee3215101
1 changed files with 0 additions and 40 deletions
Download patch file Download diff file Expand all files Collapse all files

40
Project/CalculateLaplace_mult_rot.m
View file

@ -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