• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Dynamics of magnetic microbubble transport in blood vessels

    2023-10-11 07:55:38JieChen陳杰ChenghuiWang王成會andRunyangMo莫潤陽
    Chinese Physics B 2023年9期
    關鍵詞:王成陳杰

    Jie Chen(陳杰), Chenghui Wang(王成會), and Runyang Mo(莫潤陽)

    Institute of Shaanxi Key Laboratory of Ultrasonics,Shaanxi Normal University,Xi’an 710119,China

    Keywords: magnetic microbubbles,targeted drug delivery,ultrasound,magnetic field

    1.Introduction

    Ultrasound combined with microbubble contrast agents has been widely demonstrated as a method of improving the delivery of drugs and genesin vitroandin vivo.However,the technique is limited by the relatively low delivery efficiencies achieved.Encapsulated microbubbles coupled with magnetic nanoparticles are referred to as magnetic microbubbles(MMBs), and are regarded as a recent development with the potential to greatly improve the efficacy of drug delivery and procedure sonoporation.[1–4]The shell elasticity allows volume oscillations in an acoustic field,which can be more easily controlled by an external magnetic field and demonstrate great potential as a drug delivery vehicle and multi-modal imaging tool.[5–7]

    Gene transfection enhanced by MMBs was first demonstrated by Strideet al.,[8]and then confirmed by Mulvanaet al.[9]and Vlaskouet al.[10]One possible mechanism is that both the concentration and proximity of the MMBs to the target cells are improved by the magnetic field.[11]However,this assumption needs to be further confirmed[12–14]and more experiments are needed.The effects of nanoparticles on the dynamic and acoustic responses of MMBs have been characterized.[15,16]The harmonic response and cavitation activity have been investigated experimentally.[17,18]Zhaoet al.[19]found experimentally that the magnetic susceptibility is proportional to the square of MMB radii.Meanwhile,MMBs also show powerful potential in tumor therapy and targeted drug delivery.[20,21]Both magnetic nanoparticles and nanodrugs can be released from oscillating MMBs.[12,14]Zhanget al.[2]reported an application in sonothrombolysis and confirmed that MMBs can be retained by a rotational magnetic field against blood flow.Recently, a series of sonodynamic therapies induced by oscillating MMBs in external fields were further verified.[22,23]Several specific devices were designed to maximize the superiority for tumor treatments.[13,24,25]

    The kinetics and regulation of MMBs are the basis of its application.To describe the volume oscillations exposed to ultrasound, Mulvanaet al.[16]conducted a preliminary modelling study.A fixed volume fraction of solid nanoparticles was used to examine the influence of an additional magnetic layer.Recently, Zhaoet al.[26]developed a multilayer magnetic shell model and the nonlinear oscillating characteristics were numerically evaluated.However,these studies have their limitations, because they ignore the influence of an applied magnetic field and non-spherical pulsation.Previous numerical studies on bubble dynamics in a magnetic field often were replaced by the dynamic of one or more non-magnetic bubbles in magnetic fluids.In this way, Lind[27]simulated the dynamics of non-spherical MMBs and showed that the magnetic susceptibility difference between bubbles and the ambient fluid can provide an additional means of flow control.Malvaret al.[28]formulated a new version of the Rayleigh–Plesset equation by considering the magnetic stress contribution and proved that the bubble oscillating motion can be controlled and its collapse avoided when a high magnetic field is applied.Chenet al.[29]expanded this work and found that the vibration amplitude could be enhanced by the magnetic field.Unfortunately, the coupling effect between vibration and translation has not been investigated and this incomplete understanding of the transport mechanism of MMBs makes some phenomena observed from experiments difficult to explain.

    Mathematical study acts an important role in exploring this mechanism by offering detailed mechanical processes that would be difficult or impossible to seein vitroorin vivo.In this work, theoretical models are established to simulate the dynamics of MMB transport in blood vessels under acoustic pressure and magnetic fields.In addition,some parameters associated with drug release are also simulated, especially the shear stress generated by bubble pulsation, which increases with the intensity of the applied magnetic field.The study attempts to provide insight into the details of MMB motion during targeted transport and support some experimental conclusions.

    2.Theory and method

    2.1.Translation in blood vessels

    MMBs transported in blood vessels are capable of targeted delivery and drug release when driven by ultrasound and magnetic fields.As shown in Fig.1, typical processes were summarized as translation, pulsation and drug release.MMBs can be retained in the tumor region under the guidance of a magnetic field, and perform drug release behavior in ultrasound.

    In a steady vessel flow, a single MMB is governed by three forces,expressed as

    wheremeffis the effective mass of the MMB including the added mass,and the three terms on the right side of Eq.(1)are the magnetic force,buoyancy and drag force,respectively.

    In addition,meffis small enough relative to other factors,which means the terminal velocity can be rapidly obtained.The process of MMB transport in the target region is shown in Fig.2.When MMBs arrive at the targeted region,FMis generated by a linear Halbach array of magnets,which is used in the experiment work of Barnsleyet al.[24]and expressed by

    wheremis the total magnetic moment of the magnetic nanoparticles,χbis the magnetic susceptibility of the bubbles,μ0is the permeability of vacuum space,andBis the magnetic induction.The expressionsχb=(6×10-6m)R2, dB/dz=900e(-0.75zmag/mm)T/mandB= 1.1e(-0.63zmag/mm)Tare cited from the literature.[19,24]

    The vessel is assumed to be a laminar cylindrical channel in thex-axis direction,[24]and the liquid flow ratevFis given by

    wherevAis the mean velocity andRCis the channel radius.

    Taking blood as a Newtonian fluid, the drag forceFDis followed by Stokes’law and given by viscous drag

    whereηLis the liquid viscosity,Ris the temporary radius andvis the velocity vector of the transport MMBs.Only the two dominant forcesFMandFDare considered in this model;one can obtain the result by combining Eqs.(1)–(4),given as

    IfFMis uniformly distributed in the local region,the component ofvinz-directionvzis approximated as a constant.This method has been verified by simulations and experiments in the literature.[24]

    Fig.1.Sketch of MMBs used for targeted drug delivery.

    Fig.2.Targeted delivery of MMBs in a flowing channel.The channel radius is RC with a flow in the x-direction, and each bubble is trapped in the target area marked by the shaded region from x0 to xM.Linear Halbach array magnets with thickness dmag are fixed at a distance zmag from the channel center.The dotted lines show the bubble trajectories under the attraction of the magnets.

    Primarily, our model is based on the one-dimensional cylindrical channel flow along thex-direction according to Eq.(3), and the magnetic force is applied in thez-direction.There is no driving force in they-direction, so the displacement in they-direction is ignored.Further, the coordinate componentsx(t) andz(t) are obtained by integrating the velocity components in Eq.(5).Taking the original coordinate position as (x0,y0,z0),x(t) andz(t) can be represented as follows:[24]

    Due to the strong magnetic force,MMBs are assumed to have no slippery contact with the channel wall.Thereby,the capture efficiency is expressed as[24,25]

    2.2.Dynamic model of coupled motion

    Previous studies have shown that magnetic nanoparticles have little effect on acoustic response, and MMBs are also not obviously affected by magnetic fields.[16]Donikove’s model[30]describes the vibration of a thin-shell bubble in an infinite liquid,represented as

    whereRandR0are the transient radius and the initial radius,respectively.˙Ris the bubble wall velocity, the liquid density and viscosity areρLandηL,respectively,σis the surface tension,Pg0is the equilibrium pressure,γis the polytrophic exponent,Pac(t)=Pasin(2π ft)is the acoustic driving pressure,andSrepresents the role of the encapsulated shell and is described by

    whereκ(˙R/R)=κ0+κ1˙R/Ris the shell surface viscosity (in kg·s-1) andζis the shell surface elasticity (in N·m-1).The values of shell parametersκ0,κ1andζare obtained from the literature,[30]and some observations[17,21,30]show that the stiffness and viscosity of the shell are related to magnetic particle concentration.When MMBs move in the vessel,[31,32]the kinetic energy of the system is contributed by pulsation and translation,which is defined as

    Considering the work done by pulsation and external forcesFexacting on MMBs, the potential energy[33]of the system is

    From the right terms of Eq.(8), ΔPincludes an additional magnetic pressurePM,andPMis given by Lind:[27,34]

    whereμ0is the vacuum permeability, andB0is the magnetic field strength and generally taken as 1–2 T.[8,16,27]The magnetic susceptibility difference Δχis a dimensionless value between the MMBs and the ambient fluid.Fexin Eq.(11) denotes the external forces acting on the MMBs,and can be written as

    The right terms of Eq.(14) correspond to Levich viscous drag,[35]magnetic force and secondary Bjerknes force,[36]and the bubble volume isVb.

    In addition,the Lagrangian quantity of the system is

    If other dissipative factors are ignored, the Lagrange equation is expressed as

    takingR,x,zand ˙R, ˙x, ˙zas generalized coordinates and velocities,respectively.This operation yields the equations of MMB dynamics

    Equations (17)–(19) are the coupled formulas describing the radial pulsation and translation of MMBs in the vessel.The contributions of the magnetic field and vessel wall are considered, such as the magnetic pressure, magnetic force, Bjerknes force,velocity gradient and the system kinetic energy adjusted by the vessel wall, which are the main factors in describing the boundary effect of fluid–solid interaction.This model has a similar form to the non-magnetic bubbles in a magnetic fluid,[28,29]and it also agrees with the reported models of gas bubbles in a channel and near a curved wall.[31,32]For more complex cases, the reflection and absorption of the vessel wall, even for non-spherical jets, remain to be further clarified.

    Assuming MMBs vibrating at small amplitudes, the eigenfrequency[37]is obtained by Eq.(17)as

    wheref*0andf0correspond to the vessel and infinite environment,respectively.

    2.3.Drug release dynamics

    Safely increasing transfection efficiencies in sonoporation(the amount of drug uptake by cells)remains an important issue and the shear stress caused by microstreaming has been demonstrated as an effective means.When a drug particle is detached from a MMB, the force balance relationship[12]is satisfied:

    The microstreaming, generated by vibrating bubbles when they approach the wall, can be described byψ(r,θ) in spherical coordinates.[38]The dimensionless velocityu*is a function ofθand given by Eq.(24).

    whereAis the dimensionless source strength depending on the oscillation mode of the bubble.To simplify the description,a dimensionless radial functionf(r)=is introduced to demonstrate the streamline distribution.

    According to Nyborg’s theory,[39,40]if the pulsating bubble is close to the vessel wall,thexcomponent of shear stress on the cells can be quantified by

    whereξmis the pulsating amplitude andδ=x/ld.

    3.Numerical simulation

    The parameters used in numerical analysis are as follows:R0=2 μm,P0=105Pa,ρL=103kg·m-3,σ=0.085 N·m-1,ηL= 10-3Pa·s,Pa= 60 kPa, Δχ= 0.1, andf= 1 MHz.Typical parameters of the shells of MMBs areκ0=1.45×10-8kg·s-1,κ1=6×10-15kg andζ=0.2 N·m-1.Moreover, the drug particles are simulated withdP= 100 nm,ρP=3×103kg·m-3andCA=0.5.

    3.1.Transport in gradient magnetic field

    3.1.1.Magnetic field analysis

    The strength and gradient of the magnetic field are essential for MMB localization.Figure 3(a) plots the fitting curve ofBas a function ofzmaggiven asB=1.1e(-0.63zmag/mm) T.It matches the experimental data[24](blue dots in Fig.3(a))very well.The gradient of the magnetic field depends onzmagand is expected to be dB/dz=900e(-0.75zmag/mm) T/m (inset plot in Fig.3(a)).The values ofB·dB/dzare 18 T2·m-1and 53.2 T2·m-1in the literature,[8,11]which agree with our data of 18.1 T2·m-1and 54.6 T2·m-1atzmag= 2.9 mm and 2.1 mm, respectively.The magnetic forceFMacting on MMBs is related tozmagin Fig.3(b).A semilog plot is presented to emphasize the exponential properties ofFMwithzmag.The marked data of 100 μm radius atzmag=4 mm is 2.142×10-7N, which agrees with experimental observations.[19]TheFMmarked as 8.568×10-11N for a radius of 2 μm is about 2500 times for 100 μm.This variation follows the ruleFM~R2due to magnetic nanoparticle numbers.When MMBs are attracted by magnets, their position and acting forces are changed with time.Larger bubbles carry more magnetic particles and produce stronger magnetic force.In Fig.3(c),a bubble with 2 μm radius moves fromz=3 mm toz=0 mm at the latest and takes 0.09 s,which is much longer than a larger one.This suggests that the size sorting of MMBs could be assisted by an external force,[41]and a better operation on aggregation would be realized by the attached magnetic force.Figure 3(d)shows the competition between forces on a bubble with 2 μm radius atz0=3 mm.When the bubble is close to the magnets, the velocity is gradually accelerated by the magnetic force,and the drag force is also increased until terminal equilibrium is reached.As expected, the values of both the added mass force and buoyancy force are smaller than magnetic and drag forces; this observation supports the simplification of Eq.(5).

    Fig.3.Responses of MMBs to the linear Halbach array magnets.(a) Character of gradient magnetic field; (b) relationship of FM ~zmag;(c)transport in z-direction guided by magnets;(d)forces acting on a bubble of R=2 μm at z0=3 mm.

    3.1.2.Transport trajectory

    In a vessel ofRC=0.25 mm,MMBs moving fromx0=0 mm toxM=2 mm along thexaxial direction are manipulated by the linear Halbach array magnets.For mean flow velocities ofvA= 20 mm·s-1and 40 mm·s-1, trajectories are shown in Figs.4(a) and 4(b), respectively, and each horizontal line corresponds to a starting position.As an example for comparison, the black line is retained atx= 2 mm for 40 mm·s-1andx=1 mm forvA=20 mm·s-1.Higher velocityvAproduces a horizontal displacement with a longer distance.Smaller distance to the magnets leads to better capture of MMBs with a shorter horizontal displacement, while higher flow velocity weakens the capture efficiency.The red line in Fig.4(c)is the fitting curve from Eq.(7)that matches the calculated results(statistics from a large number of tracks),indicating that capture efficiency decreases asvAincreases,and the capture behavior observed in experiments can also be quantified by Eq.(7).[24,25]In other words, Eq.(7) can serve as a simple method to predict the capture efficiency without calculating tracks.The distance to the magnets is also an essential parameter in MMB capture.The relationship between translation velocityvzand separation distancezmagis shown in Fig.4(d), Eq.(5) shows an exponential relationship between them,and the marked data represent the speed at several specific distances.Above all, the trajectories of MMBs can be established by the model shown in Eq.(6), and the capture efficiency between calculation and the analytical formula has been well verified.

    3.1.3.Capture efficiency

    MMBs can be localized at the target region when they are attracted by a gradient magnetic field.Capture efficiency is studied by Eq.(5) using the control variable method.These parameters include the component of transport velocity in thezdirectionvz, mean velocityvA, vessel radiusRCand magnet lengthxM;the results are shown in Figs.5(a)–5(d).WhenRC=0.25 mm,the efficiency increases with improvedvzandxM(Figs.5(a)and 5(b))and decreases withvA(Fig.5(c)),so a higher moving velocity and longer magnets can be more efficient,whereas a stronger flow velocity has the opposite effect.However, ifvAandxMare constant, smallerRCand largervzare more effective as shown in Fig.5(d).

    Fig.4.Transport properties of MMBs in a blood vessel.(a)and(b)The trajectories in the x–z plane at zmag=3.25 mm,vA=10 mm·s-1 and 20 mm·s-1,respectively;(c)capture efficiency related to vA;(d)vz determined by zmag.

    Fig.5.Capture efficiency of MMBs.(a) RC =0.25 mm, xM =2 mm; (b)vz =3 mm·s-1, RC =0.25 mm; (c)xM =2 mm, RC =0.25 mm;(d)vA=20 mm·s-1,xM=2 mm.

    3.2.Pulsation in uniform magnetic field

    3.2.1.Pulsation in unbound case

    To study MMBs pulsating in a uniform magnetic field,first, the unbound case without the vessel wall constraint was simulated.TheR–tcurves with code testing are shown in Fig.6(a); the parameters used here areR0= 1.4 μm,f= 3 MHz andPa= 100 kPa.The prediction of the proposed model (green line) matches well with the data (red circles) in the literature.[30]To highlight the effect of the shell,an unshelled bubble was used for comparison data(blue line).The shell effect reduces the amplitude and increases the harmonics.[42]Figure 6(b) shows the vibration of MMBs(R=2 μm,Pa=60 kPa andf=1 MHz) in magnetic fields ofB0=1 T,1.5 T and 2 T.The vibration amplitude increases withB0and one possible mechanism is that magnetic field promotes acoustic response.Similar results were found in Lind’s theoretical analysis of non-spherical bubbles in a magnetic liquid,[27]and it is also supported by violent bubble oscillation manipulated by external fields in a magnetic fluid[28,29]and electric medium.[43]Moreover, after removing the ultrasonic field at 6.5 μs, damping oscillations of MMBs are also displayed, and the final static radius is increased with the intensity of magnetic field.

    As shown in Fig.6(c), the oscillation amplitude near the resonance frequency increases significantly with magnetic field intensity.ConsideringB0=0 T,the maximum amplitude corresponding to the frequency is 2.6 MHz,which agrees with the 2.63 MHz calculated by Eq.(24) usingζ=0.2 N·m-1.Meanwhile,the resonant peak shifts from 2.6 MHz to 1.4 MHz after applyingB0=2 T, and the resonance frequency drops because the magnetic field increases the final static radius.Figure 6(d) depicts the curve of magnetic pressurePMfor Δχ=0.01, 0.05 and 0.1; the larger the susceptibility difference is, the fasterPMincreases withB0improvement.Conversely, bubbles in liquid oxygen with susceptibility 0.0035 have a smaller magnetic pressure of 1.4 kPa atB0=1 T,[34]which directly supports the idea that magnetic pressure is determined by the magnetic field and susceptibility difference between MMBs and the ambient fluid.

    Fig.6.MMBs vibrating in a uniform magnetic field.(a) Code testing and MMB analysis in B0 =0 T (f =3 MHz, Pa =100 kPa); (b) R–t curves adjusted by the magnetic field(f =1 MHz,Pa=60 kPa);(c)oscillating amplitude in B0=0 T,1 T and 2 T(Pa=60 kPa),Rmax (solid line)and Rmin (dashed line);(d)PM determined by B0 and Δχ.

    3.2.2.Resonance frequency

    First,the shell effect on resonance frequencyf0is shown as an unbound case in Fig.7(a).The hybrid shell has higher surface elasticity compared with an uncoated bubble and yields a higherf0.The prediction of the green line withζ=0.3 N·m-1aligns with the experimental data,[44]and the discrepancy decreases with increasing bubble size.Moreover,f0can also be used to evaluate the shell elasticity.[45]As shown in Figs.7(b) and 7(c), a largerRCand further distanceldyield a lower resonance frequency.Moreover,f*0/f0decreases greatly with increasingRand the observed rule is consistent with the literature,[46,47]indicating the reliability of this model.

    Fig.7.Resonance frequency of single MMBs.(a)Unbound;(b)in the vessel center;(c)close to the vessel wall.

    3.2.3.Pulsation near vessel wall

    Oscillating MMBs near the wall were used to verify the validity of the model.The red circles in Fig.8(a) represent the data carried by Garbin[48]and the green line is the predicted result of this model.The two results almost coincide and the deviation may be attributed to uncertainldor nonspherical pulsation.In contrast, a bubble far away from the wall is seen as a free one, shown as a blue line; a shorter period and larger amplitude imply that boundary constraints suppress oscillations.[49,50]As shown in Fig.8(b),the boundary constraint causes a resonance frequency decrease,and the magnetic field obviously alters the oscillatory response when compared with the unconstrained case as shown in Fig.6(c).

    Fig.8.Pulsation of MMBs near the vessel wall.(a) Compared with Garbin’s experimental data;[48] (b) oscillating amplitude in various magnetic fields,Rmax (solid line)and Rmin (dashed line).

    3.3.Coupled motion

    The coupled motion of MMBs near the vessel wall is shown in Fig.9(a).The three cases have similar trends.There is a much stronger response in the case of both fields than in the other cases,which could be ascribed to the synergistic manipulation.On the other hand, the pulsation in both fields is consistent with the non-spherical jet behavior near the vessel wall,[27]indicating that the volume vibration is improved,and then non-spherical deformation and jet generation.However,further studies on jet behavior go beyond our spherical hypothesis of MMBs.If the MMBs are far from the vessel wall andld=20 μm, Fig.9(b) displays the translational motion.The motion toward the wall is weaker in the acoustic field than in the magnetic field, which can be described by the Bjerknes force decays with (ld)2.The difference between the red and blue lines is attributed to the fact that the Levich drag in both fields is twice the Stokes drag in the magnetic field.Therefore,a rapid translation of the red line is observed at the initial stage and then the blue line approaches quickly as the MMBs approach the wall.This indicates that the translation of MMBs away from the vessel wall is mainly controlled by the magnetic force, while near the wall it is the Bjerknes force.Moreover,the combined application of magnetic and ultrasonic fields is shown with a shorter horizontal distancex(t).The capture of MMBs can be improved by the synergistic effect between magnetic force and Bjerknes force.

    Fig.9.Coupled motion of MMBs approaching the vessel wall.(a)Close to wall;(b)far from wall.

    3.4.Transport distance of drug particles

    MMBs carrying drugs and magnetic nanoparticles have the potential to improve localized drug delivery and tumor therapy.[11–13]It is assumed that the slight deformation of the vessel wall is replaced by a rigid wall.In Fig.10(a),the evolution of transport distance vs.time by solving Eq.(18) and Eq.(23)is compared with the reported trajectory.[12,14]Overall,the similar trend in both results confirms that the model can predict the drug release well.Figures 10(b),10(c),10(d)show the influences ofB0,PaandR0on the transport distance, respectively.A stronger magnetic field,higher acoustic pressure and suitable bubble radius provide a further transport distance and allow better drug release.The ultrasound triggers MMBs with violent vibrations and liquid microstreaming improves the drugs separated from MMBs.For the latter two graphics,the displacement of released drug particles corresponds well with Guillaume’s observations,[14]and the existing difference may be explained by the simplified velocity distribution and non-spherical oscillation.

    Fig.10.Transport distance of drug particles.(a)Compared with experimental data in the literature,[25] (b)effects of intensity B0,(c)effects of driving pressure at B0=0 T,(d)effects of MMBs radius at B0=0 T.

    3.5.Acoustic streaming

    Fig.11.Mechanical effect associated with vibrating MMBs near the vessel wall.(a) Streamline distribution; (b) spatial dimensionless velocity;(c)shear stress associated with MMB position(f =1 MHz);(d)shear stress determined by acoustic driving frequency f and B0 at x/R0=0.4,ld=R0.

    Acoustic streaming pattern and direction are usually related to the acoustic driving frequency and amplitude.Streamlines help clarify the trajectory of drug release.[14]Simulated results are shown in Fig.11(a).Compared to Marmottant’s pattern,[51]our calculation provides a surprisingly similar streamline distribution away from the bubble center.It is beneficial to understand the process of drug particle release from oscillating MMBs.The non-dimensional release velocityu*as a function of spatial angleθis shown in Fig.11(b).The release velocityu*is strongly dependent on the streaming direction and distribution;its maximum value is obtained atθ=0.Figure 11(c)is shear stress distribution on vessel wall when the bubble oscillates at a distanceld,the negative stress means an opposite direction.At the position ofx/R0=0.3,the variation in shear stress with driving frequency is shown in Fig.11(d);the maximum value corresponds to the resonance frequency being driven.This indicates that the combined magnetic and ultrasonic fields could improve shear stress,which may generate sonoporation to enhance transfection,[8,27]and the sonothrombolysis test implies that a better lysis rate is achieved by more cavitation nuclei retained in magnetic targeting region.[11]

    3.6.Therapy effect analysis

    As noted above, the behavior of MMBs driven by external fields was numerically analyzed, and the results are consistent with the experimental data in the literature.Experiments[8,11,12]reported that the differences in aggregated drug particles in the tumor area caused by external fields eventually lead to different therapeutic effects.[12,13]Lind[27]suggested that a stronger jet formed by the combined application of acoustic–magnetic fields and the generated highspeed collapse is most likely to maximize cell poration and gene/drug delivery.Our observations indicate that the magnitude of shear stress near the vessel wall was increased under applied external fields.In contrast, a uniform magnetic field is verified to weaken the oscillating gas bubble without coated magnetic particles, and is supported by inhibition of sonoluminescence[52–54]and reduction of the blood–brain barrier opening volume.[55]These facts could serve as contrary evidence to approve the model in this work; fiercer pulsation of MMBs is achieved by combined acoustic and magnetic fields.A higher shear stress caused oscillating MMBs in the combined acoustic and magnetic fields, which aligns with the results[11–13]in sonothrombolysis and tumor therapy.The magnetic force makes more MMBs stay in the target area and magnetic pressure improves the pulsation amplitude.The higher shear stress on the vessel wall is expected to enhance drug delivery efficiency.

    4.Conclusions

    An analytical model of single MMBs, which are transported in blood vessel and exposed to ultrasound and magnetic fields, was established.Our simulation includes translation,vibration and coupled motion of both.Vibration streaming,shear stress and their effects related to particle capture are also predicted.The results suggest that both transport and pulsation are well manipulated by the magnetic field due to the magnetic force and magnetic pressure.The effects of the vessel wall and external fields on the resonance frequency are further evaluated.Acoustic streaming near the vessel wall is demonstrated;some parameters related to drug release are quantified and consistent with the experimental data in the literature.In summary,a combined application of ultrasound and magnetic fields not only improves the capture of MMBs in the targeted region but also increases the shear stress near the vessel wall,which may enhance the permeability of the tissue,thus leading to better drug uptake.These results provide better understanding of the mechanism of MMBs in clinical applications and can be used to optimize the delivery strategy.

    Acknowledgement

    Project supported by the National Natural Science Foundation of China (Grant Nos.12074238, 11974232, and 11727813).

    猜你喜歡
    王成陳杰
    陳杰
    牡丹(2023年7期)2023-04-13 13:44:38
    封面攝影
    牡丹(2022年9期)2022-05-17 17:32:03
    Nonlinear oscillation characteristics of magnetic microbubbles under acoustic and magnetic fields
    Experimental study on plasma actuation characteristics of nanosecond pulsed dielectric barrier discharge
    陳杰
    奇妙的大自然
    只往壞處想
    偏偏只往壞處想
    上海故事(2018年9期)2018-09-30 02:57:16
    放風箏
    紀實
    河南電力(2017年7期)2017-11-30 07:09:16
    色哟哟·www| 午夜福利在线观看免费完整高清在| 久久久午夜欧美精品| 国产日韩欧美亚洲二区| 中文精品一卡2卡3卡4更新| 高清午夜精品一区二区三区| 纯流量卡能插随身wifi吗| 国产精品一区二区在线不卡| 全区人妻精品视频| 亚洲丝袜综合中文字幕| 成年免费大片在线观看| 午夜老司机福利剧场| 嫩草影院入口| 国产亚洲最大av| 我的老师免费观看完整版| 亚洲丝袜综合中文字幕| 亚洲精品成人av观看孕妇| 观看av在线不卡| 色婷婷久久久亚洲欧美| 自拍欧美九色日韩亚洲蝌蚪91 | 国产欧美亚洲国产| 伊人久久国产一区二区| 黄色视频在线播放观看不卡| 久久97久久精品| 国产国拍精品亚洲av在线观看| 久久99精品国语久久久| 日韩一本色道免费dvd| 亚洲经典国产精华液单| 不卡视频在线观看欧美| 日本黄色片子视频| 久久人人爽av亚洲精品天堂 | 欧美成人精品欧美一级黄| 久久精品久久久久久久性| 国产淫片久久久久久久久| 亚洲婷婷狠狠爱综合网| 天堂8中文在线网| 又粗又硬又长又爽又黄的视频| 亚洲精品久久久久久婷婷小说| 亚洲经典国产精华液单| 免费黄网站久久成人精品| 久久 成人 亚洲| 久久久久久伊人网av| 成人亚洲欧美一区二区av| 欧美老熟妇乱子伦牲交| 中文字幕人妻熟人妻熟丝袜美| 日韩一区二区三区影片| 色视频www国产| av在线播放精品| 最近最新中文字幕免费大全7| 国产 一区精品| 日韩一本色道免费dvd| 国产高潮美女av| 欧美最新免费一区二区三区| 亚洲欧美成人精品一区二区| 涩涩av久久男人的天堂| 如何舔出高潮| av国产久精品久网站免费入址| 在线观看免费高清a一片| 女人十人毛片免费观看3o分钟| 国产精品人妻久久久影院| 国产精品伦人一区二区| 国产老妇伦熟女老妇高清| 一区二区三区乱码不卡18| 欧美成人a在线观看| a级毛片免费高清观看在线播放| 丰满人妻一区二区三区视频av| 精品久久久久久电影网| 免费观看av网站的网址| 精品久久久噜噜| 国产精品久久久久久久久免| 欧美3d第一页| 中文字幕亚洲精品专区| 少妇 在线观看| 午夜福利影视在线免费观看| 精品少妇久久久久久888优播| 九色成人免费人妻av| 久久久精品94久久精品| a级毛片免费高清观看在线播放| 少妇人妻久久综合中文| 久久99热这里只有精品18| 久久99热这里只有精品18| 国产成人91sexporn| 婷婷色麻豆天堂久久| 亚洲精品亚洲一区二区| 亚洲成色77777| 欧美+日韩+精品| 国产欧美日韩精品一区二区| 精品99又大又爽又粗少妇毛片| 在线免费观看不下载黄p国产| 人妻系列 视频| 国精品久久久久久国模美| 国产精品三级大全| 婷婷色综合大香蕉| 制服丝袜香蕉在线| 黄色日韩在线| .国产精品久久| 欧美变态另类bdsm刘玥| 精品久久久噜噜| 国产精品久久久久久精品古装| 天天躁日日操中文字幕| 亚洲精品国产色婷婷电影| av播播在线观看一区| 欧美日韩一区二区视频在线观看视频在线| 国产精品人妻久久久久久| 亚洲精品乱码久久久v下载方式| 97热精品久久久久久| 久久午夜福利片| 精品久久国产蜜桃| 久久精品久久久久久噜噜老黄| 国内揄拍国产精品人妻在线| 黄色一级大片看看| 青春草视频在线免费观看| 尤物成人国产欧美一区二区三区| 在现免费观看毛片| 18禁在线播放成人免费| 欧美bdsm另类| 日韩欧美 国产精品| 免费大片黄手机在线观看| 亚洲国产毛片av蜜桃av| 国产乱人偷精品视频| 亚洲成人一二三区av| 亚洲四区av| 亚洲,一卡二卡三卡| 伦理电影免费视频| 久久久久久伊人网av| 在线观看免费日韩欧美大片 | 黄片wwwwww| 午夜福利在线在线| 国产精品久久久久久久电影| 欧美激情国产日韩精品一区| 在线观看人妻少妇| 欧美丝袜亚洲另类| 一区二区三区乱码不卡18| 午夜免费鲁丝| 国产综合精华液| av国产精品久久久久影院| 国产精品久久久久久av不卡| 国产高清国产精品国产三级 | 91精品国产九色| 亚洲精品中文字幕在线视频 | 午夜免费鲁丝| 黑人高潮一二区| 91久久精品国产一区二区三区| 亚洲国产欧美在线一区| 免费播放大片免费观看视频在线观看| 18禁在线播放成人免费| 麻豆成人午夜福利视频| 成年免费大片在线观看| 国产久久久一区二区三区| 国产免费又黄又爽又色| 91午夜精品亚洲一区二区三区| 男人舔奶头视频| 日韩av不卡免费在线播放| 日本vs欧美在线观看视频 | 国产 一区精品| 美女国产视频在线观看| 亚洲成人中文字幕在线播放| 欧美精品一区二区大全| 亚洲成人av在线免费| 亚洲色图av天堂| 欧美区成人在线视频| av在线观看视频网站免费| 久久99精品国语久久久| 国产精品成人在线| 在线观看一区二区三区| 午夜福利在线在线| 午夜福利在线在线| 中文欧美无线码| 观看av在线不卡| 国产一区亚洲一区在线观看| 最近中文字幕高清免费大全6| 插阴视频在线观看视频| 午夜福利影视在线免费观看| 干丝袜人妻中文字幕| 女的被弄到高潮叫床怎么办| 男女无遮挡免费网站观看| h日本视频在线播放| 亚洲av中文av极速乱| 国产精品一区二区性色av| 国产精品99久久久久久久久| 亚洲精品一区蜜桃| 老熟女久久久| 亚洲欧美日韩另类电影网站 | 人妻少妇偷人精品九色| 国内精品宾馆在线| a级毛色黄片| 久热这里只有精品99| 91狼人影院| 人人妻人人澡人人爽人人夜夜| 寂寞人妻少妇视频99o| 免费高清在线观看视频在线观看| 精品久久久久久久久亚洲| 免费观看av网站的网址| 我的老师免费观看完整版| 纵有疾风起免费观看全集完整版| 人体艺术视频欧美日本| 欧美三级亚洲精品| 亚洲av欧美aⅴ国产| 又大又黄又爽视频免费| 久久久国产一区二区| 久久久国产一区二区| 丰满少妇做爰视频| 97在线人人人人妻| 精品国产露脸久久av麻豆| 国产午夜精品久久久久久一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 搡老乐熟女国产| 热99国产精品久久久久久7| av在线老鸭窝| 伦理电影大哥的女人| 日韩av不卡免费在线播放| 国产精品秋霞免费鲁丝片| 久久久久久久久大av| 大码成人一级视频| 国产精品一区二区在线不卡| 国产精品一区二区在线不卡| 一本久久精品| 国产v大片淫在线免费观看| 国产乱人偷精品视频| 国产av一区二区精品久久 | 精品人妻一区二区三区麻豆| 如何舔出高潮| 美女视频免费永久观看网站| 久久久久性生活片| 成年人午夜在线观看视频| 国产在视频线精品| 国产免费视频播放在线视频| 18禁在线播放成人免费| 我要看日韩黄色一级片| 王馨瑶露胸无遮挡在线观看| 免费在线观看成人毛片| 男人狂女人下面高潮的视频| 中文在线观看免费www的网站| 熟妇人妻不卡中文字幕| 久久久久久久久久久丰满| 亚洲人成网站在线观看播放| 久久久久久久久久人人人人人人| 亚洲激情五月婷婷啪啪| 精品久久久久久久久亚洲| xxx大片免费视频| 高清av免费在线| 成人影院久久| 国产成人免费无遮挡视频| 女性被躁到高潮视频| 亚洲四区av| 国产精品久久久久久精品电影小说 | 日韩人妻高清精品专区| 亚洲一区二区三区欧美精品| 亚洲精品日韩av片在线观看| 亚洲av成人精品一二三区| 国产高清国产精品国产三级 | 黄色欧美视频在线观看| 最近中文字幕2019免费版| 国产免费福利视频在线观看| 国产一区二区在线观看日韩| 搡老乐熟女国产| 美女高潮的动态| 人人妻人人看人人澡| 熟女电影av网| 我的女老师完整版在线观看| 观看美女的网站| 99久久精品国产国产毛片| 免费看日本二区| 岛国毛片在线播放| 成人黄色视频免费在线看| 日日摸夜夜添夜夜爱| 欧美 日韩 精品 国产| 国产精品.久久久| 精品久久久久久电影网| 能在线免费看毛片的网站| 久久婷婷青草| 一区二区三区精品91| 美女内射精品一级片tv| 午夜免费观看性视频| 深夜a级毛片| 好男人视频免费观看在线| 免费人成在线观看视频色| 在线观看免费日韩欧美大片 | 久久国产精品大桥未久av | 国产真实伦视频高清在线观看| 少妇人妻 视频| 美女xxoo啪啪120秒动态图| 高清欧美精品videossex| 熟妇人妻不卡中文字幕| 熟女av电影| 亚洲精品aⅴ在线观看| www.色视频.com| 一二三四中文在线观看免费高清| 欧美3d第一页| 97超碰精品成人国产| 三级国产精品欧美在线观看| 美女国产视频在线观看| 嫩草影院新地址| 青春草国产在线视频| 成人综合一区亚洲| .国产精品久久| av国产免费在线观看| 国产精品精品国产色婷婷| 国产老妇伦熟女老妇高清| 亚洲国产精品国产精品| 久久精品国产亚洲av涩爱| freevideosex欧美| 五月玫瑰六月丁香| 美女内射精品一级片tv| 超碰97精品在线观看| 午夜视频国产福利| 国产av码专区亚洲av| 亚洲精品日韩av片在线观看| 欧美少妇被猛烈插入视频| 99久国产av精品国产电影| 久久精品国产a三级三级三级| 日韩,欧美,国产一区二区三区| 亚洲欧美一区二区三区黑人 | 亚洲美女搞黄在线观看| 国产精品国产av在线观看| 大话2 男鬼变身卡| 亚洲综合色惰| 日日摸夜夜添夜夜爱| 日韩一区二区三区影片| 80岁老熟妇乱子伦牲交| 国产成人精品久久久久久| 99热网站在线观看| 啦啦啦视频在线资源免费观看| 亚洲精品中文字幕在线视频 | 啦啦啦啦在线视频资源| 国产欧美亚洲国产| 国产一区二区三区av在线| 国国产精品蜜臀av免费| 色婷婷久久久亚洲欧美| 亚洲av在线观看美女高潮| 国产欧美亚洲国产| 中文欧美无线码| 成人18禁高潮啪啪吃奶动态图 | 国产大屁股一区二区在线视频| 国产精品一区二区三区四区免费观看| 天天躁夜夜躁狠狠久久av| 久久人人爽人人爽人人片va| 亚洲国产最新在线播放| 干丝袜人妻中文字幕| kizo精华| 国产乱人视频| 日本wwww免费看| 国产精品一区二区在线不卡| 下体分泌物呈黄色| 中文字幕久久专区| 午夜福利在线观看免费完整高清在| 日韩制服骚丝袜av| 狂野欧美激情性xxxx在线观看| 简卡轻食公司| 免费人成在线观看视频色| 插逼视频在线观看| 亚洲av成人精品一二三区| 91在线精品国自产拍蜜月| 国产乱来视频区| 国产精品国产av在线观看| av网站免费在线观看视频| 久久99热这里只有精品18| 国产av国产精品国产| 亚洲av综合色区一区| kizo精华| 自拍偷自拍亚洲精品老妇| 国产精品爽爽va在线观看网站| a 毛片基地| 亚洲第一av免费看| 99视频精品全部免费 在线| 22中文网久久字幕| 免费人妻精品一区二区三区视频| 啦啦啦在线观看免费高清www| 99久久精品国产国产毛片| 直男gayav资源| 伦理电影大哥的女人| 日韩中字成人| 18禁在线播放成人免费| 秋霞在线观看毛片| av专区在线播放| 日本-黄色视频高清免费观看| 精品一区二区免费观看| 国产黄片视频在线免费观看| 久久99热6这里只有精品| 赤兔流量卡办理| 午夜福利高清视频| 十分钟在线观看高清视频www | 精品人妻偷拍中文字幕| 色吧在线观看| 一级毛片我不卡| 亚洲国产精品国产精品| 黄色配什么色好看| 亚洲美女黄色视频免费看| 久久久久人妻精品一区果冻| 国产成人aa在线观看| 80岁老熟妇乱子伦牲交| 亚洲精品aⅴ在线观看| 欧美精品国产亚洲| 春色校园在线视频观看| 熟妇人妻不卡中文字幕| 国产在线免费精品| 免费观看a级毛片全部| 午夜免费鲁丝| 成人二区视频| 亚洲久久久国产精品| 国产真实伦视频高清在线观看| 偷拍熟女少妇极品色| 免费看日本二区| 国产精品一区二区三区四区免费观看| 精品久久久久久久久亚洲| 一区二区三区免费毛片| 午夜精品国产一区二区电影| 亚洲成色77777| 十分钟在线观看高清视频www | 看非洲黑人一级黄片| 久久国产亚洲av麻豆专区| 久久精品国产亚洲av天美| 黑丝袜美女国产一区| 99热网站在线观看| 日韩亚洲欧美综合| 日韩一区二区视频免费看| 纵有疾风起免费观看全集完整版| 97在线视频观看| 狂野欧美激情性bbbbbb| 中国美白少妇内射xxxbb| 五月伊人婷婷丁香| 国模一区二区三区四区视频| av在线老鸭窝| 午夜福利高清视频| 舔av片在线| 女的被弄到高潮叫床怎么办| 波野结衣二区三区在线| 国产伦精品一区二区三区视频9| 日韩欧美精品免费久久| 久久久a久久爽久久v久久| 永久免费av网站大全| 丰满乱子伦码专区| 亚洲av综合色区一区| 亚洲经典国产精华液单| 十分钟在线观看高清视频www | 爱豆传媒免费全集在线观看| 国产精品久久久久久精品电影小说 | 色综合色国产| 久久久久久久久大av| 大香蕉97超碰在线| tube8黄色片| 国产极品天堂在线| 一级毛片久久久久久久久女| 特大巨黑吊av在线直播| 国产伦在线观看视频一区| 美女主播在线视频| av不卡在线播放| 成人午夜精彩视频在线观看| 国产精品成人在线| 高清不卡的av网站| 高清在线视频一区二区三区| 日本黄大片高清| 男女边吃奶边做爰视频| 99热6这里只有精品| 99国产精品免费福利视频| 搡女人真爽免费视频火全软件| 亚洲国产最新在线播放| 黄色欧美视频在线观看| 欧美精品一区二区大全| 成人漫画全彩无遮挡| 中文字幕亚洲精品专区| 九九在线视频观看精品| 精品久久久久久久久av| 超碰97精品在线观看| 人人妻人人澡人人爽人人夜夜| 日韩不卡一区二区三区视频在线| 国产精品久久久久成人av| 国产伦精品一区二区三区四那| 国产精品无大码| 美女内射精品一级片tv| 国产色爽女视频免费观看| 国产欧美另类精品又又久久亚洲欧美| 欧美日韩国产mv在线观看视频 | 最新中文字幕久久久久| 日本一二三区视频观看| 99久久精品一区二区三区| 日本av免费视频播放| 一区二区三区免费毛片| 一级毛片aaaaaa免费看小| 亚洲精品日韩在线中文字幕| 国精品久久久久久国模美| tube8黄色片| 午夜日本视频在线| 我的老师免费观看完整版| 天堂8中文在线网| 超碰97精品在线观看| 最近中文字幕2019免费版| 男女边吃奶边做爰视频| 国产精品蜜桃在线观看| 亚洲av日韩在线播放| 超碰av人人做人人爽久久| 日日啪夜夜爽| 男女啪啪激烈高潮av片| 免费观看无遮挡的男女| 一级毛片aaaaaa免费看小| 美女主播在线视频| 欧美xxxx性猛交bbbb| 黄色配什么色好看| 夫妻午夜视频| 国产久久久一区二区三区| 国产乱人视频| 亚洲婷婷狠狠爱综合网| 亚洲精品一二三| 亚洲伊人久久精品综合| 国产精品一二三区在线看| 国产精品久久久久久久电影| 久久ye,这里只有精品| 久久精品国产a三级三级三级| 亚洲av国产av综合av卡| 我要看黄色一级片免费的| 99re6热这里在线精品视频| xxx大片免费视频| 亚洲av不卡在线观看| 亚洲av日韩在线播放| 亚洲色图综合在线观看| 国产在线男女| 美女国产视频在线观看| 亚洲av不卡在线观看| 久久精品人妻少妇| 亚洲成人手机| 中国美白少妇内射xxxbb| 亚洲综合色惰| 美女高潮的动态| 亚洲精品aⅴ在线观看| 午夜福利视频精品| 亚洲国产精品999| 国产国拍精品亚洲av在线观看| 国产欧美日韩一区二区三区在线 | 国产免费一区二区三区四区乱码| 国产在线男女| 十分钟在线观看高清视频www | 91久久精品国产一区二区三区| 国产日韩欧美亚洲二区| av免费观看日本| 久久久久性生活片| 亚洲欧美精品自产自拍| 亚洲aⅴ乱码一区二区在线播放| 少妇的逼水好多| 亚州av有码| 亚洲精品一区蜜桃| 观看av在线不卡| 秋霞在线观看毛片| 永久网站在线| videossex国产| 国产一区二区三区综合在线观看 | 一本—道久久a久久精品蜜桃钙片| 嘟嘟电影网在线观看| 人妻夜夜爽99麻豆av| 亚洲va在线va天堂va国产| 亚洲久久久国产精品| 大片电影免费在线观看免费| 人妻系列 视频| 永久免费av网站大全| 2021少妇久久久久久久久久久| 网址你懂的国产日韩在线| 免费大片18禁| 91精品一卡2卡3卡4卡| 777米奇影视久久| 日日啪夜夜爽| 国产伦精品一区二区三区视频9| 干丝袜人妻中文字幕| 男女啪啪激烈高潮av片| 国产精品久久久久久精品电影小说 | 亚洲真实伦在线观看| 97精品久久久久久久久久精品| 插逼视频在线观看| 亚洲欧美日韩卡通动漫| 亚洲精品一二三| 国产亚洲最大av| 国产日韩欧美在线精品| 国产成人免费观看mmmm| 久久久久久人妻| 99热国产这里只有精品6| 免费看日本二区| 2022亚洲国产成人精品| 亚洲精品成人av观看孕妇| 日韩欧美 国产精品| 啦啦啦中文免费视频观看日本| 国产伦理片在线播放av一区| 激情 狠狠 欧美| 国产在线男女| 一级a做视频免费观看| 在线观看人妻少妇| 国产成人91sexporn| 在线亚洲精品国产二区图片欧美 | 男人爽女人下面视频在线观看| 久久人妻熟女aⅴ| 国产在视频线精品| 大陆偷拍与自拍| videossex国产| 人体艺术视频欧美日本| 91久久精品国产一区二区成人| 91狼人影院| 欧美少妇被猛烈插入视频| 亚洲人成网站在线观看播放| 午夜激情久久久久久久| 一个人免费看片子| 亚洲国产精品专区欧美| av黄色大香蕉| 亚洲av男天堂| 中文字幕人妻熟人妻熟丝袜美| 99热网站在线观看| 精品人妻一区二区三区麻豆| 99久久精品热视频| 直男gayav资源| 久久久久久伊人网av| 精品亚洲乱码少妇综合久久| 国产伦精品一区二区三区视频9| 亚洲av在线观看美女高潮| 亚洲怡红院男人天堂| 高清毛片免费看| 十分钟在线观看高清视频www | 久久精品久久精品一区二区三区| 久久av网站| 国产综合精华液|