{'fraction_non_alphanumeric': 0.04340050897585804, 'fraction_numerical': 0.021167205447417292, 'mean_word_length': 4.456140350877193, 'pattern_counts': {'":': 0, '<': 23, '<?xml version=': 0, '>': 21, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 4, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Blood Viscoelasticity Measurement Using Interface Variations in Coflowing Streams under Pulsatile Blood Flows

Blood flows in microcirculation are determined by the mechanical properties of blood samples, which have been used to screen the status or progress of diseases. To achieve this, it is necessary to measure the viscoelasticity of blood samples under a pulsatile blood condition. In this study, viscoelasticity measurement is demonstrated by quantifying interface variations in coflowing streams. To demonstrate the present method, a T-shaped microfluidic device is designed to have two inlets (a, b), one outlet (a), two guiding channels (blood sample channel, reference fluid channel), and one coflowing channel. Two syringe pumps are employed to infuse a blood sample at a sinusoidal flow rate. The reference fluid is supplied at a constant flow rate. Using a discrete fluidic circuit model, a first-order linear differential equation for the interface is derived by including two approximate factors (F1 = 1.094, F2 = 1.1087). The viscosity and compliance are derived analytically as viscoelasticity. The experimental results showed that compliance is influenced substantially by the period. The hematocrit and diluent contributed to the varying viscosity and compliance. The viscoelasticity varied substantially for red blood cells fixed with higher concentrations of glutaraldehyde solution. The experimental results showed that the present method has the ability to monitor the viscoelasticity of blood samples under a sinusoidal flow-rate pattern.


Introduction
Cardiovascular diseases (CVDs) occur without any symptoms and can lead to unexpected death [1]. In other words, blood clotting or abnormal blood flow contributes to vasculature blockages. Currently, biochemical properties (i.e., biomarkers [2,3] or DNA [4]) are used to diagnose CVDs. However, the biochemical approach has not been considered an effective tool for early detection of CVDs, because it does not provide information on blood flows or blood clotting. Instead of the biochemical approach, a biophysical approach is required to quantify abnormal blood flows in narrow-sized vessels. Blood samples collected from patients with CVDs or disorders exhibit changes in cells (i.e., red blood cells (RBCs) [5] or platelets [6]) or plasma proteins [3,7]. To detect CVDs effectively, it is necessary to quantify the contributions of cells or plasma proteins. Blood flows in microcirculation are determined by the mechanical properties of blood samples. In addition, blood vessel walls (i.e., shape and size) contribute to varying blood flows substantially. These properties include viscosity, elasticity, RBC aggregation, and RBC deformability. According to recent reports, a blood sample collected from a patient with CVD showed significantly different biophysical properties when compared with a normal blood sample [2,8,9]. For this reason, the mechanical properties of blood samples have been used to monitor the status or progress of CVDs. Additionally, the mechanical properties have been quantified under dynamic blood flows [10]. Under ex vivo closed-circuit conditions, blood flows vary periodically over time. When a flow regulator is integrated into the fluidic circuit, a periodic flow pattern is regulated to a constant flow pattern [11]. The viscosity of the blood sample is then quantified by using reverse flow-switching phenomena under a constant blood flow. Additionally, blood viscosity is obtained by monitoring the interface in coflowing streams, while the flow rates of the blood sample and reference fluid remain constant at the same flow rate with two syringe pumps [12,13]. However, because the blood sample includes viscoelasticity (viscosity and elasticity) under a periodic flow condition, it is necessary to quantify the viscoelasticity of the blood sample without a flow regulator (i.e., periodic flow rate).
Recently, a microfluidic platform has been suggested for effectively manipulating a small volume of blood sample in microfluidic channels. A microfluidic channel has been used to investigate hemorheological properties of blood samples. Under a microfluidic platform, several techniques have been suggested to quantify the viscoelasticity of blood samples. Guido et al. have measured membrane viscoelasticity by measuring the velocity and shape of a single RBC in converging or diverging flow [14]. Kim et al. measured RBC stretching with viscoelastic cell focusing [15]. The method was then used to characterize differences in RBC deformability. Lee et al. measured the viscosity and elasticity of blood samples by infusing steady and transient blood flows sequentially [16]. The viscosity and time constant were obtained sequentially by controlling the flow rate of the reference fluid and blood sample, respectively. Here, the time constant was quantified by monitoring temporal variation of a bridge channel filled with a human blood sample at a transient blood flow. The elasticity of the blood sample was quantified by using the linear Maxwell model (elasticity = viscosity/time constant). Kang reported that viscosity and elasticity of blood samples can be obtained sequentially under a periodic on-off blood flow condition [17]. Monitoring the interface in a coflowing stream enabled blood viscoelasticity to be obtained at an interval of a specific period.
More recently, the author suggested a viscoelasticity measurement method under a sinusoidal blood flow rate (Q B (t)) (Q B (t) = Q α + Q β ·sin (ωt), where Q α is the mean flow rate, Q β is the alternating flow rate, and ω is the angular frequency of syringe pump) [18]. While the blood sample is infused periodically, the viscoelasticity of the blood sample is obtained by monitoring the interface in coflowing streams and by calculating the pulsatility index (PI) (PI = 0.5 × (Q max -Q min )/Q ave , where Q max is the maximum flow rate, Q min is the minimum flow rate, and Q ave is the average flow rate). A first-order differential equation for coflowing streams is derived by constructing a simple discrete fluidic circuit. Here, using a conventional microelectromechanical system technique, a microfluidic device has rectangular shape with an extremely low aspect ratio (aspect ratio = depth/width = 4/250) which is devised to compensate for the boundary condition difference between a real physical model and a mathematical condition. First, the viscosity of the blood sample is calculated by averaging the equation over a single period. Second, instead of an analytical solution of the equation, the time constant is obtained from the expression of PI. Variations of velocity and interface are required simultaneously to find out time constant of interface. The elasticity is then obtained with a linear Maxwell model. However, the differential equation does not include a correction factor (CF), which makes it necessary to compensate for the boundary condition difference between the real model and mathematical model under coflowing streams [19]. Additionally, the previous study did not use the analytical expression to obtain the viscoelasticity of a blood sample. At last, the method required variations of blood velocity over time. Nonetheless, the previous method shows promise for quantifying the viscoelasticity of blood samples under a periodic blood flow-rate pattern. It is extremely difficult to obtain the viscoelasticity of a blood sample circulated under an ex vivo or in vivo condition (i.e., a real and complex situation). As a preliminary study, it is necessary to develop a simple method for measuring viscoelasticity under a single sinusoidal flow patterns with a syringe pump.
In this study, to resolve these issues, a CF is inserted while deriving a differential equation for coflowing streams. Because the differential equation includes nonlinear terms, it is difficult to find an analytical solution. Conducting computational fluid dynamics (CFD) simulation enables an approximate expression of the CF to be obtained. Then, two approximation factors (F 1 , F 2 ) are suggested and calculated to convert nonlinear terms into linear terms. Analytical expressions of the viscoelasticity of blood samples are obtained by solving the differential equation. Here, viscosity and compliance are derived analytically. To demonstrate the present method, a T-shaped microfluidic channel is used. When measuring velocity fields of blood sample, a T-shaped microfluidic channel is not required to align a microscopic image in the horizontal or vertical direction. It consists two inlets, one outlet, two guiding channels (blood sample channel, reference fluid channel), and one coflowing channel. Using two syringe pumps, a blood sample is infused into the blood sample channel with a sinusoidal flow-rate pattern. A reference fluid is supplied into the reference fluid channel at a constant flow rate. By monitoring the interface of both fluids in the coflowing channel, the viscosity and compliance are obtained at an interval of a specific period. As a performance demonstration, the present method was used to evaluate the contributions of period (T), diluents (plasma, 1x phosphate-buffered saline (PBS), and hematocrit (Hct) to viscoelasticity. The present method was then employed to quantify the viscoelasticity of a fixed blood sample prepared by adding fixed RBCs into plasma.

Blood Sample Preparation
According to the ethics committee of Chosun University Hospital (CHOSUN 2018-05-11), all experiments were conducted after ensuring that the experimental protocols were appropriate and humane.
Human concentrated RBCs and fresh frozen plasma (FFP) were purchased from the Gwangju-Chonnam blood bank (Gwangju, Korea) and were stored at 4 • C and −20 • C, respectively. Because the RBCs were preserved in citrate phosphate dextrose adenine (CPDA) as an anticoagulant solution, it was necessary to remove CPDA from the concentrated RBCs. The concentrated RBCs (~7 mL) were added into 1x PBS (pH 7.4, Gibco, Life Technologies, Carlsbad, CA, USA) (~7 mL) in a 15-mL tube. After the tube was inserted into a centrifugal separator (Allegra X-30R benchtop, Beckman Coulter, Brea, CA, USA), it was set to 4000 rpm and operated for 10 min. The diluted blood was separated into two layers: an upper layer (plasma), and a lower layer (RBCs). Pure RBCs were collected after removing liquid in the upper layer. Additionally, FFP was thawed at room temperature (25 • C). Plasma was filtered to remove cellular debris and unwanted white blood cells with a syringe filter (mesh size = 5 µm, Minisart, Sartorius, Göttingen, Germany). The RBCs and plasma were stored at 4 • C in a refrigerator before the blood test [20].
First, to evaluate the effect of the contribution of Hct and diluents (1x PBS, plasma) on the viscoelasticity of blood samples, blood samples with Hct = 30%, 40%, 50%, and 60% were prepared by adding normal RBCs into 1x PBS or plasma. Except in the experiment for evaluating the contribution of Hct, all blood samples were adjusted to Hct = 50%. Second, to fix normal RBCs chemically, three different concentrations of glutaraldehyde (GA) solution (C GA = 4, 8, and 12 µL/mL) were diluted by mixing GA solution (Grade II, 25% in H 2 0, Sigma-Aldrich, St. Louis, MO, USA) into 1x PBS. Normal RBCs were fixed for consistent measurement because RBCs needed to be unchanged over experimental time. To fix normal RBCs, normal RBCs were mixed with each concentration of GA solution for 10 min prior to washing them. Fixed RBCs were collected after a washing procedure. The fixed blood sample (Hct = 50%) was then prepared by adding the fixed RBCs into plasma. Here, to evaluate the contribution of fixed RBCs to viscoelasticity effectively, it was necessary to remain constant at a level of hematocrit (i.e., Hct = 50%).

Fabrication on a Microfluidic Device and Experimental Procedure
A T-shaped microfluidic device for measuring blood viscoelasticity consisted of two inlets (a, b), one outlet (a), two guiding channels (blood sample channel (BC), reference fluid channel (RC)), and one coflowing channel (CC), as shown in Figure 1A-a and Figure S1 (Supplementary Materials). The blood sample channel (width = 250 µm, length = 7500 µm) and reference fluid channel (width = 250 µm, length = 7500 µm) were connected to the coflowing channel (width = 250 µm, length = 9200 µm). Here, dimensions of a microfluidic channel were selected to measure velocity fields and blood viscosity accurately. First, velocity fields of blood flows were obtained accurately with microscopic images captured with at least 10× objective lens. Based on fields of view, channel width and length were selected suitably. Second, a rectangular channel with low aspect to ratio was preferred to measure blood viscosity effectively under coflowing method. The channel depth of the microfluidic device was fixed at 20 µm.

Quantification of Interface (αB), Averaged Blood Velocity (<UBC>), and Averaged Image Intensity (<IBC>)
Variations of the interface in the coflowing channel were used to quantify the viscoelasticity of the blood sample. Additionally, the image intensity and velocity fields of the blood sample flowing in the blood sample channel were obtained to evaluate the erythrocyte sedimentation rate (ESR) that occurred in the driving syringe while the blood flow rate was controlled by the syringe pump. First, to obtain the interface between the blood sample and 1x PBS in the coflowing channel, a specific ROI  Conventional microelectromechanical-system fabrication techniques (photolithography and deep reactive ion etching) were used to fabricate a master mold on a 4-inch silicon wafer. polydimethylsiloxane (PDMS) (Sylgard 184, Dow Corning, Midland, MI, USA) prepolymer and a curing agent were mixed at a ratio of 10:1. After the mold was fixed in a Petri dish, the PDMS mixture was poured on the master mold. Air bubbles in the PDMS were removed with a vacuum pump for 1 h. After curing the PDMS mixture in a convective oven (70 • C for 1 h), a PDMS block was peeled from the master mold. It was cut with a razor blade. Three ports (two inlets and one outlet) were punched with a biopsy punch (outer diameter = 1.2 mm). After oxygen-plasma treatment on the PDMS block and a glass slide with an oxygen-plasma system (CUTE-MPR, Femto Science Co., Gyeonggi, Korea), a microfluidic device was finally prepared by bonding the PDMS block on the glass slide.
As shown in Figure 1A-b, two polyethylene tubes (L 1 , inner diameter = 500 µm, thickness = 500 µm, and length = 300 mm) were connected from two disposable syringes (~1 mL) to two inlets (a and b). The other polyethylene tube (L 2 , inner diameter = 500 µm, thickness = 500 µm, and length = 200 mm) was connected from an outlet (a) to a waste collection unit. To remove air bubbles in the channels and avoid nonspecific binding of plasma proteins to the inner surfaces of the channels, all channels were filled completely with bovine serum albumin (BSA) solution of C BSA = 2 mg/mL through outlet (a) with a disposable syringe. After 10 min, all channels were rinsed then filled with 1x PBS. After two disposable syringes (~1 mL) were filled with blood sample (~1 mL) and 1x PBS (~1 mL), they were installed into two syringe pumps (neMESYS, Cetoni GmbH, Germany). The blood sample was supplied into inlet (a) at a sinusoidal flow rate (Q B (t) = Q 0 + Q 1 ·sin (2πt/T)). Q 0 and Q 1 are the average and amplitude of the sinusoidal flow rate. Additionally, T represents the period. Reference fluid (1x PBS) was supplied into inlet (b) at a constant flow rate (Q R (t) = Q 0 ).
As shown in Figure 1A-c, the microfluidic device was positioned on an optical microscope (BX51, Olympus, Tokyo, Japan) equipped with a 10× objective lens (NA = 0.25). A high-speed camera (FASTCAM MINI, Photron, USA) was used to capture microscopic images of the blood sample and 1x PBS flowing in microfluidic channels. It offered a spatial resolution of 1280 × 1000 pixels. Each pixel corresponded to 10 µm. With a function generator (WF1944B, NF Corporation, Yokohama, Japan), a pulse signal with a period of 1 s triggered the high-speed camera. Microscopic images were sequentially captured at a frame rate of 5 kHz. All experiments and blood sample preparations were conducted at a room temperature of 25 • C.

Quantification of Interface (α B ), Averaged Blood Velocity (<U BC >), and Averaged Image Intensity (<I BC >)
Variations of the interface in the coflowing channel were used to quantify the viscoelasticity of the blood sample. Additionally, the image intensity and velocity fields of the blood sample flowing in the blood sample channel were obtained to evaluate the erythrocyte sedimentation rate (ESR) that occurred in the driving syringe while the blood flow rate was controlled by the syringe pump. First, to obtain the interface between the blood sample and 1x PBS in the coflowing channel, a specific ROI of 150 × 750 pixels was selected in the coflowing channel, as shown in Figure 1B-a. To obtain the interface in the coflowing channel effectively, a gray-scale image was converted into a binary-scale image by adopting Otsu's method [21]. As shown in Figure 1B-b, the blood-filled width (W B ) over the ROI was calculated by using a commercial software package (MATLAB 2019, MathWorks, Natick, MA, USA). The interface (α B ) was obtained as α B = W B /W. Here, W is the channel width of the coflowing channel. Second, to monitor variations of blood flow rate supplied from the syringe pump, velocity fields of the blood sample flowing in the blood sample channel were obtained by conducting a time-resolved micro particle image velocimetry (micro-PIV) technique. A specific ROI of 300 × 150 pixels was selected in the blood sample channel. The size of the interrogation window was 32 × 32 pixels. The window overlap was 50%. The obtained velocity fields were validated with a median filter. The averaged velocity (<U BC >) was calculated as an arithmetic average of U BC distributed over the ROI. Third, to evaluate the ESR that occurred in the driving syringe, it was necessary to quantify the microscopic image intensity of the blood sample flowing in the blood sample channel. A specific ROI with 300 × 150 pixels was selected in the blood sample channel. The image intensity of the blood sample flowing in the blood sample channel was obtained by conducting digital image processing with MATLAB. An averaged image intensity (<I BC >) was obtained by averaging variations of I BC distributed over the ROI.

Discrete Fluidic Circuit for Representing Viscoelasticity of Blood Sample
Blood samples were assumed to be Newtonian fluids. To evaluate the viscoelasticity of the blood sample, a simple mathematical model was constructed with discrete fluidic circuit elements. As shown in Figure 1B-c, two fluids flowing in the coflowing channel (i.e., blood sample stream, reference fluid stream) are represented with individual discrete fluidic circuit elements. The fluid circuit model is composed of a flow-rate element (Q B , Q R ), resistance element (R B and R R ), and compliance element (C B ). Q B and Q R are the flow rates of the blood sample and reference fluid, respectively. Ground represented zero value of pressure (P = 0).
The blood stream for representing the viscoelasticity of the blood sample was modeled as a resistance element (R B ) and compliance element (C B ) combined in parallel. The C B was included to account for the compliance effect of the RBCs, the microfluidic channel, and a connected tube. However, because the reference fluid flowed at a constant flow rate in the coflowing channel, the compliance effect of a microfluidic channel and tube filled with reference fluid was negligible. For this reason, the reference stream was simply modeled only as a single resistance element (R R ). The coflowing channel was partially filled with the blood sample stream (W B ) and reference fluid stream (W -W B ). Because ground represented zero value of pressure, both streams had the same pressure drop (P B = P R = ∆P = P). Based on mass conservation for the blood sample stream and reference fluid stream in the coflowing channel, two equations were derived: for the blood sample stream, and for the reference fluid stream. By inserting Equation (2) into Equation (1), a first-order ordinary differential equation was derived: Because a rectangular-shaped channel (width = w, depth = h, and length = l) with a lower aspect ratio (AR) (i.e., AR = depth/width = 20/250) was filled with fluid (viscosity = µ), the resistance element was modeled approximately as [21]: Based on the analytical expression of a rectangular channel, the corresponding resistance element for each stream was derived as for the blood sample stream and for the reference fluid stream. In Equation (6), the CF was included to compensate for the boundary condition difference between the real physical model and approximate circuit model [17]. When inserting Equations (5) and (6) into Equation (4), the following equation was derived.
Micromachines 2020, 11, 245 7 of 19 The first part in the right side of Equation (7) was expressed again in a different form.
In Equation (8), R WB , which assumed that the coflowing channel was filled with the blood sample, was given as R WB = 12µ B l cc Wh 3 . According to a previous study, blood viscosity remained constant with respect to the interface [18]. Because R WB and µ R µ B were independent of time, Equation (8) became a simple expression of Equation (9).
In Equation (10), CF was varied depending on the interface (α B ) (i.e., CF = CF (α B )). Here, the CF was obtained by conducting numerical simulation. Because Equation (10) had nonlinear terms, it was substantially difficult to find an analytical solution. For this reason, it was necessary to approximate the nonlinear Equation (10) as a simple linear equation. Equation (10) was modified as a simple form.
In Equation (11), F 1 and F 2 were obtained by obtaining the weighted average of CF (i.e., within a specific value of the interface. As the interface was relocated periodically within a specific range (0.1 < α B < 0.9), two approximate factors (F 1 and F 2 ) with constant values were calculated from Equations (12) and (13).
By dividing Equation (11) with F 2 , Equation (11) was changed to a simple linear differential equation.
In Equation (14), β B and time constant (τ B ) were derived as In this study, the flow rates of the blood sample and reference fluid were controlled as Q B (t) = Q 0 + Q 1 sin(ωt) and Q R (t) = Q 0 , respectively. Here, ω was given as ω = 2π T . The particular solution of Equation (14) was then derived as In Equation (15), β 0 and β 1 were given as Equations (16) and (17).
Micromachines 2020, 11, 245 8 of 19 and Additionally, time delay (ϕ) was given as ϕ = ωλ B . From Equation (16), the blood viscosity (µ B ) was given as Additionally, from Equation (17), the time constant (λ B ) was derived as According to λ B = R WB C B F 1 F 2 , the analytical expression of compliance (C B ) was derived as In other words, if β 0 and β 1 could be obtained from periodic variations of the interface (α B ) in the coflowing channel, µ B and C B as blood viscoelasticity could be evaluated from Equations (18) and (20), respectively.
As shown in Figure 1C, as a preliminary study, a blood sample (Hct = 50%) was prepared by adding normal RBCs into 1x PBS. The blood sample was supplied into inlet (a) at a sinusoidal flow rate (Q B (t) = 1 + 0.5 sin (2πt/360) mL/h). Simultaneously, 1x PBS was supplied into inlet (b) at a constant flow rate of Q R (t) = 1 mL/h. Figure 1C-a showed temporal variations of <U BC > and α B over time.

Correction Factor (CF) and Approximate Factors (F 1 and F 2 ) via Numerical Simulation
To find two approximate factors (F 1 , F 2 ) in Equation (11), it was necessary to obtain the CF by conducting numerical simulation with CFD software (CFD-ACE+, ESI Group, Paris, France). For convenience, the flow rate of the blood sample was assumed to be Q B = 1 mL/h. The viscosities of both fluids (blood sample, reference fluid) were assumed as µ B = 1 cP and µ R = 1 cP, respectively. Figure 2A showed variations of the interface (α B ) through numerical simulation with respect to the flow-rate ratio (

Correction Factor (CF) and Approximate Factors (F1 and F2) via Numerical Simulation
To find two approximate factors (F1, F2) in Equation (11), it was necessary to obtain the CF by conducting numerical simulation with CFD software (CFD-ACE+, ESI Group, Paris, France). For convenience, the flow rate of the blood sample was assumed to be QB = 1 mL/h. The viscosities of both fluids (blood sample, reference fluid) were assumed as μB = 1 cP and μR = 1 cP, respectively. Figure    According to the coflowing-streams method [22,23], the viscosity ratio of the blood sample to the reference fluid (µ B /µ R ) was obtained as 1−α B by quantifying the interface (α B ) in coflowing streams at the same flow-rate condition (Q B = Q R ). The normalized viscosity of the blood sample was obtained by dividing the estimated viscosity (µ est = µ R ·α B /(1 − α B )) by the given viscosity (µ given = 1 cP) (µ n = µ est /µ given ). As shown in Figure 2B, variations of µ n were obtained by varying the interface (0.1 < α B < 0.9). When the interface was located at the center line (α B = 0.5), µ n was given as µ n = 1. In other words, the blood viscosity could be measured accurately when the interface was located at the center line [24]. However, when α B was relocated from the center line to both walls, µ n tended to decrease by approximately 0.8. Because the viscosity of the blood sample was given as µ B = 1 cP, the coflowing-streams method exhibited a large measurement error of approximately 20% when compared with the given viscosity of the blood sample. The reason could be explained by the boundary condition difference between the real physical model and simple mathematical model. Instead of a real and complex model, to construct a simple model of coflowing streams, the coflowing-streams method assumed that the interface of two streams was a virtual wall. In other words, it assumed that the interface was a virtual-wall boundary. Because the CF varied by channel dimensions (width and depth), the CF was calculated by referring to the general procedure discussed in previous studies [17,19]. The CF was then estimated by the reciprocating µ n obtained at a specific interface (CF·(α n ) = µ n −1 for α n ). As shown in Figure 2B, the CF was obtained as CF = 1 at the center line (α B = 0.5). However, the CF tended to increase gradually when the interface moved to both walls. According to regression analysis, the CF was obtained as CF = 4.2162 α B 4 -8.4325 α B 3 + 7.0763 α B 2 − 2.86 α B + 1.4599 (R 2 = 0.9912). By inserting the correction factor into the coflowing-streams method, the viscosity of the blood sample was estimated with µ B = µ R CF(α B ) α B 1−α B . The coflowing method with correction factor could be used to measure the viscosity of blood with a specific hematocrit. Based on Equation (12), the approximate factor (F 1 ) was obtained as F 1 = 1.094. As shown in Figure 2C, variations of CF(α B ) 1 1−α B and F 1 1 1−α B were obtained with respect to α B . The normalized difference (ND) between both terms exhibited its maximum value at the center and both walls. The normalized difference was less than 10%. Additionally, according to Equation (13), the approximate factor (F 2 ) was obtained as F 2 = 1.1087. As shown in Figure 2D, variations of CF(α B ) α B 1−α B and F 2 α B 1−α B were obtained with respect to α B . The maximum value of normalized difference was estimated to be approximately 11%. The simulation study showed that the two approximate factors (F 1 =1.094, F 2 =1.1087) could give consistent results when compared with the original expression. Using F 1 and F 2 , the nonlinear Equation (10) was converted into the simple linear Equation (11) for consistency.

Effect of Period (T) on Viscoelasticity of Blood Sample
To verify the contribution of the period (T) to the viscoelasticity of the blood sample (Equations [18] and [20]), the viscosity and compliance were evaluated by varying the period (T = 120, 240, 360, and 480 s). The blood sample (Hct = 50%) was prepared by adding normal RBCs into 1x PBS. Q 0 and Q 1 of the two syringe pumps were controlled at Q 0 = 1 and Q 1 = 0.5 mL/h, respectively. For a rectangular channel with a lower aspect ratio, the shear rate ( . γ) was derived as . γ = 6Q wh 2 [19]. Based on the shear rate formula, the shear rates of the corresponding flow rate were estimated as . γ = 8333 s −1 for Q B = 0.5 mL/h and . γ = 25, 000 s −1 for Q B = 1.5 mL/h. Because the shear rate ( . γ) was much greater than 1000 s -1 , it was reasonable that the blood sample behaves as a Newtonian fluid. In other words, the blood viscosity (µ B ) remained constant within the specific flow rates of the blood sample.
As shown in Figure Figure 3B-a showed variations of β 0 and β 1 with respect to T, where β 0 and β 1 fluctuated at a shorter period (T = 120, 360 s). However, they remained stable at a longer period (T = 360, 480 s). Figure 3B-b showed variations of λ B with respect to T. The λ B tended to increase linearly for up to T = 360 s. However, the slope of λ B tended to decrease between T= 360 and T = 480 s. According to Equation (19), the λ B was linearly proportional to the period (i.e., λ B~T ). The experimental results showed appropriately consistent variations of λ B with respect to T. Using Equations (18) and (20), the blood viscosity (µ B ) and compliance (C B ) were obtained with respect to T. As shown in Figure 3B-c, µ B did not exhibit a linear dependency of T. It fluctuated at a shorter period. However, it remained constant at a longer period (T = 360, 480 s). The results agreed with Equation (18), which did not relate to the period. It was necessary, however, to set a longer period for consistently measuring the viscosity of the blood sample. Compliance (C B ) tended to increase linearly for up to T = 360 s. The slope of C B tended to decrease between T = 360 and T = 480 s. According to the mathematical relation, C B was linearly proportional to λ B . The experimental results indicated that blood viscosity was independent of period. However, compliance varied linearly depending on the period.
From the results, for consistent measurement of viscoelasticity (blood viscosity and compliance), the period of the syringe pump was set to a longer period of T = 360 s throughout all experiments for convenience. T = 480 s. According to the mathematical relation, CB was linearly proportional to λB. The experimental results indicated that blood viscosity was independent of period. However, compliance varied linearly depending on the period.
From the results, for consistent measurement of viscoelasticity (blood viscosity and compliance), the period of the syringe pump was set to a longer period of T = 360 s throughout all experiments for convenience.

Quantification of the Effect of Hematocrit on Blood Viscoelasticity
According to a previous study, hematocrit caused an increase in blood viscosity and elasticity [16]. In addition, a blood sample with a low hematocrit (Hct = 30%) exhibited a continuous ESR occurring in a driving syringe [25]. According to the previous study, to increase the ESR significantly, a blood sample was prepared by adding normal RBCs into various concentrations of dextran solution (C dex = 2, 5, 8, and 10 mg/mL). Because the hematocrit of the blood sample flowing in the microfluidic channel varied continuously over time, RBC aggregation or blood viscosity tended to vary continuously [19]. In this study, under blood perfusion with a sinusoidal flow-rate pattern, the contribution of hematocrit to blood viscoelasticity and ESR was quantified by varying the hematocrit. To induce the ESR in a driving syringe, plasma was used as the diluent. In other words, blood samples (Hct = 30%, 40%, 50%, and 60%) were prepared by adding normal RBCs into the plasma. Figure 4A showed temporal variations of β B with respect to hematocrit ((a) Hct = 30%, (b) Hct = 40%, (c) Hct = 50%, and (d) Hct = 60%). Using temporal variations of β B , blood viscosity and compliance were obtained at an interval of a specific period. As shown in Figure 4B, temporal variations of µ B and C B were obtained with respect to hematocrit ((a) Hct = 30%, (b) Hct = 40%, (c) Hct = 50%, and (d) Hct = 60%). Both parameters were obtained as mean ± standard deviation at a specific time (t = 290, 650, 1010, 1370, and 1730 s).
contribution of hematocrit to blood viscoelasticity and ESR was quantified by varying the hematocrit. To induce the ESR in a driving syringe, plasma was used as the diluent. In other words, blood samples (Hct = 30%, 40%, 50%, and 60%) were prepared by adding normal RBCs into the plasma. Figure 4A showed temporal variations of βB with respect to hematocrit ((a) Hct = 30%, (b) Hct = 40%, (c) Hct = 50%, and (d) Hct = 60%). Using temporal variations of βB, blood viscosity and compliance were obtained at an interval of a specific period. As shown in Figure 4B, temporal variations of μB and CB were obtained with respect to hematocrit ((a) Hct = 30%, (b) Hct = 40%, (c) Hct = 50%, and (d) Hct = 60%). Both parameters were obtained as mean ± standard deviation at a specific time (t = 290, 650, 1010, 1370, and 1730 s). For the blood sample with Hct = 30%, compliance (CB) fluctuated greatly over time. It tended to decrease after t = 290 s, but it tended to increase between t = 650 and t = 1370 s. After t = 1370 s, it For the blood sample with Hct = 30%, compliance (C B ) fluctuated greatly over time. It tended to decrease after t = 290 s, but it tended to increase between t = 650 and t = 1370 s. After t = 1370 s, it remained constant over time. Blood viscosity tended to decrease after t = 650 s. In other words, the ESR of RBCs in a driving syringe accelerated over time. In other words, after a certain time, the hematocrit of the blood sample tended to decrease. Thus, µ B tended to decrease, and C B tended to increase with an elapse of time. For the blood sample with Hct = 40%, blood viscosity tended to decrease after t = 1010 s. Compliance remained constant for up to t = 1010 s. After that, it fluctuated over time. When the hematocrit increased from Hct = 30% to Hct = 40%, the time when C B had the minimum value tended to increase substantially. However, for the blood samples with high hematocrit (Hct = 50%, 60%), blood viscosity and compliance remained constant over time. The blood viscosity tended to increase at a higher hematocrit. The compliance tended to decrease at a higher hematocrit. According to a previous study, blood elasticity tended to increase with respect to hematocrit [16]. When compared with the previous result, compliance tended to decrease with respect to the hematocrit. The experimental results can be considered reasonable, because compliance had the reciprocal of elasticity (i.e., compliance1 /elasticity).
To quantify the ESR that occurred in a driving syringe, the microscopic image intensities of the blood sample were obtained over time. A driving syringe was installed horizontally. Because of the ESR in the driving syringe, the hematocrit of the blood sample flowing in a microfluidic channel tended to decrease over time. Here, variations of hematocrit were used by monitoring the image intensity of the blood sample [26]. Figure 4C-a showed microscopic images of the blood sample (Hct = 30%) flowing in the blood sample channel over time (t = 0, 480, 960, 1440, and 1920 s). After t = 960 s, the number of RBCs tended to decrease significantly. The image intensity tended to increase. At t = 1920 s, the image intensity tended to decrease as RBCs tended to increase significantly. The contrast of each image was enhanced by conducting image processing with the software Image-J (NIH, Maryland, USA). A specific ROI (300 x 150 pixels) in the blood sample channel was selected to quantify the image intensity (<I BC >). To visualize the ESR in a driving syringe, a side view of a syringe filled with a blood sample (Hct = 30%) was captured sequentially with a smartphone camera (Galaxy A5, Samsung, Korea). As lower level of hematocrit exhibited higher value of <I BC >, the results of blood sample (Hct = 30%) was selected and summarized at a specific time. Figure 4C-b showed snapshots of a driving syringe over time (t) (t = 0, 480, 960, 1440, and 1920 s). Before t = 480 s, there was no indication of the ESR occurring in a driving syringe, because the RBCs were distributed uniformly. After t = 960 s, because of the continuous ESR, the blood sample inside the syringe was separated into two regions: an RBC-rich region (i.e., lower layer) and an RBC-free region (i.e., upper layer). As shown in Figure 4C-a, the ESR in the driving syringe caused a reduced number of RBCs in the microfluidic channel. To quantify the decrease in hematocrit resulting from ESR in the driving syringe, the image intensity of the blood sample (<I BC >) was obtained over time. Figure 4C-c showed temporal variations of <I BC > with respect to Hct = 30%, 40%, 50%, and 60%. For a blood sample with Hct = 30%, <I BC > tended to increase after t = 1190 s. For a blood sample with Hct = 40%, image intensity tended to increase after t = 1570 s. For blood samples with a high hematocrit (Hct = 50%, 60%), the image intensity remained constant over time. This result showed that a lower hematocrit contributed to varying image intensity (or numbers of RBCs) of blood samples flowing in a microfluidic channel. According to a specific parameter (S ESR ) suggested in a previous study [20], variations of ESR were quantified with respect to Hct. The S ESR was obtained as S ESR = t=1920 t=0 ( I BC − I min )dt. Here, I min represents the minimum value of <I BC > within a specific duration, t = 1920 s. Figure 4C-d showed variations of S ESR with respect to Hct. The S ESR tended to decrease substantially with respect to Hct. The blood sample with Hct = 30% showed the maximum value of S ESR . This result indicated that S ESR varied significantly depending on the hematocrit.
From the result, one can conclude that the viscoelasticity of the blood sample suggested by the present method can be varied with the hematocrit. In addition, it can be employed to quantify variations of the ESR occurring in the driving syringe by monitoring temporal variations of viscoelasticity.

Quantification of the Contribution of Hardened RBCs to Blood Viscoelasticity
Finally, the method was employed to quantify the contribution of hardened RBCs to the viscoelasticity of blood samples. According to previous studies [27,28], normal RBCs were hardened chemically with GA solution. The degree in rigidity increased gradually by varying concentrations of the GA solution. Normal RBCs were hardened chemically with three different concentrations of GA solution (C GA ) (C GA = 4, 8, and 12 µL/mL). The hardened blood sample (Hct = 50%) was then prepared by adding hardened RBCs into plasma. Figure 5A showed temporal variations of β B with respect to C GA . As the concentration of GA solution increased, β B tended to increase gradually. In addition, the β B exhibited steady and periodic variations over time. Figure 5B showed temporal variations of <I BC > with respect to C GA . The inset showed a microscopic image of the fixed blood sample (i.e., fixed RBCs with C GA = 12 µL/mL) captured at t = 480 s. The <I BC > tended to decrease slightly at a higher concentration of GA solution. However, it exhibited steady and periodic variations over time. The result indicated that the fixed blood sample with Hct = 50% did not induce the ESR in the driving syringe as shown in Figure 4.
GA solution (CGA) (CGA = 4, 8, and 12 μL/mL). The hardened blood sample (Hct = 50%) was then prepared by adding hardened RBCs into plasma. Figure 5A showed temporal variations of βB with respect to CGA. As the concentration of GA solution increased, βB tended to increase gradually. In addition, the βB exhibited steady and periodic variations over time. Figure 5B showed temporal variations of <IBC> with respect to CGA. The inset showed a microscopic image of the fixed blood sample (i.e., fixed RBCs with CGA = 12 μL/mL) captured at t = 480 s. The <IBC> tended to decrease slightly at a higher concentration of GA solution. However, it exhibited steady and periodic variations over time. The result indicated that the fixed blood sample with Hct = 50% did not induce the ESR in the driving syringe as shown in Figure 4.  Figure 5C showed variations of β0 and β1 with respect to CGA. Both parameters (β0 and β1) tended to increase substantially at a higher concentration of GA solution. Because the GA solution was employed to increase the rigidity of RBCs, both parameters presented distinctive variations of hardness. Based on Equations (18)- (20), variations of μB and CB were obtained at an interval of T = 360 s. As βB and <IBC> showed stable variations over time, μB and CB remained constant over the specific duration of the test as shown in Figures 5A and 5B. Thus, μB and CB were represented as average ± standard deviation (n = 5). Figure 5D showed variations of μB and CB with respect to CGA. According to previous studies [17,18], blood viscosity and elasticity tended to increase at a higher concentration of GA solution. When compared with the previous results, blood viscosity exhibited consistent variations with respect to the concentration of GA solution. Because compliance was defined as the reciprocal of elasticity, the compliance also showed consistent variations with respect to the concentration of GA solution.  Figure 5C showed variations of β 0 and β 1 with respect to C GA . Both parameters (β 0 and β 1 ) tended to increase substantially at a higher concentration of GA solution. Because the GA solution was employed to increase the rigidity of RBCs, both parameters presented distinctive variations of hardness. Based on Equations (18)- (20), variations of µ B and C B were obtained at an interval of T = 360 s. As β B and <I BC > showed stable variations over time, µ B and C B remained constant over the specific duration of the test as shown in Figure 5A,B. Thus, µ B and C B were represented as average ± standard deviation (n = 5). Figure 5D showed variations of µ B and C B with respect to C GA . According to previous studies [17,18], blood viscosity and elasticity tended to increase at a higher concentration of GA solution. When compared with the previous results, blood viscosity exhibited consistent variations with respect to the concentration of GA solution. Because compliance was defined as the reciprocal of elasticity, the compliance also showed consistent variations with respect to the concentration of GA solution.
The results lead to the conclusion that the present method could be employed to monitor variations of the viscoelasticity of blood samples while the syringe pump was set to a pulsatile flow-rate pattern.

Quantificative Comparison of Viscoelasticty Obtained with Preesent Method and Conventional Viscometer
Using a conventional viscometer, viscoelasticity of cells was modeled with the linear Maxwell model. As shown in Figure 6A-a, the viscoelasticity of each blood sample was modeled as a solid element and a fluid element connected in series. The corresponding constitutive expression of each element was given as τ = Gγ (solid element) and τ = µ . γ (fluid element), respectively. Here, G and µ represented elasticity and viscosity, respectively. τ and γ denoted shear stress and shear strain. External shear strain was excited periodically as γ(t) = γ 0 e jωt . The governing equation was then derived as γ. Here, the time constant of viscometry (l cv ) was expressed as l cv = µ/G.
The viscous effect of the viscometer was considered as negligible since the viscometer did not have an influence on relaxation time of the cells. The viscometer had operated at a wider range of frequency from ω = 0.3 rad/s to ω = 700 rad/s [29]. The previous study indicated that the time constant of whole blood was obtained as λ cv = 1.5-13.4 ms for shear flow [29], and λ cv = 114-259 ms for extensional flow [30].
Using a conventional viscometer, viscoelasticity of cells was modeled with the linear Maxwell model. As shown in Figure 6A-a, the viscoelasticity of each blood sample was modeled as a solid element and a fluid element connected in series. The corresponding constitutive expression of each element was given as τ = Gγ (solid element) and τ = μ̇ (fluid element), respectively. Here, G and μ represented elasticity and viscosity, respectively. τ and γ denoted shear stress and shear strain. External shear strain was excited periodically as γ(t) = 0 . The governing equation was then derived as ̇+ 1 =̇. Here, the time constant of viscometry (lcv) was expressed as lcv = μ/G. The viscous effect of the viscometer was considered as negligible since the viscometer did not have an influence on relaxation time of the cells. The viscometer had operated at a wider range of frequency from ω = 0.3 rad/s to ω = 700 rad/s [29]. The previous study indicated that the time constant of whole blood was obtained as λcv = 1.5-13.4 ms for shear flow [29], and λcv = 114-259 ms for extensional flow [30]. On the other hand, under periodic blood flow in the microfluidic system, Figure 6A On the other hand, under periodic blood flow in the microfluidic system, Figure 6A-b showed the simple fluidic circuit model of the microfluidic system. Based on electric circuit analysis, the governing equation of fluidic system was derived as λ B . Q + Q = Q B (t). Here, Q represented the flow rate which passed through resistance element. Time constant (λ B ) was derived as λ B = R WB ·C B . R WB and C B represented fluidic resistance of the coflowing channel filled with only blood sample and compliance, respectively. As shown in Figure 6A-b, a microfluidic system behaved as an R-C low pass filer. To effectively infuse alternating components in the periodic flow rate into a microfluidic system, the period of the excitation flow rate should be much longer than the time constant of the microfluidic system (i.e., T > λ B ). In addition, the syringe pump used in this study did not infuse blood samples during short periods. From the experimental results as shown in Figure 3B, period (T) of the sinusoidal flow rate was fixed as T = 360 s.
The λ B was determined by R WB and C B . Flexible tubing and the PDMS channels tended to vary the time constant substantially [11]. A microfluidic channel with different channel depths (H) (H = 4, 10, and 20 µm) was prepared to change R WB . Here, 1x PBS was infused into a microfluidic channel to reduce or remove the viscoelastic effects. As shown in Figure 6B-a, the corresponding time constant for each channel depth was obtained as λ B = 5.92 ± 0.81 s (H = 4 µm), λ B = 6.11 ± 0.20 s (H = 10 µm), and l B = 3.64 ± 0.40 s (H = 20 µm). In addition, Figure 6B-b showed variations of R WB and C B with respect to H. A lower channel depth contributed to increasing R WB , and decreasing C B . The result indicated that fluidic resistance (R WB ) had a strong influence on C B . To find out the contribution of compliance element (C B ), the compliance of the microfluidic system increased intentionally by securing the air cavity inside the driving syringe. As shown in Figure 6C, variations of λ B and C B were obtained with respect to V air = 0 and 0.1 mL. When the air cavity of 0.1 mL existed inside the driving syringe, λ B and C B increased considerably as λ B = 20.84 ± 3.42 s and C B = 378.59 ± 62.11 s. From the result, the air cavity tended to increase λ B and C B substantially. When compared with viscometry data, the time constant (λ B ) increased about O (10 2 ) significantly because of the compliance effect of the microfluidic system (i.e., flexible tubing, PDMS channels, and the air cavities existing in the driving syringe). In addition, minimum threshold of C B (i.e., detection limit) was estimated as 66.05 ± 7.30 µm 3 /mPa at a specific condition (i.e., V air = 0, H = 20 µm). According to a previous study [17], time constants obtained with two different systems (i.e., λ PM : microfluidic system, λ CPV : conventional viscometer) were obtained and compared with respect to C glycerin = 10%, 20%, 30% and 40%. As glycerin solution did not include viscoelasticity, both time constants remained unchanged with respect to different concentrations of glycerin solution. However, the microfluidic system had a longer time constant with O (10 0 ). The ratio of time constant between λ CPV and λ PM was obtained as λ CPV /λ PM = O (10 2 ). To quantitatively determine elasticity obtained with both methods, a scatter plot was employed by plotting G PM on the vertical axis and G CPV on the horizontal axis. According to linear regression analysis, regression coefficient (R 2 ) had a higher value (R 2 = 0.9617). This result indicated that the elasticity obtained with microfluidic system exhibited consistent variations with respect to C glycrin when compared with elasticity obtained with conventional viscometer. Thus, the microfluidic system could be employed to measure viscoelasticity effectively. The slope of 0.0022 indicated that elasticity obtained with the microfluidic system was much less than that obtained with the conventional viscometer.
According to order analysis, C B had an order of O (10 −13 ) from the analytical expression of time constant. Here, O represented order. According to experimental results as shown in Figure 5 . Thus, C B obtained for the blood sample had an order of O (10 −13 ). In order words, both approaches (i.e., analytical expression, and experimental data) exhibited the same order of O (10 −13 ). Thus, the present method could be used to monitor C B of blood samples sufficiently.
To compare the relationship between G B and C B , it was assumed that λ cv of the conventional viscometer had the same λ B as the microfluidic system (i.e., λ cv = λ B ). The time constant obtained by the microfluidic system was used to evaluate G B and C B simultaneously. The following relation was given as µ G = R WB × C B . The analytical expression indicated that G B and C B had a reciprocal relationship (i.e., G~1/C B ). Using experimental results as shown in Figure 4, G B and C B were obtained with respect to Hct = 30%, 40%, 50%, and 60%. Additionally, by referring to the previous study [18], variations of G B were represented with respect to Hct = 30%, 40%, and 50%. Here, blood samples were prepared by adding normal RBCs into plasma. As shown in Figure 7A, at Hct > 40%, G B tended to increase with respect to Hct. Inversely, C B tended to decrease with respect to Hct. When compared with the previous study, the trend of G B increased similarly with respect to Hct. It can be inferred that different microfluidic systems contributed to differences of G B between both studies. Additionally, variations of C B and G B were summarized with respect to C GA as shown in Figure 7B. Here, in the previous study [4], fixed blood samples were prepared by adding fixed RBCs into 1x PBS instead of plasma. G B of both studies tended to increase gradually with respect to C GA . From quantitative comparisons between the previous study and the present study, elasticity (G B ) and compliance (C G ) had a reciprocal relationship. Additionally, they varied significantly when the rigidity of RBCs increased substantially.
previous study [4], fixed blood samples were prepared by adding fixed RBCs into 1x PBS instead of plasma. GB of both studies tended to increase gradually with respect to CGA. From quantitative comparisons between the previous study and the present study, elasticity (GB) and compliance (CG) had a reciprocal relationship. Additionally, they varied significantly when the rigidity of RBCs increased substantially. Viscoelasticity (G) was represented as G = G1+ j G2. Here, G1 (storing modulus) and G2 (loss modulus) were expressed as 1 = G ( 2 2 1+ 2 2 ) and 2 = G ( 1+ 2 2 ). Variations of G1 and G2 were represented with respect to radial frequency (ω). G1 and G2 tended to vary depending on ·ω. However, as ·ω showed a significant difference for both systems, it was apparent that both systems exhibited different variations of G1 and G2 with respect to ω. The microfluidic system and conventional viscometer showed different characteristics in terms of angular frequency (or period) and time constants. However, according to experimental results, the microfluidic system could be used effectively to evaluate the viscoelasticity of human blood when compared with conventional viscometers.

Conclusions
In this study, a viscoelasticity measurement method for human blood samples (normal blood sample, fixed blood sample) was suggested by quantifying the interface in coflowing streams when a blood sample was infused at a sinusoidal flow rate. Using a discrete fluidic circuit model, a firstorder nonlinear differential equation for the interface (αB) in coflowing streams was derived. Two approximation factors (F1, F2) were applied to convert a nonlinear term into a linear term. The viscosity and compliance (as viscoelasticity) were derived analytically from the linear differential equation. From numerical simulation, two approximate factors were obtained: F1 = 1.094 and F2 = 1.1087. The normalized difference between the nonlinear term and linear term was less than 10%. The experimental results showed that compliance varied linearly by period (T= 120, 240, 360, and 480 s). However, the blood viscosity remained constant with respect to period. The hematocrit and diluent contributed to varying viscoelasticity. Finally, viscoelasticity varied substantially depending on the degree in rigidity of RBCs. From the experimental results, it was found that the present method had the ability to monitor variations of viscoelasticity while the syringe pump was set to a pulsatile flow-rate pattern. As a limitation, the microfluidic device used in this study had microfluidic channels that Figure 7. Quantitative compliance of G B and C B with respect to previous study [18] and present study. Viscoelasticity (G) was represented as G = G 1 + j G 2 . Here, G 1 (storing modulus) and G 2 (loss modulus) were expressed as G 1 = G λ B 2 ω 2 1+λ B 2 ω 2 and G 2 = G λ B ω 1+λ B 2 ω 2 . Variations of G 1 and G 2 were represented with respect to radial frequency (ω). G 1 and G 2 tended to vary depending on λ B ·ω. However, as λ B ·ω showed a significant difference for both systems, it was apparent that both systems exhibited different variations of G 1 and G 2 with respect to ω. The microfluidic system and conventional viscometer showed different characteristics in terms of angular frequency (or period) and time constants. However, according to experimental results, the microfluidic system could be used effectively to evaluate the viscoelasticity of human blood when compared with conventional viscometers.

Conclusions
In this study, a viscoelasticity measurement method for human blood samples (normal blood sample, fixed blood sample) was suggested by quantifying the interface in coflowing streams when a blood sample was infused at a sinusoidal flow rate. Using a discrete fluidic circuit model, a first-order nonlinear differential equation for the interface (α B ) in coflowing streams was derived. Two approximation factors (F 1 , F 2 ) were applied to convert a nonlinear term into a linear term. The viscosity and compliance (as viscoelasticity) were derived analytically from the linear differential equation. From numerical simulation, two approximate factors were obtained: F 1 = 1.094 and F 2 = 1.1087. The normalized difference between the nonlinear term and linear term was less than 10%. The experimental results showed that compliance varied linearly by period (T= 120, 240, 360, and 480 s). However, the blood viscosity remained constant with respect to period. The hematocrit and diluent contributed to varying viscoelasticity. Finally, viscoelasticity varied substantially depending on the degree in rigidity of RBCs. From the experimental results, it was found that the present method had the ability to monitor variations of viscoelasticity while the syringe pump was set to a pulsatile flow-rate pattern. As a limitation, the microfluidic device used in this study had microfluidic channels that were rectangular shape. However, human or animal blood vessels are not rigid and possess circular shape. When compared with live blood vessels, the rectangular shape of microfluidic device might contribute to varying viscoelasticity of the microfluidic channels within the blood sample. Microfluidic channels with 90 • turns can cause live mammalian cells to lyse and aggregate in corners, causing fluid flow disruptions. In our next design, we will evaluate the contributions of channel shape and dimensions to viscoelasticity characteristics of the blood samples. In the near future, the present method will be employed to obtain the viscoelasticity of a blood sample circulated under ex vivo and in vivo conditions. In other words, to measure biomechanical properties of blood circulating under an extracorporeal bypass loop or hemodialysis, the blood sample will be collected from the fluidic circuit at periodic intervals. Biomechanical properties are then evaluated with individual conventional viscometers. However, the repetitive collection tends to reduce blood volume in the circuit substantially. If the present method will be integrated into the future fluidic circuit, viscoelasticity of the blood sample will be obtained at various intervals of time, even without periodic collection of blood samples.
Blood volume loss might be reduced considerably in our microfluidic device when compared with the conventional method. Finally, the data obtained by the novel microfluidic system could be employed effectively to monitor vascular diseases and human health.