-
Notifications
You must be signed in to change notification settings - Fork 218
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Laser back scattering #3511
Comments
@cbontoiu Thanks for posting this question. It might be (unintuitively) caused by the details of your laser setup. Could you please state the following details of your setup:
According to our absorber rules, the positive y-axis absorber without moving window is only active if @psychocoderHPC I remember you explained to me, why we had this rule years ago. But I do not remember. Could you state the reason why we have this rule? |
@PrometheusPi I believe what you wrote applies for the min y boundary, while the reflection observed is at max y boundary? There is a set of .param files attached. However I can't figure which .cfg file you used @cbontoiu ? |
I think i understand The logic of that condition you mentioned, it is to properly generate the laser with initial plane y at 0. However again I don't think it's what causes the problem here |
@sbastrakov to run my simulation I use
with |
Hello @PrometheusPi and @sbastrakov and many thanks for looking into this problem. It seems that the laser pulse reaches the positive y absorber around t = 40 fs, and that corresponds to 9620 iterations. I don't know how to obtain |
@sbastrakov Sorry - you are right. I got confused with the old naming. " |
So that |
@cbontoiu sorry, it is all somewhat confusing. It is hopefully simplified later this year. So with moving window enabled. Our assumption is that the whole laser (or the "only important" part) fits the global domain in y (not counting the extra 1 GPU layer in y). So we operate as follows: the laser is first generated at near min y boundary according to With the moving window disabled, we do not do that manipulation. Since your attached |
With respect to my last point, just to be sure could you attach file |
Hello @sbastrakov he re is the file Thank you |
Thanks. So moving window is indeed off. So far I am out of ideas, will run a small simulation tomorrow with such a 1x2x1 devices and laser to investigate further |
Hello, I observed this problem again on an array of tubes as shown in the gif attached. The unusual settings from the gas-plasma PIC simulations are:
The cfg file is called Cristian |
@cbontoiu Yes, light is always scattered. IN the underdense case, only a small fraction is scattered in backward direction. In he overdense case, only backward scattering occurs. I am not sure if setting the move point to a value >1 is valid. @psychocoderHPC could you comment on that. What densities do you reach during the simulation (due to self compression of the tubes) and what longitudinal momentum do the electrons reach? |
Setting slide point to >1 will destroy the physics because the longitudinal field will hit the simulation bonderie. The default 0.9 is choosen to avoid that numerical artifacts will touch the simulation box.
Am 8. März 2021 15:35:35 MEZ schrieb Richard Pausch <notifications@github.com>:
…
@cbontoiu Yes, light is always scattered. IN the underdense case, only
a small fraction is scattered in backward direction. In he overdense
case, only backward scattering occurs.
I am not sure if setting the move point to a value >1 is valid.
@psychocoderHPC could you comment on that.
That could mean that your sliding starts after (parts of) the laser
hits the back of the simulation box.
What densities do you reach during the simulation (due to self
compression of the tubes) and what longitudinal momentum do the
electrons reach?
--
You are receiving this because you were mentioned.
Reply to this email directly or view it on GitHub:
#3511 (comment)
|
Actually, I don't know how to extract the longitudinal/transverse momentum. All I can do so far is to extract the energy and I assume it is kinetic energy expressed in J. |
Thank you. This is really important to know. I pushed the moving point beyond one just be able to to see more of the pulse tail while still keeping the simulation domain small and get the 0.25 nm mesh resolution. It all stems from the fact that I have only two GPUs on my machine. |
@cbontoiu To extract the momentum, you can use the phase space plugin or the openPMD based output. For density you can use the openPMD based output. It contains (by default) the density of each species. The openPMD-viewer allows to visualize momentum distribution, phase space and density plots. |
Hello,
I am struggling to find the "perfect" absorber power and thickness at the positive y end. I use 32 cells so for the laser
initPlaneY = 33u
and 1e-03 for the absorber power but I get a massive laser wave bouncing back. Is this because of my plasma density settings? Is it because of the short laser wavelength? am I missing anything here?I attach the model using the typical absorber as defined in
grid.param
but I get a similar effect when using PML. What are my choices?single-thick-tube-3d.zip
Many thanks.
Cristian
The text was updated successfully, but these errors were encountered: