- Open Access
Evaluation of needle movement effect on root canal irrigation using a computational fluid dynamics model
BioMedical Engineering OnLinevolume 18, Article number: 52 (2019)
Irrigation is considered to be a critical part of root canal treatment. However, little is known about the effect of needle movement on the irrigation process. Therefore, this study aimed to investigate the influence of the syringe and needle movement on root canal irrigation using a three-dimensional computational fluid dynamics (CFD) numerical model.
The CFD codes Flow-3D was adopted to simulate the root canal irrigation process with the syringe and needle moving up and down in motions at different amplitudes and frequencies. One stationary needle was adopted to allow comparison with the needles in up-and-down motions. Six cases where the needles were moving up and down with different amplitudes and frequencies were used to investigate the relationships between the motion of needle and irrigation efficacy.
The stationary needle gained relatively higher flow velocity and apical pressure all through the irrigation process, while the needles in constant up-and-down motions exhibited lower mean flow velocity and apical pressure. The larger the amplitude, the less mean flow velocity and apical pressure were developed. In addition, the needles moving with different frequencies were similar in the terms of irrigant replacement and apical pressure.
To avoid periapical extrusion accidents while obtaining adequate irrigant replacement, the needle should be moving up and down with a moderate amplitude during manual root canal irrigation; and the motion frequency was not highly relevant in terms of the irrigation efficiency.
Irrigation is considered to be a critical part of root canal treatment because it allows for the debridement and disinfection of the root canal, which cannot be achieved by root canal instrumentation alone [1, 2]. In addition, several studies have demonstrated that irrigation of the apical third of the root canal is extremely challenging not only because of the complexity of the root canal structure [3, 4], but also due to the fact that the microorganisms in this area principally communicate with the host and may cause periapical inflammation .
Among all the irrigation techniques, conventional syringe irrigation is still considered as the most widely used method [6, 7], with the irrigant being delivered into the root canal by a syringe and needle, either passively or with agitation. The latter is achieved by moving the needle up and down the canal space. Moreover, using this conventional irrigation method, the depth of needle penetration and the volume of irrigant could be comparatively better controlled [8, 9]. Many factors have been shown to be closely related to efficient syringe irrigation. For instance, high flow rates , closer approach of the needle to the apex , and larger irrigation volume . However, since the irrigant is being injected directly under positive pressure into the canal, these measures might potentially pose a threat to a well-known complication: periapical extrusion of irrigant , which may result in severe tissue damage and postoperative pain . Moreover, it would be difficult to standardize a limit for a safe irrigation protocol since irrigant extrusion in vivo is a multifactorial process and it is challenging to build a clinically relevant model to validate various factors.
Recently, by means of mathematical modeling and computer simulation, computational fluid dynamics (CFD) has been adopted to simulate the flow pattern developed within the root canal [13,14,15,16,17]. One outstanding advantage of CFD modeling is the visualization of the flow pattern developed within the root canal that experimental measurements are difficult to perform. Moreover, in contrast to experiments, various flow characteristics such as flow velocity and pressure could be easily obtained by CFD modeling, especially in microscale flow situations like the flow developed in the apical part of root canal . In addition, the reliability of CFD root canal model has also been validated by various studies by means of high-speed imaging and particle image velocimetry analysis [18, 19]. Therefore, CFD might be considered as a promising technique to investigate the efficiency of syringe irrigation and the potential risk of irrigant extrusion. However, in almost all CFD analyses, the syringes and needles remain stationary and it has been pointed out that the needle should be kept in motion during manual irrigation to facilitate decontamination and reduce the possibility of irrigant extrusion [20, 21].
Presently, little is known about the effect of needle movement on the irrigation process and there is limited evidence for the recommendation of a safe yet effective irrigation protocol during root canal treatment. As a result, the present study established a three-dimensional CFD model to simulate the root canal irrigation process with the syringe and needle moving up and down in motions at different amplitudes and frequencies. The aim of this study was to evaluate the relationships between irrigation efficacy and the amplitudes and frequencies of the up-and-down motions. Moreover, the potential for periapical extrusion under different amplitudes and frequencies was investigated and the optimal irrigation strategy was also explored.
The geometry construction for a CFD irrigation model
The geometry of the root canal and the needle was created according to a previous study  with Auto CAD 2012 software (AutoDesk, San Rafael, CA, USA) (Fig. 1). A 30-G flat needle was selected in the present study and the external diameter, internal diameter, and length of the needle were 320 µm, 196 µm, and 31 mm, respectively. All the needles were initially placed at 3.5 mm above the apical foramen to facilitate standardization among different cases. In addition, point A (1.5 mm from the apical foramen) was selected [14, 15] to analyze the variations of flow velocity under different cases (Fig. 1). In this study, 1% solution of sodium hypochlorite (NaOCl), with a density of 1.04 g/cm3, and a viscosity of 0.986 × 10−3 Pa s , was adopted and assumed to be an incompressible and homogeneous Newtonian fluid . The inlet flow velocity of the needle was set to be 8.6 m/s, which was consistent with an irrigant flow rate of 0.26 ml/s through a 30-G flat needle .
The governing equations
According to Sinaiski , the general condition for the beginning of turbulence has been established by Reynolds number (Re), which could be calculated by the formula Re = ρvD/µ (ρ was the fluid density, v was the fluid velocity, D was the domain diameter, and µ denoted the fluid viscosity). In the present study, ρ = 1.04 g/cm3, µ = 0.986 × 10−3 Pa s, v = 4.3 m/s (the mean velocity, which was half of the inlet flow velocity of the needle), D ranged from 0.63 mm to 0.81 mm for the two extreme positions of the needle tip (at 3.5 mm and 6.5 mm from the apex). Therefore, the Reynolds number in the studied case ranged from 2857 to 3674 depending on the position of the needle tip. In addition, the critical value of Re (Recr) for a pipe flow is approximately equal to 2800. Namely, the flow remains laminar until the Re stays below Recr, whereas at Re > Recr, the flow becomes turbulent . Therefore, the flow motions in the present study were modeled by solving the Reynolds-Averaged Navier–Stokes (RANS) equation and the k–ε turbulence model . Due to the adoption of the FAVOR mesh processing technique in the Flow-3D codes, the area fraction and volume fraction were introduced in the governing equations. Therefore, the governing equations for flow simulations in the present study were as follows:
where x, y, z were the coordinates; ρ was the fluid density; t represented the time; ux, uy, uz denoted the fluid velocity in the x-, y-, z- direction, respectively; Ax, Ay, Az signified the area fraction in the x-, y-, z- direction; VF meant the volume fraction; Gx, Gy, Gz were the gravitational acceleration in the x-, y-, z- direction, respectively; and fx, fy, fz denoted the viscous force due to the acceleration in the x-, y-, z- direction, respectively.
The standard k–ε turbulence model was adopted in this study to achieve the closure of the governing equations. Therefore, the k–ε equations after the introduction of the area fraction and the volume fraction could be expressed as
where kT signified the turbulence energy; εT was the dissipation rate of the turbulence energy; PT meant the parameter related to the velocity gradient; GT denoted the parameter related to buoyancy; and C1 = 1.44, C2 = 1.92, C3 = 0.2.
The VOF method was applied to track the free water surface (the interface between water and air). After considering the area fraction and volume fraction, the VOF equation was as follows :
where α denoted the volume fraction ( = 1 meant the cell was full of water).
CFD model design and mesh generation
The three-dimensional CAD canal and needle geometrical model was created and imported into the Flow-3D codes (Flow Science Inc., USA). The FAVOR technique was applied in the process of mesh generation and the computational grids was 0.02 mm, with a total of 1,647,300 hexahedral elements (Fig. 2). The fluid was ejected into the canal from the outlet of the needle and flew out through the orifice of the root canal where atmospheric pressure was imposed. The canal and needle were considered to be fully filled with irrigant. Gravity factor was directed at the direction of the negative y-axis.
The mass source was used to simulate the entrance of the needle lumen into the root canal and by changing the displacement of the mass source, the movement of the needle in clinical practice could be modeled. One needle was kept stationary at 3.5 mm from the apical foramen during the whole irrigation process (Case0). Three needles were kept in motion with amplitudes of 1, 2, 3 mm at a frequency of 1 Hz (CaseA1, CaseA2, CaseA3), which meant the needles were in constant 2, 4, 6 mm up-and-down motions every second, respectively. Another three needles were set to be moving up and down in 4 mm motions with different frequencies. To be specific, three frequencies (F = 1, 2, 3) were adopted, which meant the needle made 1,2,3 up-and-down motions in one second, respectively (CaseF1, CaseF2, CaseF3). Moreover, the movement locus of the needle was achieved by firstly moving away from the apex and then approaching to the starting point (3.5 mm from the apical foramen) in constant motions.
Boundary conditions and numeric simulation
No-slip boundary condition was adopted to the walls of the root canal and of the needle with the assumptions that the walls were rigid, smooth, and impermeable. The walls of the root canal were defined as wall boundary condition, while other boundaries were defined as pressure boundary condition, which meant the relative pressure was zero. (Moved to “The Geometry construction for a CFD irrigation model” section.)
When solving the governing equations, in order to get stable numerical solution, the time step (∆t) must be less than or at most equal to ∆x/c based on the definition of Courant number (∆x was the cell size in the direction of the velocity, c was the magnitude of the velocity through the cell) . In this study, ∆x was 2 × 10−5 m and c was 8.6 m/s. Therefore, the minimum time step for this study was 2.3 × 10−6 s, however, in order to improve the computational accuracy, a smaller time step of 1 × 10−7 s was used throughout the whole calculations. In addition, the computational time for each case was set to 10 s. When the flow developed fully and the parameters of interest reached a steady state, 1 s was adopted in order to allow one period (T) for the cases kept in motion. Computations were carried out on a Windows 7 Computer workstation with 2 CPU and 128 G RAM. A series of computations were performed where the needle was kept stationary or moved up and down in motions. The flow fields computed were compared including flow pattern, velocity magnitude, and apical pressure.
Statistical analysis was analyzed with statistical software (IBM SPSS Statistics, v22.0; IBM Corp). 1-way analysis of variance (ANOVA) was performed to identify statistically significant difference among different groups ( = 0.05).
First of all, the 3D CFD model was validated by comparing with previous numerical results [14, 15]. The mean apical pressure of the stationary case (18.04 kPa) (Table 1) generated in the present study agreed well with the mean apical pressure of previous CFD analysis using the same canal model and irrigation speed [14, 15]. Differences with other CFD studies might be attributed to different experimental settings and turbulence models . Therefore, the validity of the current model could be verified.
Comparison of the stationary needle and the needles in motions
The distribution of velocity field, velocity vector, vorticity, and pressure field at 1/2 T (T = one period) when the needle remained static (Case0) is presented in Fig. 3. In addition, the distribution of the velocity field (Fig. 4), velocity vector (Fig. 5), vorticity (Fig. 6), and pressure field (Fig. 7) of CaseA1 at certain time: 1/4 T, 1/2 T, 3/4 T, and T in the apical part of the root canal is also shown and they provide a detailed overview of the variations of fluid field of CaseA1. As shown, the value of flow parameters of CaseA1 gradually decreased as the needle moved away from the apical foramen, while increased when the needle approached the apex. In addition, the variations of velocity and pressure field of Case0 and CaseA1 within a second were also demonstrated by one video (Additional file 1: Video S1).
The variations of y-velocity at point A and apical pressure of the stationary needle (Case0) and needle in 2 mm up-and-down motions (CaseA1) within a second are presented in Figs. 8a, b and 9a, b. Like the fluid field developed within the canal (Figs. 3, 4, 5, 6, 7), the distinction between the stationary case and CaseA1 was still eminent. To be specific, Case0 exhibited relatively higher irrigant velocity and apical pressure all the time, while CaseA1 gained smaller irrigant velocity and apical pressure, especially when the needle moved away from the apex. In addition, Table 1 shows that the mean y-velocity magnitude at point A, mean apical pressure, and the percentage of time when the y-velocity magnitude at point A is higher than 0.1 m/s (P0.1) of Case0 were all larger than that of CaseA1.
Comparison of the needles in motions with different amplitudes
To investigate the influence of the amplitude on irrigation efficiency, various amplitudes of 1, 2, 3 mm (CaseA1, CaseA2, and CaseA3) were adopted. Figures 8, 9 show the variations of y-velocity at point A and apical pressure of needles with different amplitudes. Furthermore, the Vm, Pm, and P0.1 of CaseA1, CaseA2, and CaseA3 are summarized in Table 1. As shown, the flow velocity and apical pressure in different cases presented a similar trend, that is, the lager the amplitude of the needle, the more unstable and smaller the flow velocity and apical pressure. However, when comparing P0.1 of different amplitudes, CaseA1 is similar to CaseA2 (86% and 87%), while CaseA3 gained smaller P0.1 (69%).
Comparison of the needles in motions with different frequencies
The effect of frequency on irrigation efficiency was also studied and the variations of y-velocity at point A and apical pressure of needles with different frequencies (F = 1, 2, 3, CaseF1, CaseF2, and CaseF3) are presented in Fig. 10. The magnitude of flow velocity and apical pressure in different cases underwent periodic alterations and were inconsistent with the needle movement locus. Moreover, the Vm, Pm, and P0.1 are also shown in Table 1 and the results show that there were no significant differences among CaseF1, CaseF2, and CaseF3 in the magnitude of mean flow velocity (1.48, 1.58, 1.80), P0.1 (87%, 88%, 88%), and apical pressure (13.41, 14.33, 13.75) (P > 0.05).
CFD simulations have previously shown clinical relevance in a range of studies [13,14,15,16,17], but without incorporating the motion of the needles. This study described the first simulation to investigate the influence of needle movement on the root canal irrigation using the CFD model. The dynamic alteration of the flow characteristics such as flow velocity and pressure were examined when the needle was kept stationary or moved up and down in constant motions with different amplitudes and frequencies.
The stationary case (Case0) presented relative high flow velocity all the time, while the cases in motions underwent evident periodic alterations of velocity with much lower mean flow velocity. Since it has been pointed out by previous studies that the distribution of y-velocity along the longitudinal axis of the root canal (y-axis) acts as a vital role in determining the irrigant replacement and the flush effect [13,14,15,16,17] and velocities higher than 0.1 m/s were considered clinically significant for adequate irrigant replacement [14, 15], it could be speculated that the stationary case would achieve improved irrigation replacement, as depicted by P0.1 (the percentage of time when the y-velocity magnitude at point A is higher than 0.1 m/s) (Table 1). However, besides irrigation replacement, other factors should also be taken into consideration when assessing an irrigation protocol. For example, a chief complication during irrigation is the extrusion of irrigant into periapical tissues, leading to severe pain and tissue destruction [3, 12]. If the pressure at the apical foramen exceeds the backpressure of the periradicular tissue, periapical extrusion would probably happen and complications caused by cytotoxic effects of irrigant might occur . As stated before, the stationary case gained much higher pressure all through the irrigation process, hence it is more likely to force the irrigant and debris into the periapical tissue. In comparison, the needles in motions presented a distinct decrease in pressure when the needle moved away from the apex and thus gained much lower apical pressure. Such trend would pose much lower risks of injury to the periradicular tissue. Moreover, cases in motions with smaller amplitude like CaseA1 and CaseA2 could also gain adequate irrigant replacement (86% and 87%, respectively). Therefore, it is reasonable to assume that the needles in motions would produce safer irrigation effect yet adequate irrigant replacement if moving the needle up and down with smaller amplitude.
When comparing cases in motions with different amplitudes, the overall variation trend of the flow velocity and apical pressure was similar and in parallel with the periodic motion of the needles. As shown, the irrigant replacement was not sacrificed when increasing the amplitude from 1 mm (CaseA1, 86%) to 2 mm (CaseA2, 87%). However, when the amplitude was increased to 3 mm in CaseA3, the irrigant replacement became less efficient (69%) (Table 1). Therefore, it could be estimated that the irrigant replacement would be insufficient if the moving amplitude is larger than 3 mm. On the other hand, with the decrease of the moving amplitude, the apical pressure, as well as the risk of periapical extrusion would be higher [14,15,16,17]. Therefore, increasing the moving amplitude would be safer during irrigation. However, because there is no definite evidence on the minimum irrigant pressure that leads to extrusion, the risk of apical extrusion can only be estimated by comparison between different cases. Therefore, a reasonable compromise would be to move the needle with a moderate amplitude to get a relatively satisfactory flush effect and lower apical pressure meanwhile during irrigation.
Furthermore, a noteworthy finding of this study was the influence of the moving frequencies of the needle on irrigation efficiency. As indicated by the results, no significant differences were observed among different cases in the magnitude of mean flow velocity, apical pressure, and the percentage of clinically efficient velocities (P0.1). Though the mean flow velocity is a little higher when increasing the frequencies, the P0.1 and mean apical pressure in the three cases were almost the same, indicating similar irrigant replacement and risk of apical extrusion. Hence, the moving frequencies of the needles were not highly related to the irrigation efficiency and thus might not be the main concern during irrigation process. Therefore, it would be unnecessary to move the needles very fast or very slow to improve irrigation efficiency and should be just the practitioners’ personal preference in terms of motion frequency.
Nevertheless, the findings of this study should further be evaluated before any clinical suggestion is made. First of all, it is hard to say exactly which amplitude should be employed since there are no safety limits set for the apical wall pressure so far and many other additional factors such as the needle depth, needle type, and the apical foramen size and taper of canal would also have great impacts on the irrigation outcome. In addition, there are some limitations of the model applied. For example, the geometry of the root canal was simplified in order to allow comparison with previous literature [14, 15]. The walls of the root canals were assumed to be smooth and impermeable and this assumption may not be clinically relevant. Moreover, the gravity would vary according to the canal and how the patient sits on the chair was not taken into consideration in this study. Therefore, further studies should be conducted to provide a more comprehensive guidance to clinical practice.
The present study presented an initial attempt to use a CFD model to investigate the influence of the needle movement on the root canal irrigation. Within the limitations of the study, this numerical model provides some basic understanding of the effect of needle movement on the flow velocity and apical pressure that experimental measurements are difficult to perform. First of all, the stationary case gained higher irrigant replacement yet would also pose much higher risk of injury to the periradicular tissue. Therefore, the fear of irrigant extrusion has led to the suggestion of keeping the needle in motion while manually irrigating the canal. In addition, comparison of different cases in motions demonstrated that a needle moving up and down with a moderate amplitude during manual root canal irrigation would result in better flush effect while reducing the risk of periapical extrusion; and the motion frequency was not highly relevant in terms of the irrigation efficiency and should be practitioners’ personal preferences.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Haapasalo M, Shen Y, Wang Z, Gao Y. Irrigation in endodontics. Br Dent J. 2014;216:299–303.
Paqué F, Balmer M, Attin T, Peters OA. Preparation of oval-shaped root canals in mandibular molars using nickel-titanium rotary instruments: a micro-computed tomography study. J Endod. 2010;36:703–7.
Desai P, Himel V. Comparative safety of various intracanal irrigation systems. J Endod. 2009;35:545–9.
Howard RK, Kirkpatrick TC, Rutledge RE, Yaccino JM. Comparison of debris removal with three different irrigation techniques. J Endod. 2011;37:1301–5.
Ricucci D, Siqueira JF Jr. Biofilms and apical periodontitis: study of prevalence and association with clinical and histopathologic findings. J Endod. 2010;36:1277–88.
Dutner J, Mines P, Anderson A. Irrigation trends among American Association of Endodontists members: a web-based survey. J Endod. 2012;38:37–40.
Willershausen I, Wolf TG, Schmidtmann I, Berger C, Ehlers V, Willershausen B, Briseno B. Survey of root canal irrigating solutions used in dental practices within Germany. Int Endod J. 2015;48:654–60.
van der Sluis LW, Gambarini G, Wu MK, Wesselink PR. The influence of volume, type of irrigant and flushing method on removing artificially placed dentine debris from the apical root canal during passive ultrasonic irrigation. Int Endod J. 2006;39:472–6.
Haapasalo M, Shen Y, Qian W, Gao Y. Irrigation in endodontics. Den Clin N Am. 2010;54:291–312.
Sedgley CM, Nagel AC, Hall D, Applegate B. Influence of irrigant needle depth in removing bioluminescent bacteria inoculated into instrumented root canals using real-time imaging in vitro. Int Endod J. 2005;38:97–104.
Sedgley C, Applegate B, Nagel A, Hall D. Real-time imaging and quantification of bioluminescent bacteria in root canals in vitro. J Endod. 2004;30:893–8.
Hülsmann M, Rödig T, Nordmeyer S. Complications during root canal irrigation. Endod Topics. 2009;16:27–63.
Boutsioukis C, Lambrianidis T, Kastrinakis E. Irrigant flow within a prepared root canal using various flow rates: a Computational Fluid Dynamics study. Int Endod J. 2009;42:144–55.
Boutsioukis C, Lambrianidis T, Verhaagen B, Versluis M, Kastrinakis E, Wesselink PR, van der Sluis LW. The effect of needle-insertion depth on the irrigant flow in the root canal: evaluation using an unsteady computational fluid dynamics model. J Endod. 2010;36:1664–8.
Boutsioukis C, Verhaagen B, Versluis M, Kastrinakis E, van der Sluis LW. Evaluation of irrigant flow in the root canal using different needle types by an unsteady computational fluid dynamics model. J Endod. 2010;36:875–9.
Boutsioukis C, Gogos C, Verhaagen B, Kastrinakis E, Wesselink PR, Sluis LW. The effect of root canal taper on the irrigant flow: evaluation using an unsteady Computational Fluid Dynamics model. Int Endod J. 2010;43:909–16.
Shen Y, Gao Y, Qian W, Ruse ND, Zhou X, Wu H, Haapasalo M. Three-dimensional numeric simulation of root canal irrigant flow with different irrigation needles. J Endod. 2010;36:884–9.
Boutsioukis C, Verhaagen B, Versluis M, Kastrinakis E, van der Sluis LW. Irrigant flow in the root canal: experimental validation of an unsteady Computational Fluid Dynamics model using high-speed imaging. Int Endod J. 2010;43:393–403.
Gao Y, Haapasalo M, Shen Y, Kastrinakis E, Sluis LW. Development and validation of a three dimensional computational fluid dynamics model of root canalirrigation. J Endod. 2009;35:1282–7.
Druttman ACS, Stock CJR. An in vitro comparison of ultrasonic and conventional methods of irrigant replacement. Int Endod J. 1989;22:174–8.
Metzger Z, Solomonov M, Kfir A. The role of mechanical instrumentation in the cleaning of root canals. Endod Topics. 2013;29:87–109.
Guerisoli DMZ, Silva RS, Pecora JD. Evaluation of some physico-chemical properties of different concentrations of sodium hypochlorite solutions. Braz Endod J. 1998;3:21–3.
Tilton JN. Fluid and particle dynamics. In: Perry RH, Green DW, Maloney JO, editors. Perry’s chemical engineer’s handbook. 7th ed. New York: McGraw-Hill; 1999.
Boutsioukis C, Lambrianidis T, Kastrinakis E, Bekiaroglou P. Measurement of pressure and flow rates during irrigation of a root canal ex vivo with three endodontic needles. Int Endod J. 2007;40:504–13.
Sinaiski EG. Hydromechanics: theory and fundamentals. 1st ed. Weinheim: Springer; 2011. p. 124–5.
Rodi W. Turbulence models and their application in hydraulics-state-of-the-art review. 3rd ed. Rotterdam: A.A. Balkema; 1993.
Hirt CW, Nichols BD. Volume of fluid (VOF) method for the dynamics of free boundaries. J Comput Phys. 1982;39:201–25.
Anderson JD. Computational fluid dynamics: the basics with applications. 1st ed. New York: McGraw Hill; 1995. p. 153–65.
Tinaz AC, Alacam T, Uzun O, Maden M, Kayaoglu G. The effect of disruption of apical constriction on periapical extrusion. J Endod. 2005;31:533–5.
This work was supported by National Natural Science Foundation of China (81771122) and Department of Science and Technology of Sichuan Province (2018SZ0037, 19YYJC2625).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.