Comparing the thermal performance of horizontal slinky-loop and vertical slinky-loop heat exchangers

The heat pump market in the UK has grown rapidly over the last few years. Performance analyses of vertical ground-loop heat exchanger configurations have been widely carried out using both numerical modelling and experiments. However, research findings and design recommendations on horizontal slinky-loop and vertical slinky-loop heat exchangers are far fewer compared with those for vertical ground-loop heat exchanger configurations, especially where the long-term operation of the systems is concerned. The paper presents the results obtained from a numerical simulation for the horizontal slinky-loop and vertical slinky-loop heat exchangers of a ground-source heat pump system. A three-dimensional numerical heat transfer model was developed to study the thermal performance of various heat exchanger configurations. The influence of the loop pitch (loop spacing) and the depth of a vertical slinky-loop installation were investigated and the thermal performance and excavation work required for the horizontal and vertical slinky-loop heat exchangers were compared. The influence of the installation depth for vertical slinky-loop configurations was also investigated. The results of this study show that the influence of the installation depth of the vertical slinky-loop heat exchanger on the thermal performance of the system is small. The maximum difference in the thermal performance between the vertical and horizontal slinky-loop heat exchangers with the same loop diameter and loop pitch is less than 5%.

Foreword: International Space Science Institute (ISSI) workshop on the Earth’s cryosphere and sea level change

Rising sea level is perhaps the most severe consequence of climate warming, as much of the world’s population and infrastructure is located near current sea level (Lemke et al. 2007). A major rise of a metre or more would cause serious problems. Such possibilities have been suggested by Hansen and Sato (2011) who pointed out that sea level was several metres higher than now during the Holsteinian and Eemian interglacials (about 250,000 and 120,000 years ago, respectively), even though the global temperature was then only slightly higher than it is nowadays. It is consequently of the utmost importance to determine whether such a sea level rise could occur and, if so, how fast it might happen. Sea level undergoes considerable changes due to natural processes such as the wind, ocean currents and tidal motions. On longer time scales, the sea level is influenced by steric effects (sea water expansion caused by temperature and salinity changes of the ocean) and by eustatic effects caused by changes in ocean mass. Changes in the Earth’s cryosphere, such as the retreat or expansion of glaciers and land ice areas, have been the dominant cause of sea level change during the Earth’s recent history. During the glacial cycles of the last million years, the sea level varied by a large amount, of the order of 100 m. If the Earth’s cryosphere were to disappear completely, the sea level would rise by some 65 m. The scientific papers in the present volume address the different aspects of the Earth’s cryosphere and how the different changes in the cryosphere affect sea level change. It represents the outcome of the first workshop held within the new ISSI Earth Science Programme. The workshop took place from 22 to 26 March, 2010, in Bern, Switzerland, with the objective of providing an in-depth insight into the future of mountain glaciers and the large land ice areas of Antarctica and Greenland, which are exposed to natural and anthropogenic climate influences, and their effects on sea level change. The participants of the workshop are experts in different fields including meteorology, climatology, oceanography, glaciology and geodesy; they use advanced space-based observational studies and state-of-the-art numerical modelling.

On the 'downward control' of extratropical diabatic circulations by eddy-induced mean zonal forces

The situation considered is that of a zonally symmetric model of the middle atmosphere subject to a given quasi-steady zonal force F̄, conceived to be the result of irreversible angular momentum transfer due to the upward propagation and breaking of Rossby and gravity waves together with any other dissipative eddy effects that may be relevant. The model's diabatic heating is assumed to have the qualitative character of a relaxation toward some radiatively determined temperature field. To the extent that the force F̄ may be regarded as given, and the extratropical angular momentum distribution is realistic, the extratropical diabatic mass flow across a given isentropic surface may be regarded as controlled exclusively by the F̄ distribution above that surface (implying control by the eddy dissipation above that surface and not, for instance, by the frequency of tropopause folding below). This “downward control” principle expresses a critical part of the dynamical chain of cause and effect governing the average rate at which photochemical products like ozone become available for folding into, or otherwise descending into, the extratropical troposphere. The dynamical facts expressed by the principle are also relevant, for instance, to understanding the seasonal-mean rate of upwelling of water vapor to the summer mesopause, and the interhemispheric differences in stratospheric tracer transport. The robustness of the principle is examined when F̄ is time-dependent. For a global-scale, zonally symmetric diabatic circulation with a Brewer-Dobson-like horizontal structure given by the second zonally symmetric Hough mode, with Rossby height HR = 13 km in an isothermal atmosphere with density scale height H = 7 km, the vertical partitioning of the unsteady part of the mass circulation caused by fluctuations in F̄ confined to a shallow layer LF̄ is always at least 84% downward. It is 90% downward when the force fluctuates sinusoidally on twice the radiative relaxation timescale and 95% if five times slower. The time-dependent adjustment when F̄ is changed suddenly is elucidated, extending the work of Dickinson (1968), when the atmosphere is unbounded above and below. Above the forcing, the adjustment is characterized by decay of the meridional mass circulation cell at a rate proportional to the radiative relaxation rate τr−1 divided by {1 + (4H2/HR2)}. This decay is related to the boundedness of the angular momentum that can be taken up by the finite mass of air above LF̄ without causing an ever-increasing departure from thermal wind balance. Below the forcing, the meridional mass circulation cell penetrates downward at a speed τr−1 HR2/H. For the second Hough mode, the time for downward penetration through one density scale height is about 6 days if the radiative relaxation time is 20 days, the latter being representative of the lower stratosphere. At any given altitude, a steady state is approached. The effect of a rigid lower boundary on the time-dependent adjustment is also considered. If a frictional planetary boundary layer is present then a steady state is ultimately approached everywhere, with the mass circulation extending downward from LF̄ and closing via the boundary layer. Satellite observations of temperature and ozone are used in conjunction with a radiative transfer scheme to estimate the altitudes from which the lower stratospheric diabatic vertical velocity is controlled by the effective F̄ in the real atmosphere. The data appear to indicate that about 80% of the effective control is usually exerted from below 40 km but with significant exceptions up to 70 km (in the high latitude southern hemispheric winter). The implications for numerical modelling of chemical transport are noted.

The effect of model configuration on modelled hillslope–riparian interactions

The transfer of hillslope water to and through the riparian zone forms a research area of importance in hydrological investigations. Numerical modelling schemes offer a way to visualise and quantify first-order controls on catchment runoff response and mixing. We use a two-dimensional Finite Element model to assess the link between model setup decisions (e.g. zero-flux boundary definitions, soil algorithm choice) and the consequential hydrological process behaviour. A detailed understanding of the consequences of model configuration is required in order to produce reliable estimates of state variables. We demonstrate that model configuration decisions can determine effectively the presence or absence of particular hillslope flow processes and, the magnitude and direction of flux at the hillslope–riparian interface. If these consequences are not fully explored for any given scheme and application, the resulting process inference may well be misleading.

The latitudinal effect of corotating interaction regions on galactic cosmic rays

The response of galactic cosmic rays (GCRs) to an isolated enhancement of the non-axisymmetric component of the solar open magnetic field between June and November 1996 is investigated by using a combination of solar observations and numerical modelling of the interplanetary medium. The most obvious coronal hole visible from Earth at the time had little shielding effect on the flux of GCRs, as measured at Earth by neutron monitors. It is found that the evolution of the corotating interaction regions generated by a less obvious coronal hole was the principal controlling factor. Moreover, we demonstrate the imprint of the latitudinal and longitudinal evolution of that coronal hole on the variation of GCRs. The latitudinal extent of this solar minimum corotating interaction region had a determining, but local, shielding effect on GCRs, confirming previous modelling results.

Centennial changes in solar activity and the response of galactic cosmic rays

There is a growing consensus that the eleven year modulation of galactic cosmic rays (GCRs) resulting from solar activity is related to interplanetary propagating diffusive barriers (PDBs). The source of these PDBs is not well understood and numerical models describing GCR modulation simulate their effect by scaling the diffusion tensor to the interplanetary magnetic field strength (IMF). The implications of a century-scale change in solar wind speed and open solar flux, for numerical modelling of GCR modulation and the reconstruction of GCR variations over the last hundred years are discussed. The dominant role of the solar non-axisymmetric magnetic field in both forcing longitudinal solar wind speed fluctuations at solar maximum and in increasing the IMF is discussed in the context of a long-term rise in the open solar magnetic flux.

A review of urban roughness sublayer turbulence

It is becoming increasingly important that we can understand and model flow processes in urban areas. Applications such as weather forecasting, air quality and sustainable urban development rely on accurate modelling of the interface between an urban surface and the atmosphere above. This review gives an overview of current understanding of turbulence generated by an urban surface up to a few building heights, the layer called the roughness sublayer (RSL). High quality datasets are also identified which can be used in the development of suitable parameterisations of the urban RSL. Datasets derived from physical and numerical modelling, and full-scale observations in urban areas now exist across a range of urban-type morphologies (e.g. street canyons, cubes, idealised and realistic building layouts). Results show that the urban RSL depth falls within 2 – 5 times mean building height and is not easily related to morphology. Systematic perturbations away from uniform layouts (e.g. varying building heights) have a significant impact on RSL structure and depth. Considerable fetch is required to develop an overlying inertial sublayer, where turbulence is more homogeneous, and some authors have suggested that the “patchiness” of urban areas may prevent inertial sublayers from developing at all. Turbulence statistics suggest similarities between vegetation and urban canopies but key differences are emerging. There is no consensus as to suitable scaling variables, e.g. friction velocity above canopy vs. square root of maximum Reynolds stress, mean vs. maximum building height. The review includes a summary of existing modelling practices and highlights research priorities.

SHIMMER (1.0): a novel mathematical model for microbial and biogeochemical dynamics in glacier forefield ecosystems

SHIMMER (Soil biogeocHemIcal Model for Microbial Ecosystem Response) is a new numerical modelling framework designed to simulate microbial dynamics and biogeochemical cycling during initial ecosystem development in glacier forefield soils. However, it is also transferable to other extreme ecosystem types (such as desert soils or the surface of glaciers). The rationale for model development arises from decades of empirical observations in glacier forefields, and enables a quantitative and process focussed approach. Here, we provide a detailed description of SHIMMER, test its performance in two case study forefields: the Damma Glacier (Switzerland) and the Athabasca Glacier (Canada) and analyse sensitivity to identify the most sensitive and unconstrained model parameters. Results show that the accumulation of microbial biomass is highly dependent on variation in microbial growth and death rate constants, Q10 values, the active fraction of microbial biomass and the reactivity of organic matter. The model correctly predicts the rapid accumulation of microbial biomass observed during the initial stages of succession in the forefields of both the case study systems. Primary production is responsible for the initial build-up of labile substrate that subsequently supports heterotrophic growth. However, allochthonous contributions of organic matter, and nitrogen fixation, are important in sustaining this productivity. The development and application of SHIMMER also highlights aspects of these systems that require further empirical research: quantifying nutrient budgets and biogeochemical rates, exploring seasonality and microbial growth and cell death. This will lead to increased understanding of how glacier forefields contribute to global biogeochemical cycling and climate under future ice retreat.

Modelagem numérica de estruturas de contenção atirantadas em areia

A numerical study on the behavior of tied-back retaining walls in sand, using the finite element method (FEM) is presented. The analyses were performed using the software Plaxis 2D, and were focused on the development of horizontal displacements, horizontal stresses, shear forces and bending moments in the structure during the construction process. Emphasis was placed on the evaluation of wall embedment, tie-back horizontal spacing, wall thickness, and free anchor length on wall behavior. A representative soil profile of a specific region at the City of Natal, Brazil, was used in the numerical analyses. New facilities built on this region often include retaining structures of the same type studied herein. Soil behavior was modeled using the Mohr-Coulomb constitutive model, whereas the structural elements were modeled using the linear elastic model. Shear strength parameters of the soil layers were obtained from direct shear test results conducted with samples collected at the studied site. Deformation parameters were obtained from empirical correlations from SPT test results carried out on the studied site. The results of the numerical analyses revealed that the effect of wall embedment on the investigated parameters is virtually negligible. Conversely, the tie-back horizontal spacing plays an important role on the investigated parameters. The results also demonstrated that the wall thickness significantly affects the wall horizontal displacements, and the shear forces and bending moments within the retaining structure. However, wall thickness was not found to influence horizontal stresses in the structure

Mecanismos de alojamento de magmas granitóides :exemplo do plúton de Japi(RN)

This dissertation describes the igneous suites of the Japi granitoid pluton, intrusive in the Paleoproterozoic gneiss-migmatite complex of the eastern domain of the Seridó Belt, northeastern Brazil. Field relations show that the pluton is affected by strong deformation associated to the Brasiliano orogeny (known as the D3 phase) , with a NW-trending extensionalleft-hand senestral shear zone (the Japi Shear Zone, JSZ) bordering the intrusive body to the west. Four plutonic suites are found in the main pluton and as satellyte intrusions, besides Iate pegmatite and pink leucogranites. An alkaline granitoid suite, dominated by syenogranites bearing sodic augite (and subordinate hornblende), define a main elliptical intrusion. In its northern part, this intrusion is made up by concentric sheets, contrasting with a smaller rounded stock to the south. These granites display a pervasive solid-state S>L fabric developed under high T conditions, characterized by plastic deformation of quartz and feldspar. It is especially, developed along the border of the pluton, with inward dips. A regular magmatic layering is present sometimes, parallel to the tectonic foliation. The syntectonic emplacement as regards to the Brasiliano (D3) event is indicated by the common occurrence of dykes and sheets along transtensional or extensional sites of the major structure. Field relations attest to the early emplacement of the alkaline granites as regards to the other suites. A basic-to-intermediate suite occurs as a western satellyte body and occupying the southern tail of the main alkaline pluton. It comprises a wide variety of compositional terms, including primitive gabbros and gabbro-norites, differentiated to monzonitic intermediate facies containing amphibole and biotite as their main mafic phases. These rocks display transitional high-K calc-alkaline to shoshonitic affinities. Porphyritic monzogranite suítes commonly occur as dykes and minor intrusives, isolated or associated with the basic-tointermediate rocks. In the latter case, magma mingling and mixing features attest that these are contemporaneous igneous suites. These granites show K-feldspar phenocrysts and a hornblende+biotite+titanite assemblage, displaying subalkaline/monzonitic geochemical affinities. Both suites exhibit SL magmatic fabrics overprinting or transitional to solid-state D3 deformation related to the JSI. Chemical data clearly show that they are related to different parental magmas. Finally, a microgranite suite occurs along a few topographic ridges paralell to the JSI. It comprises dominantly granodiorites with a mineralogy similar to the one of the porphyritic granitoids. However, discriminant diagrams show their distinct calc-alkaline affinity. The granodiorites display an essencially magmatic fabric, even though an incipient D3 solid-state structure may be developed along the JSI. Intrusion relationships with the previous suites, as well as regards to the D3 structures, point to their Iate emplacement. All these suites are intrusive in a Paleoproterozoic, high-grade gneiss-migmatite complex affected by two previous deformation phases (D1, D2). The fabrics associated with these earlier events are folded and overprinted by the younger D3 structures along the JSZ. The younger deformation is characterized by NE-dipping foliations and N/NE-plunging stretching lineations. In the JSZ northern termination the foliation acquires an ENE orientation, containing a stretching lineation plunging to the south. Symmetric kinematic cri teria developed at this site confirms the transpressional termination of the JSZ, as also shown by orthorrombic quartz c-axis patterns. E-W-trending d extra I shear zones developed in the central part of the JSZ are interpreted as antithetic structures associated to the transtensional deformation along the JSZ. This is consistent with its extensional-transcurrent kinematics and a flat-and-ramp geometry at depth, as shown by gravimetric data. The lateral displacement of the negative residual Bouguer anomalies, as regards to the main outcropping alkaline pluton, may be modelized by other deeper-seated granite bodies. Based on numerical modelling it was possible to infer two distinct intrusion styles for the alkaline pluton. The calculated model values are consistent with an emplacement by sheeting for the northern body, as already suggested by satellyte imagery and field mapping. On the other hand, the results point to a transition towards a diapir-related style associated to the smaller. southern stock. This difference in intrusion styles may relate to intensity variations and transtensional sites of the shear deformation along the JSZ. Trace element and Sr and Nd isotopes of the alkaline granites are compatible with their derivation trom a more basic crustal source, as compared to the presently outcropping highgrade gneisses, with participation (or alternatively dominated by) of an enriched lithospheric mantle component. Like other igneous suites in the Seridó Belt, the high LlL contents and fractionated REE patterns of the basic rocks also point to an enriched mantle as the source for this kind of magmatism. Geochemical and isotope data are compatible with a lower crustal origin for the porphyritic granites. On the basis of the strong control of the JSZ on the emplacement of lower crustal (porphyritic and alkaline granites) or lithospheric mantle (basic rocks, alkaline granites or a component of them) magmas, one may infer a deep root for this structure, bearing an important role in magma extraction, transport and emplacement in the Japi region, eastern domain of the Seridó Belt

Modelagem numérica da superfície na geração da planialtimetra segundo dados de quatro equipamentos topogríficos

The aim of this present work was to compare planialtimetric charts obtained from different risings using two different theodolite types, a total station, and a precision level, used as control. Using a total station, an area was marked with clear variations of relief, following a grid, with a distance of 20 meters among stakes. After that, the stakes were read by the total station and two theodolites of different precisions. The geometric leveling was done by a precision level. The data were input in the DataGeosis software and the numerical modelling of the land was made with mesh of maximum rigidity, generating planialtimetric representation for each rising. It was verified, through comparison of the four representations that little variations occur in relation to the control. The closest representation of the control was the planialtimetry based on the data from the total station, in which the representations obtained from the theodolites were identical among themselves. It was concluded that in the process of obtaining detailed planialtimetry of small areas, submitted to the grid, it was not necessary to use composed geometric leveling, reducing work to the exclusive use of a total station or conventional theodolite.

Modelagem direta de dados de eletrorresistividade 3-D

Dentre os métodos geofísicos, o da Eletrorresistividade é um dos mais tradicionais, com o seu desenvolvimento ocorrido há mais de 80 anos. Durante esse tempo o seu uso acompanhou o avanço no poder de processamento numérico e mais recentemente, a modelagem e inversão tridimensional tornou-se uma possibilidade para o geofísico. Apresentamos, neste trabalho, a técnica de elementos finitos aplicada ao método da eletrorresistividade 3-D, através do cálculo do potencial secundário. Para o desenvolvimento da metodologia, simulamos o levantamento do método da eletrorresistividade 3-D com os arranjos Dipolo-Dipolo e Schlumberger, visando medir as variações laterais e verticais da resistividade aparente do solo. Estes arranjos consistem na injeção de corrente elétrica na superfície e de medidas de diferenças de potencial elétrico, resultante da interação da corrente elétrica com o solo. Sendo que, as fontes e receptores são localizados de acordo com os arranjos escolhidos para o levantamento. Neste trabalho, as curvas de sondagem e as pseudo-secções de resistividade aparente, são obtidas através da modelagem de eletrorresistividade 3-D, usando malha de elementos finitos regular. Para efeito de validação, os resultados são comparados com a resposta 3-D obtida a partir dos potenciais totais.

Analysis of the Stresses in Corrugated Sheets under Bending

This paper presents three different numerical models for the evaluation of the stresses in corrugated sheets under bending. Regarding the numerical simulations different approaches can be considered, i.e., a elastic linear analysis or a physical nonlinear analysis, that considers criteria to fail for the sheet material. Moreover, the construction of the finite element mesh can be used shell elements or solid elements. The choice of each finite element must be made from the consideration of their representativity before behavior to be simulated. Thus, the numerical modelling in this manuscript was performed from the three-dimensional models using the SAP2000Nonlinear software, version 7.42, which has as base the finite elements method (FEM). It was considered shell elements in the build the mesh of finite elements and an analysis of type elastic linear in this case. Five mm thick sheets were evaluated considering three different longitudinal dimensions (spans), i.e., 1100 mm, 1530 mm and 1830 mm. The applied load to the models was 2500 N/m and it was verified that the spans of support of sheets have a significant influence on the results of stresses. The sheets with larger spans present larger stresses for the same applied load. The most intense values of tension occur in the troughs (low waves) of the sheets, on the lower surface, while the most intense values of compression occur in the crests (high waves), on the upper surface of the sheet. The flanks, which are the parts among the troughs and crests of the sheets, are submitted to low levels of stresses. The numeric results of the stresses showed a good agreement with the results obtained from other researchers(3) and these results can be used to predict the behavior of corrugated sheets under bending.

Impacts of interruption of the Agulhas leakage on the tropical Atlantic in coupled ocean-atmosphere simulations

In this paper we use a coupled ocean-atmosphere model to investigate the impact of the interruption of Agulhas leakage of Indian ocean water on the tropical Atlantic, a region where strong coupled ocean-atmosphere interactions occur. The effect of a shut down of leakage of Indian ocean water is isolated from the effect of a collapse of the MOC. In our experiments, the ocean model is forced with boundary conditions in the southeastern corner of the domain that correspond to no interocean exchange of Indian ocean water into the Atlantic. The southern boundary condition is taken from the Levitus data and ensures an MOC in the Atlantic. Within this configuration, instead of warm and salty Indian ocean water temperature (cold) and salinity (fresh) anomalies of southern ocean origin propagate into the South Atlantic and eventually reach the equatorial region, mainly in the thermocline. This set up mimics the closure of the ""warm water path"" in favor of the ""cold water path"". As part of the atmospheric response, there is a northward shift of the intertropical convergence zone (ITCZ). The changes in trade winds lead to reduced Ekman pumping in the equatorial region. This leads to a freshening and warming of the surface waters along the equator. Especially in the Cold Tongue region, the cold and fresh subsurface anomalies do not reach the surface due to the reduced upwelling. The anomaly signals are transported by the equatorial undercurrent and spread away from the equator within the thermocline. Part of the anomaly eventually reaches the Tropical North Atlantic, where it affects the Guinea Dome. Surprisingly, the main effect at the surface is small on the equator and relatively large at the Guinea Dome. In the atmosphere, the northward shift of the ITCZ is associated with a band of negative precipitation anomalies and higher salinities over the Tropical South Atlantic. An important implication of these results is that the modified water characteristics due to a shut down of the Agulhas leakage remain largely unaffected when crossing the equatorial Atlantic and therefore can affect the deepwater formation in the North Atlantic. This supports the hypothesis that the Agulhas leakage is an important source region for climate change and decadal variability of the Atlantic.