![]() |
Citation: Lorenzo Azzalini, Philippe L. L’Allier, Jean-François Tanguay. Bioresorbable Scaffolds: The Revolution in Coronary Stenting?[J]. AIMS Medical Science, 2016, 3(1): 126-146. doi: 10.3934/medsci.2016.1.126
[1] | Tingting Ma, Yayun Fu, Yuehua He, Wenjie Yang . A linearly implicit energy-preserving exponential time differencing scheme for the fractional nonlinear Schrödinger equation. Networks and Heterogeneous Media, 2023, 18(3): 1105-1117. doi: 10.3934/nhm.2023048 |
[2] | Fengli Yin, Dongliang Xu, Wenjie Yang . High-order schemes for the fractional coupled nonlinear Schrödinger equation. Networks and Heterogeneous Media, 2023, 18(4): 1434-1453. doi: 10.3934/nhm.2023063 |
[3] | Michael T. Redle, Michael Herty . An asymptotic-preserving scheme for isentropic flow in pipe networks. Networks and Heterogeneous Media, 2025, 20(1): 254-285. doi: 10.3934/nhm.2025013 |
[4] | Junjie Wang, Yaping Zhang, Liangliang Zhai . Structure-preserving scheme for one dimension and two dimension fractional KGS equations. Networks and Heterogeneous Media, 2023, 18(1): 463-493. doi: 10.3934/nhm.2023019 |
[5] | Aslam Khan, Abdul Ghafoor, Emel Khan, Kamal Shah, Thabet Abdeljawad . Solving scalar reaction diffusion equations with cubic non-linearity having time-dependent coefficients by the wavelet method of lines. Networks and Heterogeneous Media, 2024, 19(2): 634-654. doi: 10.3934/nhm.2024028 |
[6] | Leqiang Zou, Yanzi Zhang . Efficient numerical schemes for variable-order mobile-immobile advection-dispersion equation. Networks and Heterogeneous Media, 2025, 20(2): 387-405. doi: 10.3934/nhm.2025018 |
[7] | Luis Almeida, Federica Bubba, Benoît Perthame, Camille Pouchol . Energy and implicit discretization of the Fokker-Planck and Keller-Segel type equations. Networks and Heterogeneous Media, 2019, 14(1): 23-41. doi: 10.3934/nhm.2019002 |
[8] | Min Li, Ju Ming, Tingting Qin, Boya Zhou . Convergence of an energy-preserving finite difference method for the nonlinear coupled space-fractional Klein-Gordon equations. Networks and Heterogeneous Media, 2023, 18(3): 957-981. doi: 10.3934/nhm.2023042 |
[9] | Raimund Bürger, Antonio García, Kenneth H. Karlsen, John D. Towers . Difference schemes, entropy solutions, and speedup impulse for an inhomogeneous kinematic traffic flow model. Networks and Heterogeneous Media, 2008, 3(1): 1-41. doi: 10.3934/nhm.2008.3.1 |
[10] | Mahmoud Saleh, Endre Kovács, Nagaraja Kallur . Adaptive step size controllers based on Runge-Kutta and linear-neighbor methods for solving the non-stationary heat conduction equation. Networks and Heterogeneous Media, 2023, 18(3): 1059-1082. doi: 10.3934/nhm.2023046 |
In this paper, we consider numerically a class of nonlinear dispersive equations [18] as follows:
$ ut−uxxt+κux+3uux=γ(2uxuxx+uuxxx), u(x,0)=u0(x), x∈R, t>0, $
|
(1.1) |
which models a variety of nonlinear dispersive phenomena depending on the parameters $ \kappa $ and $ \gamma $. If $ \kappa \ge 0 $ and $ \gamma = 1 $, it reduces to the Camassa-Holm (CH) equation:
$ ut−uxxt+κux=2uxuxx+uuxxx−3uux, u(x,0)=u0(x), x∈R, t>0, $
|
which characterizes unidirectional shallow water waves [2,3], with $ u $ representing the height of the fluid's free surface above a flat bottom (or equivalently the fluid velocity in the $ x $ direction), and $ \kappa $ denoting the critical shallow water wave speed. The CH model possesses a bi-Hamiltonian structure, is completely integrable, and has global solutions. Moreover, it has solitary waves. Compared with the Korteweg-de Vries equation, the solitary waves become peaked in the limit of $ \kappa \! \to \! 0 $ (called "peakons"), which characterize in physical context "wave-breaking" [25]. It is worth noting that the dispersive Eq (1.1), by choosing different values of $ \kappa $ and $ \gamma $ [18], may turn into other physical models, such as the Dai equation and regularized long wave equation.
The system (1.1) can be rewritten in an infinite-dimensional Hamiltonian system
$ ut=DδHδu, $
|
(1.2) |
where $ \mathcal D = (1-\partial_{xx})^{-1} \partial_x $ is a skew-adjoint operator, and $ \mathcal H $ is the Hamiltonian energy functional
$ H(t)=−12∫L0(κu2+u3+γuu2x)dx. $
|
(1.3) |
In this paper, we choose the periodic boundary condition
$ u(x,t)=u(x+L,t), x∈R, 0<t≤T. $
|
(1.4) |
Besides the energy conservation law $ \mathcal H(t) = \mathcal H(0) $, the dispersive system (1.1) further satisfies the following mass and momentum conservation laws, respectively,
$ I(t)=∫L0udx=I(0), $
|
(1.5) |
$ M(t)=12∫L0(u2+(ux)2)dx=M(0). $
|
(1.6) |
Over the past few decades, structure-preserving methods have come to the fore due to their superior capability for long-term computation over traditional numerical methods [12]. In Ref. [4], Cohen et al. first proposed multisymplectic finite difference schemes for the CH equation. Subsequently, their idea was generalized to solve generalized hyperelastic-rod wave equation [5]. In particular, the wavelet collocation method [30] and discontinuous Galerkin method [24] were further shown to be powerful tools for constructing multisymplectic schemes of the nonlinear system (1.1). In addition to the geometric structure, the nonlinear dispersive equation also satisfies three first integrals Eqs (1.3)–(1.6). It is well-known that the first integral plays a key role in the numerical analyses for conservative systems. In Ref. [18], Matsuo et al. first proposed an energy-preserving Galerkin scheme for the dispersive system (1.1), and some numerical experiments for the CH equation and Dai equation were also investigated. Moreover, numerical comparisons between the energy-preserving scheme and momentum-preserving scheme have been conducted for the CH equation [17]. Further discussions on energy-preserving Galerkin schemes can be found in Ref. [19]. Cohen et al. [5] then proposed a new energy-preserving scheme based on the discrete gradient approach. Later on, Gong et al. [9] combined the averaged vector field method and wavelet collocation method to design an energy-preserving wavelet collocation scheme, and some numerical comparisons with the multisymplectic wavelet collocation scheme [30] were also carried out. Recently, Brugnano et al. [1] developed an energy-conserving and space-time spectrally accurate method for Hamiltonian systems. However, almost all the mentioned schemes are fully implicit, which implies that the nonlinear iteration is inevitable at each time step, and thus it may be time consuming. In Ref. [8], Furihata and Matsuo first proposed a linearly implicit energy-preserving scheme for the CH model, in which only a linear system needs to be solved at each time step. Thus it is computationally much cheaper than that of the fully implicit counterparts. Another efficient strategy for constructing linearly implicit energy-preserving schemes of the nonlinear dispersive system (1.1) is based on the Kahan's method and the polarised discrete gradient method. For more details, please refer to Refs. [6,7]. Moreover, some linearly implicit momentum-preserving schemes can be founded in Refs. [13,16]. Nevertheless, to our best knowledge, the mentioned linearly implicit schemes are of only second-order accuracy in time at most.
It's worth noting that Yang et al. [27] coined the idea of invariant energy quadratization (IEQ) to develop energy stable schemes for solving the MBE model. Then Gong et al. [10] first generalized their idea to construct arbitrarily high-order energy stable algorithms for gradient flow models. Subsequently, the IEQ method [22] was further shown to be an efficient approach to design second-order and linearly implicit energy-preserving schemes for the CH model [14,15]. However, it is known that the issue with IEQ /EQ method is that it only preserves the modified energy. Recently, this has been partially addressed with a relaxation technique [29]. Moreover, for certain types of nonlinear systems, numerical schemes based on the IEQ idea have been shown to preserve the original energy law [28].
In this paper, we aim to develop a class of high-order, linearly implicit and energy-preserving schemes for the dispersive Eq (1.1). Based on the idea of the IEQ approach, the original model is firstly reformulated into an equivalent system with a quadratic energy by introducing an auxiliary variable. Inspired by the previous works [11,21], we further apply the prediction-correction strategy and the symplectic Partitioned Runge-Kutta (PRK) method [12,23] for temporal discretization, and the resulting semi-discrete system is linear, high-order accurate, and can preserve the quadratic energy of the reformulated system exactly. Finally, various numerical tests are addressed to verify the efficiency and accuracy of the proposed schemes.
In this section, we utilize the EQ idea to recast the model Eq (1.1) into an equivalent system, which possesses a modified energy function of the new variable. The new system can provide an elegant platform for developing high-order, linearly implicit and energy-preserving schemes. Here, we consider the system (1.1) in a finite domain $ [0, L]\times [0, T] $, and define the $ L^2 $ inner product as $ (f, g) = \int_0^L fg dx, \forall f, g \in L^2 ([0, L]) $.
Inspired by the EQ approach [27], we first introduce an auxiliary variable
$ q:=q(x,t)=−12(κu+u2+γu2x), $
|
and the Hamiltonian energy Eq (1.3) can be reformulated in a quadratic form rewritten as
$ H(t)=(u,q). $
|
(2.1) |
According to the energy variational principle, we further reformulate the system (1.2) to the following equivalent form
$ {ut=D(q−12κu−u2+γ(uxu)x),qt=−12κut−uut−γuxuxt, $
|
(2.2) |
with consistent initial conditions
$ u(x,0)=u0(x), q(x,0)=−12(κu0+u20+γ(u0)2x). $
|
It is readily to verify the energy conservation law. The integration by parts, together with the periodic boundary condition and Eq (2.2), yields
$ dHdt=(ut,q)+(u,qt)=(q,ut)+(u,−12κut−uut−γuxuxt)=(q−12κu−u2+γ(uxu)x,ut)=(q−12κu−u2+γ(uxu)x,D(q−12κu−u2+γ(uxu)x))=0, $
|
where the skew-symmetry of the operator $ \mathcal D $ was used.
In this section, we develop a class of high-precision, linearly implicit schemes for the reformulated system (2.2) by combing the prediction-correction approach and the PRK method, then show that the proposed schemes can exactly preserve the modified energy Eq (2.1). We here focus on temporal semi-discretization, and denote $ t_n = n\tau, n = 0, 1, 2, \cdots $, where $ \tau $ is the time step.
Inspired by the previous works in [11,21], we employ the prediction-correction strategy, together with the PRK method, to the system (2.2), and obtain the following prediction-correction method:
Scheme 3.1. Let $ b_i, \widehat b_i, a_{ij}, \widehat a_{ij} (i, j = 1, \cdots, s) $ be real numbers and let $ c_i = \sum_{j = 1}^sa_{ij}, \widehat c_i = \sum_{j = 1}^s \widehat a_{ij} $, we apply a $ s $-stage PRK method to the system (2.2) in time, and first compute the predicted values $ u_i^{n, *}, i = 1, \cdots, s $, as follows:
1. Prediction: for given $ (u^n, q^n) $, we set $ u_i^{n, 0} = u^n, q_i^{n, 0} = q^n $. Let $ M > 0 $ be a given integer, for $ m = 0 $ to $ M-1 $, we compute $ u_i^{n, m+1}, q_i^{n, m+1}, k_i^{n, m+1}, l_i^{n, m+1} $ by using
$ {un,m+1i=un+τs∑j=1aijkn,m+1j,kn,m+1i=D(qn,m+1i−12κun,m+1i−(un,mi)2+γ((un,mi)xun,mi)x),qn,m+1i=qn+τs∑j=1ˆaijln,m+1j,ln,m+1i=−12κkn,m+1i−un,mikn,mi−γ(un,mi)x(kn,mi)x,i=1,⋯,s. $
|
Then, for the predicted value, we set $ u_i^{n, *} = u_i^{n, M} $, followed by the following correction step:
2. Correction: we further compute the intermediate values $ k_i^n, l_i^n, u_i^n, q_i^n $ via
$ {uni=un+τs∑j=1aijknj, qni=qn+τs∑j=1ˆaijlnj,kni=D(qni−12κuni−un,∗iuni+γ((un,∗i)xuni)x),lni=−12κkni−un,∗ikni−γ(un,∗i)x(kni)x, $
|
(3.1) |
then $ (u^{n+1}, q^{n+1}) $ is updated via
$ un+1=un+τs∑i=1bikni, qn+1=qn+τs∑i=1ˆbilni. $
|
Theorem 3.1. If the coefficients of Scheme 3.1 satisfy
$ biˆaij+ˆbjaji=biˆbjfor i,j=1,⋯,s,bi=ˆbifor i=1,⋯,s, $
|
(3.2) |
then it preserves the following semi-discrete energy conservation law
$ En+1=En, n=0,1,⋯,N−1, $
|
(3.3) |
where $ E^n = (u^n, q^n) $.
Proof. It follows from the last two equalities of Eq (3.1) that
$ (kni,qni)+(uni,lni)=(qni,kni)+(uni,−12κkni−un,∗ikni−γ(un,∗i)x(kni)x)=(qni−12κuni−un,∗iuni+γ((un,∗i)xuni)x,kni)=(qni−12κuni−un,∗iuni+γ((un,∗i)xuni)x,D(qni−12κuni−un,∗iuni+γ((un,∗i)xuni)x))=0, i=1,⋯,s, $
|
where the skew-symmetry of $ \mathcal D $ was used. This, together with the first two equalities of Eq (3.1) and the condition (3.2), derives that
$ En+1−En=(un+1,qn+1)−(un,qn)=(un+τs∑i=1bikni,qn+τs∑i=1ˆbilni)−(un,qn)=τs∑i=1bi(kni,qn)+τs∑i=1ˆbi(un,lni)+τ2s∑i=1s∑j=1biˆbj(kni,lnj)=τs∑i=1bi(kni,qni−τs∑j=1ˆaijlnj)+τs∑i=1ˆbi(uni−τs∑j=1aijknj,lni)+τ2s∑i=1s∑j=1biˆbj(kni,lnj)=τs∑i=1bi(kni,qni)+τs∑i=1ˆbi(uni,lni)−τ2s∑i=1s∑j=1biˆaij(kni,lnj)−τ2s∑i=1s∑j=1ˆbiaij(knj,lni)+τ2s∑i=1s∑j=1biˆbj(kni,lnj)=τs∑i=1bi[(kni,qni)+(uni,lni)]+τ2s∑i=1s∑j=1(biˆbj−biˆaij−ˆbjaji)(kni,lnj)=0, $
|
which completes the proof.
Remark 3.1. If the coefficients of a partitioned Runge-Kutta method satisfy Eq (3.2), it conserves the quadratic invariant of the form $ Q(y, z) = y^{\top}Dz $ (see Ref. [12] for more details), thus the energy conservation law Eq (3.3) stands naturally. Moreover, the coefficients of PRK methods of order 4 and 6 are explicitly given in Tables 1 and 2 (please see [20]), respectively. Hereafter, the prediction-correction scheme 3.1, coupled with 4th- and 6th-order PRK methods, are denoted by 4th- and 6th-order LEQP1, respectively. Exchanging the order of Lobatto $ \rm III $A-$ \rm III $B pairs, we further derive 4th- and 6th-order LEQP2 methods.
![]() |
![]() |
Remark 3.2. In general, a numerical algorithm, which can preserve the corresponding energy exactly, is called an energy-preserving scheme. Thus, for the semi-discrete Scheme 3.1, it should be careful to choose a suitable spatial discretization, and we need to consider the following aspects: (i) It should preserve the skew-symmetry of the operator $ \mathcal D $; (ii) The accuracy of spatial discretization should be comparable to that of temporal counterpart. These, together with the periodic boundary condition (1.4), enlighten us to use the Fourier pseudo-spectral method for spatial discretization, which allows the application of FFT technique. In fact, the derived full discrete schemes also preserve the corresponding full discrete energy conservation law. The details are omitted here.
Remark 3.3. To our best knowledge, no theoretical result is found on the choice of iteration step $ M $. From our numerical experience, the $ 4th $-order LEQP methods can reach fourth-order accuracy by choosing $ M = 3 $, while the $ 6th $-order LEQP methods can reach sixth-order accuracy by choosing $ M = 5 $.
In this section, some numerical tests are carried out to investigate the accuracy, efficiency and invariants-preservation of the proposed schemes. For brevity, we only take the CH equation as a benchmark model for illustration purposes, and the corresponding prediction-correction schemes can be directly obtained by setting $ \kappa = 0, \gamma = 1 $ in Scheme 3.1. As shown above, the proposed scheme, which preserves the modified energy, could reach arbitrary high-order accuracy in time. Some comparisons would be made with the linearized Crank-Nicolson (IEQ-LCN) method [15] and energy-preserving Fourier pseudo-spectral (EPFP) method [9], where the wavelet collocation discretization in space is substituted by the standard Fourier pseudo-spectral method. In the following, the convergence rate in time is computed by the following formula
$ Rate=ln(error1/error2)ln(τ1/τ2), $
|
where $ \tau_l, error_l \ (l = 1, 2) $ are time steps and errors with time step $ \tau_l $, respectively.
Example 4.1 We first consider the CH equation as follows
$ ut−uxxt+3uux−2uxuxx−uuxxx=0, (x,t)∈Ω×(0,T],u(x,0)=u0(x), x∈Ω, $
|
(4.1) |
with the following conservation laws
$ M(t)=∫Ωudx, H(t)=−12∫Ω(u3+uu2x)dx, $
|
where $ u_0(x) = \text{ sin}(x), \Omega = [0, 2\pi] $, and $ \mathcal M, \mathcal H $ denote the mass and Hamiltonian energy, respectively. Moreover, the numerical "exact" solution is obtained by 6th-order LEQP1 scheme with small time step $ \tau = 0.001 $ and mesh size $ h = 2 \pi/128 $. Numerical results involving accuracy test are shown in Table 3. As is observed that the 4th- and 6th-order LEQP methods arrive at fourth- and sixth-order convergence rates in time, respectively. What's more, they perform obviously more accurate than the other second-order schemes, while the LEQP1 schemes seem a bit more accurate than the corresponding LEQP2 schemes for both 4th- and 6th-order cases.
$\tau=1/30$ | $\tau=1/60$ | $ \tau=1/100$ | $ \tau=1/120$ | ||
IEQ-LCN [15] | $|| e ||_{\infty, h}$ | 2.4481e-003 | 5.8420e-004 | 2.0830e-004 | 1.4438e-004 |
Rate | * | 2.06 | 2.02 | 2.01 | |
$|| e ||_h$ | 1.5831e-003 | 3.9367e-004 | 1.4180e-004 | 9.8508e-005 | |
Rate | * | 2.01 | 2.00 | 2.00 | |
EPFP [9] | $|| e ||_{\infty, h}$ | 4.8232e-004 | 1.2089e-004 | 4.3546e-004 | 3.0243e-005 |
Rate | * | 2.00 | 2.00 | 2.00 | |
$|| e ||_h$ | 6.1360e-004 | 1.5380e-004 | 5.5401e-005 | 3.8477e-005 | |
Rate | * | 2.00 | 2.00 | 2.00 | |
4th-order LEQP1 | $|| e ||_{\infty, h}$ | 1.9789e-007 | 1.2646e-008 | 1.6738e-009 | 8.1237e-010 |
Rate | * | 3.97 | 3.96 | 3.96 | |
$|| e ||_h$ | 8.6448e-008 | 5.7624e-009 | 7.7034e-010 | 3.7448e-010 | |
Rate | * | 3.91 | 3.94 | 3.96 | |
4th-order LEQP2 | $|| e ||_{\infty, h}$ | 4.8122e-007 | 3.1072e-008 | 4.0957e-009 | 1.9840e-009 |
Rate | * | 3.95 | 3.97 | 3.98 | |
$|| e ||_h$ | 2.2853e-007 | 1.4356e-008 | 1.8679e-009 | 9.0179e-010 | |
Rate | * | 3.99 | 3.99 | 3.99 | |
6th-order LEQP1 | $|| e ||_{\infty, h}$ | 2.9625e-010 | 4.1179e-012 | 1.8738e-013 | 6.2783e-014 |
Rate | * | 6.16 | 6.04 | 6.00 | |
$|| e ||_h$ | 2.3575e-010 | 2.6531e-012 | 1.1785e-013 | 3.9568e-014 | |
Rate | * | 6.47 | 6.09 | 5.99 | |
6th-order LEQP2 | $|| e ||_{\infty, h}$ | 8.9943e-010 | 1.0775e-011 | 4.5380e-013 | 1.4805e-013 |
Rate | * | 6.38 | 6.20 | 6.14 | |
$|| e ||_h$ | 4.9226e-010 | 5.0410e-012 | 1.9006e-013 | 6.0599e-014 | |
Rate | * | 6.61 | 6.42 | 6.26 |
In Figure 1a, we show the $ L^{\infty} $-norm solution error versus the execution time for different schemes, and it can be clearly seen that the proposed high-order methods perform more effective than other second-order counterparts. Subsequently, Figures 1b–d investigate the errors of invariants of different schemes in the long-time behaviour, where we choose $ \tau = 1/3000, N = 1/128 $.
As illustrated in Figure 1b, all schemes conserve the mass exactly. Though the proposed schemes can not preserve the exact Hamiltonian energy theoretically, the corresponding numerical errors can be preserved up to machine errors as demonstrated in Figure 1c. Moreover, Figure 1d illustrates that the proposed schemes preserve the quadratic energy precisely, which confirms the preceding theoretical analysis.
Example 4.2 In this example, we further apply the proposed method to simulate the three-peakon interaction of the CH equation with the initial condition [26]
$ u0(x)=ϕ1(x)+ϕ2(x)+ϕ3(x), $
|
where
$ ϕi(x)={cicosh(L/2)cosh(x−xi), |x−xi|≤L/2,cicosh(L/2)cosh(L−(x−xi)), |x−xi|>L/2, i=1,2,3. $
|
(4.2) |
The corresponding parameters are taken as $ c_1 = 2, c_2 = 1, c_3 = 0.8, x_1 = -5, x_2 = -3, x_3 = -1, L = 30 $, and the computational domain is $ \Omega = [0, L] $ with a periodic boundary condition.
In Figure 2a, we show the dynamics of three-peakon interaction computed by 4th-order LEQP1 scheme, and other LEQP schemes perform similarly. It is clearly seen that the moving peak interaction is resolved well. The taller wave overtakes the shorter counterparts, and afterwards they still keep their original shapes and velocities. Meanwhile, we study the invariants of the proposed schemes in long-time simulations. As demonstrated in Figures 2b–d, the proposed schemes preserve the mass and quadratic energy precisely, and perform more accurate than the IEQ-LCN method in terms of Hamiltonian energy.
In this paper, we develop a class of high-order, linearly implicit numerical algorithms, which can preserve the modified energy exactly, for the nonlinear dispersive equation. The proposed schemes, which can reach arbitrarily high-order accuracy in time, are easy to implement and computationally efficient. The key strategy lies in the application of the prediction-correction technique and the PRK method. Compared with the second-order linearly implicit structure-preserving schemes, the proposed methods perform more accurate and more efficient in longtime computations, and they can be directly applied to some relevant background problems, such as the Dai equation and regularized long wave equation in high dimensions.
Jin Cui's work is partially supported by the High Level Talents Research Foundation Project of Nanjing Vocational College of Information Technology (Grant No. YB20200906), the "Qinglan" Project of Jiangsu Province. Yayun Fu's work is partially supported by the National Natural Science Foundation of Henan Province (No. 222300420280).
The authors declare there is no conflict of interest.
[1] |
Grüntzig AR, Senning A, Siegenthaler WE (1979) Nonoperative dilatation of coronary-artery stenosis— percutaneous transluminal coronary angioplasty. New Engl J Med 301:61–68. doi: 10.1056/NEJM197907123010201
![]() |
[2] |
De Feyter PJ, de Jaegere PP, Serruys PW (1994) Incidence, predictors, and management of acute coronary occlusion after coronary angioplasty. Am Heart J 127: 643–651. doi: 10.1016/0002-8703(94)90675-0
![]() |
[3] | Serruys PW, Luijten HE, Beatt KJ, et al. (1988) Incidence of restenosis after successful coronary angioplasty: a time-related phenomenon. A quantitative angiographic study in 342 consecutive patients at 1, 2, 3, and 4 months. Circulation 77: 361–371. |
[4] |
Gruentzig AR, King SB 3rd, Schlumpf M, et al. (1987) Long-term follow-up after percutaneous transluminal coronary angioplasty. The early Zurich experience. New Engl J Med 316: 1127–1132. doi: 10.1056/NEJM198704303161805
![]() |
[5] |
Farooq V, Gogas BD, Serruys PW (2001) Restenosis: delineating the numerous causes of drug-eluting stent restenosis. Circ Cardiovasc Interv 4: 195–205. doi: 10.1080/14628840127766 ![]() |
[6] | Macander PJ, Agrawal SK, Roubin GS (1991) The Gianturco-Roubin balloon-expandable intracoronary flexible coil stent. J Invasive Cardiol 3: 85–94. |
[7] |
Schatz RA, Palmaz JC, Tio FO, et al. (1987) Balloon-expandable intracoronary stents in the adult dog. Circulation 76: 450–457. doi: 10.1161/01.CIR.76.2.450
![]() |
[8] |
Sigwart U, Gold S, Kaufman U, et al. (1988) Analysis of complications associated with coronary stenting (abstract). J Am Coll Cardiol 11: 66A. doi: 10.1016/0735-1097(88)90168-4
![]() |
[9] |
Sigwart U, Puel J, Mirkovitch V, et al. (1987) Intravascular stents to prevent occlusion and restenosis after transluminal angioplasty. New Engl J Med 316: 701–706. doi: 10.1056/NEJM198703193161201
![]() |
[10] |
Lansky AJ, Roubin GS, O’Shaughnessy CD, et al. (2000) Randomized comparison of GR-II stent and Palmaz-Schatz stent for elective treatment of coronary stenoses. Circulation 102: 1364–1368. doi: 10.1161/01.CIR.102.12.1364
![]() |
[11] |
Baim DS, Cutlip DE, Midei M, et al. (2001) Final results of a randomized trial comparing the MULTI-LINK stent with the Palmaz-Schatz stent for narrowings in native coronary arteries. Am J Cardiol 87: 157–162. doi: 10.1016/S0002-9149(00)01308-4
![]() |
[12] | Baim DS, Cutlip DE, O’Shaughnessy CD, et al. (2001) Final results of a randomized trial comparing the NIR stent to the Palmaz-Schatz stent for narrowings in native coronary arteries. Am J Cardiol 87: 152–156. |
[13] |
Colombo A, Hall P, Nakamura S, et al. (1995) Intracoronary stenting without anticoagulation accomplished with intravascular ultrasound guidance. Circulation 91: 1676–1688. doi: 10.1161/01.CIR.91.6.1676
![]() |
[14] | Duckers HJ, Nabel EG, Serruys PW, editors. (2007) Essentials of restenosis: for the interventional cardiologist. Totowa, NJ: Humana Press. |
[15] |
Tada T, Byrne RA, Simunovic I, et al. (2013) Risk of stent thrombosis among bare-metal stents, first-generation drug-eluting stents, and second-generation drug-eluting stents: results from a registry of 18,334 patients. JACC Cardiovasc Interv 6: 1267–1274. doi: 10.1016/j.jcin.2013.06.015
![]() |
[16] |
Daemen J, Wenaweser P, Tsuchida K, et al. (2007) Early and late coronary stent thrombosis of sirolimus-eluting and paclitaxel-eluting stents in routine clinical practice: data from a large two-institutional cohort study. Lancet 369: 667–678. doi: 10.1016/S0140-6736(07)60314-6 ![]() |
[17] | Wenaweser P, Daemen J, Zwahlen M, et al. (2008) Incidence and correlates of drug-eluting stent thrombosis in routine clinical practice. 4-year results from a large 2-institutional cohort study. J Am Coll Cardiol 52: 1134–1140. |
[18] | Sarno G, Lagerqvist B, Fröbert O, et al. (2012) Lower risk of stent thrombosis and restenosis with unrestricted use of “new-generation” drug-eluting stents: A report from the nationwide Swedish Coronary Angiography and Angioplasty Registry (SCAAR). Eur Heart J 33: 606–613. |
[19] |
Togni M, Windecker S, Cocchia R, et al. (2005) Sirolimus-eluting stents associated with paradoxic coronary vasoconstriction. J Am Coll Cardiol 46: 231–236. doi: 10.1016/j.jacc.2005.01.062
![]() |
[20] | Hofma SH, Van Der Giessen WJ, Van Dalen BM, et al. (2006) Indication of long-term endothelial dysfunction after sirolimus-eluting stent implantation. Eur Heart J 27: 166–170. |
[21] |
Onuma Y, Serruys PW (2011) Bioresorbable scaffold: the advent of a new era in percutaneous coronary and peripheral revascularization? Circulation 123: 779–797. doi: 10.1161/CIRCULATIONAHA.110.971606
![]() |
[22] | Wayangankar SA, Ellis SG (2015) Bioresorbable stents: is this where we are headed? Prog Cardiovasc Dis 26. |
[23] | Van der Giessen WJ, Slager CJ, van Beusekom HM, et al. (1992) Development of a polymer endovascular prosthesis and its implantation in porcine arteries. J Invasive Cardiol 5: 175–185. |
[24] |
Van der Giessen WJ, Lincoff AM, Schwartz RS, et al. (1996) Marked inflammatory sequelae to implantation of biodegradable and nonbiodegradable polymers in porcine coronary arteries. Circulation 94: 1690–1697. doi: 10.1161/01.CIR.94.7.1690
![]() |
[25] |
Lincoff AM, Furst JG, Ellis SG, et al. (1997) Sustained local delivery of dexamethasone by a novel intravascular eluting stent to prevent restenosis in the porcine coronary injury model. J Am Coll Cardiol 29: 808–816. doi: 10.1016/S0735-1097(96)00584-0
![]() |
[26] | Yamawaki T, Shimokawa H, Kozai T, et al. (1998) Intramural delivery of a specific tyrosine kinase inhibitor with biodegradable stent suppresses the restenotic changes of the coronary artery in pigs in vivo. J Am Coll Cardiol 32: 780–786. |
[27] | Iqbal J, Onuma Y, Ormiston J, et al. (2013) Bioresorbable scaffolds: rationale, current status, challenges, and future. Eur Heart J 35: 765–776. |
[28] |
Tamai H, Igaki K, Kyo E, et al. (2000) Initial and 6-month results of biodegradable poly-l-lactic acid coronary stents in humans. Circulation 102: 399–404. doi: 10.1161/01.CIR.102.4.399
![]() |
[29] | Serruys PW, Ormiston JA, Onuma Y, et al. (2009) A bioabsorbable everolimus-eluting coronary stent system (ABSORB): 2-year outcomes and results from multiple imaging methods. Lancet 373: 897–910. |
[30] |
Verheye S, Ormiston JA, Stewart J, et al. (2014) A next-generation bioresorbable coronary scaffold system: from bench to first clinical evaluation: 6- and 12-month clinical and multimodality imaging results. JACC Cardiovasc Interv 7: 89–99. doi: 10.1016/j.jcin.2013.07.007
![]() |
[31] | Abizaid A (2014) DESolve NX Trial clinical and imaging results. Presented at: Transcatheter Cardiovascular Therapeutics (TCT). Washington, DC; September 13-17. |
[32] | Costa RA (2012) REVA ReZolve clinical program update. Presented at: Transcatheter Cardiovascular Therapeutics (TCT). Miami Beach, FL; October 22-28. |
[33] | Muller D (2014) ReZolve 2: Bioresorbable coronary scaffold clinical program update. Presented at: EuroPCR. Paris, France; May 20-23. |
[34] | Erbel R, Di Mario C, Bartunek J, et al. (2007) Temporary scaffolding of coronary arteries with bioabsorbable magnesium stents: a prospective, non-randomised multicentre trial. Lancet 369: 1869–1875. |
[35] | Haude M, Erbel R, Erne P, et al. (2013) Safety and performance of the drug-eluting absorbable metal scaffold (DREAMS) in patients with de-novo coronary lesions: 12 month results of the prospective, multicentre, first-in-man BIOSOLVE-I trial. The Lancet 381: 836–844. |
[36] | Waksman R (2014) Long-term clinical data of the BIOSOLVE-I study with the paclitaxel-eluting absorbable magnesium scaffold (DREAMS 1st generation) and multi- modality imaging analysis. |
[37] |
Onuma Y, Dudek D, Thuesen L, et al. (2013) Five-year clinical and functional multislice computed tomography angiographic results after coronary implantation of the fully resorbable polymeric everolimus-eluting scaffold in patients with de novo coronary artery disease: the ABSORB cohort A trial. JACC Cardiovasc Interv 6: 999–1009. doi: 10.1016/j.jcin.2013.05.017
![]() |
[38] |
Serruys PW, Onuma Y, Garcia-Garcia HM, et al. (2014) Dynamics of vessel wall changes following the implantation of the Absorb everolimus-eluting bioresorbable vascular scaffold: a multi-imaging modality study at 6, 12, 24 and 36 months. EuroIntervention 9: 1271–1284. doi: 10.4244/EIJV9I11A217
![]() |
[39] | Serruys PW, Chevalier B, Dudek D, et al. (2015) A bioresorbable everolimus-eluting scaffold versus a metallic everolimus-eluting stent for ischaemic heart disease caused by de-novo native coronary artery lesions (ABSORB II): an interim 1-year analysis of clinical and procedural secondary outcomes from. Lancet 385: 43–54. |
[40] |
Ellis SG, Kereiakes DJ, Metzger DC, et al. (2015) Everolimus-eluting bioresorbable scaffolds for coronary artery disease. N Engl J Med 373: 1905–1915. doi: 10.1056/NEJMoa1509038
![]() |
[41] |
Sabaté M, Windecker S, Iñiguez A, et al. (2015) Everolimus-eluting bioresorbable stent vs. durable polymer everolimus-eluting metallic stent in patients with ST-segment elevation myocardial infarction: results of the randomized ABSORB ST-segment elevation myocardial infarction-TROFI II trial. Eur Heart J 23. |
[42] | Gao R, Yang Y, Han Y, et al. (2015) Bioresorbable vascular scaffolds versus metallic stents in patients with coronary artery disease: ABSORB China trial. J Am Coll Cardiol 25. |
[43] |
Capodanno D, Gori T, Nef H, et al. (2015) Percutaneous coronary intervention with everolimus-eluting bioresorbable vascular scaffolds in routine clinical practice: early and midterm outcomes from the European multicentre GHOST-EU registry. EuroIntervention 10: 1144–1153. doi: 10.4244/EIJY14M07_11 ![]() |
[44] |
Palmerini T, Biondi-Zoccai G, Della Riva D, et al. (2012) Stent thrombosis with drug-eluting and bare-metal stents: evidence from a comprehensive network meta-analysis. Lancet 379: 1393–1402. doi: 10.1016/S0140-6736(12)60324-9
![]() |
[45] | Abizaid A, Costa JR, Bartorelli AL, et al. (2015) The ABSORB EXTEND study: preliminary report of the twelve-month clinical outcomes in the first 512 patients enrolled. EuroIntervention 10: 1396–1401. |
[46] |
Džavík V, Colombo A (2014) The Absorb bioresorbable vascular scaffold in coronary bifurcations: insights from bench testing. JACC Cardiovasc Interv 7: 81–88. doi: 10.1016/j.jcin.2013.07.013
![]() |
[47] | Rundeken MJ, Hassell MECJ, Kraak RP, et al. (2014) Treatment of coronary bifurcation lesions with the Absorb bioresorbable vascular scaffold in combination with the Tryton dedicated coronary bifurcation stent: evaluation using two- and three-dimensional optical coherence tomography. EuroIntervention 30. |
[48] | Capranzano P, Gargiulo G, Capodanno D, et al. (2014) Treatment of coronary bifurcation lesions with bioresorbable vascular scaffolds. Minerva Cardioangiol 62: 229–234. |
[49] |
Kawamoto H, Latib A, Ruparelia N, et al. (2015) Clinical outcomes following bioresorbable scaffold implantation for bifurcation lesions: overall outcomes and comparison between provisional and planned double stenting strategy. Catheter Cardiovasc Interv 86: 644–652. doi: 10.1002/ccd.26045
![]() |
[50] | Cortese B, Silva Orrego P, Sebik R, et al. (2014) Biovascular scaffolding of distal left main trunk: experience and follow up from the multicenter prospective RAI registry (Registro Italiano Absorb). Int J Cardiol 177: 497–499. |
[51] | Panoulas VF, Miyazaki T, Sato K, et al. (2015) Procedural outcomes of patients with calcified lesions treated with bioresorbable vascular scaffolds. EuroIntervention 30. |
[52] | Vaquerizo B, Barros A, Pujadas S, et al. (2014) Bioresorbable everolimus-eluting vascular scaffold for the treatment of chronic total occlusions: CTO-ABSORB pilot study. EuroIntervention 16. |
[53] |
Diletti R, Karanasos A, Muramatsu T, et al. (2014) Everolimus-eluting bioresorbable vascular scaffolds for treatment of patients presenting with ST-segment elevation myocardial infarction: BVS STEMI first study. Eur Heart J 35: 777–786. doi: 10.1093/eurheartj/eht546
![]() |
[54] | Kocka V, Maly M, Tousek P, et al. (2014) Bioresorbable vascular scaffolds in acute ST-segment elevation myocardial infarction: a prospective multicentre study “Prague 19.” Eur Heart J 35: 787–794. |
[55] |
Ielasi A, Cortese B, Varricchio A, et al. (2015) Immediate and midterm outcomes following primary PCI with bioresorbable vascular scaffold implantation in patients with ST-segment myocardial infarction: insights from the multicentre “Registro ABSORB Italiano” (RAI registry). EuroIntervention 11: 157–162. doi: 10.4244/EIJY14M10_11
![]() |
[56] | Brugaletta S, Gori T, Low AF, et al. (2015) Absorb bioresorbable vascular scaffold versus everolimus-eluting metallic stent in ST-segment elevation myocardial infarction: 1-year results of a propensity score matching comparison: the BVS-EXAMINATION Study (bioresorbable vascular scaffold-a clinical evaluation of everolimus eluting coronary stents in the treatment of patients with ST-segment elevation myocardial infarction). JACC Cardiovasc Interv 8: 189–197. |
[57] | Ong PJ, Jafary FH, Ho HH (2013) “First-in-man” use of bioresorbable vascular scaffold in saphenous vein graft. EuroIntervention 9: 165. |
[58] | Deora S, Shah S, Patel T (2014) First-in-man implantation of bioresorbable vascular scaffold in left internal mammary artery graft. J Invasive Cardiol 26: 92–93. |
[59] |
Moscarella E, Varricchio A, Stabile E, et al. (2015) Bioresorbable vascular scaffold implantation for the treatment of coronary in-stent restenosis: results from a multicenter Italian experience. Int J Cardiol 199: 366–372. doi: 10.1016/j.ijcard.2015.07.002
![]() |
[60] |
Tamburino C, Latib A, van Geuns R-J, et al. (2015) Contemporary practice and technical aspects in coronary intervention with bioresorbable scaffolds: a European perspective. EuroIntervention 11: 45–52. doi: 10.4244/EIJY15M01_05
![]() |
[61] | Azzalini L, L’Allier PL (2015) Bioresorbable vascular scaffold thrombosis in an all-comer patient population: single center experience. J Invasive Cardiol 27: 85–92. |
[62] | Azzalini L, Al-Hawwas M, L’Allier PL (2015) Very late bioresorbable vascular scaffold thrombosis: a new clinical entity. Eurointervention 11: e1–2. |
[63] |
Karanasos A, van Geuns R-J, Zijlstra F, et al. (2014) Very late bioresorbable scaffold thrombosis after discontinuation of dual antiplatelet therapy. Eur Heart J 35: 1781. doi: 10.1093/eurheartj/ehu031
![]() |
[64] |
Karanasos A, Van Mieghem N, van Ditzhuijzen N, et al. (2015) Angiographic and optical coherence tomography insights into bioresorbable scaffold thrombosis: single-center experience. Circ Cardiovasc Interv 8: e002369. doi: 10.1161/CIRCINTERVENTIONS.114.002369
![]() |
[65] |
Nakatani S, Onuma Y, Ishibashi Y, et al. (2015) Early (before 6 months), late (6-12 months) and very late (after 12 months) angiographic scaffold restenosis in the ABSORB Cohort B trial. EuroIntervention 10: 1288–1298. doi: 10.4244/EIJV10I11A218
![]() |
[66] |
Longo G, Granata F, Capodanno D, et al. (2015) Anatomical features and management of bioresorbable vascular scaffolds failure: A case series from the GHOST registry. Catheter Cardiovasc Interv 85: 1150–1161. doi: 10.1002/ccd.25819
![]() |
[67] |
Felix C, Everaert B, Jepson N, et al. (2015) Treatment of bioresorbable scaffold failure. EuroIntervention 11: V175–180. doi: 10.4244/EIJV11SVA42
![]() |
[68] | Cassese S, Byrne RA, Ndrepepa G, et al. (2015) Everolimus-eluting bioresorbable vascular scaffolds versus everolimus-eluting metallic stents: a meta-analysis of randomised controlled trials. The Lancet 6736: 1–8. |
![]() |
![]() |
$\tau=1/30$ | $\tau=1/60$ | $ \tau=1/100$ | $ \tau=1/120$ | ||
IEQ-LCN [15] | $|| e ||_{\infty, h}$ | 2.4481e-003 | 5.8420e-004 | 2.0830e-004 | 1.4438e-004 |
Rate | * | 2.06 | 2.02 | 2.01 | |
$|| e ||_h$ | 1.5831e-003 | 3.9367e-004 | 1.4180e-004 | 9.8508e-005 | |
Rate | * | 2.01 | 2.00 | 2.00 | |
EPFP [9] | $|| e ||_{\infty, h}$ | 4.8232e-004 | 1.2089e-004 | 4.3546e-004 | 3.0243e-005 |
Rate | * | 2.00 | 2.00 | 2.00 | |
$|| e ||_h$ | 6.1360e-004 | 1.5380e-004 | 5.5401e-005 | 3.8477e-005 | |
Rate | * | 2.00 | 2.00 | 2.00 | |
4th-order LEQP1 | $|| e ||_{\infty, h}$ | 1.9789e-007 | 1.2646e-008 | 1.6738e-009 | 8.1237e-010 |
Rate | * | 3.97 | 3.96 | 3.96 | |
$|| e ||_h$ | 8.6448e-008 | 5.7624e-009 | 7.7034e-010 | 3.7448e-010 | |
Rate | * | 3.91 | 3.94 | 3.96 | |
4th-order LEQP2 | $|| e ||_{\infty, h}$ | 4.8122e-007 | 3.1072e-008 | 4.0957e-009 | 1.9840e-009 |
Rate | * | 3.95 | 3.97 | 3.98 | |
$|| e ||_h$ | 2.2853e-007 | 1.4356e-008 | 1.8679e-009 | 9.0179e-010 | |
Rate | * | 3.99 | 3.99 | 3.99 | |
6th-order LEQP1 | $|| e ||_{\infty, h}$ | 2.9625e-010 | 4.1179e-012 | 1.8738e-013 | 6.2783e-014 |
Rate | * | 6.16 | 6.04 | 6.00 | |
$|| e ||_h$ | 2.3575e-010 | 2.6531e-012 | 1.1785e-013 | 3.9568e-014 | |
Rate | * | 6.47 | 6.09 | 5.99 | |
6th-order LEQP2 | $|| e ||_{\infty, h}$ | 8.9943e-010 | 1.0775e-011 | 4.5380e-013 | 1.4805e-013 |
Rate | * | 6.38 | 6.20 | 6.14 | |
$|| e ||_h$ | 4.9226e-010 | 5.0410e-012 | 1.9006e-013 | 6.0599e-014 | |
Rate | * | 6.61 | 6.42 | 6.26 |
![]() |
![]() |
$\tau=1/30$ | $\tau=1/60$ | $ \tau=1/100$ | $ \tau=1/120$ | ||
IEQ-LCN [15] | $|| e ||_{\infty, h}$ | 2.4481e-003 | 5.8420e-004 | 2.0830e-004 | 1.4438e-004 |
Rate | * | 2.06 | 2.02 | 2.01 | |
$|| e ||_h$ | 1.5831e-003 | 3.9367e-004 | 1.4180e-004 | 9.8508e-005 | |
Rate | * | 2.01 | 2.00 | 2.00 | |
EPFP [9] | $|| e ||_{\infty, h}$ | 4.8232e-004 | 1.2089e-004 | 4.3546e-004 | 3.0243e-005 |
Rate | * | 2.00 | 2.00 | 2.00 | |
$|| e ||_h$ | 6.1360e-004 | 1.5380e-004 | 5.5401e-005 | 3.8477e-005 | |
Rate | * | 2.00 | 2.00 | 2.00 | |
4th-order LEQP1 | $|| e ||_{\infty, h}$ | 1.9789e-007 | 1.2646e-008 | 1.6738e-009 | 8.1237e-010 |
Rate | * | 3.97 | 3.96 | 3.96 | |
$|| e ||_h$ | 8.6448e-008 | 5.7624e-009 | 7.7034e-010 | 3.7448e-010 | |
Rate | * | 3.91 | 3.94 | 3.96 | |
4th-order LEQP2 | $|| e ||_{\infty, h}$ | 4.8122e-007 | 3.1072e-008 | 4.0957e-009 | 1.9840e-009 |
Rate | * | 3.95 | 3.97 | 3.98 | |
$|| e ||_h$ | 2.2853e-007 | 1.4356e-008 | 1.8679e-009 | 9.0179e-010 | |
Rate | * | 3.99 | 3.99 | 3.99 | |
6th-order LEQP1 | $|| e ||_{\infty, h}$ | 2.9625e-010 | 4.1179e-012 | 1.8738e-013 | 6.2783e-014 |
Rate | * | 6.16 | 6.04 | 6.00 | |
$|| e ||_h$ | 2.3575e-010 | 2.6531e-012 | 1.1785e-013 | 3.9568e-014 | |
Rate | * | 6.47 | 6.09 | 5.99 | |
6th-order LEQP2 | $|| e ||_{\infty, h}$ | 8.9943e-010 | 1.0775e-011 | 4.5380e-013 | 1.4805e-013 |
Rate | * | 6.38 | 6.20 | 6.14 | |
$|| e ||_h$ | 4.9226e-010 | 5.0410e-012 | 1.9006e-013 | 6.0599e-014 | |
Rate | * | 6.61 | 6.42 | 6.26 |