Expedition in Data and Harmonic Analysis on Graphs

The graph Laplacian operator is widely studied in spectral graph theory largely due to its importance in modern data analysis. Recently, the Fourier transform and other time-frequency operators have been defined on graphs using Laplacian eigenvalues and eigenvectors. We extend these results and prove that the translation operator to the i’th node is invertible if and only if all eigenvectors are nonzero on the i’th node. Because of this dependency on the support of eigenvectors we study the characteristic set of Laplacian eigenvectors. We prove that the Fiedler vector of a planar graph cannot vanish on large neighborhoods and then explicitly construct a family of non-planar graphs that do exhibit this property. We then prove original results in modern analysis on graphs. We extend results on spectral graph wavelets to create vertex-dyanamic spectral graph wavelets whose support depends on both scale and translation parameters. We prove that Spielman’s Twice-Ramanujan graph sparsifying algorithm cannot outperform his conjectured optimal sparsification constant. Finally, we present numerical results on graph conditioning, in which edges of a graph are rescaled to best approximate the complete graph and reduce average commute time.

POSITIVE FILTERED P_N METHOD FOR LINEAR TRANSPORT EQUATIONS AND THE ASSOCIATED OPTIMIZATION ALGORITHM

We propose a positive, accurate moment closure for linear kinetic transport equations based on a filtered spherical harmonic (FP_N) expansion in the angular variable. The FP_N moment equations are accurate approximations to linear kinetic equations, but they are known to suffer from the occurrence of unphysical, negative particle concentrations. The new positive filtered P_N (FP_N+) closure is developed to address this issue. The FP_N+ closure approximates the kinetic distribution by a spherical harmonic expansion that is non-negative on a finite, predetermined set of quadrature points. With an appropriate numerical PDE solver, the FP_N+ closure generates particle concentrations that are guaranteed to be non-negative. Under an additional, mild regularity assumption, we prove that as the moment order tends to infinity, the FP_N+ approximation converges, in the L2 sense, at the same rate as the FP_N approximation; numerical tests suggest that this assumption may not be necessary. By numerical experiments on the challenging line source benchmark problem, we confirm that the FP_N+ method indeed produces accurate and non-negative solutions. To apply the FP_N+ closure on problems at large temporal-spatial scales, we develop a positive asymptotic preserving (AP) numerical PDE solver. We prove that the propose AP scheme maintains stability and accuracy with standard mesh sizes at large temporal-spatial scales, while, for generic numerical schemes, excessive refinements on temporal-spatial meshes are required. We also show that the proposed scheme preserves positivity of the particle concentration, under some time step restriction. Numerical results confirm that the proposed AP scheme is capable for solving linear transport equations at large temporal-spatial scales, for which a generic scheme could fail. Constrained optimization problems are involved in the formulation of the FP_N+ closure to enforce non-negativity of the FP_N+ approximation on the set of quadrature points. These optimization problems can be written as strictly convex quadratic programs (CQPs) with a large number of inequality constraints. To efficiently solve the CQPs, we propose a constraint-reduced variant of a Mehrotra-predictor-corrector algorithm, with a novel constraint selection rule. We prove that, under appropriate assumptions, the proposed optimization algorithm converges globally to the solution at a locally q-quadratic rate. We test the algorithm on randomly generated problems, and the numerical results indicate that the combination of the proposed algorithm and the constraint selection rule outperforms other compared constraint-reduced algorithms, especially for problems with many more inequality constraints than variables.

Effects of Drive Speed Modulation on Dynamics of Slender Rotating Structures

Slender rotating structures are used in many mechanical systems. These structures can suffer from undesired vibrations that can affect the components and safety of a system. Furthermore, since some these structures can operate in a harsh environment, installation and operation of sensors that are needed for closed-loop and collocated control schemes may not be feasible. Hence, the need for an open-loop non-collocated scheme for control of the dynamics of these structures. In this work, the effects of drive speed modulation on the dynamics of slender rotating structures are studied. Slender rotating structures are a type of mechanical rotating structures, whose length to diameter ratio is large. For these structures, the torsion mode natural frequencies can be low. In particular, for isotropic structures, the first few torsion mode frequencies can be of the same order as the first few bending mode frequencies. These situations can be conducive for energy transfer amongst bending and torsion modes. Scenarios with torsional vibrations experienced by rotating structures with continuous rotor-stator contact occur in many rotating mechanical systems. Drill strings used in the oil and gas industry are an example of rotating structures whose torsional vibrations can be deleterious to the components of the drilling system. As a novel approach to mitigate undesired vibrations, the effects of adding a sinusoidal excitation to the rotation speed of a drill string are studied. A portion of the drill string located within a borewell is considered and this rotating structure has been modeled as an extended Jeffcott rotor and a sinusoidal excitation has been added to the drive speed of the rotor. After constructing a three-degree-of-freedom model to capture lateral and torsional motions, the equations of motions are reduced to a single differential equation governing torsional vibrations during continuous stator contact. An approximate solution has been obtained by making use of the Method of Direct Partition of Motions with the governing torsional equation of motion. The results showed that for a rotor undergoing forward or backward whirling, the addition of sinusoidal excitation to the drive speed can cause an increase in the equivalent torsional stiffness, smooth the discontinuous friction force at contact, and reduce the regions of negative slope in the friction coefficient variation with respect to speed. Experiments with a scaled drill string apparatus have also been conducted and the experimental results show good agreement with the numerical results obtained from the developed models. These findings suggest that the extended Jeffcott rotordynamics model can be useful for studies of rotor dynamics in situations with continuous rotor-stator contact. Furthermore, the results obtained suggest that the drive speed modulation scheme can have value for attenuating drill-string vibrations.