Summary Built around solution of conservation of energy by finite differences, HEAT is a graphically-interfaced application that allows design of a computational mesh representing rock geometry, properties, and magma body emplacement.
The original version of HEAT was a FORTRAN program described in the book "Volcanology and Geothermal Energy" by K. Wohletz and G. ... Heiken (University of California Press, 1991, 432 pages). Generally following the solution scheme described in Appendix E of that book, this HEAT version includes a number of improvements.
HEAT is a graphically interfaced application written for a PC running Microsoft WindowsTM. Designed to study 2-D and 3-D, transient thermal regimes in and around volcanic and magmatic intrusions, HEAT models a variety of geologic structures and rock properties and their effect on both conductive and convective heat flow. The graphical interface is readily used to develop and tailor the simulation to represent most geological conditions of magma intrusion and geological structure. Calculated thermal regimes are color encoded and updated in graphical display with each update stored as a file for future playback. This modeling goes well beyond that done by analytical solution of 1-D linearized expressions of thermal diffusion in that it calculates the nonlinear effects of heterogeneous media, temporal, spatial, and thermally varying properties. As with most geophysical modeling, the results of these calculations are not mathematically unique for inversion applications. However, the geological constraints applied from the first part of this study will greatly reduce the number of simulations that might satisfactorily fit observations.
This modeling goes well beyond that done by analytical solution of 1-D linearized expressions of thermal diffusion in that it calculates the nonlinear effects of heterogeneous media, temporal, spatial, and thermally varying properties. HEAT employs an explicit finite differencing scheme rather than an alternating direction implicit one in order to insure that the original differential equation solved is exactly reproduced by the finite difference equation when time and spatial steps are infinitesimal. Truncation errors that might evolve when using very short time steps are minimized by utilizing double precision. Continuous thermal gradients are assigned along the boundaries and initial conditions use a designated regional thermal gradient. Latent heats of fusion/crystallization are solved for all rocks including magma where temperatures are in that range. Convective heat transfer in the magma is calculated as a function of temperature and composition; convective heat transfer in fluid saturated rocks is approximated by Nusselt number. New magmas may be introduced into the computational mesh at any time during the simulation. All rock/magma properties are assigned by the user and they include: density, porosity (fluid saturation), heat capacity, initial temperature, spatially and thermally varying thermal conductivities, and location. As mentioned earlier, the code has been applied to several geologic areas to test its suitability. A version of heat has been adapted for laboratory rock melting experiments involving thermal variations of rocks melted by a moving hot molybdenum probe. Results thus far have shown the method predicts measured temperatures with enough accuracy to make engineering designs based on it.
Access Constraints DISCLAIMER: The development of these programs was sponsored by the United States Government. Neither the United States Government, nor the United States Department of Energy, nor the University of California, nor any of their employees makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or ... usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. References herein to any specific commercial product, process, or service by trade name, mark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof.
Use Constraints 32-bit (Windows 9x/ME/NT4/2000/XP) version 4.05.0169, 3.2 MB, 26 November 2002; self-installing file (LA-CC 99-27; Copyright ? 1998-2002 UC). When using this program remember the old adage of "garbage in garbage out!"
If you are updating KWare HEAT, you should uninstall it before reinstalling this newer version.
Old version for very old PCs: 16-bit (Windows 3.1, version 16.37.0024) 2.1 MB, 03/12/98 is an old version only to be used on old operating platforms
Name:
KEN
WOHLETZ
Phone:
(505) 667-9202
Fax:
(505) 665-3287
Email:
wohletz at lanl.gov
Email:
KWare at lanl.gov
Contact Address:
Earth and Environmental Science Division
P.O. Box 1663
Los Alamos Natl. Lab City:
Los Alamos
Province or State:
NM
Postal Code:
87545
Country:
USA
Personnel
TYLER
B.
STEVENS Role:
SERF AUTHOR
Phone:
(301) 614-6898
Fax:
301-614-5268
Email:
Tyler.B.Stevens at nasa.gov
Contact Address:
NASA Goddard Space Flight Center
Global Change Master Directory City:
Greenbelt
Province or State:
MD
Postal Code:
20771
Country:
USA
KEN
WOHLETZ Role:
TECHNICAL CONTACT
Phone:
(505) 667-9202
Fax:
(505) 665-3287
Email:
wohletz at lanl.gov
Email:
KWare at lanl.gov
Contact Address:
Earth and Environmental Science Division
P.O. Box 1663
Los Alamos Natl. Lab City:
Los Alamos
Province or State:
NM
Postal Code:
87545
Country:
USA
Links
