diff --git a/joss/paper.bib b/joss/paper.bib
index bf7e824..b5f93f5 100644
--- a/joss/paper.bib
+++ b/joss/paper.bib
@@ -75,3 +75,72 @@ @article{catmap
volume = {145},
year = {2015},
}
+
+@article{wang2013simulations,
+ author = {Wang, Hainan and Thiele, Alexander and Pilon, Laurent},
+ journal = {The Journal of Physical Chemistry C},
+ number = {36},
+ pages = {18286-18297},
+ title = {Simulations of Cyclic Voltammetry for Electric Double Layers in Asymmetric Electrolytes: A Generalized Modified {P}oisson--{N}ernst--{P}lanck Model},
+ volume = {117},
+ year = {2013},
+ doi={10.1021/jp402181e},
+ }
+}
+
+@article{zheng2023python,
+ title={Python for Electrochemistry: A Free and All-In-One Toolset},
+ author={Zheng, Weiran},
+ journal={ECS Advances},
+ volume={2},
+ number={4},
+ pages={040502},
+ year={2023},
+ publisher={IOP Publishing},
+ doi={10.1149/2754-2734/acff0b}
+}
+
+@article{sulzer2021python,
+ title={Python battery mathematical modelling (PyBaMM)},
+ author={Sulzer, Valentin and Marquis, Scott G and Timms, Robert and Robinson, Martin and Chapman, S Jon},
+ journal={Journal of Open Research Software},
+ volume={9},
+ number={1},
+ year={2021},
+ publisher={Ubiquity Press}
+}
+
+
+@article{kirby2018efficient,
+ author = {Kirby, Robert C. and Mitchell, Lawrence},
+ doi = {10.1137/17M1133208},
+ journal = {SIAM Journal on Scientific Computing},
+ number = {1},
+ pages = {C76-C98},
+ title = {Solver Composition Across the {PDE}/Linear Algebra Barrier},
+ volume = {40},
+ year = {2018},
+}
+
+@article{Mitusch2019,
+ doi = {10.21105/joss.01292},
+ year = {2019},
+ publisher = {The Open Journal},
+ volume = {4},
+ number = {38},
+ pages = {1292},
+ author = {Sebastian K. Mitusch and Simon W. Funke and Jørgen S. Dokken},
+ title = {dolfin-adjoint 2018.1: automated adjoints for {FEniCS} and {F}iredrake},
+ journal = {Journal of Open Source Software} }
+
+ @article{roy2022topology,
+ title={Topology optimization for the design of porous electrodes},
+ author={Roy, Thomas and Salazar de Troya, Miguel A and Worsley, Marcus A and Beck, Victor A},
+ journal={Structural and Multidisciplinary Optimization},
+ volume={65},
+ number={6},
+ pages={171},
+ year={2022},
+ publisher={Springer},
+ doi={10.1007/s00158-022-03249-2}
+}
\ No newline at end of file
diff --git a/joss/paper.md b/joss/paper.md
index e91c6a8..b150770 100644
--- a/joss/paper.md
+++ b/joss/paper.md
@@ -28,45 +28,80 @@ bibliography: paper.bib
The shift from fossil fuels towards renewable energy brings about a substantial increase in clean but intermittent electricity.
Thankfully, diverse electrochemical technologies, including energy storage and electrochemical manufacturing, can harness this surplus energy that would otherwise go to waste.
Managing the growing prevalence of renewable energy underscores the importance of developing and scaling up these technologies.
-Likewise, the electrification of transport creates an increasing need for efficient electrochemical energy storage devices such as batteries and supercapacitors.
-Naturally, simulation tools are required to assist in the design of efficient and industrial-scale electrochemical devices.
+Likewise, the electrification of transport creates an increasing need for energy-dense electrochemical energy storage devices.
+Naturally, simulation tools are required to assist in the design of performant and industrial-scale electrochemical devices.
Modeling and simulation are used extensively to describe the physics of the electrochemical and transport mechanisms in electrochemical devices.
These devices have various applications, such as batteries and supercapacitors, which are used for energy storage.
-Flow batteries have a similar function, but with a flowing electrolyte.
-Electrolyzers can be used to transform carbon dioxide into useful products or create hydrogen from water.
-Proton-exchange membrane fuel cells reverse this process to harness electricity from hydrogen.
+Flow batteries can have a similar function, but with a flowing electrolyte.
+Electrolyzers can be used to transform carbon dioxide into useful products or create hydrogen and oxygen from water.
+On the other hand, proton-exchange membrane fuel cells reverse this process to harness electricity from hydrogen.
While these devices vary wildly in purpose, the governing equations used to describe them are very similar.
-Transport of charged chemical species in a fluid is often modeled using the Nernst-Planck equation,
+The transport of charged chemical species in a fluid is often modeled using the Nernst-Planck equation,
which includes the usual advection and diffusion transport as well as *electromigration*, where charged species are transported by an electric field.
EchemFEM provides a high-level user interface for a finite element implementation of the Nernst-Planck equation.
-The user is simply required to provide physical parameters as well as functions describing for chemical reactions.
-All transport options are then selected using keyword arguments.
+The user is simply required to provide physical parameters as well as functions describing the chemical reactions.
+Then, the desired transport physics are selected using keyword arguments.
Ionic charge can either be modeled using the Poisson equation or the electroneutrality approximation.
-The simulated devices can be resolved electrolyte-electrode interfaces, or have homogenized porous electrodes, in which case electron conduction is also modeled.
-Finite size effects are also available to allow the use of models such as Generalized Modified Poisson-Nernst-Planck (GMPNP).
+The simulated devices can have resolved electrolyte-electrode interfaces or homogenized porous electrodes, in which case electron conduction is also modeled.
+Additionally, finite size effects are available to allow the use of models such as Generalized Modified Poisson-Nernst-Planck (GMPNP) [@wang2013simulations].
EchemFEM is based on Firedrake [@FiredrakeUserManual], an open-source finite element package,
enabling straightforward implementation of the governing equations in Python.
-Firedake has access to scalable, customizable, solvers through its interface with PETSc [@petsc-user-ref,petsc-web-page], allowing for parallelization and scalability on computing clusters.
+Firedake has access to scalable, customizable, solvers through its interface with PETSc [@petsc-user-ref; @petsc-web-page], allowing for parallelization and scalability on computing clusters.
This balance between usability and scalability permits a seamless transition from prototyping to large-scale simulation.
EchemFEM leverages Firedrake's capabilities while further increasing the ease-of-use.
-Indeed, since the governing equations are already implemented, little to no knowledge of Firedrake and the finite element method is required to use EchemFEM
+Indeed, since the governing equations are already implemented, little to no knowledge of Firedrake and the finite element method is required to use EchemFEM.
-
-
-Examples
+Examples: CO2 reduction. flow reactors
# Statement of need
+
+How simulation can help electrochemistry.
+What questions are we trying to answer here?
+Scaling up technologies: scaling up simulations.
+Higher dimensions: many codes are just 1D.
+Fluid flow and Architected electrodes: higher-dimensional effects do matter.
+
+
+At the moment, the most common codes used for electrochemistry simulations are commercial software.
+COMSOL Multiphysics®, with its detailed electrochemistry module, is the most popular, while Simcenter™ STAR-CCM+™ is also used commonly for flowing systems.
+These programs provide simple graphical user interfaces (GUI), which allow users to quickly set up new simulations.
+Additionally, other physics modules such as fluid dynamics are available and can usually be coupled with the electrochemistry simulation.
+However, there are several drawbacks to using such commercial software.
+For instance, license fees can be prohibitively expensive, therefore limiting collaboration.
+Furthermore, the closed nature of the source code limits the flexibility of the software.
+Indeed, it is not possible to implement new discretization schemes and preconditioning approaches that may be required for numerical stability or scalability, respectively.
+Finally, since everything needs to be set up through the GUI, scripting and coupling to other software are difficult tasks.
+
+There is a growing number of open-source software for electrochemistry, especially Python-based packages [@zheng2023python], many of which are specialized for specific applications, notably batteries.
+One such package, PyBaMM [@sulzer2021python], is a battery modelling code with a flexible implementation, allowing for new models and numerical methods to be tested.
+
+
+EchemFEM provides a general framework for simulating electrochemical transport: it is not specific to an application.
+Since it is based on Firedrake, any additional physics that can be implemented in a finite element framework can be coupled to EchemFEM.
+In one of the demos, the Navier-Stokes equations are solved in a reactor with an irregular surface, providing a velocity field for the transport equations.
+
+In some cases, for example for fast flows, stabilization schemes that are not offered in other software may be required.
+For continuous Galerkin (CG) elements, a streamline-upwind Petrov-Galerkin (SUPG) method for the Nernst-Planck equation is provided.
+For discontinuous Galerking (DG), a custom upwind scheme for the Nernst-Planck equation is used [@roy2023scalable].
+In both cases, the upwinding considers the combined advection-migration ``velocity''.
+
+As opposed to commercial software, custom scalable solvers are available in Firedrake.
+A plethora of solver options are available through simple PETSc keywords and custom operators for preconditioning can be defined using Firedrake [@Mitusch2019].
+In @roy2023scalable, scalable block preconditioners were developed for the electroneutral Nernst-Planck equations with DG and implemented in EchemFEM.
-[@roy2023scalable]
+Combining EchemFEM with other Python packages is rather simple.
+In @govindarajan2023coupling, multi-scale simulations for CO2; reduction in flow reactors are performed by coupling a microkinetics model from CatMAP [@catmap] with the GMPNP transport model from EchemFEM.
+The simulations are two orders of magnitude faster than a previous implementation where the transport is done in COMSOL Multiphysics®, due to the tedious interface between the commercial software and CatMAP.
-
+Another desirable feature of Firedrake is its automatic adjoint capabilities, facilitating the straightforward solution of PDE-constrained optimization problems [@Mitusch2019], which have already been used for electrochemical applications [@roy2022topology].
+Optimization problems using EchemFEM are currently being investigated.
# Acknowledgements