Efficient uncertainty quantification of a fully nonlinear and dispersive water wave model with random inputs

A major challenge in next-generation industrial applications is to improve numerical analysis by quantifying uncertainties in predictions. In this work we present a formulation of a fully nonlinear and dispersive potential flow water wave model with random inputs for the probabilistic description of the evolution of waves. The model is analyzed using random sampling techniques and non-intrusive methods based on generalized Polynomial Chaos (PC). These methods allow to accurately and efficiently estimate the probability distribution of the solution and require only the computation of the solution in different points in the parameter space, allowing for the reuse of existing simulation software. The choice of the applied methods is driven by the number of uncertain input parameters and by the fact that finding the solution of the considered model is computationally intensive. We revisit experimental benchmarks often used for validation of deterministic water wave models. Based on numerical experiments and assumed uncertainties in boundary data, our analysis reveals that some of the known discrepancies from deterministic simulation in comparison with experimental measurements could be partially explained by the variability in the model input. We finally present a synthetic experiment studying the variance based sensitivity of the wave load on an off-shore structure to a number of input uncertainties. In the numerical examples presented the PC methods have exhibited fast convergence, suggesting that the problem is amenable to being analyzed with such methods.

Similar Publications

In the standard SPH method, the interaction between two particles might be not pairwise when the support domain varies, which can result in a reduction of accuracy. To deal with this problem, a modified SPH approach is presented in this paper. First of all, a Lagrangian kernel is introduced to eliminate spurious distortions of the domain of material stability, and the gradient is corrected by a linear transformation so that linear completeness is satisfied. Read More

It is generally accepted that all models are wrong -- the difficulty is determining which are useful. Here, a useful model is considered as one that is capable of combining data and expert knowledge, through an inversion or calibration process, to adequately characterize the uncertainty in predictions of interest. This paper derives conditions that specify which simplified models are useful and how they should be calibrated. Read More

We investigate the transition to a Landau-Levich-Derjaguin film in forced dewetting using a quadtree adaptive solution to the Navier-Stokes equations with surface tension. A discretization of the capillary forces near the receding contact line is used that yields an equilibrium for a specified contact angle $\theta_\Delta$ called the numerical contact angle. Despite the well-known contact line singularity, dynamic simulations can proceed without any explicit additional numerical procedure, yielding an implicitly dynamic contact angle model. Read More

Ab initio quantum chemistry calculations for systems with large active spaces are notoriously difficult and cannot be successfully tackled by standard methods. In this letter, we generalize a Green's function QM/QM embedding method called self-energy embedding theory (SEET) that has the potential to be successfully employed to treat large active spaces. In generalized SEET, active orbitals are grouped into intersecting groups of few orbitals allowing us to perform multiple parallel calculations yielding results comparable to the full active space treatment. Read More

Drag laws for particles in fluids are often expressed in terms of the undisturbed fluid velocity, defined as the fluid velocity a particle sees before the disturbance develops in the fluid. In two-way coupled point-particle simulations the information from the undisturbed state is not available and must be approximated using the disturbed velocity field. Horwitz and Mani (2016) recently developed a procedure to estimate the undisturbed velocity for particles moving at low Reynolds number and obeying the Stokes drag law. Read More

The Wannier localization problem in quantum physics is mathematically analogous to finding a localized representation of a subspace corresponding to a nonlinear eigenvalue problem. While Wannier localization is well understood for insulating materials with isolated eigenvalues, less is known for metallic systems with entangled eigenvalues. Currently, the most widely used method for treating systems with entangled eigenvalues is to first obtain a reduced subspace (often referred to as disentanglement) and then to solve the Wannier localization problem by treating the reduced subspace as an isolated system. Read More

We provide the necessary framework for carrying out stochastic positive-P and gauge-P simulations of bosonic systems with long range interactions. In these approaches, the quantum evolution is sampled by trajectories in phase space, allowing calculation of correlations without truncation of the Hilbert space or other approximations to the quantum state. The main drawback is that the simulation time is limited by noise arising from interactions. Read More

A new highly efficient method is developed for computation of traveling periodic waves (Stokes waves) on the free surface of deep water. A convergence of numerical approximation is determined by the complex singularites above the free surface for the analytical continuation of the travelling wave into the complex plane. An auxiliary conformal mapping is introduced which moves singularities away from the free surface thus dramatically speeding up numerical convergence by adapting the numerical grid for resolving singularities while being consistent with the fluid dynamics. Read More

Machine learning has proven to be a valuable tool to approximate functions in high-dimensional spaces. Unfortunately, analysis of these models to extract the relevant physics is never as easy as applying machine learning to a large dataset in the first place. Here we present a description of atomic systems that generates machine learning representations with a direct path to physical interpretation. Read More

We developed a combined atomistic-continuum hierarchical multiscale approach to explore the effective thermal conductivity of graphene laminates. To this aim, we first performed molecular dynamics simulations in order to study the heat conduction at atomistic level. Using the non-equilibrium molecular dynamics method, we evaluated the length dependent thermal conductivity of graphene as well as the thermal contact conductance between two individual graphene sheets. Read More