{"id":"fa782b87-3be1-4f3d-b84b-953011d0d166","arxiv_id":"2411.18118","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A finite-element adjoint optimizer reconstructs thermal fields from displacement or strain measurements, outperforming temperature-only spatial interpolation in three synthetic civil-structure examples.","lead":"This paper presents a numerical method that recovers a structure's temperature field from measured displacements or strains, by repeatedly solving finite element models and minimizing the mismatch between computed and measured deformations. The authors test it on a plate, a bridge, and a dam model, and report lower errors than temperature-only spatial interpolation.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The method's forward model is mis-specified: Eq. (3)/(7) use the 3D fully-constrained coefficient αE/(1−2ν), but the plate and truss examples require αE/(1−ν) and αE respectively; if Kratos uses standard element loads, the printed equations do not describe the tested model.","rationale":"The reader's conditional verdict is appropriate. My stress-test goes one step further and locates a concrete internal inconsistency in the forward model that the reader did not identify. Eq. (3) misstates the 1D thermal stress, and Eq. (7) propagates the wrong coefficient into the FE load vector for the plane-stress and truss examples. If confirmed, this means the numerical demonstrations are not a faithful test of the derived equations, so the abstract's claim cannot be evaluated as written. This supports the CONDITIONAL verdict rather than ACCEPT: the paper needs corrected equations and/or a reproducible statement of the implemented Kratos thermal-load formulation, plus the noise/model-error tests. I do not recommend REJECT because the adjoint framework itself is standard and the error may be typographical; however, it is load-bearing and must be addressed.","tokens_in":30587,"tokens_out":14703,"duration_ms":145238,"concrete_test":"Run a single clamped triangular plane-stress element (or the Section 4.1 plate) in Kratos with a uniform ΔT and no external load, and compare the reaction forces against (a) Eq. (7) with αE/(1−2ν) and (b) the standard plane-stress thermal load αE/(1−ν). Also run one clamped truss element with ΔT and compare with αE. If the Kratos outputs match the standard coefficients, the printed equations do not describe the numerical model; if they match Eq. (7), the plate and bridge examples solve a different physical problem than the one claimed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.1 derives the thermal stress for a 'one-dimensional rod' but Eq. (3) gives σij = −αE/(1−2ν)ΔT δij, which is the triaxially constrained 3D stress, not the uniaxial rod stress (−αEΔT). The same coefficient is carried into the FE thermal load in Eq. (7). This is internally inconsistent with the examples: Section 4.1 uses plane-stress triangular elements, for which the correct thermal nodal-load coefficient is αE/(1−ν); Section 4.2 uses truss elements, for which the axial coefficient is αE. Only the 3D tetrahedral dam example matches the printed E/(1−2ν) coefficient. Since the paper claims an 'FE-method-dependent' adjoint procedure and uses Kratos for all simulations, the equations should describe the implemented forward model. If Kratos uses the standard element-consistent thermal eigenstrain loads, then Eqs. (3)-(7) are not the model actually tested, and the reported reconstructions cannot be reproduced from the paper. This is a load-bearing correctness gap in the central claim, independent of the additional noise/model-error limitations listed in Section 5.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a finite-element adjoint-based optimization procedure to reconstruct spatial temperature fields from displacement or strain measurements under known loads. The method models thermal strains as equivalent nodal loads, derives the adjoint gradient of a weighted least-squares misfit, and regularizes the ill-posed nodal temperature update with Vertex Morphing filtering. Three numerical examples are presented: a plate with a hole modeled with plane-stress triangles, a truss bridge, and a simplified Hoover Dam modeled with tetrahedra. In each case, a target temperature field is prescribed, displacement/strain data are generated synthetically, and the optimizer attempts to recover the field from sensor measurements. The results are compared with k-nearest-neighbor, ordinary kriging, and universal kriging interpolations using temperature sensors in identical configurations, with reported RMSE reductions up to 38.4%, 94%, and 40% for the three examples.","tokens_in":30884,"tokens_out":5188,"duration_ms":49626,"significance":"If the claims are borne out, the method would be a useful addition to structural health monitoring and digital twin workflows, since it would allow thermal fields to be reconstructed from deformation sensors alone, including scenarios where no temperature sensor is located in the heated region. The paper is clearly written and the optimization machinery is standard; the comparison with spatial interpolation methods is informative, and the authors are explicit about many limitations (e.g., noise, parameter uncertainty, sensor blind spots) in Section 5. The main concerns are that the printed thermal-load equations do not match the element types tested in Examples 1 and 2, and that the numerical validation uses an inverse-crime setup with noiseless synthetic data, so the stated claim of 'accurate reconstruction' is stronger than the evidence supports. These issues are fixable and do not invalidate the underlying idea.","major_comments":[{"comment":"Equation (3) gives sigma_ij = -alpha E/(1-2nu) Delta T delta_ij, which is the triaxially constrained 3D thermal stress, after the text introduces the equation as the thermal stress for a 'one-dimensional rod'. The uniaxial rod stress would be -alpha E Delta T. The same coefficient is carried into the finite-element thermal nodal load in Eq. (7). This is inconsistent with the examples: Section 4.1 uses plane-stress triangular elements, where the correct thermal nodal-load coefficient is alpha E/(1-nu), and Section 4.2 uses truss elements, where the axial coefficient is alpha E. Only the tetrahedral dam example of Section 4.3 matches the printed E/(1-2nu) coefficient. Since the paper describes an 'FE-method-dependent' procedure and all simulations were run in Kratos, Eqs. (3)-(7) should describe the implemented model; as printed, they do not, so a reader cannot reproduce the reported reconstructions from the paper. Please correct the equations or state clearly how the thermal-load coefficient is adapted for each element formulation.","section":"Section 3.1, Eqs. (3) and (7)"},{"comment":"All 'measured' data are generated by the same noiseless linear finite-element model that is used for the inversion, with perfectly known external loads, boundary conditions, and material properties. This is an inverse-crime setup. The abstract's conclusion that the approach 'can accurately reconstruct the target thermal distribution' is therefore not established for realistic measurements containing sensor noise, load uncertainty, or model error; these effects are listed as future work in Section 5. Please either add numerical experiments with additive sensor noise, perturbed material/load parameters, or a mismatched mesh/model, or qualify the central claim so that it is limited to the noiseless synthetic setting.","section":"Section 4, opening; Section 5"},{"comment":"The headline RMSE reductions (up to 38.4%, 94%, and 40%) compare the proposed displacement-based adjoint method with interpolation methods that use temperature sensors only. In configurations where no temperature sensor detects the heated region (e.g., the 6-sensor plate case in Figure 21a), the interpolation baseline returns essentially Delta T = 0 by construction, so large percentage improvements are expected and do not indicate superiority in a like-for-like information setup. The comparison is still meaningful for the claimed 'sensors not in the region of interest' scenario, but the abstract and Section 4.4.3 should state this asymmetry explicitly rather than presenting the percentages as a general accuracy comparison between reconstruction methods.","section":"Section 4.4.3, Tables 1-3"}],"minor_comments":[{"comment":"The strain misfit term in Eq. (8) uses the notation ums_ij, but the text defines the measured strain as sms_ij; please make the notation consistent.","section":"Section 3.2, Eq. (8)"},{"comment":"The symbol Js_ij is introduced to denote the displacement-strain relationship, but the cost function in Eq. (8) uses Is_ij; please align or define both symbols.","section":"After Eq. (12)"},{"comment":"The Young's modulus for the plate is given as E = 2e12 Pa, which is an order of magnitude above typical structural steel; if this value is intentional, please state so, otherwise it appears to be a typo.","section":"Section 4.1"},{"comment":"Several table captions refer to a single figure number (e.g., Table 2 refers to 'Figure 8') when the target distribution is only the bottom half of the displayed figure; please make the cross-references more precise.","section":"Tables 1-3"}],"recommendation":"major_revision","confidential_remarks":"The coefficient mismatch in Eqs. (3)-(7) is the most concrete correctness issue, but the inverse-crime validation is the main reason I cannot recommend acceptance: the numerical evidence, as presented, supports only a proof-of-concept under ideal synthetic conditions. I would encourage the authors to add at least noisy synthetic tests before resubmission, and to temper the abstract's 'accurately reconstruct' claim accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a usable engineering application of the adjoint-based inversion machinery the same group has used for weakness identification, transplanted to thermal field recovery. The paper is honest and the three test problems are well chosen. But there is a real discrepancy in the printed forward-model equations that needs to be resolved before I'd trust the paper's description of what was solved, and the validation is an inverse crime with no noise.\n\nWhat's genuinely new: the use of displacement/strain sensors to recover temperature fields in civil structures, with the demonstration that sensors don't have to sit in the heated region. The comparison against kNN and kriging is fair and useful; the RMSE tables let the reader judge, and the authors admit when the peaks come out low (the 6-sensor plate gives ~3.6°C against a 10°C target). The presentation is clear enough that a competent grad student could implement the method from the text, modulo the issue below.\n\nThe soft spots, in order of severity. First, Section 3.1 is internally inconsistent. The text introduces a one-dimensional rod, then writes the stress tensor as σ_ij = -αE/(1-2ν) ΔT δ_ij, which is the triaxially constrained 3D expression. That coefficient is carried into the FE thermal load in Eq. (7). But the plate example uses plane-stress triangles (correct coefficient αE/(1-ν)) and the bridge uses trusses (αE). Only the dam's tetrahedra match the printed formula. So either Kratos uses the standard element-consistent loads, in which case Eqs. (3)-(7) are not the model actually solved, or the simulation used the wrong coefficient for two of the three examples. Either way, the paper as written cannot be reproduced from its equations, and this touches the central claim. It's fixable—likely a presentation/derivation mistake—but it's the first thing a referee should ask about.\n\nSecond, the validation is an inverse crime: the \"measurements\" are generated by the same noiseless linear FE model that is inverted, with perfectly known loads, BCs, and material parameters. The authors list noise as future work, but for an SHM claim, a noise or model-error experiment is needed before I'd believe the accuracy statement.\n\nThird, no code or data is shipped, so none of the numbers can be independently checked.\n\nNone of this makes the method wrong, and the central idea is sensible. The paper is aimed at structural monitoring/digital-twin practitioners who want a physics-based alternative to interpolation; they'll get value even from the current form. I'd send it to peer review, but marked major revision: correct the thermal-stress coefficient, add at least one noise sensitivity study or clearly scope the claim to noiseless synthetic data, and consider releasing the input files.","headline":"A useful engineering application of standard adjoint inversion, but the paper has a serious unresolved discrepancy between the printed thermal-stress equations and the two tested element types, plus an inverse-crime validation.","tokens_in":31413,"tokens_out":4328,"would_cite":false,"duration_ms":38916,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["49K20","65K10","74F05","74S05"],"pacs":[],"model":"deepseek-v4-flash","headline":"Using only displacement or strain measurements, a finite-element adjoint optimizer reconstructs the structure's temperature field, and the sensors do not need to be inside the heated region.","keywords":["thermal field reconstruction","adjoint method","structural health monitoring","finite element method","inverse problem","vertex morphing","displacement sensors","kriging"],"falsifier":"Run the same inversion on a physical plate-with-a-hole specimen where the heated patch is known and measured with thermocouples, but feed the optimizer only strain-gauge readings that include realistic noise and unknown small load variations; if the recovered peak temperature stays far below the true value, the heated region is misplaced, or the root-mean-square error no longer beats kriging, the central claim fails.","tokens_in":30398,"feed_emoji":"🌡️","tokens_out":5842,"duration_ms":51705,"temperature":0.7,"pith_summary":"The paper proposes and tests a way to recover the temperature distribution of a structure from ordinary deformation measurements instead of temperature readings. It sets up an optimization problem: find the nodal temperature field whose thermal forces, under known mechanical loads, reproduce the displacements or strains recorded by sensors. The gradient of the mismatch is computed with an adjoint solve, and the field is updated by steepest descent with a Barzilai-Borwein step, with Vertex Morphing used to smooth ill-conditioned updates. On three synthetic examples (a plate with a hole, a truss bridge, and a simplified dam), the reconstructed fields localize the heated region even when no sensor lies inside it, and the method lowers root-mean-square error by up to 38.4%, 94%, and 40% compared with k-nearest-neighbor and kriging interpolation using temperature sensors in the same locations.","feed_headline":"Recover a structure's hidden temperature map from displacements alone","feed_subtitle":"Finite-element adjoint optimizer locates heated zones without any sensor in them, beating kriging in tests.","key_machinery":"The central object is the Lagrangian functional $\\mathcal{L}(u,\\Delta T,\\tilde{u}) = I(u,\\Delta T) + \\tilde{u}^t (f_{\\mathrm{ext}} + f_{\\Delta T} - K u)$, which couples the measurement mismatch $I$ to the finite-element equilibrium equations. The stationarity conditions produce the forward problem $K u = f_{\\mathrm{ext}} + f_{\\Delta T}$, the adjoint problem $K^t \\tilde{u} = \\frac{\\partial I}{\\partial u}$ (with $K^t = K$ for structural systems), and the gradient formula for the thermal field. Vertex Morphing acts as the regularizer: a kernel operator maps physical nodal temperatures through a control field, smoothing the gradient and suppressing high-frequency artifacts in the ill-posed inversion.","core_discovery":"The central claim is that the inverse problem of thermal field identification is tractable from displacement or strain data alone: given the load cases and the sensor readings, minimizing the weighted mismatch between measured and computed displacements yields a spatial temperature distribution that approximates the target field. The derivation shows that the gradient of the cost function with respect to the nodal temperature changes requires only the forward solve and an adjoint solve with the same stiffness matrix, since the stiffness matrix is symmetric for these structural systems, and the gradient takes the form $\\nabla I(\\Delta T) = \\frac{\\partial I}{\\partial \\Delta T} + \\tilde{u}^t \\frac{\\partial f_{\\Delta T}}{\\partial \\Delta T}$.","pith_inferences":["Inference: the method's success likely depends on the accuracy of the finite-element model and on the assumption that the measured deformation is purely thermal-mechanical; realistic sensor noise, load uncertainty, and model error are listed by the authors as future work and would probably degrade the reconstruction.","Inference: combining displacement or strain sensors with temperature sensors in the same cost function could fill the low-sensitivity 'blind spots' the paper observes near corners and edges.","Inference: since thermal loads enter the equilibrium equations as equivalent nodal forces through $f_{\\Delta T}$, the same adjoint machinery could be redirected to identify other distributed quantities, such as prestress or spatially varying material properties.","Inference: extending the cost function over time would allow recovery of transient thermal fields, at the cost of solving many more forward and adjoint problems."],"forward_implications":["Thermal fields can be reconstructed by deformation sensors alone, reducing reliance on temperature sensors in inaccessible or non-uniformly heated regions.","Sensor placement does not need to intersect the heated region; it only needs to be sensitive to the deformation that region produces.","Increasing sensor count and improving sensor placement sharpens both the magnitude and the localization of the recovered temperature field.","Vertex Morphing filtering consistently lowers reconstruction error, with root-mean-square error reductions of up to 16%, 24.4%, and 18% across the three examples.","At equal sensor counts, the proposed physics-based approach beats k-nearest-neighbor and kriging interpolation, with root-mean-square error reductions of up to 38.4%, 94%, and 40%."],"supporting_citations":[{"why":"Supplies the adjoint-based optimization machinery for identifying hidden spatial fields in structures, here adapted from weakness detection to thermal field recovery.","marker":"Airaudo et al. (2023)"},{"why":"Provides the Lagrangian and adjoint variational framework used to derive the gradient of the cost function.","marker":"Antil et al. (2018)"},{"why":"Defines Vertex Morphing, the filtering technique used to regularize the ill-conditioned inversion and smooth temperature updates.","marker":"Hojjat et al. (2014)"},{"why":"Describes the Barzilai-Borwein step-size rule used in the steepest-descent iteration.","marker":"Fletcher (2005)"},{"why":"Supplies the finite-element framework used to run all three numerical examples and generate the synthetic sensor data.","marker":"Dadvand et al. (2010)"},{"why":"Supplies the kriging implementation used as the spatial-interpolation baseline for comparison.","marker":"Murphy et al. (2024)"},{"why":"Documents that structural response in bridges is dominated by temperature effects, motivating deformation-based thermal reconstruction for structural health monitoring.","marker":"Glashier et al. (2024b)"}],"fun_headline_variants":["Heat map from displacements, no direct thermal sensors needed","Adjoint method recovers thermal fields from displacement data","Sensors can sit outside hot zones: still catch the heat map","Infer thermal map from strain: adjoint trick does it"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The demonstrations assume the displacement data come from the same noiseless linear finite-element model used in the inversion, with the loads, boundary conditions, and material properties known exactly.","fun_headline_variants_meta":{"raw":{"variants":["Heat map from displacements, no direct thermal sensors needed","Adjoint method recovers thermal fields from displacement data","Sensors can sit outside hot zones: still catch the heat map","Infer thermal map from strain: adjoint trick does it"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000896,"raw_usage":{"total_tokens":3867,"prompt_tokens":961,"completion_tokens":2906,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":577,"completion_tokens_details":{"reasoning_tokens":2838}},"tokens_in":577,"tokens_out":2906,"duration_ms":20959,"temperature":1.0,"reasoning_tokens":2838,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:30:02.153109+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same inversion on a physical plate-with-a-hole specimen where the heated patch is known and measured with thermocouples, but feed the optimizer only strain-gauge readings that include realistic noise and unknown small load variations; if the recovered peak temperature stays far below the true value, the heated region is misplaced, or the root-mean-square error no longer beats kriging, the central claim fails.","supporting_citations":[],"review_version":1}