Heat boundary conditions

Reactions in porous catalyst pellets are Invariably accompanied by thermal effects associated with the heat of reaction. Particularly In the case of exothermic reactions these may have a marked influence on the solutions, and hence on the effectiveness factor, leading to effectiveness factors greater than unity and, In certain circumstances, multiple steady state solutions with given boundary conditions [78]. These phenomena have attracted a great deal of interest and attention in recent years, and an excellent account of our present state of knowledge has been given by Arls [45]. 

In the finite element solution of the energy equation it is sometimes necessary to impose heat transfer across a section of the domain wall as a boundary condition in the process model. This type of convection (Robins) boundary condition is given as... 

BC (MAXDF) BOUNDARY CONDITIONS ARRAY VHEAT GENERATED VISCOS HEAT... 

Force field calculations often truncate the non bonded potential energy of a molecular system at some finite distance. Truncation (nonbonded cutoff) saves computing resources. Also, periodic boxes and boundary conditions require it. However, this approximation is too crude for some calculations. For example, a molecular dynamic simulation with an abruptly truncated potential produces anomalous and nonphysical behavior. One symptom is that the solute (for example, a protein) cools and the solvent (water) heats rapidly. The temperatures of system components then slowly converge until the system appears to be in equilibrium, but it is not. 

NUh2 is the Nusselt number for uniform heat flux boundary condition along the flow direction and periphery. 

When q is zero, Eq. (5-18) reduces to the famihar Laplace equation. The analytical solution of Eq. (10-18) as well as of Laplaces equation is possible for only a few boundary conditions and geometric shapes. Carslaw and Jaeger Conduction of Heat in Solids, Clarendon Press, Oxford, 1959) have presented a large number of analytical solutions of differential equations apphcable to heat-conduction problems. Generally, graphical or numerical finite-difference methods are most frequently used. Other numerical and relaxation methods may be found in the general references in the Introduction. The methods may also be extended to three-dimensional problems. 

In some convection equations, such as for turbulent pipe flow, a special correction factor is used. This factor relates to the heat transfer conditions at the flow inlet, where the flow has not reached its final velocity distribution and the boundary layer is not fully developed. In this region the heat transfer rate is better than at the region of fully developed flow. 

This is not as simple as it sounds, as the boundary conditions are very often not well known. The number and distribution of heat sources and ventilation parameters, particularly in naturally ventilated surroundings, are very often not known or even vary. The CFD engineer has to cope with this situation and can do one of the following ... 

The temperature boundary condition is in many problems of buoyant flow, e.g., heat sources (machines) or heat sinks (cold glazings), is of great importance. 

Under steady-state conditions, the temperature distribution in the wall is only spatial and not time dependent. This is the case, e.g., if the boundary conditions on both sides of the wall are kept constant over a longer time period. The time to achieve such a steady-state condition is dependent on the thickness, conductivity, and specific heat of the material. If this time is much shorter than the change in time of the boundary conditions on the wall surface, then this is termed a quasi-steady-state condition. On the contrary, if this time is longer, the temperature distribution and the heat fluxes in the wall are not constant in time, and therefore the dynamic heat transfer must be analyzed (Fig. 11.32). 

To be specific, we consider the two-dimensional growth of a pure substance from its undercooled melt in about its simplest form, where the growth is controlled by the diffusion of the latent heat of freezing. It obeys the diffusion equation and appropriate boundary conditions [95]... 

Al-Mg (5000 Series) and Al-Mg-Si (6000 Series) In the binary alloy system strength is obtained mainly by strain hardening. Stress corrosion is thought to be associated with a continuous grain boundary film of Mg,Alg which is anodic to the matrix . Air cooling prevents the immediate formation of such precipitates, but they form slowly at ambient temperatures. Thus only low Mg alloys are non-susceptible (Al-3% Mg). Widespread precipitation arising from plastic deformation with carefully controlled heat-treatment conditions can lower susceptibility. Al-5Mg alloys of relatively low susceptibility are subjected to such treatments. Mn and Cr... 

In most cases, however, heat transfer and mass transfer occur simultaneously, and the coupled equation (230) thus takes into account the most general case of the coupling effects between the various fluxes involved. To solve Eq (230) with the appropriate initial and boundary conditions one can decouple the equation by making the transformation (G3)... 

Equation 10.66 is referred to as Fick s Second Law. This also applies when up is small, corresponding to conditions where C, is always low. This equation can be solved for a number of important boundary conditions, and it should be compared with the corresponding equation for unsteady state heat transfer (equation 9.29). 

In general, the axial heat conduction in the channel wall, for conventional size channels, can be neglected because the wall is usually very thin compared to the diameter. Shah and London (1978) found that the Nusselt number for developed laminar flow in a circular tube fell between 4.36 and 3.66, corresponding to values for constant heat flux and constant temperature boundary conditions, respectively. 

The micro-channels utilized in engineering systems are frequently connected with inlet and outlet manifolds. In this case the thermal boundary condition at the inlet and outlet of the tube is not adiabatic. Heat transfer in a micro-tube under these conditions was studied by Hetsroni et al. (2004). They measured heat transfer to water flowing in a pipe of inner diameter 1.07 mm, outer diameter 1.5 mm, and 0.600 m in length, as shown in Fig. 4.2b. The pipe was divided into two sections. The development section of Lj = 0.245 m was used to obtain fully developed flow and thermal fields. The test section proper, of heating length Lh = 0.335 m, was used for collecting the experimental data. 

The analysis of the behavior of the fluid temperature and the Nusselt number performed for a circular tube at the thermal wall boundary condition 7(v = const, also reflects general features of heat transfer in micro-channels of other geometries. 

One particular characteristic of conduction heat transfer in micro-channel heat sinks is the strong three-dimensional character of the phenomenon. The smaller the hydraulic diameter, the more important the coupling between wall and bulk fluid temperatures, because the heat transfer coefficient becomes high. Even though the thermal wall boundary conditions at the inlet and outlet of the solid wall are adiabatic, for small Reynolds numbers the heat flux can become strongly non-uniform most of the flux is transferred to the fluid at the entrance of the micro-channel. Maranzana et al. (2004) analyzed this type of problem and proposed the model of channel flow heat transfer between parallel plates. The geometry shown in Fig. 4.15 corresponds to a flow between parallel plates, the uniform heat flux is imposed on the upper face of block 1 the lower face of block 0 and the side faces of both blocks... 

A variety of studies can be found in the literature for the solution of the convection heat transfer problem in micro-channels. Some of the analytical methods are very powerful, computationally very fast, and provide highly accurate results. Usually, their application is shown only for those channels and thermal boundary conditions for which solutions already exist, such as circular tube and parallel plates for constant heat flux or constant temperature thermal boundary conditions. The majority of experimental investigations are carried out under other thermal boundary conditions (e.g., experiments in rectangular and trapezoidal channels were conducted with heating only the bottom and/or the top of the channel). These experiments should be compared to solutions obtained for a given channel geometry at the same thermal boundary conditions. Results obtained in devices that are built up from a number of parallel micro-channels should account for heat flux and temperature distribution not only due to heat conduction in the streamwise direction but also conduction across the experimental set-up, and new computational models should be elaborated to compare the measurements with theory. 

The first approach developed by Hsu (1962) is widely used to determine ONE in conventional size channels and in micro-channels (Sato and Matsumura 1964 Davis and Anderson 1966 Celata et al. 1997 Qu and Mudawar 2002 Ghiaasiaan and Chedester 2002 Li and Cheng 2004 Liu et al. 2005). These models consider the behavior of a single bubble by solving the one-dimensional heat conduction equation with constant wall temperature as a boundary condition. The temperature distribution inside the surrounding liquid is the same as in the undisturbed near-wall flow, and the temperature of the embryo tip corresponds to the saturation temperature in the bubble 7s,b- The vapor temperature in the bubble can be determined from the Young-Laplace equation and the Clausius-Clapeyron equation (assuming a spherical bubble) ... 

H2. It may be noted that the HI and H2 boundary conditions for the symmetrically heated passages with no sharp corners (e.g., circular, flat, and concentric annular channels) are identical they are simply designated as H. 

Thermal boundary condition at constant wall heat flux... 

Recent revisions to the boundary conditions (ice-sheet topography and sea surface temperatures) have added uncertainty to many of the GCM calculations of the past two decades. Moreover, all of these calculations use prescriptions for at least one central component of the climate system, generally oceanic heat transport and/or sea surface temperatures. This limits the predictive benefit of the models. Nonetheless, these models are the only appropriate way to integrate physical models of diverse aspects of the Earth systems into a unified climate prediction tool. 

---