Numerical modelling of two-way propagation of non-linear dispersive waves

We describe and implement a fully discrete spectral method for the numerical solution of a class of non-linear, dispersive systems of Boussinesq type, modelling two-way propagation of long water waves of small amplitude in a channel. For three particular systems, we investigate properties of the numerically computed solutions; in particular we study the generation and interaction of approximate solitary waves.

Propagation of waves generated by a pressure disturbance moving in a channel

This paper studies the effect of ship speed and water depth on the propagation of ship generated waves. The ship is represented by a moving pressure distribution function at the free surface that is able to reproduce most of the phenomena involved in wave propagation. Results are obtained for a ship sailing along a coastal stretch made of a sloping bottom and a constant depth region. The results show that in the sloping bottom the crests of waves are bent along the slope and in the constant depth the standard Kelvin wave patterns can be found for the subcritical regime. In the critical regime the wave system is characterized by significant diverging waves and for a supercritical regime, the transverse waves disappear. © 2015 Taylor & Francis Group, London.

Equatorial waves in the lower stratosphere. II: Annual and interannual variability

The annual and interannual variability of idealized, linear, equatorial waves in the lower stratosphere is investigated using the temperature and velocity fields from the ECMWF 15-year re-analysis dataset. Peak Kelvin wave activity occurs during solstice seasons at 100 hPa, during December-February at 70 hPa and in the easterly to westerly quasi-biennial oscillation (QBO) phase transition at 50 hPa. Peak Rossby-gravity wave activity occurs during equinox seasons at 100 hPa, during June-August/September-November at 70 hPa and in the westerly to easterly QBO phase transition at 50 hPa. Although neglect of wind shear means that the results for inertio-gravity waves are likely to be less accurate, they are still qualitatively reasonable and an annual cycle is observed in these waves at 100 hPa and 70 hPa. Inertio-gravity waves with n = 1 are correlated with the QBO at 50 hPa, but the eastward inertio-gravity n = 0 wave is not, due to its very fast vertical group velocity in all background winds. The relative importance of different wave types in driving the QBO at 50 hPa is also discussed. The strongest acceleration appears to be provided by the Kelvin wave while the acceleration provided by the Rossby-gravity wave is negligible. Of the higher-frequency waves, the westward inertio-gravity n = 1 wave appears able to contribute more to the acceleration of the 50 hPa mean zonal wind than the eastward inertio-gravity n = 1 wave.

Transmission problem for waves with frictional damping April 22, 2007

In this paper we consider the transmission problem, in one space dimension, for linear dissipative waves with frictional damping. We study the wave propagation in a medium with a component with attrition and another simply elastic. We show that for this type of material, the dissipation produced by the frictional part is strong enough to produce exponential decay of the solution, no matter how small is its size. ©2007 Texas State University.

Comment on ``Piecewise potential vorticity inversion: Elementary tests''

Egger (2008) constructs some idealised experiments to test the usefulness of piecewise potential vorticity inversion (PPVI) in the diagnosis of Rossby wave dynamics and baroclinic development. He concludes that, ``PPVI does not help us to understand the dynamics of linear Rossby waves. It provides local tendencies of the streamfunction which are unrelated to the true ones. The same way, the motion of baroclinic waves in shear flow cannot be understood by using PPVI. Moreover, the effect of boundary temperatures as determined by PPVI is unrelated to the flow evolution.'' He goes further in arguing that we should not consider velocities as ``induced'' by PV anomalies defined by carving up the global domain. However, these conclusions partly reflect the limitations of his idealised experiments and the manner in which the PV components were partitioned from one another.

Planetary-wave-induced transport in the stratosphere

In the stratosphere, chemical tracers are drawn systematically from the equator to the pole. This observed Brewer–Dobson circulation is driven by wave drag, which in the stratosphere arises mainly from the breaking and dissipation of planetary-scale Rossby waves. While the overall sense of the circulation follows from fundamental physical principles, a quantitative theoretical understanding of the connection between wave drag and Lagrangian transport is limited to linear, small-amplitude waves. However, planetary waves in the stratosphere generally grow to a large amplitude and break in a strongly nonlinear fashion. This paper addresses the connection between stratospheric wave drag and Lagrangian transport in the presence of strong nonlinearity, using a mechanistic three-dimensional primitive equations model together with offline particle advection. Attention is deliberately focused on a weak forcing regime, such that sudden warmings do not occur and a quasi-steady state is reached, in order to examine this question in the cleanest possible context. Wave drag is directly linked to the transformed Eulerian mean (TEM) circulation, which is often used as a surrogate for mean Lagrangian motion. The results show that the correspondence between the TEM and mean Lagrangian velocities is quantitatively excellent in regions of linear, nonbreaking waves (i.e., outside the surf zone), where streamlines are not closed. Within the surf zone, where streamlines are closed and meridional particle displacements are large, the agreement between the vertical components of the two velocity fields is still remarkably good, especially wherever particle paths are coherent so that diabatic dispersion is minimized. However, in this region the meridional mean Lagrangian velocity bears little relation to the meridional TEM velocity, and reflects more the kinematics of mixing within and across the edges of the surf zone. The results from the mechanistic model are compared with those from the Canadian Middle Atmosphere Model to test the robustness of the conclusions.

Rossby wave propagation on potential vorticity fronts with finite width

The horizontal gradient of potential vorticity (PV) across the tropopause typically declines with lead time in global numerical weather forecasts and tends towards a steady value dependent on model resolution. This paper examines how spreading the tropopause PV contrast over a broader frontal zone affects the propagation of Rossby waves. The approach taken is to analyse Rossby waves on a PV front of finite width in a simple single-layer model. The dispersion relation for linear Rossby waves on a PV front of infinitesimal width is well known; here an approximate correction is derived for the case of a finite width front, valid in the limit that the front is narrow compared to the zonal wavelength. Broadening the front causes a decrease in both the jet speed and the ability of waves to propagate upstream. The contribution of these changes to Rossby wave phase speeds cancel at leading order. At second order the decrease in jet speed dominates, meaning phase speeds are slower on broader PV fronts. This asymptotic phase speed result is shown to hold for a wide class of single-layer dynamics with a varying range of PV inversion operators. The phase speed dependence on frontal width is verified by numerical simulations and also shown to be robust at finite wave amplitude, and estimates are made for the error in Rossby wave propagation speeds due to the PV gradient error present in numerical weather forecast models.

Análise dos resultados do tratamento da mordida cruzada posterior funcional com o expansor fixo Hyrax

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

High resolution numerical simulation of triple point collision and origin of unburned gas pockets in turbulent detonations

A key issue in pulse detonation engine development is better understanding of the detonation structure and its propagation mechanism. Thus, in the present work the turbulent structure of an irregular detonation is studied through very high resolution numerical simulations of 600 points per half reaction length. The aim is to explore the nature of the transverse waves during the collision and reflection processes of the triple point with the channel walls. Consequently the formation and consumption mechanism of unreacted gas pockets is studied. Results show that the triple point and the transverse wave collide simultaneously with the wall. The strong transverse wave switches from a primary triple point before collision to a new one after reflection. Due to simultaneous interaction of the triple point and the transverse wave with the wall in the second half of the detonation cell, a larger high-pressurised region appears on the wall. During the reflection the reaction zone detaches from the shock front and produces a pocket of unburned gas. Three mechanisms found to be of significance in the re-initiation mechanism of detonation at the end of the detonation cell; i: energy resealed via consumption of unburned pockets by turbulent mixing ii: compression waves arise due to collision of the triple point on the wall which helps the shock to jump abruptly to an overdriven detonation iii: drastic growth of the Richtmyer–Meshkov instability causing a part of the front to accelerate with respect to the neighbouring portions.

GLOBAL WELL-POSEDNESS AND NON-LINEAR STABILITY OF PERIODIC TRAVELING WAVES FOR A SCHRODINGER-BENJAMIN-ONO SYSTEM

The objective of this paper is two-fold: firstly, we develop a local and global (in time) well-posedness theory for a system describing the motion of two fluids with different densities under capillary-gravity waves in a deep water flow (namely, a Schrodinger-Benjamin-Ono system) for low-regularity initial data in both periodic and continuous cases; secondly, a family of new periodic traveling waves for the Schrodinger-Benjamin-Ono system is given: by fixing a minimal period we obtain, via the implicit function theorem, a smooth branch of periodic solutions bifurcating a Jacobian elliptic function called dnoidal, and, moreover, we prove that all these periodic traveling waves are nonlinearly stable by perturbations with the same wavelength.

Desarrollo de Algoritmos Numéricos para el Estudio de la Propagación no Lineal de Ondas en Aplicaciones Aeroespaciales y Astrofísicas. Develop to Numerical Algorithms to Study non-Linear Waves Propagation in Aerospace and Astrophysics Applications.

En este proyecto se desarrollarán algoritmos numéricos para sistemas no lineales hiperbólicos-parabólicos de ecuaciones diferenciales en derivadas parciales. Dichos sistemas tienen aplicación en propagación de ondas en ámbitos aeroespaciales y astrofísicos.Objetivos generales: 1)Desarrollo y mejora de algoritmos numéricos con la finalidad de incrementar la calidad en la simulación de propagación e interacción de ondas gasdinámicas y magnetogasdinámicas no lineales. 2)Desarrollo de códigos computacionales con la finalidad de simular flujos gasdinámicos de elevada entalpía incluyendo cambios químicos, efectos dispersivos y difusivos.3)Desarrollo de códigos computacionales con la finalidad de simular flujos magnetogasdinámicos ideales y reales.4)Aplicación de los nuevos algoritmos y códigos computacionales a la solución del flujo aerotermodinámico alrededor de cuerpos que ingresan en la atmósfera terrestre. 5)Aplicación de los nuevos algoritmos y códigos computacionales a la simulación del comportamiento dinámico no lineal de arcos magnéticos en la corona solar. 6)Desarrollo de nuevos modelos para describir el comportamiento no lineal de arcos magnéticos en la corona solar.Este proyecto presenta como objetivo principal la introducción de mejoras en algoritmos numéricos para simular la propagación e interacción de ondas no lineales en dos medios gaseosos: aquellos que no poseen carga eléctrica libre (flujos gasdinámicos) y aquellos que tienen carga eléctrica libre (flujos magnetogasdinámicos). Al mismo tiempo se desarrollarán códigos computacionales que implementen las mejoras de las técnicas numéricas.Los algoritmos numéricos se aplicarán con la finalidad de incrementar el conocimiento en tópicos de interés en la ingeniería aeroespacial como es el cálculo del flujo de calor y fuerzas aerotermodinámicas que soportan objetos que ingresan a la atmósfera terrestre y en temas de astrofísica como la propagación e interacción de ondas, tanto para la transferencia de energía como para la generación de inestabilidades en arcos magnéticos de la corona solar. Estos dos temas poseen en común las técnicas y algoritmos numéricos con los que serán tratados. Las ecuaciones gasdinámicas y magnetogasdinámicas ideales conforman sistemas hiperbólicos de ecuaciones diferenciales y pueden ser solucionados utilizando "Riemann solvers" junto con el método de volúmenes finitos (Toro 1999; Udrea 1999; LeVeque 1992 y 2005). La inclusión de efectos difusivos genera que los sistemas de ecuaciones resulten hiperbólicos-parabólicos. La contribución parabólica puede ser considerada como términos fuentes y tratada adicionalmente tanto en forma explícita como implícita (Udrea 1999; LeVeque 2005).Para analizar el flujo alrededor de cuerpos que ingresan en la atmósfera se utilizarán las ecuaciones de Navier-Stokes químicamente activas, mientras la temperatura no supere los 6000K. Para mayores temperaturas es necesario considerar efectos de ionización (Anderson, 1989). Tanto los efectos difusivos como los cambios químicos serán considerados como términos fuentes en las ecuaciones de Euler. Para tratar la propagación de ondas, transferencia de energía e inestabilidades en arcos magnéticos de la corona solar se utilizarán las ecuaciones de la magnetogasdinámica ideal y real. En este caso será también conveniente implementar términos fuente para el tratamiento de fenómenos de transporte como el flujo de calor y el de radiación. Los códigos utilizarán la técnica de volúmenes finitos, junto con esquemas "Total Variation Disminishing - TVD" sobre mallas estructuradas y no estructuradas.

Generation and nonlinear evolution of shore-oblique/transverse sand bars

The coupling between topography, waves and currents in the surf zone may selforganize to produce the formation of shore-transverse or shore-oblique sand bars on an otherwise alongshore uniform beach. In the absence of shore-parallel bars, this has been shown by previous studies of linear stability analysis, but is now extended to the finite-amplitude regime. To this end, a nonlinear model coupling wave transformation and breaking, a shallow-water equations solver, sediment transport and bed updating is developed. The sediment flux consists of a stirring factor multiplied by the depthaveraged current plus a downslope correction. It is found that the cross-shore profile of the ratio of stirring factor to water depth together with the wave incidence angle primarily determine the shape and the type of bars, either transverse or oblique to the shore. In the latter case, they can open an acute angle against the current (upcurrent oriented) or with the current (down-current oriented). At the initial stages of development, both the intensity of the instability which is responsible for the formation of the bars and the damping due to downslope transport grow at a similar rate with bar amplitude, the former being somewhat stronger. As bars keep on growing, their finite-amplitude shape either enhances downslope transport or weakens the instability mechanism so that an equilibrium between both opposing tendencies occurs, leading to a final saturated amplitude. The overall shape of the saturated bars in plan view is similar to that of the small-amplitude ones. However, the final spacings may be up to a factor of 2 larger and final celerities can also be about a factor of 2 smaller or larger. In the case of alongshore migrating bars, the asymmetry of the longshore sections, the lee being steeper than the stoss, is well reproduced. Complex dynamics with merging and splitting of individual bars sometimes occur. Finally, in the case of shore-normal incidence the rip currents in the troughs between the bars are jet-like while the onshore return flow is wider and weaker as is observed in nature.