Transport properties computer calculation

Perhaps the most time-consuming part of the detailed transport property computation is the solution of Eqs. (3.25a), which involves the inversion of approximately a 3AT x 3AT matrix. The solution of this set of equations is necessary for the evaluation of the thermal conductivity and the thermal diffusion coefficients in the mixture. However, in hydrocarbon combustion systems N may be 20 or more, and if all these components are included the cost of the calculation becomes prohibitive. One approach to overcoming the problem is to exclude from this part of the calculation minor reaction intermediates which are unlikely to contribute appreciably either to the thermal conductivity or to thermal diffusion. This reduces the mixture to seven or eight effective components which are normally H, N2 or Ar, O2, H2, CO, CO2, H2O, and the initial hydrocarbon. 

Burcat [ Thermochemical Data for Combustion Calculations, in Combustion Chemistry. (W. C. Gardiner, Jr., ed.), Chapter 8. John Wiley Sons, New York, 1984] discusses in detail the various sources of thermochemical data and their adaptation for computer usage. Examples of thermochemical data tit to polynomials for use in computer calculations are reported by McBride, B. J Gordon, S., and Reno, M. A., Coefficients for Calculating Thermodynamic and Transport Properties of Individual Species, NASA, NASA Langley, VA, NASA Technical Memorandum 4513, 1993, and by Kee, R. J., Rupley, F. M and Miller, J. A., The Chemkin Thermodynamic Data Base, Sandia National Laboratories, Livermore, CA, Sandia Technical Report SAND87-8215B, 1987. 

Computer simulation is being used increasingly in diverse areas of science in the past few years. It has also emerged to become one of the powerful means for investigating condensed matter (/). The principal tools employed in computer simulation are the Monte Carlo and the molecular dynamics methods. In these methods, properties of a collection of particles, usually between 30 and 1000 in number, interacting via a potential />(r) are obtained numerically. Reliable estimates of equilibrium and transport properties as well as microscopic properties can be obtained from such calculations. 

In a practical implementation, the domain on which the phase volume functions are specified is typically a cubic grid of Nx x Ny x Nz voxels, which corresponds to real dimensions of Lx — hNx, Ly - hNy, and Lz — hNz, where h is the voxel size. We will further call this region of real space the computational unit cell. The relationship between the unit cell and the multiphase medium of interest depends on the absolute dimensions of the medium and on the spatial resolution at which the medium is represented (feature dimensions). The unit cell can either contain the entire medium and some void space surrounding it, as in the case of virtual granules described in Section IV.D below, or be a sample of a much larger (theoretically infinite) medium, as in the case of transport properties calculation, described in Section II.E below. 

In the previous section we computed thermal transport coefficients for a water cluster whose size is reasonably similar to that of a typical globular protein. The calculation of thermal transport properties of proteins turns out not to be so simple. For one thing, there is considerable computational and experimental evidence to suggest that energy transport in proteins is non-Brownian. 

A computational design procedure of a thermoelectric power device using Functionally Graded Materials (FGM) is presented. A model of thermoelectric materials is presented for transport properties of heavily doped semiconductors, electron and phonon transport coefficients are calculated using band theory. And, a procedure of an elastic thermal stress analysis is presented on a functionally graded thermoelectric device by two-dimensional finite element technique. First, temperature distributions are calculated by two-dimensional non-linear finite element method based on expressions of thermoelectric phenomenon. Next, using temperature distributions, thermal stress distributions are computed by two-dimensional elastic finite element analysis. 

Computational approach to the design of thermoelectric FGM device was presented, which includes theoretical calculation of transport properties of thermoelectric... 

The methods developed in this book can also provide input parameters for calculations using techniques such as mean field theory and mesoscale simulations to predict the morphologies of multiphase materials (Chapter 19), and to calculations based on composite theory to predict the thermoelastic and transport properties of such materials in terms of material properties and phase morphology (Chapter 20). Material properties calculated by the correlations presented in this book can also be used as input parameters in computationally-intensive continuum mechanical simulations (for example, by finite element analysis) for the properties of composite materials and/or of finished parts with diverse sizes, shapes and configurations. The work presented in this book therefore constitutes a "bridge" from the molecular structure and fundamental material properties to the performance of finished parts. 

For a general discussion, see Chemical Reaction Equilibrium Analysis Theory and Algorithms by W. R. Smith and R. W. Missen, John Wiley Sons, New York, 1982. Also see Fortran TV Computer Program for Calculation of Thermodynamic and Transport Properties of Complex Chemical Systems by R.A. Svehla and B. J. McBride, National Aeronautics and Space Administration Technical Note D-7055, January 1973 Rand s Chemical Composition Program, Rand Corporation, Santa Monica, Calif. and others. 

In research there are many fascinating areas for speculation. At the molecular level there will have to be a small band of devoted scholars who will continue to do the painstaking, detailed derivations which only other kinetic theorists can appreciate fully. The past few years have seen a renaissance in the kinetic theory of macromolecular fluids, and many new avenues have been opened or suggested. Because of the development of high-speed computers, there is also considerable interest in developing the molecular dynamics approach to transport properties. It may turn out that transport properties will be predicted more readily by direct calculation of the molecular motions than by doing involved... 

---