An advanced computational bioheat transfer model for a human body with an embedded systemic circulation

In the present work, an elaborate one-dimensional thermofluid model for a human body is presented. By contrast to the existing pure conduction-/perfusion-based models, the proposed methodology couples the arterial fluid dynamics of a human body with a multi-segmental bioheat model of surrounding solid tissues. In the present configuration, arterial flow is included through a network of elastic vessels. More than a dozen solid segments are employed to represent the heat conduction in the surrounding tissues, and each segment is constituted by a multilayered circular cylinder. Such multi-layers allow flexible delineation of the geometry and incorporation of properties of different tissue types. The coupling of solid tissue and fluid models requires subdivision of the arterial circulation into large and small arteries. The heat exchange between tissues and arterial wall occurs by convection in large vessels and by perfusion in small arteries. The core region, including the heart, provides the inlet conditions for the fluid equations. In the proposed model, shivering, sweating, and perfusion changes constitute the basis of the thermoregulatory system. The equations governing flow and heat transfer in the circulatory system are solved using a locally conservative Galerkin approach, and the heat conduction in the surrounding tissues is solved using a standard implicit backward Euler method. To investigate the effectiveness of the proposed model, temperature field evolutions are monitored at different points of the arterial tree and in the surrounding tissue layers. To study the differences due to flow-induced convection effects on thermal balance, the results of the current model are compared against those of the widely used modelling methodologies. The results show that the convection significantly influences the temperature distribution of the solid tissues in the vicinity of the arteries. Thus, the inner convection has a more predominant role in the human body heat balance than previously thought. To demonstrate its capabilities, the proposed new model is used to study different scenarios, including thermoregulation inactivity and variation in surrounding atmospheric conditions.     

Interstitial Fluid Flow and Drug Delivery in Vascularized Tumors: A Computational Model

Interstitial fluid is a solution that bathes and surrounds the human cells and provides them with nutrients and a way of waste removal. It is generally believed that elevated tumor interstitial fluid pressure (IFP) is partly responsible for the poor penetration and distribution of therapeutic agents in solid tumors, but the complex interplay of extravasation, permeabilities, vascular heterogeneities and diffusive and convective drug transport remains poorly understood. Here we consider–with the help of a theoretical model–the tumor IFP, interstitial fluid flow (IFF) and its impact upon drug delivery within tumor depending on biophysical determinants such as vessel network morphology, permeabilities and diffusive vs. convective transport. We developed a vascular tumor growth model, including vessel co-option, regression, and angiogenesis, that we extend here by the interstitium (represented by a porous medium obeying Darcy''s law) and sources (vessels) and sinks (lymphatics) for IFF. With it we compute the spatial variation of the IFP and IFF and determine its correlation with the vascular network morphology and physiological parameters like vessel wall permeability, tissue conductivity, distribution of lymphatics etc. We find that an increased vascular wall conductivity together with a reduction of lymph function leads to increased tumor IFP, but also that the latter does not necessarily imply a decreased extravasation rate: Generally the IF flow rate is positively correlated with the various conductivities in the system. The IFF field is then used to determine the drug distribution after an injection via a convection diffusion reaction equation for intra- and extracellular concentrations with parameters guided by experimental data for the drug Doxorubicin. We observe that the interplay of convective and diffusive drug transport can lead to quite unexpected effects in the presence of a heterogeneous, compartmentalized vasculature. Finally we discuss various strategies to increase drug exposure time of tumor cells.     

Mass transport in a microchannel bioreactor with a porous wall

A two-dimensional flow model has been developed to simulate mass transport in a microchannel bioreactor with a porous wall. A two-domain approach, based on the finite volume method, was implemented. For the fluid part, the governing equation used was the Navier-Stokes equation; for the porous medium region, the generalized Darcy-Brinkman-Forchheimer extended model was used. For the porous-fluid interface, a stress jump condition was enforced with a continuity of normal stress, and the mass interfacial conditions were continuities of mass and mass flux. Two parameters were defined to characterize the mass transports in the fluid and porous regions. The porous Damkohler number is the ratio of consumption to diffusion of the substrates in the porous medium. The fluid Damkohler number is the ratio of the substrate consumption in the porous medium to the substrate convection in the fluid region. The concentration results were found to be well correlated by the use of a reaction-convection distance parameter, which incorporated the effects of axial distance, substrate consumption, and convection. The reactor efficiency reduced with reaction-convection distance parameter because of reduced reaction (or flux), and smaller local effectiveness factor due to the lower concentration in Michaelis-Menten type reactions. The reactor was more effective, and hence, more efficient with the smaller porous Damkohler number. The generalized results could find applications for the design of bioreactors with a porous wall.     

Two-dimensional standing gradient osmotic flow: a generalization of the “isotonic convection approximation”

A linear two-dimensional model of a flow system for solute and fluid transport in intercellular spaces has been obtained by using the so called isotonic convection approximation, already employed in the one-dimensional case (Segel (1970)). This is equivalent to ignoring the convective components in the relevant differential equation. The solutions found by means of the eigenfunctions of a Sturm-Liouville system have been compared with the numerical solutions of the general non-linear model, in which also the convective terms are present.The results of the linear model agree fairly well with those of the non-linear one, in the range of interest of parameters. This fact shows that the ignorance of the profile of fluid velocity does not introduce significant errors in the evaluation of solute and fluid fluxes.The differences of results between the one-dimensional and the twodimensional, the linear and non-linear models, give a way to evaluate the relative effects of the diffusive and convective terms in the differential equation and to estimate the errors introduced by the one-dimensional approximation.     

Mathematical studies of capillary-tissue exchange

A mathematical model of capillary-tissue exchange is presented and the method of solution of the resulting equations is described. The model includes the mutual interaction of fluid movement across the capillary wall and the convection and diffusion of a number of solutes. A variety of solutions for situations of physiological interest are obtained and discussed.     

The frequency-dependent Wright–Fisher model: diffusive and non-diffusive approximations

We study a class of processes that are akin to the Wright–Fisher model, with transition probabilities weighted in terms of the frequency-dependent fitness of the population types. By considering an approximate weak formulation of the discrete problem, we are able to derive a corresponding continuous weak formulation for the probability density. Therefore, we obtain a family of partial differential equations for the evolution of the probability density, and which will be an approximation of the discrete process in the joint large population, small time-steps and weak selection limit. If the fitness functions are sufficiently regular, we can recast the weak formulation in a more standard formulation, without any boundary conditions, but supplemented by a number of conservation laws. The equations in this family can be purely diffusive, purely hyperbolic or of convection–diffusion type, with frequency dependent convection. The particular outcome will depend on the assumed scalings. The diffusive equations are of the degenerate type; using a duality approach, we also obtain a frequency dependent version of the Kimura equation without any further assumptions. We also show that the convective approximation is related to the replicator dynamics and provide some estimate of how accurate is the convective approximation, with respect to the convective-diffusion approximation. In particular, we show that the mode, but not the expected value, of the probability distribution is modelled by the replicator dynamics. Some numerical simulations that illustrate the results are also presented.     

A continuum model of protrusion of pseudopod in leukocytes.

The morphology of human leukocytes, the biochemistry of actin polymerization, and the theory of continuum mechanics are used to model the pseudopod protrusion process of leukocytes. In the proposed model, the pseudopod is considered as a porous solid of F-actin network, the pores of which are full of aqueous solution. G-actin is considered as a "solute" transported by convection and diffusion in the fluid phase. The pseudopod grows as actin filaments elongate at their barbed ends at the tip of the pseudopod. The driving force of extension is hypothesized as being provided by the actin polymerization. It is assumed that elongation of actin filaments, powered by chemical energy liberated from the polymerization reaction, does mechanical work against opposing pressure on the membrane. This also gives rise to a pressure drop in the fluid phase at the tip of the pseudopod, which is formulated by an equation relating the work done by actin polymerization to the local state of pressure. The pressure gradient along the pseudopod drives the fluid filtration through the porous pseudopod according to Darcy's Law, which in turn brings more actin monomers to the growing tip. The main cell body serves as a reservoir of G-actin. A modified first-order equation is used to describe the kinetics of polymerization. The rate of pseudopod growth is modulated by regulatory proteins. A one-dimensional moving boundary problem based on the proposed mechanism has been constructed and approximate solutions have been obtained. Comparison of the solutions with experimental data shows that the model is compatible with available observations. The model is also applicable to growth of other cellular systems such as elongation of acrosomal process in sperm cells.     

Mixed Convection Flow of Viscoelastic Fluid by a Stretching Cylinder with Heat Transfer

Flow of viscoelastic fluid due to an impermeable stretching cylinder is discussed. Effects of mixed convection and variable thermal conductivity are present. Thermal conductivity is taken temperature dependent. Nonlinear partial differential system is reduced into the nonlinear ordinary differential system. Resulting nonlinear system is computed for the convergent series solutions. Numerical values of skin friction coefficient and Nusselt number are computed and discussed. The results obtained with the current method are in agreement with previous studies using other methods as well as theoretical ideas. Physical interpretation reflecting the contribution of influential parameters in the present flow is presented. It is hoped that present study serves as a stimulus for modeling further stretching flows especially in polymeric and paper production processes.     

Thermal convection over East Antarctica: Potential microorganism dispersal

The Indian Antarctic station, Maitri is in the Schirmacher oasis of east Antarctica. The oasisis covered with snow/ice; except for the local summer season when it gets deglaciated and exposes the small hilly region. During summer, minute microorganisms are observed near water bodies of this rocky terrain. In the year 1996a monostatic acoustic sounder probed planetary boundary layer dynamics over this region. From the data it is observed that the thermal convection (both free and forced) persist for 5.63% of the time in the whole year. The occurrence of free convection predominates in the local summer season around midday, while sporadic cases of forced convection have been recorded during the autumn and winter seasons. The onset of convection is mainly at 0600 hrs to 1200 hrs, while the cessation period is limited within 1400 hrs to 1900 hrs. The cessation of convection is basically governed by the katabatic wind flow around the Schirmacher region and it indicates that the interior of Antarctica or the polar cap ice becomes cooler much faster than the rocky region of Schirmacher oasis.Examples taken from the literatures on atmospheric structure and their effects on dispersal of microorganisms and their distribution by the wind are discussed. The study of monostatic acoustic sounder for thermal convections/plumes may form input for the study of dispersal, survival, metabolicactivities and dispersion model of microorganisms. The application of convection/plume to aerobiology can also lead to improvements in forecasting, monitoring and understanding the various mechanisms of organism dispersal.     

Effect of tumor shape,size, and tissue transport properties on drug delivery to solid tumors

Background

The computational methods provide condition for investigation related to the process of drug delivery, such as convection and diffusion of drug in extracellular matrices, drug extravasation from microvessels or to lymphatic vessels. The information of this process clarifies the mechanisms of drug delivery from the injection site to absorption by a solid tumor. In this study, an advanced numerical method is used to solve fluid flow and solute transport equations simultaneously to investigate the effect of tumor shape and size on drug delivery to solid tumor.

Methods

The advanced mathematical model used in our previous work is further developed by adding solute transport equation to the governing equations. After applying appropriate boundary and initial conditions on tumor and surrounding tissue geometry, the element-based finite volume method is used for solving governing equations of drug delivery in solid tumor. Also, the effects of size and shape of tumor and some of tissue transport parameters such as effective pressure and hydraulic conductivity on interstitial fluid flow and drug delivery are investigated.

Results

Sensitivity analysis shows that drug delivery in prolate shape is significantly better than other tumor shapes. Considering size effect, increasing tumor size decreases drug concentration in interstitial fluid. This study shows that dependency of drug concentration in interstitial fluid to osmotic and intravascular pressure is negligible.

Conclusions

This study shows that among diffusion and convection mechanisms of drug transport, diffusion is dominant in most different tumor shapes and sizes. In tumors in which the convection has considerable effect, the drug concentration is larger than that of other tumors at the same time post injection.

     

Influence of interstitial fluid dynamics on growth and therapy of angiogenic tumor. Analysis by mathematical model

We have developed a spatially distributed mathematical model of angiogenic tumor growth in tissue with account of interstitial fluid dynamics and bevacizumab monotherapy. In this model the process of neovascularization is initiated by tumor cells in a state of metabolic stress, vascular endothelial growth factor (VEGF) being its main mediator. The model takes into consideration the convection flows arising in dense tissue due to active proliferation and migration of tumor cells as well as interstitial fluid inflow from blood vascular system, its outflow through lymphatic system and redistribution in the area of tumor growth. The work considers the diffusive approximation of interstitial fluid dynamics in tumor and normal tissue. Numerical study of the model showed that in absence of therapy a peritumoral edema is formed due to the increase of interstitial fluid inflow from angiogenic capillaries. In the case of rapid interstitial fluid outflow through lymphatic system and its fast transport from necrotic zone to normal tissue the regimes of full growth stop are observed in case of low-invasive tumor. Under bevacizumab monotherapy the peritumoral edema vanishes and low-invasive tumor may not only decelerate its growth, but also start shrinking for a large range of parameters.     

Activation and extinction models for platelet adhesion

Adherent platelets are an important part of both thrombus formation and in certain stages of atherogenesis. Platelets can be activated by potent chemicals released from adherent platelets and adhere far more readily than unactivated ones. An analytical and numerical model is presented utilising high Peclet number for the activation and adhesion of platelets in shear flows. The model uses a similarity transformation, which characterises the relationship between convective, diffusive transport and the bulk platelet activating reaction mechanism. A first order surface reaction mechanism is used to model platelet adhesion at the wall (cell) surface. The reduced Damk?hler number, M, characterises the importance of the bulk reaction and includes both convective and diffusive terms. For a high rate of blood flow (M-->0) the activation of platelets can effectively be terminated. In contrast, for (M-->infinity) an inner layer of activated platelets exists with an infinitesimally thin reaction sheet separating activated and non-activated platelets. This characterisation by the Damk?hler number highlights results found clinically, in that thrombus forms in areas of low shear (high M) and in some cases an increased blood flow (low M) can inhibit the activation of platelets completely. The model shows the critical balance that exists between convection, diffusion and reaction.     

Fluid flow and convective transport of solutes within the intervertebral disc

Previous experimental and analytical studies of solute transport in the intervertebral disc have demonstrated that for small molecules diffusive transport alone fulfils the nutritional needs of disc cells. It has been often suggested that fluid flow into and within the disc may enhance the transport of larger molecules. The goal of the study was to predict the influence of load-induced interstitial fluid flow on mass transport in the intervertebral disc.An iterative procedure was used to predict the convective transport of physiologically relevant molecules within the disc. An axisymmetric, poroelastic finite-element structural model of the disc was developed. The diurnal loading was divided into discrete time steps. At each time step, the fluid flow within the disc due to compression or swelling was calculated. A sequentially coupled diffusion/convection model was then employed to calculate solute transport, with a constant concentration of solute being provided at the vascularised endplates and outer annulus. Loading was simulated for a complete diurnal cycle, and the relative convective and diffusive transport was compared for solutes with molecular weights ranging from 400 Da to 40 kDa.Consistent with previous studies, fluid flow did not enhance the transport of low-weight solutes. During swelling, interstitial fluid flow increased the unidirectional penetration of large solutes by approximately 100%. Due to the bi-directional temporal nature of disc loading, however, the net effect of convective transport over a full diurnal cycle was more limited (30% increase). Further study is required to determine the significance of large solutes and the timing of their delivery for disc physiology.     

Wavelengths of gyrotactic plumes in bioconvection

Bioconvection occurs as the result of the collective behaviour of many microorganisms swimming in a fluid and is realized as patterns similar to those of thermal convection which occur when a layer of fluid is heated from below. We consider the phenomenon of pattern formation due to gyrotaxis, an orientation mechanism which results from the balance of gravitational and viscous torques acting on bottom-heavy micro-organisms. The continuum model of Pedley et al. (1988, J. Fluid. Mech. 195, 223–237) is used to describe the suspension. The system is governed by the Navier-Stokes equations for an incompressible fluid coupled with a micro-organism conservation equation. These equations are solved numerically using a conservative finite-difference scheme. To examine the dependence of the horizontal pattern wavelengths on the parameters, we consider two-dimensional solutions in a wide chamber using rigid side walls. The wavelengths of the numerical computations are in good agreement with the experimental observations and we provide the first computational examples of the commonly seen ‘bottom-standing’ plumes.     

Mass transport in a microchannel enzyme reactor with a porous wall: Hydrodynamic modeling and applications

A two-dimensional flow model, incorporating mass transport, has been developed to simulate a microchannel enzyme reactor with a porous wall. A two-domain approach based on the finite volume method was implemented. Two parameters are defined to characterize the mass transports in the fluid and porous regions: the porous Damkohler number and the fluid Damkohler number. For reactions close to first-order type (enzyme reactor), the concentration results are found to be well correlated by the use of a reaction–convection distance parameter which incorporates the effects of axial distance, substrate consumption and convection. The reactor efficiency reduces with reaction–convection distance parameter because of reduced reaction (or flux) due to the lower concentration. Increased fluid convection improves the efficiency but it is limited by the diffusion in the fluid region. The correlated results can find applications for the design of enzyme reactors with a porous wall.     

Singular perturbation analysis of the nonlinear, flow-dependent compressive stress relaxation behavior of articular cartilage

The dominant mechanism giving rise to the viscoelastic response of articular cartilage during compression is the nonlinear diffusive interaction of the fluid and solid phases of the tissue as they flow relative to one another. The present study is concerned with the role of this interaction under uniaxial stress relaxation in compression. The model is a biphasic mixture of fluid and solid which incorporates the strain-dependent permeability found earlier from permeation experiments. When a ramp-displacement is imposed on the articular surface, simple, but accurate, asymptotic approximations are derived for the deformation and stress fields in the tissue for slow and moderately fast rates of compression. They are shown to agree very well with experiment and they provide a simple means for determining the material parameters. Moreover, they lead to important insights into the role of the flow-dependent viscoelastic nature of articular cartilage and other hydrated biological tissues.     

Signaling in morphogenesis: transport cues in morphogenesis

Extracellular transport processes play critical roles in morphogenesis. While diffusive transport effects on morphogenesis are well illustrated in examples like blood capillary architecture and in cell morphogenetic responses to the local extracellular protein environment, the effects of fluid convection, although important in many developing and regenerating tissues, are not well understood. Convective forces are present whenever a hydrated tissue undergoes dynamic mechanical strain, and so convection could not only dominate the transport of large molecules like proteins, but might also serve as a mechanism for mechanosensing. The complex interdependence of mechanical forces, protein transport and extracellular morphogen gradients needs to be elucidated in a comprehensive way in order for the importance of transport on morphogenesis to be fully appreciated.     

Perfused phantom models of microwave irradiated tissue

The theoretical basis, practical design considerations, and prototype testing of a perfused model suitable for simulation studies of microwave heated tissue are presented. A parallel tube heat exchanger configuration is used to simulate the internal convection effects of blood flow. The global thermal response of the phantom, on a scale of several tube spacings, is shown theoretically to be nearly identical to that predicted by Pennes' bioheat equation, which is known to give a reasonable representation of tissue under many conditions. A parametric study is provided for the relationships between the tube size, spacing and material properties and the simulated perfusion rate. A prototype with a physiologically reasonable perfusion rate was tested using a typical hyperthermia applicator. The measured thermal response of the phantom compares favorably with the numerical solution of the bioheat equation under the same irradiation conditions. This similarity sheds light on the unexpected success of the bioheat equation for modeling the thermal response of real tissue.     

A Thermal Equilibrium Analysis of Line Contact Hydrodynamic Lubrication Considering the Influences of Reynolds Number,Load and Temperature

Thermal effects such as conduction, convection and viscous dissipation are important to lubrication performance, and they vary with the friction conditions. These variations have caused some inconsistencies in the conclusions of different researchers regarding the relative contributions of these thermal effects. To reveal the relationship between the contributions of the thermal effects and the friction conditions, a steady-state THD analysis model was presented. The results indicate that the contribution of each thermal effect sharply varies with the Reynolds number and temperature. Convective effect could be dominant under certain conditions. Additionally, the accuracy of some simplified methods of thermo-hydrodynamic analysis is further discussed.     

Insect eggs at a transition between diffusion and reaction limitation: temperature, oxygen, and water. art.woods@mso.umt.edu

In diverse animal taxa, eggs and embryos are incapable of transporting oxygen by convection. In such cases, internal oxygen distributions are determined jointly by rates of oxygen consumption and diffusive transport. Here we develop a mathematical model of oxygen consumption and transport in insect eggs, with the goal of understanding-for eggs in variable-temperature environments-the interactive effects of the two processes on development. We fit the model to previously published data on development time of eggs of a sphingid moth, Manduca sexta. The fitted coefficients suggest that eggs develop at a transition point between reaction- and diffusion-limitation. We test then this conclusion with independent data on development times of eggs distributed across a set of temperatures generated by a thermal gradient bar. Finally, we develop an extension of the model that considers tradeoffs between oxygen transfer to eggs versus water loss from them. The model results provide both a rationale for why development is often mass-transfer limited and a set of new predictions about oxygen-water tradeoffs.