Search

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.

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., Comment: 35 pages, 19 figures, 2 tables, 42 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

This work presents a detailed mathematical model combined with an innovative efficient numerical model to predict heat, air and moisture transfer through porous building materials. The model considers the transient effects of air transport and its impact on the heat and moisture transfer. The achievement of the mathematical model is detailed in the continuity of Luikov's work. A system composed of two advection-diffusion differential equations plus one exclusively diffusion equation is derived. The main issue to take into account the transient air transfer arises in the very small characteristic time of the transfer, implying very fine discretisation. To circumvent these difficulties, the numerical model is based on the Du Fort-Frankel explicit and unconditionally stable scheme for the exclusively diffusion equation. It is combined with a two-step Runge-Kutta scheme in time with the Scharfetter-Gummel numerical scheme in space for the coupled advection-diffusion equations. At the end, the numerical model enables to relax the stability condition, and, therefore, to save important computational efforts. A validation case is considered to evaluate the efficiency of the model for a nonlinear problem. Results highlight a very accurate solution computed about 16 times faster than standard approaches. After this numerical validation, the reliability of the mathematical model is evaluated by comparing the numerical predictions to experimental observations. The latter is measured within a multi-layered wall submitted to a sudden increase of vapor pressure on the inner side and driven climate boundary conditions on the outer side. A very satisfactory agreement is noted between the numerical predictions and experimental observations indicating an overall good reliability of the proposed model., Comment: 46 pages, 12 figures, 3 tables, 46 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

It is well known that thermal insulation is a leading strategy for reducing energy consumption associated to heating or cooling processes in buildings. Nevertheless, building insulation can generate high expenditures so that the selection of an optimum insulation thickness requires a detailed energy simulation as well as an economic analysis. In this way, the present study proposes an innovative non-uniform adaptive method to determine the optimal insulation thickness of external walls. First, the method is compared with a reference solution to properly understand the features of the method, which can provide high accuracy with less spatial nodes. Then, the adaptive method is used to simulate the transient heat conduction through the building envelope of buildings located in Brazil, where there is a large potential of energy reduction. Simulations have been efficiently carried out for different wall and roof configurations, showing that the innovative method efficiently provides a gain of 25% on the computer run time., Comment: 36 pages, 18 figures, 6 tables, 42 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

This work presents an efficient numerical method based on spectral expansions for simulation of heat and moisture diffusive transfers through multilayered porous materials. Traditionally, by using the finite-difference approach, the problem is discretized in time and space domains (Method of lines) to obtain a large system of coupled Ordinary Differential Equations (ODEs), which is computationally expensive. To avoid such a cost, this paper proposes a reduced-order method that is faster and accurate, using a much smaller system of ODEs. To demonstrate the benefits of this approach, tree case studies are presented. The first one considers nonlinear heat and moisture transfer through one material layer. The second case - highly nonlinear - imposes a high moisture content gradient - simulating a rain like condition - over a two-layered domain, while the last one compares the numerical prediction against experimental data for validation purposes. Results show how the nonlinearities and the interface between materials are easily and naturally treated with the spectral reduced-order method. Concerning the reliability part, predictions show a good agreement with experimental results, which confirm robustness, calculation efficiency and high accuracy of the proposed approach for predicting the coupled heat and moisture transfer through porous materials., Comment: 44 pages, 17 figures, 6 tables, 50 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

This work presents an alternative view on the numerical simulation of diffusion processes applied to the heat and moisture transfer through porous building materials. Traditionally, by using the finite-difference approach, the discretization follows the Method Of Lines (MOL), when the problem is first discretized in space to obtain a large system of coupled Ordinary Differential Equations (ODEs). Thus, this paper proposes to change this viewpoint. First, we discretize in time to obtain a small system of coupled ODEs, which means instead of having a Cauchy (Initial Value) Problem (IVP), we have a Boundary Value Problem (BVP). Fortunately, BVPs can be solved efficiently today using adaptive collocation methods of high order. To demonstrate the benefits of this new approach, three case studies are presented, in which one of them is compared with experimental data. The first one considers nonlinear heat and moisture transfer through one material layer while the second one considers two material layers. Results show how the nonlinearities and the interface between materials are easily treated, by reasonably using a fourth-order adaptive method. Finally, the last case study compares numerical results with experimental measurements, showing a good agreement., Comment: 45 pages, 21 figures, 5 tables, 47 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

Comparisons of experimental observation of heat and moisture transfer through porous building materials with numerical results have been presented in numerous studies reported in literature. However, some discrepancies have been observed, highlighting underestimation of sorption process and overestimation of desorption process. Some studies intend to explain the discrepancies by analysing the importance of hysteresis effects as well as carrying out sensitivity analyses on the input parameters as convective transfer coefficients. This article intends to investigate the accuracy and efficiency of the coupled solution by adding advective transfer of both heat and moisture in the physical model. In addition, the efficient Scharfetter and Gummel numerical scheme is proposed to solve the system of advection-diffusion equations, which has the advantages of being well-balanced and asymptotically preserving. Moreover, the scheme is particularly efficient in terms of accuracy and reduction of computational time when using large spatial discretisation parameters. Several linear and non-linear cases are studied to validate the method and highlight its specific features. At the end, an experimental benchmark from the literature is considered. The numerical results are compared to the experimental data for a pure diffusive model and also for the proposed model. The latter presents better agreement with the experimental data. The influence of the hysteresis effects on the moisture capacity is also studied, by adding a third differential equation., Comment: 51 pages, 21 figures, 3 tables, 44 references. Other author's publications can be downloaded at http://www.denys-dutykh.com/

It is of great concern to produce numerically efficient methods for moisture diffusion through porous media, capable of accurately calculate moisture distribution with a reduced computational effort. In this way, model reduction methods are promising approaches to bring a solution to this issue since they do not degrade the physical model and provide a significant reduction of computational cost. Therefore, this article explores in details the capabilities of two model-reduction techniques - the Spectral Reduced-Order Model (Spectral-ROM) and the Proper Generalised Decomposition (PGD) - to numerically solve moisture diffusive transfer through porous materials. Both approaches are applied to three different problems to provide clear examples of the construction and use of these reduced-order models. The methodology of both approaches is explained extensively so that the article can be used as a numerical benchmark by anyone interested in building a reduced-order model for diffusion problems in porous materials. Linear and non-linear unsteady behaviors of unidimensional moisture diffusion are investigated. The last case focuses on solving a parametric problem in which the solution depends on space, time and the diffusivity properties. Results have highlighted that both methods provide accurate solutions and enable to reduce significantly the order of the model around ten times lower than the large original model. It also allows an efficient computation of the physical phenomena with an error lower than 10^{-2} when compared to a reference solution., Comment: 42 pages, 14 figures, 1 table, 69 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/. arXiv admin note: text overlap with arXiv:1704.07607

This paper proposes the use of a Spectral method to simulate diffusive moisture transfer through porous materials as a Reduced-Order Model (ROM). The Spectral approach is an a priori method assuming a separated representation of the solution. The method is compared with both classical Euler implicit and Crank-Nicolson schemes, considered as large original models. Their performance - in terms of accuracy, complexity reduction and CPU time reduction - are discussed for linear and nonlinear cases of moisture diffusive transfer through single and multi-layered one-dimensional domains, considering highly moisture-dependent properties. Results show that the Spectral reduced-order model approach enables to simulate accurately the field of interest. Furthermore, numerical gains become particularly interesting for nonlinear cases since the proposed method can drastically reduce the computer run time, by a factor of 100, when compared to the traditional Crank-Nicolson scheme for one-dimensional applications., Comment: 42 pages, 22 figures, 2 tables, 51 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

This paper deals with an inverse problem applied to the field of building physics to experimentally estimate three sorption isotherm coefficients of a wood fiber material. First, the mathematical model, based on convective transport of moisture, the Optimal Experiment Design (OED) and the experimental set-up are presented. Then measurements of relative humidity within the material are carried out, after searching the OED, which is based on the computation of the sensitivity functions and a priori values of the unknown parameters employed in the mathematical model. The OED enables to plan the experimental conditions in terms of sensor positioning and boundary conditions out of 20 possible designs, ensuring the best accuracy for the identification method and, thus, for the estimated parameter. Two experimental procedures were identified: i) single step of relative humidity from 10% to 75% and ii) multiple steps of relative humidity 10-75-33-75% with an 8-day duration period for each step. For both experiment designs, it has been shown that the sensor has to be placed near the impermeable boundary. After the measurements, the parameter estimation problem is solved using an interior point algorithm to minimize the cost function. Several tests are performed for the definition of the cost function, by using the L^2 or L^\infty norm and considering the experiments separately or at the same time. It has been found out that the residual between the experimental data and the numerical model is minimized when considering the discrete Euclidean norm and both experiments separately. It means that two parameters are estimated using one experiment while the third parameter is determined with the other experiment. Two cost functions are defined and minimized for this approach. Moreover, the algorithm requires less than 100 computations of the direct model to obtain the solution. In addition, the OED sensitivity functions enable to capture an approximation of the probability distribution function of the estimated parameters. The determined sorption isotherm coefficients calibrate the numerical model to fit better the experimental data. However, some discrepancies still appear since the model does not take into account the hysteresis effects on the sorption capacity. Therefore, the model is improved proposing a second differential equation for the sorption capacity to take into account the hysteresis between the main adsorption and desorption curves. The OED approach is also illustrated for the estimation of five of the coefficients involved in the hysteresis model. To conclude, the prediction of the model with hysteresis are compared with the experimental observations to illustrate the improvement of the prediction., Comment: 51 pages, 20 figures, 4 tables and 56 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

This paper presents a practical application of the concept of Optimal Experiment Design (OED) for the determination of properties of porous materials with in-situ measurements and an identification method. First, an experimental set-up was presented and used for the measurement of relative humidity within a wood fibre material submitted to single and multiple steps of relative humidity variation. Then, the application of OED enabled to plan the experimental conditions in terms of sensor positioning and boundary conditions out of 20 possible designs. The OED search was performed using the Fisher information matrix and a priori knowledge of the parameters. It ensures to provide the best accuracy of the identification method and thus the estimated parameter. Optimal design results have been found for single steps from the relative humidity phi = 10 to 75%, with one sensor located at the position X between 4 and 6 cm, for the estimation of moisture permeability coefficients, while from phi = 75% to phi = 33%, with one sensor located at X{\deg} = 3 cm, for the estimation of the advection coefficient. The OED has also been applied for the identification of couples of parameters. A sample submitted to multiple relative humidity steps (phi = 10-75-33-75%) with a sensor placed at X{\deg} = 5 cm was found as the best option for determining both properties with the same experiment. These OED parameters have then been used for the determination of moisture permeability and advection coefficients. The estimated moisture permeability coefficients are twice higher than the a priori values obtained using standard methods. The advection parameter corresponds to the mass average velocity of the order of v = 0.01 mm/s within the material and may play an important role on the simulation of moisture front., Comment: 34 pages, 14 figures, 6 tables, 43 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/. arXiv admin note: text overlap with arXiv:1704.06565

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., Comment: 38 pages, 19 figures, 3 tables, 28 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

Implicit schemes have been extensively used in building physics to compute the solution of moisture diffusion problems in porous materials for improving stability conditions. Nevertheless, these schemes require important sub-iterations when treating non-linear problems. To overcome this disadvantage, this paper explores the use of improved explicit schemes, such as Dufort-Frankel, Crank-Nicolson and hyperbolisation approaches. A first case study has been considered with the hypothesis of linear transfer. The Dufort-Frankel, Crank-Nicolson and hyperbolisation schemes were compared to the classical Euler explicit scheme and to a reference solution. Results have shown that the hyperbolisation scheme has a stability condition higher than the standard Courant-Friedrichs-Lewy (CFL) condition. The error of this schemes depends on the parameter \tau representing the hyperbolicity magnitude added into the equation. The Dufort-Frankel scheme has the advantages of being unconditionally stable and is preferable for non-linear transfer, which is the second case study. Results have shown the error is proportional to O(\Delta t). A modified Crank-Nicolson scheme has been proposed in order to avoid sub-iterations to treat the non-linearities at each time step. The main advantages of the Dufort-Frankel scheme are (i) to be twice faster than the Crank-Nicolson approach; (ii) to compute explicitly the solution at each time step; (iii) to be unconditionally stable and (iv) easier to parallelise on high-performance computer systems. Although the approach is unconditionally stable, the choice of the time discretisation $\Delta t$ remains an important issue to accurately represent the physical phenomena., Comment: 35 pages, 16 figures, 1 table, 32 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/. arXiv admin note: text overlap with arXiv:1612.07649

When comparing measurements to numerical simulations of moisture transfer through porous materials a rush of the experimental moisture front is commonly observed in several works shown in the literature, with transient models that consider only the diffusion process. Thus, to overcome the discrepancies between the experimental and the numerical models, this paper proposes to include the moisture advection transfer in the governing equation. To solve the advection-diffusion differential equation, it is first proposed two efficient numerical schemes and their efficiencies are investigated for both linear and nonlinear cases. The first scheme, Scharfetter-Gummel (SG), presents a Courant-Friedrichs-Lewy (CFL) condition but is more accurate and faster than the second scheme, the well-known Crank-Nicolson approach. Furthermore, the SG scheme has the advantages of being well-balanced and asymptotically preserved. Then, to conclude, results of the convective moisture transfer problem obtained with the SG numerical scheme are compared to experimental data from the literature. The inclusion of an advective term in the model may clearly lead to better results than purely diffusive models., Comment: 34 pages, 15 figures, 48 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

In the context of estimating material properties of porous walls based on in-site measurements and identification method, this paper presents the concept of Optimal Experiment Design (OED). It aims at searching the best experimental conditions in terms of quantity and position of sensors and boundary conditions imposed to the material. These optimal conditions ensure to provide the maximum accuracy of the identification method and thus the estimated parameters. The search of the OED is done by using the Fisher information matrix and a priori knowledge of the parameters. The methodology is applied for two case studies. The first one deals with purely conductive heat transfer. The concept of optimal experiment design is detailed and verified with 100 inverse problems for different experiment designs. The second case study combines a strong coupling between heat and moisture transfer through a porous building material. The methodology presented is based on a scientific formalism for efficient planning of experimental work that can be extended to the optimal design of experiments related to other problems in thermal and fluid sciences., Comment: 32 pages, 14 figures, 3 tables, 34 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/