{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,8,28]],"date-time":"2026-08-28T06:18:15Z","timestamp":1787897895887,"version":"build-2784847793"},"reference-count":21,"publisher":"Frontiers Media SA","license":[{"start":{"date-parts":[[2026,8,28]],"date-time":"2026-08-28T00:00:00Z","timestamp":1787875200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100000084","name":"Directorate for Engineering","doi-asserted-by":"publisher","award":["2130651"],"award-info":[{"award-number":["2130651"]}],"id":[{"id":"10.13039\/100000084","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["frontiersin.org"],"crossmark-restriction":true},"short-container-title":["Front. Bioeng. Biotechnol."],"abstract":"<jats:sec>\n                    <jats:title>Background<\/jats:title>\n                    <jats:p>High-density (HD)-inspired peripheral transcutaneous neurostimulation systems employing clusters of small-area electrode arrays may improve spatial controllability during electrical stimulation of superficial limb tissues, as used in TENS\/NMES\/FES-like wearable applications. However, the electrothermal consequences of small electrodes can affect usability due to complexities in safety considerations. Therefore it is of critical interest to understand the effects of these systems, and to determine how subdivision topology shapes localized current concentration and thermal accumulation under fixed electrode geometry conditions.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Objective<\/jats:title>\n                    <jats:p>This study investigates how high-density clusters of small-area electrodes influences current-density distribution and electrothermal behavior. This project is part of a larger study to design and optimize these high-density neurostimulation arrays to improve spatial scaling.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Methods<\/jats:title>\n                    <jats:p>A coupled electrical\u2013thermal finite-element model was developed in COMSOL Multiphysics 6.3 to compare conventional single-electrode configurations with a systematic series of HD-inspired small-electrode grid configurations ranging from 1 to 16 electrodes (grid1 through grid16). The individual electrode radius and inter-electrode spacing were held constant while the number of active electrodes was varied. Current-density distribution, temperature change, and thermal transient effects were evaluated using current-controlled pulsed stimulation with varying current amplitudes, pulse widths, and frequencies.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>\n                      Electrothermal effects did not increase linearly with the number of active electrodes. Three distinct electrothermal regimes were identified based on stimulation-side behavior: an extreme concentration regime (\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m1\">\n                          <mml:mrow>\n                            <mml:mi>N<\/mml:mi>\n                            <mml:mo>=<\/mml:mo>\n                            <mml:mn>1<\/mml:mn>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      \u20134), a transition regime (\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m2\">\n                          <mml:mrow>\n                            <mml:mi>N<\/mml:mi>\n                            <mml:mo>=<\/mml:mo>\n                            <mml:mn>5<\/mml:mn>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      \u20138), and a regime in which stimulation-side load falls below a fixed return-electrode background contribution (\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m3\">\n                          <mml:mrow>\n                            <mml:mi>N<\/mml:mi>\n                            <mml:mo>\u2265<\/mml:mo>\n                            <mml:mn>9<\/mml:mn>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      , crossover\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m4\">\n                          <mml:mrow>\n                            <mml:msup>\n                              <mml:mrow>\n                                <mml:mi>N<\/mml:mi>\n                              <\/mml:mrow>\n                              <mml:mrow>\n                                <mml:mo>*<\/mml:mo>\n                              <\/mml:mrow>\n                            <\/mml:msup>\n                            <mml:mo>=<\/mml:mo>\n                            <mml:mn>8<\/mml:mn>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      \u20139), where\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m5\">\n                          <mml:mrow>\n                            <mml:mi>N<\/mml:mi>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      is the number of active electrodes in a grid. The\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m6\">\n                          <mml:mrow>\n                            <mml:mi>N<\/mml:mi>\n                            <mml:mo>=<\/mml:mo>\n                            <mml:mn>9<\/mml:mn>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      regime showed the closest stimulation-side electrothermal behavior to the single-electrode reference among all configurations tested (\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m7\">\n                          <mml:mrow>\n                            <mml:mi mathvariant=\"normal\">\u0394<\/mml:mi>\n                            <mml:mi>T<\/mml:mi>\n                            <mml:mo>=<\/mml:mo>\n                            <mml:mn>1.340<\/mml:mn>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      \u00b0C vs. 1.370\u00a0\u00b0C, or a difference of 0.030\u00a0\u00b0C). The stimulation-side\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m8\">\n                          <mml:mrow>\n                            <mml:mi mathvariant=\"normal\">\u0394<\/mml:mi>\n                            <mml:mi>T<\/mml:mi>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      declined with electrode count across this full range (\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m9\">\n                          <mml:mrow>\n                            <mml:mi mathvariant=\"normal\">\u0394<\/mml:mi>\n                            <mml:mi>T<\/mml:mi>\n                            <mml:mo>\u2248<\/mml:mo>\n                            <mml:mn>9.69<\/mml:mn>\n                            <mml:mo>\u00d7<\/mml:mo>\n                            <mml:msup>\n                              <mml:mrow>\n                                <mml:mi>N<\/mml:mi>\n                              <\/mml:mrow>\n                              <mml:mrow>\n                                <mml:mo>\u2212<\/mml:mo>\n                                <mml:mn>0.91<\/mml:mn>\n                              <\/mml:mrow>\n                            <\/mml:msup>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      ,\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m10\">\n                          <mml:mrow>\n                            <mml:msup>\n                              <mml:mrow>\n                                <mml:mi>R<\/mml:mi>\n                              <\/mml:mrow>\n                              <mml:mrow>\n                                <mml:mn>2<\/mml:mn>\n                              <\/mml:mrow>\n                            <\/mml:msup>\n                            <mml:mo>=<\/mml:mo>\n                            <mml:mn>0.999<\/mml:mn>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      ). A dedicated area-matched control series, in which electrode radius was increased to hold total active area constant showed a decline in stimulation-side current density and a temperature rise with electrode count. A complementary, exploratory comparison holding electrode count fixed while varying the total active area also showed decline in both quantities. Together, these results indicate that electrode count and total active area each contribute to the observed electrothermal decline. However, causal effects were not isolated from one another. As expected from established current-concentration effects of small electrodes, the sparse grid1 configuration produced the highest superficial current density and temperature rise among the investigated configurations (\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m11\">\n                          <mml:mrow>\n                            <mml:mi mathvariant=\"normal\">\u0394<\/mml:mi>\n                            <mml:mi>T<\/mml:mi>\n                            <mml:mo>=<\/mml:mo>\n                            <mml:mn>10.18<\/mml:mn>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      \u00b0C, or 7.43\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m12\">\n                          <mml:mrow>\n                            <mml:mo>\u00d7<\/mml:mo>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      reference-electrode baseline). A discretized parameter-space analysis (128 sampled current amplitude\/pulse-width\/frequency combinations) showed that the\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m13\">\n                          <mml:mrow>\n                            <mml:mi mathvariant=\"normal\">\u0394<\/mml:mi>\n                            <mml:mi>T<\/mml:mi>\n                            <mml:mo>\u2265<\/mml:mo>\n                            <mml:mn>8<\/mml:mn>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      \u00b0C region occupied by each configuration declined continuously with electrode count. Specifically, from 72 combinations for grid1 (4.50\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m14\">\n                          <mml:mrow>\n                            <mml:mo>\u00d7<\/mml:mo>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      single) to 4 combinations for grid16 (0.25\n                      <jats:inline-formula>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m15\">\n                          <mml:mrow>\n                            <mml:mo>\u00d7<\/mml:mo>\n                          <\/mml:mrow>\n                        <\/mml:math>\n                      <\/jats:inline-formula>\n                      single). Therefore, electrode grid configurations with higher electrode counts distribute the current across a progressively larger area, and thus reduced per-electrode current density.\n                    <\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Conclusion<\/jats:title>\n                    <jats:p>Under the idealized modeling conditions investigated here, electrothermal behavior was jointly determined by electrode count and total active electrode area. These two variables are intrinsically coupled in the fixed-radius grid1\u2013grid16 design where both contribute to a reduction in stimulation-side electrothermal load. These findings represent comparative, geometry-specific modeling trends rather than validated thermal safety predictions. The present study provides a comparative computational framework for evaluating HD-inspired electrode subdivision strategies in compact peripheral transcutaneous neurostimulation design, where electrode count and total active area jointly determine the resulting current-distribution topology and associated localized thermal accumulation.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.3389\/fbioe.2026.1904885","type":"journal-article","created":{"date-parts":[[2026,8,28]],"date-time":"2026-08-28T05:33:43Z","timestamp":1787895223000},"update-policy":"https:\/\/doi.org\/10.3389\/crossmark-policy","source":"Crossref","is-referenced-by-count":0,"title":["Comparative electrothermal analysis of single and small-electrode grid configurations for HD-inspired peripheral transcutaneous neurostimulation"],"prefix":"10.3389","volume":"14","author":[{"given":"Yuanshan","family":"Zhong","sequence":"first","affiliation":[{"name":"Department of Biomedical Engineering, School of Engineering and Applied Sciences, University at Buffalo","place":["Buffalo, NY, United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Filip","family":"Stefanovic","sequence":"additional","affiliation":[{"name":"Department of Biomedical Engineering, School of Engineering and Applied Sciences, University at Buffalo","place":["Buffalo, NY, United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1965","published-online":{"date-parts":[[2026,8,28]]},"reference":[{"key":"B1","doi-asserted-by":"publisher","first-page":"163","DOI":"10.1088\/1741-2560\/5\/2\/007","article-title":"Transcranial current stimulation focality using disc and ring electrode configurations: FEM analysis","volume":"5","author":"Datta","year":"2008","journal-title":"J. Neural Eng."},{"key":"B2","doi-asserted-by":"publisher","first-page":"201","DOI":"10.1016\/j.brs.2009.03.005","article-title":"Gyri-precise head model of transcranial direct current stimulation: improved spatial focality using a ring electrode versus conventional rectangular pad","volume":"2","author":"Datta","year":"2009","journal-title":"Brain Stimul."},{"key":"B3","doi-asserted-by":"publisher","first-page":"046011","DOI":"10.1088\/1741-2560\/8\/4\/046011","article-title":"Optimized multi-electrode stimulation increases focality and intensity at target","volume":"8","author":"Dmochowski","year":"2011","journal-title":"J. Neural Eng."},{"key":"B5","doi-asserted-by":"publisher","first-page":"2271","DOI":"10.1088\/0031-9155\/41\/11\/003","article-title":"The dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum of tissues","volume":"41","author":"Gabriel","year":"1996","journal-title":"Phys. Med. Biol."},{"key":"B6","doi-asserted-by":"publisher","first-page":"8825","DOI":"10.1088\/1361-6560\/61\/24\/8825","article-title":"Effect of microscopic modeling of skin in electrical and thermal analysis of transcranial direct current stimulation","volume":"61","author":"Gomez-Tames","year":"2016","journal-title":"Phys. Med. Biol."},{"key":"B7","doi-asserted-by":"publisher","first-page":"036020","DOI":"10.1088\/1741-2560\/13\/3\/036020","article-title":"Optimization of focality and direction in dense electrode array transcranial direct current stimulation","volume":"13","author":"Guler","year":"2016","journal-title":"J. Neural Eng."},{"key":"B8","doi-asserted-by":"publisher","first-page":"2207251","DOI":"10.1002\/advs.202207251","article-title":"Focused epicranial brain stimulation by spatial sculpting of pulsed electric fields using high density electrode arrays","volume":"10","author":"Jain","year":"2023","journal-title":"Adv. Sci."},{"key":"B9","doi-asserted-by":"publisher","first-page":"7","DOI":"10.1177\/204946370700100103","article-title":"Transcutaneous electrical nerve stimulation: mechanisms, clinical application and evidence","volume":"1","author":"Johnson","year":"2007","journal-title":"Rev. Pain"},{"key":"B10","doi-asserted-by":"publisher","first-page":"235008","DOI":"10.1088\/1361-6560\/abb7c1","article-title":"Role of skin tissue layers and ultra-structure in transcutaneous electrical stimulation including tDCS","volume":"65","author":"Khadka","year":"2020","journal-title":"Phys. Med. Biol."},{"key":"B18","doi-asserted-by":"publisher","first-page":"349","DOI":"10.5021\/ad.2019.31.3.349","article-title":"Regional variation of human skin surface temperature","volume":"31","author":"Lee","year":"2019","journal-title":"Ann. Dermatol"},{"key":"B12","doi-asserted-by":"publisher","first-page":"4","DOI":"10.20517\/ss.2023.36","article-title":"A review: flexible devices for nerve stimulation","volume":"4","author":"Liu","year":"2024","journal-title":"Soft Sci."},{"key":"B13","doi-asserted-by":"publisher","first-page":"223","DOI":"10.1007\/s00421-010-1502-y","article-title":"Physiological and methodological considerations for the use of neuromuscular electrical stimulation","volume":"110","author":"Maffiuletti","year":"2010","journal-title":"Eur. J. Appl. Physiology"},{"key":"B14","doi-asserted-by":"publisher","first-page":"66","DOI":"10.1186\/1743-0003-9-66","article-title":"A multi-pad electrode based functional electrical stimulation system for restoration of grasp","volume":"9","author":"Male\u0161evi\u0107","year":"2012","journal-title":"J. NeuroEngineering Rehabilitation"},{"key":"B15","doi-asserted-by":"publisher","DOI":"10.1186\/s12938-026-01610-w","article-title":"Computational analysis of a high-density electrode to scale the electric field geometry and intensity for neurostimulation","author":"Mancabelli","year":"2026","journal-title":"Biomed. Eng. Online"},{"key":"B16","doi-asserted-by":"publisher","first-page":"681","DOI":"10.3390\/life11070681","article-title":"Perception of thermal comfort during skin cooling and heating","volume":"11","author":"Mekjavic","year":"2021","journal-title":"Life"},{"key":"B17","doi-asserted-by":"publisher","first-page":"501","DOI":"10.1149\/1.2424003","article-title":"Resistance for flow of current to a disk","volume":"113","author":"Newman","year":"1966","journal-title":"J. Electrochem. Soc."},{"key":"B19","doi-asserted-by":"publisher","first-page":"93","DOI":"10.1152\/jappl.1948.1.2.93","article-title":"Analysis of tissue and arterial blood temperatures in the resting human forearm","volume":"1","author":"Pennes","year":"1948","journal-title":"J. Appl. Physiology"},{"key":"B4","doi-asserted-by":"publisher","first-page":"1651","DOI":"10.1109\/TNSRE.2020.2994900","article-title":"Design of transcutaneous stimulation electrodes for wearable neuroprostheses","volume":"28","author":"RaviChandran","year":"2020","journal-title":"IEEE Trans. Neural Syst. Rehabil. Eng."},{"key":"B20","doi-asserted-by":"publisher","first-page":"2982","DOI":"10.3390\/s23062982","article-title":"Conformable electrode arrays for wearable neuroprostheses","volume":"23","author":"RaviChandran","year":"2023","journal-title":"Sensors"},{"key":"B21","doi-asserted-by":"publisher","first-page":"25","DOI":"10.1016\/j.neuroimage.2015.06.067","article-title":"On the importance of electrode parameters for shaping electric field patterns generated by tDCS","volume":"120","author":"Saturnino","year":"2015","journal-title":"NeuroImage"},{"key":"B22","doi-asserted-by":"publisher","first-page":"e50309","DOI":"10.3791\/50309","article-title":"Technique and considerations in the use of 4 \u00d7 1 ring high-definition transcranial direct current stimulation (HD-tDCS)","author":"Villamar","year":"2013","journal-title":"J. Vis. Exp."}],"container-title":["Frontiers in Bioengineering and Biotechnology"],"original-title":[],"link":[{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fbioe.2026.1904885\/full","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,8,28]],"date-time":"2026-08-28T05:33:44Z","timestamp":1787895224000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fbioe.2026.1904885\/full"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,8,28]]},"references-count":21,"alternative-id":["10.3389\/fbioe.2026.1904885"],"URL":"https:\/\/doi.org\/10.3389\/fbioe.2026.1904885","relation":{},"ISSN":["2296-4185"],"issn-type":[{"value":"2296-4185","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,8,28]]},"article-number":"1904885"}}