Unsteady undular bores in fully nonlinear shallow-water theory

We consider unsteady undular bores for a pair of coupled equations of Boussinesq-type which contain the familiar fully nonlinear dissipationless shallow-water dynamics and the leading-order fully nonlinear dispersive terms. This system contains one horizontal space dimension and time and can be systematically derived from the full Euler equations for irrotational flows with a free surface using a standard long-wave asymptotic expansion. In this context the system was first derived by Su and Gardner. It coincides with the one-dimensional flat-bottom reduction of the Green-Naghdi system and, additionally, has recently found a number of fluid dynamics applications other than the present context of shallow-water gravity waves. We then use the Whitham modulation theory for a one-phase periodic travelling wave to obtain an asymptotic analytical description of an undular bore in the Su-Gardner system for a full range of "depth" ratios across the bore. The positions of the leading and trailing edges of the undular bore and the amplitude of the leading solitary wave of the bore are found as functions of this "depth ratio". The formation of a partial undular bore with a rapidly-varying finite-amplitude trailing wave front is predicted for ``depth ratios'' across the bore exceeding 1.43. The analytical results from the modulation theory are shown to be in excellent agreement with full numerical solutions for the development of an undular bore in the Su-Gardner system.Comment: Revised version accepted for publication in Phys. Fluids, 51 pages, 9 figure

Flood propagation modelling with the Local Inertia Approximation: theoretical and numerical analysis of its physical limitations

Attention of the researchers has increased towards a simplification of the complete Shallow water Equations called the Local Inertia Approximation (LInA), which is obtained by neglecting the advection term in the momentum conservation equation. In the present paper it is demonstrated that a shock is always developed at moving wetting-drying frontiers, and this justifies the study of the Riemann problem on even and uneven beds. In particular, the general exact solution for the Riemann problem on horizontal frictionless bed is given, together with the exact solution of the non-breaking wave propagating on horizontal bed with friction, while some example solution is given for the Riemann problem on discontinuous bed. From this analysis, it follows that drying of the wet bed is forbidden in the LInA model, and that there are initial conditions for which the Riemann problem has no solution on smoothly varying bed. In addition, propagation of the flood on discontinuous sloping bed is impossible if the bed drops height have the same order of magnitude of the moving-frontier shock height. Finally, it is found that the conservation of the mechanical energy is violated. It is evident that all these findings pose a severe limit to the application of the model. The numerical analysis has proven that LInA numerical models may produce numerical solutions, which are unreliable because of mere algorithmic nature, also in the case that the LInA mathematical solutions do not exist. The applicability limits of the LInA model are discouragingly severe, even if the bed elevation varies continuously. More important, the non-existence of the LInA solution in the case of discontinuous topography and the non-existence of receding fronts radically question the viability of the LInA model in realistic cases. It is evident that classic SWE models should be preferred in the majority of the practical applications

Optical diffraction of focused spots and subwavelength structures

We have developed a numerical diffraction tool for cases in which the incident field is a focused spot and the diffracting structure is a single structure or an aperiodic surface. Our approach uses the integral formulation to solve Maxwell’s equations and is different from previously published methods in its choice of basis function. We compared numerical results with experimental measurements of the far-field intensity for a focused spot incident on an aluminum grating, and the comparison was favorable. Finally, we predict the diffraction behavior of the proposed digital video disk format for the next generation of optical disk. Our analysis shows that the reflected signal for this format has a strong dependence on the polarization of the incident light

On the power of homogeneous depth 4 arithmetic circuits

We prove exponential lower bounds on the size of homogeneous depth 4 arithmetic circuits computing an explicit polynomial in $VP$. Our results hold for the {\it Iterated Matrix Multiplication} polynomial - in particular we show that any homogeneous depth 4 circuit computing the $(1,1)$ entry in the product of $n$ generic matrices of dimension $n^{O(1)}$ must have size $n^{\Omega(\sqrt{n})}$. Our results strengthen previous works in two significant ways. Our lower bounds hold for a polynomial in $VP$. Prior to our work, Kayal et al [KLSS14] proved an exponential lower bound for homogeneous depth 4 circuits (over fields of characteristic zero) computing a poly in $VNP$. The best known lower bounds for a depth 4 homogeneous circuit computing a poly in $VP$ was the bound of $n^{\Omega(\log n)}$ by [LSS, KLSS14].Our exponential lower bounds also give the first exponential separation between general arithmetic circuits and homogeneous depth 4 arithmetic circuits. In particular they imply that the depth reduction results of Koiran [Koi12] and Tavenas [Tav13] are tight even for reductions to general homogeneous depth 4 circuits (without the restriction of bounded bottom fanin). Our lower bound holds over all fields. The lower bound of [KLSS14] worked only over fields of characteristic zero. Prior to our work, the best lower bound for homogeneous depth 4 circuits over fields of positive characteristic was $n^{\Omega(\log n)}$ [LSS, KLSS14]