0% found this document useful (0 votes)
2 views2 pages

Methodology (Box Model)

The document outlines a roadmap for running the Box Model to reproduce results related to hydrogen-combustion contrail formation, based on studies by Bier et al. (2024) and Zink et al. (2026). It details the framework for conducting 20,160 simulations using fixed baseline parameters and systematic parameter space, including calculations for trajectory-dependent dilution, effective plume area, and mixing-line slope. The final goal is to construct a systematic simulation database to evaluate the dependence of ice crystal formation on various parameters.
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
2 views2 pages

Methodology (Box Model)

The document outlines a roadmap for running the Box Model to reproduce results related to hydrogen-combustion contrail formation, based on studies by Bier et al. (2024) and Zink et al. (2026). It details the framework for conducting 20,160 simulations using fixed baseline parameters and systematic parameter space, including calculations for trajectory-dependent dilution, effective plume area, and mixing-line slope. The final goal is to construct a systematic simulation database to evaluate the dependence of ice crystal formation on various parameters.
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

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.

You might also like