Representing unsegmented vessels using available vascular data for bioheat transfer simulation

Abstract

Introduction:

A primary challenge with voxel domains generated from imaging data is associated with voxel resolution. Due to the dimensional scale of blood vessels, not all vessels are captured in a given voxel resolution, leading to discontinuous blood vessels in the segmentation. Pre-capillary vessels like arterioles, which provide the highest resistance to blood flow, are often modeled with tissue as a porous domain due to resolution limitations. This results in a loss of information that could have been modeled if these vessels were segmented and modeled distinctly from the capillary bed.

Methods:

This paper focuses on developing mathematical equations to calculate the flow resistance of unsegmented vasculature with reference to flow resistance of available segmented vascular data. A 3D vascular domain of 32 terminal vessels and five generations of bifurcation is simulated. Each generation is successively removed and substituted with the new flow resistance equations to analyze the error in heat transfer due to a lack of segmentation data.

Results:

The effect of using mathematical equations of flow resistance on bioheat transfer is analyzed. Two methods are proposed and demonstrated to show considerable error reduction in bioheat transfer.

Discussion:

Very high image resolution, which could allow modeling of pre-capillary vessels, increases the computational cost of the entire simulation domain. Instead, a mathematical representation of the pressure drop induced in these unsegmented blood vessels is used. The proposed methods show potential in reducing the error resulting from the lack of segmentation data, improving the accuracy of bioheat transfer simulations.

Introduction

Technological advancements enable visualization and modeling of the vasculature (Ivilinov Todorov et al., 2020; Silvestri et al., 2021), providing highly detailed blood vessel domains. When coupled with accurate biophysics simulations (Amare et al., 2022; Blowers et al., 2018), such realistic models can be used to illustrate, understand, and predict biological responses to different environmental conditions. Such tools can predict patient response to medical treatment, changes in blood flow distribution due to burns or clots (Ng and Chua, 2002a; Ng and Chua, 2002b; Cookson et al., 2012), drug distribution, and damage to healthy tissue during hyperthermia treatments (Bellizzi et al., 2020; Silva et al., 2020).

However, a very high resolution data maybe required depending on the scale of domain size. The resolution of the voxel domain limits the visualization of blood vessels that can be modeled. Although capillary beds can be modeled using the porous media assumption in tissue (Amare et al., 2022; Hodneland et al., 2019), pre-capillary vessels that impact heat transfer (Weinbaum and Jiji, 1985; Lemons et al., 1987; Jiji et al., 1984) also play a crucial role in vasomotion. The arterioles, which can vary flow resistance, are present in pre-capillary vessels that often get overlooked in bioheat transfer. For an accurate and thorough simulation of biophysics, the ability to simulate the effects of arterioles is important; thus, modeling the pre-capillary network is crucial.

In a voxel domain generated from imaging data, a blood vessel ends where the resolution of the voxel can no longer identify it. One option to simulate the pre-capillary vessels from this point would be the mathematical modeling of blood vessels or an algorithm to simulate vascularization. Vascularization and angiogenesis are highly complex phenomena that include chemical, physical, and biological processes (Tran et al., 2022; Stefanini et al., 2012; Corada et al., 2014; Chappell et al., 2019; Tang et al., 2014). Vascular Endothelial Growth Factor (VEGF) signaling released by the tissue cells directs the tip cell to guide vascularization. Various biocomputational models have been developed to simulate this process (Takigawa-imamura et al., 2022; Kim et al., 2012; Heck et al., 2015; Zhang et al., 2022). However, computational models developed from a biological perspective (Zhang et al., 2022) are different from models developed from an engineering perspective. One of the primary approaches engineers use to model blood vessel growth is modeling blood vessels as fractals (Merks and Glazier, 2006; Tong and Fan, 2001; Murray, 2003a; Murray, 2003b; Lorthois and Cassot, 2010; Niemeyer et al., 1984) like the Diffusion Aggregation Method proposed by Fleury (Vincent and Schwartz, 1999; Vincent and Schwartz, 2000; Witten and Sander, 1981).

Another engineering method used to simulate vascularization is called Constrained Constructive Optimization (CCO) (Schreiner, 1993; Schreiner et al., 1995; Schreiner et al., 2006; Cury et al., 2021). In CCO, the main assumption is that blood flow is equally distributed in the specific organ domain. Based on this assumption, the supply blood flow rate in the organ is equally distributed in a given number of terminals. The blood flow rate passing through a single terminal vessel and the pressure drop between the supply node and the terminal of vasculature are provided as constraints. A random point is selected within the domain and a new branch is grown towards the point. The radius of the new branch is calculated such that the total volume of the vasculature is minimized. The CCO method has undergone various modifications such as parallelizing the growth of blood vessels (Talou et al., 2021; Cury et al., 2021; Shen et al., 2021) to increase computational efficiency. An example of the application of the CCO model is found in Correa-Alfonso’s work of vascularization on mesh liver model (Correa-Alfonso et al., 2022). The minimum diameter of a blood vessel in this liver model is 100 μm. The blood vessels in this model are shunted, i.e., arteries are directly connected with veins. The shunt between arteries and veins at 100 μm cannot model the time blood spends in the capillary bed. Blood flow is the slowest in the capillary bed and is expected to have a higher energy absorption rate due to the time spent there. The challenge with using the CCO method with a voxel phantom is associated with voxel resolution. For example, in the mesh liver model, to model the blood vessels of diameter 500 μm, a voxel resolution of less than 100 μm is required. When a cuboidal voxel of size 100 μm is used to voxelize the mesh, the liver model has a total of voxels; and about 70% of hepatic arteries are lost in the voxelization of the liver mesh as they have a diameter of less than 500 μm.

Modeling every blood vessel to the capillary bed is not a feasible solution. Thus, a new method is required to model the flow resistance and heat transfer of blood vessels that exist between the capillary bed and the segmented blood vessels. In previous work (Amare et al., 2022; Hodneland et al., 2019) a multiscale, mixed-dimensional simulation framework was developed for simulating multiphysics, known as the VoM-PhyS framework. This framework uses a pressure drop parameter (Hodneland et al., 2016; Hodneland et al., 2019) to simulate the resistance of the “virtual” blood vessels that are not segmented. However, no in-vivo or empirical data are available to determine the pressure drop parameter. In previous work (Amare et al., 2023) a 2D domain was studied and the major findings were: (1) The flow resistance of unsegmented vessels can be recovered using the pressure drop parameter in the VoM-PhyS framework; (2) A constrained range was observed within which the correct pressure drop parameter value exists for a given simulation domain; and (3) The error in flow and pressure maps due to lack of segmented vessels can be reduced with the use of correct pressure drop parameter. This paper focuses on expanding on these findings for a 3D domain with the following aims: (1) Provide equations to calculate the correct pressure drop parameter for any available vascular domain; (2) Study the effect of unsegmented vessels on bioheat transfer in the VoM-PhyS framework; and (3) Provide rectification methods to reduce the temperature error.

This study aims to address three key challenges in bioheat transfer modeling: (1) Develop mathematical equations to calculate the pressure drop parameter for unsegmented vasculature using only net flow rate and flow resistance of segmented vessels. (2) Quantify the effect of unsegmented vessels on bioheat transfer simulations. (3) Propose and evaluate methods to reduce temperature errors in bioheat simulations with limited vascular data. The novelty of this work lies in: (a) The development of a mathematical framework to represent unsegmented vessels, bridging the gap between high-resolution vascular modeling and practical imaging limitations. (b) Demonstration of the impact of unsegmented vessels on bioheat transfer simulations, highlighting that correcting flow resistance alone is insufficient for accurate heat transfer modeling. (c) Proposal of two methods to reduce temperature errors: implementing an effective thermal conductivity approach and assuming equal spatial distribution from terminal vessels.

Methods

For this research, the VoM-PhyS framework was used to simulate blood flow coupled with heat transfer. The pressure drop parameter, , is an important parameter in the VoM-PhyS framework that permits modeling the unsegmented pre-capillary vessels. The ratio of blood viscosity and pressure drop parameter is an empirical representation of flow resistance offered by unsegmented blood vessels. A detailed description of this parameter can be found in (Hodneland et al., 2019; Amare et al., 2022). The goal of this work is to provide mathematical equations to calculate . To derive the equations for , a blood vessel domain generated using the CCO algorithm with the assumption of equal flow distribution in the domain was used. This domain is shown in Figure 1. In the VoM-PhyS framework, the porous tissue domain is modeled as two-compartments, arterial and venous, and coupled using the perfusion parameter, . The arterial pressure drop parameter, , controls the flow resistance of unsegmented blood vessels between arterial terminal nodes 3 and arterial compartment pressure nodes. Similarly, the venous pressure drop parameter, , controls the flow resistance between venous compartment pressure nodes and venous terminal pressure nodes 4.

FIGURE 1

Step 1: mass conservation and total flow

The blood flow rate flowing in the volume is from the arterial side, and the blood flow rate leaving the domain is , respectively. Since mass is conserved, is equal to as shown in Equation 1. In Figure 1, the arteries (red) and veins (blue) are considered to be the segmented vasculature for which the length and diameter can be calculated from imaging data.

The blood entering the tissue domain is considered to be uniformly distributed (Schreiner, 1993; Schreiner et al., 1995; Xing et al., 2022; Correa-Alfonso et al., 2022). Thus, the flow rates leaving the arterial terminal nodes are considered to be equal. With the equal distribution assumption, the blood flow rates in the terminal vessels are calculated using Equation 2, where represents the number of bifurcation generations. In the illustrated example, is equal to 2.

Step 2: define equivalent resistances and pressure drops

Using the dimensions of the vasculature, the equivalent flow resistance offered by each arterial and venous vascular tree can be calculated using parallel and series resistance methods. The equivalent resistance of arterial and venous trees are represented as and , respectively, and shown in Figure 2. Figure 2 represents the simplification of the domain from Figure 1 by calculating equivalent flow resistances of segmented arteries and veins.

FIGURE 2

Using the simplified blood vessel model, the blood flow equations in the arterial and venous trees can be written as Equation 3 and Equation 4, respectively. represents the pressure at node 3 where the segmented arteries end, and represents the pressure at node 4 where the segmented veins end.

The total pressure drop between terminal arteries and terminal veins represents the pressure drop across unsegmented vasculature. This pressure drop can be calculated using Equation 5, where and represent the pressure drop between the terminal artery and arterial compartment and venous compartment and terminal vein, respectively. represents the pressure drop across arterial and venous compartments within the tissue and is calculated using the perfusion coefficient, , and total volume, , of the tissue domain as shown in Equation 6.

Step 3: correlation between resistances and pressure drop

A correlation between the flow resistance of arterial and venous trees and the pressure drop between terminal vessels and tissue is proposed in Equation 7. The correlation states that the ratio of pressure drop between unsegmented arteries to capillary bed and unsegmented veins to capillary bed is equal to the ratio of overall flow resistance offered by segmented arteries and veins. The constand is used for derivation only and it represents the ratio of effective arterial resistance to effective venous resistance.

Using Equation 8 and substituting in Equation 5, equations for and are derived, shown in Equation 9 and Equation 10, respectively.

Step 4: solve for pressure drop parameters

The resistance equation for flow in virtual vessels can be written as Equations 11, 12 where the resistance offered by virtual arterial and venous networks is represented using blood viscosity and pressure drop parameter .

Using Equations 9–12, the equations to determine pressure drop parameter for the arterial and venous tree is derived as shown in Equation 13 and Equation 14, respectively.

Equation 13 and Equation 14 are the final forms of equations to calculate the pressure drop parameter for arterial and venous trees. It can be seen from these equations that the pressure drop parameters can be calculated using total blood flow rate in the simulation domain, pressure drop across the simulation domain, and the equivalent resistance of arterial and venous trees calculated from the segmented vasculature. The demonstration of how to use the pressure drop equations for a given biological domain is shown in Figure 3.

FIGURE 3

3D domain

To verify the applicability of the Equation 13 and Equation 14 on a 3D domain, a test domain shown in Figure 4 was generated using Rhinoceros (McNeel et al., 2010). The domain has 32 terminals for arterial and venous trees and the cuboidal tissue size encasing the vasculature was × × . The voxel dimensions are × × . The 32 terminal domain was considered as a reference solution (most segmented data available) for comparison and will be referred to as Case 1. Similar to previous work (Amare et al., 2023) on a 2D domain, the number of terminals was reduced gradually to create different simulation domains representing lack of segmented data as shown in Figure 4. The resultant domains are termed Case 2, Case 3, Case 4, and Case 5 for 16, 8, 4, and 2 terminals, respectively (Figure 4). The dimensions of the arterial and venous trees were considered same for this simulation for simplicity and are given in Table 1.

FIGURE 4

Since this was a 3D modeled domain, a reference volumetric blood flow rate had to be determined. The VoM-PhyS framework (Amare et al., 2022; Hodneland et al., 2019) was used with the flow parameters given in Table 2 to calculate the overall volumetric blood flow rate in Case 1. The pressure drop parameter for this simulation was considered as . This ensured that only the flow resistance of segmented vessels in Case 1 determine the overall blood flow rate.

The parameters used in this 3D blood flow and heat transfer simulation model are based on established physiological values and commonly accepted approximations in biomedical engineering. The blood viscosity and thermal conductivity are consistent with values reported in the IT’IS database for biological tissues (Hasgall et al., 2018). The arterial and venous permeability values (, ) are supported by previous computational models of tissue perfusion (Amare et al., 2022; Hodneland et al., 2016; Hodneland et al., 2019). The perfusion rate falls within the typical range for various tissue types (Jeong et al., 2023). The specific heat and density of the tissue are approximated to those of water, a common practice in biological heat transfer models (Duck, 2013). The ambient temperature represents a standard room temperature, while the inlet blood temperature is slightly below core body temperature, accounting for cooler peripheral blood. The metabolic heat generation rate is within the range observed in various tissues, albeit on the higher end (Wahyudi et al., 2022). The convective heat transfer coefficients for ambient and blood (, ) are typical values used in biological heat transfer simulations. While the inlet and outlet pressures are lower than physiological arterial pressures, they are representative for the specific modeling conditions and simulate a pressure drop of arround (Timothy, 2016; Li et al., 2012).

Using Equations 13, 14, the pressure drop parameters for all five cases were calculated. The blood flow was simulated with the respective pressure drop parameters and the pressure maps were compared with the reference Case 1.

Pressure drop parameter equations

The pressure drop parameters calculated using Equation 13 and Equation 14 for the 3D reference domain (Figure 4) are shown in Table 3. Nt represents the number of vascular terminals for respective cases. To understand the importance of pressure drop parameter and its effect on pressure solution, Figure 5 is shown. Here the percentage error in pressure at arterial nodes , , and of Case 2 is plotted for 1%, 2%, and 5% change in . Nodes 5 represent the last segemnted node in Case 2, Nodes 4 represent the nodes connected to Nodes 5, and Nodes 3 represent the pressure nodes connected to Nodes 4. There are many biological factors like non-newtonian behavior of blood and vasomotion which are not considered in Equations 13, 14. Thus, these equations can only provide an accurate flow resistance behavior of unsegmented vessels. A 0.6% change in pressure for 5% in pressudre drop parameter can be seen from Figure 5. This shows that an accurate pressure drop parameter will be sufficient for preliminary analysis and simulation of blood flow.

FIGURE 5

The flow equations were solved using the pressure drop parameters of respective cases and pramaters shown in Table 2. The resultant pressure values are given in Tables 4–6, for the arterial tree, venous tree, and tissue compartments, respectively. In Table 6, the average, maximum, and minimum pressures in the arterial and venous compartment of all the cases is shown. A contour map for pressure error between Cases 2 to 5 and Case 1, is shown in Figure 6, at z = 80 the location of z in the domain shown in Figure 7. The pressure difference between Case 5 and Case 1 at each tissue voxel is within . Case 5 represents the worst case possible with the least segmented vasculature available.

FIGURE 6

FIGURE 7

Bioheat transfer

To analyse the effect of lack of segmentation data on bioheat transfer, Case 1 temperature profile was considered as reference solution, and all other Cases were compared to it. The temperature errors of Case 2, Case 3, Case 4, and Case 5 at z = are shown in Figure 8. The location of z = 80 was selected as it had the maximum temperature error in the entire domain. The dimensionless temperature error was calculated using Equation 15. represents the ambient temperature used for simulation, represents the inlet blood temperature, represents the temperature of ith voxel in Case 1, and represents the temperature of ith voxel in Case c where c.

FIGURE 8

Discussion

Tables 4, 5 show that with the use of the correct value for the pressure drop parameter, the flow resistance of the unsegmented vessels can be simulated and the pressure drop in the vascular tree can be accurately determined. For Case 5, which represents the worst-case scenario with only two terminal arteries and veins, the pressure contours in the arterial and venous compartments of tissue show a maximum pressure error of . The corresponding error is less than 0.5% when compared to the arterial and venous compartment pressure values in reference Case 1.

The greatest flow resistance for a blood is at arterioles. These pre-capillary arterioles are not always captured in the segmented data. Any change in the flow resistance of these pre-capillaries would result in changes in total flow rate and pressure drop. In the example under consideration, since the pressure drop across the domain is given as boundary condition, the pressure drop parameter affects the total blood flow rate, as shown in the previous work (Amare et al., 2023). The pressure drop parameter thus controls an important aspect in the VoM-PhyS framework. A small change in pressure drop parameter changes the pressure drop across the vascular tree. A local vasomotion can be simulated by varying the pressure drop parameter of respective vessel terminal.

The pressure difference between Case 5 and Case 1 at each tissue voxel is within and negligible, the temperature difference between these two cases is within 2.5°C as shown in Figure 8. This temperature difference is substantial and cannot be ignored. Thus, correcting flow resistance alone does not guarantee reduction in bioheat simulation error. Blood vessels are considered to affect heat transfer via countercurrent flow (Weinbaum et al., 1984). When vessels are unsegmented, their effect on surrounding tissue cannot be simulated. To minimize this temperature error, two methods were considered: effective thermal conductivity and larger Sphere of Influence (SoI) radius (Amare et al., 2022; Amare et al., 2023).

Effective thermal conductivity

Between two tissue voxels, heat transfer takes place via advection and conduction. As the blood flows from one porous tissue voxel to another, it carries the heat with it resulting in advection. As the pressure map was consistent with maximum error within , there was negligible change in the flow rate among tissue voxels. However, as fewer blood vessels are segmented, they are modeled as part of tissue porosity. The porous tissue voxels now consist of larger blood vessel than the reference domain (Case 1). Though there was no change in heat transfer due to advection, thermal conduction among tissue voxels can change due to larger unsegmented vessels. In literature (Chen and Holmes, 1980; Keller and Seiler, 1971; Weinbaum et al., 1984; Weinbaum and Jiji, 1985; Roetzel and Xuan, 1997; Nakayama and Kuwahara, 2008), effective thermal conductivity is used to consider the effect of unsegmented blood vessels and counter-current heat exchange. The same concept of effective thermal conductivity was used to compensate for the effect of cross-flow heat exchange between unsegmented blood vessels. The thermal conductivity of the tissue voxels was varied between to . The resultant temperature errors between Case 5 and Case 1 for different values of tissue thermal conductivity are shown in Table 7 with the temperature contours plots at z = 80 shown in Figure 9. From Figure 9 and Table 7 it can be seen that as the tissue thermal conductivity increases, the temperature error begins to reduce till it reaches a threshold value, beyond which the temperature error seems to increase.

FIGURE 9

Larger SoI

Sphere of Influence (SoI) is a parameter introduced in previous work (Hodneland et al., 2016; Hodneland et al., 2019) and is critical in the VoM-PhyS framework (Amare et al., 2022) for coupling 1D flow with 3D flow. The SoI is a volume of sphere with origin center at a terminal vessel. The tissue voxels that fall within a given SoI are considered to exchange blood with the respective terminal vessel. The radius, , of the SoI is an empirical quantity. The effect of on blood flow and pressure drop was studied (Amare et al., 2023) and did not demonstrate any effect on pressure drop. However, has been shown to affect thermal maps in tissue (Amare et al., 2022). Thus, different values of were studied to understand its effect on reduction of temperature error due to unsegmented vessels. The temperature error values for different and different values of tissue thermal conductivity, , are given in Tables 8, 9. The temperature difference contours at z = 80 for Case 5 with as is shown for different values of . From Table 8 and Figure 10, it can be seen that as the increases, the temperature error begins to decrease and reaches a threshold value beyond which there is no change in the temperature error. This behavior is different than the effect of , where beyond the threshold, the temperature error begins to increase again.

FIGURE 10

To understand the effect of these two methods on temperature error reduction, root mean square error (RMSE) and the sum of temperature error (STE) in the entire domain were calculated. The values of RMSE and STE are shown in Tables 8, 9, and a graphical plot of RMSE and for the different values of and is shown in Figures 11, 12, respectively. In Figure 11, it can be seen that as increases, the RMSE decreases for to a threshold value. Beyond the threshold, the RMSE increases with an increase in . Using polynomial regression, six polynomial functions of the fourth order were fitted for RMSE for each value of , respectively. Similarly, six polynomial functions of fourth order were fitted for STE for each value of . Using the RMSE function, the threshold where minimum RMSE occurs was calculated, along with the corresponding STE. These values are given in Table 10. It is noteworthy that the STE is positive for all at the threshold of minimum RMSE. This denotes that the VoM-PhyS framework would result in higher temperatures than the reference for a less segmented vascular domain. Similarly, the where STE equals zero was calculated along with the corresponding RMSE. These values are given in Table 11. The values of where STE is zero are not the same where minimum RMSE occurs, as can be seen. A further detailed statistical analysis of these values could lead to greater insight into the performance of the proposed energy error reduction methods. A further analysis is needed to determine what parameters affect the value of for any given vascular data. Another noteworthy observation is that the minimum temperature for larger SoI is for , and this error increases for the same SoI when tissue thermal conductivity is varied. A SoI larger than the minimum required for 100% coverage ensures the tissue domain lies closer to the source than the periphery of the SoI. This provides more blood flow to the entire tissue domain than the simulation case when SoI is restricted to minimum for 100% coverage. As the SoI is increased, it is expected to achieve equal distribution in the entire tissue domain, and is expected when the SoI is considerably larger than the domain dimensions. This behavior ensures that we reach a plateau beyond which the temperature error cannot be decreased even if the SoI is increased.

FIGURE 11

FIGURE 12

The use of the equations to calculate correct pressure drop parameter for a given vasculature opens avenues for furthering this research field. Availability of segmentable vascular model remains a challenge, but the use of the pressure drop parameter equation provides a novel way to overcome this challenge and obtain a detailed accurate blood flow simulation. This study clearly shows the effect of unsegmented vasculature is dominant on heat transfer simulation. Further research is required to better understand methods to quantify this error and be able to provide rectification methods from available data. Similar to the equations proposed in this paper for calculating pressure drop parameter, if equations to rectify the temperature error could be derived, that would provide a major revolution in the field. Such equations could help simulate large scale domains which otherwise could not be easily modeled due to their computational memory requirement.

These findings also highlight that correcting flow resistance alone is insufficient for accurate heat transfer modeling. This insight is critical for researchers developing bioheat transfer models, as it underscores the need for additional considerations beyond flow resistance.

One of the major limitation of this work is tied to the SoI. The SoI is an empirical value and no data is available to determine its accuracy. The SoI radius is expected to vary based on the biological domain under consideration, tissue properties, vasomotion, and various other biological features. A detail study is required to better understand this parameter.

To derive the equations of pressure drop parameter, a simulation domain where segmented vasculature could be gradually removed was required. Obtaining such a simulation domain from medical imaging data still remains a challenge as discussed in the Introduction. The resolution of image required to get the desired segmentable vasculature is very high, increasing the computational overhead for simulation. The pressure drop parameter itself is the solution to model flow resistance of unsegmented blood vessels and hence an artificial 3D domain with blood vessels was developed for derivation of mathematical equations. Future work will include demonstrating the use of these equations on biological organs obtained from medical imaging scans.

While our proposed methods significantly improve bioheat transfer simulations with limited vascular data, several limitations should be noted. The pressure drop parameter equations assume steady-state flow and do not account for pulsatile effects or non-Newtonian blood behavior. The effectiveness of the error reduction methods may vary depending on the specific organ or tissue being modeled. Additionally, validation against in vivo measurements remains challenging due to the complexity of obtaining high-resolution temperature data in living tissues.

The broader implications of this work extend to various fields. In medical applications, more accurate bioheat transfer models could improve the planning and execution of thermal therapies, such as hyperthermia treatments for cancer. Our current research focuses on this and will be published in future work. In physiological research, these methods could enhance our understanding of thermoregulation in different organs. For computational biology, our approach offers a pathway to simulate large-scale domains that were previously computationally prohibitive, potentially enabling more comprehensive whole-organ or even full-body simulations.

These findings provide novel equations to calculate the pressure drop parameter for unsegmented vasculature using only the net flow rate and flow resistance of segmented vessels. This advancement allows for accurate pressure mapping in simulations with limited vascular network data, addressing a significant challenge in the field. The study provides a detailed analysis of how the lack of segmented vascular data affects temperature profiles in bioheat transfer simulations and proposes two approaches to mitigate these errors: implementing an effective thermal conductivity approach and assuming equal spatial distribution from terminal vessels in the tissue domain. These methods show significant improvements in simulation accuracy, with the effective thermal conductivity approach reducing maximum absolute temperature errors from to in the test case.

To contextualize our temperature error results, we compared them with literature benchmarks. Weinbaum and Jiji (Weinbaum and Jiji, 1985) reported temperature variations of up to - due to vascular effects in their bioheat transfer model. Our maximum absolute temperature error of in the worst-case scenario (Case 5) is consistent with this range. However, our proposed error reduction methods significantly improve upon this, with the effective thermal conductivity approach reducing the maximum error to . This improvement is particularly significant in the context of hyperthermia treatments, where temperature accuracy within 1- is crucial for treatment efficacy and safety (Kok et al., 2015).

In summary, this research contributes to the ongoing effort to develop more accurate and computationally efficient bioheat transfer models. It bridges the gap between high-resolution vascular modeling and practical limitations in medical imaging, offering a more accessible yet accurate method for simulating physiological processes. The novel approaches presented in this study have the potential to significantly impact various fields, from medical research to thermal regulation studies, paving the way for more sophisticated and realistic modeling of heat transfer in living tissues.

Conclusion

The pressure drop parameter equations derived in this work demonstrate a robust method for recovering the flow resistance of unsegmented vasculature using only the net flow rate in the simulation domain and the flow resistance of segmented vessels. This approach ensures accurate pressure mapping in simulations, even with limited vascular network data, significantly advancing our ability to model complex physiological systems. While the pressure and flow distribution can be simulated with high accuracy using these equations, the absence of segmented vessels introduces notable errors in temperature profiles during bioheat transfer simulations. This finding underscores the intricate relationship between vascular structure and heat transfer in biological tissues. To address these temperature discrepancies, two effective methods were identified: (1) Implementing an effective thermal conductivity approach (2) Assuming equal spatial distribution from terminal vessels in the tissue domain.

These methods substantially reduce simulation errors, with the optimal approach depending on specific domain characteristics. The effective thermal conductivity method shows particular promise, reducing maximum absolute temperature errors from to in our test case. The combined use of accurate pressure drop parameters and error reduction techniques represents a significant advancement in bioheat transfer modeling. This approach bridges the gap between high-resolution vascular modeling and practical limitations in medical imaging, offering a more accessible yet accurate method for simulating physiological processes. These findings open new avenues for research in computational biology, potentially enabling more comprehensive whole-organ or even full-body simulations that were previously computationally prohibitive. The methods developed in this study offer a practical solution for improving computational simulations with low-resolution data, which is particularly valuable given the ongoing challenges in obtaining high-resolution vascular imaging data.

Future work should focus on validating these methods against experimental data and exploring their applicability in diverse anatomical structures and pathological conditions. Additionally, incorporating non-Newtonian blood behavior and dynamic vasomotion effects could further enhance the physiological relevance of these simulations. In conclusion, this research contributes to the ongoing effort to develop more accurate and computationally efficient bioheat transfer models, with potential applications ranging from improving medical treatments to advancing our understanding of thermoregulation in living organisms.

Author summary

Our paper address the challenge of limted resolution in voxel domains derived from imaging data, particularly in capturing small blood vessels. We propose a mathematical representation of pressure drop in these unsegmented vessels within tisssue, reducing the need for high-resolution imaging. Since these equations can only approximate the true resistance, error is expected when compared to a detailed model. The two methods proposed to reduce the error in heat transfer show promising results. These equations can be further modified to simulate non-newtonian behavior of blood and could provide a practical solution for improving computational simulation with low-resolution data.

References

• 1

  AmareR.BahadoriA. A.EckelsS. J. (2023). “Analysis of sphere of influence (SOI) and pressure drop parameter in VoM-PhyS framework,” in Proceeding of 8th thermal and fluids engineering conference (TFEC) (Connecticut: Begellhouse), 81–88. Available at: https://dl.astfe.org/conferences/tfec2023,4dee63bd35f0bf74,34196bd44728bf18.html.10.1615/TFEC2023app045913

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AmareR.HodnelandE.RobertsJ. A.BahadoriA. A.EckelsS. (2022). Modeling a 3-D multiscale blood-flow and heat-transfer framework for realistic vascular systems. Sci. Rep.12 (1), 14610–14618. 10.1038/s41598-022-18831-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BellizziG. G.SumserK.VilasBoas-RibeiroI.CurtoS.DrizdalT.van RhoonG. C.et al (2020). Standardization of patient modeling in hyperthermia simulation studies: introducing the Erasmus Virtual Patient Repository. Int. J. Hyperth.37 (1), 608–616. 10.1080/02656736.2020.1772996

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BlowersS.MarshallI.ThrippletonM.AndrewsP.HarrisB.BethuneI.et al (2018). How does blood regulate cerebral temperatures during hypothermia?Sci. Rep.8 (1), 7877. 10.1038/s41598-018-26063-7

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  ChappellJ. C.DardenJ.PayneL. B.FinkK.BautchV. L. (2019). Blood vessel patterning on retinal astrocytes requires endothelial flt-1 (VEGFR-1). J. Dev. Biol.7 (3), 18–13. 10.3390/jdb7030018

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  ChenM. M.HolmesK. R. (1980). Microvascular contributions in tissue heat transfer. Ann. N. Y. Acad. Sci.335.1, 137–150. 10.1111/j.1749-6632.1980.tb50742.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  CooksonA. N.LeeJ.MichlerC.ChabiniokR.HydeE.NordslettenD.et al (2012). A novel porous mechanical framework for modelling the interaction between coronary perfusion and myocardial mechanics. J. Biomechanics45 (5), 850–855. 10.1016/j.jbiomech.2011.11.026

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  CoradaM.MoriniM. F.DejanaE. (2014). Signaling pathways in the specification of arteries and veins. Arteriosclerosis, Thrombosis, Vasc. Biol.34 (11), 2372–2377. 10.1161/ATVBAHA.114.303218

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  Correa-AlfonsoC. M.WithrowJ. D.DomalS. J.XingS.ShinJ.GrassbergerC.et al (2022). A mesh-based model of liver vasculature: implications for improved radiation dosimetry to liver parenchyma for radiopharmaceuticals”. In: EJNMMI Phys., 9.128. 10.1186/s40658-022-00456-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  CuryL. F. M.Maso TalouG. D.Younes-IbrahimM.BlancoP. J. (2021). Parallel generation of extensive vascular networks with application to an archetypal human kidney model. R. Soc. Open Sci.8, 12. 10.1098/rsos.210973

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  DuckF. (2013). Physical properties of tissues: a comprehensive reference book. Academic Press.

  • Google Scholar
  • View reference in article

• 12

  HasgallP. A.Di GennaroF.BaumgartnerC.NeufeldE.LloydB.GosselinM. C.et al (2018). IT’IS Database for thermal and electromagnetic parameters of biological tissues. Available at: http://www.itis.ethz.ch/virtual-population/tissue-properties/overview%0Ahttp://www.itis.ethz.ch/virtual-population/tissue-properties/overview.10.13099/VIP21000-04-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  HeckT. A. M.VaeyensM. M.Van OosterwyckH. (2015). Computational models of sprouting angiogenesis and cell migration: towards multiscale mechanochemical models of angiogenesis. Math. Model. Nat. Phenom.10 (1), 108–141. 10.1051/mmnp/201510106

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  HodnelandE.HansonE.Munthe-KaasA. Z.LundervoldA.NordbottenJ. M. (2016). Physical models for simulation and reconstruction of human tissue deformation fields in dynamic MRI. IEEE Trans. Biomed. Eng.63 (10), 2200–2210. 10.1109/TBME.2015.2514262

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  HodnelandE.HansonE.SævareidO.NævdalG.LundervoldA.ŠoltészováV.et al (2019). A new framework for assessing subject-specific whole brain circulation and perfusion using mri-based measurements and a multiscale continuous flow model. PLoS Comput. Biol.15 (6), 10070733. 10.1371/journal.pcbi.1007073

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  Ivilinov TodorovM.PaetzoldJ. C.SchoppeO.TettehG.ShitS.EfremovV.et al (2020). Machine learning analysis of whole mouse brain vasculature. Nat. Methods17 (4), 442–449. 10.1038/s41592-020-0792-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  JeongJ.-H.LeeB.HongJ.YangT. H.ParkY. H. (2023). Reproduction of human blood pressure waveform using physiology-based cardiovascular simulator. Sci. Rep.13 (1), 7856. 10.1038/s41598-023-35055-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  JijiL. M. M.WeinbaumS.LemonsD. E. E. (1984). Theory and experiment for the effect of vascular microstructure on surface tissue heat transfer—Part II: model formulation and solution. J. Biomechanical Eng.106 (4), 331–341. 10.1115/1.3138502

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  KellerK. H.SeilerL. (1971). An analysis of peripheral heat transfer in man. J. Appl. physiology30 (5), 779–786. 10.1152/jappl.1971.30.5.779

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  KimE.StamatelosS.CebullaJ.BhujwallaZ. M.PopelA. S.PathakA. P. (2012). Multiscale imaging and computational modeling of blood flow in the tumor vasculature. Ann. Biomed. Eng.40.11, 2425–2441. 10.1007/s10439-012-0585-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  KokH. P.WustP.StaufferP.BardatiF.van RhoonG.CrezeeJ. (2015). Current state of the art of regional hyperthermia treatment planning: a review. Radiat. Oncol.10, 1–14. 10.1186/s13014-015-0503-8

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  LemonsD. E.ChienS.CrawshawL. I.WeinbaumS.JijiL. M. (1987). Significance of vessel size and type in vascular heat transfer. Am. J. Physiology-Regulatory, Integr. Comp. Physiology253 (1), R128–R135. 10.1152/ajpregu.1987.253.1.r128

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  LiY.ShiZ.CaiY.FengY.MaG.ShenC.et al (2012). Impact of coronary tortuosity on coronary pressure: numerical simulation study. PLoS One7 (8), e42558. 10.1371/journal.pone.0042558

  • CrossRef
  • Google Scholar

• 24

  LorthoisS.CassotF. (2010). Fractal analysis of vascular networks: insights from morphogenesis. J. Theor. Biol.262 (4), 614–633. 10.1016/j.jtbi.2009.10.037

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  McNeelR.et al (2010). Rhinoceros 3D, version 6.0. Seattle, WA, USA: Robert McNeel and Associates.

  • Google Scholar
  • View reference in article

• 26

  MerksR. M. H.GlazierJ. A. (2006). “Dynamic mechanisms of blood vessel growth,” in Nonlinearity. 10.1088/0951-7715/19/1/000

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  MurrayJ. D. (2003a). Math. Biol. II. Spatial Models Biomed. Appl.18 (4), 476. 10.1016/jcamwa2013

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  MurrayJ. D. (2003b). On the mechanochemical theory of biological pattern formation with application to vasculogenesis. Com. Ren. Bio326, 239, 252. 10.1016/S16310691(03)00065-9

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  NakayamaA.KuwaharaF. (2008). “A general bioheat transfer model based on the theory of porous media”. In: Int. J. Heat Mass Transf.51.3190–3199 pp. 10.1016/j.ijheatmasstransfer.2007.05.030

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  NgE. Y. K.ChuaL. T. (2002a). Prediction of skin burn injury. Part 1: numerical modelling; Part 2: parametric and sensitivity analysis. Proc. Institution Mech. Eng. Part H J. Eng. Med.216 (6), 426–427. 10.1243/095441102321032229

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  NgE. Y. K.ChuaL. T. (2002b). Prediction of skin burn injury. Part 2: parametric and sensitivity analysis. Proc. Institution Mech. Eng. Part H J. Eng. Med.216 (3), 171–183. 10.1243/0954411021536388

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  NiemeyerL.PietroneroL.WiesmannH. J. (1984). “Fractal dimension of dielectric breakdown”. In: Phys. Rev. Lett.52.12 pp. 1033–1036. 10.1103/PhysRevLett.52.1033

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  RoetzelW.XuanY. (1997). Bioheat equation of the human thermal system. Tech. Rep.20, 268–276. 10.1002/ceat.270200407

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  SchreinerW.KarchR.NeumannM.NeumannF.SzawlowskiP.RoedlerS. (2006). Optimized arterial trees supplying hollow organs. Med. Eng. Phys.28.5, 416–429. 10.1016/j.medengphy.2005.07.019

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  SchreinerW. (1993). Computer generation of complex arterial tree models. J. Biomed. Eng.15 (2), 148–150. Available at: https://linkinghub.elsevier.com/retrieve/pii/0141542593900462.10.1016/0141-5425(93)90046-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  SchreinerW.NeumannF.NeumannM.EndA.RoedlerS. M.AharinejadS. (1995). The influence of optimization target selection on the structure of arterial tree models generated by constrained constructive optimization. J. General Physiology106 (4), 583–599. 10.1085/jgp.106.4.583

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  ShenJ.FaruqiA. H.JiangY.MaftoonN. (2021). Mathematical reconstruction of patient-specific vascular networks based on clinical images and global optimization. IEEE Access9, 20648–20661. 10.1109/ACCESS20213052501

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  SilvaM.FreitasB.AndradeR.Espregueira-MendesJ.SilvaF.CarvalhoÓ.et al (2020). Computational modelling of the bioheat transfer process in human skin subjected to direct heating and/or cooling sources: a systematic review. Ann. Biomed. Eng.48 (6), 1616–1639. 10.1007/s10439-020-02515-y

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  SilvestriL.MüllenbroichM. C.CostantiniI.Di GiovannaA. P.MazzamutoG.FranceschiniA.et al (2021). Universal autofocus for quantitative volumetric microscopy of whole mouse brains. Nat. Methods18.8, 953–958. 10.1038/s41592-021-01208-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  StefaniniM. O.QutubA. A.GabhannF. M.PopelA. S. (2012). Computational models of VEGF-associated angiogenic processes in cancer. Math. Med. Biol.29.1, 85–94. 10.1093/imammb/dqq025

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  Takigawa-imamuraH.HiranoS.WatanabeC.Ohtaka-MaruyamaC.EmaM.MizutaniK. (2022).Computational model exploring characteristic pattern regulation in periventricular vessels”. 12, 2069. 10.3390/life12122069

  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  TalouG. D. M.SafaeiS.HunterP. J.BlancoP. J. (2021). Adaptive constrained constructive optimisation for complex vascularisation processes. Sci. Rep.11 (1), 6180. 10.1038/s41598-021-85434-9

  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  TangL.van de VenA. L.GuoD.AndasariV.CristiniV.LiK. C.et al (2014). Computational modeling of 3D tumor growth and angiogenesis for chemotherapy evaluation. PLoS ONE9 (1), e83962. 10.1371/journal.pone.0083962

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  TimothyW. (2016). Hemodynamics. Compr. Physiol.6 (2), 975–1003. 10.1002/cphy.c150038

  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  TongS.FanY. (2001). “Numerical simulations of angiogenesis in the cornea,” in Microvascular research61 (1), 14–27. 10.1006/mvre.2000.2282

  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  TranK. A.Baldwin‐LeclairA.DeOreB. J.AntisellM.GalieP. A. (2022). Oxygen gradients dictate angiogenesis but not barriergenesis in a 3D brain microvascular model. J. Cell. Physiology237 (10), 3872–3882. 10.1002/jcp.30840

  • CrossRef
  • Google Scholar

• 47

  VincentF.SchwartzL. (1999). Diffusion limited aggregation from shear stress as a simple model of vasculogenesis. Fractals7.1, 33–39. 10.1142/S0218348X99000050

  • CrossRef
  • Google Scholar

• 48

  VincentF.SchwartzL. (2000). Modelisation of 3-D microvasculature by interlaced diffusion limited aggregation. Fractals8.3, 255–259. 10.1142/S0218348X00000317

  • CrossRef
  • Google Scholar

• 49

  WahyudiS.Hadi SetyariniP.SetyariniP. H. (2022).Effects of metabolic heat on the temperature distribution of human hands affected by sarcoma tumors given interstitial hyperthermia therapy”. Kyu. Uni. Inst. Repo.9, 262, 268. 10.5109/4793633

  • CrossRef
  • Google Scholar

• 50

  WeinbaumS.JijiL. M.LemonsD. E. (1984). Theory and experiment for the effect of vascular microstructure on surface tissue heat transfer—Part I: Anatomical foundation and model conceptualization. J. Biomechanical Eng.106 (4), 321–330. 10.1115/1.3138501

  • CrossRef
  • Google Scholar

• 51

  WeinbaumS.JijiL. M. M. M. (1985). A new simplified bioheat equation for the effect of blood flow on local average tissue temperature. J. Biomechanical Eng.107 (2), 131–139. 10.1115/1.3138533

  • CrossRef
  • Google Scholar

• 52

  WittenT. A.SanderL. M. (1981). Diffusion-limited aggregation, a kinetic critical phenomenon. Phys. Rev. Lett.47 (19), 1400–1403. 10.1103/PhysRevLett.47.1400

  • CrossRef
  • Google Scholar

• 53

  XingS.ShinJ.PursleyJ.Correa-AlfonsoC. M.DepauwN.DomalS.et al (2022). A dynamic blood flow model to compute absorbed dose to circulating blood and lymphocytes in liver external beam radiotherapy”. In: Phys. Med. Biol.67.4 p. 045010. 10.1088/1361-6560/ac4da4

  • CrossRef
  • Google Scholar

• 54

  ZhangYuWangH.OliveiraR. H. M.ZhaoC.PopelA. S. (2022). Systems biology of angiogenesis signaling: computational models and omics”. In: WIREs Mech. Dis.14.4 pp. 2692–9368. 10.1002/wsbm.1550

  • CrossRef
  • Google Scholar