Homec4science
Diffusion SpecMiCP / ReactMiCP (reactmicp/unsaturated) Browse

specmicp/README.md
fb013dc5d277reactmicp/unsaturated

Last Change
Raw File

README.md
View Options

SpecMiCP / ReactMiCP
====================
Overview
--------
SpecMiCP is a speciation solver to find the equilibrium state of a chemical
system. The system is based on a mixed complementarity formulation of the
equilibrium condition for minerals.
For a mineral with volume fraction Sm, the equilibrium condition is :
Sm >= 0, -log(IAPm/K) >= 0, -Sm*log(IAPm/K) = 0
where IAPm is the ion activity product and K the equilibrium constant.
This condition is reformulated using the penalized Fischer-Burmeister C-function and the system is solved
using a semismooth method.
ReactMiCP is a reactive transport solver built on top of specmicp.
It uses the operator splitting approach to couple transport and chemistry.
Modules
-------
SpecMiCP is not (yet) a program but a set of libraries that can be used to solve
specific problems.
The following modules are already available :
- **specmicp** : *core* of the library, set the system, and solve it
- **reactmicp** : the reactive transport solver
- **database** : manage the database, parse it from a file, swap the basis,
select the correct species, ...
- **micpsolver** : Solve a Mixed Complementarity problem
- **odeint** : integration methods for solving ordinary differential
equations (for problem involving kinetics for example)
- **dfpm** : a finite element module
The **micpsolver** and **odeint** modules can be use independantly.
Examples of use are provided in the examples directory.
Using SpecMiCP
==============
Examples
--------
To build and install SpecMiCP/ReactMiCP see the INSTALL file.
To learn how to use the API the best is to look at the examples.
The starting point should be the leaching in CO2-saturated examples :
examples/reactmicp/saturated_react/carbonation.cpp.
This file is heavily commented.
Example of use of the Python interface can be found in the tests/cython directory.
Use in your own project
-----------------------
A cmake module can be used to find SpecMiCP in your own project. It will be installed at
<prefix>/share/specmicp/cmake/FindSpecMiCP.cmake
The python module can be used directly if it is in your python path.
Documentation
-------------
The API is documented with Doxygen and the documentation can be generated with the command :
make doc
A (maybe outdated) version of the doc can be found [here][3].
Warning
-------
The code is in development and the API is not considered stable yet.
Python binding
--------------
A python binding is available at the following address :
https://bitbucket.org/specmicp/specmicppy
About
=====
SpecMiCP is developped by F. Georget (fabieng aT princeton DoT edu). It is part
of my Ph.D. work. The purpose of the PhD is to develop a reactive transport model
to model the coupling between hydration, drying and carbonation in cement paste.
This Ph.D. is funded by [Princeton university][6], [Lafarge][4] and the [CSTB][5].
References :
- F. Georget, J. H. Prévost, and R. J. Vanderbei. A speciation solver for cement paste modeling and the semismooth Newton method . Cement and Concrete Research, 68(0):139--147, 2015.
- F. Georget, J. H. Prévost and B. Huet Validation of a reactive transport solver based on a semismooth speciation solver (submitted to Computational Geosciences)
A list of the references used in the code can be found in the [documentation][7].
[2]: http://www.empa.ch/plugin/template/empa/*/62204/---/l=1
[3]: http://www.princeton.edu/~fabieng/specmicpdoc/index.html
[4]: http://www.lafarge.com/en
[5]: http://www.cstb.fr/
[6]: http://princeton.edu
[7]: http://www.princeton.edu/~fabieng/specmicpdoc/citelist.html
c4science · Help