Skip to content
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

avoid that numerical instabilities trigegring the KHI #4462

Draft
wants to merge 1 commit into
base: dev
Choose a base branch
from

Conversation

psychocoderHPC
Copy link
Member

@psychocoderHPC psychocoderHPC commented Feb 16, 2023

related to #3234

We observe since long time that our KHi is seeded by nummerical issues instead a well reproducable condition.
I came up with the idea to use the composit pusher to create maxwell conform fields over 500 steps and then start to real KHI simulation. @PrometheusPi suggested to use a random position for electrons and ions and not using the temperature seeding for electrons.

The results look quite good.
The python tool is adjusted to ignore the first 28 outputs.

./bin/picongpu -d 1 1 1 -g 192 128 128 --periodic 1 1 1 -s 2500 --fields_energy.period 20 --e_png.period 100 --e_png.axis yx --e_png.slicePoint 0.5 --e_png.folder pngElectronsYX --e_macroParticlesCount.period 100 --e_energy.period 20 --e_energy.filter all --i_energy.period 20 --i_energy.filter all

result

  • do not merge this is a work-in-progress PR to provide @PrometheusPi with the setup of my lasts tests

related to ComputationalRadiationPhysics#3234

We observe since long time that our KHi is seeded by nummerical issues
instead a well reproducable condition.
I came up with the idea to use the composit pusher to create maxwell
conform fields over 500 steps and then start to real KHI simulation.
@PrometheusPi suggested to use a random position for electrons and ions
and not using the temperature seeding for electrons.
@psychocoderHPC psychocoderHPC added the component: examples PIConGPU or PMacc examples label Feb 16, 2023
@psychocoderHPC psychocoderHPC marked this pull request as draft February 16, 2023 13:26
@PrometheusPi
Copy link
Member

@psychocoderHPC what is the status of this pull request? Do we need to discuss something? Are these just CI errors and no real code errors?

@psychocoderHPC
Copy link
Member Author

We need to discuss the physics and if my changes make the KHI better or not.
The CI can be ignored it was the day where Simeon deleted our containes^^

@PrometheusPi
Copy link
Member

Since we discussed the options offline a while ago, what is the plan with this PR?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
component: examples PIConGPU or PMacc examples
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants