Hello there!
For my bachelor’s, I was trying to find examples of simulating different antennas using bempp or bem++. I found that code examples presented in the following paper Solving Boundary Integral Problems with BEM++ https://dl.acm.org/doi/10.1145/2590830 can not be seen or downloaded anywhere on the web. I would like to get the “Exterior Dirichlet Problem” and “Scattering by a Screen” examples as presented in this paper.
Any help is greatly appreciated!
Best, Anthony.
I would greatly appreciate the code examples too, especially bolt.py which demonstrates a robin boundary implementation.
Hi @AnthonyTSV and @stuvx!
I had a copy of the example scripts from that paper, I’ve uploaded them here: https://www.dropbox.com/scl/fi/bdgciou6y1eqc98d7zlth/bempp-toms-examples.zip?rlkey=1sofquslquret6sh8wxtgk2q6&dl=0
Note that these scripts are unlikely to work any more unless you install a very old version of Bempp-legacy.
I’ll look into getting these hosted somewhere, and will then link to them from the web page, as even if they no longer work, I think it’s important that they’re available
Thanks @mscroggs - I was hoping it would shed some light on how to implement Robin boundary conditions, but I am just as confused as before
I’m not certain what the “standard” way of imposing Robin conditions is (perhaps someone else can suggest something), but you may be interested in the method of weakly imposing Robin conditions that we proposed in section 4.3 of Boundary element methods with weakly imposed boundary conditions (Betcke, Burman, Scroggs, 2019) (https://doi.org/10.1137/18M119625X, pdf at https://www.mscroggs.co.uk/papers/2019-weaklyimposed.pdf if you don’t have access).
If I remember correctly, the formulation we ended up with was tidier than many other formulations for Robin conditions that we’d seen.
Thank you - I will take a look. I am a novice in BEM and am trying to use bempp-cl replicate this paper: https://scholarcommons.sc.edu/cgi/viewcontent.cgi?article=1074&context=eche_facpub
For the Robin boundary condition, I am not sure how to obtain phi: does it get passed in as the solution to greens equation, or do I calculate it from the point coordinate? Once I have phi there is the left hand side of the BC: the derivative of phi wrt. the normal direction; and the right hand side: f(phi) = i, so do I need to define all of these terms in my RBC?