An efficient shooting algorithm for Evans function calculations in large systems

In Evans function computations of the spectra of asymptotically constant-coefficient linear operators, a basic issue is the efficient and numerically stable computation of subspaces evolving according to the associated eigenvalue ODE. For small systems, a fast, shooting algorithm may be obtained by representing subspaces as single exterior products \cite{AS,Br.1,Br.2,BrZ,BDG}. For large systems, however, the dimension of the exterior-product space quickly becomes prohibitive, growing as $\binom{n}{k}$, where $n$ is the dimension of the system written as a first-order ODE and $k$ (typically $\sim n/2$) is the dimension of the subspace. We resolve this difficulty by the introduction of a simple polar coordinate algorithm representing ``pure'' (monomial) products as scalar multiples of orthonormal bases, for which the angular equation is a numerically optimized version of the continuous orthogonalization method of Drury--Davey \cite{Da,Dr} and the radial equation is evaluable by quadrature. Notably, the polar-coordinate method preserves the important property of analyticity with respect to parameters.Comment: 21 pp., two figure

Grassmannian flows and applications to nonlinear partial differential equations

We show how solutions to a large class of partial differential equations with nonlocal Riccati-type nonlinearities can be generated from the corresponding linearized equations, from arbitrary initial data. It is well known that evolutionary matrix Riccati equations can be generated by projecting linear evolutionary flows on a Stiefel manifold onto a coordinate chart of the underlying Grassmann manifold. Our method relies on extending this idea to the infinite dimensional case. The key is an integral equation analogous to the Marchenko equation in integrable systems, that represents the coodinate chart map. We show explicitly how to generate such solutions to scalar partial differential equations of arbitrary order with nonlocal quadratic nonlinearities using our approach. We provide numerical simulations that demonstrate the generation of solutions to Fisher--Kolmogorov--Petrovskii--Piskunov equations with nonlocal nonlinearities. We also indicate how the method might extend to more general classes of nonlinear partial differential systems.Comment: 26 pages, 2 figure

Discrete mechanics and optimal control: An analysis

The optimal control of a mechanical system is of crucial importance in many application areas. Typical examples are the determination of a time-minimal path in vehicle dynamics, a minimal energy trajectory in space mission design, or optimal motion sequences in robotics and biomechanics. In most cases, some sort of discretization of the original, infinite-dimensional optimization problem has to be performed in order to make the problem amenable to computations. The approach proposed in this paper is to directly discretize the variational description of the system's motion. The resulting optimization algorithm lets the discrete solution directly inherit characteristic structural properties from the continuous one like symmetries and integrals of the motion. We show that the DMOC (Discrete Mechanics and Optimal Control) approach is equivalent to a finite difference discretization of Hamilton's equations by a symplectic partitioned Runge-Kutta scheme and employ this fact in order to give a proof of convergence. The numerical performance of DMOC and its relationship to other existing optimal control methods are investigated

Lie symmetry analysis and numerical solutions for thermo-solutal chemicallyreacting radiative micropolar flow from an inclined porous surface

Steady, laminar, incompressible thermo-solutal natural convection flow of micropolar fluid from an inclined perforated surface with convective boundary conditions is studied. Thermal radiative flux and chemical reaction effects are included to represent phenomena encountered in high-temperature materials synthesis operations. Rosseland’s diffusion approximation is used to describe the radiative heat flux in the energy equation. A Lie scaling group transformation is implemented to derive a self-similar form of the partial differential conservation equations. The resulting coupled nonlinear boundary value problem is solved with Runge-Kutta fourth order numerical quadrature (shooting technique). Validation of solutions with an optimized Adomian decomposition method algorithm is included. Verification of the accuracy of shooting is also conducted as a particular case of non-reactive micropolar flow from a vertical permeable surface. The evolution of velocity, angular velocity (micro-rotation component), temperature and concentration are examined for a variety of parameters including coupling number, plate inclination angle, suction/injection parameter, radiation-conduction parameter, Biot number and reaction parameter. Numerical results for steady state skin friction coefficient, couple stress coefficient, Nusselt number and Sherwood number are tabulated and discussed. Interesting features of the hydrodynamic, heat and mass transfer characteristics are examined

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