Denys Dutykh

Didier Clamond, Denys Dutykh

This paper describes an efficient algorithm for computing steady two-dimensional surface gravity wave in irrotational motion. The algorithm complexity is O(N log N), N being the number of Fourier modes. The algorithm allows the arbitrary precision computation of waves in arbitrary depth, i.e., it works efficiently for Stokes, cnoidal and solitary waves, even for quite large steepnesses. The method is based on conformal mapping, Babenko equation rewritten in a suitable way, pseudo-spectral method and Petviashvili's iterations. The efficiency of the algorithm is illustrated via some relevant numerical examples. The code is open source, so interested readers can easily check the claims, use and modify the algorithm.

Denys Dutykh, Didier Clamond, Angel Duran

This paper is devoted to the computation of capillary-gravity solitary waves of the irrotational incompressible Euler equations with free surface. The numerical study is a continuation of a previous work in several points: an alternative formulation of the Babenko-type equation for the wave profiles, a detailed description of both the numerical resolution and the analysis of the internal flow structure under a solitary wave. The numerical code used in this study is provided in open source for those interested readers.

Denys Dutykh, Raphaël Poncet, Frédéric Dias

A novel tool for tsunami wave modelling is presented. This tool has the potential of being used for operational purposes: indeed, the numerical code \VOLNA is able to handle the complete life-cycle of a tsunami (generation, propagation and run-up along the coast). The algorithm works on unstructured triangular meshes and thus can be run in arbitrary complex domains. This paper contains the detailed description of the finite volume scheme implemented in the code. The numerical treatment of the wet/dry transition is explained. This point is crucial for accurate run-up/run-down computations. Most existing tsunami codes use semi-empirical techniques at this stage, which are not always sufficient for tsunami hazard mitigation. Indeed the decision to evacuate inhabitants is based on inundation maps which are produced with this type of numerical tools. We present several realistic test cases that partially validate our algorithm. Comparisons with analytical solutions and experimental data are performed. Finally the main conclusions are outlined and the perspectives for future research presented.

Denys Dutykh, Dimitrios Mitsotakis, Xavier Gardeil et al.

The present study is devoted to the problem of tsunami wave generation. The main goal of this work is two-fold. First of all, we propose a simple and computationally inexpensive model for the description of the sea bed displacement during an underwater earthquake, based on the finite fault solution for the slip distribution under some assumptions on the dynamics of the rupturing process. Once the bottom motion is reconstructed, we study waves induced on the free surface of the ocean. For this purpose we consider three different models approximating the Euler equations of the water wave theory. Namely, we use the linearized Euler equations (we are in fact solving the Cauchy-Poisson problem), a Boussinesq system and a novel weakly nonlinear model. An intercomparison of these approaches is performed. The developments of the present study are illustrated on the 17 July 2006 Java event, where an underwater earthquake of magnitude 7.7 generated a tsunami that inundated the southern coast of Java.

Hayk Nersisyan, Denys Dutykh, Enrique Zuazua

We investigate the potential and limitations of the wave generation by disturbances moving at the bottom. More precisely, we assume that the wavemaker is composed of an underwater object of a given shape which can be displaced according to a prescribed trajectory. We address the practical question of computing the wavemaker shape and trajectory generating a wave with prescribed characteristics. For the sake of simplicity we model the hydrodynamics by a generalized forced Benjamin-Bona-Mahony (BBM) equation. This practical problem is reformulated as a constrained nonlinear optimization problem. Additional constraints are imposed in order to fulfill various practical design requirements. Finally, we present some numerical results in order to demonstrate the feasibility and performance of the proposed methodology.

Denys Dutykh, Theodoros Katsaounis, Dimitrios Mitsotakis

Finite volume schemes are commonly used to construct approximate solutions to conservation laws. In this study we extend the framework of the finite volume methods to dispersive water wave models, in particular to Boussinesq type systems. We focus mainly on the application of the method to bidirectional nonlinear, dispersive wave propagation in one space dimension. Special emphasis is given to important nonlinear phenomena such as solitary waves interactions, dispersive shock wave formation and the runup of breaking and non-breaking long waves.

Denys Dutykh, Frédéric Dias

The present article is devoted to the influence of sediment layers on the process of tsunami generation. The main scope here is to demonstrate and especially quantify the effect of sedimentation on vertical displacements of the seabed due to an underwater earthquake. The fault is modelled as a Volterra-type dislocation in an elastic half-space. The elastodynamics equations are integrated with a finite element method. A comparison between two cases is performed. The first one corresponds to the classical situation of an elastic homogeneous and isotropic half-space, which is traditionally used for the generation of tsunamis. The second test case takes into account the presence of a sediment layer separating the oceanic column from the hard rock. Some important differences are revealed. We conjecture that deformations in the generation region may be amplified by sedimentary deposits, at least for some parameter values. The mechanism of amplification is studied through careful numerical simulations.

Denys Dutykh, Didier Clamond, Paul Milewski et al.

After we derive the Serre system of equations of water wave theory from a generalized variational principle, we present some of its structural properties. We also propose a robust and accurate finite volume scheme to solve these equations in one horizontal dimension. The numerical discretization is validated by comparisons with analytical, experimental data or other numerical solutions obtained by a highly accurate pseudo-spectral method.

Denys Dutykh, Theodoros Katsaounis, Dimitrios Mitsotakis

We extend the framework of the finite volume method to dispersive unidirectional water wave propagation in one space dimension. In particular we consider a KdV-BBM type equation. Explicit and IMEX Runge-Kutta type methods are used for time discretizations. The fully discrete schemes are validated by direct comparisons to analytic solutions. Invariants conservation properties are also studied. Main applications include important nonlinear phenomena such as dispersive shock wave formation, solitary waves and their various interactions.

Suelen Gasparin, Julien Berger, Denys Dutykh et al.

Implicit schemes have been extensively used in building physics to compute the solution of moisture diffusion problems in porous materials for improving stability conditions. Nevertheless, these schemes require important sub-iterations when treating non-linear problems. To overcome this disadvantage, this paper explores the use of improved explicit schemes, such as Dufort-Frankel, Crank-Nicolson and hyperbolisation approaches. A first case study has been considered with the hypothesis of linear transfer. The Dufort-Frankel, Crank-Nicolson and hyperbolisation schemes were compared to the classical Euler explicit scheme and to a reference solution. Results have shown that the hyperbolisation scheme has a stability condition higher than the standard Courant-Friedrichs-Lewy (CFL) condition. The error of this schemes depends on the parameter τrepresenting the hyperbolicity magnitude added into the equation. The Dufort-Frankel scheme has the advantages of being unconditionally stable and is preferable for non-linear transfer, which is the second case study. Results have shown the error is proportional to O(Δt). A modified Crank-Nicolson scheme has been proposed in order to avoid sub-iterations to treat the non-linearities at each time step. The main advantages of the Dufort-Frankel scheme are (i) to be twice faster than the Crank-Nicolson approach; (ii) to compute explicitly the solution at each time step; (iii) to be unconditionally stable and (iv) easier to parallelise on high-performance computer systems. Although the approach is unconditionally stable, the choice of the time discretisation $Δt$ remains an important issue to accurately represent the physical phenomena.

Denys Dutykh, Marx Chhay, Francesco Fedele

Geometric discretizations that preserve certain Hamiltonian structures at the discrete level has been proven to enhance the accuracy of numerical schemes. In particular, numerous symplectic and multi-symplectic schemes have been proposed to solve numerically the celebrated Korteweg-de Vries (KdV) equation. In this work, we show that geometrical schemes are as much robust and accurate as Fourier-type pseudo-spectral methods for computing the long-time KdV dynamics, and thus more suitable to model complex nonlinear wave phenomena.

Denys Dutykh, Céline Acary-Robert, Didier Bresch

Powder-snow avalanches are violent natural disasters which represent a major risk for infrastructures and populations in mountain regions. In this study we present a novel model for the simulation of avalanches in the aerosol regime. The second scope of this study is to get more insight into the interaction process between an avalanche and a rigid obstacle. An incompressible model of two miscible fluids can be successfully employed in this type of problems. We allow for mass diffusion between two phases according to the Fick's law. The governing equations are discretized with a contemporary fully implicit finite volume scheme. The solver is able to deal with arbitrary density ratios. Several numerical results are presented. Volume fraction, velocity and pressure fields are presented and discussed. Finally we point out how this methodology can be used for practical problems.

Frédéric Dias, Denys Dutykh, Jean-Michel Ghidaglia

In the study of ocean wave impact on structures, one often uses Froude scaling since the dominant force is gravity. However the presence of trapped or entrained air in the water can significantly modify wave impacts. When air is entrained in water in the form of small bubbles, the acoustic properties in the water change dramatically and for example the speed of sound in the mixture is much smaller than in pure water, and even smaller than in pure air. While some work has been done to study small-amplitude disturbances in such mixtures, little work has been done on large disturbances in air-water mixtures. We propose a basic two-fluid model in which both fluids share the same velocities. It is shown that this model can successfully mimic water wave impacts on coastal structures. Even though this is a model without interface, waves can occur. Their dispersion relation is discussed and the formal limit of pure phases (interfacial waves) is considered. The governing equations are discretized by a second-order finite volume method. Numerical results are presented. It is shown that this basic model can be used to study violent aerated flows, especially by providing fast qualitative estimates.

Denys Dutykh, Didier Clamond

In the present study, we propose a modified version of the Nonlinear Shallow Water Equations (Saint-Venant or NSWE) for irrotational surface waves in the case when the bottom undergoes some significant variations in space and time. The model is derived from a variational principle by choosing an appropriate shallow water ansatz and imposing some constraints. Our derivation procedure does not explicitly involve any small parameter and is straightforward. The novel system is a non-dispersive non-hydrostatic extension of the classical Saint-Venant equations. A key feature of the new model is that, like the classical NSWE, it is hyperbolic and thus similar numerical methods can be used. We also propose a finite volume discretisation of the obtained hyperbolic system. Several test-cases are presented to highlight the added value of the new model. Some implications to tsunami wave modelling are also discussed.

Angel Duran, Denys Dutykh, Dimitrios Mitsotakis

Surface water waves in ideal fluids have been typically modeled by asymptotic approximations of the full Euler equations. Some of these simplified models lose relevant properties of the full water wave problem. One of them is the Galilean symmetry, which is not present in important models such as the BBM equation and the Peregrine (Classical Boussinesq) system. In this paper we propose a mechanism to modify the above mentioned classical models and derive new, Galilean invariant models. We present some properties of the new equations, with special emphasis on the computation and interaction of their solitary-wave solutions. The comparison with full Euler solutions shows the relevance of the preservation of Galilean invariance for the description of water waves.

Julien Berger, Hervé Le Meur, Denys Dutykh et al.

The reliability of a model is its accuracy in predicting the physical phenomena using the known input parameters. It also depends on the model's ability to estimate relevant parameters using observations of the physical phenomena. In this paper, the reliability of the VTT model is investigated under these two criteria for various given temperature and relative humidity constant in time. First of all, experiments are conducted on bamboo fiberboard. Using these data, five parameters of the VTT model, defining the mold vulnerability class of a material, are identified. The results highlight that the determined parameters are not within the range of the classes defined in the VTT model. In addition, the quality of the parameter estimation is not satisfactory. Then the sensitivity of the numerical results of the VTT model is analyzed by varying an input parameter. These investigations show that the VTT mathematical formulation of the physical model of mold growth is not reliable. An improved model is proposed with a new mathematical formulation. It is inspired by the logistic equation whose parameters are estimated using the experimental data obtained. The parameter estimation is very satisfactory. In the last parts of the paper, the numerical predictions of the improved model are compared to experimental data from the literature to prove its reliability.

Francesco Fedele, Denys Dutykh

Recently, Dyachenko & Zakharov (2011) have derived a compact form of the well known Zakharov integro-differential equation for the third order Hamiltonian dynamics of a potential flow of an incompressible, infinitely deep fluid with a free surface. Special traveling wave solutions of this compact equation are numerically constructed using the Petviashvili method. Their stability properties are also investigated. Further, unstable traveling waves with wedge-type singularities, viz. peakons, are numerically discovered. To gain insights into the properties of singular traveling waves, we consider the academic case of a perturbed version of the compact equation, for which analytical peakons with exponential shape are derived. Finally, by means of an accurate Fourier-type spectral scheme it is found that smooth solitary waves appear to collide elastically, suggesting the integrability of the Zakharov equation.

Yue Pu, Robert Pego, Denys Dutykh et al.

We study a regularization of the classical Saint-Venant (shallow-water) equations, recently introduced by D. Clamond and D. Dutykh (Commun. Nonl. Sci. Numer. Simulat. 55 (2018) 237-247). This regularization is non-dispersive and formally conserves mass, momentum and energy. We show that for every classical shock wave, the system admits a corresponding non-oscillatory traveling wave solution which is continuous and piecewise smooth, having a weak singularity at a single point where energy is dissipated as it is for the classical shock. The system also admits cusped solitary waves of both elevation and depression.

Julien Berger, Suelen Gasparin, Denys Dutykh et al.

Comparisons of experimental observation of heat and moisture transfer through porous building materials with numerical results have been presented in numerous studies reported in literature. However, some discrepancies have been observed, highlighting underestimation of sorption process and overestimation of desorption process. Some studies intend to explain the discrepancies by analysing the importance of hysteresis effects as well as carrying out sensitivity analyses on the input parameters as convective transfer coefficients. This article intends to investigate the accuracy and efficiency of the coupled solution by adding advective transfer of both heat and moisture in the physical model. In addition, the efficient Scharfetter and Gummel numerical scheme is proposed to solve the system of advection-diffusion equations, which has the advantages of being well-balanced and asymptotically preserving. Moreover, the scheme is particularly efficient in terms of accuracy and reduction of computational time when using large spatial discretisation parameters. Several linear and non-linear cases are studied to validate the method and highlight its specific features. At the end, an experimental benchmark from the literature is considered. The numerical results are compared to the experimental data for a pure diffusive model and also for the proposed model. The latter presents better agreement with the experimental data. The influence of the hysteresis effects on the moisture capacity is also studied, by adding a third differential equation.

Julien Berger, Denys Dutykh, Nathan Mendes

In the context of estimating material properties of porous walls based on in-site measurements and identification method, this paper presents the concept of Optimal Experiment Design (OED). It aims at searching the best experimental conditions in terms of quantity and position of sensors and boundary conditions imposed to the material. These optimal conditions ensure to provide the maximum accuracy of the identification method and thus the estimated parameters. The search of the OED is done by using the Fisher information matrix and a priori knowledge of the parameters. The methodology is applied for two case studies. The first one deals with purely conductive heat transfer. The concept of optimal experiment design is detailed and verified with 100 inverse problems for different experiment designs. The second case study combines a strong coupling between heat and moisture transfer through a porous building material. The methodology presented is based on a scientific formalism for efficient planning of experimental work that can be extended to the optimal design of experiments related to other problems in thermal and fluid sciences.

Marx Chhay, Denys Dutykh, Marguerite Gisclon et al.

In this article, we present a model of heat transfer occurring through a li\-quid film flowing down a vertical wall. This new model is formally derived using the method of asymptotic expansions by introducing appropriately chosen dimensionless variables. In our study the small parameter, known as the film parameter, is chosen as the ratio of the flow depth to the characteristic wavelength. A new Nusselt solution should be explained, taking into account the hydrodynamic free surface variations and the contributions of the higher order terms coming from temperature variation effects. Comparisons are made with numerical solutions of the full Fourier equations in a steady state frame. The flow and heat transfer are coupled through Marangoni and temperature dependent viscosity effects. Even if these effects have been considered separately before, here a fully coupled model is proposed. Another novelty consists in the asymptotic approach in contrast to the weighted residual approach which have been formerly applied to these problems.

Denys Dutykh

Water wave propagation can be attenuated by various physical mechanisms. One of the main sources of wave energy dissipation lies in boundary layers. The present work is entirely devoted to thorough analysis of the dispersion relation of the novel visco-potential formulation. Namely, in this study we relax all assumptions of the weak dependence of the wave frequency on time. As a result, we have to deal with complex integro-differential equations that describe transient behaviour of the phase and group velocities. Using numerical computations, we show several snapshots of these important quantities at different times as functions of the wave number. Good qualitative agreement with previous study [Dutykh2009] is obtained. Thus, we validate in some sense approximations made anteriorly. There is an unexpected conclusion of this study. According to our computations, the bottom boundary layer creates disintegrating modes in the group velocity. In the same time, the imaginary part of the phase velocity remains negative for all times. This result can be interpreted as a new kind of instability which is induced by the bottom boundary layer effect.

Gayaz Khakimzyanov, Denys Dutykh, Oleg Gusev

In the present manuscript, we consider the problem of dispersive wave simulation on a rotating globally spherical geometry. In this Part IV, we focus on numerical aspects while the model derivation was described in Part III. The algorithm we propose is based on the splitting approach. Namely, equations are decomposed on a uniformly elliptic equation for the dispersive pressure component and a hyperbolic part of shallow water equations (on a sphere) with source terms. This algorithm is implemented as a two-step predictor-corrector scheme. On every step, we solve separately elliptic and hyperbolic problems. Then, the performance of this algorithm is illustrated on model idealised situations with an even bottom, where we estimate the influence of sphericity and rotation effects on dispersive wave propagation. The dispersive effects are quantified depending on the propagation distance over the sphere and on the linear extent of generation region. Finally, the numerical method is applied to a couple of real-world events. Namely, we undertake simulations of the Bulgarian 2007 and Chilean 2010 tsunamis. Whenever the data is available, our computational results are confronted with real measurements.

Denys Dutykh, Theodoros Katsaounis, Dimitrios Mitsotakis

Historically the finite volume methods have been developed for the numerical integration of conservation laws. In this study we present some recent results on the application of such schemes to dispersive PDEs. Namely, we solve numerically a representative of Boussinesq type equations in view of important applications to the coastal hydrodynamics. Numerical results of the runup of a moderate wave onto a non-uniform beach are presented along with great lines of the employed numerical method (see D. Dutykh et al. (2011) for more details).

Angel Durán, Denys Dutykh, Dimitrios Mitsotakis

In 1967 D. H. Peregrine proposed a Boussinesq-type model for long waves in shallow waters of varying depth. This prominent paper turned a new leaf in coastal hydrodynamics along with contributions by F. Serre, A. E. Green \& P. M. Naghdi and many others since then. Several modern Boussinesq-type systems stem from these pioneering works. In the present work we revise the long wave model traditionally referred to as the Peregrine system. Namely, we propose a modification of the governing equations which is asymptotically similar to the initial model for weakly nonlinear waves, while preserving an additional symmetry of the complete water wave problem. This modification procedure is called the invariantization. We show that the improved system has well conditioned dispersive terms in the swash zone, hence allowing for efficient and stable run-up computations.

Julien Berger, Denys Dutykh, Nathan Mendes et al.

This work presents a detailed mathematical model combined with an innovative efficient numerical model to predict heat, air and moisture transfer through porous building materials. The model considers the transient effects of air transport and its impact on the heat and moisture transfer. The achievement of the mathematical model is detailed in the continuity of Luikov's work. A system composed of two advection-diffusion differential equations plus one exclusively diffusion equation is derived. The main issue to take into account the transient air transfer arises in the very small characteristic time of the transfer, implying very fine discretisation. To circumvent these difficulties, the numerical model is based on the Du Fort-Frankel explicit and unconditionally stable scheme for the exclusively diffusion equation. It is combined with a two-step Runge-Kutta scheme in time with the Scharfetter-Gummel numerical scheme in space for the coupled advection-diffusion equations. At the end, the numerical model enables to relax the stability condition, and, therefore, to save important computational efforts. A validation case is considered to evaluate the efficiency of the model for a nonlinear problem. Results highlight a very accurate solution computed about 16 times faster than standard approaches. After this numerical validation, the reliability of the mathematical model is evaluated by comparing the numerical predictions to experimental observations. The latter is measured within a multi-layered wall submitted to a sudden increase of vapor pressure on the inner side and driven climate boundary conditions on the outer side. A very satisfactory agreement is noted between the numerical predictions and experimental observations indicating an overall good reliability of the proposed model.

Denys Dutykh

This manuscript contains some thoughts on the discretization of the classical heat equation. Namely, we discuss the advantages and disadvantages of explicit and implicit schemes. Then, we show how to overcome some disadvantages while preserving some advantages. However, since there is no free lunch, there is a price to pay for any improvement in the numerical scheme. This price will be thoroughly discussed below. In particular, we like explicit discretizations for the ease of their implementation even for nonlinear problems. Unfortunately, when these schemes are applied to parabolic equations, severe stability limits appear for the time step magnitude making the explicit simulations prohibitively expensive. Implicit schemes remove the stability limit, but each time step requires now the solution of linear (at best) or even nonlinear systems of equations. However, there exists a number of tricks to overcome (or at least to relax) severe stability limitations of explicit schemes without going into the trouble of fully implicit ones. The purpose of this manuscript is just to inform the readers about these alternative techniques to extend the stability limits. It was not written for classical scientific publication purposes.

Suelen Gasparin, Denys Dutykh, Nathan Mendes

This work presents an efficient numerical method based on spectral expansions for simulation of heat and moisture diffusive transfers through multilayered porous materials. Traditionally, by using the finite-difference approach, the problem is discretized in time and space domains (Method of lines) to obtain a large system of coupled Ordinary Differential Equations (ODEs), which is computationally expensive. To avoid such a cost, this paper proposes a reduced-order method that is faster and accurate, using a much smaller system of ODEs. To demonstrate the benefits of this approach, tree case studies are presented. The first one considers nonlinear heat and moisture transfer through one material layer. The second case - highly nonlinear - imposes a high moisture content gradient - simulating a rain like condition - over a two-layered domain, while the last one compares the numerical prediction against experimental data for validation purposes. Results show how the nonlinearities and the interface between materials are easily and naturally treated with the spectral reduced-order method. Concerning the reliability part, predictions show a good agreement with experimental results, which confirm robustness, calculation efficiency and high accuracy of the proposed approach for predicting the coupled heat and moisture transfer through porous materials.

Denys Dutykh, Jean-Louis Verger-Gaugry

The class B of lacunary polynomials f(x) := -1 + x + x^n + x^{m_1} + x^{m_2} + ... + x^{m_s}, where s >= 0, m_1 - n >= n - 1, m_{q+1} - m_{q} >= n - 1 for 1 <= q < s, n >= 3 is studied. A polynomial having its coefficients in {0, 1} except its constant coefficient equal to -1 is called an almost Newman polynomial. A general theorem of factorization of the almost Newman polynomials of the class B is obtained. Such polynomials possess lenticular roots in the open unit disk off the unit circle in the small angular sector π/18 <= arg z <= π/18 and their nonreciprocal parts are always irreducible. The existence of lenticuli of roots is a peculiarity of the class B. By comparison with the Odlyzko - Poonen Conjecture and its variant Conjecture, an `Asymptotic Reducibility Conjecture' is formulated aiming at establishing the proportion of irreducible polynomials in this class. This proportion is conjectured to be 3/4 and estimated using Monte-Carlo methods. The numerical approximate value ~ 0.756 is obtained. The results extend those on trinomials (Selmer) and quadrinomials (Ljunggren, Mills, Finch and Jones).

Julien Berger, Thomas Busser, Denys Dutykh et al.

This paper deals with an inverse problem applied to the field of building physics to experimentally estimate three sorption isotherm coefficients of a wood fiber material. First, the mathematical model, based on convective transport of moisture, the Optimal Experiment Design (OED) and the experimental set-up are presented. Then measurements of relative humidity within the material are carried out, after searching the OED, which is based on the computation of the sensitivity functions and a priori values of the unknown parameters employed in the mathematical model. The OED enables to plan the experimental conditions in terms of sensor positioning and boundary conditions out of 20 possible designs, ensuring the best accuracy for the identification method and, thus, for the estimated parameter. Two experimental procedures were identified: i) single step of relative humidity from 10% to 75% and ii) multiple steps of relative humidity 10-75-33-75% with an 8-day duration period for each step. For both experiment designs, it has been shown that the sensor has to be placed near the impermeable boundary. After the measurements, the parameter estimation problem is solved using an interior point algorithm to minimize the cost function. Several tests are performed for the definition of the cost function, by using the L^2 or L^\infty norm and considering the experiments separately or at the same time. It has been found out that the residual between the experimental data and the numerical model is minimized when considering the discrete Euclidean norm and both experiments separately. It means that two parameters are estimated using one experiment while the third parameter is determined with the other experiment. Two cost functions are defined and minimized for this approach. Moreover, the algorithm requires less than 100 computations of the direct model to obtain the solution. In addition, the OED sensitivity functions enable to capture an approximation of the probability distribution function of the estimated parameters. The determined sorption isotherm coefficients calibrate the numerical model to fit better the experimental data. However, some discrepancies still appear since the model does not take into account the hysteresis effects on the sorption capacity. Therefore, the model is improved proposing a second differential equation for the sorption capacity to take into account the hysteresis between the main adsorption and desorption curves. The OED approach is also illustrated for the estimation of five of the coefficients involved in the hysteresis model. To conclude, the prediction of the model with hysteresis are compared with the experimental observations to illustrate the improvement of the prediction.

Julien Berger, Denys Dutykh

The fidelity of a model relies both on its accuracy to predict the physical phenomena and its capability to estimate unknown parameters using observations. This article focuses on this second aspect by analyzing the reliability of two mathematical models proposed in the literature for the simulation of heat losses through building walls. The first one, named DuFort-Frankel (DF), is the classical heat diffusion equation combined with the DuFort-Frankel numerical scheme. The second is the so-called RC lumped approach, based on a simple ordinary differential equation to compute the temperature within the wall. The reliability is evaluated following a two stages method. First, samples of observations are generated using a pseudo-spectral numerical model for the heat diffusion equation with known input parameters. The results are then modified by adding a noise to simulate experimental measurements. Then, for each sample of observation, the parameter estimation problem is solved using one of the two mathematical models. The reliability is assessed based on the accuracy of the approach to recover the unknown parameter. Three case studies are considered for the estimation of (i) the heat capacity, (ii) the thermal conductivity or (iii) the heat transfer coefficient at the interface between the wall and the ambient air. For all cases, the DF mathematical model has a very satisfactory reliability to estimate the unknown parameters without any bias. However, the RC model lacks of fidelity and reliability. The error on the estimated parameter can reach 40% for the heat capacity, 80% for the thermal conductivity and 450% for the heat transfer coefficient.

Ainagul Jumabekova, Julien Berger, Denys Dutykh et al.

The goal of this study is to propose an efficient numerical model for the predictions of capillary adsorption phenomena in a porous material. The Scharfetter-Gummel numerical scheme is proposed to solve an advection-diffusion equation with gravity flux. Its advantages such as accuracy, relaxed stability condition, and reduced computational cost are discussed along with the study of linear and nonlinear cases. The reliability of the numerical model is evaluated by comparing the numerical predictions with experimental observations of liquid uptake in bricks. A parameter estimation problem is solved to adjust the uncertain coefficients of moisture diffusivity and hydraulic conductivity.

Francesco Fedele, Denys Dutykh

We present a study on the nonlinear dynamics of a disturbance to the laminar state in non-rotating axisymmetric Poiseuille pipe flows. The associated Navier-Stokes equations are reduced to a set of coupled generalized Camassa-Holm type equations. These support singular inviscid travelling waves with wedge-type singularities, the so called peakons, which bifurcate from smooth solitary waves as their celerity increase. In physical space they correspond to localized toroidal vortices or vortexons. The inviscid vortexon is similar to the nonlinear neutral structures found by Walton (2011) and it may be a precursor to puffs and slugs observed at transition, since most likely it is unstable to non-axisymmetric disturbances.

Srijani Mukherjee, Laurent Vuillon, Liliane Bou Nassif et al.

The rapid growth of solar photovoltaic (PV) systems necessitates advanced methods for performance monitoring and anomaly detection to ensure optimal operation. In this study, we propose a novel approach leveraging Temporal Graph Neural Network (Temporal GNN) to predict solar PV output power and detect anomalies using environmental and operational parameters. The proposed model utilizes graph-based temporal relationships among key PV system parameters, including irradiance, module and ambient temperature to predict electrical power output. This study is based on data collected from an outdoor facility located on a rooftop in Lyon (France) including power measurements from a PV module and meteorological parameters.

Denys Dutykh

In a recent study [DutykhDias2007] we presented a novel visco-potential free surface flows formulation. The governing equations contain local and nonlocal dissipative terms. From physical point of view, local dissipation terms come from molecular viscosity but in practical computations, rather eddy viscosity should be used. On the other hand, nonlocal dissipative term represents a correction due to the presence of a bottom boundary layer. Using the standard procedure of Boussinesq equations derivation, we come to nonlocal long wave equations. In this article we analyse dispersion relation properties of proposed models. The effect of nonlocal term on solitary and linear progressive waves attenuation is investigated. Finally, we present some computations with viscous Boussinesq equations solved by a Fourier type spectral method.

Ahmad Deeb, Denys Dutykh, Maryam Al Zohbi

Integration of Ordinary Differential Equations (ODEs) using Backward Difference formula (BDF) methods with p backward steps achieves order p accuracy if specific conditions are met. This work extends the composition technique with complex coefficients to the implicit BDF schemes, increasing the approximation order by one without additional backward points. The imaginary part of the composed flow provides an error estimate of order p + 1. Linear stability analysis reveals that the composed schemes break the Dahlquist barrier, achieving stability up to order eight. The computational performance of the composed flow outperforms BDF schemes when using the same number of backward points, allowing for higher accuracy with lower CPU time. For non-uniform meshes, the ratio of consecutive time steps, which influences stability, appears as a parameter in the roots of algebraic equations relative to the composed flow. Having a complex root with a real positive part implies a lower bound to this ratio depending on the order. For example, the bound is 0.4506 for order three and 0.6806 for order four. Numerical tests demonstrate the effectiveness of this technique in improving the accuracy and stability compared to BDF methods.

Benoît Colange, Laurent Vuillon, Sylvain Lespinats et al.

To perform visual data exploration, many dimensionality reduction methods have been developed. These tools allow data analysts to represent multidimensional data in a 2D or 3D space, while preserving as much relevant information as possible. Yet, they cannot preserve all structures simultaneously and they induce some unavoidable distortions. Hence, many criteria have been introduced to evaluate a map's overall quality, mostly based on the preservation of neighbourhoods. Such global indicators are currently used to compare several maps, which helps to choose the most appropriate mapping method and its hyperparameters. However, those aggregated indicators tend to hide the local repartition of distortions. Thereby, they need to be supplemented by local evaluation to ensure correct interpretation of maps. In this paper, we describe a new method, called MING, for `Map Interpretation using Neighbourhood Graphs'. It offers a graphical interpretation of pairs of map quality indicators, as well as local evaluation of the distortions. This is done by displaying on the map the nearest neighbours graphs computed in the data space and in the embedding. Shared and unshared edges exhibit reliable and unreliable neighbourhood information conveyed by the mapping. By this mean, analysts may determine whether proximity (or remoteness) of points on the map faithfully represents similarity (or dissimilarity) of original data, within the meaning of a chosen map quality criteria. We apply this approach to two pairs of widespread indicators: precision/recall and trustworthiness/continuity, chosen for their wide use in the community, which will allow an easy handling by users.

Angel Durán, Denys Dutykh, Dimitrios Mitsotakis

In this paper we consider the numerical approximation of systems of Boussinesq-type to model surface wave propagation. Some theoretical properties of these systems (multi-symplectic and Hamiltonian formulations, well-posedness and existence of solitary-wave solutions) were previously analyzed by the authors in Part I. As a second part of the study, considered here is the construction of geometric schemes for the numerical integration. By using the method of lines, the geometric properties, based on the multi-symplectic and Hamiltonian structures, of different strategies for the spatial and time discretizations are discussed and illustrated.

Denys Dutykh

The topic of these notes could be easily expanded into a full one-semester course. Nevertheless, we shall try to give some flavour along with theoretical bases of spectral and pseudo-spectral methods. The main focus is made on Fourier-type discretizations, even if some indications on how to handle non-periodic problems via Tchebyshev and Legendre approaches are made as well. The applications presented here are diffusion-type problems in accordance with the topics of the PhD school.

Yannick Meyapin, Denys Dutykh, Marguerite Gisclon

In the present study we investigate analytically the process of velocity and energy relaxation in two-phase flows. We begin our exposition by considering the so-called six equations two-phase model [Ishii1975, Rovarch2006]. This model assumes each phase to possess its own velocity and energy variables. Despite recent advances, the six equations model remains computationally expensive for many practical applications. Moreover, its advection operator may be non-hyperbolic which poses additional theoretical difficulties to construct robust numerical schemes |Ghidaglia et al, 2001]. In order to simplify this system, we complete momentum and energy conservation equations by relaxation terms. When relaxation characteristic time tends to zero, velocities and energies are constrained to tend to common values for both phases. As a result, we obtain a simple two-phase model which was recently proposed for simulation of violent aerated flows [Dias et al, 2010]. The preservation of invariant regions and incompressible limit of the simplified model are also discussed. Finally, several numerical results are presented.

Denys Dutykh, Dimitrios Mitsotakis

The classical dam break problem has become the de facto standard in validating the Nonlinear Shallow Water Equations (NSWE) solvers. Moreover, the NSWE are widely used for flooding simulations. While applied mathematics community is essentially focused on developing new numerical schemes, we tried to examine the validity of the mathematical model under consideration. The main purpose of this study is to check the pertinence of the NSWE for flooding processes. From the mathematical point of view, the answer is not obvious since all derivation procedures assumes the total water depth positivity. We performed a comparison between the two-fluid Navier-Stokes simulations and the NSWE solved analytically and numerically. Several conclusions are drawn out and perspectives for future research are outlined.