Effects of geological inhomogeneity on high Rayleigh number steady state heat and mass transfer in fluid-saturated porous media heated from below

A parametric study is carried out to investigate how geological inhomogeneity affects the pore-fluid convective flow field, the temperature distribution, and the mass concentration distribution in a fluid-saturated porous medium. The related numerical results have demonstrated that (1) the effects of both medium permeability inhomogeneity and medium thermal conductivity inhomogeneity are significant on the pore-fluid convective flow and the species concentration distribution in the porous medium; (2) the effect of medium thermal conductivity inhomogeneity is dramatic on the temperature distribution in the porous medium, but the effect of medium permeability inhomogeneity on the temperature distribution may be considerable, depending on the Rayleigh number involved in the analysis; (3) if the coupling effect between pore-fluid flow and mass transport is weak, the effect of the Lewis number is negligible on the pore-fluid convective flow and temperature distribution, hut it is significant on the species concentration distribution in the medium.

Effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in deformable fluid-saturated porous media heated from below

In this paper, a solution method is presented to deal with fully coupled problems between medium deformation, pore-fluid flow and heat transfer in fluid-saturated porous media having supercritical Rayleigh numbers. To validate the present solution method, analytical solutions to a benchmark problem are derived for some special cases. After the solution method is validated, a numerical study is carried out to investigate the effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in fluid-saturated media when they are heated from below. The related numerical results have demonstrated that: (1) medium thermoelasticity has a little influence on the overall pattern of convective pore-fluid flow, but it has a considerable effect on the localization of medium deformation, pore-fluid flow, heat transfer and mineralization in a porous medium, especially when the porous medium is comprised of soft rock masses; (2) convective pore-fluid flow plays a very important role in the localization of medium deformation, heat transfer and mineralization in a porous medium. (C) 1999 Elsevier Science S.A. All rights reserved.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Simulació numèrica de fenòmens de transferència de calor i dinàmics de fluids transitoris

The work carried out during the 4 year research activity can be barely classified in two main lines. On the one hand, a considerable effort is taken to address issues related with the verification of multi-dimensional and transient solutions that are obtained by numerical simulations. Within the studied cases, we can consider cases of piston-cylinder ows within geometries similar to those of hermetic reciprocating compressors.This issue is mentioned in Part I. On the other hand, numerical simulations of different phenomena have been performed. More emphasis has been given to the natural convection ow within enclosures. This is explained in Part II. The case extensively studied has been the natural convection ow. The natural convection ow within enclosures has attracted the attention of many researchers due to its potential to model numerous applications of engineering interest, such as cooling of electronic devices, air ow in buildings, heat transfer in solar collectors, among others. The natural convection studies corresponding to the parallelepipedic enclosures can be classified into two elementary classes: i) heating from a horizontal wall (heating from below); ii) heating from a vertical wall. The characteristic example of the former case is the Rayleigh-B_enard ow, however this research is on the cavities heated from the side. This configuration is referred commonly as the differentially heated cavity.

Uusien energianormien vaikutus pientalojen rakenteisiin ympäristön kannalta

Tehokkaimpia keinoja vähentää rakennusten lämmitysenergian kulutusta ja lämmityksen aiheuttavia hiilidioksidi- ja happamoitavia päästöjä on tiukentaa rakentamismääräysten lämmöneristysvaatimuksia. Hyvin lämmöneristetyissä, tiiveissä ja ilmanvaihdoltaan optimoiduissa taloissa on pienet lämpöhäviöt. Näin ympäristöä kuormittava vaikutus saadaan paljon vähemmäksi kuin nykynormien mukaisissa asuinrakennuksissa. Johtumislämpöhäviö pienenee suoraan eristekerroksia paksuntamalla ja siihen on helpointa vaikuttaa. Mitä suurempiin eristepaksuuksiin mennään sen suuremmaksi tulee konvektion osuus kokonaislämpöhäviöstä. Tulevaisuudessa parempia ratkaisuja haetaan erityisesti konvektiosta ja säteilystä aiheutuvien lämpöhäviöiden pienentämiseksi. Eristeen osastointi ilmanpitävillä, vesihöyryä diffuusisesti läpäisevillä pystysuuntaisilla konvektiokatkoilla vähentää tehokkaasti paksun seinäeristeen kuljettumis-ilmavirtauksia. Katkoina käytetään erilaisia kalvoja ja rakennuspapereita, joilla on pieni emissiviteetti. Katkojen merkitys kasvaa, kun mennään uusien normien mukaisiin eristepaksuuksiin. Lämmöneriste voidaan toteuttaa myös kokoamalla ohuita kalvoja paketiksi, jotka jakavat ilmatilan ja siis eristeelle varatun paksuuden suljettuihin ilmaväleihin. Kun kalvoiksi valitaan pieniemissiviteettisiä pintoja, saadaan säteilylämmönsiirto lähes eliminoiduksi. Tällaisen ilmatilan lämmönjohtumisluku lähestyy paikallaan pysyvän ilman lämmönjohtumislukua, l = 0,025 W/Km, eli tällä rakennesysteemillä on mahdollista toteuttaa ohuempia rakenteita kuin perinteisillä eristeillä. Hygroskooppisen massan käyttö sisäilman kosteutta tasaavana rakenteena voi olla tulevaisuutta. Kehitystyö tuottaa uusia, kosteusteknisesti toimivia sovelluksia. Toisaalta palomääräykset tulevat kehitystyötä vastaan. Hygroskooppinen pintamateriaali on kevyt (pieni tiheys) ja paloteknisesti arka. Suoraa sähkölämmitystä ei voida pitää ympäristöystävällisenä. Sen jalostusketju on pitkä ja monivaiheinen. Millä peruspolttoaineella sähköä tuotetaan, vaikuttaa asiaan luonnollisestikin. Suoraa sähkölämmitystä voidaan suositella vain yksinäisen ihmisen taloudessa lämmitysmuotona taloudellisista syistä. Halvan polttoaineen säästöllä ei voida maksaa suuria laiteinvestointeja. Aurinkoenergian hyvä hyödyntäminen edellyttää hyvää säätöä, joka kytkee lämmityksen pois päältä silloin, kun aurinko lämmittää. Auringon hetkelliset säteilytehot ovat suuria verrattuna rakenteen lämpöhäviöihin ja huonetilojen lämmöntarpeeseen. Ratkaisu aurinkoenergian hetkellisyyteen ja paikallisuuteen on energian siirtäminen lämmöntarpeen mukaan rakennuksen eri osiin ja sen varastoiminen päivätasolla. Kun varastoivasta massasta ei ole suoraa yhteyttä ulos, voidaan kerääjäeristeeltä saatu lämpö käyttää häviöttömästi huonetilojen lämmittämiseen. Vaikka lämmitysenergian käytössä päästään 30 % vähennyksiin uudisrakennusten osalta, ei kokonaisenergian käyttö merkittävästi pienene, jos taloussähkön kulutus pysyy vakiona. Sama pätee myös CO2 -päästöihin. Saavutettava etu lämmitys-energian kulutuksessa voidaan hukata yhä suurenevaksi taloussähkön käytöksi, mikä olisi erityisen huono asia ympäristön kannalta.

Using a Double-pass solar drier for drying of bamboo shoots

Three different drying methods, a forced convection double-pass solar drier (DPSD), typical cabinet type natural convection solar drier (CD) and traditional open-sun drying (OSD) were used for draying of bamboo shoots in central Vietnam. During drying the operational parameters such as drying temperature, relative humidity, air velocity, insolation and water evaporation have been recorded hourly. The mean drying temperatures and relative humidity in the drying chamber were 55.2°C, 23.7%; 47.5°C, 37,6%; 36.2°C, 47.8% in DPSD, CD and OSD, respectively. The mean global radiation during all experimental runs was 670 Wm^−2. The result also shows that fastest drying process was occurred in DPSD where the falling-rate period was achieved after 7 hours, in change to OSD where it took 16 hours. The overall drying efficiency was 23.11%, 15.83% and 9.73% in case of DPSD, CD and OSD, respectively. Although the construction cost of DPSD was significantly higher than in CD, the drying costs per one kilogram of bamboo shoots were by 42.8% lower in case of DPSD as compared to CD. Double-pass solar drier was found to be technically and economically suitable for drying of bamboo shoots under the specific conditions in central Vietnam and in all cases, the use of this drier led to considerable reduction in drying time in comparison to traditional open-sun drying.

Thermal control of electronics: Perspectives and prospects

One of the most prominent industrial applications of heat transfer science and engineering has been electronics thermal control. Driven by the relentless increase in spatial density of microelectronic devices, integrated circuit chip powers have risen by a factor of 100 over the past twenty years, with a somewhat smaller increase in heat flux. The traditional approaches using natural convection and forced-air cooling are becoming less viable as power levels increase. This paper provides a high-level overview of the thermal management problem from the perspective of a practitioner, as well as speculation on the prospects for electronics thermal engineering in years to come.

Wind-tunnel modelling of dispersion from a scalar area source in urban-like roughness

A wind-tunnel study was conducted to investigate ventilation of scalars from urban-like geometries at neighbourhood scale by exploring two different geometries a uniform height roughness and a non-uniform height roughness, both with an equal plan and frontal density of λ p = λ f = 25%. In both configurations a sub-unit of the idealized urban surface was coated with a thin layer of naphthalene to represent area sources. The naphthalene sublimation method was used to measure directly total area-averaged transport of scalars out of the complex geometries. At the same time, naphthalene vapour concentrations controlled by the turbulent fluxes were detected using a fast Flame Ionisation Detection (FID) technique. This paper describes the novel use of a naphthalene coated surface as an area source in dispersion studies. Particular emphasis was also given to testing whether the concentration measurements were independent of Reynolds number. For low wind speeds, transfer from the naphthalene surface is determined by a combination of forced and natural convection. Compared with a propane point source release, a 25% higher free stream velocity was needed for the naphthalene area source to yield Reynolds-number-independent concentration fields. Ventilation transfer coefficients w T /U derived from the naphthalene sublimation method showed that, whilst there was enhanced vertical momentum exchange due to obstacle height variability, advection was reduced and dispersion from the source area was not enhanced. Thus, the height variability of a canopy is an important parameter when generalising urban dispersion. Fine resolution concentration measurements in the canopy showed the effect of height variability on dispersion at street scale. Rapid vertical transport in the wake of individual high-rise obstacles was found to generate elevated point-like sources. A Gaussian plume model was used to analyse differences in the downstream plumes. Intensified lateral and vertical plume spread and plume dilution with height was found for the non-uniform height roughness

Effect of mountain breeze in ventilating the city under calm and neutral atmospheric environment: Interaction of airflow structures and ventilation efficiency

Many modern cities locate in the mountainous areas, like Hong Kong, Phoenix City and Los Angles. It is confirmed in the literature that the mountain wind system developed by differential heating or cooling can be very beneficial in ventilating the city nearby and alleviating the UHI effect. However, the direct interaction of mountain wind with the natural-convection circulation due to heated urban surfaces has not been studied, to our best knowledge. This kind of unique interaction of two kinds of airflow structures under calm and neutral atmospheric environment is investigated in this paper by CFD approach. A physical model comprising a simple mountain and three long building blocks (forming two street canyons) is firstly developed. Different airflow structures are identified within the conditions of different mountain-building height ratios (R=Hm/Hb) by varying building height but fixing mountain height. It is found that the higher ventilation rate in the street canyons is expected in the cases of smaller mountain-building ratios, indicating the stronger natural convection due to increasing heated building surfaces. However, there is the highest air change rate (ACH) in the lowest-building-height case and most of the air is advective into the street canyon through the top open area, highlighting the important role played by the mountain wind. In terms of the ventilation efficiency, it is shown that the smallest R case enjoys the best air change efficiency followed by the highest R case, while the worst ventilative street canyons occur at the middle R case. In the end, a gap across the streets is introduced in the modeling. The existence of the gap can greatly channel the mountain wind and distribute the air into streets nearby. Thus the ACH can be doubled and air quality can be significantly improved.

Desenvolvimento de um programa de simulação computacional de sistemas de aquecimento solar para água

Esta Tese apresenta uma análise do comportamento térmico de um sistema de aquecimento solar operando por termossifão. Neste tipo de sistema o fluido no coletor solar é circulado por convecção natural, que acontece devido à diferença de massa específica da água ao longo circuito. Nestes sistemas a vazão mássica varia ao longo do dia e do ano, dependendo, dentre outros fatores, da irradiância solar absorvida, do perfil de temperaturas da água no sistema, da geometria, do volume e do perfil de demanda de água quente. Para uma avaliação detalhada do comportamento térmico de aquecedores solares operando por termossifão foram realizados ensaios experimentais e cálculos teóricos. Os resultados dos experimentos concordaram com aqueles apresentados na literatura e sua análise fundamentou o desenvolvimento do aplicativo TermoSim, um programa de simulação computacional do comportamento térmico de sistemas de aquecimento de água com energia solar. O tratamento matemático adotado no TermoSim compreende a modelagem dos coletores solares de acordo com a teoria de Hottel-Bliss-Whillier. O reservatório térmico é modelado com estratificação térmica, convecção e condução entre as camadas. A vazão mássica é obtida a partir do balanço da quantidade de movimento no circuito. Os modelos matemáticos empregados na construção do aplicativo TermoSim foram validados através do confronto dos resultados simulados com medidas experimentais. Foi demonstrado que a utilização destes modelos é adequada e permite reproduzir com precisão o comportamento térmico dos coletores solares e do reservatório térmico. Além do programa TermoSim, foi também desenvolvido o programa TermoDim, que é uma ferramenta para o dimensionamento de sistemas de aquecimento solar, que requer apenas o conhecimento dos parâmetros geométricos do sistema, dados meteorológicos em média mensal e informação a respeito do volume de demanda. O TermoDim é apropriado para estimar o desempenho de aquecedores solares operando por termossifão com tanques verticais e horizontais. O método de dimensionamento do TermoDim é baseado na correlação para a eficiência média mensal obtida neste trabalho a partir de um grande número de simulações.