Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aenh_v1_0.tar.gz(135 Kbytes)
Manuscript Title: A method for the accurate and smooth approximation of standard thermodynamic functions.
Authors: O. Coufal
Program title: SmoothSTF
Catalogue identifier: AENH_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 184(2013)194
Programming language: C++.
Computer: Any computer with gcc version 4.3.2 compiler.
Operating system: Debian GNU Linux 6.0. The program can be run in operating systems in which the gcc compiler can be installed, see http://gcc.gnu.org/install/specific.html.
RAM: 256 Mbytes are sufficient for the table of standard thermodynamic functions with 500 lines
Keywords: Algorithms for functional approximation, Atomic and molecular physics, Thermal plasma, Standard thermodynamic functions.
Classification: 4.9.

Nature of problem:
Standard thermodynamic functions (STF) of individual substances are given by thermal capacity at constant pressure, entropy and enthalpy. STF are continuous and smooth in every temperature interval in which no phase transformations take place. The temperature dependence of STF as expressed by the table of its values is for further application approximated by temperature functions. In the paper, a method is proposed for calculating approximation functions which, in contrast to the hitherto used approximations, are continuous and smooth in every temperature interval.

Solution method:
The approximation functions are determined by coefficients that are calculated by the least square method coupled with meeting the conditions set (calculation of minimum with equality constraints). To calculate the coefficients the values of STF derivatives with respect to temperature must be available in addition to the table of STF values. The values of the derivatives are established using cubic splines and the derivative of the Lagrange interpolation polynomial.

The program was tested and is used on the assumption that the upper bound of the temperature interval in which approximation is being performed does not exceed 50 kK. This value is sufficient in the thermal plasma field. For temperatures higher than 50 kK the calculation might entail difficulties, which could however be removed by calculating with higher than double precision.

Additional comments:
The program package includes a README file and input and output files for a test suite.

Running time:
A few seconds, if the approximation is being determined in the maximum interval, i.e. from ambient temperature to 50 kK.