[論文レビュー] On the estimation of sorption isotherm coefficients using the optimal experiment design approach

Julien Berger, Thomas Busser+2 more|arXiv (Cornell University)|Apr 15, 2017

Hygrothermal properties of building materials被引用数 1

ひとこと要約

本研究は、モデルとデータの乖離を最小限に抑えることで、木材繊維材料の吸着等温線係数3つを正確に推定するための最適実験設計（OED）手法を提案する。感度に基づくOEDと、2種の異なる湿度ステップ実験にそれぞれ別個の目的関数を用いた内点法を用いることで、パラメータの不確実性を低減し、特に水分吸着および脱着挙動の再現性を従来の研究よりも向上させる。

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 in 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 experiments. 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 allow to capture an approximation of the probability distribution function of the estimated parameters. The determined sorption isotherm coefficients enable to calibrate the numerical model and fit better the experimental data, reducing the discrepancies usually reported in the literature that underestimate the moisture adsorption and overestimate the desorption processes.

研究の動機と目的

• 建築材料、特に木材繊維における不正確な吸着等温線パラメータ推定の課題に対処すること。
• 文献において水分吸着が低く見積もられ、脱着が過大に見積もられるという乖離を低減すること。
• 感度解析を用いて、パラメータ同定の精度を最大化する最適な実験設定を開発すること。
• 最適なセンサ配置および境界条件の同定を通じて、湿気移動モデルのより正確なキャリブレーションを実現すること。
• 吸着等温線パラメータの不確実性を最小限に抑えることで、建築物理学における数値シミュレーションの信頼性を向上させること。

提案手法

• 本研究では、木材繊維内の一時的湿気拡散をシミュレートするため、対流的湿気輸送に基づく数学的モデルを採用する。
• 感度関数と事前パラメータ値を用いて、20通りの実験設計から最良の実験構成を選択する最適実験設計（OED）を適用する。
• 2つの実験手順を検証した：単一ステップの湿度変化（10% から 75%）と、8日間隔の多段階サイクル（10-75-33-75%）。
• パラメータ推定には、L₂ノルムに基づく目的関数を最小化する内点法を用いる。
• 目的関数は各実験ごとに別々に評価され、実験データとモデル予測の残差を最小化するために離散的ユークリッドノルムが使用される。
• OEDから得られた感度関数を用いて、推定パラメータの確率分布の近似が可能となり、不確実性の定量化が強化される。

実験結果

リサーチクエスチョン

• RQ1 単一ステップの湿度変化と多段階湿度サイクリングのどちらの実験設計が、吸着等温線係数の推定において最も正確か？
• RQ2 湿気移動実験において、パラメータ同定の精度を最大化するためにセンサをどこに配置すべきか？
• RQ3 目的関数ノルム（L₂ 対 L∞）および実験のグループ化の選択が、パラメータ推定性能にどのように影響するか？
• RQ4 OED手法は、従来の実験設計と比較して、推定パラメータの不確実性を低減し、モデル適合性を向上させることができるか？
• RQ5 提案手法は、木材繊維材料において一般的に報告される、吸着を低く見積もって脱着を高く見積もるという文献バイアスをどの程度低減できるか？

主な発見

• 最適実験設計により、不透湿性境界に近い位置にセンサを配置することで、感度とパラメータ精度を最大にできることが同定された。
• 単一ステップ（10% から 75%）および多段階（10-75-33-75%）の湿度実験の両方が有効であったが、各実験に対して別個の目的関数を用いることでより優れた結果が得られた。
• 離散的ユークリッドノルム（L₂）を用い、実験を別々に処理することで、モデル予測と実験データの残差が最小化された。
• 内点法は100回未満の直接モデル評価で収束したため、計算効率が優れていた。
• OEDの感度関数により、推定パラメータの確率分布の近似が可能となり、不確実性の定量化が支援された。
• キャリブレーション済みモデルは、水分吸着および脱着をよりよく捉えており、従来の研究で報告された系統的乖離が低減された。