-
Notifications
You must be signed in to change notification settings - Fork 201
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
Implementation of Bremsstrahlung collisions #5537
base: development
Are you sure you want to change the base?
Conversation
for more information, see https://pre-commit.ci
Source/Particles/Collision/BinaryCollision/Bremsstrahlung/BremsstrahlungFunc.H
Fixed
Show fixed
Hide fixed
Since the number of collision is the number of electrons there are some pairs the are not done when the number of species 2 is greater than the number of electrons.
Removed the use of the temporary sigmaC array. Now use quadratic weighting for the photon energy.
for more information, see https://pre-commit.ci
Uses when the calculated cut off is below the default value.
for more information, see https://pre-commit.ci
This gives results that more closely matches the phirad from Seltzer and Berger
Now the ion momentum is also updated
That is based on the electron energy in the center of mass frame. This avoids sqrt of negative when doing momenumt and energy conservation
collisions.collision_names = bremsstrahlung | ||
bremsstrahlung.type = bremsstrahlung | ||
bremsstrahlung.species = electrons ions | ||
bremsstrahlung.Z = 5 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is this to overwrite the ionization state or, since this is a particle-particle mode, can we also rely on the charge / ionization state of the ion we collide with?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The Z
is only used to indicate the atomic number, which is used to select which cross section data to use from the tables.
This PR adds Bremsstrahlung collisions, producing photon particles. It does pair wise collisions between electron and ion species, produce photons. The algorithm conserves momentum and energy.
It uses the cross section tables from Seltzer and Berger, but now only has the cross sections for Z = 1, 2, 5, and 6. A future PR can add individual cross section data as needed or change this to read in data from the warpx-data repo to allow any Z.
Todo: