Re: [Meep-discuss] to specify epsilon as an arbitrary function of?frequency

2012-02-16 Thread p . kinsler
hodjat hajian wrote: > Is it possible to specify epsilon as an arbitrary function of > frequency (not in Lorentzians or Drude form) in Meep? NO > If your answer is NO, please give me an alternative. > Thank you in advance for your time. You have to build a time domain response model with the

Re: [Meep-discuss] Source that only propagates in the positive?direction

2011-02-04 Thread p . kinsler
p.kins...@ic.ac.uk wrote: > You need a double source - essentially crossed electric > and magnetic dipoles (I posted something on this in > this ng sometime last year). Or drop a PML just behind > your ordinary source. http://comments.gmane.org/gmane.comp.science.electromagnetism.meep.general/35

Re: [Meep-discuss] Source that only propagates in the positive direction

2011-02-04 Thread p . kinsler
Huntoon, Nathan R wrote: > I'm trying to generate a launch source that only has a direction > of propogation in the positive Z dimension. I'm using the code > posted here on a guassian spatial profile and I specify the k > vector to be (0 0 1), but this launches fields in both the + > and - z

Re: [Meep-discuss] cylindrical FAIL at origin

2010-12-03 Thread p . kinsler
p.kins...@ic.ac.uk wrote: > Either both MEEP and my own simple FDTD code are wrong in the same > way, or an incoming cylindrical wave modelled in 1D cylindrical > coordinates (ie just r) behaves very oddly as it passes through the > origin. Having just done a Bessel function decomposition, I have

[Meep-discuss] cylindrical FAIL at origin

2010-11-30 Thread p . kinsler
Either both MEEP and my own simple FDTD code are wrong in the same way, or an incoming cylindrical wave modelled in 1D cylindrical coordinates (ie just r) behaves very oddly as it passes through the origin. I would expect that an incoming cylindrical wave (centred on r=0) would pass through the co

Re: [Meep-discuss] half wavelength

2010-10-15 Thread p . kinsler
??? wrote: > It wants to delay one source of wave of two source of waves in half > wavelength and to output it. I also want to have differently phased sources; is there a way of adding a phase factor or time-offset to sources? I can't find anything in the doco. -- ---

Re: [Meep-discuss] How to output the backwardly reflected field?pattern

2010-09-30 Thread p . kinsler
Using a linear array of point sources in 2D, the results can look pretty good: ; == ; == ; Directional source array ; ; Paul Kinsler / Imperial College London,

Re: [Meep-discuss] How to output the backwardly reflected field?pattern

2010-09-29 Thread p . kinsler
p.kins...@ic.ac.uk wrote: > Make a uni-directional light source. > [...] likewise if you want > to specify a direction the isn't along a lattice direction I think > (I haven't tried it yet) you'll get small, resolution dependent, > and remarkably persistent imperfections. Oops. Although it work

Re: [Meep-discuss] How to output the backwardly reflected field pattern

2010-09-29 Thread p . kinsler
Ziming Meng wrote: > I want to output the backwardly reflected field > pattern for normal incidence. However, it is usually superimposed > with the backward propagation of the light source. Then just using > the output-field function fails to do that. In calculating reflection > spect

Re: [Meep-discuss] Observations on dispersion, negative epsilon, ?and time averaged fields...

2010-08-22 Thread p . kinsler
Judson Ryckman wrote> When I use this definition, the simulation runs great.? However, > suppose my epsilon was now eps = -4.4 + 3.9605i.? I could arbitrarily, > change the epsilon term (eps_inf) from 0.3 to -0.2.? Strictly speaking, eps_inf should always be exactly 1 (the vacuum value); althou

Re: [Meep-discuss] about magnetostatic waves and tensors of mu

2010-08-05 Thread p . kinsler
Alexandr Sadovnikov wrote: > I tried to download the article which Dr. Paul Kinsler refer > to from the website [...] I suggest contacting the author directly -- Amiranashvili's email is visible at http://ito.wias-berlin.de/contact/staff/index.jsp?lang=1 In the other part of the thread,

Re: [Meep-discuss] about magnetostatic waves

2010-08-04 Thread p . kinsler
Alexandr Sadovnikov wrote: > And in this computational task one have to use > the relative permeability tensor (aka mu) and > the components of this tensor depend on the frequency. > And the main question is how can I do a material > with a the relative permeability tensor with the components >

Re: [Meep-discuss] 1D photonic band gap structure

2010-08-04 Thread p . kinsler
k...@vt.edu wrote: > I was reading "Photonic crystals: semiconductors of light" by Eli Yablonavitch > (Scientific American December 2001 ). I have a simple question about a figure > of 1D photonic crystal shown in that paper. > Please find the relevant figure here > http://filebox.vt.edu/users/kdas

Re: [Meep-discuss] Negative energy returned by field_energy_in_box

2009-11-19 Thread p . kinsler
Shavkat Nizamov wrote: > Surprisingly, some times the call of > field_energy_in_box(Ey, self.vol.surroundings()) > results in negative values. I wonder, what does this mean ? I suspect a non-optimal calculation of energy, see e.g. Power loss and electromagnetic energy density

Re: [Meep-discuss] Generating collapsing spherical wavefront

2009-10-15 Thread p . kinsler
ObsessiveMathsFreak wrote: > What I would like to do is generate a inwardly collapsing spherical wave > front, but meep does not seem to support sources with complex geometry. An inwardly collapsing spherical wave is a time-reversed outwardly expanding spherical wave. So why not generate an ex

Re: [Meep-discuss] magneto-electric materials

2009-08-21 Thread p . kinsler
Steven G. Johnson wrote: > On Aug 14, 2009, at 12:54 PM, p.kins...@ic.ac.uk wrote: > > I'd like to simulate magneto-electric materials (etc), > > i.e. where D depends on H, and B on E. > > As far as I can tell, meep doesn't do this. Would it > > be hard to implement? > Meep doesn't currently supp

[Meep-discuss] magneto-electric materials

2009-08-14 Thread p . kinsler
I'd like to simulate magneto-electric materials (etc), i.e. where D depends on H, and B on E. As far as I can tell, meep doesn't do this. Would it be hard to implement? -- -+- Dr. Paul Kinsler Blackett Laborator