Papers by Fayssal Benkhaldoun
Lecture Notes in Computer Science, 2023
Lecture Notes in Computer Science, 2023

Advances in Computational Mathematics, Jun 20, 2019
This paper deals with the numerical solution of an ionization wave propagation in air, described ... more This paper deals with the numerical solution of an ionization wave propagation in air, described by a coupled set of convection-diffusion-reaction equations and a Poisson equation. The standard three-species and more complex eleven-species models with simple chemistry are formulated. The PDEs are solved by a finite volume method that is theoretically second order in space and time on an unstructured adaptive grid. The upwind scheme and the diamond scheme are used for the discretization of the convective and diffusive fluxes, respectively. The Poisson equation is also discretized by the diamond scheme. The results of both models are compared in details for a test case. The influence of physically pertinent boundary conditions at electrodes is also presented. Finally, we deal with numerical accuracy study of implicit scheme in two variants for simplified standard model. It allows us in the future to compute simultaneously and efficiently a process consisting of short time discharge propagation and long-term after-discharge phase or repetitively pulsed discharge.
Journal of Elliptic and Parabolic Equations

Mathematical Methods in the Applied Sciences
We propose a class of finite volume algorithms that are both simple and efficient for solving num... more We propose a class of finite volume algorithms that are both simple and efficient for solving numerically the shallow water equations with varying densities; shallow water flows in single and two layers are considered. In these flow regimes, variable horizontal or vertical density is taken into account. The shallow water equations for the hydraulic variables are coupled with a suspended sediment transport equation for the concentration variable to construct the model. To approximate the numerical solution of the models under consideration, a generalized Rusanov method is proposed; this method avoids solving Riemann problems during the time integration process, and it is simple and accurate. The proposed method is divided into two stages: predictor and corrector; the presented finite volume approach is well balanced, conservative, and simple. Non‐oscillatory and suitable for shallow water equations when Riemann problems are challenging to solve. Several test problems for single‐layer...
Mathematical Methods in the Applied Sciences
Large-Scale Scientific Computing, 2022

Finite Volumes for Complex Applications IX - Methods, Theoretical Aspects, Examples, 2020
In this note, we establish a finite volume scheme for a model of a second order hyperbolic equati... more In this note, we establish a finite volume scheme for a model of a second order hyperbolic equation with a time delay in any space dimension. This model is considered in [10, 11] where some exponential stability estimates and oscillatory behaviour are proved. The scheme we shall present uses, as space discretization, the general class of nonconforming finite volume meshes of [5]. In addition to the proof of the existence and uniqueness of the discrete solution, we develop a new discrete a priori estimate. Thanks to this a priori estimate, we prove error estimates in discrete seminorms of \(L^\infty (H^1_0)\), \(L^\infty (L^2)\), and \(W^{1,\infty }(L^2)\). This work can be viewed as extension to the previous ones [2, 4] which dealt with the analysis of finite volume methods for respectively semilinear parabolic equations with a time delay and the wave equation.
Large-Scale Scientific Computing, 2020
SUSHI (Scheme Using Stabilization and Hybrid Interfaces) is a finite volume method has been devel... more SUSHI (Scheme Using Stabilization and Hybrid Interfaces) is a finite volume method has been developed at the first time to approximate heterogeneous and anisotropic diffusion problems. It has been applied later to approximate several types of partial differential equations. The main feature of SUSHI is that the control volumes can only be assumed to be polyhedral. Further, a consistent and stable Discrete Gradient is developed.

The European Physical Journal Plus, 2016
Abstract.In this work, we introduce a finite volume method for numerical simulation of shallow wa... more Abstract.In this work, we introduce a finite volume method for numerical simulation of shallow water equations with source terms in one and two space dimensions, and one-pressure model of two-phase flows in one space dimension. The proposed method is composed of two steps. The first, called predictor step, depends on a local parameter allowing to control the numerical diffusion. A strategy based on limiters theory enables to control this parameter. The second step recovers the conservation equation. The scheme can thus be turned to order 1 in the regions where the flow has a strong variation, and order 2 in the regions where the flow is regular. The numerical scheme is applied to several test cases in one and two space dimensions. This scheme demonstrates its well-balanced property, and that it is an efficient and accurate approach for solving shallow water equations with and without source terms, and water faucet problem.

Http Www Theses Fr, 2005
§¢ ¥¤ ) ¨£ ¢¡ ¡ ¤£ ¦¥ §¦ ¾ s ta @ yd ) » ®º ¾ s ta gc Au gi 6} ~¿ ¨ § ' ¥( )0 32 H U 2 © ¥ dU U (... more §¢ ¥¤ ) ¨£ ¢¡ ¡ ¤£ ¦¥ §¦ ¾ s ta @ yd ) » ®º ¾ s ta gc Au gi 6} ~¿ ¨ § ' ¥( )0 32 H U 2 © ¥ dU U ( @ m ( R ´U P 9 8U 0 64 G IB F2 )S 65 8U ¿ e8 &i 6Ì Ê º d |i 6a gd |h 6¿ ³d v ¥#" a gd |h Á @¹ Ô ¬¾ ¿ & C d º Ä ®Ä @Ï G b » ®º Ä ®i 6p eÉ zÊ 8v s tw ea gc Ap ³¿ gÀ º Ä ®a ¥d Tu gÊ eq 1¿ ¥v p ew ¹ º Ï Ca gp & 6v ¥s ti ¥¿ 82 dU 9 8G I2 8( 82 ( U ¥G IB F5 ´U G ! Ñ9 5 ( ' ¥( )2 HG I !G #" %$ ¥0 G 9 B tU ( '& ' ( ´U v2 © ¥ 82 ( )2 G I ¿ Á a g `º ¹ i ~h Ê º q tp ¾ ÑÀès tq ri 6w ¹ i ~Ì Ê º v ¥p ew AE p eu º 7 T7 ) ¢ % b ¤ 3¿ h y 0 º Î i 6d |d Ti 6d Tv A¿ 21 U T' (3 ' (4 ( U 0 2 8( )0 5 2 HG I 2 ± G 6' I ( ' ¥B Ð' S 65 2 " U 0 54 G I dU 9 ( Q 76 8 ´U ( P 8G 9$ I 76 Ó9 0 |' ( 12 H U 76 ¿ w ea gd &" w ed |i 6Ì Ê h ' f a g} h a gp ' Is ta gp ew Aa gp i ~w º ¿ hÒ ªs rc ) ®i ~d f¾ @Ì 6v ¥w ei ~ Ðq rÌ D c Af es tq r} TÊ ei ~d T} 6¿ ¬ TS 65 e' ¥( )2 HG I A A3 2 ( QP P I' ¥0 32 )' 4 6B tU $ U 9 ( P ´5 2 © 5 " A ( )0 35 h ( )5 0 U Q$ R ´U P 8U ¿ Á a g èc Ah º "¹ i ~h TÊ ea "w A} ¾ èès º ¹ i 6Ì Ê º AE p eu gd Tu º 7 6 T6 ¢ % y d ¤ ¿ "p ea º 3" "¿ b S) g º Ä i 6p 8 » ®º Ä i 6p eÉ zÊ 8v ¥s rw ea gc ep d¿ D ' ¥B T m 2 t ª' ¥ D$ A I' ¥B F2 U$ ' ¥( 12 HG I G #" ' Ð dU 3 V 82 ( U C ¥G IB F5 ´U Ö P 8U ´U " G I0 Ö dG I AP 8G ! ´G Q© U dU G I5 m ( U ¿ » 5 Á ¾ 7F A¿ Î D ) ¿ eÀ º eÎ i ~d Tf eq rp ³¿ eË º Ò Cd QW ¥p ei 6d ¢ 1AE w e} ¤ 3¿ D d Ta Ì 6i ~i 6w eq rp eu g} va %© h Ê Ai 8 @Ê eq td |w 1 p zh Ti 6d |p 8v h q ra gp 8v ¥s C) È è} Ta gq rc e Âa ¥p » q tp eq Dh i 5 a gs tc A Ði ~} © Qa gd Á a g ès ri Ç E¾ èès rq tÌ v Ih q ta ¥p ³¿ 8 ¥ ¥ "¿ AÀ`º ya S) " ya º Ä 05 b ¾ º Ä ®i 6d | c ew ei v p ew ¹ º AE º 5 v 6 zc ei ¥¿ XB 9 D3 2 D$ Y ´U ( QP 8G 9$ I `" G I0 @P 9 8U 0 ~4 G IB F2 H Ñ G I A U 0 a4 I' ¥( )2 HG I B D' (3 y E3 2 ( P v G I5 "0 U f( U 0 ¿ Á a g èc Ah º cb » s rc eq tw e} e ($ ¢ % y b ¤ ¿ p ea º d "¿ e b d) @ º Á Ë 8 g d ) ¡ Aº Á a gd h i 6} ~¿ z¾ º Ë Ñi ~f ec e} T} |Ì Ê ei ¥¿ Av ¥p ew `1 º 8 a ¥c e Ðq 1¿ $ U A 2 ( U @9 8U 0 ( )5 0 64 ' ¥( )2 HG I e ´U ( QP 8G 9$ ( G ( ) ¥U 0 " G I0 ( P 8U P A' ¥B B tG !3 i4 p3 ' ¥( U 0 U TS 65 A' ¥( 12 HG I A 3 2 ( P " 30 32 H ( )2 HG I ' ¥ D$ ( G 9 8G !4 © ¥0 |' 9 8P ¿ ¾ Å v ¥q ts tv ¥f es ti q tp h Ê ei Ì 6a gp e} |i 6d Å ¥v Ih q ta ¥p s tv ) èd Ti ~èd Tq rp zh } Ti ~d |Å ¥i 6d Ê h Th T`e ' $' ) ) ) º ´v Ih Ê º p h p zc º p ea ' Ì ~a gp e} Ti ~d |Å v h Tq ta gp ' ¿ ¬ % º Á Ô ') A b ¾ º bÁ Ê eq rp ep 8v È zÈ v A¿ b¾ º " A º Ô ¬i ~d Ta gc "Ç b¿ bv ¥p ew º ) "i 6u ¥c eq tp ³¿ q 3 ®U B B T4 4 ' ¥B D' ¥ d U Q$ 85 r ´U 0 32 H T' ¥B P 8U ´U s" G I0 ' 9 9 0 G ut g2 H ' ¥( )2 HG I !G #" ( QP 8U v P A' ¥B B tG v4 p3 ' ¥( U 0 ´U TS 65 A' ¥( 12 HG I A w3 2 ( QP ( G 9 8G Q© ¥0 |' (4 9 8P Yx ( P 8U 0 U G I ' ¥ d U 9 8P 8U dG ! ´U dG I ¿ E1 p h i ~d Tp 8v Ih q ta ¥p 8v ¥s ¡ a gc ed Tp ev ¥s Aa gp » q tp eq Dh i Å ga gs rc e i ¢ H g b ¤ 3¿ ¢ )i ~s ti 6Ì h d |a gp eq tÌ ¤ º g )4 6B tU ©0 |' 9 9 8G 9 0 3( w$ U V$ U A 2 ( 7 ¿ D Ê º Ë º h TÊ ei 6} |q t} ~¿ C fp eq rÅ ¥i 6d T} |q rh T{ w ei 8v ¥d |q t} 1 Ì 6i 6p h Td Ti w ¢ Öd T} v È ¥¿ & C º Á a ¥d " ¡ Á a gd h i 6} ~¿ ( P 8U G I A ( )0 35 h ( )2 HG I G #" 5 m9 D3 2 D$ P 8U ´U C" G I0 dG I 4 TU TS 65 "2 B F2 )4 60 32 5 ( )0 T' ¥ A 2 HU 8( 3 ®G !4 t9 8P A' I U & G !3 y ¿ Á a g èc Ah Ti 6d b » s tc eq rw e} $ )7 ¢ H g ¤ 3¿ d m y d) %d º AE Ä Ä ®À b h1 º AE s r ´v Ê eq H¿ » ®º Ä ®i ~p eÉ zÊ 8v ¥s tw Aa gc ep ³¿ ¥À º Ä ®a gd |u gÊ eq H¿ Iv ¥p Aw ) º À v ¥u ¥Ê 8v È g¿ ( 85 ´U 0 32 H T' ¥B z ( )5 $ ( G #" ª5 A 9 8P 8U 0 2 H T' ¥B b2 © ¥ 82 ( 12 HG I G #" ' $ ¥0 G 9 B tU ( ¿ 8v ¥Ì ~Ì 6i ~Àh { ba ¥c ed @èc ef es rq tÌ 6v h q ra gp ³¿ A ¥ g b eº AE Ä 5 + ) 1 º AE s r ´v ¥Ê Aq H¿ » ®º Ä i 6p eÉ zÊ 8v ¥s rw ea gc ep d¿ À º 5 Öq rs t} | Ði ~q ti 6d ~¿ º Ï ªs ra ¥h Ê d¿ ¾ º D yv Ih } TÌ Ê zc es rs H¿ yv ¥p ew Ë º yÎ p ei ~s H¿ ¥ 2 82 ( U A ¥G IB F5 ´U Ð 2 H v5 B D' ¥( )2 HG I G #" ' $ ¥0 G 9 B tU ( ¨& ' ( ´U ´2 © ¥ 82 ( 12 HG I §G I 5 4 ( )0 5 8 ( )5 "0 U Q$ ´U P 8U ¿ ¡ Aº Á a ¥ Ðèc "h º ¾ èès º e¹ v h TÊ º 7 43 2$ ¢ % y ¤ ¿ "p ea º ¥¿ ³ %d @ d) " ¥ "¿ ¾ ÑÀès tq ri 6w Ev ¥p Aw EÌ 6a g èc Ah v h q ra gp 8v ¥s 8h a ¥èq tÌ ~} q tp èv ¥d |h Tq Ãv ¥s hw eq V2 di ~d Ti 6p h Tq Ãv ¥s i 6 zc 8v h Tq ta gp A} ¢ HÏ Cd v " Ðv ¥w ea e¿ ³ C @ ¤ º Ï Î ) " F ) 8 º Ï v ¥s rs ta gc 3 ~h ¿ ¡ " ¹ º Î { 6d Tv ¥d Tw d¿ yv ¥p Aw º ) "i 6u ¥c eq tp ³¿ i1 G ! ´U ' 9 9 0 G ut g2 H ' ¥( U G 9$ ¥5 dG ! P 8U ´U Ó( G W G ! Ñ9 5 "( U Ó P A' ¥B B tG !3 i4 p3 ' ¥( U 0 ´U TS 65 A' ¥( 12 HG I A %3 2 ( P ( G 9 8G Q© ¥0 T' 9 8P ¿ Á a g èc Ah º b » s tc eq rw e} $ Te ¢ H g F ¤ 3¿ ep ea º b ¿ b y S) A GF º Ï ªÒ ªÔ Á A ¡ " ¹ º Ï ªÊ eq rw 8v ¥u gs rq Ãv A¿ I¾ º Ò Cc e f 8v ¥d Ta A¿ zv ¥p ew ´Ï º Ô ¬i Á a " h¿ ¨ E( P 8U Ñ 85 ´U 0 32 H T' ¥B A G IB F5 "( )2 HG I
Presented paper describes a numerical study of discharge plasma motion. This non-stationary pheno... more Presented paper describes a numerical study of discharge plasma motion. This non-stationary phenomenon with steep gradients and sharp peaks in unknowns is described as a coupled problem of convection-diffusion equation with source term for electron, ion densities and Poisson's equation for electric potential. The numerical method is 2nd order of accuracy in space and time and it uses dynamical adaptation of unstructured triangular mesh. Results of numerical studies included size of computational domain, type of boundary conditions and numerical convergence test are presented.
FVCA4, HPS, F. Benkhaldoun, D. Ouazar ( …, 2005
Flow, Turbulence and Combustion, 2006
In this work we consider a two steps finite volume scheme, recently developed to solve nonhomogen... more In this work we consider a two steps finite volume scheme, recently developed to solve nonhomogeneous systems. The first step of the scheme depends on a diffusion control parameter which we modulate, using the limiters theory. Results on Shallow water equations and two phase flows are presented. Keywords Non homogeneous systems . Finite volumes . SRNHR scheme F. Benkhaldoun ( ) • L.
Topical Problems of Fluid Mechanics 2017, 2017
Fluids
In this paper, a new rheological model for the flow of phosphate-water suspensions is proposed. T... more In this paper, a new rheological model for the flow of phosphate-water suspensions is proposed. The model’s ability to replicate the rheological characteristics of phosphate-water suspensions under different shear rate conditions is evaluated using rheometric tests, and it is found to be in good agreement with experimental data. A comprehensive methodology for obtaining the model parameters is presented. The proposed model is then incorporated into the OpenFoam numerical code. The results demonstrate that the model is capable of reproducing the rheological behavior of phosphate suspensions at both low and high concentrations by comparing it with suitable models for modeling the rheological behavior of phosphate suspensions. The proposed model can be applied to simulate and monitor phosphate slurry flows in industrial applications.

Algorithms
In this paper, an innovative methodology to handle the numerical simulation of viscoplastic flows... more In this paper, an innovative methodology to handle the numerical simulation of viscoplastic flows is proposed based on a multigrid initialization algorithm in conjunction with the SIMPLE procedure. The governing equations for incompressible flow, which consist of continuity and momentum equations, are solved on a collocated grid by combining the finite volume discretization and Rhie and chow interpolation for pressure–velocity coupling. Using the proposed solver in combination with the regularization scheme of Papanastasiou, we chose the square lid-driven cavity flow and pipe flow as test cases for validation and discussion. In doing so, we study the influence of the Bingham number and the Reynolds number on the development of rigid areas and the features of the vortices within the flow domain. Pipe flow results illustrate the flow’s response to the stress growth parameter values. We show that the representation of the yield surface and the plug zone is influenced by the chosen valu...

arXiv (Cornell University), May 24, 2021
The present paper describes a parallel unstructured-mesh Plasma simulation code based on Finite V... more The present paper describes a parallel unstructured-mesh Plasma simulation code based on Finite Volume method. The code dynamically refines and coarses mesh for accurate resolution of the different features regarding the electron density. Our purpose is to examine the performance of a new Parallel Adaptive Mesh Refinement (PAMR) procedure introduced on the ADAPT platform, which resolves of a relatively complicated system coupling the flow partial differential equations to the Poisson's equation. The implementation deals with the MUMPS parallel multi-frontal direct solver and mesh partitioning methods using METIS to improve the performance of the framework. The standard MPI is used to establish communication between processors. Performance analysis of the PAMR procedure shows the efficiency and the potential of the method for the propagation equations of ionization waves.
Springer Proceedings in Mathematics & Statistics, 2017
Adaptive unstructured finite volume methods for ionization waves are receiving increased attentio... more Adaptive unstructured finite volume methods for ionization waves are receiving increased attention mainly because of their ability to provide a flexible spatial discretization. Hence, some areas can be resolved in great detail while not over-resolving other areas. Our purpose is to examine the numerical performance of a new criteria for mesh adaptation which account only for the elliptic equation for the electric potential. The proposed adaptive finite volume method has important advantages in the discretization of the gradient fluxes and diffusion terms using unstructured grids and satisfies the conservation property. Numerical results are presented for a propagation of ionization waves in a rectangular domain.
Finite Volumes for Complex Applications IX - Methods, Theoretical Aspects, Examples, 2020
We present a finite volume model for the simulation of floods in urban areas. The model consists ... more We present a finite volume model for the simulation of floods in urban areas. The model consists of the two-dimensional shallow water equations with variable horizontal porosity which is introduced in order to reflect the effects of obstructions. An extra porosity source term appears in the momentum equations. The main advantage of this model is the significant reduction of the computational cost while preserving an acceptable level of accuracy. The finite volume method uses a modified Roe’s scheme involving the sign of the Jacobian matrix in the system for the discretization of gradient fluxes. The performance of the numerical model is demonstrated by comparing the results obtained using the proposed method to laboratory experiments for a flow problem over an array of obstacles.
Uploads
Papers by Fayssal Benkhaldoun