Write a MATLAB code to calculate the ground-level SO2 concentration for 10 km downwind of the stack and 2 km on either side of the plume centerline (a 40 km2 area), for stability classes A, C and F. Generate a contour plot of ground-level concentration for each scenario.

https://drive.google.com/file/d/1FDD1itMYzBu1vrdMfNJZ7-w1XoMI0uIl/view?usp=sharing

Please download the code from the link above.

Overview

In this project, you will use the Point Source Gaussian Plume Model to approximate concentrations of a pollutant downwind of a power plant. The Point Source Gaussian Plume Model is given by

where x is distance along the plume centerline (m), y is the horizontal distance from the plume centerline (m), z is the vertical distance from ground surface at the emission site (m), Q is the source emission rate (µg/s), H is the effective stack height (m), u_H is the average wind velocity at z=H, σ_yand σ_z are the dispersion coefficients in the y and z directions (m).

Model Inputs

The average emission rate of SO₂ from the power plant is 500 g/s, from a stack of height 100 m (assume H=120 m). The topography in the region is relatively flat so we can assume a ground surface at z=0 everywhere. The ground-level windspeed is 2m/s for stability class A, C, and F.

Model Outputs

1. Write a MATLAB code to calculate the ground-level SO₂ concentration for 10 km downwind of the stack and 2 km on either side of the plume centerline (a 40 km² area), for stability classes A, C and F. Generate a contour plot of ground-level concentration for each scenario.