Dear Developers, I would like to be able to carry out FDTD simulations on a flat (Euclidean) region of space that is finite in extent but boundless. An example of such a space would be the 2D surface of a disc, with the outer circumferential ring of the disc connected spatially in software to the outer ring of a second 2D disc of equal radius such that waves flowing outward from the first disc would smoothly cross the outer ring (and then flow inward on the second disc). Therefore, there would be no PML, absorbing, or periodic boundary conditions (there would be no boundaries at all), and the total area of the FDTD universe would be 2*pi*r^2. I would like to then set up a gradient index profile on both disc surfaces and see how waves propagate in this space.
Another example would be to carry out 2D simulations on the closed surface of a cube. It seems like because what I want to do still involves Euclidean geometry, nothing would change in the FDTD code other than being able to "connect" non-contiguous regions of space together such that waves would propagate smoothly across the stitched regions, as they should if the index is continuous. I tried to think of a way to do this with the current meep code but failed. I thought maybe I could watch fields propagate to a boundary, record somehow the fields passing the boundary, and then replay those fields in a second simulation with a custom source, but this approach failed because of a few reasons, not to mention that this would be really complicated. How difficult would it be to allow "boundary conditions" that simply connect to other structures? Yours, Aaron Danner _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
