We propose efficient Eulerian methods for approximating the finite-time Lyapunov exponent (FTLE). The idea is to compute the related flow map using the Level Set Method and the Liouville equation. There are several advantages of the proposed approach. Unlike the usual Lagrangian-type computations, the resulting method requires the velocity field defined only at discrete locations. No interpolation of the velocity field is needed. Also, the method automatically stops a particle trajectory in the case when the ray hits the boundary of the computational domain. The computational complexity of the algorithm is <i>O</i>(<i>x</i>^(<i>d</i>+1)) with <i>d</i> the dimension of the physical space. Since there are the same number of mesh points in the <i>x</i><i>t</i> space, the computational complexity of the proposed Eulerian approach is optimal in the sense that each grid point is visited for only <i>O</i>(1) time. We also extend the algorithm to compute the FTLE

We establish theoretical recovery guarantees of a family of Riemannian optimization algorithms for low rank matrix recovery, which is about recovering an mimes n rank mimes n matrix from mimes n number of linear measurements. The algorithms are first interpreted as iterative hard thresholding algorithms with subspace projections. Based on this connection, we show that provided the restricted isometry constant mimes n of the sensing operator is less than mimes n, the Riemannian gradient descent algorithm and a restarted variant of the Riemannian conjugate gradient algorithm are guaranteed to converge linearly to the underlying rank mimes n matrix if they are initialized by one step hard thresholding. Empirical evaluation shows that the algorithms are able to recover a low rank matrix from nearly the minimum number of measurements necessary.

We develop numerical methods for solving partial differential equations (PDE) defined on an evolving interface represented by the grid based particle method (GBPM) recently proposed in [S. Leung, H.K. Zhao, A grid based particle method for moving interface problems, J. Comput. Phys. 228 (2009) 77067728]. In particular, we develop implicit time discretization methods for the advectiondiffusion equation where the time step is restricted solely by the advection part of the equation. We also generalize the GBPM to solve high order geometrical flows including surface diffusion and Willmore-type flows. The resulting algorithm can be easily implemented since the method is based on meshless particles quasi-uniformly sampled on the interface. Furthermore, without any computational mesh or triangulation defined on the interface, we do not require remeshing or reparametrization in the case of highly distorted motion

Different from conventional spaceborne or airborne synthetic aperture radar (SAR) with optimal aperture length, an imaging radar with highly suboptimal aperture length acquires the data in short bursts by a geometry spreading over a large range. A polarlike or pseudopolar format grid is introduced to sample data close to the resolution, which presents the design of a separable kernel for efficient FFT implementation. The proposed imaging algorithm formulates the reflectivity image of the target scene as an interpolation-free double image series expansion with two usual approximation-induced phase error terms being taken into account, whereby more generalized application scenarios with high frequency, large bandwidth or larger aperture length for imaging a target scene located within either the far-field or the near-field of the radar aperture are processable with high accuracy. In addition, convergence

Laser cladding has been widely used in direct laser fabrication and laser refabrication for metal parts. It is a high precision and complicated process due to multiple process parameters involved. These parameters are strongly tied to each other and their effects are still far from being completely understood. In order to make the optimal choice, a feasible mathematical model is established to describe the interaction between laser and powder particles and to predict and control the cross sectional area. Additionally, analysis concentrating on the effect of the various parameters on the cross sectional area is made. Results show that laser power, powder feed rate, scan speed and radius of focused powder stream are the main factors. Correlations between main process parameters and cross sectional area of an individual clad track have been found. The model also provides predictions of range of process parameters for melting powder to form specific cross sectional area, which will guide the experiment and greatly reduce the cost of the experimental investigations.