Schemes for Smooth Discretization And Inverse Problems - Case Study on Recovery of Tsunami Source Parameters

This thesis deals with smooth discretization schemes and inverse problems, the former used in efficient yet accurate numerical solutions to forward models required in turn to solve inverse problems. The aims of the thesis include, (i) development of a stabilization techniques for a class of forward problems plagued by unphysical oscillations in the response due to the presence of jumps/shocks/high gradients, (ii) development of a smooth hybrid discretization scheme that combines certain useful features of Finite Element (FE) and Mesh-Free (MF) methods and alleviates certain destabilizing factors encountered in the construction of shape functions using the polynomial reproduction method and, (iii) a first of its kind attempt at the joint inversion of both static and dynamic source parameters of the 2004 Sumatra-Andaman earthquake using tsunami sea level anomaly data. Following the introduction in Chapter 1 that motivates and puts in perspective the work done in later chapters, the main body of the thesis may be viewed as having two parts, viz., the first part constituting the development and use of smooth discretization schemes in the possible presence of destabilizing factors (Chapters 2 and 3) and the second part involving solution to the inverse problem of tsunami source recovery (Chapter 4). In the context of stability requirements in numerical solutions of practical forward problems, Chapter 2 develops a new stabilization scheme. It is based on a stochastic representation of the discretized field variables, with a view to reduce or even eliminate unphysical oscillations in the MF numerical simulations of systems developing shocks or exhibiting localized bands of extreme plastic deformation in the response. The origin of the stabilization scheme may be traced to nonlinear stochastic filtering and, consistent with a class of such filters, gain-based additive correction terms are applied to the simulated solution of the system, herein achieved through the Element-Free Galerkin (EFG) method, in order to impose a set of constraints that help arresting the spurious oscillations. The method is numerically illustrated through its application to a gradient plasticity model whose response is often characterized by a developing shear band as the external load is gradually increased. The potential of the method in stabilized yet accurate numerical simulations of such systems involving extreme gradient variations in the response is thus brought forth. Chapter 3 develops the MF-based discretization motif by balancing this with the widespread adoption of the FE method. Thus it concentrates on developing a 'hybrid' scheme that aims at the amelioration of certain destabilizing algorithmic issues arising from the necessary condition of moment matrix invertibility en route to the generation of smooth shape functions. It sets forth the hybrid discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing approach adopted over a conventional FE-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity ( C p 1 ) and local supports of the simplex splines of degree p . In the proposed scheme, the triangles comprising the domain discretization also serve as background cells for numerical integration which here are near-aligned to the supports of the shape functions (and their intersections), thus considerably ameliorating an oft-cited source of inaccuracy in the numerical integration of MF-based weak forms. Numerical experiments establish that the proposed method can work with lower order quadrature rules for accurate evaluation of integrals in the Galerkin weak form, a feature desiderated in solving nonlinear inverse problems that demand cost-effective solvers for the forward models. Numerical demonstrations of optimal convergence rates for a few test cases are given and the hybrid method is also implemented to compute crack-tip fields in a gradient-enhanced elasticity model. Chapter 4 attempts at the joint inversion of earthquake source parameters for the 2004 Sumatra-Andaman event from the tsunami sea level anomaly signals available from satellite altimetry. Usual inversion for earthquake source parameters incorporates subjective elements, e.g. a priori constraints, posing and parameterization, trial-and-error waveform fitting etc. Noisy and possibly insufficient data leads to stability and non-uniqueness issues in common deterministic inversions. A rational accounting of both issues favours a stochastic framework which is employed here, leading naturally to a quantification of the commonly overlooked aspects of uncertainty in the solution. Confluence of some features endows the satellite altimetry for the 2004 Sumatra-Andaman tsunami event with unprecedented value for the inversion of source parameters for the entire rupture duration. A nonlinear joint inversion of the slips, rupture velocities and rise times with minimal a priori constraints is undertaken. Large and hitherto unreported variances in the parameters despite a persistently good waveform fit suggest large propagation of uncertainties and hence the pressing need for better physical models to account for the defect dynamics and massive sediment piles. Chapter 5 concludes the work with pertinent comments on the results obtained and suggestions for future exploration of some of the schemes developed here

Prospects of large scale onshore marine pearl culture along the Indian coasts

CMFRI attempted the development of a technology for cultured marine pearls in onshore tanks like any other pond systems under controlled conditions. The results of the experiment, first of its kind, are presented and discussed here

Overview of the marine fisheries research in the Lakshadweep

The Union Territory of Lakshadweep consists of 10 inhabited islands and 16 uninhabited islets. The U.T. has an area of 32 sq km and possesses 400,000 sq km of Exclusive Economic Zone. Each island except Androth has a lagoon on the western side, and the lagoon and the reef porvide an ideal coral habitat for a variety of flora and fauna

Non-continuous Froude number scaling for the closure depth of a cylindrical cavity

A long, smooth cylinder is dragged through a water surface to create a cavity with an initially cylindrical shape. This surface void then collapses due to the hydrostatic pressure, leading to a rapid and axisymmetric pinch-off in a single point. Surprisingly, the depth at which this pinch-off takes place does not follow the expected Froude$^{1/3}$ power-law. Instead, it displays two distinct scaling regimes separated by discrete jumps, both in experiment and in numerical simulations (employing a boundary integral code). We quantitatively explain the above behavior as a capillary waves effect. These waves are created when the top of the cylinder passes the water surface. Our work thus gives further evidence for the non-universality of the void collapse

A general reduction method for one-loop N-point integrals

In order to calculate cross sections with a large number of particles/jets in the final state at next-to-leading order, one has to reduce the occurring scalar and tensor one-loop integrals to a small set of known integrals. In massless theories, this reduction procedure is complicated by the presence of infrared divergences. Working in n=4-2*epsilon dimensions, it will be outlined how to achieve such a reduction for diagrams with an arbitrary number of external legs. As a result, any integral with more than four propagators and generic 4-dimensional external momenta can be reduced to box integrals.Comment: 5 pages Latex, 1 eps figure included, uses npb.sty (included). Talk presented at the conference: Loops and Legs in Quantum Field Theory, April 2000, Bastei, German

Reduction formalism for dimensionally regulated one-loop N-point integrals

We consider one-loop scalar and tensor integrals with an arbitrary number of external legs relevant for multi-parton processes in massless theories. We present a procedure to reduce N-point scalar functions with generic 4-dimensional external momenta to box integrals in (4-2\epsilon) dimensions. We derive a formula valid for arbitrary N and give an explicit expression for N=6. Further a tensor reduction method for N-point tensor integrals is presented. We prove that generically higher dimensional integrals contribute only to order \epsilon for N>=5. The tensor reduction can be solved iteratively such that any tensor integral is expressible in terms of scalar integrals. Explicit formulas are given up to N=6.Comment: 21 pages Latex, 1 eps figur

Status, prospects and management of small pelagic fisheries in India

The annual small pelagic fish production increased from 0.30 million mt during 1950-54 to 1.24 million mt during 1996 along the Indian coast. The 4 fold increase was possible due to several technological advancements. The potenrial yield from the pelagic resources of the EEZ is estimated to be 2.2 million mt. As there is no further scope for increasing the production from the inshore waters, there is need to bring the outer shelf and oceanic waters into increasing levels of exploitation

Sudden Collapse of a Granular Cluster

Single clusters in a vibro-fluidized granular gas in N connected compartments become unstable at strong shaking. They are experimentally shown to collapse very abruptly. The observed cluster lifetime (as a function of the driving intensity) is analytically calculated within a flux model, making use of the self-similarity of the process. After collapse, the cluster diffuses out into the uniform distribution in a self-similar way, with an anomalous diffusion exponent 1/3.Comment: 4 pages, 4 figures. Figure quality has been reduced in order to decrease file-siz