REVIEW 4 major objections 6 minor 42 references
Thermal diffusivity characterization of impacted composites using evaporative cryocooling excitation and inverse physics-informed neural networks
T0 review · 4 major / 6 minor · reviewed 2026-08-04 · deepseek-v4-flash
Pith's one-line read Evaporative cryocooling combined with inverse physics-informed neural networks can measure thermal diffusivity in impacted composites with acceptable accuracy, with low-energy impacts raising diffusivity and high-energy impacts lowering it.
desk verdict Worth a serious look for the cryocooling excitation idea, but the accuracy claims need to be reined in before I'd trust the diffusivity numbers. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central machinery is the inverse physics-informed neural network (IPINN), a neural network whose loss function combines data fitting with the residual of the one-dimensional heat equation and boundary conditions, so that diffusivity is a learned variable rather than a predetermined parameter. It is paired with evaporative cryocooling: a portable spray of evaporating coolant cools the front face for about two seconds, a slow excitation whose shape is ill-defined for classical Parker analysis but can be assimilated by the network. The third component is terahertz time-domain spectroscopy, which measures sample thickness point-by-point so that diffusivity changes are not confused with impac
What would settle it
A direct check: scan a single impacted sample with the evaporative-cryocooling plus IPINN setup, then measure diffusivity at the same pixel with an independent 2-D Parker analysis or a contact probe; if the two disagree by more than the claimed ~30% at the damage center, or if a steel sample of known diffusivity is predicted only when initialized near its true value, the transfer claim fails.
Extended reading notes
Core claim
The paper's central claim is that thermal diffusivity can be recovered from a long, non-impulsive cooling excitation by inverting the one-dimensional heat-conduction equation with a neural network, despite the ill-posedness that defeats the Parker flash analysis. The inverse physics-informed neural network takes time, position, and measured temperature as inputs and outputs diffusivity, with physical residuals ensuring the prediction obeys the heat equation. Validation on simulated CFRP, GFRP, steel, wood, and concrete gives relative errors below 25%, and on real impacted flax/PLA-PBAT composites the predicted diffusivities at damage and sound sites fall within roughly 30% of Parker-method v
Load-bearing premise
The method treats the impacted composite as a locally homogeneous one-dimensional medium whose thermal response obeys the same diffusion equation used in training; if lateral heat flow, anisotropy, and damage-induced thickness gradients break that equivalence, the measured diffusivity values will be biased.
Editorial extensions
If this is right
- A compact, portable cryocooling rig can replace laboratory laser-flash equipment for thermal diffusivity mapping, at least for materials within the training distribution.
- Impact-damage severity can be read from thermal diffusivity maps once thickness variation is measured: low-energy impacts (2 J and 6 J) raise diffusivity in flax/PLA-PBAT composites, while high-energy impacts produce the opposite trend.
- Correcting for thickness changes after impact changes the measured diffusivity, so previous studies that ignored thickness variation may need revisiting.
- The IPINN approach extends to any material whose thermal response is governed by the heat equation, as demonstrated on CFRP, GFRP, steel, wood, and concrete.
- Combining THz-TDS thickness maps with thermal diffusivity maps gives a way to separate geometry-induced changes from material-property changes in nondestructive evaluation.
Reading between the lines
- Inference: Because the authors note lateral diffusion was not modeled, a 2-D or 3-D IPINN should reduce the roughly 30% error band at damage edges; this is testable by comparing the current 1-D network against a finite-volume simulation on the same thickness maps.
- Inference: The same evaporative cryocooling plus IPINN pipeline could be used to monitor subtle consolidation, aging, or processing changes in polymers, not just impacts, since the method is portable and does not require optical access for a laser.
- Inference: The reported sign reversal, diffusivity increases after low-energy impacts but decreases after high-energy ones, suggests thermal diffusivity maps could serve as a severity triage tool, but this needs independent validation against CT or destructive sectioning to confirm the underlying micromechanisms.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes an inverse physics-informed neural network (IPINN) combined with evaporative cryocooling excitation to estimate the out-of-plane thermal diffusivity of impacted flax/PLA-PBAT composites. THz time-domain spectroscopy is used to map post-impact thickness, and the Parker laser-flash method serves as reference. The IPINN is first validated on synthetic one-dimensional heat-conduction data for five homogeneous materials (CFRP, GFRP, steel, wood, concrete), then applied to experimental thermograms of four composite samples with 2 J and 6 J impacts. The authors report simulation errors below 25% and experimental relative errors up to 34.59%, and conclude that low-energy impacts increase thermal diffusivity through matrix densification and improved interfacial bonding, while high-energy impacts decrease it. The central claim is that simulations and experiments demonstrate the feasibility and accuracy of the IPINN-based approach.
Significance. If validated, the method would offer a compact, portable, and low-cost excitation source for thermal diffusivity mapping, addressing a real practical limitation of flash-lamp and laser systems. Using THz-TDS to correct for impact-induced thickness variation is a meaningful improvement over prior work that ignored such thickness changes. The idea of using a physics-informed inverse network to handle non-impulsive, long-pulse excitation is interesting and potentially useful. However, the current evidence does not robustly support the accuracy claim: the simulation validation is a closed-loop inversion of the generating forward model, the experimental reference is the Parker method applied to impacted anisotropic samples, and the paper's own Table III contradicts its stated error bound. The strengths are the novel excitation modality, the explicit use of thickness mapping, and the physics-informed formulation; the weakness is a validation gap that must be closed before the feasibility/accuracy statement is credible.
major comments (4)
- [Section IV.C, Table III] Table III reports a maximum relative error of 34.59% (RTL_2_D), yet the text immediately below states 'The maximum relative error is no more than 30%.' This is a direct numerical contradiction. Several predictions also cluster near 0.3–0.5 mm2/s (e.g., RT_1_D, RT_1_S, RT_2_S, RTL_1_D), with errors of 23–28%, suggesting that the IPINN may not be resolving damage-dependent diffusivity variations. To support the stated accuracy, the authors should report per-pixel or per-region error statistics, confidence intervals, or an independent validation measurement.
- [Section IV.B, Table II and Section IV.C] The simulation validation is a closed loop: synthetic temperature curves are generated from a forward one-dimensional heat-conduction model with known alpha, and the IPINN then recovers alpha from those same curves. This tests only the inversion of the generating model. The experimental validation then applies the network, trained/configured for homogeneous isotropic 1D materials (Table I), to impacted anisotropic flax/PLA-PBAT composites with localized damage, spatially varying thickness, and possible lateral diffusion. The paper does not demonstrate that the model generalizes to this more complex setting, nor that the Parker-method 'real' values from Fig. 6 are an adequate ground truth for impacted anisotropic plates. Please provide validation on a synthetic orthotropic/impacted case or an independent experimental reference.
- [Section II.B, Eq. (7)] Equation (7) is introduced as 'the one-dimensional heat conduction problem' but contains T_xx, T_yy, and T_zz terms. This inconsistency is load-bearing because the simulation and training setup (Section III.B) uses only isotropic materials with a single diffusivity, while the experimental samples are anisotropic and have localized thickness variations. If the forward model is 1D, Eq. (7) should be T_t = alpha T_xx; if it is 3D, the simulations must specify orthotropic diffusivities and boundary conditions that match the experimental geometry. Please clarify and align the model with the experimental setup.
- [Section IV.B, initial-value sensitivity] The paper reports that the steel prediction error drops from 24.17% to 9.45% only after manually changing the initial thermal diffusivity from 1e-8 to 1e-6 m2/s, and concludes that 'it is essential to initialize the thermal diffusivity within a reasonable range.' This introduces a free parameter that materially affects the result. For impacted composites, where local diffusivity is unknown a priori, the sensitivity of the IPINN to the initialization should be quantified, and the initialization used for the experimental data in Section IV.C should be justified rather than stated as a fixed value.
minor comments (6)
- [Abstract / Section III.B] The sentence 'due to extremely small magnitude of diffusivity is extremely small' is ungrammatical; please revise.
- [Section IV.C, Fig. 8] The text says 'Fig. 8 shows the thermal diffusivity maps obtained using the laser flash method,' but the caption says 'based on laser flash method and evaporative cryocooling method.' Clarify that the Parker-method equations were applied to inverted evaporative-cryocooling data.
- [Table III and Fig. 6] The label 'Real' in Table III is misleading; these values are Parker-method estimates from Fig. 6 after denoising and averaging, not independent reference measurements. Use a more neutral term such as 'Parker-method estimate.'
- [Section III.B] 'NVIDIA 4060 Titan GPUs' appears to be a typo; please correct the hardware description.
- [Section II.B, Eqs. (12)–(16)] The notation for L_data, L_ic, L_bc and the relation to Eq. (12) is confusing; please define each loss term explicitly and consistently.
- [Section IV.C] The statement 'The training configuration mirrors the setup used in the simulation' is vague. Specify the number of experimental observation points, how the THz-TDS thickness map is supplied to the network, and how noisy data are weighted in the loss.
Circularity Check
Simulation validation is a self-consistency check; experimental validation is independent.
-
self definitional
[Section IV.B (Simulation Validation for IPINNs), Table II]
"To evaluate the feasibility and accuracy of the proposed IPINNs combined with evaporative cryocooling method, numerical simulations were first conducted. ... Table II shows the thermal diffusivity results predicted by IPINNs."
The simulated temperature curves are generated by solving the same 1D heat-conduction equation used in the IPINN physical residual (Eqs. (7) and (11)). For a curve generated with diffusivity a_real, the IPINN output a_pred minimizes f = T_t - a_pred T_xx, so a_pred equals a_real wherever T_xx is nonzero. The 'predicted' values in Table II therefore coincide with the generating values by construction; the relative errors measure only the optimizer's ability to invert the forward model, not the validity of the measurement model for impacted anisotropic composites. The experimental validation (Table III) is not similarly forced because it compares IPINN outputs to independently obtained Parker-method values from flash-lamp excitation.
full rationale
The paper's central experimental claim is not circular: IPINN outputs are obtained from evaporative-cryocooling thermograms and compared with Parker-method values derived from independent flash-lamp experiments, without training on those Parker values. The circular step is confined to the simulation validation in Section IV.B/Table II, where synthetic data are generated by the same 1D heat equation that defines the IPINN residual, so recovering the generating alpha is guaranteed by construction. This inflates the abstract's statement that 'Simulations and experimental results demonstrated the feasibility and accuracy.' The experimental portion remains independent support, although its strength is weakened by the small sample set and the paper's own internal inconsistency (Table III lists a 34.59% relative error while the text claims 'The maximum relative error is no more than 30%'); that inconsistency is a correctness/support issue, not a circularity. Self-citations (Refs. [5], [37]) are used for contrasting impact-energy regimes and are not load-bearing for the IPINN method itself.
Assumptions & free parameters
free parameters (3)
- Initial thermal diffusivity in IPINN =
1e-8 m2/s; adjusted to 1e-6 m2/s for steel
- Cooling duration for experiments =
2 s
- Loss weighting hyperparameters =
lambda formula not fully specified; boundary weight 0.1; gradient clip 1.0; learning rate 1e-4; 4000 iterations
assumptions (6)
- domain assumption One-dimensional isotropic heat conduction with constant thermal diffusivity governs the measurement (Eq. 7, Section II.B).
- domain assumption The evaporative cryocooling excitation can be represented as a prescribed surface heat flux Q(x,y,0,t) with a known time profile (Eq. 9).
- ad hoc to paper The neural network trained on homogeneous materials generalizes to impacted inhomogeneous samples.
- domain assumption The Parker laser flash method provides reliable ground-truth thermal diffusivity for impacted composites after THz thickness correction (Section IV.A, IV.C).
- domain assumption Boundary conditions include natural convection h[T - T_inf] with Grashof number Gr >> 2000 (Eq. 10, Section II.B).
- ad hoc to paper Initial thermal diffusivity is set within a reasonable range before training.
Cite this review
Pith. "Pith review of Thermal diffusivity characterization of impacted composites using evaporative cryocooling excitation and inverse physics-informed neural networks." pith.science (2026). https://pith.science/paper/N57UA4UV
@misc{pith2026250910898,
author = {Pith},
title = {Pith review of: Thermal diffusivity characterization of impacted composites using evaporative cryocooling excitation and inverse physics-informed neural networks},
year = {2026},
howpublished = {\url{https://pith.science/paper/N57UA4UV}},
note = {Machine review of arXiv:2509.10898}
}
read the original abstract
The thermal diffusivity measurement of impacted composites using pulsed methods presents an ill-posed inverse problem influenced by multiple factors such as sample thickness, cooling duration, and excitation energy. In this study, a novel excitation method, evaporative cryocooling, was introduced for measuring the thermal diffusivity of tested samples. Compared to conventional excitation modalities, evaporative cryocooling excitation is compact, portable, and low cost. However, evaporative cryocooling cannot be considered a pulsed method due to its prolonged excitation duration. In general, it is difficult to measure thermal diffusivity based on non-impulsive pulsed excitation at times commensurate with the pulse duration, often due to ill-defined pulse shape and width and the subsequent potentially complicated thermal response which may be subject to diffusive broadening. To address this challenge, inverse physics-informed neural networks (IPINNs) were introduced in this work and integrated with an evaporative cryocooling method. The Parker method combined with a photothermal method was employed as a reference. To improve the accuracy of both IPINNs and Parker methods, terahertz time-domain spectroscopy (THz-TDS) was employed for measuring the thickness of impacted composites. Simulations and experimental results demonstrated the feasibility and accuracy of the IPINN-based approach.
Reference graph
Works this paper leans on
-
[1]
In this study, a novel excitation method—evaporative cryocooling—was introduced for measuring the thermal diffusivity of tested samples
1 > REPLACE THIS LINE WITH YOUR MANUSCRIPT ID NUMBER (DOUBLE-CLICK HERE TO EDIT) < Thermal diffusivity characterization of impacted composites using evaporative cryocooling excitation and inverse physics-informed neural networks Pengfei Zhu, Graduate Student Member, IEEE, Hai Zhang*, Member, IEEE, and Stefano Sfarra, Fabrizio Sarasini, Rubén Usamentiaga, ...
2017
-
[5]
This is because low-energy impacts enhance matrix densification, improve PLA/PBAT interfacial bonding, and align fibers favorably, thereby increasing thermal diffusivity
and [37]). This is because low-energy impacts enhance matrix densification, improve PLA/PBAT interfacial bonding, and align fibers favorably, thereby increasing thermal diffusivity. In contrast, high-energy impacts induce microcracks, fiber breakage, and interfacial debonding, leading to a decrease in thermal diffusivity. The evaporative cryocooling metho...
2017
-
[6]
https://doi.org/10.1016/j.optlaseng.2021.106927
arXiv 2021
-
[10]
https://doi.org/10.1007/s10443-016-9538-8
-
[11]
https://doi.org/10.1016/j.actamat.2019.05.010
-
[13]
in 1973, who measured the thermal diffusivity in the in-plane direction of a solid sample. The TTG technique uses the interference of two pulsed laser T 2 > REPLACE THIS LINE WITH YOUR MANUSCRIPT ID NUMBER (DOUBLE-CLICK HERE TO EDIT) < beams to generate a spatially periodic thermal grating on a sample surface, the relaxation dynamics of which is monitored...
1973
-
[19]
https://doi.org/10.1016/j.ndteint.2024.103066
arXiv 2024
-
[20]
https://doi.org/10.1364/AO.53.007985
Show all 42 references
-
[21]
https://doi.org/10.1016/0032-5910(79)85016-0
-
[24]
Standard test method for thermal diffusivity by the flash method,
ASTM E1461-13, “Standard test method for thermal diffusivity by the flash method,” ASTM International: West Conshohocken, PA, USA, 2022, https://10.1520/E1461-13R22
2022 doi
-
[25]
Standard practice for thermal diffusivity by the flash method,
ASTM E2585-09, “Standard practice for thermal diffusivity by the flash method,” ASTM International: West Conshohocken, PA, USA, 2022, https://10.1520/E2585-09R22
2022 doi
-
[29]
https://doi.org/10.3390/app13106135
-
[32]
https://doi.org/10.1016/j.ndteint.2024.103133
2024
-
[33]
https://doi.org/10.1016/j.ndteint.2023.102938
2023
-
[34]
https://doi.org/10.1364/OL.30.000029
-
[35]
https://doi.org/10.1016/j.infrared.2024.105563
2024
-
[36]
https://doi.org/10.1016/j.infrared.2024.105353
2024
-
[37]
https://doi.org/10.1016/j.ndteint.2024.103083
2024
-
[38]
Pengfei Zhu (Graduate Student Member, IEEE) received the B.Eng
https://doi.org/10.1088/2631-7990/ada837. Pengfei Zhu (Graduate Student Member, IEEE) received the B.Eng. degree in engineering mechanics from North University of China, Taiyuan, China, in 2019, and the M.Eng. degree in solid mechanics from Ningbo University, Ningbo, China, in
2019 doi
-
[39]
degree in electrical engineering with Université Laval, Québec, Canada
He is currently working toward the Ph.D. degree in electrical engineering with Université Laval, Québec, Canada. His research interests include non-destructive testing, infrared thermography, deep learning, terahertz time-domain spectroscopy, and photothermal coherence tomogra...
2004
-
[41]
He has more than 90 international peer-reviewed journal papers
Since June 2016, he has been an Assistant Professor with the Faculty of Civil and Industrial Engineering, Sapienza University of Rome. He has more than 90 international peer-reviewed journal papers. His research interests include the impact response of composite materials for ...
2016
-
[1961]
https://doi.org/10.1063/1.1728417
-
[1973]
https://doi.org/10.1063/1.1662160
-
[1977]
degree of electro-mechanical engineering and the Ph.D
He received the M.Sc. degree of electro-mechanical engineering and the Ph.D. degree from the Vrije Universiteit Brussel (VUB), Ixelles, Belgium, in July 2000 and 2008, respectively. Since 2012, he has been a Full Professor with the University of Antwerp, Antwerp, Belgium, and ...
2000
-
[1979]
#!$[1+2∑(−1)%exp (&%
In this method, an incident modulated radiation heats periodically the sample surface, and an infrared camera records the periodic temperature fluctuations. By comparing the phase shift between original waveform and recorded signals, the thermal diffusivity can be extracted. P...
1980
-
[1980]
https://doi.org/10.1063/1.91395
-
[1982]
https://doi.org/10.1051/rphysap:01982001704022700
-
[1983]
https://doi.org/10.1007/BF00616610
-
[1995]
https://doi.org/10.1063/1.1146123
-
[1998]
https://doi.org/10.1063/1.1148495
-
[2005]
https://10.1007/s10765-005-8603-6
-
[2007]
https://doi.org/10.1063/1.2812542
-
[2011]
Currently, he is an Associate Professor at DIIIE-UNIVAQ and an Adjunct Professor at Université Laval, Quebec, QC, Canada
He worked as a Post-Doctoral Fellow at UNIVAQ until October 2017, where he became a Researcher in a fixed-term contract at the Department of Industrial and Information Engineering and Economics (DIIIE), UNIVAQ. Currently, he is an Associate Professor at DIIIE-UNIVAQ and an Adj...
2017
-
[2014]
https://doi.org/10.1016/j.compositesb.2014.01.034
2014 doi
-
[2015]
https://doi.org/10.1007/s11340-014-9976-x
-
[2016]
https://doi.org/10.1016/j.ndteint.2016.06.008
2016 doi
-
[2017]
https://doi.org/10.1115/1.4034605
-
[2019]
https://doi.org/10.1088/1361-6501/ab018e
-
[2022]
https://doi.org/10.1364/AO.460385
-
[2023]
https://doi.org/10.1080/17686733.2023.2215102
2023
-
[2024]
https://doi.org/10.1016/j.compositesa.2024.108220
2024
-
[2025]
https://doi.org/10.1016/j.ndteint.2025.103361
2025
Reviewed August 4, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.