Source author record

N. K. -R. Kevlahan

N. K. -R. Kevlahan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

4works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

4 published item(s)

preprint2020arXiv

Convergence analysis of a variational data assimilation scheme for bathymetry detection from surface wave observations

Accurate mapping of ocean bathymetry is a multi-faceted process, needed for safe and efficient navigation on shipping routes and for predicting tsunami waves. Currently available bathymetry data does not always provide the resolution to capture dynamics of such nonlinear waves accurately. However collection of accurate mapping data is difficult, costly, and often a dangerous affair. As an alternative, in this study we implement a variational data assimilation scheme on the one-dimensional shallow water equations to improve estimates of bathymetry, using a finite set of observations of surface wave height to optimise predictions. We show necessary conditions on system parameters for convergence, and implement a low-pass filter for increased regularity of our reconstructed bathymetry. If our objective is to use this to model tsunami propagation, we observe that a relatively higher error in the optimal reconstruction of the bathymetry still yields a highly accurate prediction of the surface wave, suggesting low sensitivity of surface waves to small scale effects in the bathymetry. These conclusions are based on numerical experiments for both Gaussian and sandbar profile bathymetry, and with different observation operators.These extension of these results to realistic models can potentially have a significant impact, as computational cost can be minimised through a priori knowledge of sufficient error tolerances needed for accurate tsunami prediction.

preprint2020arXiv

On the convergence of data assimilation for the one-dimensional shallow water equations with sparse observations

The shallow water equations (SWE) are a widely used model for the propagation of surface waves on the oceans. We consider the problem of optimally determining the initial conditions for the one-dimensional SWE in an unbounded domain from a small set of observations of the sea surface height. In the linear case we prove a theorem that gives sufficient conditions for convergence to the true initial conditions. At least two observation points must be used and at least one pair of observation points must be spaced more closely than half the effective minimum wavelength of the energy spectrum of the initial conditions. This result also applies to the linear wave equation. Our analysis is confirmed by numerical experiments for both the linear and nonlinear SWE data assimilation problems. These results show that convergence rates improve with increasing numbers of observation points and that at least three observation points are required for the practically useful results. Better results are obtained for the nonlinear equations provided more than two observation points are used. This paper is a first step in understanding the conditions for observability of the SWE for small numbers of observation points in more physically realistic settings.

preprint2020arXiv

Second order adjoint sensitivity analysis in variational data assimilation for tsunami models

We mathematically derive the sensitivity of data assimilation results for tsunami modelling, to perturbations in the observation operator. We consider results of variational data assimilation schemes on the one dimensional shallow water equations for (i) initial condition reconstruction, and (ii) bathymetry detection as presented in Kevlahan et al. (2019, 2020). We use variational methods to derive the Hessian of a cost function J representing the error between forecast solutions and observations. Using this Hessian representation and methods outlined by Shutyaev et al. (2017, 2018), we mathematically derive the sensitivity of arbitrary response functions to perturbations in observations for case (i) and (ii) respectively. Such analyses potentially substantiate results from earlier work, on sufficient conditions for convergence, and sensitivity of the propagating surface wave to errors in bathymetry reconstruction. Such sensitivity analyses would illustrate whether particular elements of the observation network are more critical than others, and help minimise extraneous costs for observation collection, and efficiency of predictive models.

preprint2012arXiv

Shock interactions, turbulence, and the origin of the stellar mass spectrum

Supersonic turbulence is an essential element in understanding how structure within interstellar gas is created and shaped. In the context of star formation, many computational studies show that the mass spectrum of density and velocity fluctuations within dense clouds, as well as the distribution of their angular momenta, trace their origin to the statistical and physical properties of gas that is lashed with shock waves. In this article, we review the observations, simulations, and theories of how turbulent-like processes can account for structures we see in molecular clouds. We then compare traditional ideas of supersonic turbulence with a simpler physical model involving the effects of multiple shock waves and their interaction in the interstellar medium. Planar intersecting shock waves produce dense filaments, and generate vortex sheets that are essential to create the broad range of density and velocity structure in clouds. As an example, the lower mass behaviour of the stellar initial mass function can be traced to the tendency of a collection of shock waves to build-up a log-normal density distribution (or column density). Vorticity - which is essential to produce velocity structure over a very broad range of length scales in shocked clouds - can also be generated by the passage of curved shocks or intersecting planar shocks through such media. Two major additional physical forces affect the structure of star forming gas - gravity and feedback processes from young stars. Both of these can produce power-law tails at the high mass end of the IMF.