Numerical simulation of fully nonlinear interaction between steep waves and 2D floating bodies using the QALE-FEM method

This paper extends the QALE-FEM (quasi arbitrary Lagrangian–Eulerian finite element method) based on a fully nonlinear potential theory, which was recently developed by the authors [Q.W. Ma, S. Yan, Quasi ALE finite element method for nonlinear water waves, J. Comput. Phys, 212 (2006) 52–72; S. Yan, Q.W. Ma, Application of QALE-FEM to the interaction between nonlinear water waves and periodic bars on the bottom, in: 20th International Workshop on Water Waves and Floating Bodies, Norway, 2005], to deal with the fully nonlinear interaction between steep waves and 2D floating bodies. In the QALE-FEM method, complex unstructured mesh is generated only once at the beginning of calculation and is moved to conform to the motion of boundaries at other time steps, avoiding the necessity of high cost remeshing. In order to tackle challenges associated with floating bodies, several new numerical techniques are developed in this paper. These include the technique for moving mesh near and on body surfaces, the scheme for estimating the velocities and accelerations of bodies as well as the forces on them, the method for evaluating the fluid velocity on the surface of bodies and the technique for shortening the transient period. Using the developed techniques and methods, various cases associated with the nonlinear interaction between waves and floating bodies are numerically simulated. For some cases, the numerical results are compared with experimental data available in the public domain and good agreement is achieved

Quasi ALE finite element method for nonlinear water waves

This paper presents a newly developed quasi arbitrary Lagrangian-Eulerian finite element method (QALE-FEM) for simulating water waves based on fully nonlinear potential theory. The main difference of this method from the conventional finite element method developed by one of authors of this paper and others (see, e.g., [11] and [22]) is that the complex mesh is generated only once at the beginning and is moved at all other time steps in order to conform to the motion of the free surface and structures. This feature allows one to use an unstructured mesh with any degree of complexity without the need of regenerating it every time step, which is generally inevitable and very costly. Due to this feature, the QALE-FEM has high potential in enhancing computational efficiency when applied to problems associated with the complex interaction between large steep waves and structures since the use of an unstructured mesh in such a case is likely to be necessary. To achieve overall high efficiency, the numerical techniques involved in the QALE-FEM are developed, including the method to move interior nodes, technique to re-distribute the nodes on the free surface, scheme to calculate velocities and so on. The model is validated by water waves generated by a wavemaker in a tank and the interaction between water waves and periodic bars on the bed of tank. Satisfactory agreement is achieved with analytical solutions, experimental data and numerical results from other methods

Improved MLPG_R method for simulating 2D interaction between violent waves and elastic structures

Interaction between violent water waves and structures is of a major concern and one of the important issues that has not been well understood in marine engineering. This paper will present first attempt to extend the Meshless Local Petrov Galerkin method with Rankine source solution (MLPG_R) for studying such interaction, which solves the Navier-stokes equations for water waves and the elastic vibration mequations for structures under wave impact. The MLPG_R method has been applied successfully to modeling various violent water waves and their interaction with rigid structures in our previous publications. To make the method robust for modeling wave elastic-structure interaction (hydroelasticity) problems concerned here, a near-strongly coupled and partitioned procedure is proposed to deal with coupling between violent waves and dynamics of structures. In addition, a novel approach is adopted to estimate pressure gradient when updating velocities and positions of fluid particles, leading to a relatively smoother pressure time history that is crucial for success in simulating problems about wavestructure interaction. The developed method is used to model several cases, covering a range from small wave to violent waves. Numerical results for them are compared with those obtained from other methods and from experiments in literature. Reasonable good agreement between them is achieved

Effects of the field modulation on the Hofstadter's spectrum

We study the effect of spatially modulated magnetic fields on the energy spectrum of a two-dimensional (2D) Bloch electron. Taking into account four kinds of modulated fields and using the method of direct diagonalization of the Hamiltonian matrix, we calculate energy spectra with varying system parameters (i.e., the kind of the modulation, the relative strength of the modulated field to the uniform background field, and the period of the modulation) to elucidate that the energy band structure sensitively depends on such parameters: Inclusion of spatially modulated fields into a uniform field leads occurrence of gap opening, gap closing, band crossing, and band broadening, resulting distinctive energy band structure from the Hofstadter's spectrum. We also discuss the effect of the field modulation on the symmetries appeared in the Hofstadter's spectrum in detail.Comment: 7 pages (in two-column), 10 figures (including 2 tables

Incompressible SPH method based on Rankine source solution for violent water wave simulation

With wide applications, the smoothed particle hydrodynamics method (abbreviated as SPH) has become an important numerical tool for solving complex flows, in particular those with a rapidly moving free surface. For such problems, the incompressible Smoothed Particle Hydrodynamics (ISPH) has been shown to yield better and more stable pressure time histories than the traditional SPH by many papers in literature. However, the existing ISPH method directly approximates the second order derivatives of the functions to be solved by using the Poisson equation. The order of accuracy of the method becomes low, especially when particles are distributed in a disorderly manner, which generally happens for modelling violent water waves. This paper introduces a new formulation using the Rankine source solution. In the new approach to the ISPH, the Poisson equation is first transformed into another form that does not include any derivative of the functions to be solved, and as a result, does not need to numerically approximate derivatives. The advantage of the new approach without need of numerical approximation of derivatives is obvious, potentially leading to a more robust numerical method. The newly formulated method is tested by simulating various water waves, and its convergent behaviours are numerically studied in this paper. Its results are compared with experimental data in some cases and reasonably good agreement is achieved. More importantly, numerical results clearly show that the newly developed method does need less number of particles and so less computational costs to achieve the similar level of accuracy, or to produce more accurate results with the same number of particles compared with the traditional SPH and existing ISPH when it is applied to modelling water waves

A hybrid method for modelling two dimensional non-breaking and breaking waves

This is the first paper to present a hybrid method coupling a Improved Meshless Local Petrov Galerkin method with Rankine source solution (IMLPG_R) based on the Navier Stokes (NS) equations, with a finite element method (FEM) based on the fully nonlinear potential flow theory (FNPT) in order to efficiently simulate the violent waves and their interaction with marine structures. The two models are strongly coupled in space and time domains using a moving overlapping zone, wherein the information from both the solvers is exchanged. In the time domain, the Runge-Kutta 2nd order method is nested with a predictor-corrector scheme. In the space domain, numerical techniques including ‘Feeding Particles’ and two-layer particle interpolation with relaxation coefficients are introduced to achieve the robust coupling of the two models. The properties and behaviours of the new hybrid model are tested by modelling a regular wave, solitary wave and Cnoidal wave including breaking and overtopping. It is validated by comparing the results of the method with analytical solutions, results from other methods and experimental data. The paper demonstrates that the method can produce satisfactory results but uses much less computational time compared with a method based on the full NS model

A hybrid model for simulating rogue waves in random seas on a large temporal and spatial scale

A hybrid model for simulating rogue waves in random seas on a large temporal and spatial scale is proposed in this paper. It is formed by combining the derived fifth order Enhanced Nonlinear Schrödinger Equation based on Fourier transform, the Enhanced Spectral Boundary Integral (ESBI) method and its simplified version. The numerical techniques and algorithm for coupling three models on time scale are suggested. Using the algorithm, the switch between the three models during the computation is triggered automatically according to wave nonlinearities. Numerical tests are carried out and the results indicate that this hybrid model could simulate rogue waves both accurately and efficiently. In some cases discussed, the hybrid model is more than 10 times faster than just using the ESBI method, and it is also much faster than other methods reported in the literature

BcB_c spectroscopy

In the framework of potential models for heavy quarkonium the mass spectrum for the system ($\bar b c$) is considered. Spin-dependent splittings, taking into account a change of a constant for effective Coulomb interaction between the quarks, and widths of radiative transitions between the ($\bar b c$) levels are calculated. In the framework of QCD sum rules, masses of the lightest vector $B_c^*$ and pseudoscalar $B_c$ states are estimated, scaling relation for leptonic constants of heavy quarkonia is derived, and the leptonic constant $f_{B_C}$ is evaluated.Comment: IHEP 94-51, LATEX, 39 page

ϒ production in p–Pb collisions at √sNN=8.16 TeV

ϒ production in p–Pb interactions is studied at the centre-of-mass energy per nucleon–nucleon collision √sNN = 8.16 TeV with the ALICE detector at the CERN LHC. The measurement is performed reconstructing bottomonium resonances via their dimuon decay channel, in the centre-of-mass rapidity intervals 2.03 < ycms < 3.53 and −4.46 < ycms < −2.96, down to zero transverse momentum. In this work, results on the ϒ(1S) production cross section as a function of rapidity and transverse momentum are presented. The corresponding nuclear modification factor shows a suppression of the ϒ(1S) yields with respect to pp collisions, both at forward and backward rapidity. This suppression is stronger in the low transverse momentum region and shows no significant dependence on the centrality of the interactions. Furthermore, the ϒ(2S) nuclear modification factor is evaluated, suggesting a suppression similar to that of the ϒ(1S). A first measurement of the ϒ(3S) has also been performed. Finally, results are compared with previous ALICE measurements in p–Pb collisions at √sNN = 5.02 TeV and with theoretical calculations.publishedVersio