Dear Steven and Meep users, I met a problem of obainning the steady state response of a cw source. For 2D case, since the data storage size is relatively small, I used the method that Steven suggested.
Quote from Steven: (To others on this list: "steady-state" response is the usual term for the exp(-iwt) response field from an exp(-iwt) source, after all transients have died away.) To get the steady-state response from a CW source, currently the best way is to simply run for a long time after smoothly turning on a source (you will want to set the fwidth or width parameter of the continuous-src to make it turn on smoothly instead of suddenly). End of Quote. http://www.mail-archive.com/[email protected]/msg00040.html For example, my strucutre is 2D complex geometry. CW locates at the -x side and wave proprogates to +x side. I want to obtain the steady state response at certain x point along y direction. I just record the time response after turning on CW source for a long time for several period at the end of the running steps into a big hdf5 file. I assume wavefront in y direction having relatively uniform phase. After loading the file into matlab, and find out the maximum value for each postion point. That is how I find out the steady state. The problem now, is, for 3D, this hdf5 is really large and unrealistic to process in matlab. Is there any easier way to solve it efficiently? One way as indicated by Steven is to find out the time point where J term has zero phase by counting integer number of periods, and record the field at that time point. For period counting method, it can only garantee the zero phase at the source. But for abitrary location, this might not be true, since the phase also changes as function of position. How to find the zero phase for arbitrary position? Hope some one can give me some guidance. Thank you. ? Lingyun
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
