Mathematical modeling of species transformations in aquatic sediments is usually based on numerical solutions to the same general one-dimensional mass-conservation equation and is likely to require substantial computation time. In this paper we present a fast numerical solution to this equation. The solution is suited for both single and multi-component models and it is based on an implicit control volume discretization of the general mass-conservation equation. The solution consists of two algorithms, one that decomposes the discretization matrix once and one that subsequently produces multiple solutions with minimal computational effort. A unique feature of these algorithms is that values of boundary conditions can vary as a simulation progresses without requiring new decompositions of the discretization matrix. This feature can reduce computation time significantly relative to commonly used procedures for modeling dynamic systems. Finally, we present four examples in which the numerical solution is applied to specific problems. From these examples guidelines are derived for the discretization in space and time required to obtain precise solutions of the general mass-conservation equation.
The Journal of Marine Research publishes peer-reviewed research articles covering a broad array of topics in physical, biological and chemical oceanography. Articles that deal with processes, as well as those that report significant observations, are welcome. In the area of biology, studies involving coupling between ecological and physical processes are preferred over those that report systematics. Authors benefit from thorough reviews of their manuscripts, where an attempt is made to maximize clarity. The time between submission and publication is kept to a minimum; there is no page charge.