Julien Berger

Suelen Gasparin, Julien Berger, Denys Dutykh et al.

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 τ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(Δ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 $Δt$ remains an important issue to accurately represent the physical phenomena.

Julien Berger, Hervé Le Meur, Denys Dutykh et al.

The reliability of a model is its accuracy in predicting the physical phenomena using the known input parameters. It also depends on the model's ability to estimate relevant parameters using observations of the physical phenomena. In this paper, the reliability of the VTT model is investigated under these two criteria for various given temperature and relative humidity constant in time. First of all, experiments are conducted on bamboo fiberboard. Using these data, five parameters of the VTT model, defining the mold vulnerability class of a material, are identified. The results highlight that the determined parameters are not within the range of the classes defined in the VTT model. In addition, the quality of the parameter estimation is not satisfactory. Then the sensitivity of the numerical results of the VTT model is analyzed by varying an input parameter. These investigations show that the VTT mathematical formulation of the physical model of mold growth is not reliable. An improved model is proposed with a new mathematical formulation. It is inspired by the logistic equation whose parameters are estimated using the experimental data obtained. The parameter estimation is very satisfactory. In the last parts of the paper, the numerical predictions of the improved model are compared to experimental data from the literature to prove its reliability.

Julien Berger, Suelen Gasparin, Denys Dutykh et al.

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.

Julien Berger, Denys Dutykh, Nathan Mendes

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.

Julien Berger, Denys Dutykh, Nathan Mendes et al.

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.

Julien Berger, Thomas Busser, Denys Dutykh et al.

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.

Julien Berger, Denys Dutykh

The fidelity of a model relies both on its accuracy to predict the physical phenomena and its capability to estimate unknown parameters using observations. This article focuses on this second aspect by analyzing the reliability of two mathematical models proposed in the literature for the simulation of heat losses through building walls. The first one, named DuFort-Frankel (DF), is the classical heat diffusion equation combined with the DuFort-Frankel numerical scheme. The second is the so-called RC lumped approach, based on a simple ordinary differential equation to compute the temperature within the wall. The reliability is evaluated following a two stages method. First, samples of observations are generated using a pseudo-spectral numerical model for the heat diffusion equation with known input parameters. The results are then modified by adding a noise to simulate experimental measurements. Then, for each sample of observation, the parameter estimation problem is solved using one of the two mathematical models. The reliability is assessed based on the accuracy of the approach to recover the unknown parameter. Three case studies are considered for the estimation of (i) the heat capacity, (ii) the thermal conductivity or (iii) the heat transfer coefficient at the interface between the wall and the ambient air. For all cases, the DF mathematical model has a very satisfactory reliability to estimate the unknown parameters without any bias. However, the RC model lacks of fidelity and reliability. The error on the estimated parameter can reach 40% for the heat capacity, 80% for the thermal conductivity and 450% for the heat transfer coefficient.

Ainagul Jumabekova, Julien Berger, Denys Dutykh et al.

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.