Author |
Message |
Jamy (neo)
Member Username: neo
Post Number: 5 Registered: 08-2006
| Posted on Wednesday, September 20, 2006 - 07:45 am: | |
Hi! Please find attached the script and the error message. It says that an internal processing error has occurred. I have not been able to resolve the problem. Please advice JM 215.7 K | Error msg.bmp "error msg" | | 4.5 K | MB.pde "script" | | |
Jamy (neo)
Member Username: neo
Post Number: 6 Registered: 08-2006
| Posted on Wednesday, September 20, 2006 - 03:59 pm: | |
It looks like the attachmets have not been uploaded properly. So I am attaching them again. Also I have noticed that I do not get the error if I use an nodelimit of 512. Please look into this probleam. Why can't I use higher values for nodelimit. With the above script I get the error after 75 minutes of running teh script in FlexPDE 5. Please advice. Regards Jamy 4.5 K | MB.pde "PDE script" | | 215.7 K | error messege.bmp "error msg" | | |
Jamy (neo)
Member Username: neo
Post Number: 7 Registered: 08-2006
| Posted on Wednesday, September 20, 2006 - 04:27 pm: | |
Attachment not working for me: So I am copying the code below. TITLE 'Spatial mode dynamics: Wide aperture edge emitting laser diode' COORDINATES CARTESIAN1 VARIABLES Er ! Real Part of Complex Amplitude {rescalsed} Eim ! Imaginary Part of Complex Amplitude {rescaled} Pr ! Real part of Complex polarization {rescaled} Pim ! Imaginary part of Complex polarization {rescaled} N ! Carrier Density {rescaled} SELECT nodelimit = 1024 !errlim = 1e-4 changelim = 0.01 DEFINITIONS Jop = 1.5e-5 b1 = -50 b2 = 50 ! Current Density Profile !J = (exp(-(-x-50)/14))*Jop*(if x>=-256 and x<b1>=b1 and x<=b2>b2 and x<=255 then 1 else 0) J = Jop*MIN(exp((50-abs(x))/14), 1) ! -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- -------------------- c_zero = 3e14 ! Velocity of light in vaccum (micron-second^-1) Diffcoef = 3.3e9 ! Diffusion Coefficient (micron^2*second^-1) qe = 1.6e-19 ! lementary electronic charge junc_thick = 60e-4 ! Active Layer Thickness (microns) gamma_E = 0.723e12 ! Field damping rate (second-1) gamma_P = 1e13 ! Polarization damping rate (second-1) gamma_N = 1/5e-9 ! Carrier Decay rate {inverse of non-radiative recombination lifetime} (pico-s-1) N_tr = 1.3e6 ! Transparency current density (microns^3) eta = 3.5 ! Junction refractive index alpha_H = 3 ! Linewidth enhancement (Henry) Factor !a_zero = 2e-6 ! Differential gain coefficient {rescaled} (microns^3*pico-second^ -1) !omega = 2.2987e3 ! Angular frequency (radians-pico-s-1) omega=2*pi*c_zero/0.8 ! Angular frequency (radians-s-1) a_zero = 2e6 ! Differential gain coefficient {rescaled} (microns^3*second^ -1) ! -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- -------------------- vs=10e12 !Surface recombination verlocity at the wings absw=30e-4 ! Absoption coefficeint at the wings gain = a_zero*(N-N_tr) ! Gain coefficient Int=Er^2+Eim^2 ! Intensity Int2=Int*100 ! Scaled intensity A = (c_zero*c_zero)/(omega*eta*eta) B = 1-(gamma_E/(2*gamma_P)) L = (1+ (alpha_H*alpha_H)) !------------------------------------------------------------------------------- -------------------------------------------------------------------------------- ---------------------- INITIAL VALUES N=5e5*MIN(exp((50-abs(x))/14), 1) !exp(-x^2/35^2) Er=1e-3*exp(-(abs(x)/(1.4142*35))^7)!*cos(10*x) Eim=1e-3*exp(-(abs(x)/(1.4142*35))^7)!*cos(10*x) Pr=1e4*exp(-(abs(x)/(1.4142*35))^7)!*cos(10*x) Pim=1e4*exp(-(abs(x)/(1.4142*35))^7)!*cos(10*x) EQUATIONS Er: 2*dt(Er) = -gamma_E*(Er - (alpha_H*Eim)) + A*dxx(Eim) + Pr Eim: 2*dt(Eim) = -gamma_E*(Eim + (alpha_H*Er)) - A*dxx(Er) + Pim Pr: dt(Pr) = - gamma_P*((Pr + (alpha_H*B*Pim)) - L*gain*Er) Pim: dt(Pim) = - gamma_P*((Pim - (alpha_H*B*Pr)) - L*gain*Eim) N: dt(N) = (J/(qe*junc_thick))-(gamma_N*N)+(Diffcoef*dxx(N)) - (Er*Pr + Eim*Pim) BOUNDARIES { The domain definition } REGION 1 { For each material region } START (-256) point load(Er)=absw*Er point load(Eim)=absw*Eim point value(Pr)=0 point value(Pim)=0 point load(N)=vs*N/Diffcoef LINE TO (-50) LINE TO (50) LINE TO(255) point load(Er)=-absw*Er point load(Eim)=-absw*Eim point value(Pr)=0 point value(Pim)=0 point load(N)=-vs*N/Diffcoef TIME 0 by 1e-17 to 50e-9 !MONITORS PLOTS { save result displays } for t = 0 by 1e-12 to 50e-9 elevation(50*Int,N,J*1e11) from (-150) to (150) as 'INTENSITY vs CARRIER DENSITY [fixed current density]' elevation(Pr,Pim) from (-256) to (255) as 'POLARIZATION [lateral variation]' Histories history(Pr,Pim) at (0) as "POLARIZATION" ! "RELAXATION OSCILLATION" history(10*N,Int2) at (0) as "RELAXATION OSCIALLTION" END |
Jamy (neo)
Member Username: neo
Post Number: 8 Registered: 08-2006
| Posted on Thursday, September 21, 2006 - 07:08 am: | |
the error msg is attached here 215.7 K | error messege.bmp "err msg" | | |
Kyle (kyle)
Member Username: kyle
Post Number: 6 Registered: 08-2006
| Posted on Friday, September 22, 2006 - 02:21 pm: | |
Hi Jamy: I tried to upload a script a few days ago and it failed to upload just like yours. There might be something wrong with the server for this user's forum. Hopefully it will get corrected soon... Regards, Kyle |
Robert G. Nelson (rgnelson)
Moderator Username: rgnelson
Post Number: 687 Registered: 06-2003
| Posted on Friday, September 22, 2006 - 03:43 pm: | |
I set NGRID=1024 in this problem to try to minimize any numerical artifacts. I then ran it to t=3e-11 on both version 5.0.11 and 5.0.12 with no sign of trouble. By this time the signal was down to 1e-9 of the initial value. So I don't know what trouble you are reporting.
|
Robert G. Nelson (rgnelson)
Moderator Username: rgnelson
Post Number: 688 Registered: 06-2003
| Posted on Friday, September 22, 2006 - 03:58 pm: | |
I suspect that the failure in uploads has to do with slow server response and timeouts on the uploading browser. It might work to set your timeouts longer. We are looking into changing our BB system, so hopefully we will get this issue resolved.
|