Re: [Meep-discuss] Field Not Converge with GaussianBeamSource

2020-11-01 Thread Ardavan Oskooi
Perhaps you are combining two or more mirror symmetries with periodic boundary conditions in 3d? This may cause field instabilities as described in: https://github.com/NanoComp/meep/issues/132. If this is the case, try removing the mirror symmetries. On 11/1/2020 13:28, Gui Max wrote:

[Meep-discuss] Field Not Converge with GaussianBeamSource

2020-11-01 Thread Gui Max
Hello I am trying to mimic a terahertz pulse (about 910GHz) signal measured from a laser spectroscopy machine. The time-domain waveform has a unique shape so I have to use CustomSource (i.e. sources = [mp.GaussianBeamSource(src=mp.CustomSource). I obtained the expression using MATLAB curve

[Meep-discuss] Field of antenna as a source

2017-06-27 Thread Jitendra
Hello users, I am a final year masters student. In an experiment, I am using a small-loop-antenna to excite field in a Loop-Gap-Resonator. The field of such antenna is given by the expression: E(r,theta)= exp(-ikr)sin(theta)/r This field is in the phi direction and other components are zero.

Re: [Meep-discuss] Field divergence in coordinate transforming (stretching)

2013-11-08 Thread Tran Quyet Thang
tran quyet thang hereiam2005@... writes: Dear professor Johnson and Meep users, It has been mentioned multiple times that a capacity similar to variable mesh size can be emulated in meep by applying the coordinate transformation technique formulated in Coordinate Transformation

[Meep-discuss] Field divergence in coordinate transforming (stretching)

2013-11-07 Thread tran quyet thang
Dear professor Johnson and Meep users, It has been mentioned multiple times that a capacity similar to variable mesh size can be emulated in meep by applying the coordinate transformation technique formulated in Coordinate Transformation Invariance in Electromagnetism. However my attempts in

[Meep-discuss] field components and propagation direction of a volume source

2012-05-28 Thread ranjeetd
dear meep users suppose i excite the volume source with Ez component with unit length in y direction then does it have cylindrical wavefront? which component of electric and magnetic fields will be present? ranjeet dwivedi IIT KANPUR ___

[Meep-discuss] field strength question

2012-01-18 Thread Marcel Spurny
Dear all, I want to calculate the intensity inside of a phc cavity. Does anyone know how a gaussian point source e-field is normalized? My idea is to normalize the field strength inside of the cavity to the source and then scale up according to the input laser. Any other ideas are welcome.

Re: [Meep-discuss] Field enhancement at dielectric tip

2011-08-03 Thread gpipc
On Tue, 2 Aug 2011 08:49:47 -0400, J. Steele Parris parri...@vcu.edu wrote: The tip is a 1x1x4 nickel block that has been cut to a pyramid point on one side and affixed to a nickel layer (1 unit thick) and a silicon substrate. There is also a thin layer of silver on the pyramid point of the

Re: [Meep-discuss] Field enhancement at dielectric tip

2011-08-03 Thread J. Steele Parris
Okay, I feel that I have solved the issue. It was multi-faceted: 1) When I was probing the h5 file with h5totxt, I was reading data from a wrong point (should've been at -x60 -y11 -z11, not -y0 -z0). 2) Also, when using output-efield, there are 3 different datasets (x, y, z) which the h5totxt was

Re: [Meep-discuss] Field enhancement at dielectric tip

2011-08-02 Thread J. Steele Parris
The tip is a 1x1x4 nickel block that has been cut to a pyramid point on one side and affixed to a nickel layer (1 unit thick) and a silicon substrate. There is also a thin layer of silver on the pyramid point of the tip. On Mon, Aug 1, 2011 at 6:26 PM, gpipc gp...@cup.uni-muenchen.de wrote: On

[Meep-discuss] Field enhancement at dielectric tip

2011-08-01 Thread J. Steele Parris
Hello fellow Meep users! I am measuring the enhancement of the electric field at the tip of a nickel nanotip with a silver coating on top. When I measure and compare to a blank field, the enhancement is 10^13. I find this improbable, and I ask the community for any advice on the matter. To

Re: [Meep-discuss] Field enhancement at dielectric tip

2011-08-01 Thread gpipc
On Mon, 1 Aug 2011 14:51:38 -0400, J. Steele Parris parri...@vcu.edu wrote: Hello fellow Meep users! I am measuring the enhancement of the electric field at the tip of a nickel nanotip with a silver coating on top. When I measure and compare to a blank field, the enhancement is 10^13. I find

Re: [Meep-discuss] field synchronization

2010-12-26 Thread Özgür Önder KARAKILINÇ
Thanks for your help Onder Message: 1 Date: Fri, 24 Dec 2010 16:13:19 -0500 From: Steven G. Johnson stevenj@gmail.com Subject: Re: [Meep-discuss] field synchronization To: meep-discuss Discuss meep-discuss@ab-initio.mit.edu Message-ID: 24baab40-fa99-4af7-9122-b5361d1a4...@fftw.org Content

[Meep-discuss] field synchronization

2010-12-24 Thread Özgür Önder KARAKILINÇ
Hi meep users, I have two sources in my structure. One of them is Ez field, the other one is Hz field. I want to get fluxes with add-flux function. In wiki page, fields are synchronized automatically in flux-in-box and field-energy-in-box routines. But I am using add-flux function to get fluxes.

Re: [Meep-discuss] field synchronization

2010-12-24 Thread Steven G. Johnson
On Dec 24, 2010, at 3:40 AM, Özgür Önder KARAKILINÇ wrote: So, should I do synchronization of fields when I use add-flux function ? No, the Fourier transforms in add-flux automatically multiply the fields by the phase for the correct times; no synchronization is needed.

Re: [Meep-discuss] field profile along the axis

2010-09-23 Thread Steven G. Johnson
On Sep 23, 2010, at 6:56 AM, Özgür Önder KARAKILINÇ wrote: I have 160x160 matrix dataset I think, row number indicate x-axis direction, column number indicate y-axis direction Is it right? Yes. (You could easily check this yourself just by making a non- square computational cell.)

[Meep-discuss] field profile along the axis

2010-09-16 Thread Özgür Önder KARAKILINÇ
I want to get electric field intensity profile along the waveguide axis (for example each point in x direction ) How can I do? in meep reference, (get-field-point c pt) function is used for single point. I need field intensity-position graphs. Maybe, I should write script to get intensity at

Re: [Meep-discuss] field profile along the axis

2010-09-16 Thread Steven G. Johnson
On Sep 16, 2010, at 4:14 AM, Özgür Önder KARAKILINÇ wrote: I want to get electric field intensity profile along the waveguide axis (for example each point in x direction ) How can I do? in meep reference, (get-field-point c pt) function is used for single point. I need field

[Meep-discuss] Field enhancement calculation

2010-07-14 Thread Carlos Aspetti
Hello Everyone, I would like to compute the near field enhancement within a plasmonic cavity (acted on by a continuous wave source). When outputting a component of the electric field (ie. E(r)) is this simply the amplitude corresponding to the real part of the electric field in arbitrary units? I

[Meep-discuss] Field does not converge

2009-07-20 Thread ryanhao
Dear Steven and other Friends, I met a problem. I calculate the transmittance of W1 PhC waveguide. I randomly change the radius of each air hole. To varify the correction, I have added a monitor at a single point to observe the feild at the output. I

[Meep-discuss] field blow up at last

2009-06-21 Thread phcgroup
Dear all, I simulate a structure made from dielectric (without active material). I use a pulse and put PML at all sides, but the field blew up at last. (set! eps-averaging? false) can't help. The version of meep is 0.10.1. Could you help me? Thank you in advance. Li

Re: [Meep-discuss] field blow up at last

2009-06-21 Thread Steven G. Johnson
On Jun 21, 2009, at 8:58 AM, phcgroup wrote: I simulate a structure made from dielectric (without active material). I use a pulse and put PML at all sides, but the field blew up at last. (set! eps-averaging? false) can't help. The version of meep is 0.10.1. Have you tried a more recent

[Meep-discuss] FIELD

2009-03-29 Thread Michal Martinů
Please, I need good advice. Let me to show my problem closer. You can see (enclosed), how my initial magnetic field from static analysis looks. It is non-continous, because the current source is present in it. I can get the magnetic intensity Hx, Hy values only in nodes. I have regular mesh

[Meep-discuss] FIELD

2009-01-26 Thread Michal Martinů
Hello, please help me. I would like to start my Meep computing with initialized magnetic field and I do not know, how to do it. I do not understand this: http://ab-initio.mit.edu/wiki/index.php/Meep_field-function_examples I can translate Hx, Hy (intenzity of magnetic field) from static 2D

Re: [Meep-discuss] Field divergence

2008-11-16 Thread jinkyu . yang
I also met this problem before. If your structure is so large and calculation time is too long, you should activate non-averaging to avoid divergence issue just like; (set! eps-averaging? false) Quoting dy zhao [EMAIL PROTECTED]: Dear Dr. Steven and meep users, Recently, I have met the

[Meep-discuss] Field blowing out

2008-09-22 Thread Seong Kyu Kim
Dear prof. Steven Johnson and meep users: I am reporting a problem of field blowing out. This may be a similar problem to the reports of Andika Asyuda in September or Jefferson Thomas in July. As I do not see their solutions, this problem may be a limitation of current meep applications. I

Re: [Meep-discuss] Field divergence in meep after many steps for a photonic crystal

2008-07-31 Thread Jefferson Thomas
Thanks. One more thing about adding the sources in case of symmetries : If I have a box of sizes 4, 4, 4 with odd symmetry on X and I want to add a source with Ex with f.add_point_source(Ex, src, vec(1, 1, 1), 10); where 10 means the amplitude then as far as I understand I have to add a

Re: [Meep-discuss] Field divergence in meep after many steps for a photonic crystal

2008-07-30 Thread Jefferson Thomas
Hi Steven, What you are saying is really very helpful. In fact I had no idea I have to supply the symmetric source. I will try it right now and see how it goes. On the other hand I dont understand the thing with the symmetric eps. Although I dont know how meep works internally with the

Re: [Meep-discuss] Field divergence in meep after many steps for a photonic crystal

2008-07-30 Thread Steven G. Johnson
On Jul 30, 2008, at 9:20 PM, Jefferson Thomas wrote: What you are saying is really very helpful. In fact I had no idea I have to supply the symmetric source. I will try it right now and see how it goes. See: http://ab-initio.mit.edu/wiki/index.php/Exploiting_symmetry_in_Meep If you don't

Re: [Meep-discuss] Field divergence in meep after many steps for a photonic crystal

2008-07-30 Thread Jefferson Thomas
You are right. When I symmetrized the eps function the divergences disappeared. I really didnt know that this one last pixel could have such a big influence. I have also some other kind of divergences in other cases. They seem not to be connected to symmetries but arise from PML layers. I am

Re: [Meep-discuss] Field divergence in meep after many steps for a photonic crystal

2008-07-29 Thread Steven G. Johnson
Your problem seems to stem from the fact that neither your structure nor your sources obey the symmetry you specified. On Jul 22, 2008, at 9:40 PM, Jefferson Thomas wrote: meep::symmetry S=meep::identity(); S = S+ meep::mirror(meep::X, v)*(-1); /// Introducing symmetries on X and Y

Re: [Meep-discuss] Field divergence in meep after many steps for a photonic crystal

2008-07-19 Thread Steven G. Johnson
On Jul 16, 2008, at 10:11 PM, Jefferson Thomas wrote: - exciting with one gaussian_src_time with freq 0.7, width 3.0 for the Try using a narrower-bandwidth source centered on your mode frequency. That will also help harminv out, in that with a narrower bandwidth source harminv has cleaner

Re: [Meep-discuss] field reaching infinite value in the beginning

2008-03-14 Thread Steven G. Johnson
On Mar 14, 2008, at 6:04 AM, shiv chawla wrote: (set! Courant 0.95) In two dimensions, the Courant factor needs to be sqrt(1/2) for the simulation to be stable, as described in the introduction chapter of the Meep manual. ___ meep-discuss

Re: [Meep-discuss] Field should decay but blows up.

2007-08-10 Thread Ruslan Sepkhanov
On 17/07/07, Steven G. Johnson [EMAIL PROTECTED] wrote: till some value and then blow up. I tried to experiment a bit, and I have figured out that the process can be made stable by disabling subpixel averaging (I hope the following table could be usefull): Yes, right now the subpixel

[Meep-discuss] Field should decay but blows up.

2007-07-07 Thread Ruslan Sepkhanov
Dear Tjeerd, Steven and other meep users, I am calculating a transmission spectrum through 2D photonic crystal and I expirienced with the same problem. For some parameters and when the structure is present I cannot get the fields to decay, while normalization run is OK. If length L is set to 16

[Meep-discuss] Field should decay but blows up.

2007-07-03 Thread Tjeerd Pinkert
Hello all, I'm simulating PhC resonators by launching pulses into the system. Like in the bend waveguide example I put a field measurement at the output to make sure that the fields at the output are decayed enough. I put the decay value to 1e-3. The fields generally tend to decay first, but

[Meep-discuss] field blowup problem.

2006-09-14 Thread K. Choi
Hi everyone. I found a strange problem with metal-dielectric structure. On my structure ctl file, the surface field blowup in pml. I set the frequency 0.55 and Wp = 1. Meep version is 0.9.== (set! geometry-lattice (make lattice (size 8 9 7)))(set! geometry (list (make block (center 0 0 0)