Exploring Heat Exchanges with the Calorimetry Simulator -- SimuF\'isica

Alex B. Siqueira; Marco P. M. de Souza

REVIEW 3 major objections 4 minor 8 references

Exploring Heat Exchanges with the Calorimetry Simulator -- SimuF\'isica

T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read The Calorimetry–SimuFísica simulator claims to reproduce standard heat-exchange equations, including phase changes, so students can check textbook calculations against it.

desk verdict A useful but flawed simulator paper: the printed model equations cannot reproduce the paper's own numbers. read the letter →

arxiv 2506.03289 v1 pith:TJL334KD submitted 2025-06-03 physics.ed-ph

Marco P. M. de Souza , Alex B. Siqueira This is my paper

classification physics.ed-ph

keywords calorimetryphysicseducationinteractivesimulationheatexchangeNewton'slawofcoolinglatentthermalequilibriumspecific

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 introduces an interactive calorimetry simulator and argues that its physics engine matches the equations of elementary calorimetry, so students can use it to verify hand calculations. In Heating mode the simulator raises a liquid's temperature linearly until boiling, then holds the temperature fixed while vaporizing liquid at a rate set by flame power and latent heat. In Thermal Equilibrium mode it models a hot metallic sphere cooling in a colder liquid via Newton's law of cooling, with temperatures approaching a common equilibrium value over time. The paper validates the tool by solving two standard problems—water heated to evaporation, and a steel sphere dropped into water—and reports that the simulator reproduces the analytic answers.

What carries the argument

The object that carries the argument is the pair of ordinary differential equations (7), Newton's law of cooling written as $c_{\rm liq}\,dT_{\rm liq}/dt = -k(T_{\rm liq}-T_{\rm esf})$ and $c_{\rm esf}\,dT_{\rm esf}/dt = k(T_{\rm liq}-T_{\rm esf})$, with $k=1$ cal/(g·°C·s) in the simulator. The analytic solutions are the exponential relaxation formulas (8), which make both temperatures approach a common limit; once the liquid reaches boiling, the piecewise equations (9) freeze the liquid temperature at the boiling point and spend the sphere's energy on vaporization. These formulas are what let the simulator display real-time temperatures, absorbed heat, and evaporated mass that can be compared directly with textbook calculations.

What would settle it

Run the Thermal Equilibrium mode with Problem 2's parameters but double the sphere's mass to 400 g at the same 600°C. Mass-weighted calorimetry predicts a final temperature near 56°C, while a literal reading of Eq. (7) predicts the same 77°C regardless of mass; the displayed final temperature distinguishes the two cases.

Watch

Extended reading notes

Core claim

The paper's central claim is that the Calorimetry–SimuFísica simulator embodies the standard calorimetry equations and therefore reproduces heat exchange faithfully, with or without phase change. In Heating mode the model is a linear temperature rise $T(t)=T_0+Pt/(m_0c)$ up to the boiling point, then a constant-temperature vaporization stage in which the evaporated mass is $P(t-t_E)/L$. In Thermal Equilibrium mode the model is the coupled Newton-cooling system (7), whose solutions (8) are exponentials relaxing to a common temperature, replaced by the piecewise boiling branch (9) once the liquid reaches its boiling point. The paper tests this claim on two worked problems and reports that the simulator's outputs—100°C with 200 g evaporated in Problem 1, and 41°C with 11.2 kcal transferred in Problem 2—match the analytic calculations.

Load-bearing premise

The whole argument depends on the simulator's internal code using each body's mass times its specific heat in the heat-exchange equations, even though the printed Eq. (7) leaves the masses out; if the code followed the printed equations literally, the paper's own Problem 2 would end near 77°C, not the reported 41°C.

Editorial extensions

If this is right

In Heating mode, a liquid at constant flame power shows a linear temperature ramp up to its boiling point, then a flat temperature plateau while the absorbed heat climbs and the liquid level falls, letting students read sensible heat and latent heat from the graph.
In Thermal Equilibrium mode, students watch the sphere cool rapidly while the water warms until both settle at the same temperature, with the final temperature given by the energy-conservation weighted average once masses are included.
The reported agreement—11.2 kcal transferred versus 11.18 kcal computed, and 41°C versus 40.97°C—means the simulator can serve as a quantitative check on hand calculations, not just a qualitative animation.
Teachers can use the predefined example configurations to set up classroom comparisons between analytical solutions and simulated graph output without losing time on parameter entry.
The simulator gives a safe, repeatable environment for exploring boiling and vaporization, processes that are difficult to control quantitatively in a school laboratory.

Reading between the lines

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

The paper leaves implicit that the approach-to-equilibrium time scale is set by the transfer constant $k$, while the final equilibrium temperature is independent of $k$; teachers could use $k$ as a free parameter to explore exponential relaxation without changing the endpoint.
A direct empirical test, which the paper lists only as future work, is to run the same two problems with real beakers, water, and a heated metal sphere; agreement there would support using the simulator in place of laboratory equipment in schools.
The printed transcription of the Thermal Equilibrium equations needs a correction: as written, Eq. (7) mixes specific heat with heat capacity, and a reader who implements it verbatim would not reproduce the reported 41°C final temperature unless masses are inserted into the denominators.
One natural extension is to include the container's heat capacity and heat losses, both neglected in the paper; classroom experiments with real containers would then reveal systematic deviations, giving students a concrete lesson in idealization.

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.

Desk Editor's Note

read the letter

Short version: this is a straightforward description of a new educational calorimetry simulator, and the worked examples are mostly right. But the central model equations in Section 2 are wrong as written: they omit masses, and with the paper's own Problem 2 values they give 77°C, not the 41°C reported. That's not cosmetic; it's the paper's main consistency claim.

What's genuinely new: the Calorimetry SimuFísica app itself — free, offline, multilingual, with two modes (heating and thermal equilibrium) and save/load. The authors walk through the interface and give two classroom-style problems with analytic solutions. The analytic calculations in Problems 1 and 2 are correct, and if the code implements the mass-weighted energy balance (which it presumably does, since the simulator shows the right numbers), then the tool does what it claims. The paper also clearly states that no pedagogical impact study has been done yet, and points to that as future work — that's honest.

The soft spots are real. Equations (7) and (8) define the thermal-equilibrium dynamics using specific heats without masses. As printed, the equilibrium is a weighted average of the specific heats, not the heat capacities. For Problem 2 that gives about 77°C, while the paper's own solution gives 40.97°C. The reported 41°C is recovered only if the code actually solves m_i c_i dT_i/dt = ... — i.e., with heat capacities. The paper never says that. So the claim that the simulator 'faithfully reproduces' the process rests on an implementation detail that is nowhere documented. This is the load-bearing issue. There's also a sign/expansion typo in Eq. (19): 12000 - 20T_f = 700T_f - 1750 should give 13750 = 720T_f, not the printed rearrangement, though the final answer happens to be right. Regardless, it's a typo.

The lack of an independent benchmark (e.g., comparing with PhET or with real lab measurements) is a limitation but not a fatal one; for a teaching tool, that can be a future step. The citation pattern is fine: self-citations are to their own platform, which is relevant.

Who's it for? Physics teachers and curriculum designers, especially in Brazil; anyone looking for a free offline calorimetry simulation.

Bottom line: this paper deserves a serious referee, not a desk rejection, but it needs a revised version fixing Eqs. (7)–(8) to include masses (or explicitly stating the code uses heat capacities), and cleaning up Eq. (19). After that, it's a serviceable contribution to physics education. I wouldn't cite it in my own work, but I'd send it to review.

Referee Report

3 major / 4 minor

Summary. This paper presents the Calorimetry simulator from the SimuFísica platform, an interactive tool for teaching heat exchange at high-school and introductory university level. The authors describe the interface, the Heating and Thermal Equilibrium modes, and the physical equations underlying them: sensible heat, latent heat, Newton's law of cooling, and a set of ordinary differential equations for the thermal contact between a sphere and a liquid. Two classroom problems are analyzed: heating 500 g of water with 148 kcal and finding the evaporated mass, and placing a 200 g steel sphere at 600°C into 700 mL of water at 25°C. In both cases the simulator outputs are shown to agree with the analytical solutions, and the paper claims that this confirms consistency between theory and simulation.

Significance. The educational context is relevant, and the paper gives explicit online links to the exact simulator configurations, which makes the two examples reproducible by readers. The analytical calculations in Problems 1 and 2 are correct, and the reported simulator values (100°C and 200 g evaporated in Problem 1; about 41°C and 11.2 kcal in Problem 2) match the textbook results. However, the central physical model as printed in Eqs. (7)-(8) is incorrect because it omits the masses, so the stated equations cannot produce the reported Problem 2 result. In addition, the comparison between theory and simulation is partly circular, since the simulator is constructed from the same model it is being checked against. These issues are fixable, but they are load-bearing for the paper's claim that the simulator 'faithfully reproduce[s] the heat exchange process.'

major comments (3)

[Section 2, Eqs. (7)-(8)] Equations (7) are written as c_liq dT_liq/dt = -k(T_liq - T_esf) and c_esf dT_esf/dt = k(T_liq - T_esf), where c_liq and c_esf are stated to be specific heats. This system has equilibrium temperature T_eq = (c_esf T_esf0 + c_liq T_liq0)/(c_esf + c_liq). Substituting the paper's own Problem 2 values (c_esf = 0.1 cal/g°C, c_liq = 1 cal/g°C, T_esf0 = 600°C, T_liq0 = 25°C) gives T_eq ≈ 77°C, not the 41°C obtained from the mass-weighted Eq. (17) and shown in Fig. 3. The reported result is recovered only if the equations use heat capacities C_i = m_i c_i, with C_esf = 20 cal/°C and C_liq = 700 cal/°C, yielding (20×600 + 700×25)/720 ≈ 40.97°C. The paper must correct Eqs. (7)-(8) to include the masses explicitly, or clearly state that c_liq and c_esf denote heat capacities in that section; as printed, the equations cannot describe the simulation data they claim to govern.
[Section 3, Problem 2 comparison] The agreement between the analytical solution (20) and the simulator output in Fig. 3 is partly by construction. The simulator was built from the same physical model (Eqs. (1)-(9)) used in the analytical treatment, so the comparison verifies internal consistency between the code and its design equations, not agreement with an independent benchmark such as laboratory data or a distinct analytical derivation. The paper should explicitly acknowledge this circularity and soften the claim that the simulator is 'consistent with classical mathematical models'; at minimum, the corrected Eqs. (7)-(8) should be checked against the independent energy-balance solution (17).
[Section 2, Eq. (9c)] The vaporized-mass formula (9c) is asserted without a validity condition. The numerator subtracts the heat needed to bring the liquid to its boiling point from the heat released by the sphere cooling to T_E; if the sphere's available heat is smaller than that required heat, the expression becomes negative and does not represent a vaporized mass. The paper should state the domain of applicability, e.g., that Eq. (9c) holds only when m_esf c_esf (T_esf - T_E) > m_liq c_liq (T_E - T_liq0), and should specify what the simulator does when this condition fails. This matters because the paper claims the Thermal Equilibrium mode handles phase change.

minor comments (4)

[Section 2, heading] The subsection heading 'B. Fundamentals' appears to be a leftover formatting artifact; it should be a numbered section heading such as '2.2 Fundamentals'.
[Section 2, Eq. (7)] The unit quoted for k, 1 cal/(g·°C·s), makes Eqs. (7) dimensionally consistent as per-gram energy rates but not as an energy balance for two bodies of different masses; the text should explain this distinction and specify the actual numerical value used in the simulator.
[Section 3, Figures 2 and 3] The screenshots in Figs. 2 and 3 are small; the temperature and heat axes are hard to read. Including enlarged graph panels or a data table would make the quantitative comparison easier for readers to verify.
[References] Reference [7] is a quantum-physics paper about a Bloch equation generator; it appears to be cited only as another example of the SimuFísica platform, but this connection is not explained in the text.

Circularity Check

1 steps flagged · score 4.0 of 10

Consistency claim is a self-consistency check: the simulator agrees with analytic solutions of the same model, so the agreement is partly by construction; an additional mass omission in Eq. (7) severs the stated model from the reported output.

self definitional [Section 2 (Fundamentals) and Section 3 (Problem 2), Eqs. (7)-(9), (20)-(21)]
"Equations (8) and (9) describe the system dynamics in Thermal Equilibrium mode and allow the simulator to faithfully reproduce the heat exchange process, with and without phase change. ... When configured in Thermal Equilibrium mode, with the sphere at 600°C and water at 25°C, the simulator shows the sphere cooling and the water heating until both reach about 41°C, as shown in Fig. 3."

The simulator is a numerical implementation of the same calorimetric model (specific heat, energy conservation, Newton's law of cooling) that is used to derive the analytic solution in Eq. (20). The paper claims Eqs. (8)-(9) describe the simulator dynamics and then compares the simulator output to the analytic solution of that same model. Hence the reported agreement is a self-consistency check: the simulator reproduces the behavior of the model from which it was built, not an independent prediction validated against external data. Because no independent benchmark is used, the observed match is guaranteed by construction up to numerical error, making the 'faithful reproduce' claim circular in its validation logic.

full rationale

The paper's central contribution is an educational simulator and its comparison with textbook calorimetry calculations. The consistency between the simulator and the analytic solutions in Problems 1 and 2 is presented as a highlight, but it is a verification that the numerical implementation matches the same physical model used in the analytic examples; both sides derive from the same equations for sensible heat, latent heat, and Newton's law of cooling. This is a self-consistency check, not an independent validation against experiment or an external benchmark (the paper itself lists laboratory comparison as future work). No parameters are fitted, and the analytic solutions are standard textbook results, so the circularity is moderate rather than severe. A separate, non-circular defect should be flagged: Eq. (7) omits masses, and with the stated constants the equilibrium of Eq. (8) is ~77°C for Problem 2, not the reported ~41°C; the reported agreement therefore requires an undocumented mass-weighted implementation, meaning the printed model cannot support the fidelity claim. This is a correctness/reporting issue, but it compounds the self-definitional character of the validation: the simulator's actual behavior is defined by an unstated model, and the analytic solution is the same unstated model, so the agreement remains by construction rather than by independent test.

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

The central model introduces one free parameter, k, and relies on standard textbook simplifications. No new physical entities are proposed. The key hidden assumption is that the simulator code implements mass-weighted heat capacities, which the written equations do not reflect.

free parameters (1)

k (heat transfer constant) = 1 cal/(g·°C·s) (adopted)
Chosen value in the simulator for the rate of heat exchange in Newton's law of cooling; not fitted to data and not validated against experiment. It affects the time evolution but not the final equilibrium temperature in the two examples.

assumptions (4)

domain assumption Thermal losses and the heat capacity of the container are neglected in both the analytic solutions and the simulator model.
Stated in Section 2 (B) before Eq (1) and repeated in Problems 1 and 2; simplification used to make the model and the worked examples tractable.
domain assumption Heat exchange between liquid and sphere follows Newton's law of cooling with a constant transfer coefficient k.
Assumed in Section 2, Eqs (7); standard cooling model but a modeling choice; k is not derived from first principles.
domain assumption Specific heats and latent heat are constant over the temperature ranges considered.
Implied throughout the worked examples; typical for introductory calorimetry exercises.
ad hoc to paper The simulator's internal numerical implementation correctly includes the masses in the energy balance, despite Eq (7) omitting them.
The paper's own Eq (7) does not contain masses, yet Problem 2's simulator output matches the mass-weighted calculation; the paper assumes, without stating, that the code uses mass times specific heat rather than bare specific heat.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Exploring Heat Exchanges with the Calorimetry Simulator -- SimuF\'isica." pith.science (2026). https://pith.science/paper/TJL334KD

@misc{pith2026250603289,
  author       = {Pith},
  title        = {Pith review of: Exploring Heat Exchanges with the Calorimetry Simulator -- SimuF\'isica},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TJL334KD}},
  note         = {Machine review of arXiv:2506.03289}
}

read the original abstract

This article presents the Calorimetry - SimuF\'isica simulator, an interactive computational tool designed for teaching heat exchange processes. The simulator enables dynamic and audiovisual exploration of phenomena such as the heating of liquids and the establishment of thermal equilibrium between bodies at different temperatures. We describe its interface and the underlying physics, based on the equations of specific heat, latent heat of vaporization, and Newton's law of cooling. Two educational application examples are analyzed, in which the simulation results are compared with calculations typically performed at the high school and early undergraduate levels, highlighting the consistency between theory and simulation and the pedagogical potential of the tool.

Figures

Figures reproduced from arXiv: 2506.03289 by the authors.

**Figure 2.** Simulator configured for Problem 1, in Heating mode and online smartphone version, showing [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

**Figure 3.** Simulator configured for Problem 2, in Thermal Equilibrium mode (desktop version). Final [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

8 extracted references · 8 canonical work pages

[1]

Rutten, W

N. Rutten, W. R. van Joolingen, and J. T. van der Veen,The learning effects of computer simulations in science education, Comput. Educ.58, 136 (2012). doi:10.1016/j.compedu.2011.07.017

work page doi:10.1016/j.compedu.2011.07.017 2012
[2]

D. H. Jonassen,Computers as Mindtools for Schools(2nd ed.), Merrill, New Jersey, USA, 2000

work page 2000
[3]

H. J. Banda and J. Nzabahimana,Effect of integrating physics education technology simulations on students’ conceptual understanding in physics: A review of literature, Phys. Rev. Phys. Educ. Res. 17, 023108 (2021). doi:10.1103/PhysRevPhysEducRes.17.023108

work page doi:10.1103/physrevphyseducres.17.023108 2021
[4]

M. P. M. de Souza, C. M. Oliveira, and R. P. P. Ara´ ujo,O simulador conserva¸ c˜ ao de energia mecˆ anica – SimuF ´ ısica®, A F ´ ıs. na Escola22, 240173 (2024). doi:10.59727/fne.v22i1.173. English version (arXiv):https://doi.org/10.48550/arXiv.2410.17951

work page Pith review arXiv doi:10.48550/arxiv.2410.17951 2024
[5]

M. P. M. de Souza, S. P. Oliveira, and V. L. Luiz,Motor el´ etrico – SimuF ´ ısica ®: um aplicativo para o ensino de eletromagnetismo, Rev. Bras. Ens. F ´ ıs.46, e20230219 (2024). doi:10.1590/1806- 9126-RBEF-2023-0219. English version (arXiv):https://doi.org/10.48550/arXiv.2410.18721

work page Pith review arXiv doi:10.48550/arxiv.2410.18721 2024
[6]

C. M. Oliveira,Ensino da tem´ atica energia com trˆ es aplicativos da plataforma SimuF ´ ısica no Novo Ensino M´ edio, Master’s thesis, Universidade Federal de Rondˆ onia – UNIR, 2024

work page 2024
[7]

M. P. M. de Souza, G. H. H. Pav˜ ao, A. A. C. de Almeida, and S. S. Vianna,Bloch Equation Generator – SimuF ´ ısica, arXiv:2506.01108 [quant-ph].https://doi.org/10.48550/arXiv.2506.01108

work page Pith review arXiv doi:10.48550/arxiv.2506.01108
[8]

W. P. Silva, J. W. Precker, C. M. D. P. S. Silva, D. D. P. S. Silva, and C. D. P. S. Silva,Medida de calor espec ´ ıfico e lei de resfriamento de Newton: um refinamento na an´ alise dos dados experimentais, Rev. Bras. Ens. F ´ ıs.25(4), 392 (2003). doi:10.1590/S0102-47442003000400015 7

work page doi:10.1590/s0102-47442003000400015 2003

Pith tools

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

[1] [1]

Rutten, W

N. Rutten, W. R. van Joolingen, and J. T. van der Veen,The learning effects of computer simulations in science education, Comput. Educ.58, 136 (2012). doi:10.1016/j.compedu.2011.07.017

work page doi:10.1016/j.compedu.2011.07.017 2012

[2] [2]

D. H. Jonassen,Computers as Mindtools for Schools(2nd ed.), Merrill, New Jersey, USA, 2000

work page 2000

[3] [3]

H. J. Banda and J. Nzabahimana,Effect of integrating physics education technology simulations on students’ conceptual understanding in physics: A review of literature, Phys. Rev. Phys. Educ. Res. 17, 023108 (2021). doi:10.1103/PhysRevPhysEducRes.17.023108

work page doi:10.1103/physrevphyseducres.17.023108 2021

[4] [4]

M. P. M. de Souza, C. M. Oliveira, and R. P. P. Ara´ ujo,O simulador conserva¸ c˜ ao de energia mecˆ anica – SimuF ´ ısica®, A F ´ ıs. na Escola22, 240173 (2024). doi:10.59727/fne.v22i1.173. English version (arXiv):https://doi.org/10.48550/arXiv.2410.17951

work page Pith review arXiv doi:10.48550/arxiv.2410.17951 2024

[5] [5]

M. P. M. de Souza, S. P. Oliveira, and V. L. Luiz,Motor el´ etrico – SimuF ´ ısica ®: um aplicativo para o ensino de eletromagnetismo, Rev. Bras. Ens. F ´ ıs.46, e20230219 (2024). doi:10.1590/1806- 9126-RBEF-2023-0219. English version (arXiv):https://doi.org/10.48550/arXiv.2410.18721

work page Pith review arXiv doi:10.48550/arxiv.2410.18721 2024

[6] [6]

C. M. Oliveira,Ensino da tem´ atica energia com trˆ es aplicativos da plataforma SimuF ´ ısica no Novo Ensino M´ edio, Master’s thesis, Universidade Federal de Rondˆ onia – UNIR, 2024

work page 2024

[7] [7]

M. P. M. de Souza, G. H. H. Pav˜ ao, A. A. C. de Almeida, and S. S. Vianna,Bloch Equation Generator – SimuF ´ ısica, arXiv:2506.01108 [quant-ph].https://doi.org/10.48550/arXiv.2506.01108

work page Pith review arXiv doi:10.48550/arxiv.2506.01108

[8] [8]

W. P. Silva, J. W. Precker, C. M. D. P. S. Silva, D. D. P. S. Silva, and C. D. P. S. Silva,Medida de calor espec ´ ıfico e lei de resfriamento de Newton: um refinamento na an´ alise dos dados experimentais, Rev. Bras. Ens. F ´ ıs.25(4), 392 (2003). doi:10.1590/S0102-47442003000400015 7

work page doi:10.1590/s0102-47442003000400015 2003