Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

eigenmodes of an acoustic shoebox using BEM in COMSOL

Please login with a confirmed email address before reporting spam

Hi all,

Has anyone ever coded a BEM formulation for the eigen frequencies and eigen modes of a shoe box ?
I'm trying, but it doesn't seem to be that easy...

Basically, I want to solve _only on the boundary_ of a 3D domain, the following equation:
Int(dG/dn(M,M').p(M').dS) + p(M)/2 = Int(G(M,M').dp/dn(M').dS)
subject to the generalized boundary condition :
i.k.rho.c.p(M) + Z*dp/dn(M) = 0
Notations:
- M is any point of the boundary
- M' is another point that lies on an elementary patch dS of the boundary, and the integration is performed such that M' (and associated surface dS) sweeps the entire boundary.
- p(M) is the pressure in M
- dp/dn(M) is the gradient of p in M normal to the boundary in M
- G(M,M') is Green's function, that is exp(i.k.r)/(4.Pi.r), where r is ||vect(MM')||
- dG/dn(M,M') = is the derivative of G(M,M') along the normal to the boundary in M'. It is computed as the derivative of G along vect(MM'), multiplied by the projection of the unit vector of that direction on the normal axis to the boundary at M'
- Z is the specific acoustic impedance of the surface

For coding the integral equations, I was inspired by this post :
www.comsol.com/community/forums/general/thread/33357

except that I used the "boundary ode" physics node, stting the time derivative coefficients to 0.

My problem : at compilation, the model fails to evaluate some expressions (the error message is rather hard to decrypt...)

I attached my example.
Is anyone geek enough to look into this ?



NB : For deriving the eauqtions, I used the direct formulation given in Kirkup's book "The BEM in Acoustics" (2007 ed., available in PDF online, linked to the website www.boundary-element-method.com ). Relevant formulas are found pp. 9, 32-35, and 107-110 of the book.



1 Reply Last Post 2 nov 2015, 02:30 GMT-5
COMSOL Moderator

Hello Sami Karkar

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.


Please login with a confirmed email address before reporting spam

Posted: 9 years ago 2 nov 2015, 02:30 GMT-5
hello,

i saw your topic in this webpage, and recently, i am working in the problem about acoustic eigenvalues Problems with the impedance boundary condition using BEM. Maybe be we can discuss some, but i use this forum in the first time, so we could send emails each other and my address is 18392622690@163.com.

hope discuss well later.
Meng, ss
hello, i saw your topic in this webpage, and recently, i am working in the problem about acoustic eigenvalues Problems with the impedance boundary condition using BEM. Maybe be we can discuss some, but i use this forum in the first time, so we could send emails each other and my address is 18392622690@163.com. hope discuss well later. Meng, ss

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.