Finite Element Simulation of the Impedance Response of a Vascular Segment as a Function of Changes in Electrode Configuration

Abstract Monitoring a biological tissue as a three dimensional (3D) model is of high importance. Both the measurement technique and the measuring electrode play substantial roles in providing accurate 3D measurements. Bioimpedance spectroscopy has proven to be a noninvasive method providing the possibility of monitoring a 3D construct in a real time manner. On the other hand, advances in electrode fabrication has made it possible to use flexible electrodes with different configurations, which makes 3D measurements possible. However, designing an experimental measurement set-up for monitoring a 3D construct can be costly and time consuming and would require many tissue models. Finite element modeling methods provide a simple alternative for studying the performance of the electrode and the measurement set-up before starting with the experimental measurements. Therefore, in this study we employed the COMSOL Multiphysics finite element modeling method for simulating the effects of changing the electrode configuration on the impedance spectroscopy measurements of a venous segment. For this purpose, the simulations were performed for models with different electrode configurations. The simulation results provided us with the possibility of finding the optimal electrode configuration including the geometry, number and dimensions of the electrodes, which can be later employed in the experimental measurement set-up.


Introduction
Importance of monitoring biological tissues as a three dimensional model Two dimensional (2D) cell cultures have been used as in vitro cell models for decades, however; a traditional two dimensional cell culture lacks the cell-cell communications in three dimensions as well as influential proteins that can be found only in a three dimensional (3D) extracellular matrix. In this way a 2D cell culture fails to mimic in vivo tissue complexity (1). Whereas, the cell-cell communication in a 3D environment in which the cells are both connected to each other and the extracellular matrix, forms a complex and dynamic system which adds a new dimension to the in vitro situations (2). 3D tissue models due to their close resemblance to an in vivo environment, make it possible to study the physiological responses of the cells in their native condition such as exchanging oxygen, nutrients and waste as well as monitoring their structural features (1,3). Moreover, the 2D and 3D cell cultures differ not only quantitatively due to their dimensions but also qualitatively due to the cell behaviors such as: cellular morphology, proliferation rates, cell-extracellular matrix interactions (4), migration (5), gene expression, differentiation, signaling (6), physiological function, and electrophysiological properties (2,7). Therefore, it is very important to monitor biological tissues as a three dimensional construct to be able to properly characterize their physiological and morphological properties and behaviors (7).
Monitoring biological tissues using bioimpedance measurement techniques A cell membrane acts as a very thin capacitor with high resistance. Therefore, when an alternating electrical signal is applied, the cell membrane acts as a capacitance. However, intracellular and extracellular fluids provide low resistivity (8)(9)(10). Consequently, the electrical impedance that the tissue produces under an alternating current is a complex quantity, which consists of a real and an imaginary part. The real part (R) is associated with the resistive pathways across the tissue such as intra and extracellular media. The imaginary part, the reactance (X), is associated with capacitive (C) and inductive (L) properties. Biological tissue and cells express capacitive current conduction through different structures as the cell membrane etc. that acts as capacitive pathways for alternating current. At the same time the inductive properties in biological tissues are typically negligibly small at frequencies below 10 MHz, when compared to their resistance and capacitance (11).
As the electrical impedance combines resistance and capacitance, it changes with variations in tissue anatomy, morphology, structure, and composition, as tissue components can demonstrate both conductive and charge storage properties.
Conductivity or current conduction in a biological tissue defines how well the tissue conducts electrical current. Electrical conductivity (σ) comprises of a non-frequency dependent component due to its ion content and their electric field induced mobility. Furthermore, a frequency dependent component associated with the energy absorbed by the dielectric relaxation process or dielectric loss factor, which defines the imaginary component of permittivity (12). On the other hand, permittivity (Ɛ) defines the extent to which the tissue can be polarized under the applied electrical excitation. Relative permittivity (Ɛr) or dielectric constant is the ratio of the permittivity of a tissue to the permittivity of vacuum (13). Both electrical conductivity and relative permittivity vary widely between different biological tissues (14,15).
The current conduction path in the tissue changes by altering the frequency of the applied signal. At low frequencies, cell membrane has a high resistance and the electric current can only travel through the extracellular fluid. However, at high frequencies, the reactance of the cell membrane decreases due to its capacitive behavior. At very high frequencies, the impedance of the membrane becomes very low and the current passes more uniformly through the tissue structure (11). Thus, measuring electrical impedance in different current conduction pathways by varying the signal frequency, makes it possible to analyze different tissue structures, composition and behaviors and study their anatomical, physiological and pathological status (9,16).
Bioimpedance spectroscopy is a technique, which calculates and analyzes the electrical impedance from the boundary voltage/current measurements of the tissue at different frequencies. By employing bioimpedance spectroscopy not only the impedance of the tissue under study but also some other complex bioelectrical properties of the tissue such as dielectric relaxation and dispersions ( , , and ) (17) can be studied. Due to the ability of high frequency AC currents to penetrate the tissue, bioimpedance spectroscopy can be considered as a low cost, noninvasive and automation-compatible, label-free method to characterize the electrical properties of 3D tissue models in real time and with high temporal resolution (3).

Influence of electrode configuration in bioimpedance spectroscopy measurements
Proper electrode configuration is substantial in establishing an optimized measurement set-up with precise performance, and includes electrode geometry, layout, dimension, structure and material (1,3). Due to the importance of the dimensions of the electrodes, in some studies, theoretical calculations are performed to design the shape and size of the electrodes, prior to the experimental measurements (3,(18)(19)(20).
As the electrodes would be in close contact with the tissue under measurement as well as the culturing medium, it is very important that the electrode material would be chosen of a chemically and physically stable material that would not cause any kind of toxicity during the experiment. Some of the commonly used materials with nontoxic properties for this purpose are gold, platinum, palladium, and titanium (21)(22)(23)(24). The other material, which is commonly used due to its low cost and stable potential is silver chloride (25). However, long term contact with biological tissue and culturing medium can cause erosion of the material and loss of silver chloride coating which lead to cytotoxicity due to change in electrode potential (22,25,26). So even by using stabilizing coatings, the silver chloride electrodes are less biocompatible than the electrodes made of gold, platinum and titanium (3).
Optimization of measurement performance by changing the electrode configuration and material has been done by many researchers before. Fosdick and Anderson (27) used the amperometric response of a flow detector for optimization of the design of its microelectrode arrays. Min and Baeumner (28) studied optimization of oxidation and reduction reactions of potassium ferro/ferrihexacyanide by changing the geometry and material of the interdigitated ultra-microelectrode arrays in a measurement system with a two electrode set-up. Sandison et al. (29) showed the dependence of the electrochemical behavior of microelectrode arrays on the center-to-center spacing of the arrays. Lempka and coworkers (30) studied optimization of the signal to noise ratio and capability of silicon based microelectrode arrays for recording neural activities as a function of changing the electrode contact size. Wang and colleagues (31) investigated variations in the sensitivity of interdigitated microelectrode arrays as the result of changes in electrode dimension as well as frequency for electric cell-substrate impedance sensing (1).
A 3D measurement set-up should employ an electrode configuration, which would provide monitoring of the electrical signal propagation in the tissue construct in a three-dimensional manner. Performance of a 3D monitoring electrode is not only dependent on the electrode layout including the shape and size, and arrangement of the electrodes and their material, but also the curvature of the sensing area. A flexible electrode substrate provides the possibility of performing 3D measurements of the tissue while this is not possible with a conventional 2D electrode platform (7).
In order to achieve a solid integration between the electrodes and the 3D tissue, both the electrodes' arrangement and dimensions as well as the mechanical properties of the electrode and its substrate should be compatible with the tissue to be monitored (32). Some researchers have reported using flexible and/or stretchable electrodes and electrode substrates, which due to their flexibility conform to the surface of the tissue to be monitored (33)(34)(35)(36)(37). Some of these flexible electrodes have been used in monitoring the electrical activities near surfaces of the heart (33-35), brain (37) and skin (36).

The role of simulation and defining a finite element model
Designing the optimal measurement set-up requires performing the experiment several times and with various set-ups, which would be costly and time-consuming (38). Simulation using Finite Element Modeling (FEM) is a tool, which would facilitate understanding the tissue behavior and electrical response by visualization of current conduction paths in different parts of the model as well as the possibility of analyzing the contribution of each part of the model to the measurement results from the whole model. It also demonstrates the dependency of impedance and its real and imaginary components on the cell compositions such as cell membrane and intra and extracellular fluids (16). COMSOL Multiphysics is one of the simulation softwares that can demonstrate current conduction paths through 3D tissue models in multi-frequency bioimpedance measurements and in this way be employed as a practical tool in finding the optimized measurement set-up for monitoring 3D tissue models (1).
In a simulation approach for studying the electrical response of a model, the concept of impedance network can be used (39), where the model under simulation is divided into its components, where each of those components have their specific electrical properties such as electrical conductivity, relative permittivity and current density, etc. (40). Then the model components would be polarized using an electric field and then numerical calculation of impedance, capacitance and all other related electrical properties of the model components would be performed (41).
As in impedance measurement set-ups usually small dimensions (mm) and low frequencies (Hz to MHz) are used, therefore, the electrical field wavelength is much larger than the dimensions of the simulated sample and thus in simulating such measurement set-ups, a quasi-static approximation would typically be used (42,43).

Aim of the study
This study has been designed to develop a tool for optimizing a measurement set-up in which a flexible electrode is employed to facilitate 3D monitoring of a venous segment and provide measurements of high precision and accuracy. Therefore, the aim of the current study is employing finite element methods for simulating different measurement set-ups through building models with different electrode configurations (electrode number and dimensions). In addition, the finite element simulation would provide information about the contribution of each part of the model to the impedance measured from the whole system. We would expect that the optimized measurement set-up is obtained in the model in which venous segment has the highest contribution to the total measured impedance.

Methods
The AC/DC module of the COMSOL Multiphysics® platform version 5.4 (COMSOL Inc., Burlington, MA, USA), was employed to perform the finite element simulations for different frequencies (44).

Defining the model geometry and assigning material properties
The 3D model (figure 1) consists of a chamber out of PMMA or Poly methyl methacrylate also known as acrylic or Plexiglas, which is filled with saline solution and two glass tubes (silica) which are mounted in the chamber by passing through the holes fixed in the walls of the chamber. A venous segment in a cylindrical form with the length of 6 cm and diameter of 1.4 cm is mounted between these two glass tubes and the whole tube construct including glass tubes and the vein segment is filled with saline. The venous segment is surrounded by a thin layer of isolating material (quartz) as the flexible electrode substrate, which is rolled around the tissue in a cuff form.
Tiny golden cylinders, which are placed on the inner surface of the cuff towards and in contact with the venous segment are considered the surface electrodes, whereas the golden cylindrical bar passing through the vein and glass tubes is considered the ground electrode in this model. Various numbers of surface electrodes (28, 32, 44 and 48) as well as different dimensions for both type of electrodes (rterminal = 1, 2, 2.5 mm and rground = 0.5, 1 mm) were chosen. These were chosen the way that it would practically conform to the experimental set-up to be used later for measurements on a venous segment with the abovementioned dimensions. Gold was chosen as the electrode material in this model due to its high electrical conductivity, biocompatibility, and low electrode polarization (21)(22)(23)(24).
The simulation system is defined by the dimensions and dielectric properties of the different components of the model (43). Due to the dependency of impedance measurements to permittivity and conductivity, the electrical properties used were frequency dependent electrical conductivity and relative permittivity (17). Therefore, the values for the electrical conductivity and relative permittivity of the different parts of the model were obtained from the inbuilt dielectric property tables for different materials in the COMSOL Multiphysics library in the way that they would match the values of the materials used for the experimental measurements later (38,45).

Defining the electrodes and electrical sources
A two-electrode configuration model with a voltage difference of 1 V was chosen. The same terminal feature was used for all the tiny surface golden electrodes were the voltage was applied, while the wire passing through the glass tubes and the venous segment was chosen as the ground electrode (figure 1) and the current flowing from the terminal electrodes to the ground electrode was calculated. Through Ohm's law, admittance (Y) was obtained and then the impedance (Z) was calculated as its inverse (Z = 1/Y). The impedance calculated this way is the impedance measured between the terminal and the ground electrode.
For our model, two different simulation methods were chosen and performed. In the first method, all the surface electrodes were considered as the terminal electrodes and the excitation was applied to all of them at the same time. While in the second method, one surface electrode at the time was chosen as terminal electrode to which the voltage was applied by turn. In this way, each surface electrode was employed as the terminal once during the simulation process.

Meshing
In this stage, an extra fine mesh (figure 2) was created. This mesh was then merged with electrical properties of the material chosen to create the model (43). Meshing is done for subdividing the model into non-overlapping elements in the three dimensional space. The response of the whole model is the result of assembly of all the elements. Mesh has different levels from coarse to extra fine. If a fine mesh is chosen, the accuracy of the model would be increased. For example, a finer mesh provides a better representation of curved edges due to the smaller elements. However, this way the simulation tool has more data to calculate, therefore more time and memory is needed (38).

Solving by using the frequency-dependent study type
Since the wavelength of the electric field in the different layers is much larger than the dimensions of the electrodes, a quasi-static electrical conduction model was applied to solve for the electric field in the saline solution (15,44). The finite element simulation using the 3D AC/DC module of COMSOL Multiphysics was then performed for impedance measurements in the model. After obtaining the solution for the electric potential, boundary integration was used to determine electric current and consequently obtaining the electric impedance. With the AC/DC module of COMSOL it is possible to obtain the solution for the electric impedance at different frequencies. Simulations were repeated for frequencies exponentially distributed from 10 Hz to 1 MHz, to characterize the frequency response and the bioimpedance readings were recorded (46). The whole process was performed for different models with different number of surface electrodes (28, 32, 44 and 48) as well as different ground and terminal electrode dimensions (rterminal = 1, 2, 2.5 mm and rground = 0.5, 1 mm).
Finally, to find the optimized electrode configuration for the measurement set-up, the contribution of the different parts of the model, especially the venous segment, to the impedance measured between the surface and the ground electrode (total impedance) was computed.
Post processing Some parameters that cannot be measured in a simple way can be studied by post processing the simulation results. These include electrical potential, current paths and current density, sensitivity, volume impedance density, volumetric energy loss density, real and imaginary parts of the impedance and consequently phase angle and impedance magnitude. Below some of these parameters are briefly explained.
Current density Current density (J) is the electric current carried by a conductor per unit cross-sectional area of the conductor. Current density of every point throughout the measured area can be quantified by current density vectors (47). The current density vector is defined as a vector whose magnitude is the electric current per cross-sectional area at a given point in space, and its direction being that of the motion of the positive charges at this point (48). Current density is proportional to the electric potential as expressed by: where J is the current density vector, E is the electric potential and σ is the electrical conductivity.
Sensitivity and volume impedance density Sensitivity (S) is the ratio of the differential change in impedance to the differential change in conductivity within a voxel (49). The point sensitivity of impedance indicates how much the resistivity of the particular point influences the total impedance. Sensitivity is given by the following equation for a two-electrode system (50): where Jec is the current density due to current passed through the model and I is the magnitude of the current, assumed constant through both the terminal and ground electrodes (38). If the sensitivity (S) is multiplied with the complex resistivity (ρ) we get the volume impedance density (z) (38).
The volume impedance density can be integrated over all points in the measured site to get the measured impedance, Z (38): Volume impedance density can be employed to find the contribution from each constituent of the model to the total impedance measured in the model (38).
Volumetric energy loss density The complex admittance Y measured between the terminal and ground electrode (the inverse of the impedance) can be defined as: = + [5] where G is the conductance, i the imaginary unit, ω the angular frequency and C the capacitance (50). This applies if the measured material can be defined as a conductance, which is proportional to the dissipated part of the electric energy and a capacitance in parallel, which is proportional to the conserved part of the electric energy and with no induction (L = 0). This way, by integrating the total electric energy (conserved and dissipated) over different parts of the model, relative contribution of each part to the total capacitance and conductance and consequently the total calculated impedance, can be calculated (51). In this model, volumetric energy loss density was computed as another simple way to find the contribution of different components of the model to the total measured impedance. This was done by integrating the power loss density over the components of interest and comparing with the overall power loss (over the whole model). It was assumed that due to the small phase angle, the dissipated energy was higher and dominant and was therefore chosen over the conserved energy.
Impedance magnitude and phase angle Impedance is a complex parameter consisting of a real and an imaginary part. The real (Zreal) and imaginary (Zimag) components of the lumped parameter of impedance (Z) which are called resistance (R) and reactance (X) respectively, can be computed by using the built-in functions and post processing methods in COMSOL simulation (43).
Complex impedance can be described by polar coordinates (figure 3) (50). The impedance magnitude |Z| signifies the length of the impedance vector (3) and is defined as: The phase shift ϕ of the impedance describes the angle between the voltage and current (3) and can be defined as: = arctan ( ) [7] Figure 3: Diagram of the complex impedance illustrating the relationship between the real and imaginary parts of the impedance, where Z real is the real or resistive component R of the impedance, whereas Z imag is the imaginary or capacitive X component. |Z| and ϕ define the impedance magnitude and phase shift respectively. Both impedance magnitude and the phase angle can be computed by post processing methods in FEM simulation.

Ethical approval
The conducted research is not related to either human or animal use.  In figures 4 and 5, parameters computed by post processing the simulation results, such as electrical potential, current density, sensitivity, and volume impedance density are demonstrated by using slice plots. These plots demonstrate the absolute values and non-vector parameters and do not provide any direction information (38). The simulations were performed both with and without the vein segment as shown in figures 4 and 5. In addition, some of these postprocessed parameters such as volume impedance density and volumetric power loss density were then employed to compute the contribution of different parts of the model to the impedance measurements in the whole model.

Figures
The simulation of the electrical response of the model visualized by the planar color plots demonstrates the effects of electrode configuration and geometry on measurement depth and sensitivity.
The detailed simulation results for different models as a function of changing the electrode configuration (number and dimension) are shown in 18 tables that due to their excessive numbers are placed in the appendix. Hereafter we refer to radius of the terminal and ground electrode as rt and rg consecutively. However, two electrode dimensions of (rt = 1 mm and rg = 0.5 mm) and (rt = 1 mm and rg = 1 mm) provided the highest impedance contribution from the vein for all number of surface electrodes. Therefore, the summary of the simulation results with these two electrode dimensions are shown in tables 1 and 2 where all the surface electrodes are employed as terminals and excited simultaneously.  In table 3, the results of two different simulations are compared where both have employed only one surface electrode as the terminal but with different dimensions of (rt = 1 mm and rg = 0.5 mm) versus (rt = 2.5 mm and rg = 0.5 mm). The terminal electrode is employed each time in different positions (middle, first end closer to the ground electrode output, and second end) to see the influence of the location of the terminal electrode on the results and this is done regardless of the number of surface electrodes. Table 3: Simulations results including the percentage of the volumetric power loss density and volume integration of volume impedance density for vein and with two different surface and ground electrode dimensions (rt = 1 mm and rg = 0.5 mm versus rt = 2.5 mm and rg = 0.5 mm), regardless of the number of surface electrodes (28, 32, 44 and 48) and with only one surface electrode excited as the terminal (middle electrode, electrode at the first end close to the ground electrode output or electrode at the second end)

Discussion
The 3D AC/DC -Quasi-statics module of COMSOL Multiphysics FEM version 5.4 was employed for simulating the electrode configuration and the measurement set-up designed for monitoring the impedance of a venous segment. The build-in graphics tool was used to draw the geometry of each component of the measurement setup to which material and relevant dielectric properties were then allocated. The geometry was meshed with extra fine density. The simulation was then performed to analyze impedance changes in a two-electrode configuration system by varying the number of the surface electrodes (28, 32, 44 and 48) and their dimensions (rt = 1, 2, and 2.5 mm) as well as the dimension of the ground electrode (rg = 0.5 and 1 mm). In this way, a complete set of simulations consisting of many models with different electrode configuration was performed, where the surface electrodes where either excited one by one or all simultaneously as the terminal electrodes.

Terminal electrodes excited one by one
Our simulation results including the data in tables 1, 2 and 3 indicate that applying excitation to the surface electrodes one by one as the terminal provides higher contribution from the vein to the total measured impedance in comparison to when all the surface electrodes are employed as terminals and are excited at the same time.
Changing the dimensions of both the terminal and the ground electrodes influences impedance measurements. The results show that for the electrode dimensions of rt = 1 mm and rg = 0.5 and 1 mm, the impedance contribution of the vein to the total measured impedance is higher in comparison to that of other electrode dimensions. This is shown in table 3 where two different electrode dimensions (rt = 1 mm and rg = 0.5 mm) and (rt = 2.5 mm and rg = 0.5 mm) are compared while one surface electrode at a time is excited as the terminal electrode. This is regardless of the number of surface electrodes as increasing and decreasing the number of electrodes while keeping the dimension of the electrode constant, only moves the electrode slightly back and forth and the results have shown that the difference it makes in the measurements is so little that it can be neglected.
In addition, to investigate the influence of location for the surface electrode as the terminal, three different simulations were performed where one surface electrode at the first end closer to the ground electrode output, middle and the second end of the venous segment was chosen as the terminal electrode. The difference in contribution from the venous segments was found to be in the range of 0.1%, when comparing the different terminal electrode positions (Table 3). This shows that the terminal electrode position along the venous segment does not have a significant influence on the obtained sensitivity and that the results in figure 5c will be similar for all the terminal electrodes when they are excited one by one.
It can also be concluded from the results that the open space between the terminal electrodes is another important factor that should be chosen carefully to avoid blind spots along segment surface with too low sensitivity.
As for the ground electrode, according to our simulation results, increasing the radius of the ground electrode from 0.5 to 1 mm gave an increase to the contribution of the venous segment to some degree. However, it should be noted that the choices for the radius of the ground electrode would be limited considering the dimensions of the different venous segments that we are aiming for in our ex-vivo use of the method. The results suggest the choice of 1 mm for the radius of the ground electrode, however, the length of the ground electrode is more important. To obtain the same symmetry between the terminal and the ground electrode, for middle and end positions of the terminal electrode, we must have a ground electrode that is significantly longer than the whole length of the venous segment.
Exciting the surface electrodes as the terminal one at the time and by turn, requires a longer measurement time.
However, each excited surface electrode provides impedance information specifically from the site proximal to the electrode. In this way, by combining the measurement results and electrode position, a high spatial resolution would be obtained. The signal to noise ratio will also be much higher for the single terminal measurement compared to simultaneous measurement with all terminals in parallel. Hence, a higher sensitivity to local changes in the venous segment will be feasible. This is especially useful when it comes to monitoring and studying the venous segment with possible local deformation or damages in its tissue layers and cellular structure, as the higher sensitivity and spatial resolution would make it easier to detect and localize damage on the venous segment.
Terminal electrodes excited all simultaneously Results from tables 1 and 2 show that when all the terminal electrodes are excited at the same time, for the electrode dimensions of rt = 1 mm and rg = 0.5 and 1 mm, the number of 28 and 32 surface electrodes provide higher impedance contribution for the vein.
By assuming that the superposition principle is valid for our electrode simulations, we can illustrate the sensitivity distribution over the complete venous segment area by simultaneous excitation on all the terminal electrodes ( figure 4). With this method, we can locate possible blind spots between the electrodes where the sensitivity will be too low to detect local changes or damage in the venous segment. By investigating the sensitivity distribution as a function of terminal electrode size and position, we can optimize our electrode setup when it comes to size, number and geometrical distribution of our electrodes. Increasing the electrode spacing while keeping the same electrode dimensions causes a significant decrease in the maximum magnitude of the electric field, as the field becomes more spread (1). On the other hand, by increasing the number of electrodes we can obtain a more even current distribution and cover more of the area between the electrodes. In the simultaneous excitation setup, each electrode will contribute with an additional current path in parallel to all the other electrodes. All these current paths will contribute equally to the measured impedance, and the signal to noise ratio for local changes will decrease with the number of electrodes. Therefore, it can be expected that the sensitivity for each electrodes will be inversely proportional to the number of electrodes excited simultaneously. Thus, although applying excitation to all the surface electrodes at the same time decreases the measurement time, it also decreases the sensitivity and takes away the spatial resolution when it comes to detecting local structural changes or damages in the venous segment wall.

General findings
In our simulations, the total impedance is obtained by multiplying sensitivity with the resistivity in each point and integrate over all points in the volume of interest. In practical use for monitoring the properties of venous segment, the sensitivity for the application typically will be defined as the capability to discover a change in the segment. Our results demonstrate that whatever configuration for the terminal electrodes excitation is chosen, improving sensitivity for the vein monitoring application is possible by changing the dimensions of both the terminal and ground electrodes. However, the dimension of the terminal electrode is way more determinant. Results for the measured impedance, volumetric power loss density and volume impedance density, indicate that for all electrode numbers employed in our simulations (28, 32, 44 and 48), the surface electrode with the radius of 1 mm and ground electrode with the radius of 1 mm provide the highest percentage of contribution to the measured impedance from the venous segment. The reason to this might be that terminal electrodes with smaller dimensions concentrate the current density more in the relatively thin venous wall, which is the site of interest rather than the saline suspension.
Improving the spatial resolution of the measurements is very much dependent on the configuration chosen for terminal electrode excitation. Aiming for detection of local changes in a venous segment, the one-by-one terminal excitation will provide a higher spatial resolution in comparison to the simultaneous excitation of all terminal electrodes.
For both terminal electrode excitation configurations, the spacing between the terminal electrodes plays an important role in determining current distribution in the venous segment wall and consequently the area to be monitored. Therefore, while designing the measurement set-up, the positioning of the electrodes should be optimized to avoid areas or sites between the electrodes with insufficiently low current density.
The results also suggest that volume impedance density is more useful than sensitivity in estimating the contributions of different model parts to the measured impedance (47). Moreover, it can be inferred from the computations of the percentage of volume impedance density and volumetric power loss density that these two factors provide almost the same value when it comes to studying the contribution of each part of the model to the total impedance measurements.
Excitation of the terminal electrodes one by one as the optimal configuration for the experimental set-up can be achieved by various solutions such as applying a multiplexer which can employ all the surface electrodes as the terminal electrode one by one during the measurement process. This would make it possible to perform a thorough scanning of the whole vein segment for the purpose of bioimpedance spectroscopy measurements.

Conclusion
The finite element simulation method was used to evaluate the feasibility of using a two-electrode bioimpedance spectroscopy set-up to monitor a three-dimensional blood vessel model. Simulations showed that we, within the assumed tissue properties, could monitor the impedance in the model with sufficient focus on the venous segment. We also found that this simulation model would be a valuable tool to design and optimize an experimental set-up with sufficiently high sensitivity and spatial resolution. Hence, the proposed simulation model will reduce the time and expenses in our project by significantly reducing the use of ex-vivo animal models. Our project will be followed up with further optimization of the measurement set-up. Electrode configurations with an optimized combination of size, number and position of the terminal electrodes will be found by repeated simulations with emphasize on even current distribution and sensitivity focus in the vein segment.

Conflict of interest
Authors state no conflict of interest.