Search

Current hygrothermal behaviour prediction models neglect the hysteresis phenomenon. This leads to a discrepancy between numerical and experimental results, and a miscalculation of buildings’ durability. In this paper, a new mathematical model of hysteresis is proposed and implemented in a hygrothermal model to reduce this discrepancy. The model is based on a symmetry property between sorption curves and uses also a homotopic transformation relative to a parameter [Formula: see text]. The advantage of this model lies in its ease of use and implementation since it could be applied with the knowledge of only one main sorption curve by considering [Formula: see text], in other words, we only use the axisymmetric property here. In the case where the other main sorption curve is known, we use this curve to incorporate the homotopy property in order to calibrate the parameter [Formula: see text].The full version of the proposed model is called Axisymmetric + Homotopic. Furthermore, it was compared not only with the experimental sorption curves of different types of materials but also with a model that is well known in the literature (CARMELIET’s model). This comparison shows that the Axisymmetric + Homotopic model reliably predicts hysteresis loops of various types of materials even with the knowledge of only one of the main sorption curves. However, the full version of Axisymmetric + Homotopic model is more reliable and covers a large range of materials. The proposed model was incorporated into the mass transfer model. The simulation results strongly match the experimental ones.

International audience; Thanks to the use of filters collected daily in the west of Paris, we confirm the passage of the lead-laden plume following the fire on the roof and spire of Notre-Dame Cathedral in April 2019. The measured concentrations on the filter (Pb = 1.4 μg m–3), scanning electron microscopy and Hysplit simulation, correspond to an estimate of a few hundred kilograms of lead that would have been volatilized and then oxidized in the form of micronic and submicronic aerosols. The concentrations found in the plume are, however, much lower than those found in the environment in the 1980s and 1990s but are 100 times higher than those found in times prior to and after the fire. The isotopic signature of the plume is almost identical to that of the fine Pb dust found inside the Notre-Dame building. It is different from the isotopic signature of the Parisian atmosphere before and after the fire, but it is similar to that of the atmospheric Pb legacy recorded by peat cores over the last 300 years in France. The presence of very fine lead-bearing particles makes them potentially transportable over long distances after large urban fires. Our study shows the value of daily aerosol sampling to retrospectively trace the plumes of air pollutants from industrial accidents but also from historical monument fires such as Notre-Dame in 2019.

Convergence zones are marked by a variety of settings that may follow each other in modern-day tectonics, including compressive phases such as subduction, obduction, collision but also extensive ones such as back-arc opening or stress-relaxation during orogenesis. Hence, the protracted evolution leading to a super-continent block amalgamation may be difficult to decipher and so may be the forcings on external enveloppes such as volcanism or erosion caused by the different phases.This question arises critically at the time of the Pan-African Orogenesis (1-0.5 Ga) assembling Gondwana, a time of supposedly dramatic and diachronical changes for external envelopes: glaciations of debated scales, deposition of various Banded Iron Formations, first (Ediacarian) fauna, replacement by Cambrian faunas. Our goal is to explore in detail the geodynamical succession leading to the amalgamation of blocks along the northern margin of the West African Craton (WAC), outcropping in the Central Anti-Atlas region, Morocco. This region is characterized by the occurrence of extended convergence-related magmatism, ophiolite emplacement and basins fillings (including BIF) during Cryogenian and Ediacaran periods.Data obtained from compilation of cartographic work, whole-rock geochemistry and datation reveals a polyphased but still poorly constrained evolution through proxies of continentality (εNd) and of crustal thickness (Sr/Y ratio). We present new data spanning metamorphic petrology, basin stratigraphy, coupled datation and trace element analysis in detrital zircons in order to better understand the evolution of the geodynamic, magmatic and drainage systems. We propose a geodynamic scenario based on these data: Development of an early oceanic arc (760-720 Ma) with juvenile magmatic signature (3

International audience; Fault Detection and Diagnosis (FDD) are important tools to perform on-going monitoring of the systems and help in their building commissioning. An innovative method is investigated based on combined data-driven and knowledge-based approaches. This article presents the method. First phase, a so-called operating map of the system is built using dimension reduction method and numerical or experimental dataset. This map is composed of several regions corresponding to nominal operation and to specific faults. The second phase focuses on the FDD. The monitored data are projected on the map. According to the position, a clear and precise fault detection and diagnosis can be carried. The method is applied to an air handling unit. The map is built using data generated with a building simulation program. The reliability of the method is proven using experimental data of nominal and fault operation generated.

Citizenship as a Commodity – of Golden Passports and the European Union “Golden passport” schemes are increasingly gaining popularity around the world. Meanwhile, this trend has also reached the European Union. It now threatens to lead to a partial commercialisation of both national citizenship and the European citizenship. This contribution examines the evolution of national citizenship law through “golden passports” and addresses the question of the compatibility of such programs with the law of the European Union. It thereby reveals the difficulty of reconciling the sovereignty of member states in matters of nationality with the principle of sincere cooperation in the EU.

The mafic-ultramafic Ougda magmatic complex is located in the west part of Tuarge Shield, in Algeria, between Tassendjanet terrane in the east and Ahnet terrane in the west. It is composed of three successive generations of magmatic rocks (Dostal et al., 1996). The first generation located in the north, includes ultramafic rocks cut by dikes of cumulate garnet-bearing mafic rocks and quartz diorite sheets. It records high-temperature metamorphic conditions, granulite facies. The second and third generation located in the south, includes undeformed cumulate and non-cumulate gabbros and intermediate to mafic dikes. The three generations record a geochemical evolution from tholeiitic to calco-alkaline magmatism with subduction-related oceanic environment (Dostal et al., 1996). The age of the first generation is around 800 Ma and the second generation is dated at 680 Ma, considered as the ages of the inception to demise of the oceanic lithosphere (Dostal et al., 1996; Caby and Monié, 2003). Here, we focus on garnet-bearing rocks that show particular interest, as they are affected by high-grade metamorphism in this area. Understanding the pressure-temperature (P-T) evolution of those garnet-bearing rocks allow a crucial constrain of the evolution of the oceanic crust in this area during the Panafrican orogeny.Petrographical investigation shows that all samples share similar mineralogical assemblages with garnet, plagioclase, amphibole, clinopyroxene, ilmenite and rutile. It is interpreted as typical of granulite facies. Garnet is the most dominate phase and show different textural types: Pokioblastic garnet with inclusions of amphibole, clinopyroxene, plagioclase, ilmenite and rutile. In some samples, garnet is very large (~2 cm), ilmenite is observed in garnet core and rutile appears with ilmenite in garnet rims. Clinopyroxene in garnet is a primary phase as it is surrounded by amphibole, which indicate a reaction with garnet. Garnet corona is around clinopyroxene and plagioclase and both are not in contact with each other. Modeling phase relationship using P-T pseudosections was calculated to constrain the P-T conditions and mineralogical evolution. For garnet growth, modal calculations with observed mineral assemblages are more consistent with a solid-state reaction where clinopyroxene and plagioclase are consumed to produce garnet. The PT path manifest with either cooling at high pressure or pressure increase stage, linked to garnet growth, 14-7 Kbar and 1000-700 °C. The P-T conditions are limited by the appearance of biotite at low temperature, solidus at high temperature and olivine at low pressure. The maximum pressure being recorded by rutile-ilmenite-bearing assemblage. This granulitisation stage is followed by a decompression in subsolidus conditions, amphibolites facies, where amphibole appears either as the product of clinopyroxene transformation or reaction between primary clinopyroxene and garnet through hydration. Lastly, hydration in low grade, greenschist facies, is recorded in garnet- and clinopyroxene-free domains with hydrous phases, chlorite, epidote and amphibole. Hence, P-T evolution recorded in garnet-bearing rocks of Ougda shows an anticlockwise PT path with granulitisation stage showing P-T peak recorded by rutile-ilmenite-bearing assemblage in garnet. Followed by a decompression in amphibolite facies with production of amphibole and ended up with late hydration in geenschist facies.

The Moroccan Anti-Atlas orogenic belt encloses several Precambrian inliers comprising two major Neoproterozoic ophiolitic complexes: the Sirwa and Bou Azzer ophiolites. These ophiolites expose crustal and mantle units, thrusting over fragments of a long-lived intra-oceanic arc system. We present a detailed geochronological and petro-geochemical study of three mafic/ultramafic units of these two ophiolites: the Khzama sequence (Sirwa ophiolite) and the Northern and Southern Aït Ahmane sequences (Bou Azzer ophiolite). The crystallization of layered metagabbros from the Bou Azzer ophiolite (North Aït Ahmane sequence) has been dated here at 759 ± 2 Ma (U-Pb on zircons). This new age for the Bou Azzer ophiolite is similar to the formation of the Sirwa ophiolite (762 Ma) and suggests that both units formed during the same spreading event. Metabasalts of the three units show tholeiitic signature but with variable subduction-related imprints marked by LILE enrichments, HFSE depletions and variable Ti contents, similar to modern back-arc basin basalts (BABB). Their back-arc origin is also supported by the geochemical signature of ultramafic units showing very low contents in major and trace incompatible elements (Al2O3: 0.12–1.53 wt%, Ti: 3.5–64.2 ppm and Nb: 0.004–0.10 ppm), attesting of a highly refractory protolith. This is in agreement with the high Cr (0.44–0.81) and low to intermediate Mg (0.25–0.73) of their constitutive Cr-spinels. Dynamic melting models suggest that these serpentinites experienced intense and polyphased hydrous melting events, strongly influenced by supra-subduction zone SSZ-fluid influx and subduction-related melt percolation. Being particularly affected by these SSZ-melt/rock interactions and closer to arc units to the south, the Sirwa ophiolite and the South Aït Ahmane unit of the Bou Azzer ophiolite likely represent an early stage of the arc-back-arc system, which has been more influenced by the magmatic products of the arc activity compared to the North Aït Ahmane unit of the Bou Azzer ophiolite., SCOPUS: ar.j, info:eu-repo/semantics/published

The design of numerical tools to model the behavior of building materials is a challenging task. The crucial point is to save computational cost and maintain high accuracy of predictions. There are two main limitations on the time scale choice, which put an obstacle to solve the above issues. First one is the numerical restriction. A number of research is dedicated to overcome this limitation and it is shown that it can be relaxed with innovative numerical schemes. The second one is the physical restriction. It is imposed by properties of a material, phenomena itself and corresponding boundary conditions. This work is focused on the study of a methodology that enables to overcome the physical restriction on the time grid. So-called Average Reduced Model (ARM) is suggested. It is based on smoothing the time-dependent boundary conditions. Besides, the approximate solution is decomposed into average and fluctuating components. The primer is obtained by integrating the equations over time, whereas the latter is an user-defined empirical model. The methodology is investigated for both heat diffusion and coupled heat and mass transfer. It is demonstrated that the signal core of the boundary conditions is preserved and the physical restriction can be relaxed. The model proved to be reliable, accurate and efficient also in comparison with the experimental data of two years. The implementation of the scarce time-step of $1 \, \sf{h}$ is justified. It is shown, that by maintaining the tolerable error it is possible to cut computational effort up to almost four times in comparison with the complete model with the same time grid.

International audience; In this paper, the spectral method is developed as a reduced-order model for the solution of parametric problems within the building refurbishment framework. We propose to use the spectral reduced-order method to solve parametric problems in an innovative way, integrating the unknown parameter as one of the coordinates of the decomposition. The residual is minimized combining the Tau–Galerkin method with the Collocation approach. The developed method is evaluated in terms of accuracy and reduction of the computational time in three different cases. The dynamic behaviour of unidimensional moisture diffusion is investigated. The cases focus on solving parametric problems in which the solution depends on space, time, diffusivity and material thickness. Results highlight that the parametric spectral reduced-order method provides accurate solutions and can reduce 10 times the degree of freedom of the solution. It allows efficient computation of the physical phenomena with a lower error when compared to traditional approaches

Within the environmental context, numerical modeling is a promising approach to assessing the energy efficiency of buildings. Resilient buildings need to be designed, and capable of adapting to future extreme heat. Simulations are required assuming a one-dimensional heat transfer problem through walls and a simulation horizon of several years (nearly 30). The computational cost associated with such modeling is quite significant and model reduction methods are worth investigating. The objective is to propose a reliable reduced-order model for such long-term simulations. For this, an alternative model reduction approach is investigated, assuming a known Proper Orthogonal Decomposition reduced basis for time, and not for space as usual. The model enables computing parametric solutions using basis interpolation on the tangent space of the \textsc{Grassmann} manifold. Three study cases are considered to verify the efficiency of the \revision{reduced-order} model. Results highlight that the model has a satisfying accuracy of $10^{\,-3}\,$ compared to reference solutions. The last case study focuses on the wall energy efficiency design under climate change according to a \revision{four-dimensional} parameter space. The latter is composed of the load material emissivity, heat capacity, thermal conductivity, and thickness insulation layer. Simulations are carried over $30$ years considering climate change. The solution minimizing the wall work rate is determined with a computational ratio of $0.1\%$ compared to standard approaches.

The estimation of wall thermal properties through an inverse problem procedure enables to increase the reliability of the model predictions for building energy efficiency. Nevertheless, it requires to define an experimental campaign to obtain on-site observations for existing buildings. The design of experiments enables to search for the optimal measurement plan that ensure the highest precision of the parameter estimation. For on-site measurement in buildings, it seeks an answer to several questions such as number and position of sensors, period and duration of the experimental campaign. In this article, the Optimal Experiment Design (OED) methodology is applied for a multi-layer building wall to resolve the aforementioned issues.

International audience; In this article, both numerical and experimental investigations were carried out on the durability of hemp concrete. For this, an accelerated aging process was performed using cycles of immersion, freezing and drying. Then, an experimental campaign was enabled to determine heat and mass transfer properties, as well as the microstructure for both aged and reference materials. Observations using a digital microscope showed the appearance of cracks at the interfaces and an increase of the porosity of about 6%. These microstructural modifications imply a non-negligible evolution of heat and mass transfer properties. Thus, a numerical model for the prediction of heat and mass transfer was developed. The prediction of physical phenomena was computed using both aged and reference material properties. It highlights the aging effects on the behaviour of the hemp concrete. The numerical simulation results showed significant discrepancies between the predicted relative humidity values for the two configurations (aged and reference) of about 18% and a maximum phase shift of 40 min, due to the amplification of the mass transfer kinetics after aging. Nevertheless, few deviations in temperature values were found. Thus, after aging, sensible heat fluxes were overestimated compared to the reference case, unlike latent heat fluxes, where an underestimation was shown.

Within the framework of building energy assessment, this article proposes to use a derivative based sensitivity analysis of heat transfer models in a building envelope. Two, global and local, estimators are obtained at low computational cost, to evaluate the influence of the parameters on the model outputs. Ranking of these estimators values allows to reduce the number of model unknown parameters by excluding non-significant parameters. A comparison with variance and regression-based methods is carried out and the results highlight the satisfactory accuracy of the continuous-based approach. Moreover, for the carried investigations the approach is $100$ times faster compared to the variance-based methods. A case study applies the method to a real-world building wall. The sensitivity of the thermal loads to local or global variations of the wall thermal is investigated. Additionally, a case study of wall with window is analyzed.

Predictions of physical phenomena in buildings are carried out by using physical models formulated as a mathematical problem and solved by means of numerical methods, aiming at evaluating, for instance, the building thermal or hygrothermal performance by calculating distributions and fluxes of heat and moisture transfer. Therefore, the choice of the numerical method is crucial since it is a compromise among (i) the solution accuracy, (ii) the computational cost to obtain the solution and (iii) the complexity of the method implementation. An efficient numerical method enables to compute an accurate solution with a minimum computational run time (CPU). On that account, this article brings an investigation on the performance of three numerical methods. The first one is the standard and widely used finite-difference approach, while the second one is the so-called RC approach, which is a particular method brought to the building physics area by means of an analogy of electric circuits. The third numerical method is the spectral one, which has been recently proposed to solve nonlinear diffusive problems in building physics. The three methods are evaluated in terms of accuracy on the assessment of the dependent variable (temperature or vapor pressure) or of density of fluxes for three different cases: i) heat diffusion through a concrete slab, ii) moisture diffusion through an aerated concrete slab and iii) heat diffusion using measured temperatures as boundary conditions. Results highlight the spectral approach as the most accurate method. The RC based model with a few number of resistances does not provide accurate results for temperature and vapor pressure distributions neither to flux densities nor conduction loads., 35 pages, 19 figures, 2 tables, 42 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

International audience; This work provides the first comprehensive petrological and multi-scale structural investigation of a suite of ultramafic bodies and enclosing metasediments from the Late Paleozoic accretionary prism of South-Central Chile. The targeted outcrops are located in the La Cabaña area and are characterized by the presence of three main types of ultramafic rocks: (1) partially serpentinized massive peridotites, (2) antigoritic schistose serpentinites, and (3) antigoritic mylonitic serpentinites. Partially serpentinized massive peridotites and schistose serpentinites recorded a common first hydration event with the formation of Fe2+ porous rim in chromites (at ~500 °C), followed by static partial transformation of olivine to lizardite (below ~300 °C) and formation of rims of Fe3+-rich chromite and magnetite around chromites. A second event of serpentinization linked to fluid passage through zones with focused deformation resulted in the partial transformation of partially lizarditized massive peridotites into antigoritic schistose serpentinites (at ~320–400 °C). This second event took place within the shallower portion of the serpentinitic subduction channel, near the base of the accretionary prism and was characterized by the formation of tubular folds made up of lenses of olivine-lizardite massive rocks embedded in a matrix of schistose serpentinites. Antigoritic mylonitic serpentinites encountered in the matrix of schistose serpentinites may locally contain metamorphic olivine (i.e., olivine-bearing mylonitic serpentinites), registering higher temperature conditions at medium pressures (~600 °C, 11 kbar) that prevailed during the initial immature stage of development of the subduction system. Moreover, some mylonitic serpentinites containing Ti-clinohumite, Ti-chondrodite and metamorphic olivine still record higher pressure conditions (>15 kbar, i.e., >50 km depth for ~500 °C), which are related with a deformational event taking place when the subduction system had already cooled. These medium- and high-P mylonitic serpentinites preserved kinematic indicators of burial, whereas their exhumation along the subduction channel towards the shallower schistose serpentinites near the base of the accretionary prism could take place through discrete shear zones possibly aided by large sheath folds. Once incorporated in the accretionary prism (c. 285 Ma ago), the serpentinites shared metamorphism and deformation with the enclosing metasediments. Subvertical shortening related to uplift processes within the accretionary system took place in the stability field of antigorite (>300 °C) and, later, a subhorizontal E-W shortening prevailed in more superficial conditions.

The goal of this study is to propose an efficient numerical model for the predictions of capillary adsorption phenomena in a porous material. The Scharfetter-Gummel numerical scheme is proposed to solve an advection-diffusion equation with gravity flux. Its advantages such as accuracy, relaxed stability condition, and reduced computational cost are discussed along with the study of linear and nonlinear cases. The reliability of the numerical model is evaluated by comparing the numerical predictions with experimental observations of liquid uptake in bricks. A parameter estimation problem is solved to adjust the uncertain coefficients of moisture diffusivity and hydraulic conductivity., Comment: 41 pages, 19 figures, 3 tables, 35 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

International audience; The 1000 km-long suture zone of the Dahomeyide belt, exposed from Southeast Ghana to South Mali, corresponds to a narrow and lithologically diverse area with symptomatic coronitic HP granulitic massifs. Based on a review of published petrological, geochemical and geochronological data along the Dahomeyide belt we propose a global scenario for the closure of the Pharusian ocean between the West African craton (WAC) and the Benino-Nigerian shield during the end of the Neoproterozoic. The onset of a long-lived oceanic subduction by 800–780 Ma is recorded by early magmatism in the Amalaoulaou intra-oceanic arc in Mali, contemporary to the Gourma and Adrar des Iforas to the North, and in Brazil to the South. The first occurrence of tonalitic plutons dated at 720 Ma, and the development of a forearc system around 650 Ma mark the onset of active margin subduction beneath the Benino-Nigerain shield and its northward prolongation in Mali. Oceanic subduction beneath the active margin ended between 640 and 630 Ma with the onset of subduction of the WAC continental margin while subduction related magmatism continues till ca. 600 Ma on the upper plate. During a short period between 620 and 610 Ma, the forearc system and the tip of the active continental margin were buried synchronously. The positive Bouguer anomaly observed to the East of the suture zone, in Benin, and not beneath the suture zone itself supports the occurrence of a massive mafic body at the base of the crust that could correspond to the underplated buried forearc. Ultimately, the exhumation and partial amphibolitization of the suture zone and a shift from a calc-alkaline magmatism to anatectic magmatism along with the onset of strike-slip faulting in the upper plate marks the transition from continental subduction to continental collision between 610 and 580 Ma. This work highlights the importance of the Pharusian suture zone s.l. To our knowledge, it represents a unique example in the world in where the forearc system is buried and partly exhumed at the transition from subduction to collision.

International audience; A suite of mafic pyroxenite xenoliths and clinopyroxene megacrysts was brought to the surface by Cenozoic nephelinites of the Jbel Saghro Volcanic Field (Anti-Atlas, Morocco). The large population of samples was subdivided into five groups: (i) clinopyroxenites sensu stricto; (ii) olivine clinopyroxenites; (iii) mica-bearing clinopyroxenites; (iv) kaersutite-bearing clinopyroxenites; (v) clinopyroxene megacrysts. These xenoliths display a cumulate texture (adcumulate, heteradcumulate with poikilitic clinopyroxene including olivine). The clinopyroxenes have the composition of augite and show an appreciable variation of MgO (7.02–14.80 wt.%), TiO2 (0.58–5.76 wt.%) and Al2O3 (2.81–12.38 wt.%) contents in grains. The clinopyroxenes are characterized by convex upward chondrite-normalized REE patterns, they display very similar trace element compositions with low contents of incompatible elements such as Rb (0−0.9 ppm), Ba (0.1–8.3 ppm), Th (0.1−0.3 ppm), U (0.01−0.04 ppm) and Nb (1.3–3.2 ppm). REE contents of the calculated melts in equilibrium with the clinopyroxene megacrysts and clinopyroxene from pyroxenite xenoliths are similar to those of the nephelinites exposed in Jbel Saghro. Crystallization temperatures of pyroxenite xenoliths and clinopyroxene megacrysts range from 950 °C to 1150 °C. Clinopyroxene barometry yielded pressure of crystallization ranging from 0.4 to 0.8 GPa for pyroxenite xenoliths and 0.3 to 0.7 GPa for clinopyroxene megacrysts. This pressure range is in agreement with pyroxenite xenoliths and clinopyroxene megacrysts being crystallized from their parental melts at the lower and upper crust.

The sorption curve is an essential feature for the modelling of heat and mass transfer in porous building materials. Several models have been proposed in the literature to represent the amount of moisture content in the material according to the water activity (or capillary pressure) level. These models are based on analytical expressions and few parameters that need to be estimated by inverse analysis. This article investigates the reliability of eight models through the accuracy of the estimated parameters. For this, experimental data for a wood fibre material are generated with special attention to the stop criterion to capture long time kinetic constants. Among five sets of measurements, the best estimate is computed. The reliability of the models is then discussed. After proving the theoretical identifiability of the unknown parameters for each model, the primary identifiability is analysed. It evaluates whether the parameters influence on the model output is sufficient to proceed the parameter estimation with accuracy. For this, a continuous derivative-based approach is adopted. Seven models have a low primary identifiability for at least one parameter. Indeed, when estimating the unknown parameters using the experimental observations, the parameters with low primary identifiability exhibit large uncertainties. Finally, an Approximation Bayesian Computation algorithm is used to simultaneously select the best model and estimate the parameters that best represent the experimental data. The thermodynamic and \textsc{Feng}--\textsc{Xing} models, together with a proposed model in this work, were the best ones selected by this algorithm.

Estimating the temperature field of a building envelope could be a time-consuming task. The use of a reduced-order method is then proposed: the Proper Generalized Decomposition method. The solution of the transient heat equation is then re-written as a function of its parameters: the boundary conditions, the initial condition, etc. To avoid a tremendous number of parameters, the initial condition is parameterized. This is usually done by using the Proper Orthogonal Decomposition method to provide an optimal basis. Building this basis requires data and a learning strategy. As an alternative, the use of orthogonal polynomials (Chebyshev, Legendre) is here proposed.

A two-dimensional model is proposed for energy efficiency assessment through the simulation of heat transfer in building envelopes, considering the influence of the surrounding environment. The model is based on the \DF ~approach that provides an explicit scheme with a relaxed stability condition. The model is first validated using an analytical solution and then compared to three other standard schemes. Results show that the proposed model offers a good compromise in terms of high accuracy and reduced computational efforts. Then, a more complex case study is investigated, considering non-uniform shading effects due to the neighboring buildings. In addition, the surface heat transfer coefficient varies with wind velocity and height, which imposes an addition non-uniform boundary condition. After showing the reliability of the model prediction, a comparison over almost $120$ cities in France is carried out between the two- and the one-dimensional approaches of the current building simulation programs. Important discrepancies are observed for regions with high magnitudes of solar radiation and wind velocity. Last, a sensitivity analysis is carried out using a derivative-based approach. It enables to assess the variability of the solution according to the modeling of the two-dimensional boundary conditions. Moreover, the proposed model computes efficiently the solution and its sensitivity to the modeling of the urban environment.

Within the environmental context, several tools based on simulations have been proposed to analyze the physical phenomena of heat and mass transfer in porous materials. However, it is still an open challenge to propose tools that do not require to perform computations to catch the dominant processes. Thus, this article proposes to explore advantages of using a dimensionless analysis by scaling the governing equations of heat and mass transfer. Proposed methodology introduces dimensionless numbers and their nonlinear distortions. The relevant investigation enables to enhance the preponderant phenomena to \emph{(i)} compare different categories of materials, \emph{(ii)} evaluate the competition between heat and mass transfer for each material or \emph{(iii)} describe the transfer in multi-layered wall configurations. It also permits to define hygrothermal kinetic, geometric and dynamic similarities among different physical materials. Equivalent systems can be characterized in the framework of experimental or wall designs. Three cases are presented for similarity studies in terms of \emph{(i)} equivalent material length, \emph{(ii)} time of heat and mass transfer and \emph{(iii)} experimental configurations. All these advantages are illustrated in the given article considering $49$ building materials separated in $7$ categories.

Innovative numerical scheme studied in this work enables to overcome two main limitations of Building Performance Simulation (BPS) programs as high computational cost and the choice of a very fine numerical grid. The method, called Super-Time-Stepping (STS), is novel to the state-of-the-art of building simulations, but has already proved to be sufficiently efficient in recent studies from anisotropic heat conduction in astrophysics (Meyer et al. 2014). The given research is focused on employment of this adopted numerical method to model drying of a rammed earth wall with an additional insulation layer. The results show considerable advantage of the STS method compared to standard Euler explicit scheme. It is possible to choose at least 100 times bigger time-steps to maintain high accuracy and to cut computational cost by more than 92% in the same time.

41 pages, 13 figures, 2 tables, 32 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/; International audience; This article proposes an efficient explicit numerical model with a relaxed stability condition for the simulation of heat, air and moisture transfer in porous material. Three innovative approaches are combined to solve the system of three partial differential equations. The Du Fort-Frankel scheme is used to solve the diffusion equation, providing an explicit scheme with an extended stability region. The two advection--diffusion equations are solved using both Scharfetter-Gummel numerical scheme for the space discretisation and the two-step Runge-Kutta method for the time variable. This combination enables to relax the stability condition by one order. The proposed numerical model is evaluated on three case studies. The first one considers quasi-linear coefficients to confirm the theoretical results by numerical computations. The stability condition is relaxed by a factor of 40 compared to the standard approach. The second case provides an analytical solution for weakly nonlinear problem. A very satisfactory accuracy is observed between the reference solution and the one provided by the numerical model. The last case study assumes more realistic application with nonlinear coefficients and Robin-type boundary conditions. The computational time is reduced 10 times by using the proposed model in comparison with the explicit Euler method.

International audience; The southwestern French Massif central in western Rouergue displays an inverted metamorphic sequence with eclogite and amphibolite facies units forming the top of the nappe stack. They are often grouped into the leptyno-amphibolite complex included, in this area, at the base of the Upper Gneiss Unit. We sampled garnet micaschists and amphibolites to investigate their metamorphic history with isochemical phase diagrams, thermobarometry and U-Pb zircon dating. Our results demonstrate that two different tectono-metamorphic units can be distinguished. The Najac unit consists of biotite-poor phengitegarnet micaschists, a basic-ultrabasic intrusion containing retrogressed eclogites and phengite orthogneisses. Pressure and temperature estimates on micaschists with syn-kinematic garnets yield a prograde with garnet growth starting at 380°C/6-7 kbar, peak pressure at 16 kbar for 570°C, followed by retrogression in the greenschist facies. The age of high pressure metamorphism has been constrained in a recent publication between ca. 383 and 369 Ma. The Laguépie unit comprises garnet-free and garnet-bearing amphibolites with isolated lenses, veins or dykes of leucotonalitic gneiss. Thermobarometry and phase diagram calculation on a garnet amphibolite yield suprasolidus peak P-T conditions at 710°C, 10 kbar followed by retrogression and deformation under greenschist and amphibolite facies conditions. New U-Pb analyses obtained on igneous zircon rims from a leucotonalitic gneiss yield an age of 363 ± 3 Ma, interpreted as the timing of zircon crystallization after incipient partial melting of the host amphibolite. The eclogitic Najac unit records the subduction of a continental margin during Upper Devonian. It is tentatively correlated to a Middle Allochthon, sandwiched between the Lower Gneiss Unit and the Upper Gneiss Unit. Such an intermediate unit is still poorly defined in the French Massif central but it can be a lateral equivalent of the Groix blueschists in the south Armorican massif. The Uppermost Devonian, amphibolite facies Laguépie unit correlates in terms of P-T-t evolution to the Upper Gneiss Unit in the Western French Massif central. This Late Devonian metamorphism is contemporaneous with active margin magmatism and confirms that the French Massif central belonged to the continental upper plate of an ocean-continent subduction system just before the stacking of Mississippian nappes.

36 pages, 23 figures, 3 tables, 44 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/; International audience; The reliability of mathematical models for heat and mass transfer in building porous material is of capital importance. A reliable model permits to carry predictions of the physical phenomenon with sufficient confidence in the results. Among the physical phenomena, the hysteresis effects on moisture sorption and moisture capacity need to be integrated in the mathematical model of transfer. This article proposes to explore the use of an smooth Bang-Bang model to simulate the hysteresis effects coupled with heat and mass transfer in porous material. This model adds two supplementary differential equations to the two classical ones for heat and mass transfer. The solution of these equations ensures smooth transitions between the main sorption and desorption curves. Two parameters are required to control the speed of transition through the intermediary curves. After the mathematical description of the model, an efficient numerical model is proposed to compute the fields with accuracy and reduced computational efforts. It is based on the Du Fort-Frankel scheme for the heat and mass balance equations. For the hysteresis numerical model, an innovative implicit-explicit approach is proposed. Then, the predictions of the numerical model are compared with experimental observations from literature for two case studies. The first one corresponds to a slow cycle of adsorption and desorption while the second is based on a fast cycling case with alternative increase and decrease of moisture content. The comparisons highlight a very satisfactory agreement between the numerical predictions and the observations. In the last Section, the reliability and efficiency of the proposed model is investigated for long term simulation cases. The importance of considering hysteresis effects in the reliability of the predictions are enhanced by comparison with classical approaches from literature.

International audience; This paper presents new petrological, geochemical, isotopic (Nd) and geochronological data on magmatic rocks from the poorly known southern Ouaddaï massif, located at the southern edge of the socalled Saharan metacraton. This area is made of greenschist to amphibolite facies metasediments intruded by large pre-to syn-tectonic batholiths of leucogranites and an association of monzonite, granodiorite and biotite granite forming a late tectonic high-K calc-alkaline suite. U-Pb zircon dating yields ages of 635 ± 3 Ma and 613 ± 8 Ma on a peraluminous biotite-leucogranite (containing numerous inherited Archean and Paleoproterozoic zircon cores) and a muscovite-leucogranite, respectively. Geochemical fingerprints are very similar to some evolved Himalayan leucogranites suggesting their parental magmas were formed after muscovite and biotite dehydration melting of metasedimentary rocks. A biotite-granite sample belonging to the late tectonic high-K to shoshonitic suite contains zircon rims that yield an age of 540 ± 5 Ma with concordant inherited cores crystallized around 1050 Ma. Given the high-Mg# (59) andesitic composition of the intermediate pyroxene-monzonite, the very similar trace-element signature between the different rock types and the unradiogenic isotopic signature for Nd, the late-kinematic high-K to shoshonitic rocks formed after melting of the enriched mantle and further differentiation in the crust. These data indicate that the southern Ouaddaï was part of the Pan-African belt. It is proposed that it represents a continental back-arc basin characterized by a high-geothermal gradient during Early Ediacaran leading to anatexis of middle to lower crustal levels. After tectonic inversion during the main Pan-African phase, late kinematic high-K to shoshonitic plutons emplaced during the final post-collisional stage.

Models for the prediction of heat and mass transfers in building porous materials have been developed since the 90’s with simulation programs such as MATCH, UMIDUS, DELPHIN and Wufi. These models are used to analyze the physical phenomena occurring and particularly the impact of moisture on buildings’ energy performance and durability. With this goal in mind, it is important to validate the representation of the physical phenomena made by the numerical models. This article reviews recent studies comparing the results obtained with numerical models and experimental data. An overall trend can be observed for most of these studies, highlighting that the experimental front always rushes faster than the simulation results. Therefore, this study analyses theses comparisons to explain these discrepancies based on the effects of (i) the type of materials, (ii) the boundary conditions, (iii) the scale of the design facility, (iv) the model used to describe the physical phenomena and (v) the influence of th...

Prediction of moisture transfer within material using a classic diffusive model may lack accuracy, since numerical simulations underestimate the adsorption process when a sample is submitted to variations of moisture level. Model equations are always established with assumptions. Consequently, some phenomena are neglected. This paper therefore investigates the impact of improving traditional diffusive models by taking into account additional phenomena that could occur in moisture transport within hygroscopic fibrous materials such as wood-based products. Two phenomena in the porous material are investigated: (1) non-equilibrium behaviour between water vapour and bound water, and (2) transport by air convection. The equations of each model are established by starting from averaging conservation equations for the different species considered within material (water vapour, bound water and air). In addition, the validity of assumptions currently used in the models is verified. Then the three models are compared with experimental data to highlight their capacity to predict both the vapour pressure and the mass of adsorbed water. This comparison tends to show a slight improvement in predictions with the new models. To increase our understanding of these models, the influence of the main parameters characterising phenomena (sorption coefficient, intrinsic permeability, Peclet number and Fourier number) is studied using local sensitivity analysis. The shape of the sensitivity coefficients shows that the first kinetics period is only impacted a little by the non-equilibrium. In other periods, the influence of the diffusion phenomenon represented by the Fourier number is much greater than that of the two other phenomena: advection and sorption. Nevertheless, the sensitivity study shows that these two phenomena have some influence on vapour pressure.

Located in the Imiter Inlier (Eastern Saghro, Anti-Atlas, Morocco), Ediacaran volcanic dykes have been studied for their petrofabric using Anisotropy of Magnetic Susceptibility (AMS) technique. Four dykes, namely TF, TD, FF and FE show andesitic compositions and are considered to belong to the same dyke swarm. They are oriented respectively N25E, N40E, N50E, and N10E and have been emplaced during a first tectonic event. The dyke FW, oriented N90E displays a composition of alkali basalt and its emplacement is attributed to a subsequent tectonic event. These rocks are propylitized under greenschist facies conditions forming a secondary paragenesis constituted by calcite, chlorite, epidote and sericite. The dykes TF, TD, FF and FE are sub-volcanic calc-alkaline, typical of post-collisional basalts/andesites, belonging to plate margin andesites. The FW dyke shows a within-plate basalt signature; alkaline affinity reflecting a different petrogenetic process. The thermomagnetic analyses show a dominantly ferromagnetic behaviour in the TF dyke core carried by single domain Ti-poor magnetite, maghemite and pyrrhotite. The dominantly paramagnetic susceptibilities in TF dyke rims and TD, FE, FF and FW dykes are controlled by ilmenite, amphibole, pyroxene and chlorite. The magnetic fabrics of the Imider dykes, determined by our AMS study, allows us to reconstitute the tectonic conditions which prevailed during the emplacement of these two generations of volcanic dykes. The first tectonic event was characterized by a roughly NE-SW compression and the second tectonic event is characterized by an E-W shortening followed by a relaxation recording the end of the Pan-African orogeny in the eastern Anti-Atlas.

The standard methods to determine the vapour permeability and the moisture sorption curve may lack of accuracy since discrepancies are observed when comparing numerical predictions to experimental data. Moreover, these properties are determined in steady state conditions while the numerical predictions are carried in transient regime. Thus, this paper presents an experimental design to estimate these properties using dynamic measurements and identification method. The experimental facility is presented, enabling to measure at the same time the relative humidity within the material and the total moisture content. The performance of the facility and protocol in terms of reproducibility, uncertainty and direction of heat and moisture transfers are checked, confirming the abilities of the set-up. Then, experimental results are used to determine the hygrothermal material properties using a trust-region algorithm. Investigations are done to analyse important issues as the choice of the observation: relative humidity and/or mass measurements, to solve the parameter estimation problem. The estimated properties are finally validated by comparing the numerical predictions with experimental data for other boundary conditions.

One possibility to improve the accuracy of building performance simulation (BPS) tools is via co-simulation techniques, where more accurate mathematical models representing particular and complex p...

A computational procedure known as co-simulation has been proposed in the literature as a possibility to extend the capabilities and improve the accuracy of building performance simulation (BPS) to...

We present the first combined U-Pb ages and Lu-Hf isotope data of zircon grains from the southern Ouaddai area in Eastern Chad in the Central Africa Orogenic Belt. Zircon xenocrysts from granites and detrital zircons from quartzites of the Ouaddai document crustal accretion from depleted mantle sources at 3.0 ± 0.5 Ga, rarely up to 3.8 Ga, and for subsequent reworking during four magmatic events, during the Neoarchean at ca. 2.6 Ga (eHft = −7.0 to 4.6), Paleoproterozoic at ca. 1.9–2.1 Ga (eHft = −17.7 to 5.3), Mesoproterozoic at ca. 1.1 Ga (eHft = −15.3 to 5.9), and Neoproterozoic/early Cambrian at 635–540 Ma (eHft = −31.1 to 1.4). Superchondritic eHft values point to new crust addition during all four events. Metamorphic zircon rims yield a Concordia age of 600 ± 2 Ma, and provide evidence that greenschist to amphibolite-facies metamorphism affecting the quartzites occurred between syn-tectonic (635 Ma) and post-kinematic (540 Ma) granite intrusions. Our new set of zircon data, combined with field relationships, also indicate that the quartzite precursor rocks were deposited between 1045 and 635 Ma, possibly after 670 Ma in a basin that subsequently became the Central Africa Orogenic Belt. The detritus was supplied from surrounding cratons, although the sources of the Mesoproterozoic zircons remain ambiguous. Comparison of age-eHft data from worldwide sources suggests sediment supply from the Irumide, Katanga and Mozambique belts, located to the south-east of the Congo Craton. Paleogeographic reconstructions further require sediment transport prior to final closure of several oceanic basins surrounding Central and Northern Africa during the Neoproterozoic, comprising the Khomas and Mozambique oceans

Implicit schemes require important sub-iterations when dealing with highly nonlinear problems such as the combined heat and moisture transfer through porous building elements. The computational cost rises significantly when the whole-building is simulated, especially when there is important coupling among the building elements themselves with neighbouring zones and with HVAC (Heating Ventilation and Air Conditioning) systems. On the other hand, the classical Euler explicit scheme is generally not used because its stability condition imposes very fine time discretisation. Hence, this paper explores the use of an improved explicit approach - the Dufort-Frankel scheme - to overcome the disadvantage of the classical explicit one and to bring benefits that cannot be obtained by implicit methods. The Dufort-Frankel approach is first compared to the classical Euler implicit and explicit schemes to compute the solution of nonlinear heat and moisture transfer through porous materials. Then, the analysis of the Dufort-Frankel unconditionally stable explicit scheme is extended to the coupled heat and moisture balances on the scale of a one- and a two-zone building models. The Dufort-Frankel scheme has the benefits of being unconditionally stable, second-order accurate in time O(dt^2) and to compute explicitly the solution at each time step, avoiding costly sub-iterations. This approach may reduce the computational cost by twenty, as well as it may enable perfect synchronism for whole-building simulation and co-simulation., 38 pages, 19 figures, 3 tables, 28 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

This work is devoted to proposing a hybrid numerical–analytical method to address the problem of heat and moisture transfer in porous soils. Several numerical and analytical models have been used to study heat and moisture transfer. The complexity of the coupled transfer in soils is such that analytical solutions exist only for limited problems, while numerical solutions can deal with more realistic ones but at a higher computational cost. Therefore, we propose to implement analytical solutions where variations of temperature and moisture content are known to be almost nonvarying, while the numerical solution is implemented in the remaining region, near the boundaries. The coupling between solutions is performed assuming the continuity of both fields and fluxes at each interface. This strategy allows assuring the physical phenomenon occurring at the interface. Numerical experiments are performed, showing the accuracy, the efficiency, and the great potential of the method regarding applications in nonlinear soil problems.

The pre-Alpine basement of the Southern Carpathians/Western Balkans contains four ophiolitic massifs dismembered by Alpine tectonics, which define the “Balkan-Carpathian Ophiolite” (BCO) for which the tectonic setting and age of formation are still debated (Precambrian or Early Devonian). In this contribution, we demonstrate that, in light of a Pre-Alpine restoration, the four massifs belonged to a unique slice of very complete, obducted oceanic lithosphere and we re-evaluate its tectonic setting. Large chromitite volumes with Al-rich spinel compositions (Cr# = 0.39–0.48), as well as major and trace geochemical results on basalts (slightly enriched N-MORBs with low negative Nb anomaly associated with calk-alkaline BABBs), point to a formation in a back-arc basin. Mantle spinel composition (Cr# = 0.49–0.51) and melting modeling indicate mean melting extents of 8.5–11% favouring intermediate spreading rate. New Sm-Nd dating on lower gabbroic rocks give a whole rock isochron, interpreted as the age of formation of the BCO crust at 409 ± 38 Ma, thus confirming an Early Devonian oceanic crust. The previous ∼563 Ma U-Pb zircon age can be interpreted as casual inheritance indicating the proximity of an old continental lithosphere. Taking into account the lithological evidences and paleocontinental affinities of the two recognized terranes separated by the BC oceanic basin (Balkans and Sredna Gora) and by analogy with other Variscan ophiolites in Western/Central Europe, we suggest that the BC ophiolite belong to the ∼400 Ma ophiolites group obducted between West and East Galatia and belonging to the southern Variscan suture. However, the BC ophiolite is the only one of this group obducted to the north and not involved in the Lower Allochthon/ophiolite/Upper Allochthon thrust pile, likely explaining its exceptional preservation. Finally, we tentatively propose a new unifying tectonic model where different terrane drift rates and highly oblique displacements create two Rheic branches, the “Rheic” and the “Galicia-Brittany-Massif Central”.