Multiphysics finite element investigation of galvanic transmission in dynamic human body communications

Abstract

Introduction:

Human body communication (HBC) utilizes the human body as a medium of communicating data. Considerable research has been done to characterize HBC channels to optimize communication techniques. However, dynamic HBC channels have been less studied.

Methods:

An approach for developing dynamic models of the human body channel for galvanic communication is presented using multiphysics finite element analysis (FEA). An analytical framework is formulated that utilizes stochastic ABCD network parameters to explore and model dynamic HBC channel segments. Channel segments were subjected to mechanical forces using the finite element method (FEM) to reveal their impact on the current density and electric field.

Results:

Linear regression modeling shows a strong relationship between applied force, current frequency, and channel response, with R² metrics exceeding 0.99. The dynamic nature of the channel reflects the need for stochastic modeling. This study examined candidate probability density functions (PDFs) to describe channel fading for the ABCD network parameters. Lognormal and Weibull distributions fit the magnitudes best while the generalized Pareto, generalized extreme value, and logistic distributions fit the phases best. Empirical modeling validated the accuracy of the lognormal distribution fits found using the FEM.

Discussion:

The dynamic channel was characterized utilizing multiphysics FEM modeling, empirical modeling, and ABCD network parameters. This information is invaluable for EM dosimetry analysis and risk assessment in body area network (BAN) device design, as well as device optimization, because stochastic HBC parameters emulate the dynamic nature of the human body channel.

1 Introduction

Human body communication (HBC) utilizes the human body as a medium for communicating data (IEEE SA, 2008). HBC can be considered as an alternative to short-range RF communication when used for medical telemetry in body area networks (BANs) (Pereira et al., 2015). Research has shown that HBC achieves greater spectral efficiency, higher data rates, better security, and better power efficiency than radio-frequency (RF)-based BANs (Bae et al., 2012b; Cho et al., 2014; Seyedi and Lai, 2014; Ali et al., 2019; Sujaya and BhanuPrashanth, 2021). However, channel models are essential for analyzing and optimizing BAN performance. These models can be leveraged to develop application-specific physical layer (PHY) technologies (including, but not limited to, automatic gain control (AGC), channel equalization, and auto-encoders) that can be used to improve the performance of BAN transceivers to achieve better error performance, better receiver sensitivity, a higher data rate, better channel resilience, and power efficiency (Ali et al., 2019; Zhang et al., 2016b; Pereira et al., 2015). Furthermore, human body state (e.g., electrode location/state, action performed by body parts, and blood analyte concentration such as glucose) correlates with the channel response (Wang et al., 2016; Zhang et al., 2016a; Maity et al., 2018; Takamatsu et al., 2021; Roopnarine and Rocke, 2024; Tang et al., 2023; Chen et al., 2016; Ahmed et al., 2019). This type of data is crucial in medical use cases that require the monitoring of patients’ medical condition or athletes’ physical state to optimize their performance (Bouazizi et al., 2017; Smith, 2011; Dhamdhere et al., 2010). Thus, channel models are essential for BAN optimization, medical state determination, and RF safety.

Extensive research has developed channel models for the two main categories of HBC communication techniques: electric field HBC (eHBC), which involves modulating electric fields to communicate through galvanic and capacitive coupling, and magnetic HBC (mHBC), which involves modulating magneto-quasistatic fields through magnetic induction (e.g., Callejon et al., 2012a; Pereira et al., 2015; Ogasawara et al., 2014; Seyedi and Lai, 2014; Rocke and Persad, 2015; Wen et al., 2022; Ahmed et al., 2019; Chen et al., 2016; Vizziello et al., 2024; Noormohammadi et al., 2021). Each communication technique has its advantages. For example, mHBC offers channel resilience since the magnetic permeability of most materials is the same; capacitive coupling has a greater frequency range of operation that offers higher data rates; both mHBC and capacitive coupling allow for electrodes to be either on or close to the body, offering more versatility with regards to use cases; galvanic coupling offers more channel resilience against the body’s environment as it requires electrodes to be in contact with the body (Rocke and Persad, 2015; Pereira et al., 2015; Seyedi and Lai, 2014; Wen et al., 2022). This paper focuses on modeling galvanic coupling because this form of communication is standardized for HBC. However, this research can be easily extended to other HBC techniques.

However, relatively limited emphasis has been placed on the effect of the dynamic state of the human body on the channel response. Analytical solutions typically explore the propagation of signals through the human body using Maxwell’s equations (Bae et al., 2012a). However, these solutions are closed-form and limited to standard geometries; they cannot model dynamic behavior that may cause irregular tissue deformations. In the literature surveyed, mostly empirical measurements were found to model dynamic channel behavior for HBC (Maity et al., 2018; Zhang et al., 2016b; Zhang et al., 2016a; Seyedi and Lai, 2014; Tang et al., 2023; Chen et al., 2016). Nonetheless, the fading characteristics found provide no insight into the electric field and current distribution of the body under such conditions. Furthermore, no study was found on the impact of externally applied forces to the skin on the body’s channel response in the HBC band (18.375–23.625 MHz). Maity et al. (2018) demonstrated channel loss and variability increases with decreasing pressure applied to electrodes through an empirical measurement study. However, that study was limited to the 1Hz-1 MHz frequency range and did not quantify the electrode pressures applied. HBC is standardized in the HBC band for galvanic communication (IEEE, 2012). Consequently, for standardized communication, HBC models must include the HBC band. Numerical solutions offer insight into complex structures by evaluating the foundation equations that defined the physics of a model proposed for experiments (Rocke and Persad, 2015).

Finite element method (FEM) analysis a type of numerical solution that allows the dynamic characterization of biological phenomena (Guo et al., 2011) and electro-quasistatic field visualization for defined HBC use cases (Callejon et al., 2014; Noormohammadi et al., 2021; Datta et al., 2021; Ahmed et al., 2019). This gives insight into the electrical transmission mechanism through human tissue, which is invaluable when considering radio frequency (RF) safety and transceiver optimization techniques. For example, Callejon et al. (2014) simulated current density and electric field strength across the tissue layers in the arm for different channel lengths and inter-electrode distances across 1 kHz to 1 MHz. However, for this study, the electric field visualizations of the cross-section of the arm shown did not include the HBC band. Ahmed et al. (2019) extended this work by considering different tissue thicknesses and bending angles, looking at the electric potential at the Rx for these dynamic conditions. Current density visualizations were provided at the 10 kHz frequency for some geometries where the thickness of the tissue composition was varied. Additionally, this study did not consider the HBC band. Datta et al. (2021) showed the electric field across and around an overly simplified geometry of the human body for capacitive coupling communication technique for frequencies under 1 MHz. They also investigated the variation of received voltages for different body postures and electrode configurations. However, there was no indication of the incorporation of the mechanical tissue properties which, when included, can affect the resultant bio-impedance of the body postures observed (Albulbul and Chan, 2012).

This paper is organized as follows. Section 2 introduces the ABCD network parameter framework, outlines FEM model characteristics, and describes the empirical modeling approach. Section 3 presents the simulated and empirical results and analyzes them. Section 4 presents the conclusions made based on the results observed. This study extends Roopnarine et al. (2024).

2 Materials and methods

2.1 System model

Galvanic coupling involves the application of a deferentially modulated signal through transmitter electrodes (Tx) directly in contact with the body (Figure 1). This induces galvanic currents, which are collected by receiver (Rx) electrodes also in direct contact with the body (Pereira et al., 2015).

FIGURE 1

Fading—the variation of channel gain—is related to bio-impedance since it affects the channel response. Several studies show that the electric field of the propagated signal is distributed through the different tissue layers (Callejon et al., 2012a; Callejon et al., 2012b, Callejon et al., 2014; Wang et al., 2016; Ahmed et al., 2019). Consequently, circuit-based analytical HBC channel models have been developed that incorporate the electrical parameters of the human body (Callejon et al., 2012a; Wang et al., 2016; Seyedi and Lai, 2014; Callejon et al., 2012b).

Bio-impedance, and by extension the channel response, is affected by the body’s state. Gabriel et al. (1996) demonstrated that the body’s circuit parameters depend on its electrical properties—tissue conductivity and permittivity. However, these properties are dependent on non-deterministic, time-varying factors such as channel geometry—for example, skin thickness, skin deformation (Albulbul and Chan, 2012; Maity et al., 2018), body posture (Datta et al., 2021) and presence of needle punctures (White et al., 2013)—and body chemical composition—for example, blood glucose concentration (Andersen et al., 2019). Bio-impedance is also different for different sections of the body and for subsections such as the dermis and epidermis skin layers (Tsai et al., 2019; Gabriel, 1996). The state of these sections also influences bio-impedance, such as if skin is wet or dry (Gabriel et al., 1996). The body’s biological condition (e.g., biological state of liver ischemia) influences bio-impedance (Tronstad and Strand-Amundsen, 2019).

The ABCD network parameters are essential for analyzing and optimizing human body communication (HBC) systems. The A and D parameters represent input and output voltage ratios that are vital in assessing signal gain and loss. At the same time, B and C describe the voltage–current relationship, thus revealing impedance and loading effects. Understanding these parameters enables the accurate modeling of signal propagation, efficient energy transfer, and enhanced reliability in dynamic environments. Integrating ABCD parameters into the channel model provides a robust framework for improving HBC performance and guiding design strategies.

ABCD network parameters are used in this study because of their effectiveness in modeling segmented systems as a transmission line. Callejon et al. (2012a) demonstrated that the human body can be modeled as a static transmission line consisting of a series of cascaded two-port networks to analyze signal propagation. Our research extends this by introducing ABCD parameters, leveraging the premise that two-port networks can model the channel response of individual tissue segments. This segmented approach allows for the cascading of segments to represent the entire channel as a lossy transmission line, as in transmission line theory (Notaros, 2010; Ulaby et al., 2010). Thus, any cascade can then be used to represent the entire length of the channel. The stochastic model is based upon the traditional use of deterministic ABCD parameters to model transmission lines (Callejon et al., 2012a; Ulaby et al., 2010; Notaros, 2010). Our approach involves extending this ABCD framework by introducing stochasticity to the ABCD parameters of each of these tissue segments, acknowledging the inherent randomness in tissue properties. This enables us to explore the resulting PDFs that describe the transmission parameters of the entire channel under dynamic conditions, thereby providing insights into dynamic channel fading.

Accordingly, the human body channel is seen as non-deterministic and time-varying. This work proposes that the human body be modeled as a lossy transmission line divided into differential segments, such that the body segment, , is depicted in Figure 2. Hence, each differential segment of the human body channel, , is modeled as an ABCD 2-port network shown in Equation 1.

FIGURE 2

This on-body segmented model allows for the representation of fading caused by random channel response variations that stem from bio-impedance fluctuations. is a 2 × 2 matrix representing the ABCD parameters for the segment. Each element in is a complex-valued random process with probability density functions (PDFs) dependent on , and (designated as , and respectively). Thus, is a random process continuous-valued, discrete-parameter random matrix, consisting of random processes , , and , with index set , where is the number of segments that aggregate to form the channel. The channel length under observation is , where is the length of the segment.

ABCD network modeling allows for the entire human body channel under observation to be represented as an aggregate two-port network equivalent to the cascade of the differential two-port networks. This forms the resultant HBC network representation under matrix multiplication, as shown in Equation 2:Thus, is the end-to-end channel ABCD matrix as a result of cascading all the ‘s in the channel.

This approach can be extended to the other HBC communication modes (i.e., capacitive coupling and galvanic coupling) through modification of the equivalent circuit to reflect the signal pathways in the channel.

2.2 Multiphysics finite element dynamic HBC model

The FEM model was derived using COMSOL multiphysics software through geometry design, material assignment, multiphysics selection, boundary condition assignment, mesh selection, domain configuration, and result processing configuration.

2.2.1 Geometry and model parameters

The geometry modeled is a multi-layered cuboid section of human tissue consisting of four tissue components—skin, fat, muscle, and a blood vessel—with Tx and Rx electrode pairs placed on-body (Figure 3). The model dimensions are captured in Table 1. These tissue component thicknesses are within the range of the true human anatomical dimensions, similar to Callejon et al. (2012a), Wegmueller et al. (2007), Persad and Rocke (2020), and DeSaix et al. (2013). The dielectric properties of the tissue components (i.e., permittivity, , and conductivity, ) were chosen based on the Cole–Cole model by Gabriel et al. (1996) to capture their frequency response when subject to electro-quasistatic fields in the HBC band (Gabriel, 1996).

FIGURE 3

The mechanical properties of the tissue components (i.e., Young’s modulus, , Poisson’s ratio, , and density, ) were chosen based on the mechanical properties of human tissue used and derived in the literature surveyed. This is captured in Table 2.

2.2.2 Model equations

The FEM model proposed was derived using COMSOL Multi-physics 5.5 through the physics of electrical currents, electrical circuits, and solid mechanics. The electrical currents and circuit physics are deployed through the AC/DC module which solves the current conservation equation based on Ohm’s law—and, by extension, Kirchhoff’s conservation laws for the voltages, currents, and charges associated with the circuit elements—using the scalar electrical potential as the dependent variable, as shown in Equations 3–5.where is current density , is an external current source , is the electrical conductivity , is an externally generated current density , is the electrical field intensity , is scalar electrical potential , is electrical displacement , and is angular frequency .

In accordance with the convention for galvanic communication, a differential signal was applied across the Tx electrodes by setting the Dirichlet boundary condition to a voltage of amplitude of frequencies in the HBC band. is specified as peak voltage (COMSOL Multiphysics, 2019). The channel was surrounded by an external layer of air (a sphere surrounding the channel of radius ) configured as electric insulation by setting the Neumann condition, , where is the surface normal. Thus, the air domain was considered an external unbounded domain. A resistor component is assigned between the Rx electrodes to read the electrical parameters necessary to derive the ABCD parameters for developing the FEM model. Separate reference nodes for the Rx and Tx pairs were configured, with Dirichlet boundary conditions of through the relevant physics, to ensure that the Tx and Rx ports were decoupled.

The solid mechanics physics is deployed through the structural mechanics module, which solves equations of motion together with a constitutive model for a solid material to compute displacements, stresses, and strains as shown in Equation 6.where is the deformation gradient, is second Piola–Kirchhoff stress , is the force per unit volume vector , is the unit tensor, and is the displacement field .

Each tissue component was configured to be linear elastic material to emulate the effects of applying practical, small-signal external stresses expected under normal operating conditions. The bottom surface of the FEM channel model, under the muscle tissue component, was configured to be a fixed constraint (i.e., its displacements are zero in all directions) to model the multi-layered tissue’s contact with bone. A boundary load, or force in , was applied on the upper surface of the FEM model, on the skin tissue component, to emulate pushes and pulls on the model’s surface from human activities (Lee and Jung, 2017).

2.2.3 Model computation

The simulation software discretizes the channel into finite elements using physics-controlled meshing sequences, which create meshes that consist of different element types and size features, to solve the applicable model equations defined in Section 2.2.2. When these model equations are solved, the relevant field parameters can be extracted to formulate the FEM model.

Recall that this FEM model focusing on the HBC band as galvanic communication is standardized at those frequencies (IEEE, 2012). The electrical length of the human body and signal frequency in the HBC band show that the communication scenario observed satisfies that of electro-quasistatics. However, as the signal frequency becomes larger, the quasistatic assumption may not hold, as propagation and inductive effects may become more significant. Thus, to avoid invalidating the quasistatic assumption, the frequencies are kept under 100 MHz. With galvanic communication, the electrical field is mostly localized to the human body (Callejon et al., 2014). The dimensions of the external layer of air were selected based on the extent of the electrical field observed by Callejon et al. (2014) and to accommodate the geometric deformations applied in this study.

2.2.4 Simulations

The channel was deformed with forces applied on the skin’s surface, ranging in magnitude from −30 N to 30 N in 1 N increments. These force magnitudes were selected based on research from typical activities involving the human body, such as pinching, gripping an object (Lee and Jung, 2017; Wells and Greig, 2001), tactile forces (Tang et al., 2015), indentations from micro-needles (Groves et al., 2012), and internal muscle activations (Steenbrink et al., 2009). All of these activities exert forces on the skin in various directions. In reality, the skin is subject to forces in multiple directions. Thus, this study applied mechanical forces to simulate pressing and pulling forces on the tissue.

For each of these deformations, the Tx electrodes were excited with of frequencies in the HBC band (18.375–23.625 MHz in 0.25 MHz increments). Each voltage was injected with two different load resistances ( and ) between the Rx electrodes to compute ABCD parameters. This data was used to determine the relationship between the network parameters, force applied, , and frequency of the voltage used, .

The electric field, current density, and voltage field distributions were simulated on the cross-sectional plane corresponding to the Tx electrodes in the plane of the channel at different channel lengths (, , , and ) over different frequencies ( and ) under different deformation forces (,-, , and ). Insights into the electric field penetration of the channel were acquired by determining the average electric field magnitude and current density magnitude for each tissue component under different deformation forces (,-, , , and ). Negative forces represent pressing forces onto the skin while positive forces represent pulling on the skin.

The channel was divided into segments, each composed of all tissue components, (Figure 4). The ABCD parameters were determined, using the currents and voltages at input and output ports recorded from the simulation for each segment as well as the whole channel between the Tx and Rx electrodes; is the reverse voltage transfer ratio (or inverse channel gain), is the short circuit resistance, is the open circuit conductance, and is the reverse current transfer ratio. These values were subsequently used to perform dispersion and correlation analyses between segments.

FIGURE 4

2.3 Empirical dynamic HBC model

The objective of this section is to empirically examine the distributions that the ABCD parameters tend toward under typical forces applied during every day dynamic activities. This analysis serves to validate the channel response distributions obtained through numerical simulations—one of the areas for further study identified by Roopnarine and Rocke (2021).

Thus, in addition to previously simulations, the candidate PDFs were examined through channel measurements of typical dynamic conditions expected in HBC use cases. The conditions studied in both simulation and empirical research were aligned to ensure consistency in the scenarios being tested. Multiple scans with averaging were used during measurements to control for experimental error, improving the reliability of the results. Transmit (Tx) and receive (Rx) electrodes, spaced 10 mm apart, were placed on a human subject in the general locations shown in Figure 5 (arm, wrist, and wrist-arm). These electrodes were connected to the Anritsu MS46122B vector network analyzer (VNA), which was used to measure the S-parameters (, , , and through reflection and transmission mode configurations) of the channel in the general locations isolated. The Tx and Rx ports were decoupled using FTB-1-1+ Baluns. The VNA connections are also shown in Figure 5. These S-parameters were measured with a human in various states: fist clenching, writing, a bicep curl without weight, a bicep curl with weight, rotating ankle, walking, and marching. In these states, forces were naturally applied to the skin, affecting channel behavior. These S-parameters were then converted to ABCD parameters in alignment with the analytical framework outlined in Section 2.1 (Blattenberger, 2023). The VNA was calibrated to offset the effect the connector probes have on the channel response. Additionally, averaging was used to minimize measurement noise. These measurements were taken across the HBC band ( and ). Additionally, the output power of the VNA was set to 0dBm to comply with the ICNIRP standard for general exposure for humans (International Commission on Non-Ionizing Radiation Protection ICNIRP, 1998).

FIGURE 5

After acquiring the data, both the empirical and simulated results were compared following simulation validation. This comparison helped identify candidate distributions for dynamic fading, capturing the stochastic nature of the dynamic HBC channel. The agreement between the simulated and measured data helped ensure that the identified PDFs could reliably represent the channel dynamics under typical HBC conditions.

3 Results and discussion

The results and analysis are broken into field-based analysis and analysis of the stochastic dynamic HBC channel model.

3.1 EM-based field analysis

Electromagnetic (EM)-based field analysis is an important aspect of electromagnetic EM dosimetry and risk assessment. In wearable and implanted system design, simulating the distribution of electric field intensity in the different tissues is essential for i) understanding the potential effects of EM radiation on human health, ii) ensuring compliance with applicable safety standards, iii) optimizing device design, and iv) facilitating personalized dosimetry analysis.

Figure 6 shows the current passing through a cross-sectional plane of the tissue components located at distances 0, 10, 20, and 30 cm (i.e., at the Rx electrodes) from the Tx electrode. These results show that the tissue current density decreases as the distance from the Tx increases, consistent across all tissue components, similar to observations made by Callejon et al. (2014) at frequencies outside of the HBC band. Additionally, the applied force significantly impacts the distribution of current density in each tissue component. Although changes in the applied pressing force cause relatively minor variations in current density, pulling force deformations notably increase the current density across all tissue components, except near the Tx electrodes. This is likely because deformation alters the geometry of the skin, increasing the cross-sectional area through which current flows, which in turn reduces resistance and impacts the current distribution.

FIGURE 6

Figure 7 shows the percentage of the electric field intensity, , for deformed tissue components of skin, blood, fat, and muscle respectively over the entire channel length. These results show that, excepting muscle tissue, increasing a downward force on all tissues decreases the percentage of across all tissue components. In contrast, increasing the downward force in the muscle increases the percentage . Furthermore, with the exception of blood, changes in the percentage remain relatively flat across frequency. Most of is concentrated in the fat, followed by muscle, skin, then blood. Knowing the E-field distribution across these different tissue layers allows for better optimization of the communication system design process as well as RF safety completeness checks (International Commission on Non-Ionizing Radiation Protection ICNIRP, 1998).

FIGURE 7

Figure 8 shows the percentage of electric current density, , for tissue components skin, blood, fat, and muscle, respectively, over the entire channel length. Similar to , with the exception of muscle tissue, the percentage current density decreases with increasing downward force on all tissues. Increasing the downward force in muscle increases the percentage . Furthermore, changes in the percentage remain relatively flat across frequency across all tissues. Most of is concentrated in muscle, followed by skin, fat, then blood, similar to the static study of Callejon et al. (2014).

FIGURE 8

3.2 Stochastic ABCD model for dynamic HBC channel

Several trends were observed from analyzing the network parameters for the 30 cm channel length as a whole and considering the segments. The network parameters, force applied, , and frequency used, exhibited a linear relationship (of general form shown in Equation 7) derived in Table 3.

The equations derived for magnitude and phase fitting for the network parameters for the 30 cm channel length were fitted using linear regression to map the channel response as a function of the frequency and force applied. These fits demonstrated an impressive goodness of fit metric. is a statistical measure used to determine how well a regression model fits with the observed data (Kutner et al., 2005; Draper and Smith, 1998). Thus, the channel response and force applied demonstrated a linear relationship.

The use of , , , and highlights how dynamic conditions impact the full range of parameters in HBC channels. , the inverse channel gain, and , representing impedance, help capture how movement and tissue deformation affect signal attenuation and impedance variations. Similarly, , which represents the transfer of current, and , the reflection coefficient, are also crucial for modeling signal behavior under dynamic conditions. Changes in and offer insights into how current transfer and signal reflection are influenced by body motion and varying tissue properties, providing a more comprehensive understanding of channel fading, signal loss, and impedance shifts in real-world HBC scenarios.

The correlation matrix for the magnitude and phases of the network parameters of all the segments are shown in Supplementary Figures S1 and S2 respectively (located in the Supplementary Data Sheet S1). The results intuitively show that closer segment’s network parameters are more correlated, with higher relative correlation coefficients, than segments spaced further apart. By definition, is the inverse channel gain. From Supplementary Figure S3 (located in Supplementary Data Sheet S1), the force applied to the skin affects the correlation of the channel response between segments such that increasing the force, either in positive or negative directions, results in a stronger correlation between closer segments. The general trend observed when a force of is applied is the same when varying force. This correlation information plays a critical role in modeling the overall channel, particularly for simulations like Monte Carlo methods. The analysis of segments provides essential insights into how different body segments interact with under varying conditions. This allows researchers to accurately model the relationship between adjacent segments, which can then be extended to model the full transmission path from the transmitter (Tx) to the receiver (Rx). By understanding the correlation between these segments, researchers can combine multiple segments in a correlated manner to derive a more realistic and dynamic channel model. This is essential for generating accurate received signal models, which are vital for performance analysis, design optimization, and ensuring robust communication in human body communication (HBC) scenarios.

The best-fit PDFs were determined for each of the parameters for each segment using the Akaike information criterion (AIC), the Bayesian information criterion (BIC) and NLogL goodness of fit metrics, and the parameters derived for each channel variation as a separate sample. With regards to the magnitude, the results show that the Weibull distribution fits and parameters best while the lognormal fits the and parameters best. With regards to the phase, the results show the generalized Pareto distribution fits best for all parameters. Figure 9 shows a sample PDF plot of the magnitude of the parameter elements for one of the segments.

FIGURE 9

When all the segments were cascaded to form , the best fit PDFs for both magnitude and phase elements were found. Regarding the magnitude, the results show that the lognormal distribution fit best for all parameters. The multiphysics simulation-based approach thus provides further support for the findings of Roopnarine and Rocke (2021), in which Monte Carlo analysis of the ABCD cascade was done for assumed distributions for the segments. Lognormal distributions typically model phenomena where a significant number of individual effects, not strictly independent, all act on a signal (Saunders and Aragón-Zavala, 2007). Regarding the phase, the results show that the generalized extreme value distribution fits the and parameters best, while the logistic distribution fits the and parameters best. Figure 10 illustrates the PDF plots of the magnitude of the parameter elements for the cascaded elements of .

FIGURE 10

The empirically derived fading characteristics from dynamic HBC scenarios investigated for can be found in Figure 11. Here, the best fit PDFs for the ABCD parameters were found by accounting for the dynamic channels investigated. Under these conditions, forces were applied to the skin through these natural positions. The results show that the lognormal distribution fits the ABCD parameters of the overall channel. This corroborates the simulated data as the top fitting distribution match.

FIGURE 11

Consequently, there is an established relationship between applied forces and the channel response in the HBC band. This was proven using the FEM framework that incorporates the ABCD network parameters in Tables 3 and 4. This same FEM model allows the current, current density, and electric field intensity across body tissues under the dynamic conditions (i.e., channel length and applied force variations) to be observed from Figures 6–8. Channel variability (from looking at these electrical parameters in the FEM model produced) for galvanic coupling in the HBC band has been observed when the channel is subject to dynamic conditions. Further, there is a correlation between adjacent surface body segments and tissue layer under these dynamic channel conditions in the results produced from the FEM model in Supplementary Figures S1 and S2. Finally, the empirical model produced strengthens the credibility of the FEM model.

4 Conclusion

A framework for multiphysics finite element analysis of dynamic HBC channels was presented, incorporating the use of stochastic ABCD network parameters. The approach was demonstrated for a galvanic communication scenario, for which tissue deformation was simulated to model the dynamic human body channel. The results showed that applied force affected the distribution of the current, current density, and electric field intensity across the different tissue layers of the channel—important elements in EM dosimetry and risk assessment studies. A linear regression model was formulated that shows the relationship between applied force, frequency of current, and the channel response. This framework can be used to consider other HBC communication modes such as capacitive and magnetic coupling as well as different dynamic conditions, such as different types of forces, environment, body posture, and activities including sitting/standing/walking states.

The dynamic nature of the channel reflected the need for stochastic modeling. Thus, the best PDFs for each of the ABCD parameters for segments and the combined channel were determined. For the segments, the results show that the Weibull distribution fits best for and magnitudes and lognormal for and magnitudes. The generalized Pareto distribution fits best for phases of all ABCD parameters. For the overall channel, the results show that the lognormal distribution fits best for the magnitudes. Generalized extreme value fits best for and phases and logistic for and phases. Based on empirical modeling, it has been observed that the lognormal distribution provides the best fit for the magnitudes , , , and in various dynamic states experienced by the human body. This finding serves to reinforce the validity and accuracy of the FEM model used in the study. Consequently, the dynamic channel was characterized by FEM modeling, empirical analysis, and ABCD network parameters.

References

• 1

  AgacheP. G.MonneurC.LevequeJ. L.De RigalJ. (1980). Mechanical properties and Young’s modulus of human skin in vivo. Arch. Dermatol Res.269, 221–232. 10.1007/bf00406415

  • CrossRef
  • Google Scholar

• 2

  AhmedD.FischerG.KirchnerJ. (2019). “Simulation-based models of the galvanic coupling intra-body communication,” in 2019 IEEE sensors applications symposium (SAS), 1–6. 10.1109/SAS.2019.8706007

  • CrossRef
  • Google Scholar

• 3

  AlbulbulA.ChanA. D. (2012). “Electrode-skin impedance changes due to an externally applied force,” in Medical measurements and applications proceedings (MeMeA), 2012 IEEE international symposium onIEEE, 1–4. 10.1109/MeMeA.2012.6226628

  • CrossRef
  • Google Scholar

• 4

  AliA.InoueK.ShalabyA.SayedM. S.AhmedS. M. (2019). Efficient autoencoder-based human body communication transceiver for wban. IEEE Access7, 117196–117205. 10.1109/ACCESS.2019.2936796

  • CrossRef
  • Google Scholar

• 5

  CallejonM. A.Naranjo-HernandezD.Reina-TosinaJ.RoaL. M. (2012a). Distributed circuit modeling of galvanic and capacitive coupling for intrabody communication. IEEE Trans. Biomed. Eng.59, 3263–3269. 10.1109/tbme.2012.2205382

  • CrossRef
  • Google Scholar

• 6

  CallejonM. A.Reina-TosinaJ.Naranjo-HernándezD.RoaL. M. (2014). Galvanic coupling transmission in intrabody communication: a finite element approach. IEEE Trans. Biomed. Eng.61, 775–783. 10.1109/TBME.2013.2289946

  • CrossRef
  • Google Scholar

• 7

  AndersenJ.-H.BjerkeO.BlakajF.FlugsrudV. M.JacobsenF. A.JonssonM.et al (2019). Bioimpedance and nir for non-invasive assessment of blood glucose. J. Electr. Bioimpedance10, 133–138. 10.2478/.joeb-2019-0019

  • CrossRef
  • Google Scholar

• 8

  BaeJ.ChoH.SongK.LeeH.YooH.-J. (2012a). The signal transmission mechanism on the surface of human body for body channel communication. IEEE Trans. Microw. Theory Tech.60, 582–593. 10.1109/TMTT.2011.2178857

  • CrossRef
  • Google Scholar

• 9

  BaeJ.SongK.LeeH.ChoH.YooH.-J. (2012b). A low-energy crystal-less double-fsk sensor node transceiver for wireless body-area network. IEEE J. Solid-State Circuits47, 2678–2692. 10.1109/jssc.2012.2211654

  • CrossRef
  • Google Scholar

• 10

  BlattenbergerK. (2023). S, h, y, z impedances - rf cafe. Available online at: https://www.rfcafe.com/references/electrical/s-h-y-z.htm.Posted: : July 25, 2023.

  • Google Scholar

• 11

  BouaziziA.ZaibiG.SametM.KachouriA. (2017). “Wireless body area network for e-health applications: overview,” in Smart, Monitored and controlled cities (SM2C), 2017 international Conference on (IEEE), 64–68.

  • Google Scholar

• 12

  CallejonM. A.RoaL. M.Reina-TosinaJ.Naranjo-HernandezD. (2012b). Study of attenuation and dispersion through the skin in intrabody communications systems. IEEE Trans. Inf. Technol. Biomed.16, 159–165. 10.1109/TITB.2011.2171702

  • CrossRef
  • Google Scholar

• 13

  ChawlaA.MukherjeeS.KarthikeyanB. (2006). Mechanical properties of soft tissues in the human chest, abdomen and upper extremities. Mech. Eng. Available at: https://web.iitd.ac.in/∼achawla/PDF%20Files/material_properties%20of%20Soft%20Tissue.pdf (Accessed Oct 2024).

  • Google Scholar

• 14

  ChenX. M.PunS. H.ZhaoJ. F.MakP. U.LiangB. D.VaiM. I. (2016). Effects of human limb gestures on galvanic coupling intra-body communication for advanced healthcare system. Biomed. Eng. OnLine15, 60. 10.1186/s12938-016-0192-z

  • CrossRef
  • Google Scholar

• 15

  ChoH.LeeH.BaeJ.YooH.-J. (2014). “A 5.2mw ieee 802.15.6 hbc standard compatible transceiver with power efficient delay-locked-loop based bpsk demodulator,” in 2014 IEEE asian solid-state circuits conference (A-SSCC), 297–300. 10.1109/ASSCC.2014.7008919

  • CrossRef
  • Google Scholar

• 16

  ClinMed International Library (1990). Young’s modulus and Poisson’s ratio for the tissues. Wilmington, DE: ClinMed International Library.

  • Google Scholar

• 17

  ComleyK.FleckN. A. (2010). A micromechanical model for the young’s modulus of adipose tissue. Int. J. Solids Struct.47, 2982–2990. 10.1016/j.ijsolstr.2010.07.001

  • CrossRef
  • Google Scholar

• 18

  COMSOL Multiphysics (2019). AC/DC module user’s guide.

  • Google Scholar

• 19

  DattaA.NathM.YangD.SenS. (2021). Advanced biophysical model to capture channel variability for eqs capacitive hbc. IEEE Trans. Biomed. Eng.68, 3435–3446. 10.1109/TBME.2021.3074138

  • CrossRef
  • Google Scholar

• 20

  DeSaixP.BettsG. J.JohnsonE.JohnsonJ. E.OksanaK.KruseD. H.et al (2013). Anatomy and physiology.

  • Google Scholar

• 21

  DhamdhereA.ChenH.KurusingalA.SivaramanV.BurdettA. (2010). “Experiments with wireless sensor networks for real-time athlete monitoring,” in Local computer networks (LCN), 2010 IEEE 35th Conference on (IEEE), 938–945.

  • Google Scholar

• 22

  DraperN. R.SmithH. (1998). Applied regression analysis. Wiley Series in Probability and Statistics, 3rd edn. 10.1002/9781118625590

  • CrossRef
  • Google Scholar

• 23

  EbrahimiA. P. (2009). Mechanical properties of normal and diseased cerebrovascular system. J. Vasc. Interv. Neurol.2, 155–162. Available at: https://pmc.ncbi.nlm.nih.gov/articles/PMC3317338/ (Accessed Oct 2024).

  • Google Scholar

• 24

  GabrielC. (1996). Compilation of the dielectric properties of body tissues at RF and microwave frequencies. Tech. Rep. U.S. Air Force Rep. AFOSR-TR-96, U.S. Air Force. 10.21236/ADA303903

  • CrossRef
  • Google Scholar

• 25

  GabrielS.LauR.GabrielC. (1996). The dielectric properties of biological tissues: iii. parametric models for the dielectric spectrum of tissues. Phys. Med. Biol.41, 2271–2293. 10.1088/0031-9155/41/11/003

  • CrossRef
  • Google Scholar

• 26

  GrovesR. B.CoulmanS. A.BirchallJ. C.EvansS. L. (2012). Quantifying the mechanical properties of human skin to optimise future microneedle device design. Comput. Methods Biomechanics Biomed. Eng.15, 73–82. PMID: 21749225. 10.1080/10255842.2011.596481

  • CrossRef
  • Google Scholar

• 27

  GuoL.-X.ZhangY.-M.ZhangM. (2011). Finite element modeling and modal analysis of the human spine vibration configuration. IEEE Trans. Biomed. Eng.58, 2987–2990. 10.1109/TBME.2011.2160061

  • CrossRef
  • Google Scholar

• 28

  IEEE (2012). Ieee standard for local and metropolitan area networks - part 15.6: wireless body area networks. IEEE Std802 (15.6-2012), 1–271doi. 10.1109/IEEESTD.2012.6161600

  • CrossRef
  • Google Scholar

• 29

  IEEES. A. (2008). IEEE P802.15 working group for wireless personal area networks (WPANs): channel Model for body area network (BAN). Rep. IEEE P802.15-08-0033-04-0006, IEEE. Available at: https://mentor.ieee.org/802.15/file/08/15-08-0780-04-0006-tg6-channel-model.pdf (Accessed Oct 2024).

  • Google Scholar

• 30

  International Commission on Non-Ionizing Radiation Protection ICNIRP (1998). Guidelines for limiting exposure to time-varying electric, magnetic, and electromagnetic fields (up to 300 GHz). International Commission on Non-Ionizing Radiation Protection. Health Phys.74, 494–522. Available at: https://www.icnirp.org/cms/upload/publications/ICNIRPemfgdl.pdf (Accessed Oct 2024).

  • Google Scholar

• 31

  ITIS Foundation (2024). Tissue properties. Zurich, Switzerland: The Foundation for Research on Information Technologies in Society (IT’IS).

  • Google Scholar

• 32

  KarimiA.SeraT.KudoS.NavidbakhshM. (2016). Experimental verification of the healthy and atherosclerotic coronary arteries incompressibility via digital image correlation. Artery Res.16, 1–7. 10.1016/j.artres.2016.08.002

  • CrossRef
  • Google Scholar

• 33

  KearneyS. P.KhanA.DaiZ.RoystonT. J. (2015). Dynamic viscoelastic models of human skin using optical elastography. Phys. Med. Biol.60, 6975–6990. 10.1088/0031-9155/60/17/6975

  • CrossRef
  • Google Scholar

• 34

  KutnerM. H.NachtsheimC. J.NeterJ.LiW. (2005). Applied linear statistical models. 5th edn. McGraw-Hill.

  • Google Scholar

• 35

  LeeK.-S.JungM.-C. (2017). “Total force of pinch and grasp by hand postures,” in Advances in physical ergonomics and human factors. Advances in intelligent systems and computingSpringer International Publishing, 208–212. 10.1007/978-3-319-60825-9_23

  • CrossRef
  • Google Scholar

• 36

  LiangX.BoppartS. A. (2010). Biomechanical properties of in vivo human skin from dynamic optical coherence elastography. IEEE Trans. Biomed. Eng.57, 953–959. 10.1109/TBME.2009.2033464

  • CrossRef
  • Google Scholar

• 37

  LimaK. M. M. E.Costa JúniorJ. F. S.PereiraW. C. d. A.OliveiraL. F. (2018). Assessment of the mechanical properties of the muscle-tendon unit by supersonic shear wave imaging elastography: a review. Ultrasonography37, 3–15. 10.14366/usg.17017

  • CrossRef
  • Google Scholar

• 38

  MaityS.MojabeK.SenS. (2018). Characterization of human body forward path loss and variability effects in voltage-mode hbc. IEEE Microw. Wirel. Components Lett.28, 266–268. 10.1109/lmwc.2018.2800529

  • CrossRef
  • Google Scholar

• 39

  NightingaleK.McAleaveyS.TraheyG. (2003). Shear-wave generation using acoustic radiation force: in vivo and ex vivo results. Ultrasound Med. and Biol.29, 1715–1723. 10.1016/j.ultrasmedbio.2003.08.008

  • CrossRef
  • Google Scholar

• 40

  NoormohammadiR.KhaleghiA.BalasinghamI. (2021). Galvanic impulse wireless communication for biomedical implants. IEEE Access9, 38602–38610. 10.1109/ACCESS.2021.3064206

  • CrossRef
  • Google Scholar

• 41

  NotarosB. M. (2010). Electromagnetics. Upper Saddle River, NJ: Pearson.

  • Google Scholar

• 42

  OgasawaraT.SasakiA.-i.FujiiK.MorimuraH. (2014). Human body communication based on magnetic coupling. IEEE Trans. Antennas Propag.62, 804–813. 10.1109/tap.2013.2292705

  • CrossRef
  • Google Scholar

• 43

  Pailler-MatteiC.BecS.ZahouaniH. (2008). In vivo measurements of the elastic mechanical properties of human skin by indentation tests. Med. Eng. and Phys.30, 599–606. 10.1016/j.medengphy.2007.06.011

  • CrossRef
  • Google Scholar

• 44

  PawlaczykM.LelonkiewiczM.WieczorowskiM. (2013). Age-dependent biomechanical properties of the skin. Adv. Dermatology Allergology5, 302–306. 10.5114/pdia.2013.38359

  • CrossRef
  • Google Scholar

• 45

  PereiraM. D.Alvarez-BoteroG. A.Rangel de SousaF. (2015). Characterization and modeling of the capacitive hbc channel. IEEE Trans. Instrum. Meas.64, 2626–2635. 10.1109/tim.2015.2420391

  • CrossRef
  • Google Scholar

• 46

  PersadJ.RockeS. (2020). Investigating the impact of deformation on a 3d-printed antenna in biomedical systems. West Indian J. Eng.42. Available at: https://sta.uwi.edu/eng/wije/vol4202_jan2020/Documents/M07-JPersadJan2020.pdf (Accessed Oct 2024).

  • Google Scholar

• 47

  RockeS.PersadJ. (2015). “Analysis of magnetically-coupled human body communications,” in 2015 COMSOL conference in Boston.

  • Google Scholar

• 48

  RoopnarineA.RockeS. (2024). “Deep learning-based bioimpedance spectroscopy using start of frame delimiter in human body communications,” in Artificial intelligence for sustainable energy. Editors MathewJ.GopalL.JuwonoF. H. (Singapore: Springer Nature Singapore), 209–219.

  • Google Scholar

• 49

  RoopnarineA.RockeS. A. (2021). Stochastic bioimpedance-based channel model of the human body for galvanic coupling. J. Electr. Bioimpedance12, 117–124. 10.2478/joeb-2021-0014

  • CrossRef
  • Google Scholar

• 50

  RoopnarineA.RockeS. A.MahadeoM.PersadJ. (2024). “Multiphysics finite element investigation of galvanic transmission for dynamic livestock-based communication,” in Proceedings of the IEEE conference on agri-food electronics (CAFE).

  • Google Scholar

• 51

  SaundersS.Aragón-ZavalaA. (2007). Antennas and propagation for wireless communication systems. John Wiley and Sons.

  • Google Scholar

• 52

  SeyediM. H.LaiD. (2014). A novel intrabody communication transceiver for biomedical applications. Melbourne, VIC: Victoria University.

  • Google Scholar

• 53

  SmithJ. M. (2011). The doctor will see you always. iEEE Spectr.48, 56–62. 10.1109/mspec.2011.6027250

  • CrossRef
  • Google Scholar

• 54

  SteenbrinkF.MeskersC. G. M.van VlietB.SlamanJ.VeegerH. E. J.De GrootJ. H. (2009). Arm load magnitude affects selective shoulder muscle activation. Med. and Biol. Eng. and Comput.47, 565–572. 10.1007/s11517-009-0482-8

  • CrossRef
  • Google Scholar

• 55

  SujayaB.BhanuPrashanthS. (2021). An efficient hardware-based human body communication transceiver architecture for wban applications. Glob. Transitions Proc.2, 152–156. 10.1016/j.gltp.2021.08.070

  • CrossRef
  • Google Scholar

• 56

  TakamatsuR.HiguchiK.MuramatsuD. (2021). “Measurement frequency evaluation for bioimpedance-based blood-glucose estimation,” in 2021 IEEE 3rd global conference on life sciences and technologies (LifeTech), 309–310. 10.1109/LifeTech52111.2021.9391845

  • CrossRef
  • Google Scholar

• 57

  TakazaM.MoermanK. M.GindreJ.LyonsG.SimmsC. K. (2013). The anisotropic mechanical behaviour of passive skeletal muscle tissue subjected to large tensile strain. J. Mech. Behav. Biomed. Mater.17, 209–220. 10.1016/j.jmbbm.2012.09.001

  • CrossRef
  • Google Scholar

• 58

  TangC.SongY.ChenY.ZhouY.ZhangX.MaY. (2023). “Dynamic channel modeling and characteristic analysis of wearable human body communication system,” in 2023 4th information communication technologies conference (ICTC), 49–53. 10.1109/ICTC57116.2023.10154827

  • CrossRef
  • Google Scholar

• 59

  TangW.ChenN.ZhangJ.ChenS.GeS.ZhuH.et al (2015). Characterization of tactile perception and optimal exploration movement. Tribol. Lett.58, 28. 10.1007/s11249-015-0507-4

  • CrossRef
  • Google Scholar

• 60

  TronstadC.Strand-AmundsenR. (2019). Possibilities in the application of machine learning on bioimpedance time-series. J. Electr. Bioimpedance10, 24–33. 10.2478/joeb-2019-0004

  • CrossRef
  • Google Scholar

• 61

  TsaiB.XueH.BirgerssonE.OllmarS.BirgerssonU. (2019). Dielectrical properties of living epidermis and dermis in the frequency range from 1 khz to 1 mhz. J. Electr. Bioimpedance10, 14–23. 10.2478/joeb-2019-0003

  • CrossRef
  • Google Scholar

• 62

  UlabyF. T.MichielssenE.RavaioliU. (2010). Fundamentals of applied electromagnetics. 6 edn. Upper Saddle River, NJ: Pearson.

  • Google Scholar

• 63

  VizzielloA.SavazziP.GuerraR. R.Dell’AcquaF. (2024). Experimental channel characterization of human body communication based on measured impulse response. IEEE Trans. Commun.72, 3970–3984. 10.1109/TCOMM.2024.3370468

  • CrossRef
  • Google Scholar

• 64

  WangH.TangX.ChoyC. S.SobelmanG. E. (2016). Cascaded network body channel model for intrabody communication. IEEE J. Biomed. Health Inf.20, 1044–1052. 10.1109/JBHI.2015.2448111

  • CrossRef
  • Google Scholar

• 65

  WangK.ChenY.HuangS.WangL.NiuW. (2021). Subcutaneous fat thickness remarkably influences contact pressure and load distribution of buttock in seated posture. J. Healthc. Eng.2021, 1–9. 10.1155/2021/4496416

  • CrossRef
  • Google Scholar

• 66

  WegmuellerM. S.KuhnA.FroehlichJ.OberleM.FelberN.KusterN.et al (2007). An attempt to model the human body as a communication channel. IEEE Trans. Biomed. Eng.54, 1851–1857. 10.1109/TBME.2007.893498

  • CrossRef
  • Google Scholar

• 67

  WellsR.GreigM. (2001). Characterizing human hand prehensile strength by force and moment wrench. Ergonomics44, 1392–1402. 10.1080/00140130110109702

  • CrossRef
  • Google Scholar

• 68

  WenE.SievenpiperD. F.MercierP. P. (2022). Channel characterization of magnetic human body communication. IEEE Trans. Biomed. Eng.69, 569–579. 10.1109/TBME.2021.3101766

  • CrossRef
  • Google Scholar

• 69

  WhiteE. A.OrazemM. E.BungeA. L. (2013). Characterization of damaged skin by impedance spectroscopy: mechanical damage. Pharm. Res.30, 2036–2049. 10.1007/s11095-013-1052-1

  • CrossRef
  • Google Scholar

• 70

  ZhangY.HeZ.LiuY.EnamoradoL. A. L.ChenX. (2016a). “Measurement and characterization on a human body communication channel,” in 2016 IEEE 27th annual international symposium on personal, indoor, and mobile radio communications (PIMRC)IEEE, 1–6.

  • Google Scholar

• 71

  ZhangY.KouB.FanD.LiuY.HeZ.ChenX. (2016b). “A dynamic pilot interval adjustment scheme for hbc channel estimation,” in Communications in China (ICCC workshops), 2016 IEEE/CIC international Conference on (IEEE), 1–5.

  • Google Scholar