Analysis of the Nonlinear Shallow Water Equations Over Nonplanar Topography

The role played by the beach bottom profile on coastal inundation phenomena is analyzed here by means of approximate analytical solutions of the nonlinear shallow water equations (NSWEs) over uneven bottoms. These are obtained by only using the assumptions of small waves at the seaward boundary and small topographic forcing. Our work, built on the Carrier and Greenspan [ 1 ] hodographic transformation and on the solution of the boundary value problem (BVP) for the NSWEs proposed by Antuono and Brocchini [ 2 ], focuses on the propagation of nonlinear non-breaking waves over quasi-planar beaches. Since the terms associated with the perturbed bottom only appear in the second-order perturbed solutions, the breaking conditions for the planar-beach bathymetry also predict well the breaking occurring on the nonplanar beaches analyzed here. The most important results, concerning the shoreline position and the near-shoreline velocity, are given for both pulse-like and periodic input waves propagating over two types of nonplanar bathymetries. The solution proposed here is a fundamental benchmark for any numerical and theoretical analyzes concerned with estimates of wave run-up on beaches of complex shape.     

Generating Boundary Conditions for the Calculation of Tsunami Propagation on Nested Grids

Some boundary conditions used to numerically simulate tsunami generation and propagation are studied. Special attention is given to generating boundary conditions thatmake it possible to simulate tsunami waves with desired characteristics (amplitude, time period and, in general, waveform). Since the water flow velocity in a propagating tsunami wave is uniquely defined by its height and ocean depth, one can simulate a wave propagating from the boundary into the simulation area. This can be done by specifying the wave height and water flow velocity on the boundary. This method is used to numerically simulate the propagation of a tsunami from the source to the coast on a sequence of refined grids. In this numerical experiment the wave parameters are transferred from the larger area to the subarea via boundary conditions. This method can also generate a wave that has certain characteristics on a specified line.     

A simple and robust boundary treatment for the forced Korteweg–de Vries equation

In this paper, we propose a simple and robust numerical method for the forced Korteweg–de Vries (fKdV) equation which models free surface waves of an incompressible and inviscid fluid flow over a bump. The fKdV equation is defined in an infinite domain. However, to solve the equation numerically we must truncate the infinite domain to a bounded domain by introducing an artificial boundary and imposing boundary conditions there. Due to unsuitable artificial boundary conditions, most wave propagation problems have numerical difficulties (e.g., the truncated computational domain must be large enough or the numerical simulation must be terminated before the wave approaches the artificial boundary for the quality of the numerical solution). To solve this boundary problem, we develop an absorbing non-reflecting boundary treatment which uses outward wave velocity. The basic idea of the proposing algorithm is that we first calculate an outward wave velocity from the solutions at the previous and present time steps and then we obtain a solution at the next time step on the artificial boundary by moving the solution at the present time step with the velocity. And then we update solutions at the next time step inside the domain using the calculated solution on the artificial boundary. Numerical experiments with various initial conditions for the KdV and fKdV equations are presented to illustrate the accuracy and efficiency of our method.     

Open boundary conditions for hyperbolic equations

A previously developed general procedure for deriving accurate difference equations to describe conditions at open boundaries for hyperbolic equations is extended and further illustrated by means of several examples of practical importance. Problems include those with both incoming and outgoing waves at the boundary, the use of locally cylindrical and spherical wave approximations at each point of the boundary, and nonlinear wave propagation. Reflected waves in all cases are minimal and less than 10^?2 of the incident wave.     

Uniform asymptotics of the boundary values of the solution in a linear problem on the run-up of waves on a shallow beach

We consider the Cauchy problem with spatially localized initial data for a two-dimensional wave equation with variable velocity in a domain Ω. The velocity is assumed to degenerate on the boundary ?Ω of the domain as the square root of the distance to ?Ω. In particular, this problems describes the run-up of tsunami waves on a shallow beach in the linear approximation. Further, the problem contains a natural small parameter (the typical source-to-basin size ratio) and hence admits analysis by asymptotic methods. It was shown in the paper “Characteristics with singularities and the boundary values of the asymptotic solution of the Cauchy problem for a degenerate wave equation” [1] that the boundary values of the asymptotic solution of this problem given by a modified Maslov canonical operator on the Lagrangian manifold formed by the nonstandard characteristics associatedwith the problemcan be expressed via the canonical operator on a Lagrangian submanifold of the cotangent bundle of the boundary. However, the problem as to how this restriction is related to the boundary values of the exact solution of the problem remained open. In the present paper, we show that if the initial perturbation is specified by a function rapidly decaying at infinity, then the restriction of such an asymptotic solution to the boundary gives the asymptotics of the boundary values of the exact solution in the uniform norm. To this end, we in particular prove a trace theorem for nonstandard Sobolev type spaces with degeneration at the boundary.     

Interacting Wave Fronts and Rarefaction Waves in a Second Order Model of Nonlinear Thermoviscous Fluids

A wave equation including nonlinear terms up to the second order for a thermoviscous Newtonian fluid is proposed. In the lossless case this equation results from an expansion to third order of the Lagrangian for the fundamental non-dissipative fluid dynamical equations. Thus it preserves the Hamiltonian structure, in contrast to the Kuznetsov equation, a model often used in nonlinear acoustics. An exact traveling wave front solution is derived from a generalized traveling wave assumption for the velocity potential. Numerical studies of the evolution of a number of arbitrary initial conditions as well as head-on colliding and confluent wave fronts exhibit several nonlinear interaction phenomena. These include wave fronts of changed velocity and amplitude along with the emergence of rarefaction waves. An analysis using the continuity of the solutions as well as the boundary conditions is proposed. The dynamics of the rarefaction wave is approximated by a collective coordinate approach in the energy balance equation.     

Long distance wave computation using nonlinear solitary waves

A recently developed method is described to propagate short wave equation pulses over indefinite distances and through regions of varying indices of refraction, including multiple reflections. The method, “Wave Confinement”, utilizes a newly developed nonlinear partial differential equation (pde) that propagates basis functions according to the wave equation. These basis functions are generated as stable solitary waves where the discretized equation can be solved without any numerical dissipation. The method can also be used to solve for harmonic waves in the high frequency (Eikonal) limit, including multiple arrivals. The solution involves discretizing the wave equation on a uniform Eulerian grid and adding a simple nonlinear “Confinement” term. This term does not change the amplitude (integrated through each point on the pulse surface) or the propagation velocity, or arrival time, and yet results in capturing the waves as thin surfaces that propagate as thin nonlinear solitary waves and remain ∼2-3 grid cells in thickness indefinitely with no numerical spreading. A new feature described in this paper involves computing scattering of short pulses from complex objects such as complete aircraft. A simple “immersed surface” approach is used, that utilizes the same uniform grid as the propagation and avoids complex, body fitted or adaptive grid schemes.The new method should be useful in areas of wave propagation, from radar scattering and long distance communications to cell phone transmission.     

Existence and uniqueness of solutions to a nonlinear elliptic boundary value problem on time scales

In this paper, we consider a boundary value problem (BVP) for second-order nonlinear partial dynamic equations on the time scale rectangles. Some explicit conditions are established that ensure existence and uniqueness of solution to the BVP under consideration.     

Two-Wave Interactions for Hyperbolic Signaling Problems and the Application to Supersonic Gas Dynamics

The signaling problem for a system of conservation laws in a single space variable is treated through the deployment of a perturbation analysis. Our method of approach involves the direct use of two nonlinear phase variables making possible the study of weakly nonlinear interacting waves arising from a boundary disturbance consisting of two wave modes. As a result of our analysis, the asymptotic solution is derived, and the class of admissible boundary disturbances is distinguished as well. An application is then made to gas dynamics in one space dimension to investigate the propagation and interaction of two sound waves for which the base state is taken to be a steady supersonic flow.     

On the interaction of deep water waves and exponential shear currents

A train of periodic deep-water waves propagating on a steady shear current with a vertical distribution of vorticity is investigated by an analytic method, namely the homotopy analysis method (HAM). The magnitude of the vorticity varies exponentially with the magnitude of the stream function, while remaining constant on a particular streamline. The so-called Dubreil–Jacotin transformation is used to transfer the original exponentially nonlinear boundary-value problem in an unknown domain into an algebraically nonlinear boundary-value problem in a known domain. Convergent series solutions are obtained not only for small amplitude water waves on a weak current but also for large amplitude waves on a strong current. The nonlinear wave-current interaction is studied in detail. It is found that an aiding shear current tends to enlarge the wave phase speed, sharpen the wave crest, but shorten the maximum wave height, while an opposing shear current has the opposite effect. Besides, the amplitude of waves and fluid velocity decay over the depth more quickly on an aiding shear current but more slowly on an opposing shear current than that of waves on still water. Furthermore, it is found that Stokes criteria of wave breaking is still valid for waves on a shear current: a train of propagating waves on a shear current breaks as the fiuid velocity at crest equals the wave phase speed. Especially, it is found that the highest waves on an opposing shear current are even higher and steeper than that of waves on still water. Mathematically, this analytic method is rather general in principle and can be employed to solve many types of nonlinear partial differential equations with variable coefficients in science, finance and engineering.     

On the interaction of deep water waves and exponential shear currents

A train of periodic deep-water waves propagating on a steady shear current with a vertical distribution of vorticity is investigated by an analytic method, namely the homotopy analysis method (HAM). The magnitude of the vorticity varies exponentially with the magnitude of the stream function, while remaining constant on a particular streamline. The so-called Dubreil–Jacotin transformation is used to transfer the original exponentially nonlinear boundary-value problem in an unknown domain into an algebraically nonlinear boundary-value problem in a known domain. Convergent series solutions are obtained not only for small amplitude water waves on a weak current but also for large amplitude waves on a strong current. The nonlinear wave-current interaction is studied in detail. It is found that an aiding shear current tends to enlarge the wave phase speed, sharpen the wave crest, but shorten the maximum wave height, while an opposing shear current has the opposite effect. Besides, the amplitude of waves and fluid velocity decay over the depth more quickly on an aiding shear current but more slowly on an opposing shear current than that of waves on still water. Furthermore, it is found that Stokes criteria of wave breaking is still valid for waves on a shear current: a train of propagating waves on a shear current breaks as the fiuid velocity at crest equals the wave phase speed. Especially, it is found that the highest waves on an opposing shear current are even higher and steeper than that of waves on still water. Mathematically, this analytic method is rather general in principle and can be employed to solve many types of nonlinear partial differential equations with variable coefficients in science, finance and engineering.     

Construction of traveling wave solutions of the (2 + 1)-dimensional modified KdV-KP equation

Traveling wave solutions have played a vital role in demonstrating the wave character of nonlinear problems emerging in the field of mathematical sciences and engineering. To depict the nature of propagation of the nonlinear waves in nature, a range of nonlinear evolution equations has been proposed and investigated in the existing literature. In this article, solitary and traveling periodic wave solutions for the (2 + 1)-dimensional modified KdV-KP equation are derived by employing an ansatz method, named the enhanced (G′/G)-expansion method. For this continued equation, abundant solitary wave solutions and nonlinear periodic wave solutions, along with some free parameters, are obtained. We have derived the exact expressions for the solitary waves that arise in the continuum-modified KdV-KP model. We study the significance of parameters numerically that arise in the obtained solutions. These parameters play an important role in the physical structure and propagation directions of the wave that characterizes the wave pattern. We discuss the relation between velocity and parameters and illustrate them graphically. Our numerical analysis suggests that the taller solitons are narrower than shorter waves and can travel faster. In addition, graphical representations of some obtained solutions along with their contour plot and wave train profiles are presented. The speed, as well as the profile of these solitary waves, is highly sensitive to the free parameters. Our results establish that the continuum-modified KdV-KP system supports solitary waves having different shapes and speeds for different values of the parameters.     

Well-Posedness of a Problem of Propagation of Nonlinear Acoustic-Gravity Waves in the Atmosphere Generated by Surface Pressure Variations

Currently there are many international microbarograph networks for high-resolution recording of wave pressure variations on the Earth’s surface. This arouses interest in wave propagation in the atmosphere generated by atmospheric pressure variations. A full system of nonlinear hydrodynamic equations for atmospheric gases with lower boundary conditions in the form of wavelike pressure variations on the Earth’s surface is considered. Since the wave amplitudes near the Earth’s surface are small, linearized equations are used in the analysis of well-posedness of the problem. With the help of a wave energy functional method, it is shown that in the non-dissipative case the solution to the boundary value problem is uniquely determined by the variable pressure field on the Earth’s surface. The corresponding dissipative problem is well-posed if, in addition to the pressure field, appropriate conditions on the velocity and temperature on the Earth’s surface are given. In the case of an isothermal atmosphere, the problem admits analytical solutions that are harmonic in the variables x and t. A good agreement between the numerical and analytical solutions is obtained. The study shows that the temperature and density can rapidly vary at the lower boundary of the boundary value problem. An example of solving the three-dimensional problem with variable pressure on the Earth’s surface taken from experimental observations is given. The developed algorithms and computer programs can be used to simulate atmospheric waves generated by pressure variations on the Earth’s surface.     

A Boundary Value Problem for Second Order Nonlinear Difference Equations on the Semi-infinite Interval

In this paper, we consider a boundary value problem (BVP) for nonlinear difference equations on the discrete semi-axis in which the left-hand side being a second order linear difference expression belongs to the so-called Weyl-Hamburger limit-circle case. The BVP is considered in the Hilbert space l 2 and is formed via boundary conditions at a starting point and at infinity. Existence and uniqueness results for solutions of the considered BVP are established.     

Rayleigh type wave dispersion in an incompressible functionally graded orthotropic half-space loaded by a thin fluid-saturated aeolotropic porous layer

This paper is concerned with the Rayleigh wave dispersion in an incompressible functionally graded orthotropic half-space loaded by a thin fluid-saturated aeolotropic porous layer under initial stress. Both the layer and half-space have subjected to the incompressible in nature. The particle motion of the Rayleigh type wave is elliptically polarized in the plane, which described by the normal to the surface and the focal point along with wave generation. The dispersion of waves refers typically to frequency dispersion, which means different wavelengths travel at a different velocity of phase. To deal with the analytical solution of displacement components of Rayleigh type waves in a layer over a half-space, we have taken the assistance of different methods like exponential, characteristic polynomial and undetermined coefficients. The dispersion relation has been derived based upon suitable boundary conditions. The finite difference scheme has been introduced to calculate the phase velocity and group velocity of the Rayleigh type waves. We also have derived the stability condition of the finite difference scheme (FDS) for the phase and group velocities. If a wave equation has to travel in the time domain, it is necessary to achieve both accuracy and stability requirements. In such cases, FDS is preferred because of its power, accuracy, reliability, rapidity, and flexibility. The effect of various parameters involved in the model like non-homogeneity, porosity, and internal pre-stress on the propagation of Rayleigh type waves have been studied in detail. Graphical representations for the effects of various parameters on the dispersion equation have been represented. Numerical results demonstrated the accuracy and versatility of the group and phase velocity depending on the stability ratio of the FDS.     

On the Analytic Solutions of a Special Boundary Value Problem for a Nonlinear Heat Equation in Polar Coordinates

The paper addresses a nonlinear heat equation (the porous medium equation) in the case of a power-law dependence of the heat conductivity coefficient on temperature. The equation is used for describing high-temperature processes, filtration of gases and fluids, groundwater infiltration, migration of biological populations, etc. The heat waves (waves of filtration) with a finite velocity of propagation over a cold background form an important class of solutions to the equation under consideration. A special boundary value problem having solutions of such type is studied. The boundary condition of the problem is given on a sufficiently smooth closed curve with variable geometry. The new theorem of existence and uniqueness of the analytic solution is proved.     

Propagation of Surface Waves at High Frequencies

An asymptotic theory is presented for the analysis of surfacewave propagation at high frequencies. The theory is developedfor scalar surface waves satisfying an impedance boundary conditionon a surface, which may be curved and, whose impedance may bevariable. A surface eikonal equation is derived for the phaseof the surface wave field, and it is shown that the wave fieldpropagates over the surface along the surface rays, which arethe characteristics of the surface eikonal equation. The wavefield in space is found by solving certain eikonal and transportequations with the aid of complex rays. The theory is then appliedto several examples: axial waves on a circular cylinder, sphericallysymmetric waves on a sphere, waves on a circular cone with avariable impedance, and waves on the plane boundary of an inhomogeneousmedium. In each case it is found that the asymptotic expansionof the exact solution agrees with the asymptotic solution.     

Love waves in a fluid-saturated porous layer under a rigid boundary and lying over an elastic half-space under gravity

In this paper, mathematical modeling of the propagation of Love waves in a fluid-saturated porous layer under a rigid boundary and lying over an elastic half-space under gravity has been considered. The equations of motion have been formulated separately for different media under suitable boundary conditions at the interface of porous layer, elastic half-space under gravity and rigid layer. Following Biot, the frequency equation has been derived which contain Whittaker’s function and its derivative that have been expanded asymptotically up to second term (for approximate result) for large argument due to small values of Biot’s gravity parameter (varying from 0 to 1). The effect of porosity and gravity of the layers in the propagation of Love waves has been studied. The effect of hydrostatic initial stress generated due to gravity in the half-space has also been shown in the phase velocity of Love waves. The phase velocity of Love waves for first two modes has been presented graphically. Frequency equations have also been derived for some particular cases, which are in perfect agreement with standard results. Subsequently the lower and upper bounds of Love wave speed have also been discussed.     

Boundary value problems in the theory of binary mixtures of thermoelastic solids

In this paper the diffusion model of the linear theory of thermoelasticity of binary mixtures is considered. The basic properties of wave numbers of the longitudinal and transverse plane waves are established. The interior boundary value problem (BVP) of steady vibration of binary mixtures of thermoelastic solids is formulated. The existence theorems of eigenfrequencies of the BVP of steady vibrations are proved. The connection between plane waves and existence of eigenfrequencies is established. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)