Researcher profile

Denys Dutykh

Denys Dutykh contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author

131works

0followers

29topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Long wave interaction with a partially immersed body. Part II: Numerical results

In this manuscript we perform an extensive numerical study of the long wave interaction problem with a fixed partially immersed body into a fluid layer. The incident wave is assumed to be an isolated solitary wave. The body in this study is assumed to be fixed with a rectangular section which is not touching the bottom of the channel. The mathematical modelling of this problem is based on Part I (Khakimzyanov and Dutykh 2020) of this series and considered models include the Nonlinear Shallow Water Equations (NSWE), fully nonlinear weakly dispersive Serre-Green-Naghdi Equations (SGN equations) (completed with appropriate compatibility conditions on solid/fluid boundaries) and the free surface irrotational full Euler equations (FEE). We study the influence of the floating body elongation, immersion depth and incident wave amplitude on the wave field before and after the obstacle. The comparison of all three models predictions and the data of small-scale laboratory experiments is performed. Moreover, in the framework of the FEE model we investigate the anomalous wave run-up behind the floating body in the close presence of a vertical wall. We demonstrate the cases where the vertical wall creates extreme wave amplitudes behind the body, but also we show the cases where the wall attenuates wave amplitudes comparing to the wave field without a wall.

preprint2022arXiv

On a class of lacunary almost Newman polynomials modulo p and density theorems

The reduction modulo $p$ of a family of lacunary integer polynomials, associated with the dynamical zeta function $ζ_β(z)$ of the $β$-shift, for $β> 1$ close to one, is investigated. We briefly recall how this family is correlated to the problem of Lehmer. A variety of questions is raised about their numbers of zeroes in $\mathbb{F}_p$ and their factorizations, via Kronecker's Average Value Theorem (viewed as an analog of classical Theorems of Uniform Distribution Theory). These questions are partially answered using results of Schinzel, revisited by Sawin, Shusterman and Stoll, and density theorems (Frobenius, Chebotarev, Serre, Rosen). These questions arise from the search for the existence of integer polynomials of Mahler measure > 1 less than the smallest Salem number 1.176280. Explicit connection with modular forms (or modular representations) of the numbers of zeroes of these polynomials in $\mathbb{F}_p$ is obtained in a few cases. In general it is expected since it must exist according to the Langlands program.

preprint2020arXiv

A Computational Simulation of Steady Natural Convection in an H-form Cavity

The simulation of natural convection problem based on the Galerkin finite-element method, with the penalty finite-element formulation of the momentum balance equation, is exploited for accurate solutions of equations describing the problem of H-Form cavity differentially heated side walls. The cavity is occupied by the air whose Prandtl number is Pr=0.71, the fluid is assumed to be steady, viscous and incompressible within thermal convection. A numerical investigation has been made for Rayleigh numbers ranging from 10 to 10^6 for three cases of total internal height aspects of H-Form cavity: 0%, 50%, and 85%. Firstly, the goal is to validate the numerical code used to resolve the equations governing the problem of this work. For that, we present a comparison between the profiles at the point (0.5, 0) for the u-component, and u-component obtained in previous work for simple square cavity. Further, a comparison of the averaged Nusselt number with previous works for simple square cavity is realized in order to ensure the numerical accuracy, and the validity of our considered numerical tool. Secondly, the objective is to investigate on the hydrodynamic effects of Rayleigh number for different total internal height aspects of H-Form cavity on the dynamics of natural convection. Shortly after, the ambition is to assess the heat transfer rate for different Rayleigh number for three cases of internal height aspects.

preprint2020arXiv

An optimal scaling to computationally tractable dimensionless models: Study of latex particles morphology formation

In modelling of chemical, physical or biological systems it may occur that the coefficients, multiplying various terms in the equation of interest, differ greatly in magnitude, if a particular system of units is used. Such is, for instance, the case of the Population Balance Equations (PBE) proposed to model the Latex Particles Morphology formation. The obvious way out of this difficulty is the use of dimensionless scaled quantities, although often the scaling procedure is not unique. In this paper, we introduce a conceptually new general approach, called Optimal Scaling (OS). The method is tested on the known examples from classical and quantum mechanics, and applied to the Latex Particles Morphology model, where it allows us to reduce the variation of the relevant coefficients from 49 to just 4 orders of magnitudes. The PBE are then solved by a novel Generalised Method Of Characteristics, and the OS is shown to help reduce numerical error, and avoid unphysical behaviour of the solution. Although inspired by a particular application, the proposed scaling algorithm is expected find application in a wide range of chemical, physical and biological problems.

preprint2020arXiv

Coupling Conditions for Water Waves at Forks

We considered the propagation of nonlinear shallow water waves in a narrow channel presenting a fork. We aimed at computing the coupling conditions for a 1D effective model, using 2D simulations and an analysis based on the conservation laws. For small amplitudes, this analysis justifies the well-known Stoker interface conditions, so that the coupling does not depend on the angle of the fork. We also find this in the numerical solution. Large amplitude solutions in a symmetric fork also tend to follow Stoker's relations, due to the symmetry constraint. For non symmetric forks, 2D effects dominate so that it is necessary to understand the flow inside the fork. However, even then, conservation laws give some insight in the dynamics.

preprint2020arXiv

Derivation of an intermediate viscous Serre-Green-Naghdi equation

In this note we present the current status of the derivation of a viscous Serre-Green-Naghdi system. For this goal, the flow domain is separated into two regions. The upper region is governed by inviscid Euler equations, while the bottom region (the so-called boundary layer) is described by Navier-Stokes equations. We consider a particular regime linking the Reynolds number and the shallowness parameter. The computations presented in this note are performed in the fully nonlinear regime. The boundary layer flow reduces to a Prantdl-like equation. Further approximations seem to be needed to obtain a tractable model.

preprint2020arXiv

Dispersive shallow water wave modelling. Part I: Model derivation on a globally flat space

In this paper, we review the history and current state-of-the-art in the modelling of long nonlinear dispersive waves. For the sake of conciseness of this review, we omit the unidirectional models and focus especially on some classical and improved Boussinesq-type and Serre-Green-Naghdi equations. Finally, we propose also a unified modelling framework which incorporates several well-known and some less known dispersive wave models. The present manuscript is the first part of a series of two papers. The second part will be devoted to the numerical discretization of a practically important model on moving adaptive grids.

preprint2020arXiv

Dispersive shallow water wave modelling. Part II: Numerical simulation on a globally flat space

In this paper, we describe a numerical method to solve numerically the weakly dispersive fully nonlinear Serre-Green-Naghdi (SGN) celebrated model. Namely, our scheme is based on reliable finite volume methods, proven to be very effective for the hyperbolic part of equations. The particularity of our study is that we develop an adaptive numerical model using moving grids. Moreover, we use a special form of the SGN equations where non-hydrostatic part of pressure is found by solving a nonlinear elliptic equation. Moreover, this form of governing equations allows determining the natural form of boundary conditions to obtain a well-posed (numerical) problem.

preprint2020arXiv

Energy of tsunami waves generated by bottom motion

In the vast literature on tsunami research, few articles have been devoted to energy issues. A theoretical investigation on the energy of waves generated by bottom motion is performed here. We start with the full incompressible Euler equations in the presence of a free surface and derive both dispersive and non-dispersive shallow-water equations with an energy equation. It is shown that dispersive effects only appear at higher order in the energy budget. Then we solve the Cauchy-Poisson problem of tsunami generation for the linearized water wave equations. Exchanges between potential and kinetic energies are clearly revealed.

preprint2020arXiv

Extreme Inundation Statistics on a Composite Beach

The runup of initial Gaussian narrow-banded and wide-banded wave fields and its statistical characteristics are investigated using direct numerical simulations, based on the nonlinear shallow water equations. The bathymetry consists of the section of a constant depth, which is matched with the beach of constant slope. To address different levels of nonlinearity, the time series with five different significant wave heights are considered. The selected wave parameters allow also seeing the effects of wave breaking on wave statistics. The total physical time of each simulated time-series is 1000 hours (~360000 wave periods). The statistics of calculated wave runup heights are discussed with respect to the wave nonlinearity, wave breaking and the bandwidth of the incoming wave field. The conditional Weibull distribution is suggested as a model for the description of extreme runup heights and assessment of extreme inundations.

preprint2020arXiv

Extreme wave runup on a vertical cliff

Wave impact and run-up onto vertical obstacles are among the most important phenomena which must be taken into account in the design of coastal structures. From linear wave theory, we know that the wave amplitude on a vertical wall is twice the incident wave amplitude with weakly nonlinear theories bringing small corrections to this result. In this present study, however, we show that certain simple wave groups may produce much higher run-ups than previously predicted, with particular incident wave frequencies resulting in run up heights exceeding the initial wave amplitude by a factor of 5, suggesting that the notion of the design wave used in coastal structure design may need to be revisited. The results presented in this study can be considered as a note of caution for practitioners, on one side, and as a challenging novel material for theoreticians who work in the field of extreme wave - coastal structure interaction.

preprint2020arXiv

Free Surface Flows in Electrohydrodynamics with a Constant Vorticity Distribution

In 1895, Korteweg and de Vries (KdV), derived their celebrated equation describing the motion of waves of long wavelength in shallow water. In doing so they made a number of quite reasonable assumptions, incompressibility of the water and irrotational fluid. The resulting equation, the celebrated KdV equation, has been shown to be a very reasonable description of real water waves. However there are other phenomena which have an impact on the shape of the wave, that of vorticity and viscosity. This paper examines how a constant vorticity affects the shape of waves in electrohydrodynamics. For constant vorticity, the vertical component of the velocity obeys a Laplace equation and also has the usual lower boundary condition. In making the vertical component of the velocity take central stage, the Burns condition can be thus bypassed.

preprint2020arXiv

How to overcome the Courant-Friedrichs-Lewy condition of explicit discretizations?

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.

preprint2020arXiv

Nonlinear waves in networks: a simple approach using the sine-Gordon equation

To study the propagation of nonlinear waves across Y- and T-type junctions, we consider the 2D sine--Gordon equation as a model and study the dynamics of kinks and breathers in such geometries. The comparison of the energies reveals that the angle of the fork plays no role. Motivated by this, we introduce a 1D effective equation whose solutions agree well with the 2D simulations for kink and breather solutions. For branches of equal width, breather crossing occurs approximately when $v > 1 - ω$, where $v$ is the breather celerity and $ω$ is its frequency. We then characterize the breathers in the two upper branches by estimating their velocity and frequency. These new breathers are slower than the initial breather and up-shifted in frequency. In perspective, this study could be generalized to more complex nonlinear waves.

preprint2020arXiv

Numerical simulation of conservation laws with moving grid nodes: Application to tsunami wave modelling

In the present article we describe a few simple and efficient finite volume type schemes on moving grids in one spatial dimension combined with appropriate predictor-corrector method to achieve higher resolution. The underlying finite volume scheme is conservative and it is accurate up to the second order in space. The main novelty consists in the motion of the grid. This new dynamic aspect can be used to resolve better the areas with large solution gradients or any other special features. No interpolation procedure is employed, thus unnecessary solution smearing is avoided, and therefore, our method enjoys excellent conservation properties. The resulting grid is completely redistributed according the choice of the so-called monitor function. Several more or less universal choices of the monitor function are provided. Finally, the performance of the proposed algorithm is illustrated on several examples stemming from the simple linear advection to the simulation of complex shallow water waves. The exact well-balanced property is proven. We believe that the techniques described in our paper can be beneficially used to model tsunami wave propagation and run-up.

preprint2020arXiv

Numerical simulation of Feller's diffusion equation

This article is devoted to Feller's diffusion equation which arises naturally in probabilities and physics (e.g. wave turbulence theory). If discretized naively, this equation may represent serious numerical difficulties since the diffusion coefficient is practically unbounded and most of its solutions are weakly divergent at the origin. In order to overcome these difficulties we reformulate this equation using some ideas from the Lagrangian fluid mechanics. This allows us to obtain a numerical scheme with a rather generous stability condition. Finally, the algorithm admits an elegant implementation and the corresponding Matlab code is provided with this article under an open source license.

preprint2020arXiv

Numerical simulation of wave impact on a rigid wall using a two-phase compressible SPH method

In this paper, an SPH method based on the SPH-ALE formulation is used for modelling two-phase flows with large density ratios and realistic sound speeds. The SPH scheme is further improved to circumvent the tensile instability that may occur in the SPH simulations. The two-phase SPH solver is then used to model a benchmark problem of liquid impact on a rigid wall. The results are compared with an incompressible Level Set solver. Furthermore, a wave impact on a rigid wall with a large entrained air pocket is modelled. The SPH simulation is initialised by the output of a fully non-linear potential flow solver. The pressure distribution, velocity field and impact pressure are then analysed.

preprint2020arXiv

Observation of the Inverse Energy Cascade in the modified Korteweg-de Vries Equation

In this Letter we demonstrate for the first time the formation of the inverse energy cascade in the focusing modified Korteweg-de Vries (mKdV) equation. We study numerically the properties of this cascade such as the dependence of the spectrum shape on the initial excitation parameter (amplitude), perturbation magnitude and the size of the spectral domain. Most importantly we found that the inverse cascade is always accompanied by the direct one and they both form a very stable quasi-stationary structure in the Fourier space in the spirit of the FPU-like reoccurrence phenomenon. The formation of this structure is intrinsically related to the development of the nonlinear stage of the Modulational Instability (MI). These results can be used in several fields such as the internal gravity water waves, ion-acoustic waves in plasmas and others.

preprint2020arXiv

Poisson and Symplectic structures, Hamiltonian action, momentum and reduction

This manuscript is essentially a collection of lecture notes which were given by the first author at the Summer School Wisl-2019, Poland and written down by the second author. As the title suggests, the material covered here includes the Poisson and symplectic structures (Poisson manifolds, Poisson bi-vectors and Poisson brackets), group actions and orbits (infinitesimal action, stabilizers and adjoint representations), moment maps, Poisson and Hamiltonian actions. Finally, the phase space reduction is also discussed. The very last section introduces the Poisson-Lie structures along with some related notions. This text represents a brief review of a well-known material citing standard references for more details. The exposition is concise, but pedagogical. The Authors believe that it will be useful as an introductory exposition for students interested in this specific topic.

preprint2020arXiv

Resonance enhancement by suitably chosen frequency detuning

In this Letter we report new effects of resonance detuning on various dynamical parameters of a generic 3-wave system. Namely, for suitably chosen values of detuning the variation range of amplitudes can be significantly wider than for exact resonance. Moreover, the range of energy variation is not symmetric with respect to the sign of the detuning. Finally, the period of the energy oscillation exhibits non-monotonic dependency on the magnitude of detuning. These results have important theoretical implications where nonlinear resonance analysis is involved, such as geophysics, plasma physics, fluid dynamics. Numerous practical applications are envisageable e.g. in energy harvesting systems.

preprint2020arXiv

The Whitham Equation with Surface Tension

The viability of the Whitham equation as a nonlocal model for capillary-gravity waves at the surface of an inviscid incompressible fluid is under study. A nonlocal Hamiltonian system of model equations is derived using the Hamiltonian structure of the free surface water wave problem and the Dirichlet-Neumann operator. The system features gravitational and capillary effects, and when restricted to one-way propagation, the system reduces to the capillary Whitham equation. It is shown numerically that in various scaling regimes the Whitham equation gives a more accurate approximation of the free-surface problem for the Euler system than other models like the KdV, and Kawahara equation. In the case of relatively strong capillarity considered here, the KdV and Kawahara equations outperform the Whitham equation with surface tension only for very long waves with negative polarity.

preprint2019arXiv

A comparative study of bi-directional Whitham systems

In 1967, Whitham proposed a simplified surface water-wave model which combined the full linear dispersion relation of the full Euler equations with a weakly linear approximation. The equation he postulated which is now called the Whitham equation has recently been extended to a system of equations allowing for bi-directional propagation of surface waves. A number of different two-way systems have been put forward, and even though they are similar from a modeling point of view, these systems have very different mathematical properties. In the current work, we review some of the existing fully dispersive systems. We use state-of-the-art numerical tools to try to understand existence and stability of solutions to the initial-value problem associated to these systems. We also put forward a new system which is Hamiltonian and semi-linear. The new system is shown to perform well both with regard to approximating the full Euler system, and with regard to well posedness properties.

preprint2019arXiv

A new model for simulating heat, air and moisture transport in porous building materials

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.

preprint2019arXiv

A spectral method for solving heat and moisture transfer through consolidated porous media

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.

preprint2019arXiv

An adaptive simulation of nonlinear heat and moisture transfer as a boundary value problem

This work presents an alternative view on the numerical simulation of diffusion processes applied to the heat and moisture transfer through porous building materials. Traditionally, by using the finite-difference approach, the discretization follows the Method Of Lines (MOL), when the problem is first discretized in space to obtain a large system of coupled Ordinary Differential Equations (ODEs). Thus, this paper proposes to change this viewpoint. First, we discretize in time to obtain a small system of coupled ODEs, which means instead of having a Cauchy (Initial Value) Problem (IVP), we have a Boundary Value Problem (BVP). Fortunately, BVPs can be solved efficiently today using adaptive collocation methods of high order. To demonstrate the benefits of this new approach, three case studies are presented, in which one of them is compared with experimental data. The first one considers nonlinear heat and moisture transfer through one material layer while the second one considers two material layers. Results show how the nonlinearities and the interface between materials are easily treated, by reasonably using a fourth-order adaptive method. Finally, the last case study compares numerical results with experimental measurements, showing a good agreement.

preprint2019arXiv

An efficient numerical method for a long-term simulation of heat and mass transfer: the case of an insulated rammed earth wall

Innovative numerical scheme studied in this work enables to overcome two main limitations of Building Performance Simulation (BPS) programs as high computational cost and the choice of a very fine numerical grid. The method, called Super-Time-Stepping (STS), is novel to the state-of-the-art of building simulations, but has already proved to be sufficiently efficient in recent studies from anisotropic heat conduction in astrophysics (Meyer et al. 2014). The given research is focused on employment of this adopted numerical method to model drying of a rammed earth wall with an additional insulation layer. The results show considerable advantage of the STS method compared to standard Euler explicit scheme. It is possible to choose at least 100 times bigger time-steps to maintain high accuracy and to cut computational cost by more than 92% in the same time.

preprint2019arXiv

An efficient numerical model for liquid water uptake in porous material and its parameter estimation

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.

preprint2019arXiv

An innovative method to determine optimum insulation thickness based on non-uniform adaptive moving grid

It is well known that thermal insulation is a leading strategy for reducing energy consumption associated to heating or cooling processes in buildings. Nevertheless, building insulation can generate high expenditures so that the selection of an optimum insulation thickness requires a detailed energy simulation as well as an economic analysis. In this way, the present study proposes an innovative non-uniform adaptive method to determine the optimal insulation thickness of external walls. First, the method is compared with a reference solution to properly understand the features of the method, which can provide high accuracy with less spatial nodes. Then, the adaptive method is used to simulate the transient heat conduction through the building envelope of buildings located in Brazil, where there is a large potential of energy reduction. Simulations have been efficiently carried out for different wall and roof configurations, showing that the innovative method efficiently provides a gain of 25% on the computer run time.

preprint2019arXiv

Dispersive and non-dispersive nonlinear long wave transformations: Numerical and experimental results

The description of gravity waves propagating on the water surface is considered from a historical point of view, with specific emphasis on the development of a theoretical framework and equations of motion for long waves in shallow water. This provides the foundation for a subsequent discussion about tsunami wave propagation and run-up on a sloping beach, and in particular the role of wave dispersion for this problem. Wave tank experiments show that wave dispersion can play a significant role for the propagation and wave transformation of wave signals that include some higher frequency components. However, the maximum run-up height is less sensitive to dispersive effects, suggesting that run-up height can be adequately calculated by use of non-dispersive model equations.

preprint2019arXiv

Dispersive effects during long wave run-up on a plane beach

Dispersive effects during long wave run-up on a plane beach are studied. We take an advantage of experimental data collection of different wave types (single pulses, sinusoidal waves, bi-harmonic waves, and frequency modulated wave trains) and simulate their run-up using two models: (i) non-dispersive nonlinear shallow water theory and (ii) dispersive Boussinesq type model based on the modified Peregrine system. It is shown, that for long positive pulses, dispersive effects are not so important and nonlinear shallow water theory can be used. However, for periodic sinusoidal and bi-harmonic pulses of the same period, the dispersive effects result in significant wave transformation during its propagation, but do not have a strong impact on its maximal run-up height. Overall, for maximum wave run-up height, we could not find a preference of dispersive model against the nondispersive one, and, therefore, suggest using nonlinear shallow water model for long wave run-up height estimation.

preprint2019arXiv

Effects of vorticity on the travelling waves of some shallow water two-component systems

In the present study we consider three two-component (integrable and non-integrable) systems which describe the propagation of shallow water waves on a constant shear current. Namely, we consider the two-component Camassa-Holm equations, the Zakharov-Ito system and the Kaup--Boussinesq equations all including constant vorticity effects. We analyze both solitary and periodic-type travelling waves using the simple and geometrically intuitive phase space analysis. We get the pulse-type solitary wave solutions and the front solitary wave solutions. For the Zakharov-Ito system we underline the occurrence of the pulse and anti-pulse solutions. The front wave solutions decay algebraically in the far field. For the Kaup-Boussinesq system, interesting analytical multi-pulsed travelling wave solutions are found.

preprint2019arXiv

Evaluation of the reliability of building energy performance models for parameter estimation

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.

preprint2019arXiv

Experimental and numerical study of the propagation of focused wave groups in the nearshore zone

The propagation of focused wave groups in intermediate water depth and the shoaling zone is experimentally and numerically considered in this paper. The experiments are carried out in a two-dimensional wave flume and wave trains derived from Pierson-Moskowitz and JONSWAP spectrum are generated. The peak frequency does not change during the wave train propagation for Pierson-Moskowitz waves; however, a downshift of this peak is observed for JONSWAP waves. An energy partitioning is performed in order to track the spatial evolution of energy. Four energy regions are defined for each spectrum type. A nonlinear energy transfer between different spectral regions as the wave train propagates is demonstrated and quantified. Numerical simulations are conducted using a modified Boussinesq model for long waves in shallow waters of varying depth. Experimental results are in satisfactory agreement with numerical predictions, especially in the case of wave trains derived from JONSWAP spectrum.

preprint2019arXiv

Hamiltonian regularisation of shallow water equations with uneven bottom

The regularisation of nonlinear hyperbolic conservation laws has been a problem of great importance for achieving uniqueness of weak solutions and also for accurate numerical simulations. In a recent work, the first two authors proposed a so-called Hamiltonian regularisation for nonlinear shallow water and isentropic Euler equations. The characteristic property of this method is that the regularisation of solutions is achieved without adding any artificial dissipation or ispersion. The regularised system possesses a Hamiltonian structure and, thus, formally preserves the corresponding energy functional. In the present article we generalise this approach to shallow water waves over general, possibly time-dependent, bottoms. The proposed system is solved numerically with continuous Galerkin method and its solutions are compared with the analogous solutions of the classical shallow water and dispersive Serre-Green-Naghdi equations. The numerical results confirm the absence of dispersive and dissipative effects in presence of bathymetry variations.

preprint2019arXiv

Interpreting Distortions in Dimensionality Reduction by Superimposing Neighbourhood Graphs

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.

preprint2019arXiv

Long wave interaction with a partially immersed body. Part I: Mathematical models

In the present article we consider the problem of wave interaction with a partially immersed, but floating body. We assume that the motion of the body is prescribed. The general mathematical formulation for this problem is presented in the framework of a hierarchy of mathematical models. Namely, in this first part we formulate the problem at every hierarchical level. The special attention is payed to fully nonlinear and weakly dispersive models since they are most likely to be used in practice. For this model we have to consider separately the inner (under the body) and outer domains. Various approached to the gluing of solutions at the boundary is discussed as well. We propose several strategies which ensure the global conservation or continuity of some important physical quantities.

preprint2019arXiv

Nonlinear deformation and run-up of single tsunami waves of positive polarity: numerical simulations and analytical predictions

The estimate of individual wave run-up is especially important for tsunami warning and risk assessment as it allows to evaluate the inundation area. Here as a model of tsunami we use the long single wave of positive polarity. The period of such wave is rather long which makes it different from the famous Korteweg-de Vries soliton. This wave nonlinearly deforms during its propagation in the ocean, what results in a steep wave front formation. Situations, when waves approach the coast with a steep front are often observed during large tsunamis, e.g. 2004 Indian Ocean and 2011 15 Tohoku tsunamis. Here we study the nonlinear deformation and run-up of long single waves of positive polarity in the conjoined water basin, which consists of the constant depth section and a plane beach. The work is performed numerically and analytically in the framework of the nonlinear shallow water theory. Analytically, wave propagation along the constant depth section and its run-up on a beach are considered independently without taking into account wave interaction with the toe of the bottom slope. The propagation along the bottom of constant depth is described by Riemann wave, while the wave 20 run-up on a plane beach is calculated using rigorous analytical solutions of the nonlinear shallow water theory following the Carrier-Greenspan approach. Numerically, we use the finite volume method with the second order UNO2 reconstruction in space and the third order Runge-Kutta scheme with locally adaptive time steps. During wave propagation along the constant depth section, the wave becomes asymmetric with a steep wave front. Shown, that the maximum run-up height depends on the front steepness of the incoming wave approaching the toe of the bottom slope. The corresponding formula for maximum 25 run-up height, which takes into account the wave front steepness, is proposed.

preprint2019arXiv

On some model equations for pulsatile flow in viscoelastic vessels

Considered here is the derivation of partial differential equations arising in pulsatile flow in pipes with viscoelastic walls. The equations are asymptotic models describing the propagation of long-crested pulses in pipes with cylindrical symmetry. Additional effects due to viscous stresses in bio-fluids are also taken into account. The effects of viscoelasticity of the vessels on the propagation of solitary and periodic waves in a vessel of constant radius are being explored numerically.

preprint2019arXiv

On the comparison of three numerical methods applied to building simulation: finite-differences, RC circuit approximation and a spectral method

Predictions of physical phenomena in buildings are carried out by using physical models formulated as a mathematical problem and solved by means of numerical methods, aiming at evaluating, for instance, the building thermal or hygrothermal performance by calculating distributions and fluxes of heat and moisture transfer. Therefore, the choice of the numerical method is crucial since it is a compromise among (i) the solution accuracy, (ii) the computational cost to obtain the solution and (iii) the complexity of the method implementation. An efficient numerical method enables to compute an accurate solution with a minimum computational run time (CPU). On that account, this article brings an investigation on the performance of three numerical methods. The first one is the standard and widely used finite-difference approach, while the second one is the so-called RC approach, which is a particular method brought to the building physics area by means of an analogy of electric circuits. The third numerical method is the spectral one, which has been recently proposed to solve nonlinear diffusive problems in building physics. The three methods are evaluated in terms of accuracy on the assessment of the dependent variable (temperature or vapor pressure) or of density of fluxes for three different cases: i) heat diffusion through a concrete slab, ii) moisture diffusion through an aerated concrete slab and iii) heat diffusion using measured temperatures as boundary conditions. Results highlight the spectral approach as the most accurate method. The RC based model with a few number of resistances does not provide accurate results for temperature and vapor pressure distributions neither to flux densities nor conduction loads.

preprint2019arXiv

On the multi-symplectic structure of Boussinesq-type systems. II: Geometric discretization

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.

preprint2019arXiv

On the Reducibility and the Lenticular Sets of Zeroes of Almost Newman Lacunary Polynomials

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).

preprint2019arXiv

On time relaxed schemes and formulations for dispersive wave equations

The numerical simulation of nonlinear dispersive waves is a central research topic of many investigations in the nonlinear wave community. Simple and robust solvers are needed for numerical studies of water waves as well. The main difficulties arise in the numerical approximation of high order derivatives and in severe stability restrictions on the time step, when explicit schemes are used. In this study we propose new relaxed system formulations which approximate the initial dispersive wave equation. However, the resulting relaxed system involves first order derivatives only and it is written in the form of an evolution problem. Thus, many standard methods can be applied to solve the relaxed problem numerically. In this article we illustrate the application of the new relaxed scheme on the classical Korteweg-de Vries equation as a prototype of stiff dispersive PDEs.

preprint2018arXiv

Accurate fast computation of steady two-dimensional surface gravity waves in arbitrary depth

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.

preprint2018arXiv

Accurate numerical simulation of moisture front in porous material

When comparing measurements to numerical simulations of moisture transfer through porous materials a rush of the experimental moisture front is commonly observed in several works shown in the literature, with transient models that consider only the diffusion process. Thus, to overcome the discrepancies between the experimental and the numerical models, this paper proposes to include the moisture advection transfer in the governing equation. To solve the advection-diffusion differential equation, it is first proposed two efficient numerical schemes and their efficiencies are investigated for both linear and nonlinear cases. The first scheme, Scharfetter-Gummel (SG), presents a Courant-Friedrichs-Lewy (CFL) condition but is more accurate and faster than the second scheme, the well-known Crank-Nicolson approach. Furthermore, the SG scheme has the advantages of being well-balanced and asymptotically preserved. Then, to conclude, results of the convective moisture transfer problem obtained with the SG numerical scheme are compared to experimental data from the literature. The inclusion of an advective term in the model may clearly lead to better results than purely diffusive models.

preprint2018arXiv

Advanced reduced-order models for moisture diffusion in porous media

It is of great concern to produce numerically efficient methods for moisture diffusion through porous media, capable of accurately calculate moisture distribution with a reduced computational effort. In this way, model reduction methods are promising approaches to bring a solution to this issue since they do not degrade the physical model and provide a significant reduction of computational cost. Therefore, this article explores in details the capabilities of two model-reduction techniques - the Spectral Reduced-Order Model (Spectral-ROM) and the Proper Generalised Decomposition (PGD) - to numerically solve moisture diffusive transfer through porous materials. Both approaches are applied to three different problems to provide clear examples of the construction and use of these reduced-order models. The methodology of both approaches is explained extensively so that the article can be used as a numerical benchmark by anyone interested in building a reduced-order model for diffusion problems in porous materials. Linear and non-linear unsteady behaviors of unidimensional moisture diffusion are investigated. The last case focuses on solving a parametric problem in which the solution depends on space, time and the diffusivity properties. Results have highlighted that both methods provide accurate solutions and enable to reduce significantly the order of the model around ten times lower than the large original model. It also allows an efficient computation of the physical phenomena with an error lower than 10^{-2} when compared to a reference solution.

preprint2018arXiv

An efficient method to estimate sorption isotherm curve coefficients

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.

preprint2018arXiv

An improved explicit scheme for whole-building hygrothermal simulation

Implicit schemes require important sub-iterations when dealing with highly nonlinear problems such as the combined heat and moisture transfer through porous building elements. The computational cost rises significantly when the whole-building is simulated, especially when there is important coupling among the building elements themselves with neighbouring zones and with HVAC (Heating Ventilation and Air Conditioning) systems. On the other hand, the classical Euler explicit scheme is generally not used because its stability condition imposes very fine time discretisation. Hence, this paper explores the use of an improved explicit approach - the Dufort-Frankel scheme - to overcome the disadvantage of the classical explicit one and to bring benefits that cannot be obtained by implicit methods. The Dufort-Frankel approach is first compared to the classical Euler implicit and explicit schemes to compute the solution of nonlinear heat and moisture transfer through porous materials. Then, the analysis of the Dufort-Frankel unconditionally stable explicit scheme is extended to the coupled heat and moisture balances on the scale of a one- and a two-zone building models. The Dufort-Frankel scheme has the benefits of being unconditionally stable, second-order accurate in time O(dt^2) and to compute explicitly the solution at each time step, avoiding costly sub-iterations. This approach may reduce the computational cost by twenty, as well as it may enable perfect synchronism for whole-building simulation and co-simulation.

preprint2018arXiv

Analysis and improvement of the VTT mold growth model: application to bamboo fiberboard

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.

preprint2018arXiv

Asymptotic nonlinear and dispersive pulsatile flow in elastic vessels with cylindrical symmetry

The asymptotic derivation of a new family of one-dimensional, weakly nonlinear and weakly dispersive equations that model the flow of an ideal fluid in an elastic vessel is presented. Dissipative effects due to the viscous nature of the fluid are also taken into account. The new models validate by asymptotic reasoning other non-dispersive systems of equations that are commonly used, and improve other nonlinear and dispersive mathematical models derived to describe the blood flow in elastic vessels. The new systems are studied analytically in terms of their basic characteristic properties such as the linear dispersion characteristics, symmetries, conservation laws and solitary waves. Unidirectional model equations are also derived and analysed in the case of vessels of constant radius. The capacity of the models to be used in practical problems is being demonstrated by employing a particular system with favourable properties to study the blood flow in a large artery. Two different cases are considered: A vessel with constant radius and a tapered vessel. Significant changes in the flow can be observed in the case of the tapered vessel.

preprint2018arXiv

Direct dynamical energy cascade in the modified KdV equation

In this study we examine the energy transfer mechanism during the nonlinear stage of the Modulational Instability (MI) in the modified Korteweg-de Vries equation. The particularity of this study consists in considering the problem essentially in the Fourier space. A dynamical energy cascade model of this process originally proposed for the focusing NLS-type equations is transposed to the mKdV setting using the existing connections between the KdV-type and NLS-type equations. The main predictions of the D-cascade model are outlined and thoroughly discussed. Finally, the obtained theoretical results are validated by direct numerical simulations of the mKdV equation using the pseudo-spectral methods. A general good agreement is reported in this study. The nonlinear stages of the MI evolution are also investigated for the mKdV equation.

preprint2018arXiv

Dispersive shallow water wave modelling. Part III: Model derivation on a globally spherical geometry

The present article is the third part of a series of papers devoted to the shallow water wave modelling. In this part, we investigate the derivation of some long wave models on a deformed sphere. We propose first a suitable for our purposes formulation of the full Euler equations on a sphere. Then, by applying the depth-averaging procedure we derive first a new fully nonlinear weakly dispersive base model. After this step, we show how to obtain some weakly nonlinear models on the sphere in the so-called Boussinesq regime. We have to say that the proposed base model contains an additional velocity variable which has to be specified by a closure relation. Physically, it represents a dispersive correction to the velocity vector. So, the main outcome of our article should be rather considered as a whole family of long wave models.

preprint2018arXiv

Dispersive shallow water wave modelling. Part IV: Numerical simulation on a globally spherical geometry

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.

preprint2018arXiv

Evaluation of tsunami wave energy generated by earthquakes in the Makran subduction zone

The Makran subduction zone, an approximate 1000 km section of the Eurasian-Arabian plate, is located offshore of Southern Iran and Pakistan. In 1945, the Makran subduction zone (MSZ) generated a tsunamigenic earthquake with a magnitude of Mw 8.1. The region has also experienced large historical earthquakes but the data regarding these events are poorly documented. Therefore, the need to investigate tsunamis in Makran must be taken into serious consideration. Using hydrodynamic numerical simulation, we evaluate the tsunami wave energy generated by bottom motion for a tsunamigenic source model distributed along the full length of the Makran subduction zone. The whole rupture of the plate boundary is divided into 20 segments with width of order of 200 km and a co-seismic slip of 10 m but with various lengths. Exchanges between kinetic and potential components of tsunami wave energy are shown. The total tsunami wave energy displays only 0.33% of the seismic energy released from the earthquake source. As a result, for every increase in magnitude by one unit, the associated tsunami wave energy becomes about 10^3 times greater.

preprint2018arXiv

Numerical modelling of surface water wave interaction with a moving wall

In the present manuscript, we consider the practical problem of wave interaction with a vertical wall. However, the novelty here consists in the fact that the wall can move horizontally due to a system of springs. The water wave evolution is described with the free surface potential flow model. Then, a semi-analytical numerical method is presented. It is based on a mapping technique and a finite difference scheme in the transformed domain. The idea is to pose the equations on a fixed domain. This method is thoroughly tested and validated in our study. By choosing specific values of spring parameters, this system can be used to damp (or in other words to extract the energy of) incident water waves.

preprint2018arXiv

On the modelling of shallow turbidity flows

In this study we investigate shallow turbidity density currents and underflows from mechanical point of view. We propose a simple hyperbolic model for such flows. On one hand, our model is based on very basic conservation principles. On the other hand, the turbulent nature of the flow is also taken into account through the energy dissipation mechanism. Moreover, the mixing with the pure water along with sediments entrainment and deposition processes are considered, which makes the problem dynamically interesting. One of the main advantages of our model is that it requires the specification of only two modeling parameters - the rate of turbulent dissipation and the rate of the pure water entrainment. Consequently, the resulting model turns out to be very simple and self-consistent. This model is validated against several experimental data and several special classes of solutions (such as travelling, self-similar and steady) are constructed. Unsteady simulations show that some special solutions are realized as asymptotic long time states of dynamic trajectories.

preprint2018arXiv

On the multi-symplectic structure of Boussinesq-type systems. I: Derivation and mathematical properties

The Boussinesq equations are known since the end of the XIXst century. However, the proliferation of various \textsc{Boussinesq}-type systems started only in the second half of the XXst century. Today they come under various flavours depending on the goals of the modeller. At the beginning of the XXIst century an effort to classify such systems, at least for even bottoms, was undertaken and developed according to both different physical regimes and mathematical properties, with special emphasis, in this last sense, on the existence of symmetry groups and their connection to conserved quantities. Of particular interest are those systems admitting a symplectic structure, with the subsequent preservation of the total energy represented by the Hamiltonian. In the present paper a family of Boussinesq-type systems with multi-symplectic structure is introduced. Some properties of the new systems are analyzed: their relation with already known Boussinesq models, the identification of those systems with additional Hamiltonian structure as well as other mathematical features like well-posedness and existence of different types of solitary-wave solutions. The consistency of multi-symplectic systems with the full Euler equations is also discussed.

preprint2018arXiv

On the solution of coupled heat and moisture transport in porous material

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.

preprint2018arXiv

On the velocity of turbidity currents over moderate slopes

In the present article we consider the problem of underwater avalanches propagating over moderate slopes. The main goal of our work is to investigate the avalanche front velocity selection mechanism when it propagates downwards. In particular, we show that the front velocity does not depend univocally on the mass of sediments. This phenomenon is investigated and explained in our study. Moreover, we derive from the first principles a depth-averaged model. Then, we assume that sediments are uniformly distributed along the slope. In this case, they can be entrained into the flow head and a self-sustained regime can be established. One of the main findings of our study is that the avalanche front velocity is not unique due to a hysteresis phenomenon. We attempt to explain this phenomenon using dynamical systems considerations.

preprint2018arXiv

Peregrine's system revisited

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.

preprint2018arXiv

Solitary wave solutions and their interactions for fully nonlinear water waves with surface tension in the generalized Serre equations

Some effects of surface tension on fully-nonlinear, long, surface water waves are studied by numerical means. The differences between various solitary waves and their interactions in subcritical and supercritical surface tension regimes are presented. Analytical expressions for new peaked travelling wave solutions are presented in the case of critical surface tension. The numerical experiments were performed using a high-accurate finite element method based on smooth cubic splines and the four-stage, classical, explicit Runge-Kutta method of order four.

preprint2018arXiv

Solving nonlinear diffusive problems in buildings by means of a Spectral Reduced-Order Model

This paper proposes the use of a Spectral method to simulate diffusive moisture transfer through porous materials as a Reduced-Order Model (ROM). The Spectral approach is an a priori method assuming a separated representation of the solution. The method is compared with both classical Euler implicit and Crank-Nicolson schemes, considered as large original models. Their performance - in terms of accuracy, complexity reduction and CPU time reduction - are discussed for linear and nonlinear cases of moisture diffusive transfer through single and multi-layered one-dimensional domains, considering highly moisture-dependent properties. Results show that the Spectral reduced-order model approach enables to simulate accurately the field of interest. Furthermore, numerical gains become particularly interesting for nonlinear cases since the proposed method can drastically reduce the computer run time, by a factor of 100, when compared to the traditional Crank-Nicolson scheme for one-dimensional applications.

preprint2018arXiv

Tsunami hazard assessment in the Makran subduction zone

The lack of offshore seismic data caused uncertainties associated to understating the behavior of future tsunamigenic earthquakes in the Makran Subduction Zone (MSZ). Future tsunamigenic events in the MSZ may trigger significant near-field tsunamis. Run-up heights in the near-field are controlled by the heterogeneity of slip over the rupture area. Considering a non-planar geometry for the Makran subduction zone, a range of random k^{-2} slip models were generated to hypothesize rupturing on the fault zone. We model tsunamis numerically and assess probabilistic tsunami hazard in the near-field for all synthetic scenarios. The main affected areas by tsunami waves are the area between Jask and Ormara along the shorelines of Iran and Pakistan and the area between Muscat and Sur along the Oman coastline. The maximum peak-run-up along the shores of Iran and Pakistan is about 16 m and about 12 m for the Oman shoreline. The slip distributions control the run-up along the Makran coastlines. The dependency of run-up to the heterogeneity of slip is higher in the most impacted areas. Those areas are more vulnerable to tsunami hazard than other areas.

preprint2018arXiv

Wave dynamics on networks: method and application to the sine-Gordon equation

We consider a scalar Hamiltonian nonlinear wave equation formulated on networks; this is a non standard problem because these domains are not locally homeomorphic to any subset of the Euclidean space. More precisely, we assume each edge to be a 1D uniform line with end points identified with graph vertices. The interface conditions at these vertices are introduced and justified using conservation laws and an homothetic argument. We present a detailed methodology based on a symplectic finite difference scheme together with a special treatment at the junctions to solve the problem and apply it to the sine-Gordon equation. Numerical results on a simple graph containing four loops show the performance of the scheme for kinks and breathers initial conditions.

preprint2018arXiv

Weakly singular shock profiles for a non-dispersive regularization of shallow-water equations

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.

preprint2017arXiv

New asymptotic heat transfer model in thin liquid films

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.

preprint2017arXiv

Non-dispersive conservative regularisation of nonlinear shallow water (and isothermal Euler) equations

A new regularisation of the shallow water (and isentropic Euler) equations is proposed. The regularised equations are non-dissipative, non-dispersive and possess a variational structure. Thus, the mass, the momentum and the energy are conserved. Hence, for instance, regularised hydraulic jumps are smooth and non-oscillatory. Another particularly interesting feature of this regularisation is that smoothed `shocks' propagates at exactly the same speed as the original discontinuous ones. The performance of the new model is illustrated numerically on some dam-break test cases, which are classical in the hyperbolic realm.

preprint2017arXiv

On supraconvergence phenomenon for second order centered finite differences on non-uniform grids

In the present study we consider an example of a boundary value problem for a simple second order ordinary differential equation, which may exhibit a boundary layer phenomenon. We show that usual central finite differences, which are second order accurate on a uniform grid, can be substantially upgraded to the fourth order by a suitable choice of the underlying non-uniform grid. This example is quite pedagogical and may give some ideas for more complex problems.

preprint2017arXiv

On the estimation of moisture permeability and advection coefficients of a wood fibre material using the optimal experiment design approach

This paper presents a practical application of the concept of Optimal Experiment Design (OED) for the determination of properties of porous materials with in-situ measurements and an identification method. First, an experimental set-up was presented and used for the measurement of relative humidity within a wood fibre material submitted to single and multiple steps of relative humidity variation. Then, the application of OED enabled to plan the experimental conditions in terms of sensor positioning and boundary conditions out of 20 possible designs. The OED search was performed using the Fisher information matrix and a priori knowledge of the parameters. It ensures to provide the best accuracy of the identification method and thus the estimated parameter. Optimal design results have been found for single steps from the relative humidity phi = 10 to 75%, with one sensor located at the position X between 4 and 6 cm, for the estimation of moisture permeability coefficients, while from phi = 75% to phi = 33%, with one sensor located at X° = 3 cm, for the estimation of the advection coefficient. The OED has also been applied for the identification of couples of parameters. A sample submitted to multiple relative humidity steps (phi = 10-75-33-75%) with a sensor placed at X° = 5 cm was found as the best option for determining both properties with the same experiment. These OED parameters have then been used for the determination of moisture permeability and advection coefficients. The estimated moisture permeability coefficients are twice higher than the a priori values obtained using standard methods. The advection parameter corresponds to the mass average velocity of the order of v = 0.01 mm/s within the material and may play an important role on the simulation of moisture front.

preprint2017arXiv

On weakly singular and fully nonlinear travelling shallow capillary-gravity waves in the critical regime

In this Letter we consider long capillary-gravity waves described by a fully nonlinear weakly dispersive model. First, using the phase space analysis methods we describe all possible types of localized travelling waves. Then, we especially focus on the critical regime, where the surface tension is exactly balanced by the gravity force. We show that our long wave model with a critical Bond number admits stable travelling wave solutions with a singular crest. These solutions are usually referred to in the literature as peakons or peaked solitary waves. They satisfy the usual speed-amplitude relation, which coincides with Scott-Russel's empirical formula for solitary waves, while their decay rate is the same regardless their amplitude. Moreover, they can be of depression or elevation type independent of their speed. The dynamics of these solutions are studied as well.

preprint2017arXiv

Some special solutions to the Hyperbolic NLS equation

The Hyperbolic Nonlinear Schrodinger equation (HypNLS) arises as a model for the dynamics of three-dimensional narrowband deep water gravity waves. In this study, the Petviashvili method is exploited to numerically compute bi-periodic time-harmonic solutions of the HypNLS equation. In physical space they represent non-localized standing waves. Non-trivial spatial patterns are revealed and an attempt is made to describe them using symbolic dynamics and the language of substitutions. Finally, the dynamics of a slightly perturbed standing wave is numerically investigated by means a highly acccurate Fourier solver.

preprint2017arXiv

Stable explicit schemes for simulation of nonlinear moisture transfer in porous materials

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.

preprint2016arXiv

Algebraic method for constructing singular steady solitary waves: A case study

This article describes the use of algebraic methods in a phase plane analysis of ordinary differential equations. The method is illustrated by the study of capillary-gravity steady surface waves propagating in shallow water. We consider the (fully nonlinear, weakly dispersive) Serre-Green-Naghdi equations with surface tension, because it provides a tractable model that, in the same time, is not too simple so the interest of the method can be emphasised. In particular, we analyse a special class of solutions, the solitary waves, which play an important role in many fields of Physics. In capillary-gravity regime, there are two kinds of localised infinitely smooth travelling wave solutions -- solitary waves of elevation and of depression. However, if we allow the solitary waves to have an angular point, the "zoology" of solutions becomes much richer and the main goal of this study is to provide a complete classification of such singular localised solutions using the methods of the effective Algebraic Geometry.

preprint2016arXiv

Conservative modified Serre-Green-Naghdi equations with improved dispersion characteristics

For surface gravity waves propagating in shallow water, we propose a variant of the fully nonlinear Serre-Green-Naghdi equations involving a free parameter that can be chosen to improve the dispersion properties. The novelty here consists in the fact that the new model conserves the energy, contrary to other modified Serre's equations found in the literature. Numerical comparisons with the Euler equations show that the new model is substantially more accurate than the classical Serre equations, specially for long time simulations and for large amplitudes.

preprint2016arXiv

Efficient computation of capillary-gravity generalized solitary waves

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.

preprint2016arXiv

Macroscopic dynamics of incoherent soliton ensembles: soliton-gas kinetics and direct numerical modeling

We undertake a detailed comparison of the results of direct numerical simulations of the integrable soliton gas dynamics with the analytical predictions inferred from the exact solutions of the relevant kinetic equation for solitons. We use the KdV soliton gas as a simplest analytically accessible model yielding major insight into the general properties of soliton gases in integrable systems. Two model problems are considered: (i) the propagation of a `trial' soliton through a one-component `cold' soliton gas consisting of randomly distributed solitons of approximately the same amplitude; and (ii) collision of two cold soliton gases of different amplitudes (soliton gas shock tube problem) leading to the formation of an incoherend dispersive shock wave. In both cases excellent agreement is observed between the analytical predictions of the soliton gas kinetics and the direct numerical simulations. Our results confirm relevance of the kinetic equation for solitons as a quantitatively accurate model for macroscopic non-equilibrium dynamics of incoherent soliton ensembles.

preprint2016arXiv

Modified Shallow Water Equations for significantly varying seabeds

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.

preprint2016arXiv

Multi-symplectic structure of fully-nonlinear weakly-dispersive internal gravity waves

In this short communication we present the multi-symplectic structure for the two-layer Serre-Green-Naghdi equations describing the evolution of large amplitude internal gravity long waves. We consider only a two-layer stratification with rigid bottom and lid for simplicity, generalisations to several layers being straightforward. This multi-symplectic formulation allows the application of various multi-symplectic integrators (such as Euler or Preissman box schemes) that preserve exactly the multi-symplecticity at the discrete level.

preprint2016arXiv

On the nonlinear dynamics of the traveling-wave solutions of the Serre system

We numerically study nonlinear phenomena related to the dynamics of traveling wave solutions of the Serre equations including the stability, the persistence, the interactions and the breaking of solitary waves. The numerical method utilizes a high-order finite-element method with smooth, periodic splines in space and explicit Runge-Kutta methods in time. Other forms of solutions such as cnoidal waves and dispersive shock waves are also considered. The differences between solutions of the Serre equations and the Euler equations are also studied.

preprint2016arXiv

On the optimal experimental design for heat and moisture parameter estimation

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.

preprint2016arXiv

Serre-type equations in deep water

This manuscript is devoted to the modelling of water waves in the deep water regime with some emphasis on the underlying variational structures. The present article should be considered as a review of some existing models and modelling approaches even if new results are presented as well. Namely, we derive the deep water analogue of the celebrated Serre-Green-Naghdi equations which have become the standard model in shallow water environments. The relation to existing models is discussed. Moreover, the multi-symplectic structure of these equations is reported as well. The results of this work can be used to develop various types of robust structure-preserving variational integrators in deep water. The methodology of constructing approximate models presented in this study can be naturally extrapolated to other physical flow regimes as well.

preprint2015arXiv

A new run-up algorithm based on local high-order analytic expansions

The practically important problem of the wave run-up is studied in this article in the framework of Nonlinear Shallow Water Equations (NSWE). The main novelty consists in the usage of high order local asymptotic analytical solutions in the vicinity of the shoreline. Namely, we use the analytical techniques introduced by S. Kovalevskaya and the analogy with the compressible gas dynamics (i.e. gas outflow problem into the vacuum). Our run-up algorithm covers all the possible cases of the wave slope on the shoreline and it incorporates the new analytical information in order to determine the shoreline motion to higher accuracy. The application of this algorithm is illustrated on several important practical examples. Finally, the simulation results are compared with the well-known analytical and experimental predictions.

preprint2015arXiv

A plethora of generalised solitary gravity-capillary water waves

The present study describes, first, an efficient algorithm for computing capillary-gravity solitary waves solutions of the irrotational Euler equations with a free surface and, second, provides numerical evidences of the existence of an infinite number of generalised solitary waves (solitary waves with undamped oscillatory wings). Using conformal mapping, the unknown fluid domain, which is to be determined, is mapped into a uniform strip of the complex plane. In the transformed domain, a Babenko-like equation is then derived and solved numerically.

preprint2015arXiv

Adaptive modeling of shallow fully nonlinear gravity waves

This paper presents an extended version of the celebrated Serre-Green-Naghdi (SGN) system. This extension is based on the well-known Bona-Smith-Nwogu trick which aims to improve the linear dispersion properties. We show that in the fully nonlinear setting it results in modifying the vertical acceleration. Even if this technique is well-known, the effect of this modification on the nonlinear properties of the model is not clear. The first goal of this study is to shed some light on the properties of solitary waves, as the most important class of nonlinear permanent solutions. Then, we propose a simple adaptive strategy to choose the optimal value of the free parameter at every instance of time. This strategy is validated by comparing the model prediction with the reference solutions of the full Euler equations and its classical counterpart. Numerical simulations show that the new adaptive model provides a much better accuracy for the same computational complexity.

preprint2015arXiv

Numerical study of the generalised Klein-Gordon equations

In this study, we discuss an approximate set of equations describing water wave propagating in deep water. These generalized Klein-Gordon (gKG) equations possess a variational formulation, as well as a canonical Hamiltonian and multi-symplectic structures. Periodic travelling wave solutions are constructed numerically to high accuracy and compared to a seventh-order Stokes expansion of the full Euler equations. Then, we propose an efficient pseudo-spectral discretisation, which allows to assess the stability of travelling waves and localised wave packets.

preprint2015arXiv

On the multi-symplectic structure of the Serre-Green-Naghdi equations

In this short note, we present a multi-symplectic structure of the Serre-Green-Naghdi (SGN) equations modelling nonlinear long surface waves in shallow water. This multi-symplectic structure allow the use of efficient finite difference or pseudo-spectral numerical schemes preserving exactly the multi-symplectic form at the discrete level.

preprint2015arXiv

Practical use of variational principles for modeling water waves

This paper describes a method for deriving approximate equations for irrotational water waves. The method is based on a 'relaxed' variational principle, i.e., on a Lagrangian involving as many variables as possible. This formulation is particularly suitable for the construction of approximate water wave models, since it allows more freedom while preserving the variational structure. The advantages of this relaxed formulation are illustrated with various examples in shallow and deep waters, as well as arbitrary depths. Using subordinate constraints (e.g., irrotationality or free surface impermeability) in various combinations, several model equations are derived, some being well-known, other being new. The models obtained are studied analytically and exact travelling wave solutions are constructed when possible.

preprint2014arXiv

Numerical simulation of a solitonic gas in KdV and KdV-BBM equations

The collective behaviour of soliton ensembles (i.e. the solitonic gas) is studied using the methods of the direct numerical simulation. Traditionally this problem was addressed in the context of integrable models such as the celebrated KdV equation. We extend this analysis to non-integrable KdV-BBM type models. Some high resolution numerical results are presented in both integrable and nonintegrable cases. Moreover, the free surface elevation probability distribution is shown to be quasi-stationary. Finally, we employ the asymptotic methods along with the Monte-Carlo simulations in order to study quantitatively the dependence of some important statistical characteristics (such as the kurtosis and skewness) on the Stokes-Ursell number (which measures the relative importance of nonlinear effects compared to the dispersion) and also on the magnitude of the BBM term.

preprint2014arXiv

The Whitham Equation as a Model for Surface Water Waves

The Whitham equation was proposed as an alternate model equation for the simplified description of uni-directional wave motion at the surface of an inviscid fluid. As the Whitham equation incorporates the full linear dispersion relation of the water wave problem, it is thought to provide a more faithful description of shorter waves of small amplitude than traditional long wave models such as the KdV equation. In this work, we identify a scaling regime in which the Whitham equation can be derived from the Hamiltonian theory of surface water waves. The Whitham equation is integrated numerically, and it is shown that the equation gives a close approximation of inviscid free surface dynamics as described by the Euler equations. The performance of the Whitham equation as a model for free surface dynamics is also compared to two standard free surface models: the KdV and the BBM equation. It is found that in a wide parameter range of amplitudes and wavelengths, the Whitham equation performs on par with or better than both the KdV and BBM equations.

preprint2014arXiv

Visco-potential flows in electrohydrodynamics

In this study we consider the problem of the interface motion under the capillary-gravity and an external electric forces. The infinitely deep fluid layer is assumed to be viscous, perfectly conducting and the flow to be incompressible. The weak viscous effects are introduced using the Helmholtz-Leray decomposition and the visco-potential flow approach. The electric charge distributions above and on the free surface are considered. Finally, we derive some linearized analytical solutions for the free surface elevation shape under the localised pressure distribution and the combined action of the forces mentioned hereinabove.

preprint2013arXiv

Boussinesq modeling of surface waves due to underwater landslides

Consideration is given to the influence of an underwater landslide on waves at the surface of a shallow body of fluid. The equations of motion which govern the evolution of the barycenter of the landslide mass include various dissipative effects due to bottom friction, internal energy dissipation, and viscous drag. The surface waves are studied in the Boussinesq scaling, with time-dependent bathymetry. A numerical model for the Boussinesq equations is introduced which is able to handle time-dependent bottom topography, and the equations of motion for the landslide and surface waves are solved simultaneously. The numerical solver for the Boussinesq equations can also be restricted to implement a shallow-water solver, and the shallow-water and Boussinesq configurations are compared. A particular bathymetry is chosen to illustrate the general method, and it is found that the Boussinesq system predicts larger wave run-up than the shallow-water theory in the example treated in this paper. It is also found that the finite fluid domain has a significant impact on the behavior of the wave run-up.

preprint2013arXiv

Derivation of dissipative Boussinesq equations using the Dirichlet-to-Neumann operator approach

The water wave theory traditionally assumes the fluid to be perfect, thus neglecting all effects of the viscosity. However, the explanation of several experimental data sets requires the explicit inclusion of dissipative effects. In order to meet these practical problems, the theory of visco-potential flows has been developed (see P.-F. Liu & A. Orfila (2004) and D. Dutykh & F. Dias (2007)). Then, usually this formulation is further simplified by developing the potential in an entire series in the vertical coordinate and by introducing thus, the long wave approximation. In the present study we propose a derivation of dissipative Boussinesq equations which is based on asymptotic expansions of the Dirichlet-to-Neumann (D2N) operator. Both employed methods yield the same system by different ways.

preprint2013arXiv

Efficient computation of steady solitary gravity waves

An efficient numerical method to compute solitary wave solutions to the free surface Euler equations is reported. It is based on the conformal mapping technique combined with an efficient Fourier pseudo-spectral method. The resulting nonlinear equation is solved via the Petviashvili iterative scheme. The computational results are compared to some existing approaches, such as the Tanaka method and Fenton's high-order asymptotic expansion. Several important integral quantities are numerically computed for a large range of amplitudes. The integral representation of the velocity and acceleration fields in the bulk of the fluid is also provided.

preprint2013arXiv

Evolution of Random Wave Fields in the Water of Finite Depth

The evolution of random wave fields on the free surface is a complex process which is not completely understood nowadays. For the sake of simplicity in this study we will restrict our attention to the 2D physical problems only (i.e. 1D wave propagation). However, the full Euler equations are solved numerically in order to predict the wave field dynamics. We will consider the most studied deep water case along with several finite depths (from deep to shallow waters) to make a comparison. For each depth we will perform a series of Monte--Carlo runs of random initial conditions in order to deduce some statistical properties of an average sea state.

preprint2013arXiv

Fast accurate computation of the fully nonlinear solitary surface gravity waves

In this short note, we present an easy to implement and fast algorithm for the computation of the steady solitary gravity wave solution of the free surface Euler equations in irrotational motion. First, the problem is reformulated in a fixed domain using the conformal mapping technique. Second, the problem is reduced to a single equation for the free surface. Third, this equation is solved using Petviashvili's iterations together with pseudo-spectral discretisation. This method has a super-linear complexity, since the most demanding operations can be performed using a FFT algorithm. Moreover, when this algorithm is combined with the multi-precision arithmetics, the results can be obtained to any arbitrary accuracy.

preprint2013arXiv

Finite volume and pseudo-spectral schemes for the fully nonlinear 1D Serre equations

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.

preprint2013arXiv

On the Galerkin / finite-element method for the Serre equations

A highly accurate numerical scheme is presented for the Serre system of partial differential equations, which models the propagation of dispersive shallow water waves in the fully-nonlinear regime. The fully-discrete scheme utilizes the Galerkin / finite-element method based on smooth periodic splines in space, and an explicit fourth-order Runge-Kutta method in time. Computations compared with exact solitary and cnoidal wave solutions show that the scheme achieves the optimal orders of accuracy in space and time. These computations also show that the stability of this scheme does not impose restrictive conditions on the temporal step size. In addition, solitary, cnoidal, and dispersive shock waves are studied in detail using this numerical scheme for the Serre system and compared with the 'classical' Boussinesq system for small-amplitude shallow water waves. The results show that the interaction of solitary waves in the Serre system is more inelastic. The efficacy of the numerical scheme for modeling dispersive shocks is shown by comparison with asymptotic results. These results have application to the modeling of shallow water waves of intermediate or large amplitude, such as occurs in the nearshore zone.

preprint2013arXiv

On the Galilean invariance of some dispersive wave equations

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.

preprint2013arXiv

Run-up amplification of transient long waves

The extreme characteristics of long wave run-up are studied in this paper. First we give a brief overview of the existing theory which is mainly based on the hodograph transformation (Carrier & Greenspan, 1958). Then, using numerical simulations, we build on the work of Stefanakis et al. (2011) for an infinite sloping beach and we find that resonant run-up amplification of monochromatic waves is robust to spectral perturbations of the incoming wave and resonant regimes do exist for certain values of the frequency. In the setting of a finite beach attached to a constant depth region, resonance can only be observed when the incoming wavelength is larger than the distance from the undisturbed shoreline to the seaward boundary. Wavefront steepness is also found to play a role in wave run-up, with steeper waves reaching higher run-up values.

preprint2012arXiv

Finite volume methods for unidirectional dispersive wave models

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.

preprint2012arXiv

On the contribution of the horizontal sea-bed displacements into the tsunami generation process

The main reason for the generation of tsunamis is the deformation of the bottom of the ocean caused by an underwater earthquake. Usually, only the vertical bottom motion is taken into account while the horizontal co-seismic displacements are neglected in the absence of landslides. In the present study we propose a methodology based on the well-known Okada solution to reconstruct in more details all components of the bottom coseismic displacements. Then, the sea-bed motion is coupled with a three-dimensional weakly nonlinear water wave solver which allows us to simulate a tsunami wave generation. We pay special attention to the evolution of kinetic and potential energies of the resulting wave while the contribution of the horizontal displacements into wave energy balance is also quantified. Such contribution of horizontal displacements to the tsunami generation has not been discussed before, and it is different from the existing approaches. The methods proposed in this study are illustrated on the July 17, 2006 Java tsunami and some more recent events.

preprint2012arXiv

On the modelling of tsunami generation and tsunami inundation

While the propagation of tsunamis is well understood and well simulated by numerical models, there are still a number of unanswered questions related to the generation of tsunamis or the subsequent inundation. We review some of the basic generation mechanisms as well as their simulation. In particular, we present a simple and computationally inexpensive model that describes the seabed displacement during an underwater earthquake. This model is based on the finite fault solution for the slip distribution under some assumptions on the kinematics of the rupturing process. We also consider an unusual source for tsunami generation: the sinking of a cruise ship. Then we review some aspects of tsunami run-up. In particular, we explain why the first wave of a tsunami is sometimes less devastating than the subsequent waves. A resonance effect can boost the waves that come later. We also look at a particular feature of the 11 March 2011 tsunami in Japan - the formation of macro-scale vortices - and show that these macro-scale vortices can be captured by the nonlinear shallow water equations.

preprint2012arXiv

Solitary wave interaction in a compact equation for deep-water gravity waves

In this study we compute numerical traveling wave solutions to a compact version of the Zakharov equation for unidirectional deep-water waves recently derived by Dyachenko & Zakharov (2011) Furthermore, by means of an accurate Fourier-type spectral scheme we find that solitary waves appear to collide elastically, suggesting the integrability of the Zakharov equation.

preprint2012arXiv

Special solutions to a compact equation for deep-water gravity waves

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.

preprint2012arXiv

The Force of a Tsunami on a Wave Energy Converter

With an increasing emphasis on renewable energy resources, wave power technology is fast becoming a realistic solution. However, the recent tsunami in Japan was a harsh reminder of the ferocity of the ocean. It is known that tsunamis are nearly undetectable in the open ocean but as the wave approaches the shore its energy is compressed creating large destructive waves. The question posed here is whether a nearshore wave energy converter (WEC) could withstand the force of an incoming tsunami. The analytical 3D model of Renzi & Dias (2012) developed within the framework of a linear theory and applied to an array of fixed plates is used. The time derivative of the velocity potential allows the hydrodynamic force to be calculated.

preprint2011arXiv

Dispersive wave runup on non-uniform shores

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).

preprint2011arXiv

Dispersive waves generated by an underwater landslide

In this work we study the generation of water waves by an underwater sliding mass. The wave dynamics are assumed to fell into the shallow water regime. However, the characteristic wavelength of the free surface motion is generally smaller than in geophysically generated tsunamis. Thus, dispersive effects need to be taken into account. In the present study the fluid layer is modeled by the Peregrine system modified appropriately and written in conservative variables. The landslide is assumed to be a quasi-deformable body of mass whose trajectory is completely determined by its barycenter motion. A differential equation modeling the landslide motion along a curvilinear bottom is obtained by projecting all the forces acting on the submerged body onto a local moving coordinate system. One of the main novelties of our approach consists in taking into account curvature effects of the sea bed.

preprint2011arXiv

Finite volume schemes for Boussinesq type equations

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.

preprint2011arXiv

Finite volume schemes for dispersive wave propagation and runup

preprint2011arXiv

Hamiltonian form and solitary waves of the spatial Dysthe equations

The spatial Dysthe equations describe the envelope evolution of the free-surface and potential of gravity waves in deep waters. Their Hamiltonian structure and new invariants are unveiled by means of a gauge transformation to a new canonical form of the evolution equations. An accurate Fourier-type spectral scheme is used to solve for the wave dynamics and validate the new conservation laws, which are satisfied up to machine precision. Traveling waves are numerically constructed using the Petviashvili method. It is shown that their collision appears inelastic, suggesting the non-integrability of the Dysthe equations.

preprint2011arXiv

Local Runup Amplification By Resonant Wave Interactions

Until now the analysis of long wave runup on a plane beach has been focused on finding its maximum value, failing to capture the existence of resonant regimes. One-dimensional numerical simulations in the framework of the Nonlinear Shallow Water Equations (NSWE) are used to investigate the Boundary Value Problem (BVP) for plane and non-trivial beaches. Monochromatic waves, as well as virtual wave-gage recordings from real tsunami simulations, are used as forcing conditions to the BVP. Resonant phenomena between the incident wavelength and the beach slope are found to occur, which result in enhanced runup of non-leading waves. The evolution of energy reveals the existence of a quasi-periodic state for the case of sinusoidal waves, the energy level of which, as well as the time required to reach that state, depend on the incident wavelength for a given beach slope. Dispersion is found to slightly reduce the value of maximum runup, but not to change the overall picture. Runup amplification occurs for both leading elevation and depression waves.

preprint2011arXiv

Long wave runup on random beaches

The estimation of the maximum wave sunup height is a problem of practical importance. Most of the analytical and numerical studies are limited to a constant slope plain shore and to the classical Nonlinear Shallow Water (NSW) equations. However, in nature the shore is characterized by some roughness. In order to take into account the effects of the bottom rugosity various ad-hoc friction terms are usually used. In this paper we study the effect of the roughness of the bottom on the maximum runup height. A stochastic model is proposed to describe the bottom irregularity and its effect is quantified using Monte-Carlo simulations. For the discretization of the NSW equations we employ modern finite volume schemes. Moreover, the results of the random bottom model are compared with the more conventional approaches.

preprint2011arXiv

On the use of finite fault solution for tsunami generation problems

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.

preprint2011arXiv

Shallow water equations for large bathymetry variations

In this study, we propose an improved version of the nonlinear shallow water (or Saint-Venant) equations. This new model is designed to take into account the effects resulting from the large spacial and/or temporal variations of the seabed. The model is derived from a variational principle by choosing the appropriate shallow water ansatz and imposing suitable constraints. Thus, the derivation procedure does not explicitly involve any small parameter.

preprint2011arXiv

The VOLNA code for the numerical modelling of tsunami waves: generation, propagation and inundation

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.

preprint2010arXiv

Mathematical modeling of powder-snow avalanche flows

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.

preprint2010arXiv

On the relevance of the dam break problem in the context of nonlinear shallow water equations

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.

preprint2010arXiv

Velocity and energy relaxation in two-phase flows

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.

preprint2009arXiv

A compressible two-fluid model for the finite volume simulation of violent aerated flows. Analytical properties and numerical results

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.

preprint2009arXiv

A two-fluid model for violent aerated flows

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. 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 and analyze some of its properties. It is shown that this model can successfully mimic water wave impacts on coastal structures. The governing equations are discretized by a second-order finite volume method. Numerical results are presented for two examples: the dam break problem and the drop test problem. It is shown that this basic model can be used to study violent aerated flows, especially by providing fast qualitative estimates.

preprint2009arXiv

Group and phase velocities in the free-surface visco-potential flow: new kind of boundary layer induced instability

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.

preprint2009arXiv

Influence of sedimentary layering on tsunami generation

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.

preprint2009arXiv

Two-fluid barotropic models for powder-snow avalanche flows

In the present study we discuss several modeling issues of powder-snow avalanche flows. We take a two-fluid modeling paradigm. For the sake of simplicity, we will restrict our attention to barotropic equations. We begin the exposition by a compressible model with two velocities for each fluid. However, this model may become non-hyperbolic and thus, represents serious challenges for numerical methods. To overcome these issues, we derive a single velocity model as a result of a relaxation process. This model can be easily shown to be hyperbolic for any reasonable equation of state. Finally, an incompressible limit of this model is derived.

preprint2008arXiv

Simulation of Free Surface Compressible Flows Via a Two Fluid Model

The purpose of this communication is to discuss the simulation of a free surface compressible flow between two fluids, typically air and water. We use a two fluid model with the same velocity, pressure and temperature for both phases. In such a numerical model, the free surface becomes a thin three dimensional zone. The present method has at least three advantages: (i) the free-surface treatment is completely implicit; (ii) it can naturally handle wave breaking and other topological changes in the flow; (iii) one can easily vary the Equation of States (EOS) of each fluid (in principle, one can even consider tabulated EOS). Moreover, our model is unconditionally hyperbolic for reasonable EOS.

preprint2008arXiv

Tsunami generation by dynamic displacement of sea bed due to dip-slip faulting

In classical tsunami-generation techniques, one neglects the dynamic sea bed displacement resulting from fracturing of a seismic fault. The present study takes into account these dynamic effects. Earth's crust is assumed to be a Kelvin-Voigt material. The seismic source is assumed to be a dislocation in a viscoelastic medium. The fluid motion is described by the classical nonlinear shallow water equations (NSWE) with time-dependent bathymetry. The viscoelastodynamic equations are solved by a finite-element method and the NSWE by a finite-volume scheme. A comparison between static and dynamic tsunami-generation approaches is performed. The results of the numerical computations show differences between the two approaches and the dynamic effects could explain the complicated shapes of tsunami wave trains.

preprint2008arXiv

Visco-potential free-surface flows and long wave modelling

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.

preprint2007arXiv

Comparison between three-dimensional linear and nonlinear tsunami generation models

The modeling of tsunami generation is an essential phase in understanding tsunamis. For tsunamis generated by underwater earthquakes, it involves the modeling of the sea bottom motion as well as the resulting motion of the water above it. A comparison between various models for three-dimensional water motion, ranging from linear theory to fully nonlinear theory, is performed. It is found that for most events the linear theory is sufficient. However, in some cases, more sophisticated theories are needed. Moreover, it is shown that the passive approach in which the seafloor deformation is simply translated to the ocean surface is not always equivalent to the active approach in which the bottom motion is taken into account, even if the deformation is supposed to be instantaneous.

preprint2007arXiv

Dissipative Boussinesq equations

The classical theory of water waves is based on the theory of inviscid flows. However it is important to include viscous effects in some applications. Two models are proposed to add dissipative effects in the context of the Boussinesq equations, which include the effects of weak dispersion and nonlinearity in a shallow water framework. The dissipative Boussinesq equations are then integrated numerically.

preprint2007arXiv

Viscous potential free-surface flows in a fluid layer of finite depth

It is shown how to model weakly dissipative free-surface flows using the classical potential flow approach. The Helmholtz-Leray decomposition is applied to the linearized 3D Navier-Stokes equations. The governing equations are treated using Fourier--Laplace transforms. We show how to express the vortical component of the velocity only in terms of the potential and free-surface elevation. A new predominant nonlocal viscous term is derived in the bottom kinematic boundary condition. The resulting formulation is simple and does not involve any correction procedure as in previous viscous potential flow theories [Joseph2004]. Corresponding long wave model equations are derived.

preprint2007arXiv

Water waves generated by a moving bottom

Tsunamis are often generated by a moving sea bottom. This paper deals with the case where the tsunami source is an earthquake. The linearized water-wave equations are solved analytically for various sea bottom motions. Numerical results based on the analytical solutions are shown for the free-surface profiles, the horizontal and vertical velocities as well as the bottom pressure.

preprint2006arXiv

Dynamics of tsunami waves

The life of a tsunami is usually divided into three phases: the generation (tsunami source), the propagation and the inundation. Each phase is complex and often described separately. A brief description of each phase is given. Model problems are identified. Their formulation is given. While some of these problems can be solved analytically, most require numerical techniques. The inundation phase is less documented than the other phases. It is shown that methods based on Smoothed Particle Hydrodynamics (SPH) are particularly well-suited for the inundation phase. Directions for future research are outlined.

preprint2006arXiv

Linear theory of wave generation by a moving bottom

The computation of long wave propagation through the ocean obviously depends on the initial condition. When the waves are generated by a moving bottom, a traditional approach consists in translating the ``frozen'' sea bed deformation to the free surface and propagating it. The present study shows the differences between the classical approach (passive generation) and the active generation where the bottom motion is included. The analytical solutions presented here exhibit some of the drawbacks of passive generation. The linearized solutions seem to be sufficient to consider the generation of water waves by a moving bottom.

Denys Dutykh

Quick read

Decide how to stay connected

How to connect with this researcher

Open a focused conversation when the fit is right

See the researcher in context

Building this graph slice

131 published item(s)

Long wave interaction with a partially immersed body. Part II: Numerical results

On a class of lacunary almost Newman polynomials modulo p and density theorems

A Computational Simulation of Steady Natural Convection in an H-form Cavity

An optimal scaling to computationally tractable dimensionless models: Study of latex particles morphology formation

Coupling Conditions for Water Waves at Forks

Derivation of an intermediate viscous Serre-Green-Naghdi equation

Dispersive shallow water wave modelling. Part I: Model derivation on a globally flat space

Dispersive shallow water wave modelling. Part II: Numerical simulation on a globally flat space

Energy of tsunami waves generated by bottom motion

Extreme Inundation Statistics on a Composite Beach

Extreme wave runup on a vertical cliff

Free Surface Flows in Electrohydrodynamics with a Constant Vorticity Distribution

How to overcome the Courant-Friedrichs-Lewy condition of explicit discretizations?

Nonlinear waves in networks: a simple approach using the sine-Gordon equation

Numerical simulation of conservation laws with moving grid nodes: Application to tsunami wave modelling

Numerical simulation of Feller&#39;s diffusion equation

Numerical simulation of wave impact on a rigid wall using a two-phase compressible SPH method

Observation of the Inverse Energy Cascade in the modified Korteweg-de Vries Equation

Poisson and Symplectic structures, Hamiltonian action, momentum and reduction

Resonance enhancement by suitably chosen frequency detuning

The Whitham Equation with Surface Tension

A comparative study of bi-directional Whitham systems

A new model for simulating heat, air and moisture transport in porous building materials

A spectral method for solving heat and moisture transfer through consolidated porous media

An adaptive simulation of nonlinear heat and moisture transfer as a boundary value problem

An efficient numerical method for a long-term simulation of heat and mass transfer: the case of an insulated rammed earth wall

An efficient numerical model for liquid water uptake in porous material and its parameter estimation

An innovative method to determine optimum insulation thickness based on non-uniform adaptive moving grid

Dispersive and non-dispersive nonlinear long wave transformations: Numerical and experimental results

Dispersive effects during long wave run-up on a plane beach

Effects of vorticity on the travelling waves of some shallow water two-component systems

Evaluation of the reliability of building energy performance models for parameter estimation

Experimental and numerical study of the propagation of focused wave groups in the nearshore zone

Hamiltonian regularisation of shallow water equations with uneven bottom

Interpreting Distortions in Dimensionality Reduction by Superimposing Neighbourhood Graphs

Long wave interaction with a partially immersed body. Part I: Mathematical models

Nonlinear deformation and run-up of single tsunami waves of positive polarity: numerical simulations and analytical predictions

On some model equations for pulsatile flow in viscoelastic vessels

On the comparison of three numerical methods applied to building simulation: finite-differences, RC circuit approximation and a spectral method

On the multi-symplectic structure of Boussinesq-type systems. II: Geometric discretization

On the Reducibility and the Lenticular Sets of Zeroes of Almost Newman Lacunary Polynomials

On time relaxed schemes and formulations for dispersive wave equations

Accurate fast computation of steady two-dimensional surface gravity waves in arbitrary depth

Accurate numerical simulation of moisture front in porous material

Advanced reduced-order models for moisture diffusion in porous media

An efficient method to estimate sorption isotherm curve coefficients

An improved explicit scheme for whole-building hygrothermal simulation

Analysis and improvement of the VTT mold growth model: application to bamboo fiberboard

Asymptotic nonlinear and dispersive pulsatile flow in elastic vessels with cylindrical symmetry

Direct dynamical energy cascade in the modified KdV equation

Dispersive shallow water wave modelling. Part III: Model derivation on a globally spherical geometry

Dispersive shallow water wave modelling. Part IV: Numerical simulation on a globally spherical geometry

Evaluation of tsunami wave energy generated by earthquakes in the Makran subduction zone

Numerical modelling of surface water wave interaction with a moving wall

On the modelling of shallow turbidity flows

On the multi-symplectic structure of Boussinesq-type systems. I: Derivation and mathematical properties

On the solution of coupled heat and moisture transport in porous material

On the velocity of turbidity currents over moderate slopes

Peregrine&#39;s system revisited

Solitary wave solutions and their interactions for fully nonlinear water waves with surface tension in the generalized Serre equations

Solving nonlinear diffusive problems in buildings by means of a Spectral Reduced-Order Model

Tsunami hazard assessment in the Makran subduction zone

Wave dynamics on networks: method and application to the sine-Gordon equation

Weakly singular shock profiles for a non-dispersive regularization of shallow-water equations

New asymptotic heat transfer model in thin liquid films

Non-dispersive conservative regularisation of nonlinear shallow water (and isothermal Euler) equations

On supraconvergence phenomenon for second order centered finite differences on non-uniform grids

On the estimation of moisture permeability and advection coefficients of a wood fibre material using the optimal experiment design approach

On weakly singular and fully nonlinear travelling shallow capillary-gravity waves in the critical regime

Some special solutions to the Hyperbolic NLS equation

Stable explicit schemes for simulation of nonlinear moisture transfer in porous materials

Algebraic method for constructing singular steady solitary waves: A case study

Numerical simulation of Feller's diffusion equation

Peregrine's system revisited