This study aimed to reduce the risk of graft occlusion by evaluating the two-phase flow of blood and LDL nanoparticles in coronary artery grafts. The study considered blood as an incompressible ...Newtonian fluid, with the addition of LDL nanoparticles, and the artery wall as a porous medium. Two scenarios were compared, with constant inlet velocity (CIV) and other with pulsatile inlet velocity (PIV), with LDL nanoparticles experiencing drag, wall-induced lift, and induced Saffman lift forces, or drag force only. The study also evaluated the concentration polarization of LDLs (CP of LDLs) near the walls, by considering the artery wall with and without permeation. To model LDL nanoparticles, the study randomly injected 100, 500, and 1000 nanoparticles in three release states at each time step, using different geometries. Numerical simulations were performed using COMSOL software, and the results were presented as relative collision of nanoparticles to the walls in tables, diagrams, and shear stress contours. The study found that a graft implantation angle of 15° had the most desirable conditions compared to larger angles, in terms of nanoparticle collision with surfaces and occlusion. The nanoparticle release modes behaved similarly in terms of collision with the surfaces. A difference was observed between CIV and PIV. Saffman lift and wall-induced lift forces having no effect, possibly due to the assumption of a porous artery wall and perpendicular outlet flow. In case of permeable artery walls, relative collision of particles with the graft wall was larger, suggesting the effect of CP of LDLs.
In this article, a simplified 3D magneto-thermal model to study the temperature distribution over the housing of a Switched Reluctance Machine (SRM) is proposed. The main objective of such model is ...to study the temperature distribution profile on the body of the machine under various working conditions in order to introduce a temperature based signature to develop condition monitoring procedures. The analysis of fluid flow and heat transfer inside the SRM supports the hypothesis that the effect of turbulent flow of the fluid and the complex heat transfer phenomenon of the rotating parts can be simplified to a heat flux boundary condition representation of the heat loss. The accuracy of the model is validated through simulation and experimental results on a prototype SRM. The developed simplified model can be used as for prediction of temperature profile on SRM surface, well suited for fault diagnosis in SRM. The proposed methodology can be modified applied in developing a thermal model for other types of electric machines.
Laminar forced convection heat transfer of water–Cu nanofluids in a microchannel was studied utilizing the lattice Boltzmann method (LBM). The entering flow was at a lower temperature compared to the ...microchannel walls. Simulations were performed for nanoparticle volume fractions of 0.00 to 0.04 and slip coefficient from 0.005 to 0.02. The model predictions were found to be in good agreement with earlier studies. The effects of wall slip velocity and temperature jump of the nanofluid were studied for the first time by using lattice Boltzmann method. Streamlines, isotherms, longitudinal variations of Nusselt number, slip velocity and temperature jump as well as velocity and temperature profiles for different cross sections were presented. The results indicate that LBM can be used to simulate forced convection for the nanofluid micro flows. Moreover, the effect of the temperature jump on the heat transfer rate is significant. Also, the results showed that decreasing the values of slip coefficient enhances the convective heat transfer coefficient and consequently the Nusselt number (Nu) but increases the wall slip velocity and temperature jump values.
Deviation from the optimal bifurcation structure causes abnormal hemodynamic stress, increasing the risk of cerebral aneurysm initiation at the arterial bifurcation apexes of the Circle of Willis. ...Although Murray’s law describes the optimal relationship between bifurcation calibers, major arterial bifurcations within the Circle of Willis show deviations from this law. This study introduces a novel scaling law that describes the optimum relationship between bifurcation characteristics based on pulsatile flow and the internal surface of vessels. The proposed scaling law applies to major intracranial arteries, such as the basilar, internal carotid and common carotid arteries, encompassing both symmetrical and asymmetrical bifurcations. One of the merits of this scaling law is its sole dependence on the Womersley number of parent vessels to determine bifurcation characteristics. The diameter ratios suggested by these relationships are in good agreement with available clinical morphometric data. Numerical simulations of pulsatile flow for several Womersley numbers indicate that the flow resistance and temperature stability of the proposed scaling law are preferable to those of Murray’s law. That might be the reason this scaling law is the optimality principle governing the major cerebral arteries, particularly those arterial blood vessels responsible for the brain’s thermoregulatory, because the brain’s thermoregulatory and temperature stability are the physiological and anatomical constraints of the human brain.
Specifying exact geometry of vessel network and its effect on temperature distribution in living tissues is one of the most complicated problems of the bioheat field. In this paper, the effects of ...blood vessels on temperature distribution in a skin tissue subjected to various thermal therapy conditions are investigated. Present model consists of counter-current multilevel vessel network embedded in a three-dimensional triple-layered skin structure. Branching angles of vessels are calculated using the physiological principle of minimum work. Length and diameter ratios are specified using length doubling rule and Cube law, respectively. By solving continuity, momentum and energy equations for blood flow and Pennes and modified Pennes bioheat equations for the tissue, temperature distributions in the tissue are measured. Effects of considering modified Pennes bioheat equation are investigated, comprehensively. It is also observed that blood has an impressive role in temperature distribution of the tissue, especially at high temperatures. The effects of different parameters such as boundary conditions, relaxation time, thermal properties of skin, metabolism and pulse heat flux on temperature distribution are investigated. Tremendous effect of boundary condition type at the lower boundary is noted. It seems that neither insulation nor constant temperature at this boundary can completely describe the real physical phenomena. It is expected that real temperature at the lower levels is somewhat between two predicted values. The effect of temperature on the thermal properties of skin tissue is considered. It is shown that considering temperature dependent values for thermal conductivity is important in the temperature distribution estimation of skin tissue; however, the effect of temperature dependent values for specific heat capacity is negligible. It is seen that considering modified Pennes equation in processes with high heat flux during low times is significant.
•Effects of blood vessels on temperature distribution of a skin are investigated.•Relaxation time, metabolism and pulse heat flux effects on temperature are studied.•Temperature dependency of thermal conductivity affects skin temperature remarkably.
The dynamics of tumor growth and associated events cover multiple time and spatial scales, generally including extracellular, cellular and intracellular modifications. The main goal of this study is ...to model the biological and physical behavior of tumor evolution in presence of normal healthy tissue, considering a variety of events involved in the process. These include hyper and hypoactivation of signaling pathways during tumor growth, vessels’ growth, intratumoral vascularization and competition of cancer cells with healthy host tissue. The work addresses two distinctive phases in tumor development—the avascular and vascular phases—and in each stage two cases are considered—with and without normal healthy cells. The tumor growth rate increases considerably as closed vessel loops (anastomoses) form around the tumor cells resulting from tumor induced vascularization. When taking into account the host tissue around the tumor, the results show that competition between normal cells and cancer cells leads to the formation of a hypoxic tumor core within a relatively short period of time. Moreover, a dense intratumoral vascular network is formed throughout the entire lesion as a sign of a high malignancy grade, which is consistent with reported experimental data for several types of solid carcinomas. In comparison with other mathematical models of tumor development, in this work we introduce a multiscale simulation that models the cellular interactions and cell behavior as a consequence of the activation of oncogenes and deactivation of gene signaling pathways within each cell. Simulating a therapy that blocks relevant signaling pathways results in the prevention of further tumor growth and leads to an expressive decrease in its size (82% in the simulation).
Molecular dynamics simulations of static argon gas at three different levels of rarefaction are conducted for a channel of 5.4
nm
height to investigate the simultaneous effect of the wall force ...field and the gas temperature on the stress distribution along the channel height. Using the interactive thermal wall model, different temperatures are applied on the channel walls to be able to investigate the effect of the wall temperature and the induced heat flux through the gas medium on the stress distribution. Considering the monoatomic neutral argon gas, the kinetic, particle-particle virial, and surface-particle virial are considered for computing the stress distribution along the channel height. The normal stress components in the bulk gas region are distributed isotropically regardless of the gas density, temperature, and induced heat flux through the domain, while an anisotropy is observed due to the presence of the surface-particle virial. As the gas becomes hotter, the velocity of the gas atoms increases, and thus the kinetic stress component also increases. Besides, the gas density in the wall force field region reduces which eventually attenuates the surface-particle and particle-particle virial stress within 1
nm
from each wall. This effect was also observed as the gas becomes cooler. It is shown that the combination of gas density, wall temperature, and induced heat flux are the main parameters which determine the distribution of stress within the gas medium especially in the wall force field region where repulsive and attractive interactions exist.
This article, examines the flow of argon inside a nanochannel with respect to the molecular dynamics (MD) in the free molecular flow regime using LAMMPS software. The nanochannel is made of copper ...featuring a square cross-section and obstacles of varying dimensions and values. In this study, the flow of argon fluid is three-dimensional. To gain a deeper understanding of the effect of solid walls within the nanochannel and their influence on flow behavior, the research is simulated in a nanochannel with all side walls for the 3D model and without side walls for the 2D model. This research assesses the effect of the obstacles’ dimensions and values on the nanochannel wall surface and areas above the wall surface. The total dimensions of all simulated two- and three-dimensional atomic structures with a square cross-section are assumed to be 60 × 60 × 100 Å3. and the presence of square obstacles (with dimensions of 8 × 8 × 8 Å3) and rectangular obstacles (with dimensions of 8 × 18 × 8 Å3) is examined. This study seeks to understand the influence on flow behavior, temperature distribution, density, heat flux, velocity, and thermal conductivity coefficient. This study is simulated using a time step of 1 fs for 10,000 time steps, involving approximately 10,000–15,000 argon and copper atoms. The results of this research indicate that obstacles with structures of P and R and larger dimensions increase the number of solid atoms exhibiting stronger attractive forces. Compared to a smooth nanochannel, the thermal exchange between fluid and solid atoms results in a density increase of 17.5 % and 17.3 %, respectively. On the other hand, in the 3D nanochannel, the sidewalls of the nanochannel have reduced the effect of the presence of R and P obstacles with larger dimensions, which comparing to a smooth nanochannel, have increased the density by 8.21 % and 7.53 %, respectively. The obstacles with different spatial positions (P and R structures) in the two-dimensional nanochannel cause a rise in the thermal conductivity coefficient. The P structure obstacles have a better effect on the thermal conductivity coefficient in the 2D nanochannel compared to the R structure. In the three-dimensional nanochannel, utilizing smaller obstacles proves to be more effective because it results in better atom distribution or temperature distribution due to increased atomic collisions in the central region compared to the wall regions.
This paper aims to study the gravity effects on the mixed convection heat transfer in a microchannel using lattice Boltzmann method. To include these effects, hydrodynamic boundary condition ...equations are modified. In this problem, cold fluid enters the microchannel and leaves it after cooling the hot walls. For a wide range of inlet Knudsen number (Kn), computations are performed, and for validation, appropriate comparisons between present and previous available results are made.
As the results, stream lines, longitudinal variations of friction coefficient, Nusselt number, slip velocity and temperature jump, and velocity and temperature profiles in different cross sections are presented. The results show that lattice Boltzmann method can be used to simulate mixed convection in a microchannel, and the effects of buoyancy forces are important for Kn < 0.05, specially for hydrodynamic properties, and thus should be included. For Kn > 0.05, these effects can be ignored. In addition, it is observed that buoyancy forces generate a rotational cell in the microchannel flow, leading to the negative slip velocity at Kn = 0.005.
► For Kn < 0.05 the effects of gravity are important for hydrodynamic properties. ► The negative slip velocity was seen, for the first time in this work. ► Mixed convection causes large variations in slip velocity and friction coefficient. ► To use LBM, the hydrodynamic boundary conditions were modified.
In this study, the structure of the fluid flow and vortices are investigated for the circular Couette flow, Taylor vortex and wavy vortex regimes by changing the Womersley number and periods of ...oscillations of the inner cylinder rotation. The rotational velocity of the inner cylinder,
ω
, increases linearly with time from zero to a maximum value and then decreases to zero in a periodic manner. The Womersley number,
Wo
, varies between
0.27
≤
W
o
≤
12.14
,
and the flow is assumed laminar. The results consist of two separate parts as the inner cylinder rotates with 1) positive and 2) negative accelerations. It is observed that by increasing
Wo
, the flow cannot follow the given boundary condition and there is a time delay. Therefore, the values of the critical Taylor number become different from those of the steady-state conditions. For
Wo
=0.27 and when the inner cylinder rotates with positive acceleration, the primary critical Taylor number is 40% higher than that of steady state, while as the inner cylinder rotates with negative acceleration, this is only 7.2% higher than the corresponding steady-state condition. It should be noted that for two values of the Womersley numbers,
Wo
=12.14 and
Wo
=8.59, no flow instability occurs and no vortex appears until the second period of the inner cylinder oscillations. The reason is that the timescale of the dynamics of flow is lower than the timescale of the flow instability; thus, the flow is circular Couette without any vortex.