【comsol学习2】-电化学模块学习与运行错误之循环伏安法

这里参考的官网案例文件仍是上一篇文章中的 电化学模块文件，当然也可以直接参考官网给的 案例文件，在实际模拟过程中，发现模拟结果与官网的实在相差太大，在此记录分析。

一维下的宏观电极上的循环伏安法

作用：记录工作电极上的电流随加载的电极电位变化
与橙子电池不同，循环伏安法只需要研究工作电极上的电流变化，那么从工作电极上发生的电极反应以及影响电极反应的方向分析，可提取两方面信息：电化学反应和传递属性
疑问

• 瞬态分析与稳态分析的含义
• 工作电极上发生的反应是氧化还是还原（这个应该是由施加电位决定的，若电极与电源负极相接，电极表面发生还原反应，反之则是发生氧化反应）
• ‌氧化，还原反应对应的BV方程是什么样的
• ‌法拉第电解定律
• ‌求解峰值电流公式：RS方程
• 双电层电容什么意思

模型定义

几何模型：由于只需要研究工作电极及其附近（半径L），因此可简化为研究工作电极上的一个点到其直线延伸距离(L,长度远大于扩散层能达到的平均厚度)

边界条件：x=0时，指电极表面上A发生氧化反应；x=L时，指反应物A本体浓度，此时生成物B浓度为零（这个应该是在工作电极与电源正极相接时的情况）
电极反应（以还原方向为正方向）：

模型建立过程

新建模型

选择物理场：电分析
变量：cA,cB

模型定义

1.导入模型参数

2.构建几何模型

物理场定义（即电分析定义）

依次定义：
电解质（扩散系数），浓度（边界条件浓度），电极表面（选择所在边界），循环伏安法设置（扫描电位范围），电极反应（化学计量数），双电层电容，初始浓度值（将两物质的浓度与电解质本体浓度及x联系起来），实验温度
注意：在定义物理场时，要特别注意边界选择，正确的选择方式是，例如，电极反应的边界应是边界1，一般我们会直接双击边界1，但实际上这并没有选择成功，应该要在边界选择那里，选择手动，然后在右边的几何图形上选择边界2（实际上，一开始我也是直接选择边界1，但运行结果不对，也试过在手动模式下，在几何图形中选择边界1，但运行结果中的物质浓度变化跟案例结果不同。因此才认为正确的操作应是手动模式下选择其相反的边界，但正确性仍需继续考究）

模型求解与结果

设置参数化扫描（指定扫描速率范围），然后进行运算。
整个模型运行下来就会逐渐清楚该模型的任务是：研究在不同扫描速率以及考虑物质扩散的情况下，工作电极上产生的电流变化。在模型中，重点在于电分析物理场中的边界选择以及边界条件的正确选择。下面展示此重点的重要性：

案例结果	实际模拟结果1	差异原因
		在定义电分析物理场时，错误选择边界，或者说没有选择边界
		同上
		同上

前面有提到边界选择的正确操作，在后面，我重新模拟了两次，比较两种不同边界选择方式产生的结果。以选择电极表面边界为例（边界1）：
方式一：在边界选择时，选择手动模式，然后在几何图形中选择边界1，此时设置及几何窗口变为：

设置窗口	图形窗口

方式二：在边界选择时，选择手动模式，然后在几何图形中选择边界2，此时设置及几何窗口变为：

设置窗口	图形窗口

下面来比较这两种边界选择方式的运行结果：

方式一	方式二	差异
		图形差异看不出来
		电流密度差异也看不出来
		在不同扫描速率下，明显A物质浓度变化不同，另外B物质变化也不同
		这里提取其中一个时间再次展示两种方式下A物质浓度变化差异

这两种方式产生的结果，原因在于边界选择刚好相反，在定义电分析物理场中，需要指定边界的有：

• 浓度（边界x=L即溶液本体处(边界2)，cA=c_bulk，cB=0）
• 电极表面(边界x=0即边界1)
  若两者定义的边界刚好相反，即
  
  这对实际电化学过程中产生的电流变化没有影响，而由于边界厚度L是基于物质的扩散来定义大小的，那么当将物质浓度与边界长度的关系转为图形时，当然是不一样的，就好比，正确的是物质A浓度与边界长度x的函数关系图，由于边界选择错误，实际上展示的是物质A浓度与L-x的函数关系图。