Computational Acoustic Modeling Laboratory - projects:fem
http://caml.music.mcgill.ca/
2026-05-07T06:35:49+00:00Computational Acoustic Modeling Laboratory
http://caml.music.mcgill.ca/
http://caml.music.mcgill.ca/lib/exe/fetch.php?media=favicon.icotext/html2025-06-02T18:11:29+00:00Anonymous (anonymous@undisclosed.example.com)flangedpipe
http://caml.music.mcgill.ca/doku.php?id=projects:fem:flangedpipe&rev=1748887889&do=diff
FEM Simulation of a Flanged Pipe
%
% COMSOL / Matlab Script
% Calculation of the acoustic radiation impedance of a flanged pipe
% with the FEM (Helmholtz equation)
%
% Antoine Lefebvre, McGill University
%
% Copyright (c) 2010 Antoine Lefebvre
%
flclear fem
rho = 1.25;
c = 343;
a = 0.01;
L = 0.03;
N = 3; % ratio of sphere radius to cylinder radius
radiation_domain = ...
block3(2*N*a,2*N*a,N*a,'base','corner','pos',[-N*a,-N*a,0]) * sphere3(N*a);
model = (cylinder3(a,L,'pos',[0,0,-L]) + ra…text/html2025-06-02T18:11:29+00:00Anonymous (anonymous@undisclosed.example.com)toneholeflangeddalmont
http://caml.music.mcgill.ca/doku.php?id=projects:fem:toneholeflangeddalmont&rev=1748887889&do=diff
FEM Simulation of the Toneholes Studied by Dalmont et al. (2002)
%
% COMSOL / Matlab script
% FEM Simulation of two toneholes with the same geometry as in
% Dalmont et al. (2002), Experimental determination of the equivalent circuit
% of an open side hole: Linear and non linear behavior. Acustica, 88, 567-575.
%
% Antoine Lefebvre, McGill University
%
% Copyright (c) 2010 Antoine Lefebvre
%
flclear fem
% COMSOL version
clear vrsn
vrsn.name = 'COMSOL 3.5';
vrsn.ext = 'a';
vrsn.major = 0;
vrsn…text/html2025-06-02T18:11:29+00:00Anonymous (anonymous@undisclosed.example.com)unflangedpipe
http://caml.music.mcgill.ca/doku.php?id=projects:fem:unflangedpipe&rev=1748887889&do=diff
FEM Simulation of an Unflanged Pipe
%
% COMSOL / Matlab script
% Calculation of the acoustic radiation impedance of an unflanged pipe
% with the FEM (Helmholtz equation)
%
% Antoine Lefebvre, McGill University
%
% Copyright (c) 2010 Antoine Lefebvre
%
flclear fem
rho = 1.25;
c = 343;
a = 0.01;
L = 0.03;
N = 5;
ka = [0.01 colon(0.1,0.1,1.0)];
k = ka/a;
block = block3(2*N*a,2*N*a,2*N*a, 'pos',[0,0,-N*a]);
sphere = sphere3(N*a)*block;
pipe = cylinder3(a,L,'pos', [0,0,-L])*block;
radiation = s…