Pith. sign in

REVIEW 2 major objections 5 minor 20 references

Accurate Computation of Activated Volume in Electromagnetic Heating

T0 review · 2 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash

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

Pith's one-line read Combining fine- and coarse-grid triangulated volumes cancels the leading error and gives higher-order-accurate isosurface volumes at no extra simulation cost.

desk verdict A useful, modest numerical trick—Richardson extrapolation on triangulated isosurface volumes—supported by clean experiments; main gap is the underspecified boundary cap in the heating test. read the letter →

arxiv 2607.25994 v1 pith:AVD45TOO submitted 2026-07-28 math-ph math.MP

Hongyun Wang , Parthiv Seetharaman , Shannon E. Foley , Hong Zhou This is my paper · ORCID
classification math-phmath.MP MSC 65D1865B05
keywords isosurfacevolumecomputationtriangulationextrapolationdivergencetheoremelectromagneticheatingactivatedrectangulargrid
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

This paper claims that the volume enclosed by an isosurface of a function known only on a rectangular grid can be computed to higher accuracy by triangulating the isosurface, computing the enclosed volume from the divergence theorem as a sum of tetrahedron volumes, and then extrapolating the fine-grid result against the volume obtained from a downsampled coarse grid. The extrapolated combination Vextrap = (4V(Δx) − V(2Δx))/3 cancels the dominant second-order error of the raw triangulation, leaving an error of order o((Δx)^2) rather than O((Δx)^2). Seven test problems spanning spheres, ellipsoids, superellipsoids, ring tori, and a model of electromagnetic skin heating support the claim. The importance is practical: the activated skin volume in millimeter-wave exposure—where temperature is numerically computed on a grid—can be assessed more accurately using information already present in the solution, without refining the simulation.

What carries the argument

The central mechanism is extrapolation applied to a volume computed from a triangulated isosurface. The volume is expressed via the divergence theorem as one third of a surface integral of (x − x0)·dS; over each planar triangle this integral reduces to one sixth of a triple scalar product, so the enclosed volume is a summed set of tetrahedron volumes. The key identity is the assumed second-order error expansion V(Δx) = V_true + C2(Δx)^2 + o((Δx)^2) for both the fine grid and the downsampled coarse grid, which makes the combination (4V(Δx) − V(2Δx))/3 cancel C2(Δx)^2.

What would settle it

For a known-volume isosurface, compute V(Δx) and V(2Δx) on a sequence of nested grids and form the ratio (V(2Δx) − V_true)/(V(Δx) − V_true). If the claimed expansion holds, this ratio should approach 4; persistent deviations from 4 across resolutions—for example, for an isosurface aligned to pass through grid vertices—would falsify the method's accuracy claim.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is the identity Vextrap(Δx) ≡ (4V(Δx) − V(2Δx))/3 = V_true + o((Δx)^2), where V(Δx) is the volume enclosed by the linear-interpolation triangulation of the isosurface on a grid of spacing Δx and V(2Δx) is the same quantity computed after downsampling the grid. Because V(Δx) has the asymptotic form V_true + C2(Δx)^2 + o((Δx)^2) and V(2Δx) has the same leading coefficient C2, the combination eliminates the leading error term. The paper's seven analytic test problems all show raw triangulation error following a C2(Δx)^2 trend and extrapolated error far below it, including for a non-symmetric activated-skin-volume problem.

Load-bearing premise

The paper assumes, without proof, that the triangulated-volume error has the form V(Δx) = V_true + C2(Δx)^2 + o((Δx)^2) and that the coefficient C2 is identical for the fine grid and the downsampled coarse grid; the expansion is inferred from the analogy with line-segment curve approximation and can break down for surfaces passing through grid vertices, non-manifold topology, or boundary-clipped domains.

Editorial extensions

If this is right

  • A single fine-grid temperature solution yields both V(Δx) and, by downsampling, V(2Δx); no extra PDE solve is needed for the accuracy gain.
  • In every test case the extrapolated volume error is substantially smaller than the raw triangulation error, so the method can replace much finer grids in practice.
  • The procedure requires no symmetry and applies to any discrete 3D function, including the non-axially-symmetric skin temperature distributions that arise with moving or oblique beams.
  • More accurate activated volumes imply better estimates of the thermal nociceptor activation region and potential thermal injury in millimeter-wave exposure assessments.
  • The observed second-order trend in all seven problems supports using the extrapolation as a cheaper accuracy boost for isosurface volume computation.

Reading between the lines

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

  • (Editorial inference) Because extrapolation can be applied recursively, combining V(Δx), V(2Δx), and V(4Δx) could push the error to higher order than the paper's single-step formula, at the cost of requiring an additional downsampling.
  • (Editorial inference) The assumed expansion may fail when the isosurface passes through grid vertices, has sharp creases, or is clipped by a boundary—conditions near the skin surface z = 0 in the application, where the paper does not specify how the boundary cap is included in the triangulation integral.
  • (Editorial inference) The same triangulation-plus-extrapolation recipe could improve other surface integrals over isosurfaces, such as integrated heat flux or surface area, but the paper does not test those.
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

2 major / 5 minor

Summary. The paper proposes a method for computing the volume of the region where a function defined on a rectangular grid exceeds a threshold. The boundary is approximated by a triangulated isosurface with vertices obtained by linear interpolation on grid lines, and the enclosed volume is computed as one sixth of the sum of signed scalar triple products (Eq. 6). The central device is Richardson extrapolation: from the fine-grid volume V(Δx) and the downsampled coarse-grid volume V(2Δx), Eq. 9 defines Vextrap = (4V(Δx) − V(2Δx))/3. This is justified by assuming the error expansion V(Δx) = Vtrue + C2(Δx)^2 + o((Δx)^2) with the same coefficient C2 for the coarse grid (Eqs. 7–8). The method is tested on seven analytic benchmark volumes — sphere, ellipsoids, superellipsoids, and elliptic ring tori, in aligned and rotated orientations — and on a simplified electromagnetic-heating model with an analytic activated volume. In all tests the error in V(Δx) is reported to scale roughly as (Δx)^2, and extrapolation markedly reduces the error.

Significance. If the extrapolation assumption holds, the method is simple, inexpensive, and practically valuable: it requires no additional PDE solves, only a downsampled coarse representation of an already computed grid function. The numerical evaluation is a genuine strength: the benchmarks are independent analytic solutions, no parameters are fitted to the target volumes, and the extrapolation formula is derived before testing rather than tuned afterward. The convergence behavior is consistent across smooth closed surfaces with nontrivial topology (tori) and across orientations. However, the central error expansion is asserted rather than proved, and the electromagnetic-heating application in §4.3 leaves the treatment of the planar boundary cap unspecified. These two issues bear directly on whether the headline claim — Vextrap = Vtrue + o((Δx)^2) — is established, and they must be addressed before the paper can be accepted.

major comments (2)
  1. [§2, Eqs. (7)–(9)] The entire extrapolation procedure rests on the asymptotic expansion V(Δx) = Vtrue + C2(Δx)^2 + o((Δx)^2) with the same coefficient C2 for the downsampled coarse grid. This is asserted, not derived; the text says only 'Thus, we expect' and draws an analogy with polygonal curve approximation. For cell-based triangulations of a level set, the per-cell error is not obviously of the form C2 times the global mesh spacing squared summed over the surface, and the coarse-grid triangulation is not the same geometric object as the fine-grid triangulation. The equality of C2 between the two grids is especially nontrivial. Since Eq. (9) is the paper's central contribution, I ask for either a local error analysis that yields Eq. (7)–(8), or an explicit statement that extrapolation is a heuristic validated only empirically. The paper should also discuss and test failure modes, such as isosurfaces pass
  2. [§4.3, Eq. (28)] The activated skin region D in the electromagnetic-heating test is bounded by the isosurface and the skin surface z = 0. Equation (6), however, is derived from the divergence theorem for a closed surface S = ∂D. The manuscript states that D is 'enclosed by the isosurface and the skin surface' but never specifies how the planar cap on z = 0 is included in the triangulation or in the tetrahedron sum. If only the isosurface triangles are used, the sum in Eq. (6) is not a well-defined volume and depends on the arbitrary point x0; the reported V(Δx) and V(2Δx) for the headline application then cannot be reproduced or verified from the text. The authors must specify the cap triangulation (or its equivalent) and the choice of x0, and clarify how the finite computational domain is handled.
minor comments (5)
  1. [§3.1] The sentence 'the extrapolation component does not require addition data or information' contains a typo: 'addition' should be 'additional'.
  2. [Figures 2–10] The convergence claims are supported only by visual inspection of log-log plots. Reporting the fitted slopes of the error curves, or a small table of errors at each Δx, would make the O((Δx)^2) statement and the improvement factor quantitative and easier to verify.
  3. [Abstract and §1] The paper is framed as computing the volume 'enclosed by an isosurface,' but in the physical skin-heating problem the activated region is truncated by the skin surface z = 0. This distinction is made in §4.3 but should appear earlier, since it affects the definition of the volume being approximated.
  4. [Eq. (25)] The analytical solution for U(z,t) has a typesetting issue: the exponentials and erfc arguments appear to have lost fraction bars and parentheses. Please check the formula carefully, since it is used to generate the benchmark volume in §4.3.
  5. [§5] The summary says both curve and surface approximations 'are expected to produce second order approximations.' This is an honest hedging, but the earlier text and the abstract state the second-order behavior as a fact. Reconcile these statements or add a caveat that Eq. (7) is an assumption to be verified.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: Eq. (9) is an algebraic consequence of the stated error expansions and is tested against independent analytic volumes; self-citations are contextual, not load-bearing.

full rationale

The central derivation chain is: exact divergence-theorem identity (2), planar-triangle approximation (3)-(5), and the tetrahedral sum V(Δx) in (6). The extrapolation formula (9) is obtained by algebraically combining the assumed expansions (7) and (8); it is not fitted to any benchmark volume. The seven test problems use independently known analytic volumes, and the exact volume is used only to evaluate errors, not to determine any parameter in the method. Thus there is no fitted-input-called-prediction structure. The unproved same-C2 quadratic error expansion is a modeling assumption and a correctness limitation, but it is not equivalent by construction to the claimed o((Δx)^2) result: it is the hypothesis under which (9) is an algebraic consequence. The §4.3 boundary-cap ambiguity is a well-posedness/implementation gap, not circularity. Self-citations [11,12,17] provide coordinate conventions and the displayed analytic solution U(z,t); they do not justify the extrapolation identity or forbid alternatives, and the method is independently supported by the other benchmarks. Under the quoted-evidence standard, no circular step is present.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The method itself has no fitted parameters; all listed constants are test-case choices. The central load-bearing assumption is the second-order error expansion with a grid-independent leading coefficient C2; if that fails, the extrapolation does not improve accuracy. The EM test additionally assumes the triangulated surface plus skin-surface cap forms a closed boundary, which is not explicitly described.

assumptions (5)
  • standard math Divergence theorem: volume(D) = (1/3) ∮_{∂D} (x−x0)·dS
    Used in Eqs. (1)–(2) to convert volume to a surface integral; exact for the true isosurface.
  • domain assumption The isosurface is well approximated by a plane-triangle triangulation with vertices obtained by linear interpolation on grid lines, and the resulting polyhedron is closed and encloses the region.
    Eqs. (3)–(6); fails near grid degeneracies or when the isosurface intersects domain boundaries in an unresolved way.
  • ad hoc to paper Volume error has the asymptotic form V(Δx) = V_true + C2(Δx)^2 + o((Δx)^2) with the same C2 for the downsampled 2Δx grid.
    Eqs. (7)–(8); this is the key assumption that makes extrapolation (9) valid and is asserted by analogy, not proved.
  • domain assumption Downsampling the fine grid yields a valid coarse discrete representation of the same function at spacing 2Δx.
    Stated after Eq. (9); requires nested grids and ignores possible aliasing of features smaller than 2Δx.
  • ad hoc to paper In the EM test, the volume of the activated region is the volume enclosed by the isosurface and the skin surface z=0, and the numerical method accounts for this boundary.
    §4.3, Eq. (28); the paper does not explain how the boundary cap enters the triangulation integral from §2.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Accurate Computation of Activated Volume in Electromagnetic Heating." pith.science (2026). https://pith.science/paper/AVD45TOO

@misc{pith2026260725994,
  author       = {Pith},
  title        = {Pith review of: Accurate Computation of Activated Volume in Electromagnetic Heating},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/AVD45TOO}},
  note         = {Machine review of arXiv:2607.25994}
}
read the original abstract

In electromagnetic heating and other applications, we need to compute the volume enclosed by an isosurface of the 3D temperature distribution that is numerically represented on a rectangular grid. This situation arises naturally when the temperature distribution is obtained by solving a partial differential equation numerically using a finite difference method. Given the temperature distribution on the 3D grid, the isosurface of a prescribed value is represented approximately by a triangulation, a collection of triangles with vertices on the grid lines. The vertices are determined by a linear interpolation to approximate the locations where the temperature is at the prescribed level. The region enclosed by the isosurface is approximated by that enclosed by the triangulation, which is a set of tetrahedrons. The enclosed volume is approximated by summing those of tetrahedrons. This volume approximation is analogous to approximating a curve using line segments and is limited in accuracy. In this study, we combine extrapolation with the triangulation approximation to develop a more accurate method for computing the volume enclosed by an isosurface of the 3D temperature distribution on a rectangular grid.

Figures

Figures reproduced from arXiv: 2607.25994 by the authors.

Figure 1
Figure 1. Top: region enclosed by the triangulation of an isosurface. Bottom: detailed view [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Left: Triangulation of the unit sphere as an isosurface of function [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Left: Triangulation of the axis-aligned ellipsoid as an isosurface of [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: Left: Triangulation of the rotated ellipsoid as an isosurface of [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: Left: Triangulation of the axis-aligned superellipsoid as an isosurface of [PITH_FULL_IMAGE:figures/full_fig_p010_5.png]
Figure 6
Figure 6. Figure 6: Left: Triangulation of the rotated superellipsoid as an isosurface of [PITH_FULL_IMAGE:figures/full_fig_p012_6.png]
Figure 7
Figure 7. Figure 7: Left: Triangulation of the axis-aligned elliptic ring torus as an isosurface of [PITH_FULL_IMAGE:figures/full_fig_p013_7.png]
Figure 8
Figure 8. Figure 8: Left: Triangulation of the rotated elliptic ring torus as an isosurface of [PITH_FULL_IMAGE:figures/full_fig_p014_8.png]
Figure 9
Figure 9. Figure 9: The coordinate system for skin and incident electromagnetic wave. [PITH_FULL_IMAGE:figures/full_fig_p015_9.png]
Figure 10
Figure 10. Figure 10: Left: Triangulation of the activated skin region as an isosurface of [PITH_FULL_IMAGE:figures/full_fig_p020_10.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references · 5 canonical work pages

  1. [1]

    Begley, A

    Romanenko, S., R. Begley, A. R. Harvey, L. Hool, and V. P. Wallace (2017) The interaction between electromagnetic fields at megahertz, gigahertz and terahertz fre- quencies with cells, tissues and organisms: risks and potentials. Journal of the Royal Society Interface 14: 20170585

    2017

  2. [2]

    Chahat, R

    Zhadobov, M., N. Chahat, R. Sauleau, C. Le Quement, and Y. Le Drean (2011) Millimeter-wave interactions with the human body: state of knowledge and recent 22 advances. International Journal of Microwave and Wireless Technologies 3:237-247

    2011

  3. [3]

    Le Quement, C., C. N. Nicolaz, D. Habauzit, M. Zhadobov, R. Sauleau, and Y. Le Drean (2014) Impact of 60-GHz millimeter waves and corresponding heat effect on endoplasmic reticulum stress sensor gene expression. Bioelectromagnetics 35:444-451

    2014

  4. [4]

    Nelson, D.A., M. T. Nelson, T. J. Walters and P. A. Mason (2000) Skin heating effects of millimeter-wave irradiation – thermal modeling results. IEEE Transactions on Microwave Theory and Techniques 48:2111-2120

    2000

  5. [5]

    Foster, K. R., H. Zhang and J. M. Osepchuk (2010) Thermal response of tissues to mil- limeter waves: implications for setting exposure guidelines. Health physics 99(6):806- 810

    2010

  6. [6]

    and O’Connor, K

    Cazares, S.M., Snyder, J.A., Belanich, J., Biddle, J.C., Buytendyk, A.M., Teng, S.H.M. and O’Connor, K. (2019) Active Denial Technology Computational Human Effects End-to-End Hypermodel for Effectiveness (ADT CHEETEH-E).Human Factors and Mechanical Engineering for Defense and Safety,3, Article No. 13. doi:10.1007/s41314- 019-0023-7

    work page doi:10.1007/s41314- 2019

  7. [7]

    and Oberhammer, J

    Topfer, F. and Oberhammer, J. (2015) Millimeter-wave Tissue Diagnosis: The most Promising Fields for Medical Applications.IEEE Microwave Magazine,16(4), 97-113

    2015

  8. [8]

    (2020) A

    Srivastave, A. (2020) A. Srivastave, Energy efficient transmission trends towards future green cognitive radio networks (5G): Progress, taxonomy and open challenges, Journal of Network and Computer Applications 168: 102760

    2020

Show all 20 references
  1. [9]

    Sheen, D.M., , D. L. McMakin and T. E. Hall (2007) Detection of explosives by millimeter-wave imaging, in Counterterrorist Detection Techniques of Explosives, edited by J. Yinon. 237-278

    2007

  2. [10]

    (2020) Robot Systems for Rail Transit Applications

    Liu, H. (2020) Robot Systems for Rail Transit Applications. Elsevier

    2020

  3. [11]

    Burgei, W.A

    Wang, H., W. Burgei, W.A. and H. Zhou (2020) A concise model and analysis for heat- induced withdrawal reflex caused by millimeter wave radiation. American Journal of Operations Research, 10, 31-81. 23

    2020

  4. [12]

    Wang, H. ,W. Burgei, W. and H. Zhou (2020) Non-Dimensional Analysis of Thermal Effect on Skin Exposure to an Electromagnetic Beam. American Journal of Operations Research, 10, 147-162

    2020

  5. [13]

    H. Wang, S. E. Foley, and H. Zhou (2025) Asymptotic Solution for Skin Heating by an Electromagnetic Beam at an Incident Angle,Electronics, vol. 14, no. 15, p. 3061

    2025

  6. [14]

    and Crouch, S.R

    Skoog, D.A., Holler, F.J. and Crouch, S.R. (2017) Principal of Instrumental Analysis. 7th Edition, Sunder College Publisher, New York

    2017

  7. [15]

    E., and Zhou, H

    Jaime-Yepez, U., Wang, H., Foley, S. E., and Zhou, H. (2024) Asymptotic solution of electromagnetic heating of skin tissue with lateral heat conduction.J Eng Math,147,

    2024

  8. [16]

    doi: doi.org/10.1007/s10665-024-10390-y

    doi

  9. [17]

    (2000) Heating and Pain Sensation Produced in Human Skin by Millimeter Waves: Comparison to a Sim- ple Thermal Model.Health Physics,78(3): 259–267

    Walters T.J., Blick D.W., Johnson L.R., Adair E.R., Foster K.R. (2000) Heating and Pain Sensation Produced in Human Skin by Millimeter Waves: Comparison to a Sim- ple Thermal Model.Health Physics,78(3): 259–267. doi:10.1097/00004032-200003000- 00003

    2000 doi

  10. [18]

    and Zhou, H

    Wang, H., Burgei, W.A. and Zhou, H. (2020) Analytical solution of one-dimensional Pennes’ bioheat equation.Open Physics,18(1), 1084-1092. doi: 10.1515/phys-2020- 0197

    2020 doi

  11. [19]

    A. S. Rodr ´ ıguez-P´ erez, H. E. Gilardi-Vel´ azquez, and S. E. Vel´ azquez-P´ erez (2025) Analysis and Simulation of Dynamic Heat Transfer and Thermal Distribution in Burns with Multilayer Models Using Finite Volumes,Dynamics, vol. 5, no. 4, Art. 41. doi:10.3390/dynamics5040041

    2025 doi

  12. [20]

    Joshi, F

    A. Joshi, F. Wang, Z. Kang, B. Yang, and D. Zhao (2022) A three-dimensional thermoregulatory model for predicting human thermophysiological responses in var- ious thermal environments,Building and Environment, vol. 207, Art. 108506. doi:10.1016/j.buildenv.2021.108506. 24

    2022

Pith tools

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