高强度聚焦超声间歇式治疗中焦域温度分布的 仿真研究*

探究不同治疗深度、组织类型和治疗模式对高强度聚焦超声焦域温度场的影响。采用有限元法建立高强度聚焦超声辐照组织的二维轴对称仿真模型,通过Westervelt方程和Pennes生物热传导方程计算高强度聚焦超声焦域的声场和温度场分布。仿真结果表明:随着治疗深度的增加,焦域内温度逐渐降低,有效温升面积减小;不同组织类型在相同条件下的焦域温度不同,但焦域形状差别不大;单次治疗时间长、治疗时间间隔短的模式焦域温升速率快,有效温升面积区域大。焦域温度是高强度聚焦超声肿瘤治疗中判断治疗有效性和安全性的重要因素之一,该文通过数值仿真法,实现了预测高强度聚焦超声间歇式治疗模式下焦域组织温度场分布,有望为制定安全、有效的高强度聚焦超声术前治疗方案提供理论参考。     

剪切波对HIFU经颅聚焦形成温度场影响的数值仿真研究

在高强度聚焦超声经颅治疗时,既有纵波又有剪切波,为了保障该治疗方法的安全有效性,有必要分析剪切波对HIFU治疗温度场的影响。该文基于人体头颅CT数据和曲率半径为150 mm的256阵元的半球相控换能器建立三维高强度聚焦超声经颅声波传播模型,利用时域有限差分法结合Westervelt声波非线性传播方程、动量方程、质量守恒方程和Pennes生物热传导方程数值仿真其形成温度场,研究在相同输入功率、不同聚焦角度条件下对应阵元数进行激励时,剪切波对换能器形成温度场的影响。结果表明,随换能器聚焦角度减小,在几何焦点处形成的焦域面积逐渐增大,考虑剪切波形成的温度场达到65?C所需时间逐渐延长,焦点前移程度越大;在相同聚焦角度条件下,考虑剪切波的温度场达到65?C所需时间更短,旁瓣更少,在颅骨处的温度更高,对焦点前移几乎没有影响;随换能器聚焦角度减小,考虑剪切波的模型形成的焦域面积变化范围更大;幂指数函数形式对不同聚焦角度下焦域面积大小的拟合优度高,可预测不同聚焦角度换能器形成的焦域面积。     

基于电阻抗层析成像的高强度聚焦超声温度监测技术

作为一种对正常组织无损伤且不易引起癌细胞转移的非入侵肿瘤治疗手段,高强度聚焦超声(HIFU)治疗过程中焦域的温度监测是实现剂量精准控制的关键.本文基于生物组织的温度-电阻抗的关系,将电阻抗层析成像(EIT)和HIFU治疗相结合,提出了一种利用组织焦平面的表面电压实现电阻抗重构的检测技术.建立了HIFU治疗和EIT综合系统模型,在考虑组织的声吸收条件下,对三维Helmholtz方程在柱坐标下的声场计算进行了二维简化,并引入Pennes生物热传导方程来计算HIFU焦域的声压和温升分布特性;引入生物组织的温度-电阻抗关系,基于麦克斯韦电磁场理论,建立了具有温度分布HIFU焦域的电流和电压计算模型,利用恒流注入的边界条件实现电场计算,获得焦平面的表面电压分布.在数值计算中,利用实验聚焦换能器参数,模拟了在固定声功率下组织焦域的声场和温度场分布,以及中心和偏心聚焦条件下不同治疗时刻的电导率分布;然后通过对称电极的循环电流注入,计算了组织模型焦平面内的电流密度和电势分布,获得了焦平面圆周分布的表面电极电压;进一步采用修正的牛顿-拉夫逊算法,利用32×32的表面电极电压实现了焦平面内电导率分布的重建.结果表明,基于温度-电阻抗关系的EIT电导率重建技术不但能准确定位HIFU焦域中心,还能恢复HIFU治疗中焦域的温度分布,证明了EIT用于HIFU治疗中温度监测的可行性,为其疗效评估和剂量控制提供了一种无创电阻抗测量和成像新方法.     

基于换能器驱动信号特征的高强度聚焦超声焦域损伤实时监测

高强度聚焦超声(high intensity focused ultrasound,HIFU)焦域的实时监测是聚焦超声临床治疗面临的关键问题,目前临床常采用B超图像强回声的变化实现焦域组织损伤的监测,而B超图像出现的强回声大多与焦域处的空化及沸腾气泡相关,无法准确、实时地监测治疗状态.HIFU治疗中焦域组织会伴随温度升高、空化、沸腾和组织特性等变化,换能器表面的声学负载也在不断变化,针对该问题,本文构建了换能器电压电流实时检测平台,通过测量换能器电学参数来感知焦域组织的状态变化.以离体牛肝组织作为HIFU辐照对象,并将相位差变化的结果与离体牛肝组织损伤的结果进行了对照,实验结果表明,在HIFU辐照过程中,换能器电压与电流的相位会出现由相对平稳到大幅波动的过程,此时停止辐照可见焦域出现明显损伤,而此时B超图像灰度无明显变化;此外,当焦域出现空化时,其波动幅度与范围将较之更大.此方法可为HIFU焦域组织损伤监测提供一种新的研究方案和手段.     

HIFU治疗中换能器驱动电功率变化机制及规律*

高强度聚焦超声(HIFU)治疗中的驱动电功率对治疗效率起着非常关键的作用,驱动电功率控制的精准性势必会影响治疗的效率和安全性。前期研究表明：HIFU治疗过程中焦域瞬态物理特性的变化会导致换能器的负载阻抗发生变化,进而影响HIFU驱动电功率,但驱动电功率与焦域瞬态物理特性之间的影响关系及规律尚不明确。该文基于电压、电流传感器、空化检测探头和温度传感器等器件,构建了一种HIFU治疗中驱动电功率实时监测及焦域声空化、温度检测系统。基于该实验研究系统,以离体牛心组织作为HIFU辐照对象,分别研究了HIFU焦域温度变化、声空化及组织损伤与驱动电功率之间的变化关系及规律。研究结果表明：当焦域温度升高时,驱动电功率缓慢上升,驱动电功率与温度变化有良好的相关性;当空化产生时,驱动电功率出现明显的波动;当组织出现损伤时,驱动电功率呈陡然下降的变化。三种情景下,驱动电功率变化有明显区别,这有望为区分HIFU治疗过程中焦域处发生损伤和空化以及实时监测靶组织损伤程度提供一种新的解决方案。     

HIFU温度场可视化测量初探

本实验用蓖麻油模拟人体环境、热敏打印纸作为记录载体,测量高强度聚焦超声(HIFU)治疗系统焦点附近的温度分布情况。实验中热敏打印纸显影所指示的温度已经达到HIFU治疗所需要的温度,实验的结果对进一步研究HIFU温度场分布有借鉴作用。     

微泡对高强度聚焦超声焦域影响的研究

微泡对高强度聚焦超声(HIFU)治疗具有增效作用,而HIFU治疗中不同声学条件下微泡对HIFU治疗焦域的影响尚不清楚。本文基于声传播方程、Yang-Church气泡运动方程、生物热传导方程、时域有限差分法(FDTD)、龙格-库塔(RK)法数值仿真研究输入功率、激励频率和气泡初始半径对HIFU在含气泡体模中形成焦域的影响,并利用含Sono Vue造影剂的仿组织体模研究进行实验验证。结果表明,增大输入功率、气泡初始半径和升高激励频率均可增大焦域,随着输入功率的增大,焦域形状可能发生变化,而随着激励频率升高和气泡初始半径的增大,焦域会向远离换能器的方向移动。     

连续激光辐照锗材料损伤的数值模拟研究

 以热传导理论和热弹应力理论为基础,建立了连续平顶激光辐照锗材料的二维非稳态物理模型。对受辐照锗材料的热-力场分布进行数值模拟,比较了激光辐照时间对靶材热-力场分布的影响,探讨了温度场及应力场空间梯度随半径的分布以及激光辐照时间对热-力场空间梯度的影响。计算结果表明,靶材中心温升效应最大,而应力最大值出现在光斑半径外侧;光斑边缘处温度场和应力场梯度绝对值均达到最大,随着辐照时间的延长,材料同一点温度场和应力场梯度均增大。在给定条件下,计算了靶材的损伤时间,应力损伤域值最先达到。     

1 319 nm连续激光辐照皮肤组织的光热效应实验研究

实验测定了离体人皮肤组织在1 319 nm连续Nd:YAG激光辐照下径向和轴向不同位置的温度分布,研究了激光辐照功率密度和作用时间对皮肤组织温升分布的影响.结果表明,皮肤组织中的温度分布与探测位置、激光辐照功率密度以及作用时间密切相关;激光输出功率密度15.80 W/cm2时安全辐照时间阈值为120 s,功率密度48.92 W/cm2时安全辐照时间阈值为6 s.     

微泡对高强度聚焦超声声压场影响的仿真研究*

微泡对高强度聚焦超声(HIFU)治疗焦域具有增效作用,而HIFU治疗中不同声学条件下微泡对HIFU形成声压场的影响尚不清楚。本文基于气液混合声波传播方程、Keller气泡运动方程、时域有限差分(FDTD)法和龙格-库塔(RK)法数值仿真研究输入声压、激励频率、气泡初始空隙率和气泡初始半径对HIFU形成声压场的影响。研究结果表明,随着输入声压的增大,焦点处声压升高但焦点处最大声压与输入声压的比值减小,焦点位置几乎不变;随着激励频率和气泡初始半径的增大,焦点处声压升高且焦点位置向远离换能器方向移动;随着气泡初始空隙率的增大,焦点处声压降低且焦点位置向换能器方向移动。     

Chalcogenide Fibers: An Overview of Selected Applications

Chalcogenide glass fibers have many unique properties that make them attractive for applications that are not possible with silica fibers. Selected applications of chalcogenide fibers in research, commerce, and the military are reviewed.     

Some mathematical aspects of Anderson localization: boundary effect,multimodality, and bifurcation

Anderson localization is a famous wave phenomenon that describes the absence of diffusion of waves in a disordered medium. Here we generalize the landscape theory of Anderson localization to general elliptic operators and complex boundary conditions using a probabilistic approach, and further investigate some mathematical aspects of Anderson localization that are rarely discussed before. First, we observe that under the Neumann boundary condition, the low energy quantum states are localized on the boundary of the domain with high probability. We provide a detailed explanation of this phenomenon using the concept of extended subregions and obtain an analytical expression of this probability in the one-dimensional case. Second, we find that the quantum states may be localized in multiple different subregions with high probability in the one-dimensional case and we derive an explicit expression of this probability for various boundary conditions. Finally, we examine a bifurcation phenomenon of the localization subregion as the strength of disorder varies. The critical threshold of bifurcation is analytically computed based on a toy model and the dependence of the critical threshold on model parameters is analyzed.     

辣椒碱与β-环糊精包合物的核磁共振研究

利用核磁共振方法研究了辣椒碱β-环糊精包合物的化学计量比、空间结构信息及其在水溶液中的自扩散系数. 通过测定不同浓度比的辣椒碱和β-环糊精混合溶液的 ¹H NMR数据,绘制Job's曲线,辣椒碱和β-环糊精的Job's曲线均在r=0.5处出现拐点. 同时测定了该包合物的2D ROESY和DOSY谱图,ROESY谱图中NOE交叉信号出现在辣椒碱的H-1~H-8和β-环糊精的H-3′、H-5′、H-6′之间,DOSY测得β-环糊精和辣椒碱形成包合物前后的表观自扩散系数. 结果表明,辣椒碱β-环糊精包合物的主客体分子的化学计量比为1∶1,辣椒碱分子的异丙基端从β-环糊精的宽口端进入疏水腔,其中H-1~H-8部分在空腔内部,包合物的自扩散系数为2.95×10^-10m²/s.     

The ρ-π Puzzle of J/φ and φ′ Decays

The recent BES Collaboration data on φ →PV, particularly the isospin violating mode φ′ → π^oω^o and finding of a finite number for B(φ′ → K^*OK^O), enable us now to deal more precisely about the challenges to theory concerning this extraordinary and remarkableso-called ρ-π puzzle of J/φ and φ′ decays. In terms of the existing data and deploying the simplest phenomenology, measurement of φ′ → π⁺π^- and whether a finite number for the K^*+K^- mode might require a significantly large accumulation of data remain interesting questions.     

Interaction of Wave Trains with Defects

The evolution and transition of planar wave trains propagating through defects(obstacles) in an excitable medium are studied. When the frequency of the planar wave trains is increased, three different dynamical regimes,namely fusion, "V" waves, and spiral waves, are observed in turn and the underlying mechanism is discussed. The dynamics is concerned with the shapes of the defects. Circle, triangle, and rectangle defects with different sizes are considered. The increase of pacing frequency broadens the fan-shaped broken region in the behind of a rectangle defect.The increase of width of a triangle defect leads to breakup of wave trains easier while the change of height shows opposite effect, which is presented in a phase diagram. Dynamical comparison on defects with different shapes indicates that the decrease of the defect width along the propagation of wave trains makes the fan-shaped region and the minimal frequency for breakup of spiral both increased.     

环位形下鱼骨模驱动快离子输运的数值模拟研究

考虑包含动理学效应的鱼骨模结构,使用导心轨道程序 ORBIT,在磁面坐标下研究了不同的扰动模 幅度、频率对快离子再分布的影响,并分析了粒子与扰动发生共振的条件。模拟得出,鱼骨模扰动会使快离子在 实空间、相空间中发生再分布,芯部( ψ_p/ψ _w≤0.2 ,ψ_p 为极向磁通,ψ_w 为最后一个闭合磁面的磁通)快离子密 度下降约 20%,中间磁面位置上( ψ_p/ψ _w≥0.2,ψ_p/ψ _w≤ 0.6)的快离子密度增加约 7%;通过扫描频率发现,相空间中快 离子的再分布对模式频率敏感,并分析了快离子与鱼骨模扰动共振的条件。     

Darboux Transformation and N-soliton Solution for Extended Form of Modified Kadomtsev-Petviashvili Equation with Variable-Coefficient

The extended form of modified Kadomtsev-Petviashvili equation with variable-coefficient is investigated in the framework of Painlevé analysis. The Lax pairs are obtained by analysing two Painlevé branches of this equation. Starting with the Lax pair, the N-times Darboux transformation is constructed and the N-soliton solution formula is given, which contains 2n free parameters and two arbitrary functions. Furthermore, with different combinations of the parameters, several types of soliton solutions are calculated from the first order to the third order. The regularity conditions are discussed in order to avoid the singularity of the solutions. Moreover, we construct the generalized Darboux transformation matrix by considering a special limiting process and find a rational-type solution for this equation.     

HL-2M 装置时序控制与数据采集系统设计

基于 FPGA 和 IEEE1588PTPv2 协议设计了 HL-2M 装置时序控制系统,用于在精准时刻为测控系统 提供触发。数据采集系统采用客户端/服务端双模式管理。在 HL-2M 装置初始等离子体放电实验中测试结果表明, 时序控制和采集系统具有纳秒级精度触发,实时性强,数据传输速度快,稳定可靠,方便部署和管理的特点,满 足 HL-2M 等离子体放电实验的应用需求。     

Thermodynamics of 4D dS/AdS Gauss–Bonnet black holes according to consistent gravity theory in the presence of a cloud of strings

By looking at the Lovelock theorem one can infer that the gravity model given by [1] cannot be applicable for all types of 4D Einstein–Gauss–Bonnet (EGB) curved space-time. The reason for this is that in 4D space-time, the Gauss–Bonnet invariant is a total derivative and hence it does not contribute to gravitational dynamics. Hence, the authors of [2] presented an alternative consistent EGB gravity model instead of [1] by applying a break-of-diffeomorphism property. In this work, we use the alternative model to produce a de Sitter (dS)/Anti-de Sitter (AdS) black hole metric and then investigate its thermodynamic behavior in the presence of a cloud of Nambu–Goto strings. Mathematical derivations show that the resulting diagrams of pressure vs specific volume at a constant temperature are similar to that for a van der Waals gas/fluid in an ordinary thermodynamic system in the dS sector but not in the AdS background. From this, we infer that the black hole participates in the small-to-large black hole phase transition in the dS background, while it exhibits a Hawking–Page phase transition in the AdS background. In the latter case, an evaporating black hole eventually reaches an AdS vacuum space because of its instability.     

Bidirectional Quantum Secure Direct Communication Network Protocol with Hyperentanglement

We propose a bidirectional quantum secure direct communication (QSDC) network protocol with the hyperentanglment in both the spatial-mode ad the polarization degrees of freedom of photon pairs which can in principle be produced with a beta barium borate crystal. The secret message can be encoded on the photon pairs with unitary operations in these two degrees of freedom independently. Compared with other QSDC network protocols, our QSDC network protocol has a higher capacity as each photon pair can carry 4 bits of information. Also, we discuss the security of our QSDC network protocol and its feasibility with current techniques.