A generic reaction-based biogeochemical simulator

Yilin Fang, Steve B. Yabusaki, Gour Tsyh Yeh

Research output: Contribution to journalArticlepeer-review

Abstract

This paper presents a generic biogeochemical simulator, BIOGEOCHEM. The simulator can read a thermodynamic database based on the EQ3/EQ6 database. It can also read user-specified equilibrium and kinetic reactions (reactions not defined in the format of that in EQ3/EQ6 database) symbolically. BIOGEOCHEM is developed with a general paradigm. It overcomes the requirement in most available reaction-based models that reactions and rate laws be specified in a limited number of canonical forms. The simulator interprets the reactions and rate laws of virtually any type for input to the MAPLE symbolic mathematical software package. MAPLE then generates Fortran code for the analytical Jacobian matrix used in the Newton-Raphson technique, which is copiled and linked into the BIOGEOCHEM executable. With this feature, the users are exempted from recoding the simulator to accept new equilibrium expressions or kinetic rate laws. Two examples are used to demonstrate the new features of the simulator.

Original languageEnglish
Pages (from-to)869-878
Number of pages10
JournalDevelopments in Water Science
Volume55
Issue numberPART 1
DOIs
StatePublished - 2004

Fingerprint

Dive into the research topics of 'A generic reaction-based biogeochemical simulator'. Together they form a unique fingerprint.

Cite this