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

Delete Project/CalculateLaplace_mult.m

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

45
Project/CalculateLaplace_mult.m
View file

@ -1,45 +0,0 @@
function [K, F] = CalculateLaplace_mult(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);
% --- build region array manually ---
regions = zeros(Nelems,1);
wallElems = findElements(mesh,'region','Face',1);
fluidElems = findElements(mesh,'region','Face',2);
airElems = findElements(mesh,'region','Face',3);
regions(wallElems) = 1;
regions(fluidElems) = 2;
regions(airElems) = 3;
for e = 1:Nelems
vert = elements(:,e);
x = nodes(1,vert);
y = nodes(2,vert);
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)];
% conductivity by material
switch regions(e)
case 1
lambda = lambda_wall;
case 2
lambda = lambda_fluid;
case 3
lambda = lambda_air;
otherwise
error('Element %d has no material assignment', e);
end
Ke = (lambda/(4*Ae)) * (b*b.' + c*c.');
K(vert,vert) = K(vert,vert) + Ke;
end
end