Numerical simulations of the primary Bjerknes force experienced by bubbles in a standing ultrasonic field: Nonlinear vs. linear

The primary Bjerknes force experienced by a population of multiple bubbles in a liquid set in a nonlinear ultrasonic standing field and their translation are calculated and analyzed by numerical simulations. The force field is evaluated by considering the nonlinear bubble oscillations as well as the nonlinear character of the ultrasonic pressure field (both variables are unknown in the coupled nonlinear differential system). The results at small amplitudes agree with the classical theory on bubble translation, depending on the driving frequency in relation to the bubble resonance. It is shown that, when amplitudes are raised, the force field exhibits important modifications that strongly affect the motion of the bubbles and the way they form agglomerates. An analysis is performed on the importance of the terms in the differential system that provoke (a) the nonlinearity of the bubble oscillations and (b) the nonlinearity of the acoustic wave. This study reveals that both features should be considered to better approximate the primary Bjerknes force field. Simulations of the nonlinear ultrasonic field after the bubbles form agglomerates under the influence of this force are also performed.     

Forced oscillations of two gas bubbles in a fluid in the vicinity of bubble contact

For a theoretical derivation of bubble coalescence conditions, nonlinear forced oscillations of two closely spaced spherical bubbles subjected to the action of a periodic external pressure field are considered. The equations, asymptotic with respect to a small distance between the bubble surfaces, are derived to describe the approach of the bubbles under the action of (i) the Bjerknes attraction force averaged over the oscillation period and (ii) the viscous drag. It is shown that due to nonlinear interaction of the viscous drag with the radial and translational oscillations of the bubbles a unidirectional repulsive force is generated, which prevents the approach of the bubbles. The coalescence of the bubbles is possible when the nondimensional parameter combined from the amplitude and frequency of the external pressure field, the bubble radius, and the fluid viscosity is greater than a certain critical value. The obtained coalescence condition is qualitatively confirmed by experiments.     

非线性周期波及破碎波对杆件的作用力

本文采用十五阶 Stokes 波的 Pade 逼近,获得了与实验较为一致的流场,并且利用已有的破碎波的速度、加速度场,计算了非线性波和破碎波对各种杆件作用力,比较了它们的主要特征,为海洋工程设计提供依据。     

Primary pulse transmission in a strongly nonlinear periodic system

The aim of this work is to study the transmission of stress waves in an impulsively forced semi-infinite repetitive system of linear layers which are coupled by means of strongly nonlinear coupling elements. Only primary pulse transmission and reflection at each nonlinear element is considered. This permits the reduction of the problem to an infinite set of first-order strongly nonlinear ordinary differential equations. A subset of these equations is solved both analytically and numerically. For a system possessing clearance nonlinearities it is found that the primary transmitted pulse propagates to only a finite number of layers, and that further transmission of energy to additional layers can occur only through time-delayed secondary pulses or does not occur at all. Hence, clearance nonlinearities in a periodic layered system can lead to energy entrapment in the leading layers. An alternative continuum approximation methodology is also outlined which reduces the problem of primary pulse transmission to the solution of a single strongly nonlinear partial differential equation. The use of the continuum approximation for studying maximum primary pulse penetration in the system with clearance nonlinearities is discussed.     

三维气泡与刚性壁面的相互作用研究

基于势流理论建立水下爆炸气泡运动三维模型,采用边界积分法求解拉普拉斯方程,得到气泡的变形及位置,并在计算过程中引入弹性网格技术,避免了因网格扭曲而导致的数值发散,进而模拟了刚性壁面附近三维气泡的动态特性。在数值模拟过程中,将本文计算值与实验数据进行对比分析,结果表明,计算值与实验数据吻合良好。在此基础上,分别模拟了弱浮力、强Bjerknes力,强浮力、弱Bjerknes力以及浮力与Bjerknes力相当时壁面附近气泡的运动特征,并将各种工况的计算结果与基于开尔文冲量理论(Kelvin Impulse)的Blake准则进行对比分析与讨论,得到了不同参数下气泡的运动特征。     

Finite Element Analysis of Unsteady Natural Convection

In this study, a new finite element method (the MSR-method) is proposed for unsteady three-dimensional thermal-fluid analyses. This method is a combination of a modified Galerkin method (MGM) and the SIMPLER formulation. In the MSR-method the velocity and pressure are computed using the SIMPLER procedure and the approximate velocity and the energy equation are solved using the MGM. In the MGM, the inertia term and the pressure term are considered explicitly, so only the symmetrical matrixes appear. Then an artificial viscosity is introduced through an error analysis approach to improve its accuracy and stability. In this paper, the natural convection problems in a three-dimensional cavity are simulated up to the Rayleigh number of 10⁸, and converged solutions are obtained. Authors confirmed that our proposed method gives reasonable results for these problems comparing with other research works.     

波浪作用下直立结构物附近强湍动掺气流体运动的数值模拟

波浪破碎卷入气体易对建筑物受力产生压力振荡,了解波浪作用下建筑物附近掺气水流的运动特性是精确计算建筑物受力的前提.基于OpenFOAM开源程序包和修正速度入口造波方法建立三维数值波浪水槽,模型采用S-A IDDES湍流模型进行湍流封闭,并采用修正的VOF方法捕捉自由液面,数值模拟了规则波在1:10的光滑斜坡上与直立结构物的相互作用过程,重点分析了结构物附近的水动力和掺气水流运动特性.结果表明,建立的数值模型能精确地捕捉波浪作用下直立结构物附近的自由液面的变化以及气泡输运过程,较好地描述气体卷入所形成的气腔形态以及多气腔之间的融合、分裂等过程;波浪与直立结构物相互作用产生强湍动掺气水流,其运动过程十分复杂;掺气流体输运过程中水气界面周围一直伴随着涡的存在,其中,气泡的分裂与周围正负涡量剪切作用密切相关,且其输运轨迹主要受周围流场的影响;研究揭示了结构物附近湍动能与掺气特性的关系,发现波浪作用下直立结构物附近湍动能的分布与掺气水流特征参数(气泡数量、空隙率)整体呈现一定的线性关系.     

Invariant submodels and exact solutions of Khokhlov–Zabolotskaya–Kuznetsov model of nonlinear hydroacoustics with dissipation

We study three-dimensional Khokhlov–Zabolotskaya–Kuznetsov (KZK) model of the nonlinear hydroacoustics with dissipation. This model is described by third order quasilinear partial differential equation of the (KZK). We obtained that the (KZK) equation admits an infinite Lie group of the transformations, depending on the three arbitrary functions. This is due to the fact that in the (KZK) model the main direction of the wave’s propagation is singled out. The submodels of the (KZK) model.are described by the invariant solutions of the (KZK) equation. We studied essentially distinct, not linked by means of the point transformations, invariant solutions of rank 0 and 1 of this equation. Also we considered the invariant solutions of rank 2 and 3. The invariant solutions of rank 0 and 1 are found either explicitly, or their search is reduced to the solution of the nonlinear integro-differential equations. For example, we obtained the invariant solutions that we called by “Ultrasonic knife” and “Ultrasonic destroyer”. The submodel “Ultrasonic knife” have the following property: at each fixed moment of the time in the field of the existence of the solution near a some plane the pressure increases indefinitely and becomes infinite on this plane. The submodel “Ultrasonic destroyer” contains a countable number of “Ultrasonic knives”. The presence of the arbitrary constants in the integro-differential equations, that determine invariant solutions of rank 1 provides a new opportunities for analytical and numerical study of the boundary value problems for the received submodels, and, thus, for the original (KZK) model. With a help of these invariant solutions we researched a propagation of the intensive acoustic waves (one-dimensional, axisymmetric and planar) for which the acoustic pressure, speed and acceleration of its change, or the acoustic pressure , speed and acceleration of its change in the radial direction, or the acoustic pressure, speed and acceleration of its change in the direction of one of the axes are specified at the initial moment of the time at a fixed point. Under the certain additional conditions, we established the existence and the uniqueness of the solutions of boundary value problems, describing these wave processes. Mechanical relevance of the obtained solutions is as follows: (1) these solutions describe nonlinear and diffraction effects in ultrasonic fields of a special kind, (2) these solutions can be used as a test solutions in the numerical calculations performed in studies of ultrasonic fields generated by powerful emitters. Application of the obtained formula generating the new solutions for the found solutions gives families of the solutions containing three arbitrary functions.     

Invariant solutions of the Westervelt model of nonlinear hydroacoustics without dissipation

We study three-dimensional Westervelt model of a nonlinear hydroacoustics without dissipation. We received all of its invariant submodels. We studied all invariant submodels described by the invariant solutions of rank 0 and 1. All invariant solutions of rank 0 and 1 are found either explicitly, or their search is reduced to the solution of the nonlinear integral equations. With a help of these invariant solutions we researched: (1) a propagation of the intensive acoustic waves (self-similar, axisymmetric, planar and one-dimensional) for which the acoustic pressure and a speed of its change, or the acoustic pressure and its derivative in the direction of one of the axes are specified at the initial moment of the time at a fixed point , (2) a spherically symmetric ultrasonic field for which the acoustic pressure and a speed of its change, or the acoustic pressure and its radial derivative are specified at the initial moment of the time at a fixed point. Solving of the boundary value problems describing these processes is reduced to the solving of nonlinear integral equations. We are established the existence and uniqueness of solutions of these boundary value problems under some additional conditions. Mechanical relevance of the obtained solutions is as follows: (1) these solutions describe nonlinear and diffraction effects in ultrasonic fields of a special kind, (2) these solutions can be used as a test solutions in the numerical calculations performed in studies of ultrasonic fields generated by powerful emitters. We found all the conservation laws of the first order for the Westerveld equation written in dimensionless variables.     

Axisymmetric spreading of a thin power-law fluid under gravity on a horizontal plane

Separable solutions admitted by a nonlinear partial differential equation modeling the axisymmetric spreading under gravity of a thin power-law fluid on a horizontal plane are investigated. The model equation is reduced to a highly nonlinear second-order ordinary differential equation for the spatial variable. Using the techniques of Lie group analysis, the nonlinear ordinary differential equation is linearized and solved. As a consequence of this linearization, new results are obtained.     

Numerical investigations into effects of three-dimensional wake patterns on unsteady aerodynamic characteristics of a circular cylinder at Re=1.3×105

In contrast with a wide range of applications concerning flows around a circular cylinder at upper subcritical Reynolds numbers (Re), there is no systematic understanding about the fundamentals of so-called random flow patterns, and their effects on intermittent modulations in the time history of pressure or force, and the decrease in their spanwise correlations. This paper employed the large-eddy simulation (LES) technique to predict flows past a circular cylinder at Re=1.3×10⁵ and to provide images based on flow visualization that can clarify the physical mechanism responsible for these outcomes. A reasonably sufficient spanwise length was adopted for the numerical model by taking into consideration the effect of aspect ratios (the spanwise length to the diameter). We found that even at such high Res, a three-dimensional pattern of vortical field is present in the wake resulting in total force modulation and weak spanwise correlation, e.g., obvious oblique shedding. The whole development process of the three-dimensional wake is exhibited as a universal. The results revealed that local phase variations in primary vortex shedding are the starting points of three-dimensional wake patterns, which are induced by the “irregular” streamwise vortex. The three-dimensional near wake following local phase variations is associated with a successive evolution composed of certain stages in order. Quantitative analyses based on the time series of sectional lift coefficients show that intermittent increase in primary shedding periods and sectional lift streak divisions are closely related to local phase variations and vortex division in the development process of the three-dimensional pattern. In addition to the phase difference along the span, the three-dimensional pattern also weakens vortex shedding in cross sections perpendicular to the axis of the cylinder, resulting in modulation of the sectional lift coefficient.     

Nonlinear Pressure Diffusion in Flow of Compressible Liquids Through Porous Media

In the flow of liquids through porous media, nonlinear effects arise from the dependence of the fluid density, porosity, and permeability on pore pressure, which are commonly approximated by simple exponential functions. The resulting flow equation contains a squared gradient term and an exponential dependence of the hydraulic diffusivity on pressure. In the limiting case where the porosity and permeability moduli are comparable, the diffusivity is constant, and the squared gradient term can be removed by introducing a new variable y, depending exponentially on pressure. The published transformations that have been used for this purpose are shown to be special cases of the Cole–Hopf transformation, differing in the choice of integration constants. Application of Laplace transformation to the linear diffusion equation satisfied by y is considered, with particular reference to the effects of the transformation on the boundary conditions. The minimum fluid compressibilities at which nonlinear effects become significant are determined for steady flow between parallel planes and cylinders at constant pressure. Calculations show that the liquid densities obtained from the simple compressibility equation of state agree to within 1% with those obtained from the highly accurate Wagner-Pru?  equation of state at pressures to 20 MPa and temperatures approaching 600 K, suggesting possible applications to some geothermal systems.     

N表面张力对近固壁二空化泡影响的数值研究

在忽略浮力下，用边界积分方法数值模拟了表面张力对固壁之上且靠近固壁的二轴对称空化泡生长和溃灭的影响，发现在下空泡最大等效半径为上空泡一半情形，若固壁对下空泡的Bjerknes力大于上空泡对下空泡的Bjerknes力，则表面张力的作用将使下空泡溃灭加速，使其向下的液体射流变强变宽；若固壁对下空泡的Bjerknes力小于上空泡对下空泡的Bjerknes力，则表面张力的作用将使下空泡溃灭变慢，使其向上射流变弱变细长；若这两个Bjerknes力近于相等，则表面张力将会对下空泡溃灭有重大作用，如改变下空泡射流的方向甚至形式(如由环状变向下或由向上变环状)，当上空泡等于或小于下空泡时，表面张力将不会对这两个空泡的行为产生显著影响，定性地分析了表面张力作用的机理。     

A Conceptual Study of Cavity Aeroacoustics Control Using Porous Media Inserts

In this study, an integrated flow simulation and aeroacoustics prediction methodology is applied to testing a sound control technique using porous inserts in an open cavity. Large eddy simulation (LES) combined with a three-dimensional Ffowcs Williams–Hawkings (FW–H) acoustic analogy is employed to predict the flow field, the acoustic sources and the sound radiation. The Darcy pressure – velocity law is applied to conceptually mimic the effect of porous media placed on the cavity floor and/or rear wall. Consequently, flow in the cavity could locally move in or out through these porous walls, depending on the local pressure differences. LES with “standard” subgrid-scale models for compressible flow is carried out to simulate the flow field covering the sound source and near fields, and the fully three-dimensional FW–H acoustic analogy is used to predict the sound field. The numerical results show that applying the conceptual porous media on cavity floor and/or rear wall could decrease the pressure fluctuations in the cavity and the sound pressure level in the far field. The amplitudes of the dominant oscillations (Rossiter modes) are suppressed and their frequencies are slightly modified. The dominant sound source is the transverse dipole term, which is significantly reduced due to the porous walls. As a result, the sound pressure in the far field is also suppressed. The preliminary study reveals that using porous-inserts is a promising technology for flow and sound radiation control.     

Secondary Acoustic Waves in a Polydisperse Bubbly Medium

Wave interaction of bubbles in a cavitation cluster is studied theoretically and experimentally. It is found that two groups of bubbles with different phases of oscillations and collapse time are formed at the early stage of bubblecluster formation. It is shown that small bubbles accumulating energy in the initial pulsed acoustic wave collapse in the field of internal positive pressure of the cluster. A secondary compression wave generated by inertial expansion of large bubbles is determined by means of numerical simulation and registered experimentally. The amplitude of this wave is comparable to the amplitude of the initial pulse. A decrease in the oscillation period and strengthening of the microbubble collapse are observed in the cluster under pulsed compression.     

Propagation of surface (seismic) waves: Ordinary differential equations with strongly nonlinear damping

Second-order ordinary differential equations (ODEs) with strongly nonlinear damping (cubic nonlinearities) govern surface wave motions that entail nonlinear surface seismic motions. They apply to dynamic crack propagation and nonlinear oscillation problems in physics and nonlinear mechanics. It is shown that the nonlinear surface seismic wave equation (Rayleigh equation) admits several functional transformations and it is possible to reduce it to an equivalent first-order Abel ODE of the second kind in normal form. Based on a recently developed methodology concerning the construction of exact analytic solutions for the type of Abel equations under consideration, exact solutions are obtained for the nonlinear seismic wave (NLSW) equation for initial conditions of the physical problem. The method employed is general and can be applied to a large class of relevant ODEs in mathematical physics and nonlinear mechanics.     

基于数据驱动的流场控制方程的稀疏识别

利用有限数据建立系统的非线性动力学模型是具有挑战性的重要课题. 数据驱动的稀疏识别方法是近年来发展的从数据识别动力系统控制方程的有效方法. 本文基于数据驱动稀疏识别方法对不同流场的控制方程进行了识别. 采用非线性动力学偏微分方程函数识别(partial differential equations functional identification of nonlinear dynamics, PDE-FIND)方法和最小绝对收缩和选择算子(least absolute shrinkage and selection operator, LASSO)方法对二维圆柱绕流、顶盖驱动方腔流、Rayleigh-Bénard (RB)对流和三维槽道湍流的控制方程进行了识别. 在稀疏识别过程中, 采用直接数值模拟得到的流场数据来计算过完备候选库中的每一项, 候选库中变量最高保留到二次, 变量导数最高保留到二阶, 非线性项最高保留到四阶. 结果发现PDE-FIND方法和LASSO方法对于不含有非线性项的控制方程, 如涡量输运方程、热输运方程和连续性方程, 都能准确识别. 对于含有强非线性项的控制方程, 如Navier-Stokes方程的识别, PDE-FIND方法正确地识别出了控制方程及流场的Rayleigh数和Reynolds数, 而LASSO方法识别结果不正确, 这是因为候选库中的项之间存在分组效应, LASSO方法通常只取分组中的一项. 本文还发现选择流动结构丰富的区域的数据进行控制方程的稀疏识别可以提高识别的准确性.      

Effect of damping and wave parameters on offshore structure under random excitation

This paper describes the linearized and nonlinear dynamic response of a tension leg platform (TLP) to random waves and current forces. The forcing term of the equation of motion is inherently nonlinear due to the nonlinear drag force. Two analysis procedures are used: nonlinear time domain analysis and linear frequency domain analysis. For the nonlinear analysis, the random wave particle velocities and accelerations are simulated for a given wave spectrum. The nonlinear equation of motion is then integrated directly to obtain the system response statistics. For the linear frequency domain analysis, the nonlinear drag force is linearized through an introduction of linearization coefficients. The main objective of this paper is to investigate the effect of the structural damping and wave parameters on both nonlinear and linear dynamic response of the TLP by parametric studies. The results of stochastic nonlinear and linear dynamic response of the TLP, with and without the presence of current, are presented and compared.     

Free Convection Boundary Layer Flow Past a Vertical Surface in a Porous Medium with Temperature-Dependent Properties

Numerical solutions for the free convection heat transfer in a viscous fluid at a permeable surface embedded in a saturated porous medium, in the presence of viscous dissipation with temperature-dependent variable fluid properties, are obtained. The governing equations for the problem are derived using the Darcy model and the Boussinesq approximation (with nonlinear density temperature variation in the buoyancy force term). The coupled non-linearities arising from the temperature-dependent density, viscosity, thermal conductivity, and viscous dissipation are included. The partial differential equations of the model are reduced to ordinary differential equations by a similarity transformation and the resulting coupled, nonlinear ordinary differential equations are solved numerically by a second order finite difference scheme for several sets of values of the parameters. Also, asymptotic results are obtained for large values of | f _w|. Moreover, the numerical results for the velocity, the temperature, and the wall-temperature gradient are presented through graphs and tables, and are discussed. It is observed that by increasing the fluid variable viscosity parameter, one could reduce the velocity and thermal boundary layer thickness. However, quite the opposite is true with the non-linear density temperature variation parameter.