氢氧自由基化学在空气质量模型创建数值稳定性问题?

Question

大气混合比率的氢氧自由基相对较短的寿命(微秒的顺序)。建模时空气质量或天气,时间步通常远远大于羟基的半衰期。这是如何计算数值稳定?是数值稳定的问题简单地忽略了?

Answer 1

氢氧自由基化学不是空气质量模型中显式地建模,所以数值的稳定性不是问题。相反,哦在pseudo-steady-state举行。例如,看到CMAQ的文档州:

CMAQ的气相,哦是假定为pseudo-steady状态,它不是一个运输的物种。这是因为氢氧自由基反应往往是催化(例如,消费和生产)。在水相化学,哦,是被云吸收水和在一个开放的云模型(即。设计方法目前使用),吸收物种(例如,哦)将通过gas-to-cloud分区补充。然而,由于运营商分裂、水、气相化学不同时解决。占这个和其他不确定性预测哦水相浓度(例如,忽视生产反应(过氧化氢+高压- > 2 OH)在aqchem目前未实现),一个稳态假设的哦,采用常规水化学。

通常,如果你想模型类似的氢氧自由基,可以使用盒模型如麦加,而不是化学传输模式。

Answer 2

事实上,有许多化学运输模型解决激进的化学不稳态近似。方程非常僵硬,但是可以集成与合适的方法,通常是隐含的。讨论了一些方法这篇评论。的动态预处理器在评论中提到,包括一组解决硬系统和在大气化学变得有点受欢迎。一些模型,使用生成的代码(默认或作为一个选项)包括GEOS-Chem和WRF-Chem。我使用了KPP实现CB4化学机制空气质量模型。

这样的方法是可能的因为算子分裂,这意味着流程包含在模型解决了非耦合的操作序列。例如,时间步模型可以包含调用显式方案求解平流,然后Crank-Nicolson垂直扩散,然后一个。对化学动力学解算器。