代写Air Pollution Engineering (MCHE 4400/6400), Spring 2024 Homework Assignment - 4代写Matlab编程

日期：2024-03-25 03:12

Homework Assignment - 4

Air Pollution Engineering (MCHE 4400/6400), Spring 2024

Due on: Friday, March 29 at 11:59 pm

Submit to eLC (follow instructions that Anita posted oneLC)

Problem statement:

You are tasked with building a box model that simulates some key processes of smog formation over a small city during a subsidence inversion episode for the period between 8:00AM and 6:00PM. The city has dimensions of 10 km by 10 km, and the mixing height can be considered to be 200 m. The wind speed is 1.5 m/s. The temperature is 20 °C and the relative humidity is 60%. Your model should account for the following chemical reactions:

i)          O3   photolysis leading to OH production. Use the same rate constants for the relevant reactions as those given in class (assume that the photolysis rate of O3  does not change overtime).

ii)         Reaction of OH with volatile organic compounds (VOCs) according to: VOCs + OH à Products (k = 2 x 10-12  (cm3/molec.s)

The emission rates and background concentrations (which you can also assume to be initial concentrations in the box at 8:00AM) are:

Pollutant	Emission rate	Background concentration
VOCs	50 kmol/hour	10 ppb
O3	120 kmol/hour*	10 ppb
OH	0	106 molec/cm3
O1D	0	0

*Note: This is a simplification. As you know, O3 is a secondary pollutant that is not emitted directly, but forms in the atmosphere. You can think of the emission rate here as a ‘formation rate.’

1.   Present the governing box model differential equations and chemical reactions. (30 pts)

2.   Implement the resulting system of equations numerically in MATLAB in order to calculate the evolution of the concentrations of each of the pollutants between 8:00AM and 6:00PM. Use a time step of 0.1 seconds. Generate plots of (i) O3 and VOCs mixing ratios versus time on the same plot, and (ii) OH molecular concentration versus time on a separate plot. (70 pts)