首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
2.
Although left ventricular assist devices (LVADs) have had success in supporting severe heart failure patients, thrombus formation within these devices still limits their long term use. Research has shown that thrombosis in the Penn State pulsatile LVAD, on a polyurethane blood sac, is largely a function of the underlying fluid mechanics and may be correlated to wall shear rates below 500 s(-1). Given the large range of heart rate and systolic durations employed, in vivo it is useful to study the fluid mechanics of pulsatile LVADs under these conditions. Particle image velocimetry (PIV) was used to capture planar flow in the pump body of a Penn State 50 cubic centimeters (cc) LVAD for heart rates of 75-150 bpm and respective systolic durations of 38-50%. Shear rates were calculated along the lower device wall with attention given to the uncertainty of the shear rate measurement as a function of pixel magnification. Spatial and temporal shear rate changes associated with data collection frequency were also investigated. The accuracy of the shear rate calculation improved by approximately 40% as the resolution increased from 35 to 12 μm/pixel. In addition, data collection in 10 ms, rather than 50 ms, intervals was found to be preferable. Increasing heart rate and systolic duration showed little change in wall shear rate patterns, with wall shear rate magnitude scaling by approximately the kinematic viscosity divided by the square of the average inlet velocity, which is essentially half the friction coefficient. Changes in in vivo operating conditions strongly influence wall shear rates within our device, and likely play a significant role in thrombus deposition. Refinement of PIV techniques at higher magnifications can be useful in moving towards better prediction of thrombosis in LVADs.  相似文献   

3.
In order to bridge the gap of existing artificial heart technology to the diverse needs of the patient population, we have been investigating the viability of a scaled-down design of the current 70 cc Penn State artificial heart. The issues of clot formation and hemolysis may become magnified within a 50 cc chamber compared to the existing 70 cc one. Particle image velocimetry (PIV) was employed to map the entire 50 cc Penn State artificial heart chamber. Flow fields constructed from PIV data indicate a rotational flow pattern that provides washout during diastole. In addition, shear rate maps were constructed for the inner walls of the heart chamber The lateral walls of the mitral and aortic ports experience high shear rates while the upper and bottom walls undergo low shear rates, with sufficiently long exposure times to potentially induce platelet activation or thrombus formation. In this study, we have demonstrated that PIV may adequately map the flow fields accurately in a reasonable amount of time. Therefore, the potential exists of employing PIV as a design tool.  相似文献   

4.
Fu Y  Kunz R  Wu J  Dong C 《PloS one》2012,7(2):e30721
Tumor cell adhesion to the endothelium under shear flow conditions is a critical step that results in circulation-mediated tumor metastasis. This study presents experimental and computational techniques for studying the local hydrodynamic environment around adherent cells and how local shear conditions affect cell-cell interactions on the endothelium in tumor cell adhesion. To study the local hydrodynamic profile around heterotypic adherent cells, a side-view flow chamber assay coupled with micro particle imaging velocimetry (μPIV) technique was developed, where interactions between leukocytes and tumor cells in the near-endothelial wall region and the local shear flow environment were characterized. Computational fluid dynamics (CFD) simulations were also used to obtain quantitative flow properties around those adherent cells. Results showed that cell dimension and relative cell-cell positions had strong influence on local shear rates. The velocity profile above leukocytes and tumor cells displayed very different patterns. Larger cell deformations led to less disturbance to the flow. Local shear rates above smaller cells were observed to be more affected by relative positions between two cells.  相似文献   

5.
Experimental models that mimic the flow conditions in microcapillaries have suggested that the local shear stresses and shear rates can mediate tumor cell and leukocyte arrest on the endothelium and subsequent sustained adhesion. However, further investigation has been limited by the lack of experimental models that allow quantitative measurement of the hydrodynamic environment over adherent cells. The purpose of this study was to develop a system capable of acquiring quantitative flow profiles over adherent cells. By combining the techniques of side-view imaging and particle image velocimetry (PIV), an in vitro model was constructed that is capable of obtaining quantitative flow data over cells adhering to the endothelium. The velocity over an adherent leukocyte was measured and the shear rate was calculated under low and high upstream wall shear. The microcapillary channel was modeled using computational fluid dynamics (CFD) and the calculated velocity profiles over cells under the low and high shear rates were compared to experimental results. The drag force applied to each cell by the fluid was then computed. This system provides a means for future study of the forces underlying adhesion by permitting characterization of the local hydrodynamic conditions over adherent cells.  相似文献   

6.
OBJECTIVE: To describe the local hemodynamics and pressure losses of crural bypass anastomoses using instantaneous velocity fields acquired by particle image velocimetry (PIV). METHODS: Silastic models of a Taylor patch, a Miller cuff and a femoro-crural patch prosthesis (FCPP) were attached to a circuit driven by a Berlin Heart, providing a pulsatile flow with an amplitude of 450 to 25 ml/min (mean 200 ml/min). An outflow resistance of 0.5 mmHg/ml/min (peripheral resistance units, PRU) was modeled using small silastic tubes providing a phase shift of -12 degrees between flow and pressure curves. The working fluid consisted of a glycerine/water mixture with a viscosity of 4 mPas. Hollow glass spheres with a mean size of 9-13 microm were used as tracer particles. Instantaneous velocity fields were obtained by means of PIV and shear rates as well as shear stresses were calculated. Triggered by the flowmeter signal, 10 measurements at 100 ms intervals per cardiac cycle were obtained. The pressures were measured on the inflow and at both distal outflows. The resulting mean pressure losses due to flow separation and distal fluid acceleration were calculated. RESULTS: Inside the Taylor patch anastomosis a large flow separation at the hood containing a clockwise rotating vortex was found. Additionally a smaller flow separation at the heel and a flow stagnation zone on the floor of the recipient artery were observed. Conversely, inside the Miller cuff a counterclockwise rotating vortex was seen inside a large heel flow separation. The FCPP also showed typical separation areas at the hood and heel of the anastomosis, although these were smaller compared to the other anastomoses. Inside the FCPP anastomosis no vortex creation was observed throughout the cardiac cycle. The mainstream velocities at the inlet levels were comparable for the three anastomoses. A significant fluid acceleration was present at the antegrade as well as the retrograde outlets of the Taylor and Miller cuff, while the fluid acceleration at the antegrade outflow of the FCPP was small, which was attributed to the end-to-end configuration of the antegrade FCPP leg. The calculated normalized antegrade and retrograde pressure losses for the Taylor form were 0.90 and 0.88, for the Miller cuff 0.89 and 0.86 and for the FCPP 0.94 and 0.86, respectively. The shear stresses inside the flow separations of the three anastomoses were significantly lower than normal wall shear stresses. High shear stress levels were found inside the transition zones between flow separation and high velocity mainstream. CONCLUSIONS: The flow pattern inside cuffed or funnel shaped anastomoses consists of large flow separation zones, which are thought to be associated with intimal hyperplasia development. In addition, fluid accelerations at the distal outlets result in pressure losses, which may contribute to impaired crural perfusion.  相似文献   

7.
This study examines the effect of aneurysmal wall elasticity on the structure of flow within an elastic aneurysm during pulsatile flow. We visualized flow structure in a model of an elastic saccular aneurysm located at the bifurcation of the anterior cerebral artery and extending to the anterior communicating artery, and measured changes in the diameter of the aneurysm wall during pulsatile flow using particle imaging velocimetry (PIV). We similarly measured these features during steady flow by PIV and found that dilation of the aneurysmal wall absorbed the dynamic energy within the aneurysm. Accordingly, aneurysm wall elasticity functions as a biocompatible reaction that relieves wall shear stress acting on the vascular wall during pulsatile flow, and should thus inhibit the development and rupture of an aneurysm.  相似文献   

8.
BACKGROUND: Computational fluid dynamics (CFD) simulations using medical-image-based anatomical vascular geometry are now gaining clinical relevance. This study aimed at validating the CFD methodology for studying cerebral aneurysms by using particle image velocimetry (PIV) measurements, with a focus on the effects of small geometric variations in aneurysm models on the flow dynamics obtained with CFD. METHOD OF APPROACH: An experimental phantom was fabricated out of silicone elastomer to best mimic a spherical aneurysm model. PIV measurements were obtained from the phantom and compared with the CFD results from an ideal spherical aneurysm model (S1). These measurements were also compared with CFD results, based on the geometry reconstructed from three-dimensional images of the experimental phantom. We further performed CFD analysis on two geometric variations, S2 and S3, of the phantom to investigate the effects of small geometric variations on the aneurysmal flow field. Results. We found poor agreement between the CFD results from the ideal spherical aneurysm model and the PIV measurements from the phantom, including inconsistent secondary flow patterns. The CFD results based on the actual phantom geometry, however, matched well with the PIV measurements. CFD of models S2 and S3 produced qualitatively similar flow fields to that of the phantom but quantitatively significant changes in key hemodynamic parameters such as vorticity, positive circulation, and wall shear stress. CONCLUSION: CFD simulation results can closely match experimental measurements as long as both are performed on the same model geometry. Small geometric variations on the aneurysm model can significantly alter the flow-field and key hemodynamic parameters. Since medical images are subjected to geometric uncertainties, image-based patient-specific CFD results must be carefully scrutinized before providing clinical feedback.  相似文献   

9.
This study employed particle image velocimetry (PIV) to validate a numerical model in a complementary approach to quantify hemodynamic factors in distal coronary anastomoses and to gain more insights on their relationship with anastomotic geometry. Instantaneous flow fields and wall shear stresses (WSS) were obtained from PIV measurement in a modified life-size silastic anastomosis model adapted from a conventional geometry by incorporating a smooth graft-artery transition. The results were compared with those predicted by a concurrent numerical model. The numerical method was then used to calculate cycle-averaged WSS (WSS(cyc)) and spatial wall shear stress gradient (SWSSG), two critical hemodynamic factors in the pathogenesis of intimal thickening (IT), to compare the conventional and modified geometries. Excellent qualitative agreement and satisfactory quantitative agreement with averaged normalized error in WSS between 0.8% and 8.9% were achieved between the PIV experiment and numerical model. Compared to the conventional geometry, the modified geometry produces a more uniform WSS(cyc) distribution eliminating both high and low WSS(cyc) around the toe, critical in avoiding IT. Peak SWSSG on the artery floor of the modified model is less than one-half that in the conventional case, and high SWSSG at the toe is eliminated. The validated numerical model is useful for modeling unsteady coronary anastomotic flows and elucidating the significance of geometry regulated hemodynamics. The results suggest the clinical relevance of constructing smooth graft-artery transition in distal coronary anastomoses to improve their hemodynamic performance.  相似文献   

10.
Blood flow patterns in the human left ventricle (LV) have shown relation to cardiac health. However, most studies in the literature are limited to a few patients and results are hard to generalize. This study aims to provide a new framework to generate more generalized insights into LV blood flow as a function of changes in anatomy and wall motion. In this framework, we studied the four-dimensional blood flow in LV via computational fluid dynamics (CFD) in conjunction with a statistical shape model (SSM), built from segmented LV shapes of 150 subjects. We validated results in an in-vitro dynamic phantom via time-resolved optical particle image velocimetry (PIV) measurements. This combination of CFD and the SSM may be useful for systematically assessing blood flow patterns in the LV as a function of varying anatomy and has the potential to provide valuable data for diagnosis of LV functionality.  相似文献   

11.
Endovascular coiling is an acceptable treatment of intracranial aneurysms, yet long term follow-ups suggest that endovascular coiling fails to achieve complete aneurysm occlusions particularly in wide-neck and giant aneurysms. Placing of a stentlike device across the aneurysm neck may be sufficient to occlude the aneurysm by promoting intra-aneurysmal thrombosis; however, conclusive evidence of its efficacy is still lacking. In this study, we investigate in vitro the efficacy of custom designed flow divertors that will be subsequently implanted in a large cohort of animals. The aim of this study is to provide a detailed database against which in vivo results can be analyzed. Six custom designed flow divertors were fabricated and tested in vitro. The design matrix included three different porosities (75%, 70%, and 65%). For each porosity, there were two divertors with one having a nominal pore density double than that of the other. To quantify efficacy, the divertors were implanted in a compliant elastomeric model of an elastase-induced aneurysm model in rabbit and intra-aneurysmal flow changes were evaluated by particle image velocimetry (PIV). PIV results indicate a marked reduction in intra-aneurysmal flow activity after divertor implantation in the innominate artery across the aneurysm neck. The mean hydrodynamic circulation after divertor implantation was reduced to 14% or less of the mean circulation in the control and the mean intra-aneurysmal kinetic energy was reduced to 29% or less of its value in the control. The intra-aneurysmal wall shear rate in this model is low and implantation of the flow divertor did not change the wall shear rate magnitude appreciably. This in vitro experiment evaluates the characteristics of local flow phenomena such as hydrodynamic circulation, kinetic energy, wall shear rate, perforator flow, and changes of these parameters as a result of implantation of stentlike flow divertors in an elastomeric replica of elastase-induced saccular aneurysm model in rabbit. These initial findings offer a database for evaluation of in vivo implantations of such devices in the animal model and help in further development of cerebral aneurysm bypass devices.  相似文献   

12.
Simulated pathline visualization of computed periodic blood flow patterns   总被引:4,自引:0,他引:4  
Improvements in computer hardware and software have made it possible to model pulsatile blood flow in realistic arterial geometries. Such studies produce enormous amounts of velocity data, which are often difficult to interpret and communicate using traditional contour and/or vector field plots. Inspired by in vitro flow visualization techniques such as particle image velocimetry (PIV), we describe a simple and effective method for visualizing periodic three-dimensional velocity data, based on the subdivision and sequential display of computed particle trajectories. Analogous to a PIV experiment, the length and spacing of such simulated particle pathlines are controlled by user-specified shutter-speed and frame rate variables. Strategies for color-coding pathlines to highlight important hemodynamic features such as recirculation zones and branch flow division are presented.  相似文献   

13.
Blood-perfused tissue models are joining the emerging field of tumor engineering because they provide new avenues for modulation of the tumor microenvironment and preclinical evaluation of the therapeutic potential of new treatments. The characterization of fluid flow parameters in such in-vitro perfused tissue models is a critical step towards better understanding and manipulating the tumor microenvironment. However, traditional optical flow measurement methods are inapplicable because of the opacity of blood and the thickness of the tissue sample. In order to overcome the limitations of optical method we demonstrate the feasibility of using phase-contrast x-ray imaging to perform microscale particle image velocimetry (PIV) measurements of flow in blood perfused hydrated tissue-representative microvessels. However, phase contrast x-ray images significantly depart from the traditional PIV image paradigm, as they have high intensity background, very low signal-to-noise ratio, and volume integration effects. Hence, in order to achieve accurate measurements special attention must be paid to the image processing and PIV cross-correlation methodologies. Therefore we develop and demonstrate a methodology that incorporates image preprocessing as well as advanced PIV cross-correlation methods to result in measured velocities within experimental uncertainty.  相似文献   

14.
This study is part of a FDA-sponsored project to evaluate the use and limitations of computational fluid dynamics (CFD) in assessing blood flow parameters related to medical device safety. In an interlaboratory study, fluid velocities and pressures were measured in a nozzle model to provide experimental validation for a companion round-robin CFD study. The simple benchmark nozzle model, which mimicked the flow fields in several medical devices, consisted of a gradual flow constriction, a narrow throat region, and a sudden expansion region where a fluid jet exited the center of the nozzle with recirculation zones near the model walls. Measurements of mean velocity and turbulent flow quantities were made in the benchmark device at three independent laboratories using particle image velocimetry (PIV). Flow measurements were performed over a range of nozzle throat Reynolds numbers (Re(throat)) from 500 to 6500, covering the laminar, transitional, and turbulent flow regimes. A standard operating procedure was developed for performing experiments under controlled temperature and flow conditions and for minimizing systematic errors during PIV image acquisition and processing. For laminar (Re(throat)=500) and turbulent flow conditions (Re(throat)≥3500), the velocities measured by the three laboratories were similar with an interlaboratory uncertainty of ~10% at most of the locations. However, for the transitional flow case (Re(throat)=2000), the uncertainty in the size and the velocity of the jet at the nozzle exit increased to ~60% and was very sensitive to the flow conditions. An error analysis showed that by minimizing the variability in the experimental parameters such as flow rate and fluid viscosity to less than 5% and by matching the inlet turbulence level between the laboratories, the uncertainties in the velocities of the transitional flow case could be reduced to ~15%. The experimental procedure and flow results from this interlaboratory study (available at http://fdacfd.nci.nih.gov) will be useful in validating CFD simulations of the benchmark nozzle model and in performing PIV studies on other medical device models.  相似文献   

15.
In most applications of biotechnology and downstream processing proteins are exposed to fluid stresses in various flow configurations which often lead to the formation of unwanted protein aggregates. In this paper we present physical degradation experiments for proteins under well-defined flow conditions in a four-roll apparatus. The flow field was characterized numerically by computational fluid dynamics (CFD) and experimentally by particle image velocimetry (PIV). The local shear strain rate as well as the local shear and elongation rate was used to characterize the hydrodynamic stress environment acting on the proteins. Lysozyme was used as a model protein and subjected to well-defined fluid stresses in high and low stress environment. By using in situ turbidity measurements during stressing the aggregate formation was monitored directly in the fluid flow. An increase in absorbance at 350 nm was attributed to a higher content of visible particles (>1 μm). In addition to lysozyme, the formation of aggregates was confirmed for two larger proteins (bovine serum albumin and alcohol dehydrogenase). Thus, the presented experimental setup is a helpful tool to monitor flow-induced protein aggregation with high reproducibility. For instance, screening experiments for formulation development of biopharmaceuticals for fill and finish operations can be performed in the lab-scale in a short time-period if the stress distributions in the application are transferred and applied in the four-roll mill.  相似文献   

16.
The transport of mass, momentum, and energy in fluid flows is ultimately determined by spatiotemporal distributions of the fluid velocity field.1 Consequently, a prerequisite for understanding, predicting, and controlling fluid flows is the capability to measure the velocity field with adequate spatial and temporal resolution.2 For velocity measurements in optically opaque fluids or through optically opaque geometries, echo particle image velocimetry (EPIV) is an attractive diagnostic technique to generate "instantaneous" two-dimensional fields of velocity.3,4,5,6 In this paper, the operating protocol for an EPIV system built by integrating a commercial medical ultrasound machine7 with a PC running commercial particle image velocimetry (PIV) software8 is described, and validation measurements in Hagen-Poiseuille (i.e., laminar pipe) flow are reported.For the EPIV measurements, a phased array probe connected to the medical ultrasound machine is used to generate a two-dimensional ultrasound image by pulsing the piezoelectric probe elements at different times. Each probe element transmits an ultrasound pulse into the fluid, and tracer particles in the fluid (either naturally occurring or seeded) reflect ultrasound echoes back to the probe where they are recorded. The amplitude of the reflected ultrasound waves and their time delay relative to transmission are used to create what is known as B-mode (brightness mode) two-dimensional ultrasound images. Specifically, the time delay is used to determine the position of the scatterer in the fluid and the amplitude is used to assign intensity to the scatterer. The time required to obtain a single B-mode image, t, is determined by the time it take to pulse all the elements of the phased array probe. For acquiring multiple B-mode images, the frame rate of the system in frames per second (fps) = 1/δt. (See 9 for a review of ultrasound imaging.)For a typical EPIV experiment, the frame rate is between 20-60 fps, depending on flow conditions, and 100-1000 B-mode images of the spatial distribution of the tracer particles in the flow are acquired. Once acquired, the B-mode ultrasound images are transmitted via an ethernet connection to the PC running the PIV commercial software. Using the PIV software, tracer particle displacement fields, D(x,y)[pixels], (where x and y denote horizontal and vertical spatial position in the ultrasound image, respectively) are acquired by applying cross correlation algorithms to successive ultrasound B-mode images.10 The velocity fields, u(x,y)[m/s], are determined from the displacements fields, knowing the time step between image pairs, ΔT[s], and the image magnification, M[meter/pixel], i.e., u(x,y) = MD(x,y)/ΔT. The time step between images ΔT = 1/fps + D(x,y)/B, where B[pixels/s] is the time it takes for the ultrasound probe to sweep across the image width. In the present study, M = 77[μm/pixel], fps = 49.5[1/s], and B = 25,047[pixels/s]. Once acquired, the velocity fields can be analyzed to compute flow quantities of interest.  相似文献   

17.
Emphysema is a progressive lung disease that involves permanent destruction of the alveolar walls. Fluid mechanics in the pulmonary region and how they are altered with the presence of emphysema are not well understood. Much of our understanding of the flow fields occurring in the healthy pulmonary region is based on idealized geometries, and little attention has been paid to emphysemic geometries. The goal of this research was to utilize actual replica lung geometries to gain a better understanding of the mechanisms that govern fluid motion and particle transport in the most distal regions of the lung and to compare the differences that exist between healthy and emphysematous lungs. Excised human healthy and emphysemic lungs were cast, scanned, graphically reconstructed, and used to fabricate clear, hollow, compliant models. Three dimensional flow fields were obtained experimentally using stereoscopic particle image velocimetry techniques for healthy and emphysematic breathing conditions. Measured alveolar velocities ranged over two orders of magnitude from the duct entrance to the wall in both models. Recirculating flow was not found in either the healthy or the emphysematic model, while the average flow rate was three times larger in emphysema as compared to healthy. Diffusion dominated particle flow, which is characteristic in the pulmonary region of the healthy lung, was not seen for emphysema, except for very small particle sizes. Flow speeds dissipated quickly in the healthy lung (60% reduction in 0.25 mm) but not in the emphysematic lung (only 8% reduction 0.25 mm). Alveolar ventilation per unit volume was 30% smaller in emphysema compared to healthy. Destruction of the alveolar walls in emphysema leads to significant differences in flow fields between the healthy and emphysemic lung. Models based on replica geometry provide a useful means to quantify these differences and could ultimately improve our understanding of disease progression.  相似文献   

18.
Thoracic endovascular aortic repair (TEVAR) has been introduced as a less invasive approach to the treatment of thoracic aortic aneurysm (TAA). However, the effectiveness of TEVAR in the treatment of TAA is often limited due to the complex anatomy of aortic arch. Flow preservation at the three supra-aortic branches further increases the overall technical difficulty. This study proposes a novel stent graft design with slit perforations that can positively alter the hemodynamics at the aortic arch while maintaining blood flow to supra-aortic branches. We carried out a computational fluid dynamic (CFD) analysis to evaluate flow characteristics near stented aortic arch in simplified TAA models, followed by in-vitro experiments using particle image velocimetry (PIV) in a mock circulatory loop. The hemodynamics result was studied in terms of time-averaged wall shear stress (TAWSS), oscillating shear index (OSI), and endothelial cell action potential (ECAP). The results showed that the stent graft with slit perforations can reduce the disturbed flow region considerably. Furthermore, the effect of the slits on flow preservation to the supra-aortic branches was simulated and compared with experimental results. The effectiveness of the stent graft with slit perforations in preserving flow to the branches was demonstrated by both simulated and experimental results. Low TAWSS and elevated ECAP were observed in the aortic arch aneurysm after the placement of the stent graft with slits, implying the potential of thrombus formation in the aneurysm. On the other hand, the effects of the stent grafts with full-slit design and half-slit design on the shear stress did not differ significantly. The present analysis indicated that not only could the stent graft with slit perforations shield the aneurysm from rupture, but also it resulted in a favorable environment for thrombus that can contribute to the shrinkage of the aneurysm.  相似文献   

19.
It is desired to understand the effect of alginic acid sodium salt from brown algae (alginate) as a viscosity modifier on the behavior of blood in vitro using a micro-particle image velocimetry (µPIV) system. The effect of alginate on the shape of the velocity profile, the flow rate and the maximum velocity achieved in rectangular microchannels channels are measured. The channels were constructed of polydimethylsiloxane (PDMS), a biocompatible silicone. Porcine blood cells suspended in saline was used as the working fluid at twenty percent hematocrit (H = 20). While alginate was only found to have minimal effect on the maximum velocity and the flow rate achieved, it was found to significantly affect the shear rate at the wall by between eight to a hundred percent.  相似文献   

20.
Velocity measurements were performed for the flow field generated by tethered krill Euphausia pacifica. The particle image velocimetry (PIV) technique was used to measure the velocity field in vertical planes aligned with the krill body axis. The krill generates a narrow jet-like flow behind and below the pleopods (roughly 25° below horizontal). The volume of fluid moving at greater than 10% of the maximum velocity near the pleopods is roughly 18 times larger than the volume of the krill. Thus, the hydrodynamic disturbance occupies a significantly larger region than the animal body. Other krill, sensing the flow disturbance, may take advantage of the flow induced by a neighbor to locate a mate or to draft for efficient propulsion.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号