Skip to content

Develop surface reaction #17

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 11 commits into from
May 27, 2023
Merged

Develop surface reaction #17

merged 11 commits into from
May 27, 2023

Conversation

mattldawson
Copy link
Collaborator

Adds a surface reaction type for gas phase species:

X -> Y + Z + ...

where X, Y, and Z are gas-phase species and the reaction rate depends on total particle surface area for a given aerosol phase.

@mattldawson mattldawson requested review from K20shores and jcurtis2 May 13, 2023 16:40
@K20shores K20shores requested a review from jtgasparik May 15, 2023 14:22
end do
end do

allocate(this%condensed_data_int(NUM_INT_PROP_ + 2 + 3 * products%size() + &
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Could you leave a comment here and for the following allocate statement on what the constants are (+2 + 3 *)?

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

done


/** \brief Update reaction data for new environmental conditions
*
* For surface reactions this only involves recalculating the rate
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This doesn't appear to be true for this.

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

updated the comment

double *float_data = rxn_float_data;
double *env_data = model_data->grid_cell_env;

// save the mean free path [m] for calculating condensation rates
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks like this is mean speed and not mean free path.

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

updated the comment

src/util.h Outdated
@@ -150,4 +161,58 @@ static inline double d_gas_aerosol_rxn_rate_constant_d_radius(
mean_free_path__m, radius__m, alpha));
}

/** Calculate the gas-aerosol reaction rate for the continuum regime
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Can we add some sort of citation for the source of this equation?

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

added citation to Tie et al. 2003

@mattldawson mattldawson merged commit 6855309 into main May 27, 2023
Copy link
Collaborator

@jtgasparik jtgasparik left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

none!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants