Filipe Fonseca Saboia, Cauane Soares Pereira, Joao Gabriel Pinheiro Sousa, Samuel Jorge Marques Cartaxo
{"title":"Inverse physics-informed neural network for CaCO3 crystallization fouling: Integrating initial rate and asymptotic resistance across multi-source data","authors":"Filipe Fonseca Saboia, Cauane Soares Pereira, Joao Gabriel Pinheiro Sousa, Samuel Jorge Marques Cartaxo","doi":"10.1016/j.cherd.2026.09.041","DOIUrl":null,"url":null,"abstract":"<div><div>Crystallization fouling of CaCO<sub>3</sub> in cooling water systems degrades heat transfer efficiency, yet the transient initial fouling rate and the asymptotic fouling resistance are traditionally evaluated using separate empirical correlations developed on distinct experimental rigs. Reconciling both descriptors within a single data-driven formulation is challenging because experimental datasets are scarce, heterogeneous across flow geometries, and sensitive to reporting conventions. This study proposes a compact Inverse Physics-Informed Neural Network (<span>CaCO<sub>3</sub>-PINN</span>) with two outputs and 9703 trainable parameters, informed by dimensionless supersaturation, temperature, and wall shear stress. The network is trained on 25 experimental points compiled from three independent sources (Watkinson 1983; Choi, via Liu 1999; Arsenyeva 2013) with a physics-constrained loss coupling rate and resistance via the Arsenyeva relationship <span><math><mrow><msubsup><mrow><mi>R</mi></mrow><mrow><mi>f</mi></mrow><mrow><mo>∗</mo></mrow></msubsup><mo>=</mo><mi>b</mi><mo>/</mo><msubsup><mrow><mi>τ</mi></mrow><mrow><mi>w</mi></mrow><mrow><mi>α</mi></mrow></msubsup></mrow></math></span>. Across five deterministic seeds, the model achieves <span><math><mrow><msup><mrow><mi>R</mi></mrow><mrow><mn>2</mn></mrow></msup><mo>=</mo><mn>0.9959</mn></mrow></math></span>–0.9985 for asymptotic resistance and <span><math><mrow><msup><mrow><mi>R</mi></mrow><mrow><mn>2</mn></mrow></msup><mo>=</mo><mn>0.853</mn></mrow></math></span>–0.910 for initial rate. The inversely recovered coefficient <span><math><mrow><mi>b</mi><mo>=</mo><mn>1.30</mn><mo>×</mo><mn>1</mn><msup><mrow><mn>0</mn></mrow><mrow><mo>−</mo><mn>4</mn></mrow></msup></mrow></math></span> K s/m matches the published plate-geometry value (<span><math><mrow><mn>1.45</mn><mo>×</mo><mn>1</mn><msup><mrow><mn>0</mn></mrow><mrow><mo>−</mo><mn>4</mn></mrow></msup></mrow></math></span> K s/m) within −10%, while the effective shear exponent (<span><math><mrow><mi>α</mi><mo>≈</mo><mn>1.43</mn></mrow></math></span>) reflects the localized operating envelope rather than a universal physical revision. Leave-One-Out cross-source validation shows that the physics constraint enables partial structural transfer across disparate geometries (<span><math><mrow><msup><mrow><mi>R</mi></mrow><mrow><mn>2</mn></mrow></msup><mo>></mo><mn>0.78</mn></mrow></math></span> on retained sources), though predictive skill remains constrained by domain boundaries. Operational scenario evaluations demonstrate that the model functions reliably as an engineering screening tool for cooling-water heat exchangers, with demonstrated validity strictly confined to CaCO<sub>3</sub> precipitation under the calibrated operational envelope.</div></div>","PeriodicalId":10019,"journal":{"name":"Chemical Engineering Research & Design","volume":"234 ","pages":"Pages 861-874"},"PeriodicalIF":4.0000,"publicationDate":"2026-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Chemical Engineering Research & Design","FirstCategoryId":"5","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0263876226006386","RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"2026/9/25 0:00:00","PubModel":"Epub","JCR":"Q2","JCRName":"ENGINEERING, CHEMICAL","Score":null,"Total":0}
引用次数: 0
Abstract
Crystallization fouling of CaCO3 in cooling water systems degrades heat transfer efficiency, yet the transient initial fouling rate and the asymptotic fouling resistance are traditionally evaluated using separate empirical correlations developed on distinct experimental rigs. Reconciling both descriptors within a single data-driven formulation is challenging because experimental datasets are scarce, heterogeneous across flow geometries, and sensitive to reporting conventions. This study proposes a compact Inverse Physics-Informed Neural Network (CaCO3-PINN) with two outputs and 9703 trainable parameters, informed by dimensionless supersaturation, temperature, and wall shear stress. The network is trained on 25 experimental points compiled from three independent sources (Watkinson 1983; Choi, via Liu 1999; Arsenyeva 2013) with a physics-constrained loss coupling rate and resistance via the Arsenyeva relationship . Across five deterministic seeds, the model achieves –0.9985 for asymptotic resistance and –0.910 for initial rate. The inversely recovered coefficient K s/m matches the published plate-geometry value ( K s/m) within −10%, while the effective shear exponent () reflects the localized operating envelope rather than a universal physical revision. Leave-One-Out cross-source validation shows that the physics constraint enables partial structural transfer across disparate geometries ( on retained sources), though predictive skill remains constrained by domain boundaries. Operational scenario evaluations demonstrate that the model functions reliably as an engineering screening tool for cooling-water heat exchangers, with demonstrated validity strictly confined to CaCO3 precipitation under the calibrated operational envelope.
CaCO3在冷却水系统中的结晶结垢降低了传热效率,然而,瞬态初始结垢率和渐近结垢阻力传统上是通过在不同的实验装置上开发的单独的经验关联来评估的。在单一数据驱动的公式中协调这两种描述符是具有挑战性的,因为实验数据集是稀缺的,跨流体几何形状的异质性,并且对报告惯例很敏感。本研究提出了一个紧凑的逆物理信息神经网络(CaCO3-PINN),具有两个输出和9703个可训练参数,由无因次过饱和度、温度和壁面剪切应力通知。该网络在来自三个独立来源(Watkinson 1983; Choi, via Liu 1999; Arsenyeva 2013)的25个实验点上进行训练,通过Arsenyeva关系Rf∗=b/τwα,具有物理约束的损耗耦合率和电阻。在5个确定性种子中,渐近抗性的R2=0.9959 ~ 0.9985,初始速率的R2=0.853 ~ 0.910。逆恢复系数b=1.30×10−4 K s/m与公布的板几何值(1.45×10−4 K s/m)在−10%以内匹配,而有效剪切指数(α≈1.43)反映的是局部操作包线,而不是普遍的物理修正。Leave-One-Out交叉源验证表明,物理约束允许跨不同几何形状的部分结构转移(在保留的源上R2>;0.78),尽管预测技能仍然受到域边界的限制。运行情景评估表明,该模型作为冷却水热交换器的工程筛选工具可靠地运行,其有效性严格限于校准操作包线下的CaCO3沉淀。
期刊介绍:
ChERD aims to be the principal international journal for publication of high quality, original papers in chemical engineering.
Papers showing how research results can be used in chemical engineering design, and accounts of experimental or theoretical research work bringing new perspectives to established principles, highlighting unsolved problems or indicating directions for future research, are particularly welcome. Contributions that deal with new developments in plant or processes and that can be given quantitative expression are encouraged. The journal is especially interested in papers that extend the boundaries of traditional chemical engineering.