homogeneous nucleation of liquid droplets & expansion chamber example
Copy/pasting the notes @slayoo added from our conversation with @AgnieszkaZaba this morning.
- [x] the laboratory paper uses a Normal distribution rather than Lognormal - perhaps worth adding to PySDM
- [x] a new
ExpansionChamberenvironment (featuring perfect gas expansion, but not featuring hydrostatics) instead of usingParcelwould help in making the code readable - [x] new
adiabatic_exponentformulae - [x] a
SuperParticleSpawningDynamicbase-class could be best to derive from two dynamics:SeedingandHomogeneousLiquidNucleation - [x] rename
SeedingMethodsto sharedSpawningMethods - [ ] could we construct a simple adaptive substepping logic by comparing the nucleation rate and the timestep?
- [x] test idea: J_homogeneous table from S&P
- [ ] test idea: long vs. short timestep
- [ ] test idea: S_max for short vs. long timestep
- [ ] test idea: conservation of total water
- [ ] test idea: check droplet number matches expectation from multiplicity of nucleated SDs
The expansion chamber example is aimed to mimic the laboratory experiments described here: https://arxiv.org/abs/2501.01467
@AgnieszkaZaba @slayoo - Any status update on the super class NewParticle and code restructure to inherit Seeding dynamic from there? Probably with that structure in place I can easily extend to make the HomogeneousLiquidNucleation dynamic. Let me know if you want to call to discuss further!
Codecov Report
Attention: Patch coverage is 88.30409% with 20 lines in your changes missing coverage. Please review.
Project coverage is 85.24%. Comparing base (
17cf24a) to head (7150f84).
Additional details and impacted files
@@ Coverage Diff @@
## main #1492 +/- ##
========================================
Coverage 85.24% 85.24%
========================================
Files 379 389 +10
Lines 9297 9436 +139
========================================
+ Hits 7925 8044 +119
- Misses 1372 1392 +20
:umbrella: View full report in Codecov by Sentry.
:loudspeaker: Have feedback on the report? Share it here.