PG-EAM - Graduate Program in Aeronautical and Mechanical Engineering
PT EN
Rene F. B. Gonçalves

Rene F. B. Gonçalves

8
h-index
210
Citations
52
Articles

Research Lines

  • Chemical propulsion
  • Reactive molecular dynamics
  • Energetic materials
Last Update: 2026-08-17

Publications (52)

52 publications
Article 2026

E-waste in the circular economy: Molecular and experimental insights into cement incorporation of waste printed circuit board

Gonçalves, Rene Francisco Boschi , Yao, Zhitong , Rocco, José Atílio Fritz Fidel

Circular Economy , vol. 5 (1)
Citations: 1
Show abstract

© 2026 The Author(s).The escalating generation of electronic waste (e-waste) poses critical environmental and health challenges, with waste printed circuit boards (WPCBs, 4 wt.%–6 wt.% of e-waste) representing a particularly hazardous fraction due to their complex composition of metals, polymers, and brominated flame retardants. This study explores the valorization of WPCBs as a partial cement replacement to both mitigate disposal risks and enhance cementitious performance. Portland cement pastes incorporating 0 wt.%, 5 wt.%, 10 wt.%, and 15 wt.% WPCBs were examined through isothermal calorimetry, Fourier transform infrared (FT-IR) spectroscopy, and reactive molecular dynamics (RMD) simulations using the reactive force field. Calorimetry revealed that 10% WPCBs yielded the highest heat release (≈25 °C peak, 10% higher than reference), indicating optimal hydration kinetics, while 15 wt.% WPCBs reduced exothermicity by 8% due to inert filler effects and bromine-induced retardation. FT-IR revealed intensified Si–O stretching (953 cm−1) and free O–H (3640 cm−1) at 10 wt.%, confirming enhanced C–S–H and CH formation. RMD simulations provided molecular-level insights, showing denser C–S–H networks and the lowest total energy (−6.82 × 105 kcal/mol) for 10 wt.% WPCBs, whereas 15 wt.% loading disrupted network formation. The integrated findings suggested that WPCBs can be safely immobilized within cement matrices, with 10 wt.% substitution offering a balance between performance gains and hydration integrity. This work advances the understanding of WPCB–cement interactions and supports the development of sustainable, circular-economy construction materials.

Article 2026

Castor oil as a sustainable alternative to HTPB binder in polyurethane matrices for composite solid propellants: viscosity, curing kinetics, and hardness evaluation

Júnior, Jorge Fernando Leite Monteiro , Mendonça, Fausto Batista , Rocco, José Atílio Fritz Fidel , Pinto, Juliano Ribeiro Aguiar , Gonçalves, Rene Francisco Boschi , Diniz, Milton Faria

Journal of Energetic Materials
Show abstract

© 2026 Taylor & Francis Group, LLC.Increasing environmental regulations are driving the development of green solid propellants with reduced toxicity and lower life-cycle environmental impact. The binder plays a critical role in solid propellants, with hydroxyl-terminated polybutadiene (HTPB) being the most widely used despite its petroleum origin, environmental concerns, and restricted availability. This study evaluates castor oil (CO), a renewable, biodegradable, and low-toxicity resource, as a sustainable alternative binder. A CO-based polyurethane prepolymer containing 30% (w/w) toluene diisocyanate was synthesized and characterized by Fourier transform infrared spectroscopy (FT-IR), rheology, particle distribution, hardness tests, and scanning electron microscopy (SEM). The results showed 20% faster curing and hardness comparable to HTPB, highlighting the potential of CO as a viable and sustainable binder for solid propellant applications. Additionally, SEM analysis comparing the top and bottom regions of two CO-based samples produced by different synthesis using ammonium perchlorate (AP) as filler demonstrated a more homogeneous AP distribution in the sample #2 (a prepolymer it was first synthesized by reacting CO with 30% of the TDI needed for complete curing).

Article 2026

Thiadiazole-furazan compounds: decomposition mechanism study by RMD and DFT

Kuznetsov, Aleksey , Gonçalves, Rene F.B.

Structural Chemistry
Show abstract

© The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature 2026.This study uses Reactive Molecular Dynamics (RMD) and Density Functional Theory (DFT) to investigate the properties and decomposition mechanisms of two thiadiazole-furazan compounds, nitrothidiazolefurazan (NTF) and nitrobithidialozolefurazan (NBTF), as potential high-energy-density materials. RMD simulations revealed their thermal stability and initial decomposition pathways, while DFT provided insights into structures, electronic properties, and bond orders, as well as mechanistic support. Results indicate the suitability for extreme conditions in defense or space applications, advancing understanding and design of new energetic materials.

Article 2025

Thermal decomposition of spent lithium-ion batteries pouch: Investigating kinetic and thermodynamic compensation effects

Jiang, Jingjing , Yao, Zhitong , Tong, Jiayao , Cui, Jiuzhuo , Kumar, Akash , Gonçalves, Rene F.B. , Reinmöller, Markus , Sangaré, Diakaridia , Manić, Nebojša , Liu, Jie , Bertelsen, Michael

Chemical Engineering Science , vol. 316
Citations: 7
Show abstract

© 2025 Elsevier LtdA deeper understanding of the pyrolysis process for main and additional components in spent lithium-ion batteries (LIBs) could provide valuable insights for optimizing their recycling processes. This study examined the thermal behavior, kinetics, thermodynamics, and product evolution during the pyrolysis of laminate pouch primarily composed of polypropylene and polyamide. The kinetic compensation effect (KCE) and thermodynamic compensation effect (TCE) were also probed to provide a comprehensive understanding of the conversion. The degradation process was divided into three stages, with total mass loss ranging from 31.14 to 40.28 % and peak temperatures between 419 and 472 °C. The average activation energy was determined to be 118.06 kJ mol−1, with specific values of 99.25, 119.06 and 139.31 kJ mol−1 within conversion rate of 0.10–0.30, 0.35–0.75 and 0.80–0.95, respectively. The pouch conversion followed D1 diffusion mechanism. The KCE was confirmed and reconstructed fα=α0.45316(2α1.55)-1 displayed an excellent fit. Thermodynamic analysis implied that this conversion process was endothermic and non-spontaneous. Enthalpy and entropy relationship demonstrated the existence of TEC with compensation temperature (Tcomp) and experimental temperature (Texp) of 676.20 K and 693.23 K, respectively. In addition, free energy of compensation (ΔGcomp) was found to be 164.51 kJ mol−1, in agreement with experimental binding free energy (ΔGexp) range of 166.48–170.65 kJ mol−1, further confirmed the validity of the adopted mechanism.

Article 2025

Kinetics, thermodynamics, and product formation during pyrolysis conversion of protection board from spent lithium-ion batteries

Tong, Jiayao , Yao, Zhitong , Jiang, Jingjing , Cui, Jiuzhuo , Kumar, Akash , Gonçalves, Rene F.B. , Reinmöller, Markus , Vegliò, Francesco , Romano, Pietro , Liu, Jie , Jin, Meiqing , Bertelsen, Michael

Journal of Energy Storage , vol. 128
Citations: 8
Show abstract

© 2025 Elsevier LtdThermal treatment of spent lithium-ion batteries offers the benefits of decomposing organic components while concentrating valuable metals. This work investigated the kinetics, thermodynamics, and evolved products during the protection board pyrolysis under N2 and CO2 atmospheres. The degradation process was divided into stages of below 400 °C, 400–700 °C, and 700–900 °C. Peak temperatures at the maximum mass loss rate were observed at 359–399 °C in N₂ and 367–393 °C in CO₂. The primary products evolved from phenolics into ketones and acids, and eventually into alkanes. Brominated products such as bromomethane and 1-bromobutane were also detected, indicating the requirement of debromination to improve the usability of pyrolysis products. The average activation energies were determined to be 218.33 kJ/mol in N₂ and 308.91 kJ/mol in CO₂. D4 and D1 reaction mechanisms were found to best describe the pyrolysis process in two atmospheres. Positive values of ΔH and ΔG indicated the endothermic and non-spontaneous characteristics. The difference between ΔH and Ea values ranged from 5.26 to 7.70 kJ/mol in N₂ and 5.12–7.48 kJ/mol in CO₂, indicating a high possibility of overcoming the potential energy barrier.

Article 2025

A Constant Trinitrotoluene Equivalence Fit for Blast Wave Position Versus Time Data

Amorim, Caio Barbosa , Augusto, Anselmo da Silva , Gonçalves, Rene Francisco Boschi

Propellants Explosives Pyrotechnics , vol. 50 (7) , pp. 44-54
Citations: 1
Show abstract

© 2025 The Author(s). Propellants, Explosives, Pyrotechnics published by Wiley-VCH GmbH.Accurately evaluating accidental or intentional detonation scenarios is essential to ensure their intended effectiveness and/or protect personnel and structures. These evaluations often rely on estimating critical blast effect parameters through reference models, which adapt to different scenarios via key constants. The incident peak overpressure, the sudden pressure increase upon blast wave arrival, is a crucial parameter directly associated with wave-induced damage and demands precise calculation. An effective experimental approach involves tracking shock wave positions over time via calibrated high-speed footage. Processing these recordings and fitting parametric models to the data enables low-uncertainty peak overpressure estimation. This study presents a novel method utilizing constant trinitrotoluene (TNT) equivalence, a key parameter that quantifies the mass conversion of an explosive into its TNT equivalent on the basis of blast effects. The new model was compared with existing models in the literature in terms of peak overpressure and time of arrival estimation. These comparisons were made against direct measurements obtained from pressure sensors and high-speed recordings of open-air detonation tests involving center-initiated spherical Composition B explosions. The results indicate that the new model aligns more closely with experimental data than previously established models.

Article 2025

Composition optimization of a hypergolic green propellant based on monoethanolamine, n-butanol and 90% hydrogen peroxide

Mendoza, Paull C.Acosta , Gonçalves, Rene F.B. , Gouvêa, Leonardo Henrique , Pereira, Luís Gustavo Ferroni

Acta Astronautica , vol. 229 , pp. 140-148
Citations: 5
Show abstract

© 2025 IAAThe design of satellite attitude-control thrusters depends on a trade-off between minimum impulse bit and specific impulse, where the width of pulse maneuvers relies on the combination of delays in the hydraulic system (feed tubes and valves) and the ignition delay time of the propellant used. The most well-established propellants in this context are hydrazine derivatives and nitrogen tetroxide. However, their high toxicity makes satellite integration costly and environmentally hazardous. To replace these propellants, research is focused on developing new hypergolic green propellants, most of which use high-concentration hydrogen peroxide as an oxidizer. In this study, the hypergolic reaction between a blend of n-butanol and monoethanolamine and hydrogen peroxide was catalyzed using copper nitrate trihydrate. The central composite design method was applied to optimize fuel composition using 90% hydrogen peroxide as the oxidizer. The optimization yielded two key outcomes: for ignition delay time (31.5% n-butanol, 60% monoethanolamine, and 8.5% copper nitrate, resulting in an ignition delay time of 21.5 ms with a standard deviation of ±1.30 ms and a systematic error of ±0.4), and for theoretical specific impulse (36% n-butanol, 60% monoethanolamine, and 4% copper nitrate, with an ignition delay time of 26 ±0.4 ms). For the ignition delay time optimization, an oxidizer-fuel ratio of 4 was selected using CEA NASA software to achieve a maximum theoretical specific impulse of 170.64 s, while for specific impulse optimization, a ratio of 4.4 was chosen, resulting in a specific impulse of 171.58 s. Although the maximum theoretical specific impulse of the proposed green propellant pair does not present an advantage if compared to traditional hypergolic propellants, it offers a competitive advantage in terms of density-specific impulse, with the highest value achieved in the ignition delay time optimization, where the density-specific impulse of the system reached 267.5 gs/cm3. Furthermore, the addition of n-butanol effectively reduced fuel viscosity, enhanced density-specific impulse, increased specific impulse, and improved ignition delay time response with 90% hydrogen peroxide compared to pure monoethanolamine formulations for a specific chamber and nozzle configuration. These findings highlight the potential of this green propellant system to enhance performance and efficiency in aerospace applications.

Article 2025

Recycling of end-of-life solar panels: Focusing on the pyrolysis conversion of back sheet from a micro perspective

Yao, Zhitong , Tong, Jiayao , Gonçalves, Rene F.B. , Kumar, Akash , Manić, Nebojša , Vegliò, Francesco , Romano, Pietro , Jiang, Jingjing , Cui, Jiuzhuo , Liu, Jie , Qi, Wei

Journal of Cleaner Production , vol. 490
Citations: 21
Show abstract

© 2025 Elsevier LtdThe accelerated deployment of solar photovoltaic (PV) systems will inevitably result in an increasing volume of end-of-life PV panels, which will pose significant environmental challenges and could potentially hinder the growth of renewable energy systems. This study provided a comprehensive examination of the pyrolysis behavior, kinetics, thermodynamics, and evolved products of typical back sheet PVDF/PET/fluorine film (KPF). In addition, an analysis of the enthalpy-entropy compensation (EEC) was performed. Reactive force field molecular dynamics (ReaxFF-MD) simulations were employed to identify atomic-level intermediates and investigate the reaction pathway. The decomposition of KPF sample occurred in three stages, characterized by temperature ranges of below 473.15 K, 473.15–923.15 K, and 923.15–1173.15 K. The corresponding mass losses were found to be 0.61–0.90, 78.30–82.06, and 1.38–2.56 wt%, respectively. The predominant products identified included benzoic acid and its derivatives, which corroborated the strong presence of the C=O group in the FTIR analysis. ReaxFF-MD simulations revealed the formation of C7H4O2, C7H4O and C7H5O2 species, and the decomposition process was found to involve random scission, decarboxylation and decarbonylation reactions. Activation energies from the FWO, KAS, and Friedman methods exhibited a declining trend, decreasing from 72.12 to 40.92 kJ mol−1. The master-plot analysis indicated that the P2 mechanism provided a more accurate description of KPF pyrolysis. The positive ΔH and ΔG values confirmed that KPF decomposition was an endothermic and non-spontaneous process. The ΔH-ΔS relationship indicated the presence of an EEC, with a compensation temperature of 687.49 K and an experimental temperature of 766.60 K.

Conference Paper 2025

Impact of Ethanol Droplet Size and Methane Addition on the Combustion Performance of Ethanol/LOx Bipropellant Systems: Insights from Reactive Molecular Dynamics Simulations

Gontijo, Mauricio Sa , Goncalves, Rene F.B. , Rossetto, Vinicius , Rocco, Jose A.F.F. , Rocco, Bruno T. , Rocco, Leopoldo

Proceedings of the International Astronautical Congress Iac , vol. 2-F219594 , pp. 988-990
Show abstract

Copyright ©2025 by the International Astronautical Federation (IAF). All rights reserved.Tins study investigates the influence of ethanol droplet size on the combustion performance of an ethanol/hqmd oxygen (LOx) brpropellant system using reactive molecular dynamics (RMD) simulations, with an additional focus on the effect of methane addition to larger ethanol droplets to enhance combustion efficiency. Liquid propulsion systems utilizing ethanol as a fuel offer advantages such as lower toxicity and environmental impact compared to traditional hypergohc propellants. However, the combustion efficiency of ethanol/LOx systems is highly dependent on the atomization and vaporization processes of the fuel, which are directly influenced by droplet size. This work employs RMD to model the molecular level interactions and reaction dynamics between ethanol droplets of varymg sizes (ranging from 20 to 100 Angstroms) and LOx under high-pressure and high-temperature conditions representative of rocket combustion chambers. Additionally, the study explores the incorporation of methane into larger ethanol droplets (50 and 100 Angstroms) to improve then combustion performance. The simulations reveal that smaller ethanol droplets (e.g., 10 A) exhibit faster vaporization rates and more efficient mixing with LOx, leadmg to enhanced combustion performance. Specifically, the ignition delay time for 10 A droplets was found to be 300/c shorter compared to 100 A droplets, while the heat release rate increased by approximately 250/c. For larger droplets, the addition of methane (10% by mass) resulted in a 15% reduction in ignition delay time and a 20% increase in heat release rate, demonstrating its potential to compensate for the lower efficiency of larger droplets. The presence of methane also altered the reaction pathways, promoting the formation of key intermediate species that enhanced overall combustion kinetics. These findings underscore the importance of precise fuel atomization in liquid propulsion systems and suggest that the strategic addition of methane can significantly improve the performance of ethanol/LOx-based engines, particularly when droplet srze minimization is challenging. Future work will focus on validating these results through experimental studies.

Conference Paper 2025

TNT Synthesis Byproducts in Solid Propellants

Gonçalves, Rene F.B. , Rossetto, Vinicius , Rocco, José A.F.F. , Rocco, Bruno T. , Rocco, Leopoldo

Proceedings of the International Astronautical Congress Iac , vol. 1-F219594 , pp. 291-297
Show abstract

Copyright ©2025 by the International Astronautical Federation (IAF). All rights reserved.The synthesis of trinitrotoluene (TNT) produces various byproducts and co-products, often regarded as waste that poses significant environmental and disposal challenges. This study explores the feasibility of repurposing these TNT synthesis co-products, such as nitrated aromatic compounds and nitro derivatives, as partial substitutes for conventional oxidants like ammonium Perchlorate or ammonium nitrate in solid propellant formulations. The research seeks to diminish dependence on traditional oxidants, enhance waste management practices, and foster sustainability in energy material production. The co-products were analyzed for their chemical composition, thermal stability, and energy characteristics, subsequently integrated into propellant mixtures at different concentrations to replace a portion of the primary oxidant. The performance of these modified propellants was assessed through theoretical modeling and experimental evaluations, including RMD simulations conducted across 1000 K to 3000 K. Results reveal that TNT co-products effectively contribute to the oxidizer balance, with the triple-base plus TNT system demonstrating enhanced combustion properties compared to the triple-base alone. Notably, the addition of TNT byproducts increased total energy release by up to ~100000 kcal/mol at 3000 K and sustained higher species counts of intermediates like N02 and CHO, reflecting improved reactivity. The apparent activation energy rose from 198.2 kcal/mol for the triple-base system to 246.4 kcal/mol with TNT byproducts, indicating a higher energy barrier but also greater exothermic potential. Experimentally, incorporating 10% of a nitroaromatic co-product boosted the burning rate by 8% (from 12 mm/s to 13 mm/s) and improved specific impulse to 250 s from 245 s in the standard formulation. Additionally, this approach reduced environmental disposal risks by 15% and offered a potential cost reduction in propellant production of up to 12%. Nevertheless, the study emphasizes the importance of optimizing co-product levels to preserve mechanical and ballistic properties while mitigating risks of sensitivity and instability, as suggested by the increased nitrogen oxide production observed in the simulations.