Pith. sign in

REVIEW 3 major objections 5 minor 27 references

Towards a Physics Engine to Simulate Robotic Laser Surgery: Finite Element Modeling of Thermal Laser-Tissue Interactions

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper claims that a finite-element solution of the heat equation with a laser absorption source predicts laser-irradiated tissue temperature within about 2 °C in most benchtop conditions, enough to bring surgical lasers into robot…

desk verdict A sound, standard FEM model for laser-tissue thermal simulation that fills a real gap in surgical robotics, but the headline <2°C accuracy is partly in-sample because the two tissue parameters are tuned on the same data used to compute the error. read the letter →

arxiv 2411.14249 v1 pith:7XK7EIJB submitted 2024-11-21 cs.RO

Nicholas E. Pacheco , Kang Zhang , Ashley S. Reyes , Christopher J. Pacheco , Lucas Burstein , Loris Fichera This is my paper · ORCID
classification cs.RO MSC 65M6080A1992C50
keywords surgicalroboticslaser-tissueinteractionfiniteelementmethodthermalsimulationCO2laserheatconductionconvectionboundaryconditiontissuephantom
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper is trying to establish that a finite-element model of heat conduction, driven by a laser absorption source and cooled by convection, can simulate the temperature of tissue during surgical laser exposure accurately enough for robot simulators. The authors solve the heat equation on a cuboid mesh, with the laser power entering as a Lambert-Beer absorption term and surface heat loss handled through Newton's law of cooling. They verify the solver against infrared-camera measurements of CO2 laser irradiation on agar phantoms and ex-vivo chicken muscle, reporting average root-mean-square temperature error below 2 °C in most conditions and a worst observed error of 3.31 °C. If accurate, this gives surgical roboticists the missing thermal dynamics needed to train and test laser procedures in simulation rather than on tissue.

What carries the argument

The load-bearing machinery is the finite-element discretization of the heat-conduction initial/boundary-value problem: eight-node trilinear hexahedral elements assembled into the semidiscrete system $M\dot{d}+Kd=F$, advanced in time with the Crank-Nicolson method. The laser enters through an absorption source term proportional to the laser intensity, modeled with a Lambert-Beer law with a Gaussian beam profile, so focusing changes the local heating rate. The boundary physics enters through Newton's law of cooling $q_n=h(T_\infty-u)$ on air-exposed surfaces and a fixed-temperature Dirichlet heat sink on the bottom surface, with the natural-convection variant used in the verification scaling $h$ by $(T-T_\infty)^{1/4}$.

What would settle it

Run a hold-out calibration: tune the two coefficients on, say, the first twenty of the thirty reported runs, then compute RMSE on the remaining ten without retuning; if the held-out RMSE exceeds the reported 3.31 °C ceiling or grows with peak temperature, the current accuracy does not generalize beyond the tuning data.

Watch

Extended reading notes

Core claim

The central claim is that the thermal response of laser-irradiated tissue can be predicted by solving the heat-conduction initial/boundary-value problem $c_v \frac{\partial T}{\partial t} = \nabla\cdot(\kappa\nabla T)+S$ with a volumetric source $S=\mu_a I$, bottom-surface Dirichlet heat-sink boundary conditions, and Newton convection $q_n = h(T_\infty - T)$ on the air-exposed surfaces. In the authors' implementation, the laser enters through an absorption source term modeled with a Lambert-Beer law and a Gaussian beam profile, the conduction problem is discretized with eight-node trilinear hexahedral elements, and Crank-Nicolson time stepping advances the resulting linear system $M\dot{d}+Kd=F$. Comparing simulated surface temperatures to infrared-camera data from 30 experimental runs (two tissue types, three laser focus distances), the average root-mean-square error (RMSE) stayed below 2 °C for most conditions, with the largest observed RMSE 3.31 °C at the tightest focus. The authors state that the absorption coefficient $\mu_a$ and the convection coefficient $h$ were manually tuned to reduce modeling error on these experiments.

Load-bearing premise

The load-bearing premise is that the absorption coefficient and convection coefficient, which the authors tuned by hand to reduce error on the same experiments used for evaluation, keep their values on other tissue samples and beam conditions; if those values are sample-specific, the reported accuracy is a fit rather than a prediction.

Editorial extensions

If this is right

  • Surgical robot simulators can add a laser end effector whose thermal effects, not just geometry, are rendered, allowing training and reinforcement-learning agents to react to temperature.
  • Closed-loop laser control policies, such as regulating power or focal distance to hold surface temperature below a threshold, can be developed in simulation before tissue experiments.
  • The same finite-element solver can be applied to other wavelengths and tissue types by swapping the tissue parameters $\mu_a$, $\kappa$, $c_v$, and $h$.
  • With the authors' planned extensions, perfusion cooling and an Arrhenius damage model, the simulation would move from temperature histories to predictions of coagulation and cutting.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The reported accuracy is likely to degrade on a truly unseen tissue type or focus geometry, because the two coefficients that most control the fit were tuned on the same data used for evaluation; a held-out validation would show how much of the sub-2 °C error is predictive.
  • The model's own stated simplifications—no radiation, no temperature-dependent tissue properties, no scattering—imply that errors should grow at high temperatures and tight focus, so charting when that growth becomes clinically relevant is a natural follow-up.
  • Pairing the temperature history with an Arrhenius integral would let the simulator output a coagulation boundary, which is one step from predicting laser cutting depth for task planning.
  • If real-time performance can be demonstrated, a parameter-based finite-element model could replace the data-hungry machine-learning temperature predictors previously tried for laser surgery.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The manuscript proposes a finite-element model of thermal laser-tissue interaction as a component for surgical-robot simulators. Heat conduction with a Beer-Lambert laser source is discretized with trilinear hexahedral elements and Crank-Nicolson time stepping; boundary conditions include a bottom Dirichlet heat sink and natural convection on the remaining surfaces. The model is evaluated against benchtop CO2 laser experiments on agar phantoms and ex-vivo chicken muscle at three focal distances, reporting average RMSE mostly below 2 °C. The paper acknowledges several simplifications (no perfusion, scattering, radiation, or temperature-dependent properties), and the authors are transparent about the two manually tuned parameters.

Significance. If validated, the model would address a genuine gap, because common surgical robot simulators lack support for lasers and energy-based end effectors. The FEM derivation is standard and generally correct, and the paper includes a clear experimental comparison with multiple tissue types and focal distances. The main weakness is that the headline accuracy is an in-sample calibration result rather than a predictive validation, since the absorption coefficient and convection coefficient are tuned on the same experiments used to compute the reported RMSE. The authors are honest about this tuning and about the model's simplifications, but the central generalizability claim is therefore not yet established.

major comments (3)
  1. [Section III-B.2 and Tables II-III] The statement in Section III-B.2 that the absorption coefficient µa and the heat transfer coefficient h were 'manually tuned to reduce modeling error' makes the RMSE values in Tables II and III in-sample calibration errors. Because the same experimental runs were used both to select these parameters and to compute the reported accuracy, the abstract's claim that the model 'can accurately predict' tissue temperature is not supported by the current evidence. Please add a genuine out-of-sample validation, for example tune on some conditions and test on held-out focal distances or tissue types; alternatively, provide a sensitivity analysis showing that RMSE remains below the claimed threshold for a physically plausible range of µa and h. At minimum, the abstract and conclusion should distinguish calibration accuracy from predictive accuracy.
  2. [Section III-B.3, Eq. (27)] Equation (27) is dimensionally inconsistent as written. The factor 2P/(πw(z)) has units of W/cm, not W/cm^2, and the exponent -2(x^2+y^2)/w(z) is not dimensionless. A Gaussian beam intensity with the radially integrated power equal to P should be I(r,z) = 2P/(πw(z)^2) exp[-2r^2/w(z)^2] exp[-µa z]. This is load-bearing because the heat source S = µaI drives the entire thermal response used in the validation. Please correct the formula and state explicitly whether the implementation used the printed expression; if the implementation also has this form, the numerical results need to be recomputed with the corrected source.
  3. [Section III-A and III-B] Several inputs required to reproduce the verification are omitted: the laser power P, the beam waist w0, the laser wavelength λ, the fixed temperature assigned to the Dirichlet bottom boundary, and the exact expression used for the temperature-dependent convection coefficient (e.g., whether the implemented flux is q_n = h(T-T∞)^(5/4) or another form). These quantities enter Eqs. (24) and (27) and the boundary conditions, and without them the reported RMSE cannot be independently assessed. Please report all of these values or provide a code/data availability statement containing them.
minor comments (5)
  1. [Table I] The unit for µa is listed as 'cm'; the correct unit for an absorption coefficient is cm^-1. This should be corrected to avoid ambiguity with Eq. (27).
  2. [Fig. 4 and Tables II-III] Please specify explicitly whether the RMSE is computed per trial and then averaged, or computed on the trial-averaged temperature curve. The reported standard deviations suggest the former, but this should be stated in the text.
  3. [Section II-F and Fig. 2] The caption of Fig. 2 refers to a 'left-handed global frame'; please define the axes precisely and reconcile this with the z-axis of the optical beam used in Eq. (27), since boundary normal directions and the heat source both depend on the coordinate convention.
  4. [Section II-E] The time step Δt used in the Crank-Nicolson scheme is not reported. A brief statement of Δt and a short mesh/time-step convergence check would strengthen confidence in the numerical discretization.
  5. [General] Consider adding a code/data availability statement. The implementation is described as a C++ MEX file, but no repository or experimental dataset is referenced, which limits reproducibility.

Circularity Check

1 steps flagged · score 6.0 of 10

The reported <2°C 'prediction' accuracy is in-sample calibration: µa and h are manually tuned on the same experiments whose RMSE is reported.

  1. fitted input called prediction [Section III-B.2 (Tissue Physical Properties), with RMSE results in Section III-C and Tables II-III]
    ""Our approach for the selection of these two parameters was to manually tune them to reduce modeling error." ... "Average temperature tracking errors (root-mean-square error, RMSE) for each experimental condition are reported in Tables II and III.""

    The model's two tissue-specific free parameters, absorption coefficient µa and convection coefficient h, are manually tuned to reduce modeling error against the benchtop temperature measurements. The same set of 30 experimental runs (six conditions times five replicates) is then used to compute the RMSE values in Tables II and III that support the abstract's claim of sub-2°C prediction accuracy. Because no data are held out during tuning and no cross-validation is reported, the RMSE quantifies goodness of fit to the calibration data rather than predictive generalization; the 'prediction' claim reduces to the fitted parameters by construction.

full rationale

The FEM formulation itself is standard and self-contained: the heat equation, weak form, element assembly, Crank-Nicolson stepping, and boundary-condition treatment are derived from external references and are not defined in terms of the experimental outcomes. No load-bearing self-citation chain is present; the group's prior MCmatlab paper and the laser review are background, not justification for the model's validity. The circularity is confined to the experimental verification step. The authors state explicitly that µa and h were manually tuned to reduce modeling error, yet the abstract, discussion, and conclusion phrase the resulting RMSE as the model's ability to 'predict' temperature. Since the same experiments supply both the tuning objective and the reported error, the central accuracy claim is an in-sample fit. This is a substantial circularity concern for the headline claim, but it is not a definitional equivalence of the whole derivation; a proper hold-out or sensitivity analysis could restore the predictive interpretation. Score 6 reflects that the key predictive claim reduces, in part, to fitting the calibration constants.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The benchmark is the FEM solver itself; the only fitted quantities are the two tissue coefficients. No new physical entities are postulated.

free parameters (2)
  • absorption coefficient µa = 31 /cm (agar), 26 /cm (chicken)
    Tuned manually in Section III-B.2 to reduce modeling error.
  • heat transfer coefficient h = 0.022 W/cm2/degC (agar), 0.029 W/cm2/degC (chicken)
    Tuned manually together with µa in Section III-B.2.
assumptions (4)
  • standard math Bioheat equation with source S = µa I (Eqs. 1-2)
    Taken from Niemz [18]; this is the governing PDE for the entire model.
  • domain assumption Lambert-Beer intensity profile (Eq. 27)
    Assumes exponential attenuation without scattering; the paper acknowledges scattering is neglected.
  • domain assumption Newton's law of cooling with natural convection scaling (Eq. 23 and Section III-B.1)
    Convection boundary condition used for all free surfaces; scaling by (T - T∞)^{1/4} is a chosen model.
  • domain assumption Tissue thermal properties constant during exposure
    Acknowledged as limitation in Section IV; real properties vary with temperature.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Towards a Physics Engine to Simulate Robotic Laser Surgery: Finite Element Modeling of Thermal Laser-Tissue Interactions." pith.science (2026). https://pith.science/paper/7XK7EIJB

@misc{pith2026241114249,
  author       = {Pith},
  title        = {Pith review of: Towards a Physics Engine to Simulate Robotic Laser Surgery: Finite Element Modeling of Thermal Laser-Tissue Interactions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7XK7EIJB}},
  note         = {Machine review of arXiv:2411.14249}
}
read the original abstract

This paper presents a computational model, based on the Finite Element Method (FEM), that simulates the thermal response of laser-irradiated tissue. This model addresses a gap in the current ecosystem of surgical robot simulators, which generally lack support for lasers and other energy-based end effectors. In the proposed model, the thermal dynamics of the tissue are calculated as the solution to a heat conduction problem with appropriate boundary conditions. The FEM formulation allows the model to capture complex phenomena, such as convection, which is crucial for creating realistic simulations. The accuracy of the model was verified via benchtop laser-tissue interaction experiments using agar tissue phantoms and ex-vivo chicken muscle. The results revealed an average root-mean-square error (RMSE) of less than 2 degrees Celsius across most experimental conditions.

Figures

Figures reproduced from arXiv: 2411.14249 by the authors.

Figure 1
Figure 1. Graphical representation of an arbitrary two-dimensional (2D) [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Domain discretization and bi-unit domain. Without loss of [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Experiments used a surgical CO2 laser whose beam is delivered through an articulated (passive) arm. The tissue surface temperature was monitored with an infrared thermal camera at a rate of 20 frames per second (fps), and spatial resolution of 70 pixel cm . The distance between the beam’s focal point and the tissue surface (df ) was controlled with a robotic arm. the A655sc (Teledyne FLIR, Oregon, USA), and compared… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Simulated temperature results and experimental temperature results on agar tissue phantoms (top) and chicken muscle (bottom). The experimental [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

27 extracted references · 26 canonical work pages

  1. [1]

    Autonomy in surgical robotics,

    A. Attanasio, B. Scaglioni, E. De Momi, P. Fiorini, and P. Valdastri, “Autonomy in surgical robotics,” Annual Review of Control, Robotics, and Autonomous Systems , vol. 4, no. V olume 4, 2021, pp. 651–679, 2021

    work page 2021

  2. [2]

    Medical robotics—regulatory, ethical, and legal considerations for increasing levels of autonomy,

    G.-Z. Yang, J. Cambias, K. Cleary, E. Daimler, J. Drake, P. E. Dupont, N. Hata, P. Kazanzides, S. Martel, R. V . Patel, et al. , “Medical robotics—regulatory, ethical, and legal considerations for increasing levels of autonomy,” Science Robotics , vol. 2, no. 4, p. eaam8638, 2017

    work page 2017

  3. [3]

    Surgical gym: A high-performance gpu-based platform for reinforcement learning with surgical robots,

    S. Schmidgall, A. Krieger, and J. Eshraghian, “Surgical gym: A high-performance gpu-based platform for reinforcement learning with surgical robots,” in 2024 IEEE International Conference on Robotics and Automation (ICRA) . IEEE, 2024, pp. 13 354–13 361

    work page 2024

  4. [4]

    Orbit-surgical: An open-simulation framework for learning surgical augmented dexterity,

    Q. Yu, M. Moghani, K. Dharmarajan, V . Schorp, W. C.-H. Panitch, J. Liu, K. Hari, H. Huang, M. Mittal, K. Goldberg, and A. Garg, “Orbit-surgical: An open-simulation framework for learning surgical augmented dexterity,” in 2024 IEEE International Conference on Robotics and Automation (ICRA) , 2024, pp. 15 509–15 516

    work page 2024

  5. [5]

    Lapgym - an open source framework for reinforcement learning in robot-assisted laparoscopic surgery,

    P. M. Scheikl, B. Gyenes, R. Younis, C. Haas, G. Neumann, M. Wag- ner, and F. Mathis-Ullrich, “Lapgym - an open source framework for reinforcement learning in robot-assisted laparoscopic surgery,”Journal of Machine Learning Research , vol. 24, no. 368, pp. 1–42, 2023

    work page 2023

  6. [6]

    Ambf-rl: A real-time simulation based reinforce- ment learning toolkit for medical robotics,

    V . M. Varier, D. K. Rajamani, F. Tavakkolmoghaddam, A. Munawar, and G. S. Fischer, “Ambf-rl: A real-time simulation based reinforce- ment learning toolkit for medical robotics,” in 2022 International Symposium on Medical Robotics (ISMR) , 2022, pp. 1–8

    work page 2022

  7. [7]

    Open simulation environment for learning and practice of robot- assisted surgical suturing,

    A. Munawar, J. Y . Wu, G. S. Fischer, R. H. Taylor, and P. Kazanzides, “Open simulation environment for learning and practice of robot- assisted surgical suturing,” IEEE Robotics and Automation Letters , vol. 7, no. 2, pp. 3843–3850, 2022

    work page 2022

  8. [8]

    Surrol: An open-source reinforcement learning centered and dvrk compatible platform for surgical robot learning,

    J. Xu, B. Li, B. Lu, Y .-H. Liu, Q. Dou, and P.-A. Heng, “Surrol: An open-source reinforcement learning centered and dvrk compatible platform for surgical robot learning,” in 2021 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS) . IEEE, 2021, pp. 1821–1828

    work page 2021

Show all 27 references
  1. [9]

    Soft tissue simulation environment to learn manipulation tasks in autonomous robotic surgery,

    E. Tagliabue, A. Pore, D. Dall’Alba, E. Magnabosco, M. Piccinelli, and P. Fiorini, “Soft tissue simulation environment to learn manipulation tasks in autonomous robotic surgery,” in 2020 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), 2020, pp. 3261– 3266

    2020

  2. [10]

    An asynchronous multi-body sim- ulation framework for real-time dynamics, haptics and learning with application to surgical robots,

    A. Munawar and G. S. Fischer, “An asynchronous multi-body sim- ulation framework for real-time dynamics, haptics and learning with application to surgical robots,” in 2019 IEEE/RSJ International Con- ference on Intelligent Robots and Systems (IROS) , 2019, pp. 6268– 6275

    2019

  3. [11]

    A v-rep simulator for the da vinci research kit robotic platform,

    G. A. Fontanelli, M. Selvaggio, M. Ferro, F. Ficuciello, M. Vendittelli, and B. Siciliano, “A v-rep simulator for the da vinci research kit robotic platform,” in 2018 7th IEEE International Conference on Biomedical Robotics and Biomechatronics (Biorob) , 2018, pp. 1056– 1061

    2018

  4. [12]

    Faure, C

    F. Faure, C. Duriez, H. Delingette, J. Allard, B. Gilles, S. Marchesseau, H. Talbot, H. Courtecuisse, G. Bousquet, I. Peterlik, and S. Cotin, SOFA: A Multi-Model Framework for Interactive Physical Simulation. Berlin, Heidelberg: Springer Berlin Heidelberg, 2012, pp. 283–321

    2012

  5. [13]

    Pybullet, a python module for physics sim- ulation for games, robotics and machine learning,

    E. Coumans and Y . Bai, “Pybullet, a python module for physics sim- ulation for games, robotics and machine learning,” http://pybullet.org, 2016–2021

    2016

  6. [14]

    GitHub - NVIDIAGameWorks/PhysX: NVIDIA PhysX SDK — github.com,

    “GitHub - NVIDIAGameWorks/PhysX: NVIDIA PhysX SDK — github.com,” https://github.com/NVIDIAGameWorks/PhysX, [Ac- cessed 29-10-2024]

    2024

  7. [15]

    When the end effector is a laser: A review of robotics in laser surgery,

    H. C. Lee, N. E. Pacheco, L. Fichera, and S. Russo, “When the end effector is a laser: A review of robotics in laser surgery,” Advanced Intelligent Systems, vol. 4, no. 10, p. 2200130, 2022

    2022

  8. [16]

    µralp and beyond: Micro-technologies and systems for robot-assisted endoscopic laser microsurgery,

    L. S. Mattos, A. Acemoglu, A. Geraldes, A. Laborai, A. Schoob, B. Tamadazte, B. Davies, B. Wacogne, C. Pieralli, C. Barbalata, et al., “µralp and beyond: Micro-technologies and systems for robot-assisted endoscopic laser microsurgery,” Frontiers in Robotics and AI , vol. 8, p....

    2021

  9. [17]

    Bringing the light inside the body to perform better surgery,

    L. Fichera, “Bringing the light inside the body to perform better surgery,” Science Robotics, vol. 6, no. 50, p. eabf1523, 2021

    2021

  10. [18]

    M. H. Niemz, Laser-Tissue Interactions: Fundamentals and Applica- tions. Springer International Publishing, 2019

    2019

  11. [19]

    Learning Tempera- ture Dynamics on Agar-Based Phantom Tissue Surface During Single Point CO2 Laser Exposure,

    D. Pardo, L. Fichera, D. Caldwell, and L. Mattos, “Learning Tempera- ture Dynamics on Agar-Based Phantom Tissue Surface During Single Point CO2 Laser Exposure,” Neural Processing Letters, vol. 42, no. 1, 2015

    2015

  12. [20]

    Thermal supervision during robotic laser microsurgery,

    D. Pardo, L. Fichera, D. G. Caldwell, and L. S. Mattos, “Thermal supervision during robotic laser microsurgery,” in 5th IEEE RAS/EMBS International Conference on Biomedical Robotics and Biomechatronics . IEEE, Aug. 2014, pp. 363–368. [Online]. Available: https://ieeexplore.iee...

    2014

  13. [21]

    Identification of tissue optical properties during thermal laser-tissue interactions: An ensemble kalman filter- based approach,

    A. Arnold and L. Fichera, “Identification of tissue optical properties during thermal laser-tissue interactions: An ensemble kalman filter- based approach,” International Journal for Numerical Methods in Biomedical Engineering, vol. 38, no. 4, p. e3574, 2022

    2022

  14. [22]

    MCmatlab: an open-source, user-friendly, MATLAB-integrated three- dimensional Monte Carlo light transport solver with heat diffusion and tissue damage,

    D. Marti, R. N. Aasbjerg, P. E. Andersen, and A. K. Hansen, “MCmatlab: an open-source, user-friendly, MATLAB-integrated three- dimensional Monte Carlo light transport solver with heat diffusion and tissue damage,” Journal of Biomedical Optics , vol. 23, no. 12, p. 1, Dec. 2018...

    2018

  15. [23]

    Fish and T

    J. Fish and T. Belytschko, A First Course in Finite Elements . Chich- ester, England: John Wiley & Sons, Apr. 2007

    2007

  16. [24]

    Newton’s law of cooling—A critical assessment,

    C. T. O’Sullivan, “Newton’s law of cooling—A critical assessment,” American Journal of Physics , vol. 58, no. 10, pp. 956–960, Oct

  17. [25]

    Laser ablation of discs of agar gel,

    S. Rastegar, M. J. C. v. Gemert, A. J. Welch, and L. J. Hayes, “Laser ablation of discs of agar gel,” Physics in Medicine and Biology, vol. 33, no. 1, p. 133–140, Jan. 1988. [Online]. Available: http://dx.doi.org/10.1088/0031-9155/33/1/012

    1988 doi

  18. [26]

    P. A. Hasgall, F. Di Gennaro, C. Baumgartner, E. Neufeld, B. Lloyd, M. C. Gosselin, D. Payne, A. Klingenbock, and N. Kuster. (2022, Feb) IT’IS database for thermal and electromagnetic parameters of biological tissues. [Online]. Available: itis.swiss/database

    2022

  19. [1990]

    Available: https://pubs.aip.org/ajp/article/58/10/956/ 1053575/Newton-s-law-of-cooling-A-critical-assessment

    [Online]. Available: https://pubs.aip.org/ajp/article/58/10/956/ 1053575/Newton-s-law-of-cooling-A-critical-assessment

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.