Asymptotic analysis of vertical geothermal boreholes in the limit of slowly varyng heat injection rates

Theoretical models for the thermal response of vertical geothermal boreholes often assume that the characteristic time of variation of the heat injection rate is much larger than the characteristic diffusion time across the borehole. In this case, heat transfer inside the borehole and in its immediate surroundings is quasi-steady in the first approximation, while unsteady effects enter only in the far field. Previous studies have exploited this disparity of time scales, incorporating approximate matching conditions to couple the near-borehole region with the outer unsteady temperatura field. In the present work matched asymptotic expansion techniques are used to analyze the heat transfer problem, delivering a rigorous derivation of the true matching condition between the two regions and of the correct definition of the network of thermal resistances that represents the quasi-steady solution near the borehole. Additionally, an apparent temperature due to the unsteady far field is identified that needs to be taken into account by the near-borehole region for the correct computation of the heat injection rate. This temperature differs from the usual mean borehole temperature employed in the literatura.

Difference schemes for numerical solutions of lagging models of heat conduction

Non-Fourier models of heat conduction are increasingly being considered in the modeling of microscale heat transfer in engineering and biomedical heat transfer problems. The dual-phase-lagging model, incorporating time lags in the heat flux and the temperature gradient, and some of its particular cases and approximations, result in heat conduction modeling equations in the form of delayed or hyperbolic partial differential equations. In this work, the application of difference schemes for the numerical solution of lagging models of heat conduction is considered. Numerical schemes for some DPL approximations are developed, characterizing their properties of convergence and stability. Examples of numerical computations are included.

Difference schemes for time-dependent heat conduction models with delay

Different non-Fourier models of heat conduction, that incorporate time lags in the heat flux and/or the temperature gradient, have been increasingly considered in the last years to model microscale heat transfer problems in engineering. Numerical schemes to obtain approximate solutions of constant coefficients lagging models of heat conduction have already been proposed. In this work, an explicit finite difference scheme for a model with coefficients variable in time is developed, and their properties of convergence and stability are studied. Numerical computations showing examples of applications of the scheme are presented.

Retrofit of heat exchanger networks with pressure recovery of process streams at sub-ambient conditions

This paper presents a new mathematical programming model for the retrofit of heat exchanger networks (HENs), wherein the pressure recovery of process streams is conducted to enhance heat integration. Particularly applied to cryogenic processes, HENs retrofit with combined heat and work integration is mainly aimed at reducing the use of expensive cold services. The proposed multi-stage superstructure allows the increment of the existing heat transfer area, as well as the use of new equipment for both heat exchange and pressure manipulation. The pressure recovery of streams is carried out simultaneously with the HEN design, such that the process conditions (streams pressure and temperature) are variables of optimization. The mathematical model is formulated using generalized disjunctive programming (GDP) and is optimized via mixed-integer nonlinear programming (MINLP), through the minimization of the retrofit total annualized cost, considering the turbine and compressor coupling with a helper motor. Three case studies are performed to assess the accuracy of the developed approach, including a real industrial example related to liquefied natural gas (LNG) production. The results show that the pressure recovery of streams is efficient for energy savings and, consequently, for decreasing the HEN retrofit total cost especially in sub-ambient processes.

THE TRANSPORT OF HEAT BY FLOWING GROUNDWATER IN A DISCRETE FRACTURE

In typical theoretical or experimental studies of heat migration in discrete fractures, conduction and thermal dispersion are commonly neglected from the fracture heat transport equation, assuming heat conduction into the matrix is predominant. In this study analytical and numerical models are used to investigate the significance of conduction and thermal dispersion in the plane of the fracture for a point and line sources geometries. The analytical models account for advective, conductive and dispersive heat transport in both the longitudinal and transverse directions in the fracture. The heat transport in the fracture is coupled with a matrix equation in which heat is conducted in the direction perpendicular to the fracture. In the numerical model, the governing heat transport processes are the same as the analytical models; however, the matrix conduction is considered in both longitudinal and transverse directions. Firstly, we demonstrate that longitudinal conduction and dispersion are critical processes that affect heat transport in fractured rock environments, especially for small apertures (eg. 100 μm or less), high flow rate conditions (eg. velocity greater than 50 m/day) and early time (eg. less than 10 days). Secondly, transverse thermal dispersion in the fracture plane is also observed to be an important transport process leading to retardation of the migrating heat front particularly at late time (eg. after 40 days of hot water injection). Solutions which neglect dispersion in the transverse direction underestimate the locations of heat fronts at late time. Finally, this study also suggests that the geometry of the heat sources has significant effects on the heat transport in the system. For example, the effects of dispersion in the fracture are observed to decrease when the width of the heat source expands.

Heat flow in the Black and Mediterranean seas

Results of heat-flow measurements in the western Black Sea and within the African-Sicilian region of the Mediterranean Sea are presented. New data reveal a well-defined correlation with the deep-seated geologic structure of the regions examined. A heat-flow anomaly within the Pantelleria graben, in the authors' opinion, is a consequence of connective heat transfer along channels of disjunctive faults.

Heat flow determinations of ODP Hole 111-504B

Hole 504B in the eastern equatorial Pacific has been the focus of five scientific drilling expeditions since it was first drilled in 1979. During these five legs, a series of temperature logs has been obtained over a time span of almost 8 yr, documenting the geothermal and hydrologic state of the oceanic crust in this region. Immediately following reentry at the onset of ODP Leg 111 operations, a high-resolution temperature probe was lowered into the borehole and a precise record of temperature vs. depth in Hole 504B was recorded down to 1300 mbsf. As was observed during previous legs, the temperature gradient in the upper 400 m was reduced, indicating that downhole flow of cool ocean waters through the casing continued, though at a diminished rate. As subhydrostatic pressures in the upper basement have gradually diminished, the volume of flow has decayed from an estimated 6000-7000 L/hr in late 1979 to about 80 L/hr during Leg 111. At depths below 480 mbsf, a predominantly conductive heat transfer environment enabled the temperature gradient log to be analyzed with respect to lithology on both fine and broad scales. Anomalies in the gradient log in the cased section through the sedimentary column were found to correspond to biostratigraphic age markers and/or sharp changes in sediment composition and texture. Broad variations in temperature gradient within the basement correlated with large-scale porosity trends. Conductive heatflow estimates depict a systematic reduction with depth, ranging from approximately 196 mW/m**2 in the sediments to 120 ± 17 mW/m**2 at 1300 mbsf. Possible causes for this observation were examined from several perspectives, but none was suitably convincing. A fluid instability analysis indicated the likely existence of convection cells within the borehole and substantiated the hypothesis of mixing within the borehole postulated from isotopic and chemical studies of borehole waters. However, such mixing of borehole fluids does not provide an adequate explanation for the heatflow variations, and the disparity between surficial and deep values of heat flow remains unresolved.

Study of plutonium oxide heat generation in a wood moderated container /

The generation of heat due to the radioactive decay of plutonium in the form of plutonium dioxide stored in wooden containers is considered. The results of two experiments in which plutonium dioxide is placed inside wooden blocks and the temperature rise measured at various points is reported. Heat transfer calculations are then performed to show that the data obtained are those which would be expected, i.e., that Fourier's Law adequately describes the situation. The heat build up in a proposed wood moderated shipping container is calculated, using some of the assumptions verified in the described experiments.

Gasdynamic theory of heat

Cover title: Gasdynamic theory of heat transfer.

Patch and whole of surface mass-transfer rate measurements on a stepped rotating cylinder electrode

A recently developed whole of surface electroplating technique was used to obtain mass-transfer rates in the separated flow region of a stepped rotating cylinder electrode. These data are compared with previously reported mass-transfer rates obtained with a patch electrode. It was found that the two methods yield different results, where at lower Reynolds numbers, the mass-transfer rate enhancement was noticeably higher for the whole of the surface electrode than for the patch electrode. The location of the peak mass transfer behind the step, as measured with a patch electrode, was reported to be independent of the Reynolds number in previous studies, whereas the whole of the surface electrode shows a definite Reynolds number dependence. Large eddy simulation results for the recirculating region behind a step are used in this work to show that this difference in behavior is related to the existence of a much thinner fluid layer at the wall for which the velocity is a linear junction of distance from the wall. Consequently, the diffusion layer no longer lies well within a laminar sublayer. It is concluded that the patch electrode responds to the wall shear stress for smooth wall flow as well as for the disturbed flow region behind the step. When the whole of the surface is electro-active, the response is to mass transfer even when this is not a sole function of wall shear stress. The results demonstrate that the choice of the mass-transfer measurement technique in corrosion studies can have a significant effect on the results obtained from empirical data.

Velocity distributions in a plate heat exchanger channel

The development of a Laser Doppler Anemometer technique to measure the velocity distribution in a commercial plate heat exchanger is described. Detailed velocity profiles are presented and a preliminary investigation is reported on flow behaviour through a single cell in the channel matrix. The objective of the study was to extend previous investigations of plate heat exchanger flow patterns in the laminar range with the eventual aim of establishing the effect of flow patterns on heat transfer performance, thus leading to improved plate heat exchanger design and design methods. Accurate point velocities were obtained by Laser Anemometry in a perspex replica of the metal channel. Oil was used as a circulating liquid with a refractive index matched to that of the perspex so that the laser beams were not distorted. Cell-by-cell velocity measurements over a range of Reynolds number up to ten showed significant liquid mal-distribution. Local cell velocities were found to be as high as twenty seven times average velocity, contrary to the previously held belief of four times. The degree of mal-distribution varied across the channel as well as in the vertical direction, and depended on the upward or downward direction of flow. At Reynolds numbers less than one, flow zig-zagged from one side of the channel to the other in wave form, but increases in Reynolds number improved liquid distribution. A detailed examination of selected cells showed velocity variations in different directions, together with variation within individual cells. Experimental results are also reported on the flow split when passing through a single cell in a section of a channel . These observations were used to explain mal-distribution in the perspex channel itself.

Local shell-side coefficients in shell and tube exchanges the use of a mass transfer technique

A diffusion-controlled electrochemical mass transfer technique has been employed in making local measurements of shell-side coefficients in segmentally baffled shell and tube heat exchangers. Corresponding heat transfer data are predicted through the Chilton and Colburn heat and mass transfer analogy. Mass transfer coefficients were measured for baffle spacing lengths of individual tubes in an internal baffle compartment. Shell-side pressure measurements were also made. Baffle compartment average coefficients derived from individual tube coefficients are shown to be in good agreement with reported experimental bundle average heat transfer data for a heat exchanger model of similar geometry. Mass transfer coefficients of individual tubes compare favourably with those obtained previously by another mass transfer technique. Experimental data are reported for a variety of segmental baffle configurations over the shell-side Reynolds number range 100 to 42 000. Baffles with zero clearances were studied at three baffle cuts and two baffle spacings. Baffle geometry is shown to have a large effect on the distribution of tube coefficients within the baffle compartment. Fluid "jetting" is identified with some baffle configurations. No simple characteristic velocity is found to correlate zonal or baffle compartment average mass transfer data for the effect of both baffle cut and baffle spacing. Experiments with baffle clearances typical of commercial heat exchangers are also reported. The effect of leakage streams associated with these baffles is identified. Investigations were extended to double segmental baffles for which no data had previously been published. The similarity in the shell-side characteristics of this baffle arrangement and two parallel single segmental baffle arrangements is demonstrated. A general relationship between the shell-side mass transfer performance and pressure drop was indicated by the data for all the baffle configurations examined.

CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors. Part B: Heat, momentum and mass transport in bubbling fluidised beds

The fluid–particle interaction inside a 150 g/h fluidised bed reactor is modelled. The biomass particle is injected into the fluidised bed and the heat, momentum and mass transport from the fluidising gas and fluidised sand is modelled. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. Heat transfer from the bubbling bed to the discrete biomass particle, as well as biomass reaction kinetics are modelled according to the literature. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase. FLUENT 6.2 has been used as the modelling framework of the simulations with the whole pyrolysis model incorporated in the form of user-defined function (UDF). The study completes the fast pyrolysis modelling in bubbling fluidised bed reactors.

Flow and pressure drop on the shellside of cylindrical heat exchangers

This work presents pressure distributions and fluid flow patterns on the shellside of a cylindrical shell-and-tube heat exchanger. The apparatus used was constructed from glass enabling direct observation of the flow using a dye release technique and had ten traversable pressure instrumented tubes permitting detailed pressure distributions to be obtained. The `exchanger' had a large tube bundle (278 tubes) and main flow areas typical of practical designs. Six geometries were studied: three baffle spacings both with and without baffle leakage. Results are also presented of three-dimensional modelling of shellside flows using the Harwell Laboratory's FLOW3D code. Flow visualisation provided flow patterns in the central plane of the bundle and adjacent to the shell wall. Comparison of these high-lighted significant radial flow variations. In particular, separated regions, originating from the baffle tips, were observed. The size of these regions was small in the bundle central plane but large adjacent to the shell wall and extended into the bypass lane. This appeared to reduce the bypass flow area and hence the bypass flow fraction. The three-dimensional flow modelling results were presented as velocity vector and isobar maps. The vector maps illustrated regions of high and low velocity which could be prone to tube vibration and fouling. Separated regions were also in evidence. A non-uniform crossflow was discovered with, in general, higher velocities in the central plane of the bundle than near the shell wall._The form of the isobar maps calculated by FLOW3D was in good agreement with experimental results. In particular, larger pressure drops occurred across the inlet than outlet of a crossflow region and were higher near the upstream than downstream baffle face. The effect of baffle spacing and baffle leakage on crossflow and window pressure drop measurements was identified. Agreement between the current measurements, previously obtained data and commonly used design correlations/models was, in general, poor. This was explained in terms of the increased understanding of shellside flow. The bulk of previous data, which dervies from small-scale rigs with few tubes, have been shown to be unrepresentative of typical commerical units. The Heat Transfer and Fluid Flow Service design program TASC provided the best predictions of the current pressure drop results. However, a number of simple one-dimensional models in TASC are, individually, questionable. Some revised models have been proposed.

Transfer coefficients over inclined tube banks

The thesis presents experimental results for shell-side transfer coefficients and pressure drops across four different tube banks, using small-scale models, with yawed tubes, as found in many types of heat exchangers, boilers and nuclear reactors. The tube banks investigated have a staggered tube layout on a rotated square pitch, with a 1.25 pitch-to-diameter ratio. The angle of attack was varied between 45o and 90o. An extensive range of Reynolds number, i.e. 0.5. to 12,600, covering so-called laminar, transition and turbulent flows, was investigated. A diffusion-controlled electrochemical mass transfer technique has been employed to measure mass transfer coefficients. The heat transfer coefficients may be then readily obtained from the mass transfer values by applying the well-established Chilton-Colburn analogy. The results for the normal tube bank, which forms the base case for the study on inclined tube banks, show close agreement with previous work. The transfer coefficients and pressure drops of the inclined tube banks are compared with results from the ideal normal tube bank to examine the effect of inclination angle on heat transfer and pressure drop variations. The variation of the transfer coefficients row-by-row and the entrance and exit effects have also been investigated. An auxilary investigation has been carried out on the role of natural convection. A preliminary correlation of transfer coefficients and pressure drops against the variation in the yaw angle has been attempted. The results are discussed in the light of the few existing theoretical treatments and experimental data for these situations, and recommendations made for future work.