a 65-year - old male patient who presented with right arm weakness revealed an acute focal cortical ischemic change in the borderzone type on the diffusion - weighted image ( not shown ) .
driver ( medtronic ireland , minneapolis , mn ) with 2.25 mm diameter and 8 mm length were deployed to 2.35 mm diameter at 12 atmosphere of balloon pressure .
our institutional review board approved this study , and we obtained written informed consent from the patient .
three - dimensional ( 3d ) rotational angiography was obtained in axiom artis zee ( siemens , erlangen , germany ) and then the images were transferred to syngo workstation ( version vb 15d ) to generate 3d angiography with reducing fov and adjusting the target vessel .
the final image was transformed to stl format in 256  256 pixels . after smoothing of vessel surface and trimming of branch vessels in magics rp ( materialise , belgium ) ,
, we carried out cfd and fully coupled fsi simulations with an assumption of rigid wall and compliant wall , respectively .
blood flow was assumed to be newtonian , homogeneous flow with density 1050 kg / m , viscosity of 3.5  10 pas .
the si physical unit of dynamic viscosity is the pascalsecond ( pas ) , ( equivalent to ns / m , or kg/(ms ) ) . to achieve true boundary conditions corresponding with patient - specific , flow velocity before and after installed stent
were obtained by transcranial doppler ultrasounds method from gated phase contrast mr angiography in an age - matched male .
pulsatile flow rates waveform based on in vivo measurement data was specified at the posterior communicating artery , the anterior cerebral artery and the middle cerebral artery for pre- and post - stent simulations respectively and time - varying pressure condition according to perktold et al .  
was imposed in the internal carotid artery inlet . in simulations with compliant wall model ,
intracranial arterial wall and plaque components were assumed to be isotropic , linear , elastic with a constant young modulus e = 2.6  10 pa , density of p = 1000 kg / m and a poisson 's ratio v = 0.49 .
to simplify the complex stent geometry , stent was modeled as shell type tube with mechanical property of density e = 2.2  10pa , p = 8420.39 kg / m and possion 's ratio v = 0.234 ( co - cr alloy , astm f563 ) .
. computational analysis of blood flow in intracranial blood vessels was performed using the commercial finite element software adina version 8.7.2 ( adina r & d , inc . ,
fluid domain was divided into 525,659 tetrahedral elements and elastic wall discretization yielded 77,264 shell elements for pre - stent case . because of rigorously complicated flow field through stenosis , dense mesh with gradual ratio change from normal to lesion region was required and especially computational mesh had the highest density around stenosis region . in case of post - stent model , 354,009 tetrahedral fluid elements and 55,359 shell elements were used for fluid and structure domain respectively . in order to prevent torsion and rotation of vessel , constraints at intracranial ends were placed and rotation motions around global coordinate axis system were fixed .
we compared difference of wss before and after stenting regarding to maximum wss at systole and enddiastole .
a 65-year - old male patient who presented with right arm weakness revealed an acute focal cortical ischemic change in the borderzone type on the diffusion - weighted image ( not shown ) .
driver ( medtronic ireland , minneapolis , mn ) with 2.25 mm diameter and 8 mm length were deployed to 2.35 mm diameter at 12 atmosphere of balloon pressure .
our institutional review board approved this study , and we obtained written informed consent from the patient .
three - dimensional ( 3d ) rotational angiography was obtained in axiom artis zee ( siemens , erlangen , germany ) and then the images were transferred to syngo workstation ( version vb 15d ) to generate 3d angiography with reducing fov and adjusting the target vessel .
the final image was transformed to stl format in 256  256 pixels . after smoothing of vessel surface and trimming of branch vessels in magics rp ( materialise , belgium ) ,
in the current study , we carried out cfd and fully coupled fsi simulations with an assumption of rigid wall and compliant wall , respectively .
blood flow was assumed to be newtonian , homogeneous flow with density 1050 kg / m , viscosity of 3.5  10 pas .
the si physical unit of dynamic viscosity is the pascalsecond ( pas ) , ( equivalent to ns / m , or kg/(ms ) ) . to achieve true boundary conditions corresponding with patient - specific , flow velocity before and after installed stent
were obtained by transcranial doppler ultrasounds method from gated phase contrast mr angiography in an age - matched male .
pulsatile flow rates waveform based on in vivo measurement data was specified at the posterior communicating artery , the anterior cerebral artery and the middle cerebral artery for pre- and post - stent simulations respectively and time - varying pressure condition according to perktold et al .  
was imposed in the internal carotid artery inlet . in simulations with compliant wall model ,
intracranial arterial wall and plaque components were assumed to be isotropic , linear , elastic with a constant young modulus e = 2.6  10 pa , density of p = 1000 kg / m and a poisson 's ratio v = 0.49 . to simplify the complex stent geometry , stent
was modeled as shell type tube with mechanical property of density e = 2.2  10pa , p = 8420.39 kg / m and possion 's ratio v = 0.234 ( co - cr alloy , astm f563 ) .
computational analysis of blood flow in intracranial blood vessels was performed using the commercial finite element software adina version 8.7.2 ( adina r & d , inc . ,
fluid domain was divided into 525,659 tetrahedral elements and elastic wall discretization yielded 77,264 shell elements for pre - stent case . because of rigorously complicated flow field through stenosis , dense mesh with gradual ratio change from normal to lesion region was required and especially computational mesh had the highest density around stenosis region . in case of post - stent model , 354,009 tetrahedral fluid elements and 55,359 shell elements were used for fluid and structure domain respectively . in order to prevent torsion and rotation of vessel , constraints at intracranial ends were placed and rotation motions around global coordinate axis system were fixed .
we compared difference of wss before and after stenting regarding to maximum wss at systole and enddiastole .
the wall shear stress distributions from cfd simulation with rigid wall assumption as well as fsi simulation before and after stenting were compared in fig .
overall , the difference of wss between rigid wall and compliant wall model both in pre- and post - stent case is only minor except at the stenosis region .
in the case of pre - stenting , fsi simulation created noticeable reduction of maximum wss at the stenosis compared to the simulation with rigid wall assumption . before stent implantation , maximum wsss at the stenosis were calculates as 167 pa at systole and 109 pa at end - diastole in cfd simulation with rigid wall assumption and as 138 pa at systole and 104 pa at end - diastole in fsi simulation with considering wall compliance ( fig .
these wss values were greatly reduced after stenting to 15~20 pa at systole and 3~5 pa at end - diastole in cfd simulation , which are similar in fsi simulations .
this revealed that the structural stress at the stenosis region is not rather high compared with the one on normal vessel , but in similar range even before implanting stent .
this implies that wss may play more important role in plaque rupture by fissuring of thin fibrous cap than structural stress on plaque .
fluid - structure interaction ( fsi ) analysis has been employed as a powerful tool because it combines blood flow simulation through cfd with wall tensile stress ( wts ) analysis in the plaque region by the finite element method .
blood flow through a compliant artery requires appropriate fluid - structure coupling in order to account for the interaction between the flowing blood and the deforming arterial wall . with anatomically realistic plaque geometry ,
fsi simulation is able to provide detailed stress analysis for the plaque which has been demonstrated with individual case studies .
as we could demonstrated the difference of fsi compared to cfd without considering the fsi , ultimately it is hoped that biomechanical simulations based on patient - specific data will serve as a reliable and robust clinical tool in the characterization of lesions and choice of treatment plan .
although various arterial sites have been analyzed by fsi such as the carotid arterial bifurcation , an abdominal aortic aneurysm model due to relative ease of imaging and simpler physical boundary conditions , and a stented aortic aneurysm model , our results revealed that the application of cfd and fsi can also be possible in the intracranial vessels which has small diameter in the range of 2 - 3 mm [ 13 , 19 , 20 ] .
risk assessment based solely on the degree of luminal stenosis will tend to underestimate the severity of the atherosclerotic lesion .
for these reasons , there is a need for more sensitive methods to risk stratify patients with intracranial atherosclerotic disease .
currently , patient - specific 3d fsi simulations suffer from two major bottle necks in terms of the time required to get from medical images to mechanical stress and strain predictions .
the first challenge is the generation of a suitable computational mesh that can accurately represent the components of the diseased artery and can provide meaningful numerical results .
the second challenge relates to the time and computational resources required to solve realistic finite element models on a patient - specific basis .
the two challenges are inter - related , and reduction of one generally means an increase in the other .
to estimate stress levels in the fibrous cap , fsi analysis has emerged as a tool combining blood - flow simulation through cfd with finite element analysis of the corresponding stress levels in the surrounding tissues .
this study has some limitations . intrinsically , vessel wall is nonlinear hyperelastic material and even its properties vary with locations and the components of vessel , in particular , at the stenosis with plaque .
we do note , however , that our finding that wss would be essential in plaque rupture compared to structural stress in the perspective of hemodynamics would hold because the narrower and thicker vessel in stenosis is , the less structural stress , but higher wss would be produced .
however , this would be proper assumption in modeling because real stent would be relatively rigid after implantation with outward expanding narrowed vessel .
the assumption of newtonian rheology made in this study have also shown to be of relatively minor influence on hemodynamic characteristics in large vessel .
 purposeimage - based computational models with fluid - structure interaction ( fsi ) can be used to perform plaque mechanical analysis in intracranial artery stenosis . 
 we described a process in fsi study applied to symptomatic severe intracranial ( m1 ) stenosis before and after stenting.materials and methodsreconstructed 3d angiography in stl format was transferred to magics for smoothing of vessel surface and trimming of branch vessels and to hypermesh for generating tetra volume mesh from triangular surface - meshed 3d angiogram . computational analysis of blood flow in the blood vessels was performed using the commercial finite element software adina ver 8.5 . 
 the distribution of wall shear stress ( wss ) , peak velocity and pressure was analyzed before and after intracranial stenting.resultsthe wall shear stress distributions from computational fluid dynamics ( cfd ) simulation with rigid wall assumption as well as fsi simulation before and after stenting could be compared . 
 the difference of wss between rigid wall and compliant wall model both in pre- and post - stent case is only minor except at the stenosis region . 
 these wss values were greatly reduced after stenting to 15~20 pa at systole and 3~5 pa at end - diastole in cfd simulation , which are similar in fsi simulations.conclusionour study revealed that fsi simulation before and after intracranial stenting was feasible despite of limited vessel wall dimension and could reveal change of wss as well as flow velocity and wall pressure . 
MATERIALS AND METHODS
Subject-specific Geometry Reconstruction
Computational Fluid Dynamics
RESULTS
DISCUSSION