PFLOTRAN
From WaterWiki
For details on PFLOTRAN, please visit the project website at LANL
A number of environmental problems of national importance, e.g. geologic sequestration of CO2, require detailed modeling of reactive flows through geologic media. Typical simulations are very computationally demanding, involving 10 or more chemical degrees of freedom on a grid of millions of nodes.
A massively parallel 3-D reservoir simulator, PFLOTRAN, that can model multiphase reactive flows in geologic formations based on continuum scale mass and energy conservation equations has been developed at LANL/ORNL. It employs the PETSc toolkit's modular and efficient Newton-Krylov solver framework.
For details, check out the SciDAC Groundwater project involving PFLOTRAN.
Contents |
[edit] Performance Analysis and Optimization on Cray XT4 (Jaguar) at ORNL
There are a lot of additional results obtained exploring several environmental optimizations, I/O options including HDF options on Cray XT4 (Jaguar at ORNL). If interested, please email me.
[edit] Performance Analysis on Neptune (Opteron Cluster)
[edit] Pflow
The performance of Pflow (example problem - tough) has been analyzed from 1 to 16 processors.
- Baseline Performance Results and Profiling Overhead
- TAU Profiling Outline
- TAU Profiling Function Distribution
- TAU Profiling User events
Problem Description (CO2 Injection):
This problem involves injection of supercritical CO2 into a reservoir at depth.
[edit] Ptran
The performance of Ptran (example problem - copper_leaching) has been analyzed from 1 to 3 processors.
- Baseline Performance Results and Profiling overhead
- TAU Profiling Outline
- TAU Profiling Function Distribution
- TAU Profiling User events
Problem Description (Copper_leaching):
The copper leaching problem describes leaching of a copper ore body through injection of sulfuric acid. A quarter symmetry element of a five-spot well pattern is represented with injection of sulfuric acid at the center of the five-spot pattern and extraction from a corner well. A 2D planar geometry is used with 30 x 30 nodes. A system of 12 partial differential equations are solved simultaneously for 12 primary chemical species. Reactions take place between primary and secondary minerals in the ore body (porphyry copper deposit) and an aqueous solution. The primary ore mineral is chrysocolla which is dissolved and replaced with secondary minerals consisting of gypsum, amorphous silica, jurbanite alunite, and others. Alunite precipitates at an extremely sharp front which is difficult to resolve with a fixed grid.
