A Justification of the Modulation Approximation to the 3D Full Water Wave Problem

We consider small amplitude wave packet-like solutions to the 3D inviscid incompressible irrotational infinite depth water wave problem neglecting surface tension. Formal multiscale calculations suggest that the modulation of such a solution is described by a profile traveling at group velocity and governed by a hyperbolic cubic nonlinear Schr\"odinger equation. In this paper we show that, given wave packet initial data, the corresponding solution exists and retains the form of a wave packet on natural NLS time scales. Moreover, we give rigorous error estimates between the true and formal solutions on the appropriate time scale in Sobolev spaces using the energy method. The proof proceeds by directly applying modulational analysis to the formulation of the 3D water wave problem developed by Sijue Wu

Refraction-corrected ray-based inversion for three-dimensional ultrasound tomography of the breast

Ultrasound Tomography has seen a revival of interest in the past decade, especially for breast imaging, due to improvements in both ultrasound and computing hardware. In particular, three-dimensional ultrasound tomography, a fully tomographic method in which the medium to be imaged is surrounded by ultrasound transducers, has become feasible. In this paper, a comprehensive derivation and study of a robust framework for large-scale bent-ray ultrasound tomography in 3D for a hemispherical detector array is presented. Two ray-tracing approaches are derived and compared. More significantly, the problem of linking the rays between emitters and receivers, which is challenging in 3D due to the high number of degrees of freedom for the trajectory of rays, is analysed both as a minimisation and as a root-finding problem. The ray-linking problem is parameterised for a convex detection surface and three robust, accurate, and efficient ray-linking algorithms are formulated and demonstrated. To stabilise these methods, novel adaptive-smoothing approaches are proposed that control the conditioning of the update matrices to ensure accurate linking. The nonlinear UST problem of estimating the sound speed was recast as a series of linearised subproblems, each solved using the above algorithms and within a steepest descent scheme. The whole imaging algorithm was demonstrated to be robust and accurate on realistic data simulated using a full-wave acoustic model and an anatomical breast phantom, and incorporating the errors due to time-of-flight picking that would be present with measured data. This method can used to provide a low-artefact, quantitatively accurate, 3D sound speed maps. In addition to being useful in their own right, such 3D sound speed maps can be used to initialise full-wave inversion methods, or as an input to photoacoustic tomography reconstructions

A Rigorous Finite-Element Domain Decomposition Method for Electromagnetic Near Field Simulations

Rigorous computer simulations of propagating electromagnetic fields have become an important tool for optical metrology and design of nanostructured optical components. A vectorial finite element method (FEM) is a good choice for an accurate modeling of complicated geometrical features. However, from a numerical point of view solving the arising system of linear equations is very demanding even for medium sized 3D domains. In numerics, a domain decomposition method is a commonly used strategy to overcome this problem. Within this approach the overall computational domain is split up into smaller domains and interface conditions are used to assure continuity of the electromagnetic field. Unfortunately, standard implementations of the domain decomposition method as developed for electrostatic problems are not appropriate for wave propagation problems. In an earlier paper we therefore proposed a domain decomposition method adapted to electromagnetic field wave propagation problems. In this paper we apply this method to 3D mask simulation.Comment: 9 pages, 7 figures, SPIE conference Advanced Lithography / Optical Microlithography XXI (2008

An Unstaggered Constrained Transport Method for the 3D Ideal Magnetohydrodynamic Equations

Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with nonlinear numerical instabilities. One approach for controlling the discrete divergence is through a so-called constrained transport method, which is based on first predicting a magnetic field through a standard finite volume solver, and then correcting this field through the appropriate use of a magnetic vector potential. In this work we develop a constrained transport method for the 3D ideal MHD equations that is based on a high-resolution wave propagation scheme. Our proposed scheme is the 3D extension of the 2D scheme developed by Rossmanith [SIAM J. Sci. Comp. 28, 1766 (2006)], and is based on the high-resolution wave propagation method of Langseth and LeVeque [J. Comp. Phys. 165, 126 (2000)]. In particular, in our extension we take great care to maintain the three most important properties of the 2D scheme: (1) all quantities, including all components of the magnetic field and magnetic potential, are treated as cell-centered; (2) we develop a high-resolution wave propagation scheme for evolving the magnetic potential; and (3) we develop a wave limiting approach that is applied during the vector potential evolution, which controls unphysical oscillations in the magnetic field. One of the key numerical difficulties that is novel to 3D is that the transport equation that must be solved for the magnetic vector potential is only weakly hyperbolic. In presenting our numerical algorithm we describe how to numerically handle this problem of weak hyperbolicity, as well as how to choose an appropriate gauge condition. The resulting scheme is applied to several numerical test cases.Comment: 46 pages, 12 figure