BOX MODEL - ROADMAP
Objective:
• To run Box Model and Reproduce results for hydrogen-combustion contrail formation
Base studies:
1. Bier et al. (2024)
2. Zink et al. (2026)
Assumption:
1. No soot particles
2. No sulfur particles
3. Ice crystals form on entrained ambient aerosol particles
Framework
1. Set Fixed baseline parameters (Zink Table 1)
EIv, Q, η, U∞, TE, AE, Ujet,E, σaer
2. Define systematic parameter space (Zink Table 2)
Ta, pa, RHice,a, naer, rd,aer, κaer
Generation of all possible parameter combinations using parameter space at prescribed levels
7 × 4 × 4 × 5 × 9 × 4 = 20,160 simulations
Note:
For each simulation case, use 1000 merged FLUDILES trajectories because each trajectory represents one part of
the plume. (For trajectory dataset, I already send email to corresponding authors).
3. For each FLUDILES trajectory
• Calculate trajectory-dependent dilution Dk(t) [Bier et al. (2024), Eq. (2)]
(To Determine aerosol entrainment from the dilution state)
• Calculate effective plume area [Bier et al. (2024), Eq. (9)]
• Thermodynamic evolution [Bier et al. (2022), Eq. (9) & Eq. (10)]
(Plume temp and water vapor equation)
4. Calculate mixing-line slope [Zink et al. (2026), Eq. (2)]
• To determine plume supersaturation
Note:
Represent aerosols using SIPs because simulation particles represent many physical particles
5. Apply κ-Köhler activation theory [Bier et al. (2024), Eq. (14)]
(Aerosols activate into liquid droplets if supersaturation is sufficient)
6. Calculate condensational droplet growth [Bier et al. (2022), Eq. (2)]
7. Apply homogeneous freezing parameterization [Bier et al. (2024), Appendix B]
For each timestep:
• Calculate homogeneous freezing temperature Tfrz using droplet water volume, cooling rate and water activity.
• If plume/droplet temperature reaches or falls below Tfrz,
Then Liquid droplet SIP is converted into ice crystal SIP else it continues condensational growth
Single Trajectory Output: Nice,f,k
Note:
We have to repeat this for all 1000 trajectories
Then Combine trajectory results
All Extensive quantities that include Nice,f,k are summed
Mass-weighted averaged for all Intensive quantities
Simulation Output:
Nice,f
Again, repeat same procedure for all parameter combinations to construct systematic simulation database and
evaluate dependence of N_ice,f on all six variable parameters.