S Dumont

Imen Manoubi
email: imen.manoubi@yahoo.fr

Numerical analysis of a water wave model with a nonlocal viscous dispersive term using the diffusive approach

Keywords: decay rate, diffusive realization, fractional derivatives, nonlocal viscous model, time splitting, water waves

In this paper, we numerically study the water wave model with a nonlocal viscous term

where

ds is the Riemann-Liouville half-order derivative in time. We propose and compare different numerical schemes based on the diffusive realizations of fractional operators.

INTRODUCTION

State of the art

Modeling the effects of viscosity in the propagation of long waves is a challenging issue. In the last decade, Liu-Orfila 1 and Dutykh-Dias 2 have independently derived mathematical models in 2D and in 3D where the viscosity is modeled with nonlocal in time operators. In the case of a mono-dimension wave, the model is reduced as follows:

u t + u x + 𝛽u xxx + √ 𝜈 √ 𝜋 ∫ t 0 u t (s) √ t -s ds + uu x = 𝛼u xx , (1)
where We denote as 𝛽, 𝜈,and𝛼 the parameters dedicated to balance or unbalance the effects of viscosity and dispersion against nonlinear effects. Some numerical results of the decay rate of a solution of (1) are performed in Chen et al [START_REF] Chen | Decay of solutions to a water wave model with a nonlocal viscous dispersive term[END_REF] and Dumont and Duval. [START_REF] Dumont | Numerical investigation of asymptotical properties of solutions to models for water waves with non local viscosity[END_REF] Also, in the recent work, [START_REF] Manoubi | Theoretical and numerical analysis of the decay rate of solutions to a water wave model with a nonlocal viscous dispersive term with Riemann-Liouville half derivative[END_REF] the author has considered the model (1) where she used the Riemann-Liouville half derivative instead of that of Caputo, namely,

u t + u x + 𝛽u xxx + √ 𝜈 √ 𝜋 𝜕 𝜕t ∫ t 0 u(s) √ t -s ds + uu x = 𝛼u xx . ( 2
)
1

Comparing with the problem (1), she improved the estimates on the decay rates of the solution. Particularly, she proved that the local and global existence result and decay estimates for the integro-differential equation 2 when 𝛽 = 0, 𝜈 = 𝛼 = 1 supplemented with the initial condition u 0 ∈ L 1 (R)∩L 2 (R). Precisely, she stated the following theorem.

Theorem 1. Let u 0 ∈ L 2 (R), [START_REF] Manoubi | Theoretical and numerical analysis of the decay rate of solutions to a water wave model with a nonlocal viscous dispersive term with Riemann-Liouville half derivative[END_REF] then there exists a unique local solution u ∈ C([0, T); L 2 x (R)) of (2). Moreover, for u 0 ∈ L 1 (R)∩L 2 (R), there exists a positive constant C 0 > 0 that depends on u 0 such that if || u 0 || L 1 (R) is small enough, there exists a unique global solution u ∈ C(R + ; L 2

x (R)) ∩ C 1∕2 (R + ; H -2 x (R)) of (2) given by u(t, x)=[K RL (t, •) ⋆ u 0 ](x)-N ⊛ u 2 (t, x), (3) where K RL and N are given by

K RL (t, x)= 1 2 √ 𝜋t e -x 2 4t e -x - ( 1 - 1 2 ∫ +∞ 0 e -𝜇 2 4t -𝜇|x| 2t -𝜇 2 d𝜇
) and

N(t, x)= 1 4 √ 𝜋t 𝜕 x ( e -x 2 4t e -x - ( 1 - 1 2 ∫ +∞ 0 e -𝜇 2 4t -𝜇|x| 2 -𝜇 2 d𝜇
))

with x -= |x|-x 2 = max(-x, 0); ⋆ represents the usual convolution product and ⊛ is the time-space convolution product defined by v ⊛ w(t, x)= ∫ t 0 ∫ R v(ts, x -𝑦)w(s,𝑦)dsd𝑦 whenever the integrals make sense. In addition, we have the following estimate:

max(t 1∕4 , t 3∕4 )|| u(t, •)|| L 2 x (R) + max(t 1∕2 , t)|| u(t, •)|| L ∞ x (R) ≤ C 0 . ( 4
)
The proof of this theorem is presented in Manoubi. [START_REF] Manoubi | Theoretical and numerical analysis of the decay rate of solutions to a water wave model with a nonlocal viscous dispersive term with Riemann-Liouville half derivative[END_REF] The approximation of time-fractional operators has received a lot of interest during last decades for their wide application in fluid and solid mechanics and visco-elasticity. The formulation of a numerical stable scheme is crucial but also a difficult question because of the nonlocal feature of such operators. The classical methods used in the literature [START_REF] Galucio | The G 𝛼 -scheme for approximation of fractional derivatives: application to the dynamics of dissipative systems[END_REF][START_REF] Galucio | An adaptation of the Gear scheme for fractional derivatives[END_REF][START_REF] Koh | Application of fractional derivatives to seismic analysis of base-isolated models[END_REF][START_REF] Li | A fast time stepping method for evaluating fractional integrals[END_REF][START_REF] Shokooh | A comparison of numerical methods applied to a fractional model of damping materials[END_REF][START_REF] Zhang | Numerical algorithm for dynamic problems involving fractional operators[END_REF] consist on the approximation of these fractional operators using either convolution integrals or the so-called Gear scheme for fractional operators. However, these methods are expensive in terms of run time evaluations and of the memory required to store the function values, especially for large time simulations. Also, Zang and Xu proposed in Zhang and Xu [START_REF] Zhang | Finite difference/spectral approximations to a water wave model with a nonlocal viscous term[END_REF] numerical schemes based on finite difference/spectral approximations to a water wave model. In particular, they used the known (2 -𝛼) order scheme for the 𝛼-order fractional derivative and a mixed linearization for the nonlinear term. They proved the unconditionally stability and some convergence rates for the schemes. Moreover, in their work, [START_REF] Liu | Finite difference/finite element method for a nonlinear time-fractional fourth-order reaction-diffusion problem[END_REF] Liu et al constructed finite difference/finite element method for a nonlinear time-fractional fourth-order reaction-diffusion problem. They proved unconditionally stability of the schemes and some a priori estimates of L 2 norm with optimal order of convergence. Independently, a number of authors suggested a method for the direct approximation of the fractional operators. This method, called diffusive representation, [START_REF] Montseny | Représentation Diffusive[END_REF] is an operator theory approach and is devoted to sech pseudo-differential equations and time nonlocal problems. This approach was developed by Montseny, [START_REF] Montseny | Diffusion monodimensionnelle et intégration d'ordre 1∕2[END_REF] Montseny et al, [START_REF] Montseny | Modèle simple d'amortisseur viscoélastique. Application à une corde vibrante[END_REF][START_REF] Montseny | Optimal models of fractional integrators and application to systems with fading memory[END_REF] and Staffans. [START_REF] Staffans | Well-posedness and stabilizability of a viscoelastic equation in energy space[END_REF] Different applications of this approach can be found in other studies. [START_REF] Audounet | A threshold phenomenon in the propagation of a point-source initiated flame[END_REF][START_REF] Helie | Diffusive representations for the analysis and simulation of flared acoustic pipes with visco-thermal losse[END_REF][START_REF] Matignon | Asymptotic stability of linear conservative systems when coupled with diffusive systems[END_REF][START_REF] Montseny | Diffusive representation for pseudo-differentially damped nonlinear systems[END_REF] The main idea of this method is to replace the nonlocal operator by a linear differential equation. The resulting diffusive model is infinite dimensional, but local in time, which has many advantages from the mathematical and numerical points of view.

Diffusive formulation of the model

In the literature, there are several diffusive realizations of the Riemann-Liouville half-order derivative D 1/2 u(t). We recall in the following some of these formulations. First, the diagonal form of the diffusive realization of

D 1/2 u(t) is given for all t > 0by ⎧ ⎪ ⎨ ⎪ ⎩ 𝜕 t 𝜓(t,𝜎)=-𝜎𝜓(t,𝜎)+u(t),𝜓 (0,𝜎)=0,𝜎 ≥ 0, D 1∕2 u(t)= ∫ +∞ 0 1 𝜋 √ 𝜎 𝜕 t 𝜓(t,𝜎)d𝜎. ( 5
)
Second, the partial differential equation (PDE) form of the diffusive realization of D 1/2 u(t) is given for all t > 0by ⎧ ⎪ ⎨ ⎪ ⎩ 𝜕 t Φ(t,𝑦)=Φ 𝑦𝑦 (t,𝑦)+u(t) ⊗𝛿 𝑦=0 , Φ(0,𝑦)=0,𝑦 ∈ R,

D 1∕2 u(t)=2 <𝛿 𝑦=0 ,𝜕 t Φ(t,𝑦)>  ′ , = 2 d dt Φ(t, 0), (6)
where 𝛿 y=0 is the Dirac delta function at y = 0andu(t) ⊗𝛿 y=0 is the tensorial product in the distributions sense of the applications t  → u(t) and y  → 𝛿 y=0 . Finally the diffusive realization of D 1/2 u(t) used in the remaining of this article is given for all t > 0by

⎧ ⎪ ⎨ ⎪ ⎩ 𝜕 t 𝜙(t,𝜎)=-𝜎 2 𝜙(t,𝜎)+ 2 𝜋 u(t),𝜙 (0,𝜎)=0,𝜎 ≥ 0, D 1∕2 u(t)= ∫ +∞ 0 ( 2
𝜋 u(t)-𝜎 2 𝜙(t,𝜎) ) d𝜎. ( 7
)
We note that the author has used the diffusive realizations ( 6) and (7) in her PhD thesis to study mathematically the integro-differential equation 2 when 𝛽 = 0, 𝜈 = 𝛼 = 1 using 2 different methods. For more details, we refer the readers to Manoubi. [START_REF] Manoubi | Modèle visqueux asymptotique pour la propagation d'une onde dans un canal[END_REF] In this paper, we study numerically the integro-differential equation 2 using the diffusive approach. Hence, we extend the diffusive realization (7) of the Riemann-Liouville half-order derivative as follows:

⎧ ⎪ ⎨ ⎪ ⎩ 𝜕 t 𝜓(t, x,𝜎)=-𝜎 2 𝜓(t, x,𝜎)+ 2 𝜋 u(t, x),𝜓 (0, x,𝜎)=0, ∀𝜎 ≥ 0, D 1∕2 u(t, x)= ∫ +∞ 0 ( 2
𝜋 u(t, x)-𝜎 2 𝜓(t, x,𝜎) ) d𝜎. ( 8
)
Then, the integro-differential equation 2 can be written as a PDE-ordinary differential equation (ODE) coupled system as follows:

⎧ ⎪ ⎨ ⎪ ⎩ u t (t, x)+ √ 𝜈 ∫ +∞ 0 ( 2 𝜋 u(t, x)-𝜎 2 𝜓(t, x,𝜎) ) d𝜎 = 𝛼u xx (t, x) -u x (t, x)-𝛽u xxx (t, x)-𝛾u(t, x)u x (t, x), t > 0, x ∈ R, 𝜕 t 𝜓(t, x,𝜎)=-𝜎 2 𝜓(t, x,𝜎)+ 2 𝜋 u(t, x), t > 0, x ∈ R,𝜎 ≥ 0, (9)
supplemented with the initial conditions

∀x ∈ R, ∀𝜎 ≥ 𝜓(0, x,𝜎)=0, ∀x ∈ R, u(0, x)=u 0 (x).
The remaining of this article is as follows. Different quadrature methods are presented in Section 2 to approximate the Riemann-Liouville half-order derivative. In Section 3, we present a first numerical scheme associated to the system (9) where we used the quadrature methods followed by several numerical simulations and discussion. In Section 4, we present a second numerical scheme for (9) using a splitting procedure. Numerical results are proposed in order to check the efficiency of the numerical scheme. A dispersion analysis is presented in Section 5.

QUADRATURE METHODS

To approximate the Riemann-Liouville half-order derivative in (2), we need to approximate the generalized integral in (9).

To this end, we use a quadrature formula with N points.W enotebyw i the weights and by 𝜎 i the nodes (or abscissae) of the appropriate quadrature method used in the approximation. We get

D 1∕2 u(t, x)≃ N ∑ i=1 w i ( 2 𝜋 u(t, x)-𝜎 2 i 𝜓(t, x,𝜎 i ) ) = N ∑ i=1 w i ( 2 𝜋 u(t, x)-𝜎 2 i 𝜓 i (t, x)
) ,

where 𝜓 i verifies the ODE (8). Hence, system ( 9) is approximated as follows:

⎧ ⎪ ⎪ ⎨ ⎪ ⎪ ⎩ u t (t, x)+ √ 𝜈 N ∑ i=1 w i ( 2 𝜋 u(t, x)-𝜎 2 i 𝜓 i (t, x) ) = 𝛼u xx (t, x) -u x (t, x)-𝛽u xxx (t, x)-𝛾u(t, x)u x (t, x), t > 0, x ∈ R, 𝜕 t 𝜓 i (t, x)=-𝜎 2 i 𝜓 i (t, x)+ 2 𝜋 u(t, x), t > 0, x ∈ R, i = 0, •••, N, (10)
supplemented with the initial conditions

∀x ∈ R, ∀i = 0, •••, N 𝜓 i (0, x)=0, ∀x ∈ R, u(0, x)=u 0 (x).
In the sequel, we are interested in the choice of the 2N coefficients w i and 𝜎 i of the diffusive representation given in (10). These coefficients aim to approach the improper integrals in the form

∫ +∞ 0 𝜓(𝜎)d𝜎 ≃ N ∑ i=1 w i 𝜓(𝜎 i ). (11)
For seek of convenience, we drop here the variables t and x. The choice of the coefficients w i and 𝜎 i is an important issue, because it effects directly the accuracy of the method and the efficiency of the approximation. We note that many methods developed in the literature are especially based on the Gauss quadrature. In the sequel, we present some quadrature methods that will be used in the remaining of this article.

The Gauss-Laguerre quadrature

A first method consists in approximating the generalized integral (11) using the standard Gauss-Laguerre quadrature over the interval [0, +∞[ with the weight function w(𝜎)=e -𝜎 (see Abramowitz and Stegun [START_REF] Abramowitz | Handbook of mathematical functions with formulas, graphs and mathematical tables[END_REF] )

∫ +∞ 0 e -𝜎 𝜓(𝜎)d𝜎 ≃ N ∑ i=0 𝜇 i 𝜓(z i ). (12)
Particularly, this method is more suitable in the case of an integrand with an exponential decay. The coefficients 𝜇 i (resp z i ) are the weights (resp the nodes) of the Gauss-Laguerre quadrature formula. Hence, (12) implies

∫ +∞ 0 𝜓(𝜎) d𝜎 ≃ N ∑ i=0 𝜇 i e z i 𝜓(z i ),
As a consequence, the coefficients in (11) are given by

w i = 𝜇 i e z i ,𝜎 i = z i .
For illustrative purposes, we present in Table 1 the weights and the nodes of this quadrature method when N = 8. We note that this quadrature method was proposed by Yuan and Agrawal 24 in order to approximate the Caputo fractional derivative of order 𝛼 where 0 <𝛼<1.

The Gauss-Jacobi quadrature method

In this paragraph, a second method is proposed, approximating the generalized integral (11) with the Gauss-Jacobi quadrature method. For that purpose, it is necessary to transform the domain of integration from [0, +∞[ to [-1, 1], then to apply the Gauss-Jacobi quadrature to the resulting integral. Hence, we choose the following change of variables:

𝜎 = 1 -z 1 + z then d𝜎 dz = -2 (1 + z) 2 .

It follows that

∫ +∞ 0 𝜓(𝜎)d𝜎 = ∫ 1 -1 2 (1 + z) 2 𝜓 ( 1 -z 1 + z ) dz.
ψ(z)= 2 (1 + z) 2 𝜓 ( 1 -z 1 + z ) ,
we get

∫ +∞ 0 𝜓(𝜎)d𝜎 = ∫ 1 -1 ψ(z)dz ≃ N ∑ i=0 𝜇 i ψ(z i ),
where 𝜇 i (resp z i ) are the weights (resp the nodes) of the standard Gauss-Jacobi quadrature formula over [-1, 1].F o r illustrative purposes, we present in the Table 1 the weights and the nodes of this quadrature formula for N = 8. By identification, we deduce that the quadrature coefficients in (11) are given by

w i = 2 (1 + z i ) 2 𝜇 i ,𝜎 i = 1 -z i 1 + z i .
We note that this strategy was proposed by Diethelm [START_REF] Diethelm | An investigation of some nonclassical methods for the numerical approximation of Caputo-type fractional derivatives[END_REF] for the approximation of Caputo fractional-order derivative.

A FIRST NUMERICAL SCHEME

In this section, we first present the numerical scheme. After that, we present different numerical simulations with both Gauss-Laguerre and Gauss-Jacobi quadrature methods. Finally, we compare these results in order to choose the appropriate quadrature that will be used in the second numerical scheme. In the sequel, we take t ∈[ 0, T], x ∈[ 0, A] and we consider periodic boundary conditions, ie, u| x=0 = u| x=A . Since we run the numerics with an initial data that provides a wave that moves from the left to the right, we expect our computations to be physically relevant until this wave reaches the right boundary.

Presentation of the model

Let Δt > 0b eaf i x e dt i m es t e pa n dh > 0 be a fixed space step. We set for all m = 0, •••, N p , t m = mΔt and for all j = 1, •••, K, x j = jh. We use the following notations:

∀i = 1, •••, N, u m (x)≈u(mΔt, x),𝜓 m i (x)≈𝜓(mΔt, x,𝜎 i ).
For the PDE in (10), we suggest a semi-implicit scheme, where the half derivative and the nonlinear term are explicitly approximated. We use the Crank-Nicolson approximation in time. Thus, the semi-discretized scheme associated to (10) is given by

u m+1 -u m Δt + √ 𝜈 N ∑ i=1 w i ( 2 𝜋 u m -𝜎 2 i 𝜓 m i ) = 𝛼 ( u m+1 + u m 2 ) xx - ( u m+1 + u m 2 ) x -𝛽 ( u m+1 + u m 2 ) xxx -𝛾 ( (u m ) 2 2 ) x . ( 13
)
For the ODE in (10), we use the following implicit-explicit scheme

𝜓 m+1 i -𝜓 m i Δt =-𝜎 2 i 𝜓 m+1 i + 2 𝜋 u m . ( 14
)
Then, we apply the Fourier transform in space to ( 13) and ( 14). We get

ûm+1 (1 + Δt 2 (𝛼𝜉 2 + i𝜉 -i𝜉 3 𝛽)) = ûm (1 - Δt 2 (𝛼𝜉 2 + i𝜉 -i𝜉 3 𝛽)) - √ 𝜈Δt N ∑ i=1 w i ( 2 𝜋 ûm -𝜎 2 i ψm i )-i𝛾Δt𝜉 (u m ) 2 2 ( 15
)
and

ψm+1 i (1 + 𝜎 2 i Δt)= ψm i + 2 𝜋 Δt ûm , ( 16
)
where the initial conditions are approximated as follows:

∀i = 1, •••, N, ψ0 i = 0, û0 ≈ û0 .
Hence, our scheme is a 2-step scheme: Let the approximated solution (u m ,𝜓 m ) be given, we first calculate an approximate solution 𝜓 m+1 of ( 16). Then, we use this solution to approximate the solution u m+1 of ( 15), and so on.

Numerical results and discussion

SimilartoChenetal 3 and Manoubi, [START_REF] Manoubi | Theoretical and numerical analysis of the decay rate of solutions to a water wave model with a nonlocal viscous dispersive term with Riemann-Liouville half derivative[END_REF] we choose an initial datum u 0 , which provides a small amplitude and long-wave Korteweg-de Vries (KdV) soliton for 𝛼 = 𝜈 = 0, 𝛽 = 1, and 𝛾 = 6. We denote by x 0 the middle of the space interval, then the initial datum reads u 0 (x)=0.32 * sech 2 (0.4 * (x -x 0 )).

All numerical simulations are realized with the time step Δt = 0.1, the space step h = 0.1, and on the space interval [0, A]=[ 0, 500]. In a first step, we study the convergence of the numerical schemes ( 15) and ( 16) using the quadrature coefficients of Gauss-Laguerre and Gauss-Jacobi. In a second step, we compare these numerical solutions with those constructed using the scheme (68) developed in Manoubi. [START_REF] Manoubi | Theoretical and numerical analysis of the decay rate of solutions to a water wave model with a nonlocal viscous dispersive term with Riemann-Liouville half derivative[END_REF] In this paragraph, we take the following parameters values 𝛼 = 𝜈 = 1, 𝛾 = 𝛽 = 0, and T = 100, which correspond to the nonlinear case without geometric dispersion. In Figure 1, we present the numerical solution associated to the schemes ( 15) and ( 16) increasing the number of Gauss-Laguerre quadrature nodes; more precisely, we take N = 30, N = 40, N = 50, N = 60, and N = 80. Moreover, we present the L 2 and L ∞ decay rates of the solution with N = 40, N = 50, N = 60, and N = 80. We observe that the numerical solutions converge for N ≥ 60. However, we do not observe the qualitative convergence of the decay rates even for N = 80, which represents an important number of points. Also, we observed this slow convergence using other parameters values of the scheme. This may be due to the poor quality of the algorithm of the Gauss-Laguerre quadrature; this fact has been observed by other authors. [START_REF] Lu | Wave field simulation for heterogeneous porous media with singular memory drag force[END_REF][START_REF] Schmidt | On a critique of a numerical scheme for the calculation of fractionally damped dynamical systems[END_REF] Later, this method was mathematically analyzed by Diethelm. [START_REF] Diethelm | An investigation of some nonclassical methods for the numerical approximation of Caputo-type fractional derivatives[END_REF] Diethelm identified that this unsatisfactory behavior is due to the asymptotic characteristics of the integrand function at the end points of the interval [0, +∞).On the basis of this result, Diethelm proposed an alternative quadrature formula, namely, the Gauss-Jacobi quadrature, that gives rise to a much better overall performance of the scheme and that is more promising for integrands with an algebraic decay.

In Figure 2, we present the approximate solution of the schemes ( 15) and ( 16) with the Gauss-Jacobi quadrature coefficients and the L 2 and L ∞ decay rates of the solution. We observe that the numerical solutions and the decay rates converge for a number of nodes greater than N = 15. We have constructed different numerical simulations using different parameters values, and we get the same observation. In Figure 3, we compare the numerical solutions and decay rates resulting from the schemes ( 15) and ( 16) using the Gauss-Laguerre and Gauss-Jacobi rules with those arising from the numerical scheme (68) developed in Manoubi. [START_REF] Manoubi | Theoretical and numerical analysis of the decay rate of solutions to a water wave model with a nonlocal viscous dispersive term with Riemann-Liouville half derivative[END_REF] We use 80 points for Gauss-Laguerre quadrature and 15 points for Gauss-Jacobi quadrature. We observe that the 3 numerical solutions overlap each other. Moreover, it is sufficient to use only N = 15 points with the Gauss-Jacobi rule to obtain satisfying numerical results for the numerical solution. However, using Gauss-Laguerre quadrature, we need to use N = 50 points or more if we need to get more precise results. Moreover, using Gauss-Jacobi rule, we observe that the decay rates of the solution are smoother than both numerical decay rates resulting from the scheme (68) developed in Manoubi [START_REF] Manoubi | Theoretical and numerical analysis of the decay rate of solutions to a water wave model with a nonlocal viscous dispersive term with Riemann-Liouville half derivative[END_REF] and from Gauss-Laguerre rule. Taking in account these observations, we choose in the sequel the Gauss-Jacobi quadrature to all numerical simulations. Now, we simulate the numerical schemes ( 15) and ( 16) in the nonlinear case. To this end, we take the following parameters values: 𝛼 = 2, 𝜈 = 0.2, 𝛾 = 1, 𝛽 = 0.5, and T = 100.

In Figure 4, we present the numerical solution and L 2 and L ∞ decay rates of the schemes ( 15)and ( 16) using N = 15 points in the Gauss-Jacobi rule. We compare these results with those constructed from the numerical scheme (68) developed in Manoubi. [START_REF] Manoubi | Theoretical and numerical analysis of the decay rate of solutions to a water wave model with a nonlocal viscous dispersive term with Riemann-Liouville half derivative[END_REF] We observe that the numerical solutions overlap each other. Moreover, the decay rates in L 2 and L ∞ norms are very close. Then, we can deduce that both schemes converge. Also, in Table 2, we compare the computation time of solutions in linear and nonlinear cases (in seconds) between the numerical schemes ( 15) and ( 16) for different Gauss quadrature rules and the numerical scheme (68) developed in Manoubi. [START_REF] Manoubi | Theoretical and numerical analysis of the decay rate of solutions to a water wave model with a nonlocal viscous dispersive term with Riemann-Liouville half derivative[END_REF] We observe that the numerical scheme (68) developed in Manoubi 5 is relatively expensive in computation time. This is due to the numerical method used to approximate the half-order derivative.

A SECOND NUMERICAL SCHEME

To improve the precision and the efficiency of the numerical scheme used before, we construct in this subsection a second numerical scheme associated to (8) following the method proposed in Lombard and Mercier. [START_REF] Lombard | Numerical modeling of nonlinear acoustic waves in a tube connected with Helmholtz[END_REF]

Presentation of the model

The PDE-ODE system given by ( 9) can be written as

⎧ ⎪ ⎨ ⎪ ⎩ 𝜕 t u(t, x)+𝜕 x (u + 𝛾 2 u 2 )=- √ 𝜈 ∫ ∞ 0 ( 2 𝜋 u(t, x)-𝜎 2 𝜓(t, x,𝜎) ) d𝜎 + 𝛼u xx (t, x)-𝛽u xxx (t, x), t > 0, x ∈ R 𝜕 t 𝜓(t, x,𝜎)=-𝜎 2 𝜓(t, x,𝜎)+ 2 𝜋 u(t, x), t > 0, x ∈ R,𝜎 ≥ 0, (18)
supplemented with the initial conditions To approximate the integral in (18), we use the Gauss-Jacobi quadrature with N nodes. We denote by w i and 𝜎 i for i = 1, •••, N the weights and abscissas, respectively. Hence, the system ( 18) is written as a first-order system as follows:

∀x ∈ R,∀𝜎 ≥ 0,𝜓 (0, x,𝜎)=0, ∀x ∈ R, u(0, x)=u 0 (x).
⎧ ⎪ ⎪ ⎨ ⎪ ⎪ ⎩ 𝜕 t u(t, x)+ 𝜕 x (u + 𝛾 2 u 2 )=- √ 𝜈 N ∑ i=1 w i ( 2 𝜋 u(t, x)-𝜎 2 i 𝜓 i (t, x)
)

+𝛼u xx (t, x)-𝛽u xxx (t, x), t > 0, x ∈ R, 𝜕 t 𝜓 i (t, x)=-𝜎 2 i 𝜓 i (t, x)+ 2 𝜋 u(t, x), t > 0, x ∈ R, i = 1, •••, N, (19)
endowed with the initial conditions Now, we note by

∀x ∈ R, ∀i = 1, •••, N,𝜓 i (0, x)=0, ∀x ∈ R, u(0, x)=u 0 (x).
U =(u,𝜓 1 , •••,𝜓 N ) T ,
the vector of (N + 1) unknowns, by

 (U)=(u + 𝛾 2 u 2 , 0, •••, 0) T ,
and finally

S = ⎛ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎝ - √ 𝜈 N ∑ i=1 w i √ 𝜈w 1 𝜎 2 1 ••• ••• √ 𝜈w N 𝜎 2 N 2 𝜋 -𝜎 2 1 0 ••• 0 2 𝜋 0 -𝜎 2 2 ⋱ 0 ⋮⋮ ⋱ 2 𝜋 0 ••• ••• -𝜎 2 N ⎞ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎠ .
Thus, the problem ( 19) is written in the following form:

𝜕 t U + 𝜕 x  (U)=S(U)+G 1 𝜕 2 x U + G 2 𝜕 3 x U, (20)
where G 1 and G 2 are diagonal matrices of order N + 1. In the sequel, we introduce the so-called splitting scheme.

The splitting method

Let Δt > 0, for all m = 0, •••, N p ,wenotebyt m = mΔt and

U m (x)≈U(mΔt, x).
From Equation 20, we consider the propagation equation

𝜕 t U + 𝜕 x  (U)=G 1 𝜕 2 x U + G 2 𝜕 3 x U, (21)
and the diffusive equation

𝜕 t U = S(U). (22)
We note by H a (respectively, H b ) the discrete operator of the solution of (21) (respectively, the solution of Equation 22). Then, a Strang splitting method [START_REF] Holden | Operator splitting for the KDV equation[END_REF][START_REF] Leveque | Numerical Methods for Conservation Laws[END_REF] of order 2 between t n and t n+1 is used to solve respectively ( 21) and ( 22) as follows:

U (1) = H b ( Δt 2 )U m ,
U [START_REF] Dutykh | Viscous potential free surface flows in a fluid layer of finite depth[END_REF] = H a (Δt)U [START_REF] Liu | Viscous effects on transient long wave propagation[END_REF] ,

U m+1 = H b ( Δ 2
)U [START_REF] Dutykh | Viscous potential free surface flows in a fluid layer of finite depth[END_REF] .

(

) 23
Here, the constructed operators H a and H b are stable and of order 2. Then, the scheme (23) provides an approximation of order 2 in time to the problem (21). In the sequel, we present the discretization of ( 21) and (22).

Propagation equation 21

Here, u is a solution of the KdV equation. The discretization in time of this equation is provided by a semi-discrete in time scheme, using the Crank-Nicolson scheme for the linear part and Adams-Bashforth scheme (see Kalisch and Bona 32 )for the nonlinear part. The space approximation of the solutions are performed by standard Fourier methods. First, we note that the approximate solution û1 is provided by a fixed-point method that verifies the semi-discrete scheme (of order 2):

û1 -û0 Δt = û1 + û0 2 (-𝛼𝜉 2 -i𝜉 + i𝛽𝜉 3 )- i𝛾𝜉 2 (u 1 ) 2 + (u 0 ) 2 4 . ( 24
)
Then, for m = 1, •••, N p , we consider the following semi-discrete scheme:

ûm+1 -ûm Δt = ûm+1 + ûm 2 (-𝛼𝜉 2 -i𝜉 + i𝛽𝜉 3 )- i𝛾𝜉 2 3 (u m ) 2 -( u m-1 ) 2 4 . ( 25
)
This scheme has local truncation error of order (Δt) 2 , and a second-order convergence is observed (for more details, see Kalisch and Bona 32 ).

Diffusion equation 22

We can solve mathematically (22) as follows:

H b ( Δt 2 )U = e S Δt 2 U. ( 26
)
In addition, for all N > 0, the exponential of the matrix S Δt 2 is calculated numerically using the "scaling and squaring" method with a (6, 6) Padé approximation [START_REF] Moler | Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later[END_REF] that corresponds, with Matlab, to the function "expm(S Δt 2 )." We recall that the (p, q) Padé approximation of e A is given by

R 𝑝q (A)=(D 𝑝q (A)) -1 N 𝑝q (A),
where

N 𝑝q (A)= 𝑝 ∑ 𝑗=0 (𝑝 + q -𝑗)!𝑝! (𝑝 + q)!𝑗!(𝑝 -𝑗)! A 𝑗 , D 𝑝q (A)= 𝑝 ∑ 𝑗=0 (𝑝 + q -𝑗)!q! (𝑝 + q)!𝑗!(q -𝑗)! (-A) 𝑗 .
Hence, the numerical scheme ( 23) is a 3-step scheme. Since the scheme (25) as well as the Strang splitting method are of order 2 in time, it follows that ( 23) is of order 2 in time.

Numerical results and discussion

In this paragraph, we present some numerical results from the numerical scheme (23) with similar parameters and initial datum as in the case of the first schemes 15 and (16). We analyze the convergence in time, and we study the effects of the number of quadrature points on the numerical solution. To this end, we choose 𝛼 = 2, 𝜈 = 1, 𝛾 = 3, 𝛽 = 1, T = 100 and h = 0.05.

Convergence in time

To study the convergence in time of the numerical solution of (23), we take N = 40 Gauss-Jacobi quadrature points. In Figure 5, we present different numerical solutions for Δt = 0.05, Δt = 0.1, Δt = 0.2, Δt = 0.3, and Δt = 0.4, as well as the decay rates in norms L 2 and L ∞ in this case. We observe that the difference between different numerical solutions and decay rates are negligible.

To evaluate the error due to the time discretization of the scheme ( 23 we illustrate the L 2 -error E m in terms of the time step Δt in the log-log scale. We observe that the measured values are close to a straight line with slope 2. This numerically confirms that the numerical scheme ( 23) is of order 2 in time.

Influence of the number of Gauss-Jacobi quadrature points

To analyze the influence of the number of the grid points N on the numerical solution of ( 23), we denote by u m Ref the reference solution, which is the numerical solution for Δt = 0.05 and N = 40. We denote by u m N the numerical solution for different values of

N.W edefinetheerrorE m S𝑝 (N)=|| u m Ref -u m N || 2 in norm L 2 .
In Figure 7, we present E m S𝑝 in terms of N. We observe that as expected the error decreases as the number of N increases. Also, the measured values are close to a straight line with slope -5, which let us conclude that E m S𝑝 ≈ C(

DISPERSION ANALYSIS

We recall the definition of the half-order derivative of Riemann-Liouville

D 1∕2 u(t)= 1 √ 𝜋 𝜕 t ∫ +∞ 0 u(s) √ t -s ds. ( 27
)
In this section, we present the dispersion analysis of the model (2) in the linear case (𝛾 = 0) using 2 methods. First, let us define the Laplace and Fourier transforms in time and space, respectively. The first method consists in applying Laplace-Fourier transform to (2) in order to derive the dispersion relation between 𝜏 and k. Then, we define the symbol of the half-order derivative ( 27)

𝜒(𝜏)= √ i𝜏 = √ 𝜏 ( 1 √ 2
+ i 1 √ 2
) .

We get from (2) the following identity:

(i𝜏 + √ 𝜈 √ i𝜏) û +(ik -i𝛽k 3 ) û =-𝛼k 2 û.
We set 𝛼(𝜏,k)=-i𝛽k 3 + 𝛼k 2 + ik + i𝜏.

Hence, the dispersion relation is given by 𝛼(𝜏,k)+ √ 𝜈𝜒(𝜏)=0. (28) We note that we get the same dispersion relation as in Chen et al [START_REF] Chen | Decay of solutions to a water wave model with a nonlocal viscous dispersive term[END_REF] where the authors studied the linear dispersion of (2) when the nonlocal diffusive-dispersive term is described by the Caputo half-order derivative with an initial condition u(0)=0.

The second method is devoted to the linear dispersion relation of the system (20) (we take 𝛾 = 0). To this end, we apply Fourier transforms in time and space to (20), we get a set of (N + 1) equations whose determinant must be vanishing. Then, the matrix form of the resulting system is given by

A Û = 0, where A = ⎛ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎝ 𝛼(𝜏,k)+ 2 𝜋 √ 𝜈 N ∑ i=1 w i - √ 𝜈w 1 𝜎 2 1 ••• ••• - √ 𝜈w N 𝜎 2 N 2 𝜋 -𝜎 2 1 -i𝜏 0 ••• 0 2 𝜋 0 -𝜎 2 2 -i𝜏 ⋱ 0 ⋮⋮ ⋱ 2 𝜋 0 ••• ••• -𝜎 2 N -i𝜏 ⎞ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎠ .
Thus, the dispersion relation between 𝜏 and k is given by

𝛼(𝜏,k)+ 2 𝜋 √ 𝜈 ( N ∑ 𝑗=1 w 𝑗 - N ∑ 𝑗=1 w 𝑗 𝜎 2 𝑗 𝜎 2 𝑗 + i𝜏 ) = 0. ( 29
)
We denote by

χ(𝜏)= 2 𝜋 ( N ∑ 𝑗=1 w 𝑗 - N ∑ 𝑗=1 w 𝑗 𝜎 2 𝑗 𝜎 2 𝑗 + i𝜏 ) ,
then the dispersion relation ( 29) can be written in a similar way as ( 28)

𝛼(𝜏,k)+ √ 𝜈χ(𝜏)=0. ( 30
)
It is worth to note that when we consider the main Equation 2 without nonlocal diffusive-dispersive term, ie, when 𝜈 = 0, the dispersion relations ( 28) and (30) reduced to the KdV-dispersion relation, namely, 𝛼(𝜏,k)=0.

FIGURE 1

 1 FIGURE 1 Numerical solution and decay rates using Gauss-Laguerre rule for different N when 𝛼 = 𝜈 = 1, 𝛽 = 𝛾 = 0 [Colour figure can be viewed at wileyonlinelibrary.com]

FIGURE 2

 2 FIGURE 2 Numerical solution and decay rates using the Gauss-Jacobi for different N with 𝛼 = 𝜈 = 1, 𝛽 = 𝛾 = 0 [Colour figure can be viewed at wileyonlinelibrary.com]

FIGURE 3

 3 Scheme (68) from[19] Gauss-Laguerre N=80 Gauss-Jacobi N=15

FIGURE 4

 4 Scheme (68) from[20] Gauss Jacobi N=15

4 FIGURE 5

 45 FIGURE 5 Comparison of the splitting scheme solutions and decay rates in L ∞ and L 2 norms where 𝛼 = 2, 𝜈 = 1, 𝛾 = 3, 𝛽 = 1 for different time steps Δt [Colour figure can be viewed at wileyonlinelibrary.com]

e

  t (𝑓 )(𝜏)= f (𝜏)= ∫ +∞ 0 𝑓 (t)e -𝜏t dt,  x (𝑓 )(k)= f (k)= ∫ +∞ -∞ 𝑓 (x)e -ikx dx,where 𝜏 is the angular frequency and k is the wavenumber.The Laplace-Fourier transform of a function f readsf (𝜏,k)= ∫ -ixk 𝑓 (t, x)dx ) dt.

  ds represents the Caputo half derivative in time. Here, u is the horizontal velocity of the fluid, -𝛼u xx is the usual diffusion, 𝛽u xxx is the geometric dispersion, and[START_REF] Liu | Viscous effects on transient long wave propagation[END_REF]

	1 √ 𝜋	∫	t 0	u t (s) √ t-s
				√ 𝜋	∫	t 0	u t (s) √ t-s	ds stands for the nonlocal diffusive-dispersive term.

TABLE 1

 1 Weights and abscissae of the Gauss-Laguerre and Gauss-Jacobi quadrature formulas for N = 8

	Gauss-Laguerre		Gauss-Jacobi
	w i	𝜎 i	w i	z i	𝜎 i
	0.4377	0.1703	128.3897	-0.9603	49.3650
	1.0339	0.9037	10.7575	-0.7967	8.8361
	1.6697	2.2511	2.7870	-0.5255	3.2153
	2.3769	4.2667	1.0879	-0.1834	1.4493
	3.2085	7.0459	0.5179	0.1834	0.6899
	4.2686	10.7585	0.2696	0.5255	0.3110
	5.8181	15.7407	0.1378	0.7967	0.1132
	8.9062	22.8631	0.0527	0.9603	0.0203
	Moreover, letting

TABLE 2

 2 Computation time of the numerical solution

		Linear case Nonlinear case
	Scheme (68) in Manoubi 5 104.26	125.81
	Gauss-Laguerre (N = 80)	61.195	61.45
	Gauss-Jacobi (N = 15)	40.44	40.24

  Error of the time discretization [Colour figure can be viewed at wileyonlinelibrary.com] Error of the splitting scheme in terms of N [Colour figure can be viewed at wileyonlinelibrary.com]

						|| Err( t) || 2
		10 -4				line of slope 2
	Time Error	10 -5
		10 -6
		0.02	0.04	0.06 0.08 0.1	0.2	0.4
				t
	15 FIGURE 6 N L 2 Error 10 -7 10 -6 10 -5 10 -4	20	25 || Err(N) || 2 line with slope -5	30
	FIGURE 7

1 N )

[START_REF] Manoubi | Theoretical and numerical analysis of the decay rate of solutions to a water wave model with a nonlocal viscous dispersive term with Riemann-Liouville half derivative[END_REF]

where C is a constant.

CONCLUSIONIn this article, we propose different numerical schemes to approximate the solution to an asymptotical water wave model where the nonlocal viscosity is described by the Riemann-Liouville half-order derivative. We compare different numerical results of different numerical schemes where the half derivative is described by a diffusive representation. Moreover, we compare our numerical results with those given in Minoubi.[START_REF] Manoubi | Theoretical and numerical analysis of the decay rate of solutions to a water wave model with a nonlocal viscous dispersive term with Riemann-Liouville half derivative[END_REF]

ACKNOWLEDGEMENTS

This work was supported by PHC Utique ASEO and exchange program CNRS/DGRS "modèles non locaux en dynamique des fluides." A part of this article was performed during the stay of the second author in the LAMFA with the support of the exchange program "SSHN2015." ORCID I. Manoubi http://orcid.org/0000-0003-2243-1330


Pascal Floquet

Xavier Joulia

Alain Vacher

Martin Gainville

Michel Pons

Pressure-Driven Steady-State Simulation of Oilfield Infrastructure

Numerical strategy in Pressure-driven simulation Pressure-driven problem is a particular case of design problem Variables associated to pressure constraints are the well flowrates other variables Two types of problem are treated : Flowrates/Pressure problems Well flowrates and riser top pressure are known Action variables are chokes openings or well pressures Pressures/Pressure problems Well pressures and riser top pressure are known Action variables are well flowrates

Purpose

  IFP research collaborative project A platform for integrated multiphase flow simulations From reservoir to process facilities For flow assurance application Based on an open software architecture with INDISS as reference simulator IFP research collaborative project A platform for integrated multiphase flow simulations From reservoir to process facilities For flow assurance application Based on an open software architecture with INDISS as reference simulator IFP research collaborative project A platform for integrated multiphase flow simulations From reservoir to process facilities For flow assurance application Based on an open software architecture with INDISS as reference simulator IFP research collaborative project A platform for integrated multiphase flow simulations From reservoir to process facilities For flow assurance application Based on an open software architecture with INDISS as reference simulator IFP research collaborative project A platform for integrated multiphase flow simulations From reservoir to process facilities For flow assurance application Based on an open software architecture with INDISS as reference simulator IFP research collaborative project A platform for integrated multiphase flow simulations From reservoir to process facilities For flow assurance application Based on an open software architecture with INDISS as reference simulator IFP research collaborative project A platform for integrated multiphase flow simulations From reservoir to process facilities For flow assurance application Based on an open software architecture with INDISS as reference simulator Extend simultaneous modular strategy for solving steady-state pressure-driven simulation and design problems in Oil and Gas production networks Interoperate via CAPE-OPEN Extend INDISS simulator capabilities by integration of ProSim CO-SPEC module Extend INDISS simulator capabilities by integration of IFP TOTAL pipeline multiphase flow module Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 study Extend simultaneous modular strategy for solving steady-state pressure-driven simulation and design problems in Oil and Gas production networks Interoperate via CAPE-OPEN Extend INDISS simulator capabilities by integration of ProSim CO-SPEC module Extend INDISS simulator capabilities by integration of IFP TOTAL pipeline multiphase flow module Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 for solving steady-state pressure-driven simulation and design problems in Oil and Gas production networks Interoperate via CAPE-OPEN Extend INDISS simulator capabilities by integration of ProSim CO-SPEC module Extend INDISS simulator capabilities by integration of IFP TOTAL pipeline multiphase flow module Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 for solving steady-state pressure-driven simulation and design problems in Oil and Gas production networks Interoperate via CAPE-OPEN Extend INDISS simulator capabilities by integration of ProSim CO-SPEC module Extend INDISS simulator capabilities by integration of IFP TOTAL pipeline multiphase flow module Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 for solving steady-state pressure-driven simulation and design problems in Oil and Gas production networks Interoperate via CAPE-OPEN Extend INDISS simulator capabilities by integration of ProSim CO-SPEC module Extend INDISS simulator capabilities by integration of IFP TOTAL pipeline multiphase flow module Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 for solving steady-state pressure-driven simulation and design problems in Oil and Gas production networks Interoperate via CAPE-OPEN Extend INDISS simulator capabilities by integration of ProSim CO-SPEC module Extend INDISS simulator capabilities by integration of IFP TOTAL pipeline multiphase flow module Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 for solving steady-state pressure-driven simulation and design problems in Oil and Gas production networks Interoperate via CAPE-OPEN Extend INDISS simulator capabilities by integration of ProSim CO-SPEC module Extend INDISS simulator capabilities by integration of IFP TOTAL pipeline multiphase flow module Base Case Base case description Two subsea production clusters, Two and three subsea wells controlled by choking wellhead valves, Two flowlines connected to a riser and a basic surface process Part II Pressure-Driven Steady-State Simulation Pressure-driven vs classical simulation Sequential Modular Simulation X 0 :Temperature T, Pressure P, Composition z and total flowrate Q d : operating and design parameters of the modules are the standard input of a pure simulation case Pressure-driven Simulation Characterize Oil and Gas upstream operations For example, pipes connected to the same manifold must operate at the same pressure n w Pressure Equality Constraints vs classical simulation Sequential Modular Simulation X 0 :Temperature T, Pressure P, Composition z and total flowrate Q d : operating and design parameters of the modules are the standard input of a pure simulation case Pressure-driven Simulation Characterize Oil and Gas upstream operations For example, pipes connected to the same manifold must operate at the same pressure n w Pressure Equality Constraints We are interested to examine the ability of convergence of Sequential Modular Simulator in Pressure Driven Problem Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 pressure at riser top pressures, flowrates and temperature are fixed (base case) Action variables : 5 pressure drops of the 5 chokes Quasi-Newton Strategy used Results obtained in 4 iterations and 11 flowsheet simulations, for 15 bar Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 pressure at riser top pressures, flowrates and temperature are fixed (base case) Action variables : 5 pressure drops of the 5 chokes Quasi-Newton Strategy used Results obtained in 4 iterations and 11 flowsheet simulations, for 15 bar Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 pressure at riser top Pressures drops of the 5 chokes, well pressures and temperature are fixed Action variables : 5 well flowrates Results obtained in 5 iterations and 12 flowsheet simulations, for 15 bar ...but it depends on initialization ! pressure at riser top Pressures drops of the 5 chokes, well pressures and temperature are fixed Action variables : 5 well flowrates Results obtained in 5 iterations and 12 flowsheet simulations, for 15 bar ...but it depends on initialization ! Specification on the pressure at riser top Action variables : pressure drops of 4 chokes and flowrate of gas-lift Tear stream : compressor output (15 + 2 (T, P) iterative variables) Less eruptive choke completly open Results obtained in 5 iterations and 12 flowsheet simulations, for 30 bar Recycle rate versus Riser top pressure specification

Pascal Floquet 1 ,

Pressure-Driven Steady-State Simulation of Oilfield Infrastructure Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 , Martin Gainville 3 , Michel Pons 4 Without recycle With recycle Flowrates/Pressure Problem Specification on the pressure at riser top Action variables : pressure drops of 4 chokes and flowrate of gas-lift Tear stream : compressor output (15 + 2 (T, P) iterative variables) Less eruptive choke completly open Results obtained in 5 iterations and 12 flowsheet simulations, for 30 bar Recycle rate versus Riser top pressure specification Pressure-Driven Steady-State Simulation of Oilfield Infrastructure Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 , Martin Gainville 3 , Michel Pons 4 CO Results Part IV CAPE-OPEN Integration Pressure-Driven Steady-State Simulation of Oilfield Infrastructure Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 , Martin Gainville 3 , Michel Pons 4 CO Results Interoperability ProSimPlus is reference simulator chosen to adjust the approach IFP TOTAL pipeline multiphase flow modules are specialized upstream Oil an Gas modules ProSim CO-SPEC (CAPE-OPEN Unit Operation 1.0) and IFP TOTAL pipeline multiphase flow modules integrated in INDISS-TINA environment Pressure-Driven Steady-State Simulation of Oilfield Infrastructure Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 , Martin Gainville 3 , Michel Pons 4 CO Results Interoperability ProSimPlus is reference simulator chosen to adjust the approach IFP TOTAL pipeline multiphase flow modules are specialized upstream Oil an Gas modules ProSim CO-SPEC (CAPE-OPEN Unit Operation 1.0) and IFP TOTAL pipeline multiphase flow modules integrated in INDISS-TINA environment Pressure-Driven Steady-State Simulation of Oilfield Infrastructure Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 , Martin Gainville 3 , Michel Pons 4 CO Results Interoperability ProSimPlus is reference simulator chosen to adjust the approach IFP TOTAL pipeline multiphase flow modules are specialized upstream Oil an Gas modules ProSim CO-SPEC (CAPE-OPEN Unit Operation 1.0) and IFP TOTAL pipeline multiphase flow modules integrated in INDISS-TINA environment Pressure-Driven Steady-State Simulation of Oilfield Infrastructure Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 , Martin Gainville 3 , Michel Pons 4 CO Results Interoperability ProSimPlus is reference simulator chosen to adjust the approach IFP TOTAL pipeline multiphase flow modules are specialized upstream Oil an Gas modules ProSim CO-SPEC (CAPE-OPEN Unit Operation 1.0) and IFP TOTAL pipeline multiphase flow modules integrated in INDISS-TINA environment Pressure-Driven Steady-State Simulation of Oilfield Infrastructure Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 , Martin Gainville 3 , Michel Pons 4 Part V Conclusion and future work Pressure-Driven Steady-State Simulation of Oilfield Infrastructure Pascal Floquet 1 , Xavier Joulia 1 , Alain Vacher 2 , Martin Gainville 3 , Michel Pons 4

To conclude...

CAPE tools such as ProSimPlus are able to solve pressure-driven simulation CAPE-OPEN standards are the best way to plug-and-play CAPE software components Future Works Multi period optimization in an Oil and Gas context