Conduction, heat numerical integration

Lotkin (L10) gives a scheme for numerical integration of the heat conduction equation in a finite ablating slab, using unequal subdivisions in both space and time variables. Near the melting surface it is advantageous to choose rather small integration steps. Stability characteristics of the method are established. 

Unlike the radiant loss from an optically thin flame, conductive or convective losses never can be consistent exactly with the plane-flame assumption that has been employed in our development. Loss analyses must consider non-one-dimensional heat transfer and should also take flame shapes into account if high accuracy is to be achieved. This is difficult to accomplish by methods other than numerical integration of partial differential equations. Therefore, extinction formulas that in principle can be used with an accuracy as great as that of equation (21) for radiant loss are unavailable for convective or conductive loss. The most convenient approach in accounting for convective or conductive losses appears to be to employ equation (24) with L(7 ) estimated from an approximate analysis. The accuracy of the extinction prediction then depends mainly on the accuracy of the heat-loss estimate. Rough heat-loss estimates are readily obtained from overall balances. 

For creeping flow (0 < Re < 1), the solutions of the conduction-convection equation with flow held of the Hadamard-Rybczynski or Stokes are given by numerical integration [1]. The numerical results show that the concentration contours are not symmetrical (Figure 5.1 and Figure 5.2) and that the how inside and outside the sphere largely inhuences heat or mass transfer. In the case of a sphere with weak viscosity ratio, the heat or mass transfer is facilitated. 

A Tian-Calvet heat flux calorimeter was used in the measurements described in ref. [2]. This type of calorimeter is also called isothermal [4, 5], in contrast to other kinds of calorimeter. A tutorial [6] on heat-conduction calorimetry gives a good account of the technique. Peak integration of the heat flux against time may be performed by a numerical integration method, such as Simpson s method, on a personal computer interfaced to the calorimeter [7]. 

For a given inlet gas composition, pressure, and temperature, numerical integration was carried out throughout the reactor. The rate constant k2 in the Temkin-Pyzhev equation was adjusted so that the calculated temperature profiles matched the measured ones. Radial gradients and axial heat conduction were ignored. The catalyst particles were assumed to have the same temperature... 

Heat conductivity has been studied by placing the end particles in contact with two thermal reservoirs at different temperatures (see (Casati et al, 2005) for details)and then integrating the equations of motion. Numerical results (Casati et al, 2005) demonstrated that, in the small uj regime, the heat conductivity is system size dependent, while at large uj, when the system becomes almost fully chaotic, the heat conductivity becomes independent of the system size (if the size is large enough). This means that Fourier law is obeyed in the chaotic regime. 

The first type, which includes, for example, the problem of strong explosion or propagation of heat in a medium with nonlinear thermal conductivity [3], is characterized by the fact that the exponents are found from physical considerations, from the conservation laws and their dimensionality. In addition, the exponents turn out to be rational numbers. The task of the calculation is to find the dimensionless functions by integration of ordinary differential equations. After this the problem is completely solved, since the numerical constants are determined by normalizing the solution to the conserved quantity (the total energy released in these examples). 

Example 5.3 The Semi-infinite Solid with Variable Thermophysical Properties and a Step Change in Surface Temperature Approximate Analytical Solution We have stated before that the thermophysical properties (k, p, Cp) of polymers are generally temperature dependent. Hence, the governing differential equation (Eq. 5.3-1) is nonlinear. Unfortunately, few analytical solutions for nonlinear heat conduction exist (5) therefore, numerical solutions (finite difference and finite element) are frequently applied. There are, however, a number of useful approximate analytical methods available, including the integral method reported by Goodman (6). We present the results of Goodman s approximate treatment for the problem posed in Example 5.2, for comparison purposes. 

Temperature profiles can be determined from the transient heat conduction equation or, in integral models, by assuming some functional form of the temperature profile a priori. With the former, numerical solution of partial differential equations is required. With the latter, the problem is reduced to a set of coupled ordinary differential equations, but numerical solution is still required. The following equations embody a simple heat transfer limited pyrolysis model for a noncharring polymer that is opaque to thermal radiation and has a density that does not depend on temperature. For simplicity, surface regression (which gives rise to convective terms) is not explicitly included. 

Conduction is treated from both the analytical and the numerical viewpoint, so that the reader is afforded the insight which is gained from analytical solutions as well as the important tools of numerical analysis which must often be used in practice. A similar procedure is followed in the presentation of convection heat transfer. An integral analysis of both free- and forced-convection boundary layers is used to present a physical picture of the convection process. From this physical description inferences may be drawn which naturally lead to the presentation of empirical and practical relations for calculating convection heat-transfer coefficients. Because it provides an easier instruction vehicle than other methods, the radiation-network method is used extensively in the introduction of analysis of radiation systems, while a more generalized formulation is given later. 

Problems of physical limits for miniaturization of electronic devices in integrated circuits are discussed. The quantization of both electrical and thermal conductance in nanostructures is considered and estimated numerically. Problems of heat exchange in nanostructures are also discussed. 

The relevance of contact resistances in SPS process has been simulated and confirmed, with the simulation shown in Fig. 6.15 as an example [4]. The system simulated is a Model 1050-Sumitomo SPS, where a solid graphitic cylinder is inserted into the die. The 2D cylindrical coordinate system of coupled thermal and electrical problems is numerically solved by using Abaqus (FEM). The heat losses due to radiation from all exposed surfaces, except those on the ends of the rams, have been considered, where a constant temperature of 25 °C is used for the simulation. Thermophysical parameters of all materials are available in that study. A proportional feedback controller based on the outer surface temperature of the die is modeled, in order to determine the voltage drop applied at two ends of the rams. This controller is used to imitate the actual proportional integral derivative (PID), which is observed in real SPS facilities. It is used to apply electric power input to the system when experiments are conducted in terms of temperature controlling. 

Equation (6.109) is known as Fourier s law of heat conduction. As noted above, there will be a potential contribution as well. Solution to the integral equation, Eq. (6.95), for the function A w) gives numerical values of the thermal conductivity from Eq. (6.110)." ... 

On the basis of the theory of numerical methods and mathematical modeling the problem of the calculation and forecast of the distribution of the temperature field in a two-phase nanocomposite environment is solved. The mathematical statement of the problem is formulated as the integral equation of thermal balance with a heat flux taken into account, which changes according to Fourier s law. Jumps of enthalpy and heat conductivity coefficient are considered. Various numerical schemes and methods are examined and the best one is selected - the method of control volume. Calculation of the dynamics of the temperature field in the nanostructure is hold using the software. 

---