旋转地球流体的一个非线性奇异孤立波精确解

本文求出了旋转地球流体的f平面浅水波方程组的一个非线性精确解.这是一个处处连续、分片光滑的凹形孤立波,在波谷处为一尖点（奇异点）,其移速比长波速度gh^1/2略小,为（1+a/h）（gh（1+a/h））^1/2,其中a<0.其水平尺度与水深几乎无关.     

Runup of Tsunami Waves in U-Shaped Bays

The problem of tsunami wave shoaling and runup in U-shaped bays (such as fjords) and underwater canyons is studied in the framework of 1D shallow water theory with the use of an assumption of the uniform current on the cross-section. The wave shoaling in bays, when the depth varies smoothly along the channel axis, is studied with the use of asymptotic approach. In this case a weak reflection provides significant shoaling effects. The existence of traveling (progressive) waves, propagating in bays, when the water depth changes significantly along the channel axis, is studied within rigorous solutions of the shallow water theory. It is shown that traveling waves do exist for certain bay bathymetry configurations and may propagate over large distances without reflection. The tsunami runup in such bays is significantly larger than for a plane beach.     

Convection in rotating annular channels heated from below. Part 3: Experimental boundary conditions

Results are presented from both linear stability analysis and numerical simulations of three-dimensional nonlinear convection in a Boussinesq fluid in an annular channel, under experimental boundary conditions, rotating about a vertical axis uniformly heated from below. The focus is placed on the Prandtl number Pr = 7.0, representing liquid water at room temperature. The linear analysis shows that, when the aspect ratio is sufficiently small, there exists only one stationary mode that occupies the whole fluid container. When the aspect ratio is moderate or large, however, there exist three different linear solutions: (i) the outer sidewall-localized traveling wave propagating against the sense of rotation; (ii) the inner sidewall-localized traveling wave propagating in the same sense as rotation; and (iii) both the counter-traveling waves occurring simultaneously. Guided by the result of the linear stability analysis, fully three-dimensional simulations are then performed for a channel with a moderate aspect ratio. It is found that neither the prograde nor the retrograde mode is physically realizable near threshold and beyond. The dynamics of nonlinear convection in a rotating channel are chiefly characterized by the interaction between the sidewall-localized waves and the interior convection cells/rolls, producing an interesting and unusual nonlinear phenomenon. In order to compare with the classical Rayleigh–Bénard problem without vertical sidewalls, we also study linear and nonlinear convection at exactly the same parameters but in an infinitely extended layer with periodic horizontal conditions. This reveals that both the linear instability and nonlinear convection in a rotating channel are characteristically different from those in a rotating layer with periodic horizontal conditions.     

Discontinuous Steady‐State Analytical Solutions of the Boussinesq Equation and Their Numerical Representation by Modflow

Known analytical solutions of groundwater flow equations are routinely used for verification of computer codes. However, these analytical solutions (e.g., the Dupuit solution for the steady‐state unconfined unidirectional flow in a uniform aquifer with a flat bottom) represent smooth and continuous water table configurations, simulating which does not pose any significant problems for the numerical groundwater flow models, like MODFLOW. One of the most challenging numerical cases for MODFLOW arises from drying‐rewetting problems often associated with abrupt changes in the elevations of impervious base of a thin unconfined aquifer. Numerical solutions of groundwater flow equations cannot be rigorously verified for such cases due to the lack of corresponding exact analytical solutions. Analytical solutions of the steady‐state Boussinesq equation, associated with the discontinuous water table configurations over a stairway impervious base, are presented in this article. Conditions resulting in such configurations are analyzed and discussed. These solutions appear to be well suited for testing and verification of computer codes. Numerical solutions, obtained by the latest versions of MODFLOW (MODFLOW‐2005 and MODFLOW‐NWT), are compared with the presented discontinuous analytical solutions. It is shown that standard MODFLOW‐2005 code (as well as MODFLOW‐2000 and older versions) has significant convergence problems simulating such cases. The problems manifest themselves either in a total convergence failure or erroneous results. Alternatively, MODFLOW‐NWT, providing a good match to the presented discontinuous analytical solutions, appears to be a more reliable and appropriate code for simulating abrupt changes in water table elevations.     

雅可比椭圆函数在大气和海洋动力学中的应用:二维非线性Rossby波研究

提出了扩展雅可比椭圆函数方法,来求得Petviashvili方程的精确解析解.Petviashvili方程被视为正压准地转位涡度方程的非地转扩展,应用该方法可以得到很多二维非线性Rossby波的周期波解,在取极限情况下,也可以得到二维Rossby孤立子解.     

High-order finite volume WENO schemes for the shallow water equations with dry states

The shallow water equations are used to model flows in rivers and coastal areas, and have wide applications in ocean, hydraulic engineering, and atmospheric modeling. These equations have still water steady state solutions in which the flux gradients are balanced by the source term. It is desirable to develop numerical methods which preserve exactly these steady state solutions. Another main difficulty usually arising from the simulation of dam breaks and flood waves flows is the appearance of dry areas where no water is present. If no special attention is paid, standard numerical methods may fail near dry/wet front and produce non-physical negative water height. A high-order accurate finite volume weighted essentially non-oscillatory (WENO) scheme is proposed in this paper to address these difficulties and to provide an efficient and robust method for solving the shallow water equations. A simple, easy-to-implement positivity-preserving limiter is introduced. One- and two-dimensional numerical examples are provided to verify the positivity-preserving property, well-balanced property, high-order accuracy, and good resolution for smooth and discontinuous solutions.     

Weakly nonlinear shallow water magnetohydrodynamic waves

We consider an electrically conducting rotating fluid governed by the shallow water magnetohydrodynamic equations with no diffusion. We use an a priori asymptotic technique (the method of geometric optics or ray method) to study weakly nonlinear hydromagnetic waves. These waves are intermediate in length in the following sense: they are much longer than the fluid depth but much shorter than the radius of the earth. The time scale for the waves is much longer than that of the free surface oscillations and the approximation varies on an even longer timescale. The waves we are considering are studied in the beta plane approximation for an ambient magnetic field parallel to the equator which varies in the direction perpendicular to the equator. The leading order approximation gives a dispersion relation for the waves, which are generally found to be confined to bands about the equator as well as in bands at higher and lower latitudes. At the next order of approximation, a conservation law is found for the wave amplitude. We also obtain an equation governing the behavior of the leading order mean azimuthal velocity which is forced to grow linearly with time.     

Fission of a weakly nonlinear interfacial solitary wave at a step

The transformation of a weakly nonlinear interfacial solitary wave in an ideal two-layer flow over a step is studied. In the vicinity of the step the wave transformation is described in the framework of the linear theory of long interfacial waves, and the coefficients of wave reflection and transmission are calculated. A strong transformation arises for propagation into shallower water, but a weak transformation for propagation into deeper water. Far from the step, the wave dynamics is described by the Korteweg-de Vries equation which is fully integrable. In the vicinity of the step, the reflected and transmitted waves have soliton-like shapes, but their parameters do not satisfy the steady-state soliton solutions. Using the inverse scattering technique it is shown that the reflected wave evolves into a single soliton and dispersing radiation if the wave propagates from deep to shallow water, and only dispersing radiation if the wave propagates from shallow to deep water. The dynamics of the transmitted wave is more complicated. In particular, if the coefficient of the nonlinear quadratic term in the Korteweg-de Vries equation is not changed in sign in the region after the step, the transmitted wave evolves into a group of solitons and radiation, a process similar to soliton fission for surface gravity waves at a step. But if the coefficient of the nonlinear term changes sign, the soliton is destroyed completely and transforms into radiation. The effects of cubic nonlinearity are studied in the framework of the extended Korteweg-de Vries (Gardner) equation which is also integrable. The higher-order nonlinear effects influence the amplitudes of the generated solitons if the amplitude of the transformed wave is comparable with the thickness of lower layer, but otherwise the process of soliton fission is qualitatively the same as in the framework of the Korteweg-de Vries equation.     

Linear Magnetoconvection in Rotating Fluid Spheres Permeated by a Uniform Axial Magnetic Field

A linear analysis of thermally driven magnetoconvection is carried out with emphasis on its application to convection in the Earth's core. We consider a rotating and self-gravitating fluid sphere (or spherical shell) permeated by a uniform magnetic field parallel to the spin axis. In rapidly rotating cases, we find that five different convective modes appear as the uniform field is increased; namely, geostrophic, polar convective, magneto-geostrophic, fast magnetostrophic and slow magnetostrophic modes. The polar convective (P) and magneto-geostrophic (E) modes seem to be of geophysical interest. The P mode is characterized by such an axisymmetric meridional circulation that the fluid penetrates the equatorial plane, suggesting that generation of quadrapole from dipole fields could be explained by a linear process. The E mode is characterized by a few axially aligned columnar rolls which are almost two-dimensional due to a modified Proudman-Taylor theorem.     

切变基本纬向流中非线性赤道Rossby长波

为了解决观测和理论研究中的一些问题以及更好地了解热带大气动力学 ,有必要进一步研究基本气流的变化对大气中赤道Rossby波动的影响 .本文研究分析基本气流对赤道Rossby长波的影响 ,利用一个简单赤道 β平面浅水模式和摄动法 ,研究纬向基本气流切变中非线性赤道Rossby波 ,推导出在切变基本纬向流中赤道Rossby长波振幅演变所满足的非线性KdV方程并得到其孤立波解 .分析表明 ,孤立波存在的必要条件是基本气流有切变 ,而且基流切变不能太强 ,否则将产生正压不稳定 .     

切变基本纬向流中非线性赤道Rossby长波

为了解决观测和理论研究中的一些问题以及更好地了解热带大气动力学 ,有必要进一步研究基本气流的变化对大气中赤道Rossby波动的影响 .本文研究分析基本气流对赤道Rossby长波的影响 ,利用一个简单赤道 β平面浅水模式和摄动法 ,研究纬向基本气流切变中非线性赤道Rossby波 ,推导出在切变基本纬向流中赤道Rossby长波振幅演变所满足的非线性KdV方程并得到其孤立波解 .分析表明 ,孤立波存在的必要条件是基本气流有切变 ,而且基流切变不能太强 ,否则将产生正压不稳定 .     

EVOLUTION EQUATIONS FOR NONLINEAR WAVES IN TWO-LAYER FLUID SYSTEM OF WATER AND SOFT MUD

LINTRODUCTIONSoftmudhaVebeenfoundinmanycoasts,rivers,andeStheriesallovertheworld.ThemutUaleffectSofpropagatingwavesandsoftmudseedbedshavebeendiscussedbymostinveshgatorsbasedonlinearwavetheory.Buttheinteractionbetweenwaterwaveandsoftmudisnonlinear.OneaPProachisthePe~ationmethod,usedbyPeregrine(1967).AlternativeformofBoussunesqeqUationsisderiVedbyusingthevelocityatanarbitrarydistancefromthestillwaterlevelasthevelocityvariableinsteadofthecommonlyuseddepth-averagedvelocitybyNWogU(1993).…     

Statistical distribution of nonlinear random wave height in shallow water

Here we present a statistical model of random wave,using Stokes wave theory of water wave dynamics,as well as a new nonlinear probability distribution function of wave height in shallow water.It is more physically logical to use the wave steepness of shallow water and the factor of shallow water as the parameters in the wave height distribution.The results indicate that the two parameters not only could be parameters of the distribution function of wave height but also could reflect the degree of wave heigh...     

Validation and Verification of a Numerical Model for Tsunami Propagation and Runup

A robust numerical model to simulate propagation and runup of tsunami waves in the framework of non-linear shallow water theory is developed. The numerical code adopts a staggered leapfrog finite-difference scheme to solve the shallow water equations formulated for depth-averaged water fluxes in spherical coordinates. A temporal position of the shoreline is calculated using a free-surface moving boundary algorithm. For large scale problems, the developed algorithm is efficiently parallelized employing a domain decomposition technique. The developed numerical model is benchmarked in an exhaustive series of tests suggested by NOAA. We conducted analytical and laboratory benchmarking for the cases of solitary wave runup on simple beaches, runup of a solitary wave on a conically-shaped island, and the runup in the Monai Valley, Okushiri Island, Japan, during the 1993 Hokkaido-Nansei-Oki tsunami. In all conducted tests the calculated numerical solution is within an accuracy recommended by NOAA standards. We summarize results of numerical benchmarking of the model, its strengths and limits with regards to reproduction of fundamental features of coastal inundation, and also illustrate some possible improvements.     

A reduced model for nonlinear dispersive waves in a rotating environment

Abstract

The simplest model for geophysical flows is one layer of a constant density fluid with a free surface, where the fluid motions occur on a scale in which the Coriolis force is significant. In the linear shallow water limit, there are non-dispersive Kelvin waves, localized near a boundary or near the equator, and a large family of dispersive waves. We study weakly nonlinear and finite depth corrections to these waves, and derive a reduced system of equations governing the flow. For this system we find approximate solitary Kelvin waves, both for waves traveling along a boundary and along the equator. These waves induce jets perpendicular to their direction of propagation, which may have a role in mixing. We also derive an equivalent reduced system for the evolution of perturbations to a mean geostrophic flow.     

The instability of precessing flow

Abstract

An explanation is put forward for the instability observed within a precessing, rotating spheroidal container. The constant vorticity solution for the flow suggested by Poincaré is found to be inertially unstable through the parametric coupling of two inertial waves by the underlying constant strain field. Such resonant couplings are due either to the elliptical or shearing strains present which elliptically distort the circular streamlines and shear their centres respectively. For the precessing Earth's outer core, the shearing of the streamlines and the ensuing shearing instability are the dominant features. The instability of some exact, linear solutions for finite precessional rates is established and used to corroborate the asymptotic analysis. A complementary unbounded analysis of a precessing, rotating fluid is also presented and used to deduce a likely upperbound on the growth rate of a small disturbance. Connection is made with past experimental studies.     

Numerical modeling of 3D fully nonlinear potential periodic waves

A simple and exact numerical scheme for long-term simulations of 3D potential fully nonlinear periodic gravity waves is suggested. The scheme is based on the surface-following nonorthogonal curvilinear coordinate system. Velocity potential is represented as a sum of analytical and nonlinear components. The Poisson equation for the nonlinear component of velocity potential is solved iteratively. Fourier transform method, the second-order accuracy approximation of vertical derivatives on a stretched vertical grid and the fourth-order Runge–Kutta time stepping are used. The scheme is validated by simulation of steep Stokes waves. A one-processor version of the model for PC allows us to simulate evolution of a wave field with thousands degrees of freedom for hundreds of wave periods. The scheme is designed for investigation of nonlinear 2D surface waves, generation of extreme waves, and direct calculations of nonlinear interactions.     

Space discontinuous Galerkin method for shallow water flows—kinetic and HLLC flux,and potential vorticity generation

In this paper, a second order space discontinuous Galerkin (DG) method is presented for the numerical solution of inviscid shallow water flows over varying bottom topography. Novel in the implementation is the use of HLLC and kinetic numerical fluxes¹ in combination with a dissipation operator, applied only locally around discontinuities to limit spurious numerical oscillations. Numerical solutions over (non-)uniform meshes are verified against exact solutions; the numerical error in the L₂-norm and the convergence of the solution are computed. Bore–vortex interactions are studied analytically and numerically to validate the model; these include bores as “breaking waves” in a channel and a bore traveling over a conical and Gaussian hump. In these complex numerical test cases, we correctly predict the generation of potential vorticity by non-uniform bores. Finally, we successfully validate the numerical model against measurements of steady oblique hydraulic jumps in a channel with a contraction. In the latter case, the kinetic flux is shown to be more robust.     

Baroclinic,Kelvin and inertia-gravity waves in the barostrat instability experiment

The differentially heated rotating annulus is a laboratory experiment historically designed for modelling large-scale features of the mid-latitude atmosphere. In the present study, we investigate a modified version of the classic baroclinic experiment in which a juxtaposition of convective and motionless stratified layers is created by introducing a vertical salt stratification. The thermal convective motions are suppressed in a central region at mid-depth of the rotating tank, therefore double-diffusive convection rolls can develop only in thin layers located at top and bottom, where the salt stratification is weakest. For high enough rotation rates, the baroclinic instability destabilises the flow in the top and the bottom shallow convective layers, generating cyclonic and anticyclonic eddies separated by the stable stratified layer. Thanks to this alternation of layers resembling the convective and radiative layers of stars, the planetary’s atmospheric troposphere and stratosphere or turbulent layers at the sea surface above stratified waters, this new laboratory setup is of interest for both astrophysics and geophysical sciences. More specifically, it allows to study the exchange of momentum and energy between the layers, primarily by the propagation of internal gravity waves (IGW). PIV velocity maps are used to describe the wavy flow pattern at different heights. Using a co-rotating laser and camera, the wave field is well resolved and different wave types can be found: baroclinic waves, Kelvin and Poincaré type waves. The signature of small-scale IGW can also be observed attached to the baroclinic jet. The baroclinic waves occur at the thin convectively active layer at the surface and the bottom of the tank, though decoupled they show different manifestation of nonlinear interactions. The inertial Kelvin and Poincaré waves seem to be mechanically forced. The small-scale wave trains attached to the meandering jet point to an imbalance of the large-scale flow. For the first time, the simultaneous occurrence of different wave types is reported in detail for a differentially heated rotating annulus experiment.     

Magnetic effects of sea tides

The tidal motion of sea water across the earth's magnetic field is known to induce small electric currents to flow in the oceans and the surrounding land. There has been recently a great deal of interest in this phenomenon. Here we consider a model ocean in the form of an infinitely long channel with a tidal wave passing along it. It is shown that for the lunar tidal frequency of 12.45 hours the solution of this problem for typical ocean depths exhibits characteristics which are essentially due to the high frequency. The solution therefore differs appreciably from low-frequency solutions which are therefore only applicable as models of the shallow seas. In our model we observe concentrations of electric current near the coast and phase changes relative to the oceanic tide. A new “equivalent field exclusion principle” is presented and used to explain some of the results we have obtained. The method used is simple and applies readily to waves of any frequency, and results are given for waves of higher frequencies. The amplitude and phase lag of the magnetic field of a channel simulating the Atlantic Ocean are discussed.