Next Article in Journal
Automatic Illumination Control Method for Indoor Luminaires Based on Multichromatic Quantum Dot Light-Emitting Diodes
Next Article in Special Issue
Fabrication and Evaluation of Tubule-on-a-Chip with RPTEC/HUVEC Co-Culture Using Injection-Molded Polycarbonate Chips
Previous Article in Journal
Deep Learning Based Feature Selection Algorithm for Small Targets Based on mRMR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical and Experimental Analysis of Shear Stress Influence on Cellular Viability in Serpentine Vascular Channels

1
Department of Biomedical Engineering, National Institute of Technology, Raipur 492010, India
2
Biomedical Engineering, Florida Tech, Melbourne, FL 32901, USA
*
Authors to whom correspondence should be addressed.
Micromachines 2022, 13(10), 1766; https://doi.org/10.3390/mi13101766
Submission received: 11 August 2022 / Revised: 19 September 2022 / Accepted: 14 October 2022 / Published: 18 October 2022

Abstract

:
3D bioprinting has emerged as a tool for developing in vitro tissue models for studying disease progression and drug development. The objective of the current study was to evaluate the influence of flow driven shear stress on the viability of cultured cells inside the luminal wall of a serpentine network. Fluid–structure interaction was modeled using COMSOL Multiphysics for representing the elasticity of the serpentine wall. Experimental analysis of the serpentine model was performed on the basis of a desirable inlet flow boundary condition for which the most homogeneously distributed wall shear stress had been obtained from numerical study. A blend of Gelatin-methacryloyl (GelMA) and PEGDA200 PhotoInk was used as a bioink for printing the serpentine network, while facilitating cell growth within the pores of the gelatin substrate. Human umbilical vein endothelial cells were seeded into the channels of the network to simulate the blood vessels. A Live-Dead assay was performed over a period of 14 days to observe the cellular viability in the printed vascular channels. It was observed that cell viability increases when the seeded cells were exposed to the evenly distributed shear stresses at an input flow rate of 4.62 mm/min of the culture media, similar to that predicted in the numerical model with the same inlet boundary condition. It leads to recruitment of a large number of focal adhesion point nodes on cellular membrane, emphasizing the influence of such phenomena on promoting cellular morphologies.

1. Introduction

Small diameter blood vessels followed by upstream vasoconstriction is often prone to secondary stenosis [1]. Hence, regenerated blood vessels will be an appropriate solution for replacing these small diameter vessels (which were affected with secondary stenosis). Various techniques, such as extrusion, electrospinning, thermal-induced phase separation, braiding, hydrogel tubing, and gas foaming had been used for developing artificial blood vessels [2,3,4,5]. Serpentine vascular geometry having a rigid wall has been used in vascular regeneration [6]. Unlike a straight channel blood vessel, both bifurcated and arch-type blood vessels do not facilitate constant and unidirectional blood flow [7]. Instead, the flow becomes turbulent with larger eddies and torsion. Such phenomena encourages formulation and analysis of the flow driven mechanical force distributions within the serpentine vascular network.
Bioprinting methodologies commonly use hydrogels as a cell-laden material as they are able to come into contact with the cells without damaging their viability [8,9,10,11]. Moreover, hydrogels can readily be mixed with cells, while simultaneously allowing for high cell density and homogenous cell distribution throughout the scaffold [12]. The source of the hydrogel is an important factor because of varying chemical and mechanical properties of synthetic and natural polymers [13]. Naturally derived hydrogels, such as GelMA, consist of inherent signaling molecules that promote cell adhesion. On the other hand, hydrogels derived from other organisms such, as alginate, lack these binding sites for cell adhesion and attachment [14]. Natural biomaterials regulate interactions between the cells and the extracellular matrix and provide an excellent environment for the cells to grow, differentiate, and proliferate [13,15]. Since GelMa is a gelatin based bioink, it can also resolve intricate vascular networks and channels that offer endothelial cells with the essential properties of their native environments. Rheological modifiers also called rheological additives are mixed with gelatin based bioinks to enhance the rheology of the resulting bioink, especially with regard to viscosity and yield stress.
The pulsatile and sinusoidal flows have been reported as the simplified form of physiological blood flow in circulatory system [16,17,18]. These inlet velocity profiles represent simulated biomimetic flow conditions for the serpentine vessels [19,20,21,22]. Further, the flow facilitates the seeded endothelial cells with the physiological stress at the luminal wall of the serpentine structure [23]. Limited study can be found dealing with both physiological and sinusoidal flow to generate hemodynamic stress on the annular surface of the straight and bifurcated substrates [24,25,26]. However, such inlet boundary conditions were not considered in a serpentine network with an elastic wall. Previous studies had considered serpentine models independent of the upstream stenosis (or secondary stenosis). The physiological flow of blood had been found to reproduce adequate stress for the seeded endothelial cells on the luminal wall, while the rigid model of the wall compensates the viability of such cells under a shear stress condition [27,28]. This stress was generated due to the flow of culture media through the bioprinted channels. Due to the positioning of endothelial cells, it is more sensitive and responds quickly to the fluid flow driven stress [29]. It had also been reported that generation of vascular growth factors, conservation of blood vessels and migration of endothelial cells were being improved in the presence of modulating wall shear stress [30,31]. Shear stress derived from flow structures were also found to influence interconnected points between the cytoskeleton and extracellular structures of endothelial cells [32,33,34]. However, it was found that elasticity of the vessel wall promotes migratory phenomena of endothelial cells in the presence of lateral stress in the periodic form of lymphatics. Hence, it is very vital to consider the bioprinted serpentine structure with the elastic wall for enhancing the cell viability of seeded cells in the presence of physiological and sinusoidal flow conditions.
Most of the previous studies deal with the systemic development of a 3D scaffold for formulation of straight channel blood vessels. In addition, a perfusion bioreactor has traditionally been used because of its mechanical advantages rather than a rotatory bioreactor. The present study deals with the use of a 3D printer along with a rotatory bioreactor to fabricate a serpentine vascular structure. The use of a rotatory bioreactor facilitates adequate amounts of shear stress for stimulating the growth of human endothelial cells on the vessel walls. Numerical analysis was conducted to enumerate the relationship between flow driven different mechanical parameters and different inlet and wall boundary conditions. Variation of pressure, stress, and shear rate were found influenced by the inlet and outlet boundary conditions as well as by the viscosity of the acting fluid. There is limited data on cellular viability of vascular channels bioprinted using GelMA bioink mixed with PEGDA photoink. Cellular viability of GelMa hydrogels in a 3D cell culture model demonstrated that the hydrogel scaffolds provide a cell promoting environments for mesenchymal stem cells [35]. Experimental characterization of the bioprinted channels has been performed to validate the optimum flow-driven mechanical parameters, while maintaining higher cellular viability.

2. Methodology

A 3D model of the serpentine blood vessel was implemented using COMSOL Multiphysics. Results of the numerical models were further used to optimize the media flow parameters: pressure, shear rate, stress, and axial velocity for different inlet flow boundary conditions. These optimized parameters were further utilized for maintaining the printed serpentine channels in a rotary bioreactor. Finally, a microfluidic based chip system was developed to fabricate the serpentine vascular channels.

2.1. Numerical Analysis

The serpentine channel model is shown in Figure 1. Fluid flow (with Newtonian viscosity model for blood equivalent) was evaluated for different wall boundary conditions [slip (S) and no-slip (NS)]. To ensure that the vascular serpentine channel model reflects the realistic simulation of in vivo conditions, appropriate relevant boundary conditions were used for numerical analysis. The wall of the serpentine vascular channel was taken as a rigid body, with a no-slip boundary condition [36]. Fluid–structure interaction (FSI) physics was used for simulating the elastic behavior of the serpentine channel model in another study. For the simulation, the wall of the serpentine vascular structure was considered as a moving mesh. The following properties of blood were used for numerical analysis: density = 1066 kg/m3, dynamic viscosity = 3.5 × 10−3 Pa.s. Young’s Modulus = 0.4 × 106 Pa, Poisson’s ratio = 0.5. GelMA material properties (density = 1020 kg/m3, dynamic viscosity = 4.2 × 10−3 Pa.s, Young Modulus 3.18 KPa) were used for the modeling of the serpentine wall. Simulations were performed at a reference temperature of 293.15 K. Navier–Stokes equation was used to represent the momentum of the fluid model. Time-dependent partial differential equation of incompressible Navier–Stokes (Equation (1)) was used for the simulation.
ρ u t μ 2 u + ρ u × u + P = 0
× u = 0
where u is the blood velocity and P is the blood pressure, ρ is the blood density and is the blood viscosity. Fluid was considered to be incompressible and hence a continuity equation of incompressible form was used. Linear elastic material properties of solid mechanics and fluid properties of laminar flow were considered for the study. An isotropic solid model was used for both slip and no-slip conditions. A fully coupled option was selected in COMSOL Multiphysics to establish a possible connection between solid mechanics physics and laminar flow physics. In FSI analysis, all of the meshes created were physics controlled and automatically generated at a normal element size. The statistics of mesh for slip and no slip boundary conditions are listed in Table 1.
A transient model was used to realize the transport phenomena of the fluid in a 2D serpentine channel. A MUMPS based direct linear solver was used to solve the model with a damping factor of unity and a recovery damping factor on 0.75. Velocities corresponding to mean arterial pressure: 80 mmHg, 85 mmHg, and 90 mmHg were considered for the study in Microchannel [37,38,39,40].
Three different flow rate based inlet velocity models of sinusoidal flow (SF) [V1 = 4.48 mm/min, V2 = 4.62 mm/min and V3 = 4.76 mm/min] and physiological flow-PF [V4 = 4.48 mm/min, V5 = 4.62 mm/min and V6 = 4.76 mm/min] were considered in the study. The physiological functions (waveform) were derived from 4D-Laser Doppler data of blood flow across a cross-section of sub-clavicular artery of healthy subject. The study was approved by Institute Ethical Committee (NITRR/IEC/2021/12). Obtained tempero-spatial functions were interpolated thereafter using Curve-Fitting MATLAB toolbox version 2018a, Mathworks, India. Thereafter, the interpolated function with lowest RMSE value was selected for the study. The functions corresponding to inlet velocity models for V1 to V6 are given in Table 2. One complete cycle of both sinusoidal and physiological flow corresponds to one complete cardiac cycle (represented by T). Different study conditions were formulated based on the above boundary conditions: no-slip (wall) and sinusoidal inlet flow, slip (wall) and sinusoidal inlet flow, no-slip (wall) and physiological inlet flow, slip (wall) and physiological inlet flow (as tabulated in Table 3). A no-slip boundary condition will exclusively provide the influence of inlet flow type in the development of flow physics of media inside the serpentine. On the other hand, slip wall boundary condition will simulate the elasticity of the serpentine wall (similar to the natural blood vessel wall).
Three different locations corresponding to positions P1 (4, −2.5), P2 (5.8, 2.4) and P3 (7.7, −2.5) on the neck, abdominal, and rear region, respectively, were identified on the serpentine channel. These three regions are the downstream region immediately after curvilinear orientation of ascending and descending phases of the serpentine structure, which are more prone to Coriolis forces within the serpentine model. Flow parameters (pressure, shear rate, stress) over the entire cardiac cycle were evaluated at an interval of T/3, T/2, and T cycles.

Grid Convergence Test

A grid convergence test was conducted for optimizing the grid size of the flow domain. Grid Convergence Index (GCI) provides a uniform measure of convergence for grid refinement study [41]. The discretization of sinusoidal flow for velocity V1 with no slip condition had been selected for three different element sizes. For all element sizes—pressure, shear rate and velocity at probe points P1 and P2 (as marked in Figure 1) were measured for the entire time cycle. In the current study, the grid refinement ratio r, which is equivalent to mean refinement ratio r m e a n of r 12 and r 23 were calculated using Equation (3a)
r = r m e a n = r 12 + r 23 2
where,
r 12 = element value of Fine discretization element value of Normal discretization
r 23 = element value of Normal discretization element value of coarse discretization
Richardson extrapolation [42] introducing the p -th order method, as shown by Equation (3b):
p = l n f 3 f 2 f 2 f 1 / l n r
where f 1 , f 2 and f 3 are magnitude of derived parameters (e.g., velocity) at point P1 and P2 at different grid resolution at time T.
With the different parameters listed in Table 4, values of r and p were evaluated from Equation (3a,b), respectively, and the calculated grid convergence magnitude for two probe points is listed in Table 5. Further, the Grid Convergence Index ( G C I ) was evaluated using the Richardson extrapolation method based on estimated fractional error [41], as shown in Equation (4a):
G C I = F s e r p 1 × 100 %
The safety factor (Fs) selected for this study was considered to be 1.25.
It was observed that GCI for the fine grid ( G C I 12 ) is relatively low if compared to the coarse grid ( G C I 23 ), which demonstrates that the numerical simulation dependency on the cell size was minimum with a decrease in grid size. Figure 2 shows the plot between extrapolated value and the Richardson extrapolation estimation for centerline velocity. The y-axis of Figure 2 represents ƹ i + 1 , i , which is the difference in the values of a particular parameter at two different grid resolution and given by:
ƹ i + 1 , i = f i + 1 f i
It was observed that the extrapolated values for the velocity at point P1 and P2 did not change significantly on further decreasing the mesh size to represent finer grid resolution.

2.2. Experimental Analysis

A co-axial nozzle configuration of a 3D printer was used to print the serpentine channel. The printing bed had been located inside closed chamber housing, as shown in Figure 3. The construct was printed using a blend of GelMA photoInk purchased from CELLINK, Boston USA and PEGDA200 photoInk (purchased from CELLINK, Boston, MA, USA). This photo-ink blend equipped the final construct with cell attachment sites, high shape fidelity, and minimal construct swelling. CaCl2 was used as a chemical cross-linker during printing process of the serpentine. The printing environment parameters are presented in Table 6.

2.2.1. Endothelization of Vascular Network

A total of 50–500 µL of HUVEC suspensions (107 cells per mL) was injected via a micro-pipette to fill the channel. Simultaneously, the inlet and the outlet of the construct was sealed with a pinch-clamp. Channels were injected with HUVECs stained with CellTrace CFSE dye (ThermoFisher, Waltham, MA, USA). The device had been centrifuged at 40 rcf for 1 min with a slow acceleration and deceleration to let the cells settle to the bottom. Incubation at 37 °C facilitated the adhesion of cells to form the innermost layer of the vessels. After incubation for 30 min, the construct was placed over 180° for cell adhesion on the other side of the vessel. Finally, the cells were incubated for 5 h at 37 °C on a shaker.
A blood-equivalent working fluid was put in a 10 mL BD Plastipak syringe. The filled syringe was loaded in the programable Shenchen syringe pump and was operated in infusion mode. A syringe pump was connected to microfluidic accessories (one-way luer lock check valve, microfluidic fitting female luer lock adopter kit, barbed to male adopter and syringe tube) in a series to complete the flow driven system. The other end of the tube was connected to a fabricated microfluidic channel placed in the petri dish via a barbed male adopter, microfluidic fitting female luer lock kit, male luer lock kit, fitting adopter male luer to male thread, and mini luer to luer adopter. All the microfluidic components were procured from DARWIN microfluidics.
HUVECS media was fed to the reservoir for cell culture over the printed layers. Regulatory control of each component, such as media pressure, temperature, valve, pH and CO2 and N2 delivery, etc., were achieved. Induction of media to the housing was made by using a pulsatile pump, as shown in Figure 4. It creates pulsatile hydrostatic pressure to the printed bioink construct, similar to the optimized inlet velocity profiles finalized from numerical study. The printed constructs were imaged using fluorescence microscopy to evaluate cell attachment and cell distribution within the channel.

2.2.2. In-Vitro Testing

Microfluidic channel had been formulated by combining two casted segments. The upper part and the lower part of the printed structure of the model were merged to form a microfluidic system to support a serpentine structure, for housing the printed construct. Media had been transfused through the serpentine structure from the inlet of the structure in the form of sinusoidal waveform and physiological waveform using a programmable syringe pump (DARWIN microfluidic Shenchen Syringe pump, Paris, France), as shown in Figure 4.

2.2.3. Uncertainty Analysis

Uncertainty analysis was performed to validate the model characteristics and to determine the shear stress sensitivity. The uncertainty equation of shear stress on the fabricated channel model was represented by Equation (5). Shear stress was generated when the printed serpentine model came into contact with the flowing HUVECS media.
Δ τ = Δ τ Δ P σ P 2 + Δ τ Δ Y σ Y 2 + Δ τ Δ V σ V 2
Here, σ P , σ V and σ Y are the standard deviations of the pressure, velocity, and the strain rate standard deviation, respectively. From Equation (5), the uncertainty value of stress equivalent was evaluated as ±179.36501 N/m2.

2.2.4. Cellular Viability Assessment

The viability of the bioprinted vascular channel was evaluated over a period of 14 days by performing LIVE-DEAD assay. The channels were stained with the LIVE-DEAD cell imaging kit (Fisher Scientific, Waltham, MA, USA). Live cells were stained with Calcein-AM (2 μM), the fluorescence was measured at 494–517 nm, and dead cells were stained with ethidium homodimer-1 (4 μM), which was measured at 528–617 nm. The channels were incubated with LIVE-DEAD stain for 30 min and the samples were washed with PBS. Imaging was then performed using a Nikon C1 Confocal fluorescence microscope. Image J 1.53k open source software developed by Wayne Rasband and contributors National institute of Health, USA was used to evaluate the cell viability.

2.2.5. Sensitivity Analysis

The objective of the sensitivity analysis was to obtain the change of output with respect to change in input parameters. Input parameters may be a material property, variation in geometry and loading condition etc. For a proposed serpentine geometry, sensitivity analysis was performed for a no-slip wall condition of physiological flow corresponding to velocity model V6. The sensitivity of the serpentine model was calculated using Equation (6)
p = v 2 ρ 2
where p is the calculated pressure, v is the maximum velocity (4.76 mm/min), ρ is the blood density 1066 kg/m3. Change in the pressure due to change in density can be expressed by the Equation (7) (modified version of Equation (6))
p = v 2 ρ 2
Sensitivity = p ρ = v 2 2
The analytical value of sensitivity was calculated by Equation (8)
p ρ = 4.76 × 4.76 2
Sensitivity(analytical) = 11.32
The computed sensitivity was calculated from pressure obtained at blood density (1066 Kg/m3) and 6% of blood density (997 kg/m3) at a point P1, respectively, as given in Equation (9)
Sensitivity ( computed ) = p ρ
Sensitivity ( computed ) = 12050 11300 1066 997
Sensitivity(computed) = 10.86
The error of sensitivity was calculated between analytical sensitivity and computed sensitivity by below Equation (10)
Error ( sensitivity ) = Sensitivity computed Sensitivity analytical Sensitivity analytical
Error ( sensitivity ) = 11.32 10.86 11.32
Error (sensitivity) = 4.06%
At 6% of change in input parameter (density), it produces a minimum error of sensitivity [43]. From the above analysis it was concluded that an independent material property (density) of the flowing fluid influenced the dependent parameters (i.e., pressure) with a 6% change in density of fluid, as shown in Figure 5.

3. Results

Numerical simulation was performed for evaluating flow induced mechanical parameters (stress, pressure, shear rate, and velocity) on the inner lumen of the serpentine model. For condition I (slip wall with sinusoidal flow, refer to Table 3), significant variation of pressure, shear rate, velocity, and stress were observed across the entire model, as shown in Figure 6. For velocity models V1 and V3, there was no variation in the pressure parameter over time. Maximum pressure generation for velocity models V1, V2 and V3 were recorded as 2–5 × 105 Pa. Similarly, minimum pressure generation for different velocity models V1, V2, and V3 were obtained as 0.2 × 105, −1 × 105, and −1 × 105 Pa, respectively. In the velocity model V2, pressure generation on point P1 at a time T was maximum, compared to other points, as shown in Figure 6a. It has been observed that the shear rate profile was globally constant for all velocity models, ranging from 0.5 × 104 to 4.5 × 104 1/s. The shear rate results for all velocity models, around points P1 and P3 at time T was maximum, as shown in the Figure 6b. For all velocity models, it can be concluded that the flow development was linearly correlated with the progression of time within the T cycle, as shown in Figure 6c. The minimum and maximum axial velocity of 2 mm/min at neck region (P1) and 20 mm/min in between the abdominal (P2) and rear regions (P3), respectively, were observed for condition I. From Figure 6d, it had been inferred that the local maximum stress having a magnitude of 18 × 10−10 N/m2 had been generated for the case of velocity model V2 at the end of the cycle in between the abdominal (P2) and rear regions (P3).
For the sinusoidal flow with no-slip wall boundary condition (condition II), local minimal variation in pressure, shear rate, and axial velocity were analyzed for velocity models V1, V2, and V3 (Figure 7). A significant variation of pressure distribution was recorded for velocity model V1 compared to the two other velocity models. For velocity models V2 and V3, the local maximum pressure equal to 4 × 105 Pa was observed in the neck region (P1) of the channel at end of cycle time T, as shown in Figure 7a. Maximum pressure was generated at the time T/2 sec in the periphery of the rear region (P2). It was noticed that the value of axial velocity was similar for V1, V2, and V3 velocity models. The maximum velocity was generated at the end of the T cycle for all inlet velocity models. The minimum and maximum axial velocities were observed as 0 and 20 mm/min, respectively, in the neck region (P1) and periphery of the rear region (P3), as shown in Figure 7c. Negligible variation had been observed in the shear rate for different velocity models in all regions of the serpentine at T/3, T/2, and T sec. The shear rate profiles for velocity models V1, V2, and V3 were similar and its range of magnitude was 0.5–4.5 × 105 1/s, as shown in Figure 7b.
The contour plot of pressure for velocity models V4, V5, and V6 at a measured time period T/3, T/2, and T sec was evaluated similarly for physiological flow with no slip wall boundary condition (condition III). For each velocity model, it was observed that local maximum pressure (3.5–4 × 105 Pa) developed at the T/3 and T/2 sec at the neck region of the channel, as shown in Figure 8a. For the same condition, the local maximum (at neck region) and minimum (rear and abdominal regions) shear rate value were predicted as 4.5 × 105 1/s and 0.5 × 105 1/s, respectively (Figure 8b). The range of local minimum and maximum axial velocity was calculated between 0 to 18 mm/min, respectively (Figure 8c). The velocity contour was found to be higher for periods T/3 and T/2 sec. Maximum velocity had been developed at the neck region of the serpentine channel.
The contour plot of physiological flow with slip wall boundary condition (condition IV) depicts that for velocity model V4, the pressure magnitude at T/3 and T/2 sec was minimum (0.2 × 1010 Pa). The maximum pressure (1 × 1010 Pa) was developed at approximately the periphery of the neck region and the rear region at the end of the T cycle, as shown in Figure 9a. At T/3 and T/2 secs, local maximum pressure was developed near the periphery of P1, as shown in Figure 9a. The contour plot of the shear rate was found similar for all velocity models of the physiological flow condition. The minimum shear rate corresponding to V4, V5, and V6 were obtained as 1 × 107, 0.2 × 105, and 0 × 105, respectively. Similarly, the maximum shear rate for V4, V5, and V6 were 8 × 107, 1.8 × 105, and 2.5 × 105, respectively, as shown in Figure 9b. A maximum axial velocity of 200 mm/min was developed on all measured points at T sec, which was extremely high as compared to other conditions shown in Figure 9c. For velocity models V5 and V6, a maximum pressure of 2.5 × 105 Pa had been observed on the entire channel at T/3 and T/2 sec at the inlet region. For the same period, the local maximum velocity for model V6 were found to be higher than that of V5 in the neck as well as the periphery of rear region. For the results obtained in Figure 9d, it was observed that the stress decreases with an increase in velocity. Global values of stress for all velocities fluctuated between 2 and 25 × 10−10 N/m2 in the rear region of the serpentine channel.
After evaluating the magnitude of various mechanical parameters, the printed serpentine channels were imaged every alternate day using fluorescence microscopy to assess cell attachment and cell distribution within the channel. A sample channel image stained with CFSE is shown in Figure 10a, but the CFSE fluorescent dye may still be producing fluorescence even if the cells are dead. So, for longer time periods samples were imaged without CFSE dye. Sample image of the channel and a cross-section of the channel after 14 days is shown in Figure 10b. Analysis of the LIVE-DEAD assay demonstrated that on average 98.5% of the cells were viable after 14 days in the incubator, as shown in Figure 10c, corresponding to the inlet flow boundary condition of 4.62 mm/min for the luminal fluid (or media). Figure 10d, show a live dead assay fluorescence image for an inlet boundary condition of 4.48 mm/min and 4.76 mm/min, respectively. It is evident from Figure 10d, that cell density in the serpentine channel is less than that corresponding to optimized mechanical parameters, as presented in Figure 10c.

4. Discussion

The variation in magnitude of stress, pressure, shear rate, and axial velocity of the fluid flowing through the serpentine vascular channels depends on the material properties, orientation of serpentine geometry, and the flow architecture of working fluid at points P1, P2, and P3. For condition I (as given in Table 3), the proposed architecture follows the Hagen–Poiseuille model for a Newtonian fluid where pressure and velocities are proportional to each other [44]. Variation in the axial velocity influences the shear rate distribution in a closed channel flow. The magnitude of the shear rate describes the flow behavior of working fluid inside the serpentine channel. A minimal variation in shear rate has been found due to a nominal variation in the axial velocity with respect to the radius of the serpentine because of considered rigid wall configuration [45]. Variation in axial velocities also influences the downstream generated pressure in serpentine structures. Localized maximal pressure was generated due to chaotic motion of fluid particles. When the axial velocity of a fluid increases, some of the energy used by the random motion particles to follow the fluid direction develops a lower downstream pressure [46]. The developing downstream pressure further gives rise to localized stress. It was also correlated that the velocity enhancement at the downstream region leads to formation of maximal stress. At higher axial velocity, fluid flow moved slowly near the wall due to diffusion and dispersion of fluid particles [47]. Hence, the inlet velocity V2 was found to create higher stress than its counterpart velocity V3 = Performing transient analysis, the maximum variation in the axial velocity profile was obtained at the end of the full cycle in condition II (see Table 3) due to positive acceleration of fluid particles [48].
For the no-slip wall boundary condition with physiological flow at the inlet of serpentine (as given in Table 3), a small variation compared to the maximum magnitude of the shear rate was obtained near the wall due to rigidity of the wall of the serpentine channel as well as due to negligible variation of axial velocity magnitude [49]. There was minimum fluid velocity near the wall, which produces a non-significant variation of axial velocity near the wall. According to Hagen–Poiseuille model, maximum pressure magnitude was responsible for producing maximum axial velocity in a serpentine vascular model [50]. For a no-slip condition, when fluid comes into contact with the wall of the serpentine channel, there is no relative movement between them, which decreases the stress magnitude [51].
For the slip boundary wall condition with physiological flow at the inlet of the serpentine (as presented in Table 3), the abrupt changes in axial velocity were contoured by the inlet velocity profile (V4), which occurred due to slip conditions at the end of the T cycle. When a fluid flow comes into contact with the curvature section of the serpentine, it creates Coriolis force, which causes a transversal slope in the flowing fluid. Resultant interaction between Coriolis force and transversal slope develops a secondary force on the flow cross-section. This secondary force had disseminated to the curvature section, producing a higher axial velocity magnitude [52]. Maximum axial velocity and presence of curvature of proposed serpentine structure were the effective parameters for developing maximum pressure near the neck and rear region of the serpentine structure at the end of T cycle.
The temporo-spatial deviations in flow-derived parameters were found to affect cell viability and functionality, as given in Table 7. For condition I, the overall maximum and minimum deviation was obtained for pressure and shear rate parameters, respectively. The minimum deviation was obtained due to small changes in the measured point’s value. Serpentine structure and different axial velocity magnitude were responsible for producing maximum pressure deviation [53]. The maximum velocity deviation was obtained due to development of Coriolis force by the serpentine structure, which had caused a radial pressure gradient [54]. Increased pressure gradient radially had further influenced the radial flow of the media, thereby enhancing the deviation in shear stress of adjacent regions [55].
For condition II (refer to Table 3), minimum deviation was obtained for the axial velocity on point P2 at the end of T/3 cycle. Similarly, the minimum deviation was reported for the pressure on point P3 at the end of the T cycle. As the flow approached near the wall, it was converted into a transitional flow due to elasticity of the wall of the serpentine. Hence, the velocity of flowing fluid increases, producing a maximum deviation of velocities [54].
The laminar flow used in the model and the absence of resistivity to the flow of fluid particles were responsible for producing a slight deviation in the velocity profile [56,57]. For condition III, it was observed that the deviation of all parameters had increased with respect to time. The variation in measured parameters had occurred due to the serpentine structure of the model, responsible for inducing the heterogeneity in the nature of the flow [58]. The maximum and minimum deviation was obtained at point P3 in boundary condition IV (refer to Table 3) due to the physiological relevance flow in the presence of the elastic wall. Such a flow condition brought a non-stationary axial velocity profile over the period of time [59].
The numerical analysis aided in the selection of optimized values of different biofluid dynamics features, such as axial velocity, shear rate, pressure, and stress of the proposed serpentine blood vessel model for producing an ideal condition for cell viability. High cellular viability was observed even after 14 days of seeding of HUVECs in the serpentine channel maintained in the bioreactor. Cell attachment was almost uniform across the channel as can be observed from the fluorescent microscopy and confocal images.

5. Conclusions

A 3D bioprinting system using an extrusion-based technique was developed to fabricate vascular serpentine channels. The optimized perfusion rate of the media through the printed serpentine channels provided a biomimetic micro-niche for proliferation of seeded HUVECs. Dynamic cultured media at the extracellular surface of seeded cells imposed external forces on the membrane in the form of shear stress and its time derivatives on serpentine walls. Such a distributed pattern of extracellular forces induced faster and synchronous mechanotransduction to the focal points of adhesion of these cells. Hence, a variation in cellular proliferative activities was observed in response to the different perfusion rate of media. Therefore, it can be concluded that flow dynamics features, such as axial velocity, shear rate, stress, and pressure played a vital role in reproducing the physiological behavior and protruding adhesion properties of cells under in vitro conditions on the surface of the serpentine plane. Longitudinal stress generated by the flowing fluid also affects the HUVECS cell’s functionality. It was also observed that the modeling of the serpentine with an elastic wall and consideration of physiological flow conditions at the inlet of the model brought maximal dynamicity in the axial velocity. On the other hand, identified local maxima and minima of flow parameters at designated regions of the serpentine were found to be a useful tool for effective seeding of cells for their proliferation over the surface of the serpentine. The minimum deviation of pressure, shear rate, and velocity were obtained for the sinusoidal flow with a slip wall condition. Result of Live-Dead assay showed the viability of HUVECs after 14 days of printing the channels, demonstrating that the obtained value of fluid dynamics parameters in the serpentine channel was helpful in maintaining the cell proliferation inside the bioreactor.

Author Contributions

Conceptualization, A.B. and K.M.; Research Execution, K.D.; Funding, S.G., A.B. and K.M.; Result Analysis, K.D., A.B., S.G. and K.M.; Manuscript preparation, K.D. and A.B.; Manuscript Review, A.B. and K.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received funding from NASA Florida Space Grant Program and “SPARC [SPARC/2018-2019/140/SL (IN)]”.

Institutional Review Board Statement

The study was approved by Institute Ethical Committee (NITRR/IEC/2021/12).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be made available based on request made.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mount, C.; Dusserrre, N.; McAllister, T.; L’Heureux, N. Tissue-engineered cardiovascular grafts and novel applications of tissue engineering by self-assembly (TESATM). In Cardiac Regeneration and Repair: Biomaterials and Tissue Engineering; Woodhead Publishing Ltd.: Cambridge, UK, 2014; pp. 410–451. [Google Scholar] [CrossRef]
  2. Ng, H.Y.; Lee, K.-X.A.; Kuo, C.-N.; Shen, Y.-F. Bioprinting of artificial blood vessels. Int. J. Bioprinting 2018, 4, 140. [Google Scholar] [CrossRef] [PubMed]
  3. Fazal, F.; Raghav, S.; Callanan, A.; Koutsos, V.; Radacsi, N. Recent advancements in the bioprinting of vascular grafts. Biofabrication 2021, 13, 032003. [Google Scholar] [CrossRef] [PubMed]
  4. Jin, Q.; Fu, Y.; Zhang, G.; Xu, L.; Jin, G.; Tang, L.; Ju, J.; Zhao, W.; Hou, R. Nanofiber electrospinning combined with rotary bioprinting for fabricating small-diameter vessels with endothelium and smooth muscle. Compos. Part B Eng. 2022, 234, 109691. [Google Scholar] [CrossRef]
  5. Forget, A.; Derme, T.; Mitterberger, D.; Heiny, M.; Sweeney, C.; Mudili, L.; Dargaville, T.R.; Shastri, V.P. Architecture-inspired paradigm for 3D bioprinting of vessel-like structures using extrudable carboxylated agarose hydrogels. Emergent Mater. 2019, 2, 233–243. [Google Scholar] [CrossRef]
  6. Liu, H.; Zhou, H.; Lan, H.; Liu, T.; Liu, X.; Yu, H. 3D printing of artificial blood vessel: Study on multi-parameter optimization design for vascular molding effect in alginate and gelatin. Micromachines 2017, 8, 237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Nagy, J.A.; Chang, S.-H.; Dvorak, A.M.; Dvorak, H.F. Why are tumour blood vessels abnormal and why is it important to know? Br. J. Cancer 2009, 100, 865–869. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Schuurman, W.; Levett, P.A.; Pot, M.W.; van Weeren, P.R.; Dhert, W.J.A.; Hutmacher, D.W.; Melchels, F.P.W.; Klein, T.J.; Malda, J. Gelatin-methacrylamide hydrogels as potential biomaterials for fabrication of tissue-engineered cartilage constructs. Macromol. Biosci. 2013, 13, 551–561. [Google Scholar] [CrossRef]
  9. Askari, M.; Naniz, M.A.; Kouhi, M.; Saberi, A.; Zolfagharian, A.; Bodaghi, M. Recent progress in extrusion 3D bioprinting of hydrogel biomaterials for tissue regeneration: A comprehensive review with focus on advanced fabrication techniques. Biomater. Sci. 2021, 9, 535–573. [Google Scholar] [CrossRef]
  10. Kim, S.H.; Seo, Y.B.; Yeon, Y.K.; Lee, Y.J.; Park, H.S.; Sultan, T.; Lee, J.M.; Lee, J.S.; Lee, O.J.; Hong, H.; et al. 4D-bioprinted silk hydrogels for tissue engineering. Biomaterials 2020, 260, 120281. [Google Scholar] [CrossRef]
  11. Sasmal, P.; Datta, P.; Wu, Y.; Ozbolat, I.T. 3D bioprinting for modelling vasculature. Microphysiol. Syst. 2018, 2, 9. [Google Scholar] [CrossRef]
  12. Wang, X.; Liu, C. Fibrin hydrogels for endothelialized liver tissue engineering with a predesigned vascular network. Polymers 2018, 10, 1048. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Ahmed, E.M. Hydrogel: Preparation, characterization, and applications: A review. J. Adv. Res. 2015, 6, 105–121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Gasperini, L.; Mano, J.F.; Reis, R.L. Natural polymers for the microencapsulation of cells. J. R. Soc. Interface 2014, 11, 20140817. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Sun, W.; Gregory, D.A.; Tomeh, M.A.; Zhao, X. Silk fibroin as a functional biomaterial for tissue engineering. Int. J. Mol. Sci. 2021, 22, 1499. [Google Scholar] [CrossRef]
  16. Azadani, A.N. In vitro experimental methods for assessment of prosthetic heart valves. In Principles of Heart Valve Engineering; Academic Press: Cambridge, MA, USA, 2019. [Google Scholar] [CrossRef]
  17. Neofytou, P.; Drikakis, D. Effects of blood models on flows through a stenosis. Int. J. Numer. Methods Fluids 2003, 43, 597–635. [Google Scholar] [CrossRef]
  18. Morbiducci, U.; Ponzini, R.; Rizzo, G.; Cadioli, M.; Esposito, A.; Montevecchi, F.M.; Redaelli, A. Mechanistic insight into the physiological relevance of helical blood flow in the human aorta: An in vivo study. Biomech. Model. Mechanobiol. 2011, 10, 339–355. [Google Scholar] [CrossRef]
  19. Wang, Y.; Nunna, B.B.; Talukder, N.; Lee, E.S. Microfluidic-Based Novel Optical Quantification of Red Blood Cell Concentration in Blood Flow. Bioengineering 2022, 9, 247. [Google Scholar] [CrossRef]
  20. Stalmans, I.; Vandewalle, E.; Anderson, D.R.; Costa, V.P.; Frenkel, R.E.P.; Garhofer, G.; Grunwald, J.; Gugleta, K.; Harris, A.; Hudson, C.; et al. Use of colour Doppler imaging in ocular blood flow research. Acta Ophthalmol. 2011, 89, e609–e630. [Google Scholar] [CrossRef]
  21. Mittal, N.; Zhou, Y.; Linares, C.; Ung, S.; Kaimovitz, B.; Molloi, S.; Kassab, G.S. Analysis of blood flow in the entire coronary arterial tree. Am. J. Physiol.-Heart Circ. Physiol. 2005, 289, H439–H446. [Google Scholar] [CrossRef]
  22. Miao, P.; Rege, A.; Li, N.; Thakor, N.V.; Tong, S. High resolution cerebral blood flow imaging by registered laser speckle contrast analysis. IEEE Trans. Biomed. Eng. 2010, 57, 1152–1157. [Google Scholar] [CrossRef]
  23. Mandal, S.; Mukhopadhyayy, S.; Layek, G. Pulsatile flow of shear-dependent fluid in a stenosed artery. Theor. Appl. Mech. 2012, 39, 209–231. [Google Scholar] [CrossRef]
  24. Chen, H.; Yu, Z.; Bai, S.; Lu, H.; Xu, D.; Chen, C.; Liu, D.; Zhu, Y. Microfluidic models of physiological or pathological flow shear stress for cell biology, disease modeling and drug development. TrAC Trends Anal. Chem. 2019, 117, 186–199. [Google Scholar] [CrossRef]
  25. Ali, N.; Sajid, M.; Javed, T.; Abbas, Z. Heat transfer analysis of peristaltic flow in a curved channel. Int. J. Heat Mass Transf. 2010, 53, 3319–3325. [Google Scholar] [CrossRef]
  26. Lamberti, G.; Soroush, F.; Smith, A.; Kiani, M.F.; Prabhakarpandian, B.; Pant, K. Adhesion patterns in the microvasculature are dependent on bifurcation angle. Microvasc. Res. 2015, 99, 19–25. [Google Scholar] [CrossRef] [Green Version]
  27. Charonko, J.; Karri, S.; Schmieg, J.; Prabhu, S.; Vlachos, P. In vitro, time-resolved PIV comparison of the effect of stent design on wall shear stress. Ann. Biomed. Eng. 2009, 37, 1310–1321. [Google Scholar] [CrossRef]
  28. Zorlutuna, P.; Rong, Z.; Vadgama, P.; Hasirci, V. Influence of nanopatterns on endothelial cell adhesion: Enhanced cell retention under shear stress. Acta Biomater. 2009, 5, 2451–2459. [Google Scholar] [CrossRef]
  29. Buchanan, C.F.; Verbridge, S.S.; Vlachos, P.P.; Rylander, M.N. Flow shear stress regulates endothelial barrier function and expression of angiogenic factors in a 3D microfluidic tumor vascular model. Cell Adhes. Migr. 2014, 8, 517–524. [Google Scholar] [CrossRef] [Green Version]
  30. Metaxa, E.; Meng, H.; Kaluvala, S.R.; Szymanski, M.P.; Paluch, R.A.; Kolega, J. Nitric oxide-dependent stimulation of endothelial cell proliferation by sustained high flow. Am. J. Physiol. Heart Circ. Physiol. 2008, 295, H736–H742. [Google Scholar] [CrossRef]
  31. Chakraborty, A.; Chakraborty, S.; Jala, V.R.; Haribabu, B.; Sharp, M.K.; Berson, R.E. Effects of biaxial oscillatory shear stress on endothelial cell proliferation and morphology. Biotechnol. Bioeng. 2012, 109, 695–707. [Google Scholar] [CrossRef]
  32. Dick, M.; MacDonald, K.; Tardif, J.-C.; Leask, R.L. The effect of simvastatin treatment on endothelial cell response to shear stress and tumor necrosis factor alpha stimulation. Biomed. Eng. Online 2015, 14, 58. [Google Scholar] [CrossRef] [Green Version]
  33. Franzoni, M.; Cattaneo, I.; Ene-Iordache, B.; Oldani, A.; Righettini, P.; Remuzzi, A. Design of a cone-and-plate device for controlled realistic shear stress stimulation on endothelial cell monolayers. Cytotechnology 2016, 68, 1885–1896. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Dolan, J.M.; Sim, F.J.; Meng, H.; Kolega, J. Endothelial cells express a unique transcriptional profile under very high wall shear stress known to induce expansive arterial remodeling. Am. J. Physiol. Physiol. 2012, 302, C1109–C1118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Pepelanova, I.; Kruppa, K.; Scheper, T.; Lavrentieva, A. Gelatin-methacryloyl (GelMA) hydrogels with defined degree of functionalization as a versatile toolkit for 3D cell culture and extrusion bioprinting. Bioengineering 2018, 5, 55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Karanasiou, G.S.; Gatsios, D.A.; Lykissas, M.G.; Stefanou, K.A.; Rigas, G.A.; Lagaris, I.E.; Kostas-Agnantis, I.P.; Gkiatas, I.; Beris, A.E.; Fotiadis, D.I. Modeling of blood flow through sutured micro-vascular anastomoses. In Proceeding of the 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Milano, Italy, 25–29 August 2015; pp. 1877–1880. [Google Scholar]
  37. Nilsson, A.; Zhang, Q.; Styf, J. Using the amplitude of pulse-synchronous intramuscular pressure oscillations when diagnosing chronic anterior compartment syndrome. Orthop. J. Sport. Med. 2014, 2, 2325967114556443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Leone, M.; Asfar, P.; Radermacher, P.; Vincent, J.L.; Martin, C. Optimizing mean arterial pressure in septic shock: A critical reappraisal of the literature. Crit. Care 2015, 19, 1–7. [Google Scholar] [CrossRef] [Green Version]
  39. Kundu, R.N.; Biswas, S.; Das, M. Mean Arterial Pressure Classification: A Better Tool for Statistical Interpretation of Blood Pressure Related Risk Covariates. Cardiol. Angiol. Int. J. 2017, 6, 1–7. [Google Scholar] [CrossRef]
  40. Papaioannou, T.G.; Protogerou, A.D.; Vrachatis, D.; Konstantonis, G.; Aissopou, E.; Argyris, A.; Nasothimiou, E.; Gialafos, E.J.; Karamanou, M.; Tousoulis, D.; et al. Mean arterial pressure values calculated using seven different methods and their associations with target organ deterioration in a single-center study of 1878 individuals. Hypertens. Res. 2016, 39, 640–647. [Google Scholar] [CrossRef] [PubMed]
  41. Bit, A.; Chattopadhyay, H. Numerical investigations of pulsatile flow in stenosed artery. Acta Bioeng. Biomech. 2014, 16, 33–44. [Google Scholar] [CrossRef]
  42. Roache, P.J. A method for uniform reporting of grid refinement studies. ASME 1993, 158, 109. [Google Scholar] [CrossRef]
  43. Brown, A.L.; Glaze, D.J.; Pierce, F. Sensitivity analysis and verification of a 1-D surface solid combustion model for a fire CFD boundary condition. In Proceeding of the 11th AIAA/ASME Joint Thermophysics and Heat Transfer Conference, Atlanta, GA, USA, 16–20 June 2014; p. 2256. [Google Scholar]
  44. Debnath, S.; Saha, A.K.; Mazumder, B.; Roy, A.K. Hydrodynamic dispersion of reactive solute in a Hagen–Poiseuille flow of a layered liquid. Chin. J. Chem. Eng. 2017, 25, 862–873. [Google Scholar] [CrossRef]
  45. Abed, W.M.; Whalley, R.D.; Dennis, D.J.; Poole, R.J. Experimental investigation of the impact of elastic turbulence on heat transfer in a serpentine channel. J. Non-Newtonian Fluid Mech. 2016, 231, 68–78. [Google Scholar] [CrossRef] [Green Version]
  46. Pecanha, R.P. Fluid Particles: A Review. J. Chem. Eng. Process Technol. 2015, 6, 2–6. [Google Scholar] [CrossRef] [Green Version]
  47. Sun, Q.; Wu, Y.; Liu, L.; Wiwatanapataphee, B. Study of a newtonian fluid through circular channels with slip boundary taking into account electrokinetic effect. Abstr. Appl. Anal. 2013, 2013, 718603. [Google Scholar] [CrossRef] [Green Version]
  48. Jha, B.K.; Malgwi, P.B. Computational analysis on unsteady hydromagnetic Couette flow of fluid—Particle suspension in an accelerated porous channel. Partial Differ. Equ. Appl. Math. 2022, 5, 100370. [Google Scholar] [CrossRef]
  49. Raj, A.; Rajak, D.K.; Gautam, S.; Guria, C.; Pathak, A.K. Shear rate estimation: A detailed review. In Proceedings of the Annual Offshore Technology Conference, Houston, TX, USA, 2–5 May 2016. [Google Scholar] [CrossRef]
  50. Qin, R.; Duan, C. The principle and applications of Bernoulli equation. J. Phys. Conf. Ser. 2017, 916, 012038. [Google Scholar] [CrossRef]
  51. Rao, I.J.; Rajagopal, K.R. Effect of the slip boundary condition on the flow of fluids in a channel. Acta Mech. 1999, 135, 113–126. [Google Scholar] [CrossRef]
  52. Pradhan, A.; Khatua, K.K.; Sankalp, S. Variation of Velocity Distribution in Rough Meandering Channels. Adv. Civ. Eng. 2018, 2018, 1569271. [Google Scholar] [CrossRef]
  53. Dutta, P.; Nandi, N. Effect of bend curvature on velocity & pressure distribution from straight to a 90° pipe bend-A Numerical Study. REST J. Emerg. Trends Model. Manuf. 2016, 2, 103–108. [Google Scholar]
  54. Amaratunga, M.; Rabenjafimanantsoa, H.A.; Time, R.W. Estimation of shear rate change in vertically oscillating non-Newtonian fluids: Predictions on particle settling. J. Non-Newton. Fluid Mech. 2020, 277, 104236. [Google Scholar] [CrossRef]
  55. Swearingen, C. >Velocity-Profile Deviations Influence Flowmeter Performance. Cole-Parmer. 2018. Available online: https://www.coleparmer.com/tech-article/selecting-the-right-flowmeter-for-the-liquid (accessed on 30 November 2021).
  56. Rogers, D.F. Laminar Flow Analysis; Cambridge University Press: Cambridge, UK, 1992; ISBN 0-521-41152-1. [Google Scholar]
  57. Nave, R. “Laminar Flow”. In HyperPhysics; Georgia State University: Atlanta, GA, USA, 2005. [Google Scholar]
  58. Ansari, S.; Rashid, A.I.; Waghmare, P.R.; Nobes, D.S. Measurement of the flow behavior index of Newtonian and shear-thinning fluids via analysis of the flow velocity characteristics in a mini-channel. SN Appl. Sci. 2020, 2, 1–15. [Google Scholar] [CrossRef]
  59. Wolański, W.; Gzik-Zroska, B.; Joszko, K.; Kawlewska, E.; Sobkowiak, M.; Gzik, M.; Kaspera, W. Impact of vessel mechanical properties on hemodynamic parameters of blood flow. Adv. Intell. Syst. Comput. 2017, 623, 271–278. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the proposed serpentine model.
Figure 1. Schematic diagram of the proposed serpentine model.
Micromachines 13 01766 g001
Figure 2. Comparison of velocity variables and extrapolated values between two grid solutions and the Richardson extrapolation estimation.
Figure 2. Comparison of velocity variables and extrapolated values between two grid solutions and the Richardson extrapolation estimation.
Micromachines 13 01766 g002
Figure 3. Schematic of 3D printer system for printing of serpentine channels.
Figure 3. Schematic of 3D printer system for printing of serpentine channels.
Micromachines 13 01766 g003
Figure 4. Experimental set up details (a) syringe pump setup, (b) magnifying image of microfluidic chip.
Figure 4. Experimental set up details (a) syringe pump setup, (b) magnifying image of microfluidic chip.
Micromachines 13 01766 g004
Figure 5. Sensitivity analysis on pressure parameters for two different values of densities.
Figure 5. Sensitivity analysis on pressure parameters for two different values of densities.
Micromachines 13 01766 g005
Figure 6. Contour Plot of (a) pressure, (b) shear rate, (c) axial velocity, (d) wall stress for sinusoidal flow, slip wall boundary condition for velocity models V1, V2 and V3.
Figure 6. Contour Plot of (a) pressure, (b) shear rate, (c) axial velocity, (d) wall stress for sinusoidal flow, slip wall boundary condition for velocity models V1, V2 and V3.
Micromachines 13 01766 g006
Figure 7. Contour Plot of (a) pressure, (b) shear rate, (c) axial velocity for sinusoidal flow, no-slip wall boundary condition for velocity models V1, V2, and V3.
Figure 7. Contour Plot of (a) pressure, (b) shear rate, (c) axial velocity for sinusoidal flow, no-slip wall boundary condition for velocity models V1, V2, and V3.
Micromachines 13 01766 g007
Figure 8. Contour Plot of (a) pressure, (b) shear rate, (c) axial velocity for physiological flow, no-slip wall boundary condition for velocity models V4, V5 and V6.
Figure 8. Contour Plot of (a) pressure, (b) shear rate, (c) axial velocity for physiological flow, no-slip wall boundary condition for velocity models V4, V5 and V6.
Micromachines 13 01766 g008
Figure 9. Contour Plot of (a) pressure, (b) shear rate, (c) axial velocity, (d) wall stress for physiological flow, slip wall boundary condition for velocity models V4, V5 and V6.
Figure 9. Contour Plot of (a) pressure, (b) shear rate, (c) axial velocity, (d) wall stress for physiological flow, slip wall boundary condition for velocity models V4, V5 and V6.
Micromachines 13 01766 g009
Figure 10. (a) Sample bioprinted tissue construct with embedded serpentine channel, (b) fluorescent imaging after 14 days, (c) confocal imaging after 14 days by staining with LIVE-DEAD assay of serpentine channel seeded with HUVECs at inlet flow velocity = 4.62 mm/min, (d) 4.48 mm/min, (e) 4.76 mm/min.υ.
Figure 10. (a) Sample bioprinted tissue construct with embedded serpentine channel, (b) fluorescent imaging after 14 days, (c) confocal imaging after 14 days by staining with LIVE-DEAD assay of serpentine channel seeded with HUVECs at inlet flow velocity = 4.62 mm/min, (d) 4.48 mm/min, (e) 4.76 mm/min.υ.
Micromachines 13 01766 g010aMicromachines 13 01766 g010b
Table 1. Mesh statistics.
Table 1. Mesh statistics.
ParametersSlipNo Slip
Mesh vertices22201562
Element typeAll elements
Triangles26081578
Quads712608
Edge elements762328
Vertex elements4038
Number of elements33202186
Minimum element quality0.27580.3608
Average element quality0.76630.7733
Element area ratio0.07610.04016
Mesh area20.31 mm221.33 mm2
Table 2. Velocity function for V1 to V6.
Table 2. Velocity function for V1 to V6.
Inlet VelocityModelParameters ValueFunction
Sinusoidal Flow (SF)V1f = 1.25 Hz, t = 0 to 1 sV1(t) = 4.48 × sin 2 × π × f × t
V2V2(t) = 4.62 × sin 2 × π × f × t
V3V3(t) = 4.76 × sin 2 × π × f × t
Physiological Flow (PF)
Micromachines 13 01766 i001
V4t = 0 to 1 s, a1 = 3.04, b1 = 0.8634, c1 = −0.07247, a2 = 1.24, b2 = 14.02, c2 = 2.164, a3 = 2.619, b3 = 5.81, c3 = −0.5393, a4 = 0.4497, b4 = 24.25, c4 = 1.474, a5 = 0.3847, b5 = 32.5, c5 = −2.181, a6 = 0.4739, b6 = 20.77, c6 = 0.4904, a7 = 0.2708, b7 = 37.65, c7 = 2.495, a8 = 0.1386, b8 = 52.47, c8 = −0.7126 V4(t) = a 1 × s i n b 1 × t + c 1 + a 2 × s i n b 2 × t + c 2 + a 3 × s i n b 3 × t + c 3 + a 4 × s i n b 4 × t + c 4 + a 5 × s i n b 5 × t + c 5 + a 6 × s i n b 6 × t + c 6 + a 7 × s i n b 7 × t + c 7 + a 8 × s i n b 8 × t + c 8
V5t = 0 to 1 s, a1 = 2.721, b1 = 1.588, c1 = −0.4734, a2 = 1.364, b2 = 13.48, c2 = 2.332, a3 = 3.081, b3 = 5.052, c3 = −0.1051, a4 = 0.3837, b4 = 24.37, c4 = 1.263, a5 = 0.4042, b5 = 32.44, c5 = −2.133, a6 = 0.4619, b6 = 19.567, c6 = 0.9404, a7 = 0.2744, b7 = 37.67, c7 = 2.482, a8 = 0.1435, b8 = 52.48, c8 = −0.7228V5(t) = a 1 × s i n b 1 × t + c 1 + a 2 × s i n b 2 × t + c 2 + a 3 × s i n b 3 × t + c 3 + a 4 × s i n b 4 × t + c 4 + a 5 × s i n b 5 × t + c 5 + a 6 × s i n b 6 × t + c 6 + a 7 × s i n b 7 × t + c 7 + a 8 × s i n b 8 × t + c 8
V6t = 0 to 1 s, a1 = 2.864, b1 = 1.405, c1 = −0.3165, a2 = 1.322, b2 = 13.96, c2 = 2.122, a3 = 3.008, b3 = 5.506, c3 = −0.3209, a4 = 0.7938, b4 = 22.63, c4 = 2.259, a5 = 0.4209, b5 = 31.95, c5 = −1.982, a6 = 0.8909, b6 = 20.65, c6 = 0.3462, a7 = 0.3141, b7 = 37.4, c7 = 2.577, a8 = 0.1474, b8 = 52.51, c8 = −0.7326V6(t) = a 1 × s i n b 1 × t + c 1 + a 2 × s i n b 2 × t + c 2 + a 3 × s i n b 3 × t + c 3 + a 4 × s i n b 4 × t + c 4 + a 5 × s i n b 5 × t + c 5 + a 6 × s i n b 6 × t + c 6 + a 7 × s i n b 7 × t + c 7 + a 8 × s i n b 8 × t + c 8
Table 3. Different boundary conditions for numerical analysis.
Table 3. Different boundary conditions for numerical analysis.
ConditionWall Boundary ConditionInlet Boundary ConditionInlet Velocity
ISlip (S)Sinusoidal Flow (SF)V1
V2
V3
IINo-Slip (NS)Sinusoidal Flow (SF)V1
V2
V3
IIINo-Slip (NS)Physiological Flow (PF)V4
V5
V6
IVSlip (S)Physiological Flow (PF)V4
V5
V6
Table 4. Generated data for pressure, shear rate, and velocity from simulation.
Table 4. Generated data for pressure, shear rate, and velocity from simulation.
ProbeGridNumber of ElementsPressure at Time TThe Shear Rate at Time TVelocity at Time T
P1Fine2366210,698.842,091.8112.64
Normal1578215,248.6043,115.8212.71
Coarse1106230,802.244,968.6312.99
P2Fine2366125,274.6910,222.5312.66
Normal1578123,677.9310,938.7212.25
Coarse1106136,541.8416,772.0712.09
Table 5. Grid convergence index for different variables.
Table 5. Grid convergence index for different variables.
VariablerpFs G C I 23 ( % ) G C I 12 ( % )
Pressure (P1)1.2096.481.253.721.11
Shear Rate (P1)3.126.63.7
Velocity (P1)6.411.130.3
Pressure (P2)1022.90.28
Shear Rate (P2)9.3913.481.77
Velocity (P2)270.0030.24
Table 6. Parameters used for 3D printing of serpentine channel.
Table 6. Parameters used for 3D printing of serpentine channel.
ParametersSpecification
Printing house temperature37 °C
Extrusion pressure of the coaxial nozzle0.1 to 1 Mpa
Nozzle diameter1 mm
Needle size20 gauge
Extrusion rate5 mm/s
Dispensing speed4 mm/s
Printing bed temperature and material2 °C to 40 °C and Aluminum/Glass
Pulsatile pump speed, pressure and flow rate5 to 36 rpm, 80 to 120 mmHg; and 1–4 mL/min
pH of media7.2 to 7.4
Table 7. Calculated deviation for pressure, shear rate and velocity.
Table 7. Calculated deviation for pressure, shear rate and velocity.
ConfigurationParameterPointTimeDeviation (%)
S_SF_V1, V2, V3PressureP3 (V2) percent [V2, V3]0.5Min0.24937
P2 (V2) percent [V2, V3]0.3Max231.1239
Shear rateP2 (V2) percent [V2, V3]0.5Min0.19267
P2 (V1) percent [V2, V1]0.3Max19.17148
VelocityP2 (V3) percent [V2, V3]0.5Min0.324655
P1(V3) percent [V1, V2]0.5Max36.753
NS_SF_ V1, V2, V3PressureP1 (V3) percent [V2, V3]0.3Min3.603224
P3 (V1) percent [V3, V1]1Max14.14604
Shear rateP1 (V3) percent [V2, V3]0.3Min1.854657
P2 (V1) percent [V3, V1]0.5Max9.254331
VelocityP2 (V3) percent [V2, V3]0.3Min1.775892
P1 (V1) percent [V3, V1]1Max6.57622
NS_PF_V4, V5, V6PressureP1 (V6) percent [V6, V5]1Min3.075083
P3 (V4) percent [V6, V4]0.5Max14.18746
Shear rateP2 (V5) percent [V4, V5]1Min2.104332
P2 (V4) percent [V6, V4]0.5Max8.072614
VelocityP3 (V6) percent [V6, V4]1Min2.813501
P1 (V4) percent [V6, V4]0.5Max6.5163
S_PF_ V4, V5, V6PressureP3 (V6) percent [V6, V5]0.3Min2.485286
P3 (V6) percent [V6, V4]1Max−3.9 × 108
Shear rateP3 (V6) percent [V6, V4]1Min2.959357
P1 (V6) percent [V6, V5]0.3Max3229923
VelocityP3 (V4) percent [V4, V5]0.5Min0.7726
P3 (V5) percent [V4, V5]1Max291,407
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Deshmukh, K.; Gupta, S.; Mitra, K.; Bit, A. Numerical and Experimental Analysis of Shear Stress Influence on Cellular Viability in Serpentine Vascular Channels. Micromachines 2022, 13, 1766. https://doi.org/10.3390/mi13101766

AMA Style

Deshmukh K, Gupta S, Mitra K, Bit A. Numerical and Experimental Analysis of Shear Stress Influence on Cellular Viability in Serpentine Vascular Channels. Micromachines. 2022; 13(10):1766. https://doi.org/10.3390/mi13101766

Chicago/Turabian Style

Deshmukh, Khemraj, Saurabh Gupta, Kunal Mitra, and Arindam Bit. 2022. "Numerical and Experimental Analysis of Shear Stress Influence on Cellular Viability in Serpentine Vascular Channels" Micromachines 13, no. 10: 1766. https://doi.org/10.3390/mi13101766

APA Style

Deshmukh, K., Gupta, S., Mitra, K., & Bit, A. (2022). Numerical and Experimental Analysis of Shear Stress Influence on Cellular Viability in Serpentine Vascular Channels. Micromachines, 13(10), 1766. https://doi.org/10.3390/mi13101766

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop