Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDispersive shallow water wave modelling. Part IV: Numerical simulation on a globally spherical geometry
53
0
0.0
(
0
)
Authors
Gayaz Khakimzyanov
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
The multiscale patch scheme is built from given small micro-scale simulations of complicated physical processes to empower large macro-scale simulations. By coupling small patches of simulations over unsimulated spatial gaps, large savings in computational time are possible. Here we discuss generalising the patch scheme to the case of wave systems on staggered grids in 2D space. Classic macro-scale interpolation provides a generic coupling between patches that achieves arbitrarily high order consistency between the emergent macro-scale simulation and the underlying micro-scale dynamics. Eigen-analysis indicates that the resultant scheme empowers feasible computation of large macro-scale simulations of wave systems even with complicated underlying physics. As examples we use the scheme to simulate some wave scenarios via a turbulent shallow water model.
We present a parareal in time algorithm for the simulation of neutron diffusion transient model. The method is made efficient by means of a coarse solver defined with large time steps and steady control rods model. Using finite element for the space discretization, our implementation provides a good scalability of the algorithm. Numerical results show the efficiency of the parareal method on large light water reactor transient model corresponding to the Langenbuch-Maurer-Werner (LMW) benchmark [1].
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.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
