Proton Radiography Module #895
Conversation
Ignoring particles that will not ever pass through the grid as a speed improvement.
|
Check out this pull request on Review Jupyter notebook visual diffs & provide feedback on notebooks. Powered by ReviewNB |
pheuer
added
the
plasmapy.diagnostics
Codecov Report
@@ Coverage Diff @@
## master #895 +/- ##
==========================================
+ Coverage 96.62% 97.02% +0.40%
==========================================
Files 64 65 +1
Lines 6249 6559 +310
==========================================
+ Hits 6038 6364 +326
+ Misses 211 195 -16
Continue to review full report at Codecov.
|
Co-authored-by: Dominik Stańczak <[email protected]>
Co-authored-by: Dominik Stańczak <[email protected]>
Co-authored-by: Dominik Stańczak <[email protected]>
Co-authored-by: Dominik Stańczak <[email protected]>
Co-authored-by: Dominik Stańczak <[email protected]>
Co-authored-by: Dominik Stańczak <[email protected]>
…/PlasmaPy into proton_radiography_module
Co-authored-by: Dominik Stańczak <[email protected]>
Co-authored-by: Dominik Stańczak <[email protected]>
Co-authored-by: Dominik Stańczak <[email protected]>
Co-authored-by: Dominik Stańczak <[email protected]>
…nthetic_radiograph
|
This feature is large enough that hey, might warrant 0.6 in a few weeks' time (I'm swamped until then)! |
|
@StanczakDominik YAY! clinks virtual champagne Thanks for all your help! Getting this functionality onto the release branch sounds great to me: I have lots of people I'm planning on sharing it with, but it's definitely easier for people to install from release! |
|
Sure, but unless you wanna try handling the release... hope you can wait! I should be able to get it out between March 6 and 14. I'll raise the topic at the telecon tomorrow. |
|
Hahaha no no no sometime in March sounds good to me! |
Proton radiography is an increasingly common diagnostic technique used to interrogate the electric and magnetic fields inside high energy density plasmas. The plasma of interest is positioned between a bright pulsed source of protons and a detector plane. Electric and magnetic fields in the plasma deflect the protons, producing patterns on the detector. Since this represents a non-linear and line-integrated measurement of the fields, the interpretation of these "proton radiographs" is a complicated business. The goal of this module is to provide tools for planning these experiments and analyzing their results.
Eventually, this new module should contain two type different types of functions:
Functions that construct synthetic proton radiographs given simulated or calculated fields, either for planning an experiment or analyzing the results of one.
Functions that use assumptions about symmetries etc. to "invert" experimental radiographs to recover the magnitudes of
the fields that created them.
The SyntheticProtonRadiograph object contained in this initial PR belongs to the first category. The SyntheticProtonRadiograph object is initialized with a
plasmapy.plasma.grid.AbstractGridsubclass containing electric and magnetic fields, as well as the locations of the proton source and detector. A particle tracing algorithm is then run to determine the trajectories of the protons, and the resulting positions in the detector plane are used to create a synthetic radiograph (2D histogram).This function can currently support field data on either regular or non-uniform (although the former is much faster) to support complicated simulation meshes. Time steps are computed in an adaptive manner based on both the fields currently experienced by the particles and the grid resolution. Source and detector positions can be specified anywhere in Cartesian, cylindrical, or spherical coordinates systems to re-create geometric effects.
The example fields generated within the function for testing and demonstration are intentionally simple, but this technique can re-create complicated structures, such as this off-axis radiograph of a cylindrical, radial electric field: