Jin Zhang,Alexander F. Chadwick,Peter W. Voorhees
Abstract:Modeling microstructure evolution in electrochemical systems is vital for understanding the mechanism of various electrochemical processes. In this work, we propose a general phase field framework that is fully variational and thus guarantees that the energy decreases upon evolution in an isothermal system. The bulk and interface free energies are decoupled using a grand potential formulation to enhance numerical efficiency. The variational definition of the overpotential is used, and the reaction kinetics is incorporated into the evolution equation for the phase field to correctly capture capillary effects and eliminate additional model parameter calibrations. A higher-order kinetic correction is derived to accurately reproduce general reaction models such as the Butler-Volmer, Marcus, and Marcus-Hush-Chidsey models. Electrostatic potentials in the electrode and the electrolyte are considered separately as independent variables, providing additional freedom to capture the interfacial potential jump. To handle realistic materials and processing parameters for practical applications, a driving force extension method is used to enhance the grid size by three orders of magnitude. Finally, we comprehensively verify our phase field model using classical electrochemical theory.
What problem does this paper attempt to address?
This paper attempts to solve the problem of microstructure evolution modeling in electrochemical systems, especially in order to better understand the mechanisms of various electrochemical processes. Specifically, the author proposes a general phase - field framework, which is fully variational and ensures that the energy decreases with evolution in isothermal systems. This framework addresses some limitations of existing models in the following ways:
1. **Energy Decrease Guarantee**: The model is constructed by the full variational method to ensure that the energy of the system always decreases during the evolution process, which is thermodynamically consistent.
2. **Improvement of Numerical Efficiency**: The large - potential formula is used to decouple the bulk free energy and the interface free energy, improving the efficiency of numerical calculations.
3. **Accurate Description of Reaction Kinetics**: The overpotential is defined variationally, and the reaction kinetics is incorporated into the evolution equation of the phase - field variable, so as to accurately capture the capillary effect and avoid additional model parameter calibration.
4. **Higher - Order Kinetic Correction**: The higher - order kinetic correction terms are derived to accurately reproduce the kinetic behaviors of general reaction models (such as Butler - Volmer, Marcus, Marcus - Hush - Chidsey models).
5. **Treatment of Potential Jumps**: The electrostatic potentials in the electrode and the electrolyte are considered as independent variables respectively, providing additional degrees of freedom to capture the interface potential jumps and accurately capture the potential jumps varying along the interface.
6. **Mesh Size Expansion**: The driving - force expansion method is used to increase the mesh size by three orders of magnitude, so as to handle actual materials and processing parameters, which is suitable for micron - scale simulations.
7. **Model Verification**: The proposed phase - field model is comprehensively verified, including equilibrium behavior, the influence of interface energy or capillarity, the recovery ability of general reaction kinetics, the coupling of ion diffusion and migration, and time - dependent behavior.
In summary, this paper aims to solve the problem of microstructure evolution modeling in electrochemical systems by proposing a general and efficient phase - field framework, especially providing theoretical support for performance optimization and material design in practical applications such as batteries.