基于心电脉搏信号的实时动态血管建模方法

当前血管微创介入手术三维图像导航系统中血管模型是由病人的CT或MRI二维扫描图像经三维重建生成的,为静态模型,不能准确反映术中血管的形态。针对以上问题,研究了基于实时心电脉搏信号的动态血管模型建模方法,主要包括编程建立参数化动态血管模型、实时采集处理心电脉搏信号获得模型运动参数、基于参数的动态模型建立三个部分。通过实验验证了动态模型的实时性和准确性。研究可以为血管微创介入手术图像导航系统提供较准确的实时动态血管模型,提高手术导航系统的视觉临场感,为医生提供重要的视觉参考。     

Pulsatile flow of viscous and viscoelastic fluids in constricted tubes

The unsteady flow of blood through stenosed artery, driven by an oscillatory pressure gradient, is studied. An appropriate shape of the time-dependent stenoses which are overlapped in the realm of the formation of arterial narrowing is constructed mathematically. A msathematical model is developed by treating blood as a non-Newtonian fluid characterized by the Oldroyd-B and Cross models. A numerical scheme has been used to solve the unsteady nonlinear Navier-stokes equations in cylindrical coordinate system governing flow, assuming axial symmetry under laminar flow condition so that the problem effectively becomes two-dimensional. Finite difference technique was used to investigate the effects of parameters such as pulsatility, non-Newtonian properties and the flow time on the velocity components, the rate of flow, and the wall shear stress through their graphical representations quantitatively at the end of the paper in order to validate the applicability of the present improved mathematical model under consideration.     

Numerical investigation of non-uniform magnetic field effects on the blood velocity and magnetic nanoparticles concentration inside the vessel

Among other applications, Magnetic nanoparticles (MNPs) are used as drug carriers for guided drug delivery purposes. The aim of the current study is to numerically investigate the effect of non-uniform magnetic field on the blood flow and MNPs concentration inside a vessel. The coupled governing equations under the influence of non-uniform magnetic field are utilized and non-dimensionalized. Two dimensionless parameters, Magnetic force number (MFN) and Magnetic drift number (MDN), are introduced and their effects on the blood flow velocity and the MNPs concentration are studied. The results show for the same magnetic force number, the MNPs concentration on the vessel wall increases as the magnetic drift number increases. While under the same magnetic drift number, magnetic nanoparticles decreases as magnetic force number increases. This implies that the MNPs concentration at the entrance is low while MNPs radius is large. In addition, the greater the magnetic force and the magnetic drift numbers is the higher the deformation of streamlines becomes.     

Magnetic diagnostics for equilibrium reconstructions in the presence of nonaxisymmetric eddy current distributions in tokamaks (invited)

The lithium tokamak experiment (LTX) is a modest-sized spherical tokamak (R(0)=0.4?m and a=0.26?m) designed to investigate the low-recycling lithium wall operating regime for magnetically confined plasmas. LTX will reach this regime through a lithium-coated shell internal to the vacuum vessel, conformal to the plasma last-closed-flux surface, and heated to 300-400?°C. This structure is highly conductive and not axisymmetric. The three-dimensional nature of the shell causes the eddy currents and magnetic fields to be three-dimensional as well. In order to analyze the plasma equilibrium in the presence of three-dimensional eddy currents, an extensive array of unique magnetic diagnostics has been implemented. Sensors are designed to survive high temperatures and incidental contact with lithium and provide data on toroidal asymmetries as well as full coverage of the poloidal cross-section. The magnetic array has been utilized to determine the effects of nonaxisymmetric eddy currents and to model the start-up phase of LTX. Measurements from the magnetic array, coupled with two-dimensional field component modeling, have allowed a suitable field null and initial plasma current to be produced. For full magnetic reconstructions, a three-dimensional electromagnetic model of the vacuum vessel and shell is under development.     

Influence of different bell-shaped stenoses on the progression of the disease,atherosclerosis

In the present paper, the effect of symmetrical and asymmetrical bell-shaped stenoses on wall pressure drop, streamline contour, and rise in wall shear stress for the progression of the disease, atherosclerosis has been investigated numerically. The governing equations have been solved by finite volume method. Both steady and pulsatile flow at inlet is considered in our study. It is revealed from the study that the impact of wall pressure and peak wall shear stress on progression of disease are always high for asymmetrical shaped stenosis for both steady and pulsatile flow. The impact of asymmetrical shape on plaque deposition zone is less, if the aggravation changes the shape of stenosis due to change in stricture length only keeping percentage of restriction same. Whereas, the impact of asymmetrical shape on plaque deposition zone will be high, if shape of stenosis changes by increasing both stricture length and percentage of restriction for both steady and pulsatile flow. Impact of pulsatile nature of flow on the aggravation of disease is higher at some timesteps in comparison to steady flow.     

Thermophoretic transport and deposition of particles in vertical tube flow with variable wall temperature and thermal radiation

A particulate two phase flow with variable wall temperature has been studied for examining the deposition of particles in the thermal radiation and mixed convection flow associated with the manufacture of optical fiber preforms. The two-dimensional governing equations of continuity, momentum and energy have been solved numerically including the effects of thermal radiation and buoyancy (upward or downward flow) in the vertical tube flow. A particle trajectory model has been adopted to predict the particle transport, and P-1 approximation has been used to evaluate the radiation heat transfer. In the upward flow case, a high deposition efficiency is obtained and the deposition zone of the downward flow is broader than that of the upward flow. Thermal radiation makes the deposition zone broader and the deposition efficiency smaller.     

Numerical simulation of fully developed flow and heat transfer characteristics in a curved tube with pulsating pressure gradient

Characteristics of fluid flow and convective heat transfer of a pulsating flow in a curved tube have been investigated numerically. The tube wall is assumed to be maintained at a uniform temperature peripherally in a fully developed pulsating flow region. The temperature and flow distributions over a cross-section of a curved tube with the associated velocity field need to be studied in detail. This problem is of particular interest in the design of Stirling engine heat exchangers and in understanding the blood flow in the aorta. The time-dependent, elliptic governing equations are solved, employing finite volume technique. The periodic steady state results are obtained for various governing dimensionless parameters, such as Womersley number, pulsation amplitude ration, curvature ratio and Reynolds number. The numerical results indicate that the phase difference between the pressure gradient and averaged axial velocity increases gradually up to π/2 as Womersley number increases. However, this phase difference is almost independent of the amplitude ratio of pulsation. It is also found that the secondary flow patterns are strongly affected by the curvature ratio and Reynolds number. These, in turn, give a strong influence on the convective heat transfer from the pipe wall to the pulsating flow. The results obtained lead to a better understanding of the underlying physical process and also provide input that may be used to design the relevant system. The numerical approach is discussed in detail, and the aspects that must be included for an accurate simulation are discussed.     

Pulsatile blood flows through a bileaflet mechanical heart valve with different approach methods of numerical analysis; pulsatile flows with fixed leaflets and interacted with moving leaflets

Many researchers have investigated the blood flow characteristics through bileaflet mechanical heart valves using computational fluid dynamics (CFD) models. Their numerical approach methods can be classified into three types; steady flow analysis, pulsatile flow analysis with fixed leaflets, and pulsatile flow analysis with moving leaflets. The first and second methods have been generally employed for two-dimensional and three-dimensional calculations. The pulsatile flow analysis interacted with moving leaflets has been recently introduced and tried only in two-dimensional analysis because this approach method has difficulty in considering simultaneously two physics of blood flow and leaflet behavior interacted with blood flow. In this publication, numerical calculation for pulsatile flow with moving leaflets using a fluid-structure interaction method has been performed in a three-dimensional geometry. Also, pulsatile flow with fixed leaflets has been analyzed for comparison with the case with moving leaflets. The calculated results using the fluid-structure interaction model have shown good agreements with results visualized by previous experiments. In peak systole, calculations with the two approach methods have predicted similar flow fields. However, the model with fixed leaflets has not been able to predict the flow fields during opening and closing phases. Therefore, the model with moving leaflets is rigorously required for advanced analysis of flow fields.     

Thermal behaviour modelling of tapered optical fibres for scanning near-field microscopy

A simple model allowing the calculation of the thermal field inside a metal-coated fibre tip is presented. The approach has been based on previous temperature measurements which operated in steady state and periodic rate. The modelhas been inspired from the general theory of heat transfer inside fins, after having divided the taper into a set of layers. The advantage of the method is the possibility to consider any taper shapes. Moreover, any kind of coating thickness and external heat transfer distributions can be considered. As a mean of comparison with some previous works, results obtained for simple configurations are presented. Then, a study of the main governing parameters provides the basic thermal behaviour analysis of optical tips and a comparison with experience is given in order to confirm the validity of our approach.     

Spectral element modeling and analysis of the blood flow through a human blood vessel

As the blood flow characteristics are closely related to various cardiovascular diseases, it is very important to predict them accurate enough in an efficient way. Thus, this paper proposes a one-dimensional spectral element model for human blood vessels with varying cross-sections. The spectral element model is formulated by using the variational approach of finite element formulation. The wave solutions analytically solved in the frequency-domain to satisfy governing equations are used to determine the frequency-dependent interpolation functions. The spectral finite element model is then applied to an example blood vessel to investigate the blood flow rate and blood pressure through the blood vessel.     

彩色视网膜眼底图像血管自动检测方法

为了给视网膜图像配准、光照校正及视网膜内部病理学检测等问题提供有效依据,本文提出一种有效检测及识别彩色视网膜眼底图像血管的全自动方法。针对视网膜可见血管呈长条型管状、局部具有较好直线型结构的形态特点,本文采取适用于条状结构的组合移位滤波响应模型进行特征提取。针对血管和血管末端特征的不同,分别配置对称和非对称的两种滤波模型进行跟踪,利用组合移位滤波模型(对称和非对称)获取到的响应及G通道像素灰度值共同构建特征向量库,采用AdaBoost分类器对各个像素点进行分类判定。基于国际公共数据库DRIVE与STARE的实验结果表明,该方法针对两个标准数据库的分割结果(DRIVE:Accuracy=0.948 9,Sensitivity=0.765 7,Specificity=0.980 9;STARE:Accuracy=0.956 7,Sensitivity=0.771 7,Specificity=0.976 6)均优于已有方法,适用于彩色视网膜眼底图像的计算机辅助定量分析,可作为临床借鉴。     

A venous pulse Doppler catheter-tip flowmeter for measuring arterial blood velocity, flow, and diameter in deep arteries.

Numerous schemes have been used for measuring hemodynamic properties of deeply lying arteries; however, all have their limitations. This paper describes a new relatively nontraumatic intravenous approach that uses a catheter in connection with a pulsed ultrasonic Doppler velocity meter (PUDVM) and an echo track. The catheter was initially tested in a hydraulic model system for calibration of velocity and flow parameters. Lately, the catheter has permittted measurements of local instantaneous blood velocity, flow, and wall motion characteristics in adult Beagle dogs in the abdominal aorta and iliac artery. Evaluation studies have been conducted to compare the catheter-tip recordings with an independent method for measuring blood flow and wall motion. Coupling of this catheter-tip device with the PUDVM and echo track provides chronic measurements of hemodynamic parameters in these deep vessels which were virtually impossible to obtain previously. This technique may prove useful in monitoring vessel pathology longitudinally as well as in basic experimental situations requiring flow and arterial wall mechanical properties.     

A computational study of annular cascade flows using aq-ω turbulence model with a wall function

A computational code has been developed for steady viscous flows in three dimensional annular cascades. This code solves a special form of the thin-layer Navier-Stokes equations with a two-equationq-ω turbulence model in curvilinear coordinates using a time asymptotic method for steady state solutions. It employs a scalar implicit approximate factorization in time and a finite volume formulation with second-order upwind-differencing in space. A wall function treatment is implemented at solid boundaries for turbulence equations instead of integration to the wall to relieve gridding requirements. In order to validate the effectiveness of this code, computational studies have been made to access modeling capability for complex turbulent flow fields in three dimensional annular cascade geometries which typically include laminar-turbulent boundary layer transition. The results have been compared with both the computational studies with integration to the wall and the experimental studies. The wall function treatment was found to be reliable by predicting secondary flows and loss contours reasonably well.     

基于PMAC开放式数控系统软件开发

基于PC和多轴运动控制器的开放式数控系统是理想的开放式数控系统。介绍了基于PMAC的开放式数控系统结构形式,PMAC的差补、位置控制、伺服功能、以PMAC和PC机为硬件平台搭建了数控系统,并对其硬件构成和软件设计结构进行了分析。着重从软件设计的角度,介绍了PTALK控件的功能和作用,对数控系统软件构成进行了详细的阐述。并设计出了友好的用户界面,在实际应用中具有重要意义。     

Computational study of the mixed cooling effects on the in-vessel retention of a molten pool in a nuclear reactor

The retention of a molten pool vessel cooled by internal vessel reflooding and/or external vessel reactor cavity flooding has been considered as one of severe accident management strategies. The present numerical study investigates the effect of both internal and external vessel mixed cooling on an internally heated molten pool. The molten pool is confined in a hemispherical vessel with reference to the thermal behavior of the vessel wall. In this study, our numerical model used a scaled-down reactor vessel of a KSNP (Korea Standard Nuclear Power) reactor design of 1000 MWe (a Pressurized Water Reactor with a large and dry containment). Well-known temperature-dependent boiling heat transfer curves are applied to the internal and external vessel cooling boundaries. Radiative heat transfer has been considered in the case of dry internal vessel boundary condition. Computational results show that the external cooling vessel boundary conditions have better effectiveness than internal vessel cooling in the retention of the melt pool vessel failure.     

Determination of the optimal axial length of the electrode in an electrical capacitance tomography sensor

Electrical capacitance tomography (ECT) is one of the most successful process tomography techniques, which has been developed for more than a decade. Thus far, most of the research on ECT has been based on two-dimensional modeling, which simplifies the ECT sensor structure as two dimensional. However, certain key issues such as the three-dimensional effect of the ECT sensor are not considered at all. In this paper, the finite element method is used to solve the partial differential equation governing the ECT sensing zone, which is actually three dimensional. The sensitivity distributions between different electrode pairs are calculated and compared with those based on the two-dimensional modeling. The impacts of the length of both the measurement electrode and the guard electrode on the sensitivity distribution are investigated. The effect of the length of the measurement electrode on the normalized capacitance is also discussed. According to the numerical results presented in this paper, it is suggested that the optimal length of the measurement electrode for electrical capacitance tomography is the same as the width or diameter of the sensor.     

Contact and friction between catheter and blood vessel

Therapeutic vascular catheterization techniques are sometimes hampered by the frictional forces between the blood vessel and the catheter, when contact points of the vessel are changing and deforming during the movement of the catheter. The goal of the present study was to characterize frictional interactions between the blood vessel wall and the catheter using experimental and numerical analysis. First, the frictional force was measured with an experimental apparatus that uses a ball and flattened porcine aorta to simulate frictional forces between the catheter and the vessel. Second, catheter motion was characterized by two-dimensional numerical calculations based on the experimental results. Experimental analysis demonstrated that slip occurred and that friction coefficient between the vessel and the catheter and the deformation of the specimen were small when the contact between the ball and the aorta occurred at a small angle. The compliance of the specimen in the normal direction obtained by the experiment was by far larger than that calculated according to the Hertzian contact theory. Numerical analysis shows that this difference of the parameter of the vessel, which must be determined accurately in surgical simulator, could affect the trajectory of the catheter.     

Film condensation on horizontal tube with wall suction effects

This study performs a theoretical investigation into the problem of two-dimensional steady filmwise condensation flow on a horizontal tube with suction effects at the tube surface. An effective suction function is introduced to model the effect of the wall suction on the thickness of the liquid condensate film. The local condensate film thickness and the local Nusselt number are then derived using a simple series numerical method. The results show that the Nusselt number varies as a function of the Jakob number Ja, the Rayleigh number Ra, and the suction parameter Sw. It is found that the wall suction effect has a significant influence on the heat transfer performance. An analytical solution is derived for the mean Nusselt number for the case in which the wall suction effect is ignored. Finally, a closed-form correlation is presented for the mean Nusselt number subject to a wall suction effect.     

Calibration of electrodiffusion probes for turbulent flow measurements

Calibration of electrodiffusion (ED) probes with respect to the wall velocity gradient measurements has been experimentally studied in a fully developed turbulent channel flow over the Reynolds number range 14000–23000. In steady state conditions, direct calibration concerning the mean wall velocity gradient can be provided using simultaneous transient and steady state diffusion limiting current measurements. Indirect calibration of the ED probes regarding the turbulent fluctuation has been undertaken using a spectral analysis performed with probes of different size and geometry. This indirect calibration method has been supported using the measurements of thermo-dependence of molecular diffusivity. In this case, only one ED probe can be used for calibration. Dynamic calibration provides necessary information about critical longitudinal dimension of ED probes. When the longitudinal dimension of the probes exceeds the critical length, the probes becomes sensitive to the normal fluctuation as well. This fact has been confirmed using different orientation of the rectangular ED probe with respect to the flow direction.