Studies of Scattering, Reflectivity, and Transmitivity in WBAN Channel: Feasibility of Using UWB

The Wireless Personal Area Network (WPAN) is one of the fledging paradigms that the next generation of wireless systems is sprouting towards. Among them, a more specific category is the Wireless Body Area Network (WBAN) used for health monitoring. On the other hand, Ultra-Wideband (UWB) comes with a number of desirable features at the physical layer for wireless communications. One big challenge in adoption of UWB in WBAN is the fact that signals get attenuated exponentially. Due to the intrinsic structural complexity in human body, electromagnetic waves show a profound variation during propagation through it. The reflection and transmission coefficients of human body are highly dependent upon the dielectric constants as well as upon the frequency. The difference in structural materials such as fat, muscles and blood essentially makes electromagnetic wave attenuation to be different along the way. Thus, a complete characterization of body channel is a challenging task. The connection between attenuation and frequency of the signal makes the investigation of UWB in WBAN an interesting proposition. In this paper, we study analytically the impact of body channels on electromagnetic signal propagation with reference to UWB. In the process, scattering, reflectivity and transmitivity have been addressed with analysis of approximate layer-wise modeling, and with numerical depictions. Pulses with Gaussian profile have been employed in our analysis. It shows that, under reasonable practical approximations, the human body channel can be modeled in layers so as to have the effects of total reflections or total transmissions in certain frequency bands. This could help decide such design issues as antenna characteristics of implant devices for WBAN employing UWB.


Wave Propagation through Biological Media
From the appendix A, if the incident wave is a linearly polarized uniform plane wave travelling along the z-direction, then, for E and H, Equation (A.9) is of the form:

Scattering, Reflection, and Anti-Reflection
Sensor nodes find a human body, when they are placed, to be layered media. Fat, muscles etc. are such independent layers. Signals from wireless sensors in human body experience scattering, reflection, and diffraction by those layers. In Appendix B, we started with acoustic wave for mathematical formulation and then used these formulae for electromagnetic wave. We have established the basic equations for the wave propagation through the layered media. Now we will use

Variation of Penetration Depth with Frequency
Muscle Fat those mathematical equations for calculating some wave properties such as scattering, reflection, and transmission while passing through layered media.

Scattering by a Single Interface
In this section we consider the case in which two homogeneous half-spaces are separated by an interface at z = 0 ( Figure 3): The goal of this section is to analyze the scattering problem in terms of right-and left-going modes (Appendix B).
We introduce the local velocities = and impedances = and the right-and left-going modes defined by: For j = 0, 1, the pairs (A j , B j ) satisfy the following system in their respective half-spaces: which means that ( , ) is a function of − only, and ( , ) is a function of + only. We assume that a right-going wave with the time profile f is incoming from the left and is partly reflected by the interface. We also assume a radiation condition in the right half-space so that no wave is coming from the right. Assume that f is completely supported in 0, ∞ . We next introduce two ways to define proper boundary conditions: (I) We can consider an initial value problem with initial conditions given at some time 0 < 0 by: As shown in the Appendix B, these initial conditions generate a pure right-going wave whose support at t = t 0 is in the interval ∈ −∞, 0 0 , which lies in the left half-space.
(II) We can consider a point source located at some point 0 < 0 and generating a forcing term of the from: As seen in the Appendix B, this point source generates two waves. The left-going wave is propagating into the negative z-direction and will never interact with the interface, so we will ignore it. The right-going wave first propagates in the homogeneous left half-space and it eventually interacts with the interface z = 0.
In terms of the right-and left-going waves, these two formulations give the same descriptions. We have 0 , = − 0 < 0, and B 1 (t, z) = 0 for > 0, and consequently, at the interface z = 0: Note that the delays introduced in the initial conditions (10) and in the forcing term (11) have been chosen so that the boundary conditions (12) have a very simple form.
The pairs (A 0 , B 0 ) and (A 1 , B 1 ) are coupled by the jump conditions at z = 0 corresponding to the continuity of the velocity and pressure fields: which gives: with (±) = 1 2 1 0 ± 0 1 . Note that (+) 2 − (−) 2 = 1. The matrix can be interpreted as a propagator, since it -propagates‖ the rightand left-going modes from the left side of the interface to the right side. Such a propagator matrix will be called interface propagator in the following.
Taking into account the boundary conditions (12) yields: and solving this equation gives: where ℛ and are the reflection and transmission coefficients of the interface: These coefficients satisfy the energy-conservation relation: meaning that the sum of the energies of the reflected and transmitted waves is equal to the energy of the incoming waves. Finally, the complete solution for < 0 in terms of the right-and left-going modes is: and for > 0: Using (7-8) we can obtain the pressure and velocity fields ( Figure 4).

Single-Layer Case: Scattering
In this section, we consider the case of a homogeneous slab with thickness L embedded between two homogeneous half-spaces ( Figure 5). Three regions can be described as follows: with j = 0 for z < 0, j = 1 for ∈ 0, , and j = 2 for z = L. The boundary conditions correspond to an impinging pulse at the interface z = 0 and a radiation condition at z = L 2 : The propagation equations (9) in each homogeneous region show that is a function of − only and is a function of + only. The waves inside the slab [0, L] are therefore of the form: 1 , = 1 ( − 1 ), 1 , = 1 + 1 while the reflected wave for z <0 is of the form: and the transmitted wave for z > L is of the form: We want to indentify the functions b 0 and a 2 , which give the shapes of the reflected and transmitted waves. with 1 (±) = 1 2 2 1 ± 1 2 . We can write these relations in terms of the functions a j , b j as: The second equation of (14) can be written as: The syplectic matrix 1 is a propagator in the frequency domain. It propagates the right-and left-going modes from the right side of the interface 0 to the right side of the interface 1, and it depends on the layer thickness L. Finally, combining the first equation of (14) and (15), we obtain the relation: where the frequency-dependent syplectic matrix: is the overall propagator of the slab. Equation (16) shows that K 0 propagates the right-and left-going modes from the left side of the interface 0 to the right side of the interface 1. We find explicitly: By solving equation (16), whose unknowns are 2 and 0 and using the expressions of (±) , we obtain: where the frequency-dependent reflection and transmission coefficients are: which means that the individual energies of the frequency components of the incoming pulse are preserved by the scattering process. The main qualitative difference between the scattering by a single interface and the scattering by single layer is that the reflection and transmission coefficients in the layer case are frequency-dependent. This frequency dependence originates from interference effects between the waves that are scattered back and forth by the two interfaces of the layer.

Reflection
Let us consider a layer embedded between two homogeneous half-spaces that have the same material properties, i.e., the situation in which ρ 2 = ρ 0 and K 2 = K 0 . We then have R 1 = -R 0 and T 1 = T 0 , which implies that the global reflectivity of the layer can be written as: The reflectivity is periodic with respect to the angular frequency ω with the period ω c = πc 1 /L. As a function of the angular frequency the reflectivity goes from the minimal value: to the maximal value: This shows that for any value of the reflection coefficient R 0 of a single interface, there exist frequencies that are fully transmitted or fully reflected by the layer. If we consider the case of strong scattering 0 2 ≪ 1, then the transmitted frequency bands have a width of the order of

Anti-Reflection
The total transmission phenomenon is also encountered in situations in which the two half-spaces are different. Indeed, consideration of human body part as an ideal, fully-transmitting layer is certainly beyond perfection. In a microscopic or constituent-wise sense, a human body-part, striated muscle for example, iscomposed of water (70.09%), ether-soluble extract (6.60%), crude protein (21.94%), etc. [35]. When one is faced with the task of modeling portions of human body as a medium for electromagnetic (signal) propagation, there are inevitable practical assumptions and approximations to be made. Thinner layers can be considered with more homogeneous characteristics, while for thicker setting with internal variability, the aggregate behavior sums up by and large. Thus a certain part of a human body like muscle or fat, when considered in macroscopic perspective, can be regarded as a continuum, and hence the idea of a plannerlayered-medium assumption of human body tissues contextually holds for all practical modeling considerations. Such model could greatly affect related system design pertaining to crucial parameters (antenna and others); for instance, in addition to the frequency and bandwidth employed for the signal, the measurement of the depth of the body at which the transceiver of an implant device would function with acceptable accuracy has a lot to do with such parameters as permittivity, permeability, and impedance of the intermittent layers of body tissues.  shows a grossly approximated propagation system consisting of three layers: air, human body channel, and a transceiver. Admittedly, finding fully-transmitting or fully-reflecting layers in body-parts is unrealistic, but it is widely adopted practice to model systems using phantoms comprised of components having somewhat homogeneous characteristics, yet closely resembles body-parts in aggregate behavior. In some recent works ( [33] and [34]) UWB antenna impedance matching has been studied in the context of biomedical implants. We assume that the two homogeneous half spaces ( Figure 7) have different impedances ζ 0 ≠ ζ 2 , then it is possible to choose the thickness L and the impedance ζ 1 of the layer so that a given frequency ω will be fully transmitted from one half-space to the other one, which would not be the case in absence of such a layer. From the analysis of the reflectivity function: One can show that a necessary and sufficient condition for ℛ 2 to be zero is that In the case ζ 0 ≠ ζ 2 this in turn enforces one to choose the impedance of the layer to be 1 = 0 2 (so that R 0 = R 1 ) and the thickness L to be chosen so that ωL/(πc 1 ) is half an integer (so that cos(2 ωL/c 1 ) = -1). Usually the thickness is chosen to be equal to a quarter of the wavelength, meaning ωL/(πc 1 ) = ½.

Path Loss in the Human Body (Near Field Far Field Consideration)
When EM RF waves propagate in freespace, the power received decreases at a rate of 1 , n being the coefficient of pathloss. Other kinds of losses would be fading of signals due to multipath propagation. However, for propagation of EM waves in a lossy medium like human tissue, the losses would be mainly due to absorption of power in the tissue, where it is dissipated as heat. As the tissue medium is lossy and mostly consists of water, the EM waves are attenuated considerably before they reach the receiver. The Specific Absorption Rate (SAR) is useful in determining the amount of power lost due to heat dissipation. SAR is defined as power absorbed per unit mass of the tissue [29]. SAR is a standard measure of how much power is absorbed in the tissue and depends upon E-and H-field strengths. By determining the average SAR over the entire mass of the tissue between the transmitter and the receiver, we are able to compute the total power lost. SAR in the near field of the transmitting antenna depends mainly on the H-field, whereas the SAR in the far field of the transmitting antenna depends mainly on the E-field. We use Maxwell's E-and H-fields equations for lossy medium to obtain the average SAR of the medium between the transmitting and the receiving antenna in the far field and near field, respectively. WBAN applications involve wireless communications between implanted biosensor nodes inside human body.
These nodes exchange data among themselves and also with the base-station. In general, the system model consists of numerous biosensor nodes placed inside the various parts of the human body surrounded by tissues. In particular, for the development of this model, we consider only one transmitting and one receiving antenna separated by a distance d. An elemental short dipole (dipole length_wavelength) in a lossy human tissue medium is considered for this purpose [28] and is shown in Figure 9. A small area of tissue surrounding the antenna is considered for our analysis. Thus we can safely assume the human tissue under consideration to be a homogeneous medium with no sharp edges, no rough surfaces and having uniform electric and magnetic properties. The received power is assumed to be due only to the power from the transmitter and not from any other source. The space around the radiating antenna is divided into near field and far field regions as shown in Figure 10.  The region of space immediately surrounding the antenna is known as the near field region. The extent of the near field in the case of short dipoles is given by d 0 = λ/2, where λ is the wavelength [30]. In the near field, the E-and H-field strengths vary rapidly with the distance from the antenna. The far field is the entire region beyond the near field. In the far field region, the E-and H-field exhibit a plane wave behavior. Power absorbed between the transmitting and receiving antennas can be considered as the sum of power absorbed in near field (P NF ) and far field (P FF ) regions. The total power absorbed between the two antennas is computed by numerical integration.
Consider an elemental oscillating electric dipole in a lossy medium of conductivity σ (S/m), permittivity є (F/m), permeability μ (H/m), complex propagation constant γ, complex intrinsic impedance = + [28] at frequency ω, as shown in Figure 8. The dipole consists of a short conducting wire of length , terminated in two small conductive spheres or disks. Assume that the current I is uniform and varies sinusoidally with time [28]. The electromagnetic field at a distance ‗R' for an Hertzian dipole is derived from the vector potential A, given by [28]: is the unit vector in the z-direction, is the propagation constant, given by = + ; attenuation constant and phase constant is given by as [28]:

, a R A R + a θ A θ + a φ A φ ) are given by A R = A z cosθ, A θ = -A z singθ and
A φ = 0. The magnetic field intensity H and the electric field intensity E is given by [28]: Solving the above magnetic and electric field equations for lossy medium and expressing in terms of complex impedance we get:

Power Absorbed in the Near Field
The SAR in the near field is given by [31]: where is the density of the medium and c corr is the correction factor to take into account the changed reflection properties for small distances R of the antenna from the scatterer. Since we assume both the transmitting and receiving antennae are in a same homogeneous medium, the plane wave reflection coefficient is zero. By substituting = 0 and RMS value of the H-field, the above equation reduces to: which gives the value of SAR at a point at distance ‗R' and angle ‗θ' from the dipole. Power at infinitely small volume = 2 is: The power absorbed in the near field of the lossy tissue can be obtained by computing the average SAR over the entire tissue mass in the near field, which is obtained by integrating ∆ over the entire mass in the near field region, i.e., from the surface of the antenna (R = r) to the end of the near-field region (R = d 0 ): where: The antenna dimensions depend on the wavelength of the wave in the medium given by = 2 [28].

Power Absorbed in the Far Field
Neglecting 1 2 , 1 3 … .. terms from field Equations (21)-(23) for the far field, we have: In the far field the specific absorption rate depends only on the value which is given by [29]: The power absorbed in the infinitely small volume = 2 in the far field, at a distance R and angle from the dipole can again be obtained from (24): The total power absorbed in the far field of the lossy tissue between the source and destination antennas can be obtained by computing the average SAR over the entire tissue mass in the far field from distance d 0 to d (d 0 is the point where the far field starts). This is obtained by integrating ∆ over the mass in the far field between the two antennas: Power received The effective radiated power (ERP) is obtained by subtracting the loss in the near field (P NF ) and far field (P FF between the transmitting and receiving antennas) from the transmitted power P T (i.e., (P T -P Loss )G t ), where P Loss = P NF +P FF is obtained from (25) and (26). The power density (P e , Power per unit area) at a distance 'd' is different in near field and far field regions: (a) P R in the Near Field: There is no general formula for the estimation of field strength in the near field zone [30]. Only measurements can provide a simple means of field evaluation. However, reasonable calculations can be made for antennas like dipole or monopole. When the receiving antenna is in the near field region of the transmitting antenna, the power density does not necessarily depend on the distance from the antenna, but varies rapidly with distance, and may exhibit oscillatory behavior. The magnitude of on-axis (main beam) power density varies according to the location in the near field and its maximum value is approximated by [32] = 16 / 2 , where L is the largest dimension of the antenna, P is P T -P NF , and is the aperture efficiency (typically 0.5-0.75) [32]. It can be approximated as = / (A e is the effective aperture and A is the physical area of the antenna). The power received by the receiving antenna in the near field can be approximated by: (b) P R in the Far Field: On the other hand when the receiving antenna is in the far field region of the transmitting antenna, the power density is dependent on the distance d and is given by: The power received by the receiving antenna in the far field is = , where the receiving antenna aperture A e is given by = Here, G t and G r are the gain of the transmitting and receiving antenna, respectively. Thus the received power is: and a total phase change of − is involved during the propagation of the wave. Thus, PMBA can be used for calculating the propagation loss using the two Equations (25) and (26).

Numerical Results
The permittivity of biological tissues depends on the type of tissues (e.g. skin, fat, or muscle), water content, temperature, and frequency. However, the permittivity and frequency may also determine how far the EM wave penetrates into the body. The term depth of penetration (D p ) usually quantifies this. It is observed from Equations (1) and (2) that the wave gets attenuated as it propagates in the biological material along the z-axis. As shown in Figure 1, variation of radiation power density has been compiled for four different frequencies (27 MHz,100 MHz, 433 MHz, and 1,500 MHz) with respect to the depth in muscle. At a given depth, usage of lower frequency results in a higher power density as illustrated in Figure 1. We discovered a distinct feature (demonstrated by Figure 2) which states that attenuation is more in fat than that in muscle with respect to frequency. Expressions for left-and right-going modes of a pulse have been derived in appendix B. Equations (19) and (20) express the reflectivity and transmitivity of a layer respectively. Equations (7) and (8) are the expressions of the left-and right-going modes respectively of a pulse when scattered by a single interface between two layers. Figure4. shows the scattering of a Gaussian pulse by an interface separating two homogeneous half-spaces (c 0 , ζ 0 , z < 0) and (c 1 , ζ 1 , z > 0). The spatial profiles of the velocity field (a) and of the pressure field (b) are plotted at times t = −4, t = −3,..., t = 6. Reflectivity ℛ 2 versus frequency has been depicted by Figure 6. We can see here that, at a certain period ω c = πc 1 /L, the reflectivity maximum which is almost 0.04 for a layer with R 0 = −R 1 = 0.1 (a) and almost 1.0 for a layer with R 0 = −R 1 = 0.9 (b). In Figure 8, Transmitivity 2 versus frequency curves have been drawn.
Transmitivity has been found to vary periodically with a certain frequency. The period of this frequency band depends upon the choice of the layer thickness L. For a layer with R 0 = −R 1 = 0.1, transmitivity is about 1.0 (a) and for a layer with R 0 = −R 1 = 0.9, transmitivity is about 0.9 (b). Here the period is ω c = πc 1 /2L. From Equation (20), we can see that by the proper choice of the impedance of the layer to be 1 = 0 2 (so that R 0 = R 1 ) and the thickness L to be chosen so that ωL/(πc 1 ) is half an integer (so that cos(2 ωL/c 1 ) = -1), we can form a fully transmitting layer. Usually the thickness is chosen to be equal to a quarter of the wavelength, meaning ωL/(πc 1 ) = ½. Therefore, from the results shown above, we can infer that a layer can either fully reflect or fully transmit any incoming wave at a certain frequency or frequency band. We can use these results to UWB by proper choice of impedance and the thickness L. Power loss in near filed and far field due to absorption has also been analyzed. A propagation loss model (PMBA) for homogeneous tissue bodies has been presented, which compares PMBA with the freespace propagation model ( = 4 ; with loss coefficient n = 3). A frequency range of (900 MHz to 3 GHz) has been considered here. We have been able to make a conclusion that, compared to freespace, there is an additional 30-35 dB of attenuation at small distances the far field ( Figure 11). This loss increases further with the distance and frequency. It is argued that the human body cannot be considered as layered a media. However, we have been able to homogenize the human body channel, which is shown in Appendix C.

Conclusions
Employing UWB in WBAN involves a lot of promise, just as there are a number of relevant challenges. We studied the technical feasibility in this regard with a concentration in electromagnetic propagation of the signal across human body. Unlike conventional wireless channels, human body comes with a great deal of structural complexity requiring significantly different design considerations. The reflection and transmission coefficients of human body are heavily dependent upon the dielectric constants as well as upon the frequency. In this work, we investigated a layer-wise model for electromagnetic propagation across the components of the body in regard to such key aspects as scattering, reflectivity, and transmitivity. Naturally, the segmentation in precise layers are not what we come across in a body. But, the approximate model employing homogenization could help assess the aggregate behavior of the wireless communication involving implant devices, thus guiding the potential design issues for antenna characteristics, for instance. We also presented numerical depictions of some of the pertinent signal characteristics. From here on, one could expect to further improve the model in terms of suitable layering and other parameters of approximations.

Appendix B
In this appendix, we use a number of essential transformations of the wave equation that are specific to layered media. We consider the particular case in which the parameters of the medium vary in a piecewise-constant manner; in other words, we consider a stack of layers made of homogeneous media. We study the propagation of a normally incident plane wave, which enables us to reduce the problem to the one-dimensional acoustic wave equations [9]. We will see that the problem can be recast as a product of matrices corresponding to the scattering of the wave by the successive interfaces between the layers. This is a classical setup for waves propagating in this particular type of layered media, and it is extremely useful for direct numerical simulations. The equations for the three-dimensional velocity u and pressure p are: where ρ is the density of the medium and K the bulk modulus of the medium. These two equations correspond respectively to conservation of momentum and mass. The density and bulk modulus are assumed to be spatially varying along the z-coordinate. If the initial conditions correspond to a plane wave that is normally incident to the layered medium, then the solution of the equations remains independent of the transverse variables, and the transverse velocity is zero. The system can then be reduced to the one-dimensional wave equations. It should be noted that more general conditions, corresponding in particular to point source, require a more general three-dimensional framework. In this appendix, we focus our attention to the one-dimensional case.
In one-dimensional medium the equations for the velocity u and pressure p are: In this representation the material parameters ρ and K may either vary or remain constant (special case) with respect to the space coordinate z.

Right-and Left-Going Waves
We consider the special case with a homogeneous medium in which the coefficients ρ and K are constant. Consequently the speed of sound c and the impedance ζ are constant, and the system can be written: and the waves can be written , = − and , = + for some wave-shape functions a and b. Thus, in the constant medium case we have decomposed the wave into the right-and left-going waves A and B, which do not interact.
To fully specify the problem we have to prescribe initial conditions, for instance the velocity and pressure profiles at time t = 0: = 0, = 0 , = 0, = 0 We then translate these initial conditions for u and p into initial conditions for the modes A and B: