两层密度分层流体毛细重力波三阶Stokes波解

以小振幅波理论为基础,利用摄动方法研究了两层密度成层状态下的毛细重力波,求得了两层密度成层状态下各层流体速度势的三阶解及毛细重力波波面位移的三阶Stokes波解.结果表明:三阶方程的解均受到表面张力的影响.三阶Stokes波解描述了毛细重力波的三阶非线性修正,波速不仅取决于波数和各层流体的厚度,而且还与波幅及表面张力有关.     

On Stabilizing the Strongly Nonlinear Internal Wave Model

A strongly nonlinear asymptotic model describing the evolution of large amplitude internal waves in a two-layer system is studied numerically. While the steady model has been demonstrated to capture correctly the characteristics of large amplitude internal solitary waves, a local stability analysis shows that the time-dependent inviscid model suffers from the Kelvin–Helmholtz instability due to a tangential velocity discontinuity across the interface accompanied by the interfacial deformation. An attempt to represent the viscous effect that is missing in the model is made with eddy viscosity, but this simple ad hoc model is shown to fail to suppress unstable short waves. Alternatively, when a smooth low-pass Fourier filter is applied, it is found that a large amplitude internal solitary wave propagates stably without change of form, and mass and energy are conserved well. The head-on collision of two counter-propagating solitary waves is studied using the filtered strongly nonlinear model and its numerical solution is compared with the weakly nonlinear asymptotic solution.     

The method of successive approximations for solving the problem of the decay of a discontinuity and its application to the equations of two-layer shallow water theory

The method of successive approximations for solving the problem of the decay of a small amplitude discontinuity is proposed for hyperbolic systems of conservation laws. In the linear approximation, a Cauchy problem for a linear hyperbolic system is obtained. Its solution represents lines of discontinuity separating the regions in which the solution is constant. Most attention is paid to the first and second approximations, within the limits of which the discontinuities obtained in the first approximation decay into stable shock waves and rarefaction waves. An analysis of the qualitatively different flow conditions that arise when solving the problem of the failure of a dam for a two-layer shallow water model with a free boundary is presented as an example.     

有流存在下两层密度分层流体毛细重力波的三阶Stokes波解

以小振幅波理论为基础,利用奇异摄动方法研究了有背景流存在下两层密度成层状态下的毛细重力波,求得了两层密度成层状态下各层流体速度势的三阶解及毛细重力波波面位移的三阶Stokes波解,并讨论了毛细重力波的kelvin-Helmholtz的不稳定性.结果表明在有流存在的情况下,两层密度成层流体毛细重力波的一阶渐近解、频散关系,二阶渐近解及三阶渐近解不仅依赖于各层流体的厚度和密度,也依赖于表面张力和各层流体的背景流流场;毛细重力波的三阶解描述了背景流场与毛细重力波之问的三阶非线性相互作用.对于给定的波数k(实数)毛细重力波可能出现kelvin-Helmholtz不稳定性.     

Effect of interconvertibility of electromagnetic and gravitational waves in strong external electromagnetic fields and the propagation of waves in the field of a charged “black hole” : PMM vol. 38, n 6, 1974, pp. 1122–1129

It is shown that the behavior of an arbitrary wave propagating in the field of a nonrotating charged black hole is defined (with the use of quadratures) by four functions. Each of these functions obeys its second order equation of the wave kind. Short electromagnetic waves falling onto a black hole are reflected by its field in the form of gravitational and electromagnetic waves whose amplitude was explicitly determined. In the case of the wave carrying rays winding around the limit cycle the reflection and transmission coefficients were obtained in the form of analytic expressions.Various physical processes taking place inside, as well as outside a collapsing star, may induce perturbations of the gravitational, electromagnetic and other fields, and lead to the appearance in the surrounding space of waves of various kinds which propagate over a distorted background and are dissipated along its inhomogeneities.In the absence of rotation and charge in a star, the analysis of small perturbations of the gravitational fields is based on the system of Einstein equations linearized around the Schwarzschild solution. In [1, 2] this system of equations, after expansion of perturbations in spherical harmonics and Fourier transformation with respect to time, was reduced to two independent linear ordinary differential equations of second order of the form of the stationary Schrödinger equation for a particle in a potential force field. Each of these equations defines one of two possible independent perturbation kinds: “even” and “odd” (the different behavior of spherical tensor harmonics at coordinate inversion is the deciding factor in the determination of the kind of perturbation [1, 2]). Although these equations were derived with the superposition on the perturbations of the metric of specific coordinate conditions, they define, as shown in [4], the behavior of invariants of the perturbed gravitational field, which imparts to the potential barriers appearing in these equations an invariant meaning.The system of Maxwell equations on the background of Schwarzschild solution also reduces to similar equations, which differ from the above only by the form of potential barriers appearing in these [5].In the presence in the unperturbed solution of a strong electromagnetic field the gravitational and electromagnetic waves interact with each other, and transmutation takes place. The train of short periodic electromagnetic waves generates the accompanying train of gravitational waves. This phenomenon was first analyzed in [6] on and arbitrary background. It was shown in [7, 8] that dense stars surrounded by hot plasma may acquire a charge owing to splitting of charges by radiation pressure and the “sweeping out” of positrons nascent in vapors in strong electrostatic fields. The interaction of waves becomes particularly clearly evident in the neighborhood of black holes which may serve as “valves” by maintaining equilibrium between the relict electromagnetic and gravitational radiation in the Universe. Rotation of black holes intensifies this effect [6].If a nonrotating star possesses an electrostatic charge, the definition of perturbations of the electromagnetic and gravitational fields must be based on the complete system of Einstein-Maxwell equations linearized around the Nordström-Reissner solution. (Small perturbations of electromagnetic field outside a charged black hole were considered in [9, 10] on the basis of the system of Maxwell equations on a “rigid” background of the Nordström-Reissner solution, without taking into account the interconvertibility of gravitational and electromagnetic waves, which materially affects their behavior in the neighborhood of a charged black hole). Here this system of equations which define the interacting gravitational and electromagnetic perturbations are reduced to four independent second order differential equations, two for each kind of perturbations (an importsnt part is played here by the coordinate conditions imposed on the perturbations of the metric, proposed by the authors in [4]). Perturbation components of the metric and of the electromagnetic field are determined in quadratures by the solutions of these equations. If the charge of a star tends to vanish, two of the derived equations convert to equations for gravitational waves on the background of the Schwarzschild solution [1, 2], while the twoothers become equations which are equivalent to Maxwell solutions on the same background. The short-wave asymptotics of derived equations is determined throughout including the neighborhood of the limit cycle for the wave carrying rays. These solutions far away from the point of turn coincide with those obtained in [6] for any arbitrary background. Approximation of geometric optics does not provide correct asymptotics for impact parameters of rays which are close to critical for which the Isotropie and geodesic parameters wind around the limit cycle. This case is investigated below.A similar situation in the Schwarzschild field was analyzed in [11], where analytic expressions for the wave reflection and transmission coefficients were determined, and the integral radiation stream trapped by a black hole produced by another radiation component of the dual system was calculated.     

两层流体中具有β变化和地形影响的Rossby波

本文研究了两层流体中具有变化的Rossby参数和地形Rossby波的问题.利用行波法和摄动的方法,获得了Rossby波振幅满足齐次KdV方程和齐次mKdV方程,推广了Rossby参数和地形对Rossby孤立波的影响.     

Solitary wave solutions for a general Boussinesq type fluid model

The possible solitary wave solutions for a general Boussinesq (GBQ) type fluid model are studied analytically. After proving the non-Painlevé integrability of the model, the first type of exact explicit travelling solitary wave with a special velocity selection is found by the truncated Painlevé expansion. The general solitary waves with different travelling velocities can be studied by casting the problems to the Newtonian quasi-particles moving in some proper one dimensional potential fields. For some special velocity selections, the solitary waves possess different shapes, say, the left moving solitary waves may possess different shapes and/or amplitudes with those of the right moving solitons. For some other velocities, the solitary waves are completely prohibited. There are three types of GBQ systems (GBQSs) according to the different selections of the model parameters. For the first type of GBQS, both the faster moving and lower moving solitary waves allowed but the solitary waves with“middle” velocities are prohibit. For the second type of GBQS all the slower moving solitary waves are completely prohibit while for the third type of GBQS only the slower moving solitary waves are allowed. Only the solitary waves with the almost unit velocities meet the weak non-linearity conditions.     

On linear and nonlinear Rossby waves in an ocean

This paper deals with recent developments of linear and nonlinear Rossby waves in an ocean. Included are also linear Poincaré, Rossby, and Kelvin waves in an ocean. The dispersion diagrams for Poincaré, Kelvin and Rossby waves are presented. Special attention is given to the nonlinear Rossby waves on a β-plane ocean. Based on the perturbation analysis, it is shown that the nonlinear evolution equation for the wave amplitude satisfies a modified nonlinear Schrödinger equation. The solution of this equation represents solitary waves in a dispersive medium. In other words, the envelope of the amplitude of the waves has a soliton structure and these envelope solitons propagate with the group velocity of the Rossby waves. Finally, a nonlinear analytical model is presented for long Rossby waves in a meridional channel with weak shear. A new nonlinear wave equation for the amplitude of large Rossby waves is derived in a region where fluid flows over the recirculation core. It is shown that the governing amplitude equations for the inner and outer zones are both KdV type, where weak nonlinearity is balanced by weak dispersion. In the inner zone, the nonlinear amplitude equation has a new term proportional to the 3/2 power of the difference between the wave amplitude and the critical amplitude, and this term occurs to account for a nonlinearity due to the flow over the vortex core. The solution of the amplitude equations with the linear shear flow represents the solitary waves. The present study deals with the lowest mode (n=1) analysis. An extension of the higher modes (n?2) of this work will be made in a subsequent paper.     

缓变深度分层流体中的准周期波和准孤立波

本文讨论具缓变深度二流体系统中的非线性波,该系统由一不规则底部与一水平固壁间的两层常密度无粘流体所组成.文中用约化摄动法导出了所考虑模型的变系数Korteweg-de Vries方程,并用多重尺度法求出了该方程的近似解,发现底部固壁的不规则变化将产生所谓准周期波和准孤立波.它们的周期、波速和波形将发生缓慢变化,文中给出了准周期波的周期随深度的变化关系式以及准孤立波波幅、波速随深度的变化关系式,底部水平情形和单层流体情形可看成是本文的特例.     

Waves in an inhomogeneous fluid in the presence of a dock : PMM vol. 39, n 6, 1975, pp. 1140–1142

We investigate the propagation of waves generated by oscillations of a section of the bottom of a tank through a two-layer fluid, in the presence of a dock. Wave motions in an inhomogeneous fluid generated by displacement of a section of the bottom of a tank were studied in [1] where the upper surface of the fluid was assumed either to be completely free, or completely covered with ice. In the present paper we use the method given in [2] to investigate a similar problem under the assumption that the fluid surface is partly covered with an immovable rigid plate. The expressions obtained for the velocity potential are used to determine the form of the free surface and of the interface. We show that when the fluid is inhomogeneous, the wave amplitude on the free surface increases, while the presence of a plate reduces the amplitude of the surface waves, as well as of the internal waves in the region between the plate and the oscillating section of the bottom.     

Inhibiting Shear Instability Induced by Large Amplitude Internal Solitary Waves in Two-Layer Flows with a Free Surface

We consider a strongly nonlinear long wave model for large amplitude internal waves in two-layer flows with the top free surface. It is shown that the model suffers from the Kelvin–Helmholtz (KH) instability so that any given shear (even if arbitrarily small) between the layers makes short waves unstable. Because a jump in tangential velocity is induced when the interface is deformed, the applicability of the model to describe the dynamics of internal waves is expected to remain rather limited. To overcome this major difficulty, the model is written in terms of the horizontal velocities at the bottom and the interface, instead of the depth-averaged velocities, which makes the system linearly stable for perturbations of arbitrary wavelengths as long as the shear does not exceed a certain critical value.     

Stokes waves with constant vorticity: I. Numerical computation

Periodic traveling waves are numerically computed in a constant vorticity flow subject to the force of gravity. The Stokes wave problem is formulated via a conformal mapping as a nonlinear pseudodifferential equation, involving a periodic Hilbert transform for a strip, and solved by the Newton‐GMRES method. For strong positive vorticity, in the finite or infinite depth, overhanging profiles are found as the amplitude increases and tend to a touching wave, whose surface contacts itself at the trough line, enclosing an air bubble; numerical solutions become unphysical as the amplitude increases further and make a gap in the wave speed versus amplitude plane; another touching wave takes over and physical solutions follow along the fold in the wave speed versus amplitude plane until they ultimately tend to an extreme wave, which exhibits a corner at the crest. Touching waves connected to zero amplitude are found to approach the limiting Crapper wave as the strength of positive vorticity increases unboundedly, while touching waves connected to the extreme waves approach the rigid body rotation of a fluid disk.     

Explosive instability in nonlinear waves in media with negative viscosity : PMM vol. 38, n 6, 1974, pp. 991–995

The propagation of a wave of a finite amplitude in a medium with a nonlinearity of the second degree and negative viscosity, is examined. It is shown that in a finite time singularities appear in the solution. The exact solution of the Cauchy problem is given for a specific case. Recently the effects of negative viscosity which cause an increase in the energy of the wave motion have been studied intensively in electrodynamics, plasma physics, the Earth's atmosphere, in the theory of the circulation of the oceans and of flow in open channels [1–4], Wave amplification caused by an energy transfer from turbulent to regular motions, is possible in any medium having space-time fluctuations, provided the correlation time is sufficiently small [5, 6]. As the wave amplitude increases, nonlinear effects become important; they have been taken into account in cases where the interaction of a finite number of harmonics [2, 4] and the structure of steady motions have been examined [3].It is shown in this paper that in a medium with negative viscosity and a second degree dynamic nonlinearity, a solution of the Cauchy problem for an arbitrary “good” form of the initial perturbation, exists over a finite time interval. An example of such a solution is given.     

Extreme wave phenomena in down-stream running modulated waves

Modulational, Benjamin-Feir, instability is studied for the down-stream evolution of surface gravity waves. An explicit solution, the soliton on finite background, of the NLS equation in physical space is used to study various phenomena in detail. It is shown that for sufficiently long modulation lengths, at a unique position where the largest waves appear, phase singularities are present in the time signal. These singularities are related to wave dislocations and lead to a discrimination between successive ‘extreme’ waves and much smaller intermittent waves. Energy flow in opposite directions through successive dislocations at which waves merge and split, causes the large amplitude difference. The envelope of the time signal at that point is shown to have a simple phase plane representation, and will be described by a symmetry breaking unfolding of the steady state solutions of NLS. The results are used together with the maximal temporal amplitude MTA, to design a strategy for the generation of extreme (freak, rogue) waves in hydrodynamic laboratories.     

Internal waves in an unbounded non-Boussinesq flow

Weakly nonlinear internal waves in an unbounded non-Boussinesq flow with uniform stratification are treated with a Laurent-type expansion. The expansion eliminates the problem encountered with a traditional expansion in wave amplitude where higher harmonics grow exponentially faster with higher order. The results show that the second-order wave correction to the linear estimate of the wave speed of internal waves in an unbounded layer is always negative, meaning that higher amplitude waves travel slower.     

Influence of Finite Amplitude Disturbances on the Nonstationary Modes of a Compressible Boundary Layer Flow

A weakly nonlinear stability analysis is performed to search for the effects of compressibility on a mode of instability of the three-dimensional boundary layer flow due to a rotating disk. The motivation is to extend the stationary work of [ 1 ] (hereafter referred to as S90) to incorporate into the nonstationary mode so that it will be investigated whether the finite amplitude destabilization of the boundary layer is owing to this mode or the mode of S90. Therefore, the basic compressible flow obtained in the large Reynolds number limit is perturbed by disturbances that are nonlinear and also time dependent. In this connection, the effects of nonlinearity are explored allowing the finite amplitude growth of a disturbance close to the neutral location and thus, a finite amplitude equation governing the evolution of the nonlinear lower branch modes is obtained. The coefficients of this evolution equation clearly demonstrate that the nonlinearity is destabilizing for all the modes, the effect of which is higher for the nonstationary waves as compared to the stationary waves. Some modes particularly having positive frequency, regardless of the adiabatic or wall heating/cooling conditions, are always found to be unstable, which are apparently more important than those stationary modes determined in S90. The solution of the asymptotic amplitude equation reveals that compressibility as the local Mach number increases, has the influence of stabilization by requiring smaller initial amplitude of the disturbance for the laminar rotating disk boundary layer flow to become unstable. Apart from the already unstable positive frequency waves, perturbations with positive frequency are always seen to compete to lead the solution to unstable state before the negative frequency waves do. Also, cooling the surface of the disk will be apparently ineffective to suppress the instability mechanisms operating in this boundary layer flow.     

Forced Vibration of a Prestretched Two-Layer Slab on a Rigid Foundation

Within the framework of a piecewise homogeneous body model, with the use of the three-dimensional linearized theory of elastic waves in initially stressed bodies, the forced vibration of a prestretched two-layer slab resting on a rigid foundation is studied. To the upper plane of the slab, a harmonic point force is applied. It is assumed that the layer materials are incompressible, and their elastic properties are characterized by Treloar’s potential. Numerical results are presented for the case where the material stiffness of the lower layer is greater than that of the upper one. The influence of prestretching the layers on the frequency dependences of the normal stresses operating on the interface between the layers and between the slab and the rigid foundation are analyzed.__________Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 41, No. 3, pp. 335–350, May–June, 2005.     

Multiplicity of eigenvalues in certain boundary value problems for the Helmholtz equation

Boundary value problems for two-layer cylindrical waveguides (namely, for a circular two-layer closed waveguide and a dielectric waveguide in an unbounded medium) are considered. It is shown that the same eigenvalues determining the wave numbers of waves in the waveguides can correspond to different solutions of the boundary value problem.     

CDS simulation and pattern formation of phase separation

Several properties of the generation and evolution of phase separating patterns for binary material studied by CDS model are proposed. The main conclusions are (1) for alloys spinodal decomposition, the conceptions of “macro-pattern” and “micropattern” are posed by “black-and- white graph” and “gray-scale graph” respectively. We find that though the four forms of map f that represent the self-evolution of order parameter in a cell (lattice) are similar to each other in “macro-pattern”, there are evident differences in their micro-pattern, e.g., some different fine netted structures in the black domain and the white domain are found by the micro-pattern, so that distinct mechanical and physical behaviors shall be obtained. (2) If the two constitutions of block copolymers are not symmetric (i.e. r ≠ 0.5), a pattern called “grain-strip cross pattern” is discovered, in the 0.43 <r <0.45.     

Boundary-value problems for the heat conduction equation with a fractional derivative in the boundary conditions. Difference methods for numerical realization of these problems

Boundary-value problems for the heat conduction equation are considered in the case where the boundary conditions contain a fractional derivative. Problems of this type arise when the heat processes are simulated by a nonstationary heat flow by using the one-dimensional thermal model of a two-layer system (coating — base). It is proved that the problem under consideration is correct. A one-parameter family of difference schemes is constructed; it is shown that these schemes are stable and convergent in the uniform metric.Translated from Ukrainskii Matematicheskii Zhurnal, Vol. 45, No. 9, pp. 1289–1398, September, 1993.