Natural sloshing frequencies and modes in a rectangular tank with a slat-type screen

Being installed in tanks, screens play the role of slosh-suppressing devices which may strongly change resonant sloshing frequencies and yield an extra nonlinear damping due to cross-flow resulting in either flow separation or jet flow. Employing the linear sloshing theory and domain decomposition method, we construct an accurate analytical approximation of the natural sloshing modes in a rectangular tank with a slat-type screen at the tank middle. Two-dimensional irrotational flow of an ideal incompressible liquid is assumed. Because the considered flow model does not account for flow separation and jet flow at the screen, the velocity field is locally singular at the sharp edges. The constructed solution captures this singularity. Analyzing this solution establishes a complex dependence of the natural sloshing frequencies on the solidity ratio, the number of submerged screen gaps, the liquid depth, and the position of perforated openings relative to the mean free surface. Results are compared with experimental data. Natural surface wave profiles are discussed in the context of a jump of the velocity potential at the screen and the local inflow component to the screen.     

Low-gravity liquid nonlinear sloshing analysis in a tank under pitching excitation

Under pitch excitation, the sloshing of liquid in circular cylindrical tank includes planar motion, rotary motion and rotary motion inside planar motion. The boundaries between stable motion and unstable motion depend on the radius of the tank, the liquid height, the gravitational intension, the surface tensor and the sloshing damping. In this article, the differential equations of nonlinear sloshing are built first. And by variational principle, the Lagrange function of liquid pressure is constructed in volume integration form. Then the velocity potential function is expanded in series by wave height function at the free surface. The nonlinear equations with kinematics and dynamics free surface boundary conditions through variation are derived. At last, these equations are solved by multiple-scales method. The influence of Bond number on the global stable response of nonlinear liquid sloshing in circular cylinder tank is analyzed in detail. The result indicates that the system's amplitude–frequency response changes from a ‘soft-spring’ to a ‘hard-spring’ in the planar motion with the decreasing of the Bond number, while it changes from a ‘hard-spring’ to a ‘soft-spring’ in the rotary motion. At the same time, jump, lag and other nonlinear phenomena of liquid sloshing are discovered.     

A numerical study of the effects of the vertical baffle on liquid sloshing in two-dimensional rectangular tank

The liquid sloshing in a moving partially filled rectangular tank with a vertical baffle is investigated. A numerical algorithm based on the volume of fluid (VOF) technique is used to study the nonlinear behavior of liquid sloshing. The numerical model solves the complete Navier–Stokes equations in primitive variables by using of finite difference approximations with the moving coordinate system. The ratio of baffle height to the initial liquid depth has been changed in the range of 0≤h_B/h≤1.2. The critical baffle height to reach the roof of the tank and the baffle height beyond the liquid does not get over the baffle anymore have been investigated. The vortex originated by the flow separation from the baffle tip became weaker with increasing the baffle height. In order to assess the accuracy of the method used, some results with baffle height are compared with the experimental results. Comparisons show good agreement for slosh loads in the cases investigated. The free surface elevation and the time variations of pressures have been also presented.     

运行VOF法模拟LNG液舱的晃荡特性

液舱晃荡引起的载荷己成为航行中载液船舶安全性评估的重要内容之一、以及船舶结构动力学的一个热点.文中首先运用VOF法建立LNG晃荡液舱模型,给出对LNG液舱晃荡特性数值模拟方法,运用模拟结果讨论不同工况(如装载率、频率、横摇角度幅值)下晃荡动压特性,在此基础上给出了LNG液舱的晃荡载荷特性与抑制晃荡的措施.     

Frequency response of sloshing in an annular cylindrical tank subjected to pitching excitation

Frequency responses of stable planar and rotary motions in a partially filled annular cylindrical tank, subjected to a pitching excitation at a frequency in the neighborhood of the lowest resonant frequency, are investigated. The nonlinearity of the liquid surface oscillation and the nonlinear coupling between the dominant modes and other modes (e.g., an axisymmetric mode) are considered in the response analysis of the sloshing motion. The basic equations of the liquid motion are derived by using the variational principle and the nonlinear equations of motion of the liquid surface displacement are formulated. The characteristics of the liquid motion in an annular cylindrical tank are discussed. The equations governing the amplitude of the stable planar and rotary liquid motions are derived and the stability of each motion is analyzed. An experiment was carried out using a model tank. It is shown that the nonlinear characteristic of the liquid motion in an annular cylindrical tank is more complicated than that in a circular cylindrical tank. Furthermore, it is shown that the nonlinear analysis is important for estimating the sloshing responses.     

Large-amplitude nonlinear vibrations of a Mooney–Rivlin rectangular membrane

An analysis of the linear and nonlinear vibration response and stability of a pre-stretched hyperelastic rectangular membrane under harmonic lateral pressure and finite initial deformations is presented in this paper. Geometric nonlinearity due to finite deformations and material nonlinearity associated with the hyperelastic constitutive law are taken into account. The membrane is assumed to be made of an isotropic, homogeneous, and incompressible Mooney–Rivlin material. The results for a neo-Hookean material are obtained as a particular case and a comparison of these two constitutive models is carried out. First, the exact solution of the membrane under a biaxial stretch is obtained, being this initial stress state responsible for the membrane stiffness. The equations of motion of the pre-stretched membrane are then derived. From the linearized equations, the natural frequencies and mode shapes of the membrane are analytically obtained for both materials. The natural modes are then used to approximate the nonlinear deformation field using the Galerkin method. A detailed parametric analysis shows the strong influence of the stretching ratios and material parameters on the linear and nonlinear oscillations of the membrane. Frequency–amplitude relations, resonance curves, and bifurcation diagrams, are used to illustrate the nonlinear dynamics of the membrane. The present results are compared favorably with the results evaluated for the same membrane using a nonlinear finite element formulation.     

不同截面水槽流体参数晃动的SPH模拟

用光滑粒子流体动力学方法(SPH)对几种不同截面水槽内的流体参数晃动问题进行数值模拟.以二维矩形截面水槽为例,对其一阶正对称参数晃动过程进行分析,验证这一非线性晃动问题特征:包括水体波动特性,频率关系.对U形和圆形截面的前二阶反对称与正对称参数晃动过程进行数值模拟.结果表明:SPH方法可以有效地实现槽内流体在竖向谐振作用下的参数晃动过程.     

一维Kerr非线性光子晶体中的场分布

研究具有Kerr非线性的一维光子晶体中场的分布.计算光子晶体中加入Kerr非线性材料后,入射光强对光子晶体中场分布极大值位置的调制,发现随着光强的增加,峰值位置向入射端移动.同时,由于非线性的作用,随着光强的变化,线性光量子阱中的电磁本征模将进入线性带隙,带隙中将出现没有线性本征模对应项的新孤子类电磁模.     

Wave turbulent statistics in non-weak wave turbulence

In wave turbulence, which is made by nonlinear interactions among waves, it has been believed that statistical properties are well described by the weak turbulence theory, where separation of linear and nonlinear time scales derived from weak nonlinearity is assumed. However, the separation of the time scales is often violated. To get rid of this inconsistency, closed equations are derived in wave turbulence without assuming the weak nonlinearity according to Direct-Interaction Approximation (DIA), which has been successful in Navier-Stokes turbulence. The DIA equations is a natural extension of the conventional kinetic equation to not-necessarily-weak wave turbulence.     

A numerical study of three-dimensional liquid sloshing in tanks

A numerical model NEWTANK (Numerical Wave TANK) has been developed to study three-dimensional (3-D) non-linear liquid sloshing with broken free surfaces. The numerical model solves the spatially averaged Navier–Stokes equations, which are constructed on a non-inertial reference frame having arbitrary six degree-of-freedom (DOF) of motions, for two-phase flows. The large-eddy-simulation (LES) approach is adopted to model the turbulence effect by using the Smagorinsky sub-grid scale (SGS) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate volume-of-fluid (VOF) method is used to track the distorted and broken free surface. Laboratory experiments are conducted for both 2-D and 3-D non-linear liquid sloshing in a rectangular tank. A linear analytical solution of 3-D liquid sloshing under the coupled surge and sway excitation is also developed in this study. The numerical model is first validated against the available analytical solution and experimental data for 2-D liquid sloshing of both inviscid and viscous fluids. The validation is further extended to 3-D liquid sloshing. The numerical results match with the analytical solution when the excitation amplitude is small. When the excitation amplitude is large where sloshing becomes highly non-linear, large discrepancies are developed between the numerical results and the analytical solutions, the former of which, however, agree well with the experimental data. Finally, as a demonstration, a violent liquid sloshing with broken free surfaces under six DOF excitations is simulated and discussed.     

非线性对大气介质中阵列聚焦声场分布影响的研究

基于时域有限差分算法将大气中近似到二阶微小项的非线性声波波动方程进行离散化,得到了模拟采用的差分波动方程.在此基础上,数值模拟了初始声压强弱不同的5个点声源组成的线阵列垂直或斜向辐射的连续正弦波在大气中传播时二维声场的分布情况.将线性条件下的模拟结果与非线性条件下的模拟结果进行比较后发现:弱非线性会对声场的分布和阵列聚焦增益产生一定的影响,使声场分布波形比线性条件下的声场分布波形更加靠近阵列,聚焦效果变差;强非线性会使波形发生更严重畸变,这是由于产生了基频以外的其他频率声波引起的;非线性对斜向传播时声场分布的影响与垂直传播时的影响效果基本相同,但由于斜向辐射时的声波几何扩展造成的轴向声压衰减要大于垂直辐射时的轴向声压衰减,因此聚焦增益和强非线性的影响都将小于垂直辐射时的情况.     

Energy spectrum of the ensemble of weakly nonlinear gravity-capillary waves on a fluid surface

We consider nonlinear gravity-capillary waves with the nonlinearity parameter ? ～ 0.1–0.25. For this nonlinearity, time scale separation does not occur and the kinetic wave equation does not hold. An energy cascade in this case is built at the dynamic time scale (D-cascade) and is computed by the increment chain equation method first introduced in [15]. We for the first time compute an analytic expression for the energy spectrum of nonlinear gravity-capillary waves as an explicit function of the ratio of surface tension to the gravity acceleration. We show that its two limits—pure capillary and pure gravity waves on a fluid surface—coincide with the previously obtained results. We also discuss relations of the D-cascade model with a few known models used in the theory of nonlinear waves such as Zakharov’s equation, resonance of modes with nonlinear Stokes-corrected frequencies, and the Benjamin-Feir index. These connections are crucial in understanding and forecasting specifics of the energy transport in a variety of multicomponent wave dynamics, from oceanography to optics, from plasma physics to acoustics.     

Recent advances in vibro-impact dynamics and collision of ocean vessels

The treatment of ship impacts and collisions takes different approaches depending on the emphasis of each discipline. For example, dynamicists, physicist, and mathematicians are dealing with developing analytical models and mappings of vibro-impact systems. On the other hand, naval architects and ship designers are interested in developing design codes and structural assessments due to slamming loads, liquid sloshing impact loads in liquefied natural gas tanks and ship grounding accidents. The purpose of this review is to highlight the main differences of the two disciplines. It begins with a brief account of the theory of vibro-impact dynamics based on modeling and mapping of systems experiencing discontinuous changes in their state of motion due to collision. The main techniques used in modeling include power-law phenomenological modeling, Hertzian modeling, and non-smooth coordinate transformations originally developed by Zhuravlev and Ivanov. In view of their effectiveness, both Zhuravlev and Ivanov non-smooth coordinate transformations will be described and assessed for the case of ship roll dynamics experiencing impact with rigid barriers. These transformations have the advantage of converting the vibro-impact oscillator into an oscillator without barriers such that the corresponding equation of motion does not contain any impact term. One of the recent results dealing with the coefficient of restitution is that its value monotonically decreases with the impact velocity and not unique but random in nature. Slamming loads and grounding events of ocean waves acting on the bottom of high speed vessels will be assessed with reference to the ship structural damage. It will be noticed that naval architects and marine engineers are treating these problems using different approaches from those used by dynamicists. The problem of sloshing impact in liquefied natural gas cargo and related problems will be assessed based on the numerical and experimental results. It is important for vessel designers to determine the capacity of ships to resist random slamming loads, sloshing loading impact, grounding accidents and ships collisions.     

Absence of wave packet diffusion in disordered nonlinear systems

We study the spreading of an initially localized wave packet in two nonlinear chains (discrete nonlinear Schr?dinger and quartic Klein-Gordon) with disorder. Previous studies suggest that there are many initial conditions such that the second moment of the norm and energy density distributions diverges with time. We find that the participation number of a wave packet does not diverge simultaneously. We prove this result analytically for norm-conserving models and strong enough nonlinearity. After long times we find a distribution of nondecaying yet interacting normal modes. The Fourier spectrum shows quasiperiodic dynamics. Assuming this result holds for any initially localized wave packet, we rule out the possibility of slow energy diffusion. The dynamical state could approach a quasiperiodic solution (Kolmogorov-Arnold-Moser torus) in the long time limit.     

Analytical prediction of break-out noise from a reactive rectangular plenum with four flexible walls

This paper describes an analytical calculation of break-out noise from a rectangular plenum with four flexible walls by incorporating three-dimensional effects along with the acoustical and structural wave coupling phenomena. The breakout noise from rectangular plenums is important and the coupling between acoustic waves within the plenum and structural waves in the flexible plenum walls plays a critical role in prediction of the transverse transmission loss. The first step in breakout noise prediction is to calculate the inside plenum pressure field and the normal flexible plenum wall vibration by using an impedance-mobility approach, which results in a compact matrix formulation. In the impedance-mobility compact matrix (IMCM) approach, it is presumed that the coupled response can be described in terms of finite sets of the uncoupled acoustic subsystem and the structural subsystem. The flexible walls of the plenum are modeled as an unfolded plate to calculate natural frequencies and mode shapes of the uncoupled structural subsystem. The second step is to calculate the radiated sound power from the flexible walls using Kirchhoff-Helmholtz (KH) integral formulation. Analytical results are validated with finite element and boundary element (FEM-BEM) numerical models.     

Range of applicability of the linear fluid slosh theory for predicting transient lateral slosh and roll stability of tank vehicles

An analytical model is developed to study the transient lateral sloshing in horizontal cylindrical containers assuming inviscid, incompressible and irrotational flows. The model is derived by implementing the linearized free-surface boundary condition and bipolar coordinate transformation, resulting in a truncated system of linear ordinary differential equations, which is numerically solved to determine the fluid velocity potentials followed by the hydrodynamic forces and moment. The model results are compared with those obtained from the multimodal solution. The free-surface elevation and hydrodynamic coefficients are also compared with the reported experimental and analytical data as well as numerical simulations to establish validity of the model. The capability of the model for predicting non-resonant slosh is also evaluated using the critical free-surface amplitude. The model validity is further illustrated by comparing the transient liquid slosh responses of a partially filled tank subject to steady lateral acceleration characterizing a vehicle turning maneuver with those obtained from fully nonlinear CFD simulations and pendulum models. It is shown that the linear slosh model yields more accurate prediction of dynamic slosh than the pendulum models and it is significantly more computationally efficient than the nonlinear CFD model. The slosh model is subsequently applied to roll plane model of a suspended tank vehicle to study the effect of dynamic liquid slosh on steady-turning roll stability limit of the vehicle under constant and variable axle load conditions. The results suggest that the roll moment arising from the dynamic fluid slosh yields considerably lower roll stability limit of the partly-filled tank vehicle compared to that predicted from the widely reported quasi-static fluid slosh model.     

Profiles of initially sinusoidal waves propagating in nonlinear microinhomogeneous materials

The asymptotic analytical theory predicting acoustic wave profiles in microinhomogeneous materials with hysteretic quadratic nonlinearity and attenuation proportional to an even power of frequency is developed. The theory predicts that the influence on the nonlinear wave of the Rayleigh scattering of acoustic waves, which is proportional to the forth power of frequency, results in the net diminishing of wave attenuation. This is due to the suppression (diminishing) by scattering of the nonlinear hysteretic losses which is more important than direct increase in linear losses added by scattering.     

K2_SPH方法及二维破碎波的模拟

通过对一些关键技术的改进,特别是自由表面的判断方法和改进的固壁边界条件,使K2_SPH方法能够模拟破碎波问题,完整实现波浪的爬升、翻转、破碎等过程.对大幅液舱晃荡问题的模拟表明,K2_SPH方法较传统SPH方法在波形和压力分布都有明显改进.     

金属和Kerr非线性介质界面上表面等离子体激元的色散关系

从麦克斯韦方程组出发,结合边界条件,分别得到TM波和TE波在金属和Kerr非线性介质界面上表面等离子体激元的色散关系.由于非线性的存在,TM波的色散关系变得复杂,与光强、非线性系数有关.和线性情况一样,此界面不存在TE波.     

Nonlinear wave mechanics and particulate self-focusing

A previous model for treating electromagnetic nonlinear wave systems is examined in the context of wave mechanics. It is shown that nonlinear wave mechanics implies harmonic generation of new quasiparticle wave functions, which are absent in linear systems. The phenomenon is interpreted in terms of pair (and higher order ensembles) coherence of the interacting particles. The implications are far-reaching, and the present approach might contribute toward a common basis for diverse physical phenomena involving nonlinearity. An intimate relationship connecting coherence, nonlocal interaction, and nonlinearity has been previously noticed in the physics of superconductivity. It is shown here that all these ingredients are consistently contained in the present formalism. The present theory may contribute to elucidate a controversial theory proposed by Panarella, who claims to have measured high-energy photons due to high-intensity laser radiation, which cannot be predicted on the basis of linear quantum theory. Panarella explains the new phenomena by stipulating a nonlinear intensity-dependent photon energy. It is argued here that nonlinearity, manifested in the presence of high intensity, may give rise to high- and low-energy photons, the so-called effective and tired photons, respectively. However, the present explanation does not involve ad hoc assumptions regarding the foundations of quantum theory. In analogy with the electrodynamic model, the present theory leads to particulate self-focusing in high-density streams of particles. Since such particulate beams are currently under consideration in connection with fusion reactions, this might be of future interest.On leave of absence from the Department of Electrical Engineering, Ben-Gurion University of the Negev, Beer Sheva, Israel.