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

We present in this work a new reconstruction scheme, so-called MUSCL-THINC-BVD scheme, to solve the five-equation model for interfacial two phase flows. This scheme employs the traditional shock capturing MUSCL (Monotone Upstream-centered Schemes for Conservation Law) scheme as well as the interface sharpening THINC (Tangent of Hyperbola for INterface Capturing) scheme as two building-blocks of spatial reconstruction using the BVD (boundary variation diminishing) principle that minimizes the variations (jumps) of the reconstructed variables at cell boundaries, and thus effectively reduces the numerical dissipations in numerical solutions. The MUSCL-THINC-BVD scheme is implemented to all state variables and volume fraction, which realizes the consistency among volume fraction and other physical variables. Read More

Laterally periodic nanostructures were investigated with grazing incidence small angle X-ray scattering (GISAXS) by using the diffraction patterns to reconstruct the surface shape. To model visible light scattering, rigorous calculations of the near and far field by numerically solving Maxwell's equations with a finite-element method are well established. The application of this technique to X-rays is still challenging, due to the discrepancy between incident wavelength and finite-element size. Read More

This paper proposes the use of the Spectral method to simulate diffusive moisture transfer through porous materials, which can be strongly nonlinear and can significantly affect sensible and latent heat transfer. An alternative way for computing solutions by considering a separated representation is presented, which can be applied to both linear and nonlinear diffusive problems, considering highly moisture-dependent properties. The Spectral method is compared with the classical implicit Euler and Crank-Nicolson schemes. Read More

We have developed the first laser damage simulation algorithm capable of determining crater and surface modification morphology from microscopic physics. Rapid progress in the field of high intensity ultrafast lasers has enabled its utility in a myriad of applications. Simulation plays an important role in this research by allowing for closer analysis of the physical mechanisms involved, but current techniques struggle to meet both the spatial scope or resolution requirements for modeling such dynamics, typically specializing in one or the other. Read More

A multi-scale framework was recently proposed for more realistic molecular dynamics simulations in continuum solvent models by coupling a molecular mechanics treatment of solute with a fluid mechanics treatment of solvent, where we formulated the physical model and developed a numerical fluid dynamics integrator. In this study, we incorporated the fluid dynamics integrator with the Amber simulation engine to conduct atomistic simulations of biomolecules. At this stage of the development, only nonelectrostatic interactions, i. Read More

A self-consistent procedure for ion-induced soft error rate calculation in space environment taking into account Error Correcting Codes is proposed. The method is based on partitioning of multiple cell events into groups with different multiplicities. Read More

A procedure of space proton-induced Soft Error Rate calculation based on heavy ion testing data is proposed. The approach relies on Geant4-assisted Monte Carlo simulation of the secondary particle LET spectra produced by proton interactions. Read More

We recalculate the leading relativistic corrections for the ground electronic state of the hydrogen molecule using variational method with explicitly correlated functions which satisfy the interelectronic cusp condition. The new computational approach allowed for the control of the numerical precision which reached about 8 significant digits. More importantly, the updated theoretical energies became discrepant with the known experimental values and we conclude that the yet unknown relativistic recoil corrections might be larger than previously anticipated. Read More

We evaluate the practical performance of the zero-multiple summation method (ZMM), a method for approximately calculating electrostatic interactions in molecular dynamics simulations. The performance of the ZMM is compared with that of the smooth particle mesh Ewald method (SPME). Even though the ZMM uses a larger cutoff distance than the SPME does, the performance of the ZMM is found to be comparable to or better than that of the SPME. Read More

There is considerable interest in collective effects in hybrid systems formed by molecular or atomic ensembles strongly coupled by an electromagnetic resonance. For analyzing such collective effects, we develop an efficient and general theoretical formalism based on the natural modes of the resonator. The main strength of our approach is its generality and the high level of analyticity enabled by modal analysis, which allows one to model complex hybrid systems without any restriction on the resonator shapes or material properties, and to perform statistical computations to predict general properties that are robust to spatial and polarization disorders. Read More