An efficient method to estimate sorption isotherm curve coefficients

DOI10.1080/17415977.2018.1495720zbMATH Open1483.65148arXiv1704.06565OpenAlexW3099235953WikidataQ129492403 ScholiaQ129492403MaRDI QIDQ4990724FDOQ4990724

Authors: Julien Berger, Thomas Busser, Denys Dutykh, N. Mendes 

Publication date: 31 May 2021

Published in: Inverse Problems in Science and Engineering (Search for Journal in Brave)

Abstract: 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.

Full work available at URL: https://arxiv.org/abs/1704.06565

Recommendations

• Numerical determination of the sorption isotherm for nonisothermal sorption
• Numerical methods for solving some inverse problems in sorbtion dynamics
• Mathematical modeling numerical determination of coefficients in some mathematical models of nonisothermal sorption dynamics
• Numerical study for the parameter estimation of the moisture transfer coefficient : 2D case
• Numerical solution for the parameter estimation of the moisture transfer coefficient

zbMATH Keywords

hysteresisparameter estimation problemsensitivity functionsconvective moisture transferoptimal experiment design (OED)sorption moisture coefficients

Mathematics Subject Classification ID

Numerical optimization and variational techniques (65K10) Sensitivity (robustness) (93B35) Ill-posed problems for PDEs (35R25) Inverse problems for PDEs (35R30) Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations (65L06) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Numerical methods for ill-posed problems for initial value and initial-boundary value problems involving PDEs (65M30) Numerical methods for inverse problems for initial value and initial-boundary value problems involving PDEs (65M32) System identification (93B30) PDEs in connection with classical thermodynamics and heat transfer (35Q79) Inverse problems in thermodynamics and heat transfer (80A23)

Cites Work

• The MATLAB ODE Suite
• Optimal Measurement Methods for Distributed Parameter System Identification
• A trust region method based on interior point techniques for nonlinear programming.
• Definitions and examples of inverse and ill-posed problems
• A gradient descent method for solving an inverse coefficient heat conduction problem
• Inverse and ill-posed problems. Theory and applications.
• Extreme methods for solving ill-posed problems with applications to inverse heat transfer problems. Transl. from the Russian
• Title not available (Why is that?)
• Qualitative and quantitative experiment design for phenomenological models - a survey
• Heat and mass transfer in capillary-porous bodies. Transl. from the Russian
• Conservative modelling of the moisture and heat transfer in building components under atmospheric excitation
• A method for predicting heat and moisture transfer through multilayered walls based on temperature and moisture content gradients
• Global approaches to identifiability testing for linear and nonlinear state space models
• Multitridiagonal-Matrix Algorithm for Coupled Heat Transfer in Porous Media: Stability Analysis and Computational Performance
• Analysis of parameter sensitivity and experimental design for a class of nonlinear partial differential equations
• Estimation of thermophysical properties of moist materials under different drying conditions

Cited In (3)

• Numerical study for the parameter estimation of the moisture transfer coefficient : 2D case
• Numerical solution for the parameter estimation of the moisture transfer coefficient
• Modelling the water sorption isotherms of yoghurt powder spray

Uses Software

• Matlab
• ode23
• ode45
• Ode15s
• ode113
• ode23s
• MATLAB ODE suite
• IBPT
• WUFI

This page was built for publication: An efficient method to estimate sorption isotherm curve coefficients

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4990724)