No abstract uploaded!

We introduce a new adaptive method for analyzing nonlinear and nonstationary data. This method is inspired by the empirical mode decomposition (EMD) method and the recently developed compressed sensing theory. The main idea is to look for the sparsest representation of multiscale data within the largest possible dictionary consisting of intrinsic mode functions of the form {a(t) cos(θ(t))}, where a ≥ 0 is assumed to be smoother than cos(θ(t)) and θ is a piecewise smooth increasing function. We formulate this problem as a nonlinear L1 optimization problem. Further, we propose an iterative algorithm to solve this nonlinear optimization problem recursively. We also introduce an adaptive filter method to decompose data with noise. Numerical examples are given to demonstrate the robustness of our method and comparison is made with the EMD method. One advantage of performing such a decomposition is to preserve some intrinsic physical property of the signal, such as trend and instantaneous frequency. Our method shares many important properties of the original EMD method. Because our method is based on a solid mathematical formulation, its performance does not depend on numerical parameters such as the number of shifting or stop criterion, which seem to have a major effect on the original EMD method. Our method is also less sensitive to noise perturbation and the end effect compared with the original EMD method.

In this article, we present a second-order in time implicit-explicit (IMEX) local discontinuous Galerkin (LDG) method for computing the Cahn-Hilliard equation, which describes the phase separation phenomenon. It is well-known that the Cahn-Hilliard equation has a nonlinear stability property, i.e., the free-energy functional decreases with respect to time. The discretized Cahn-Hilliard system modeled by the IMEX LDG method can inherit the nonlinear stability of the continuous model. We apply a stabilization technique and prove unconditional energy stability of our scheme. Numerical experiments are performed to validate the analysis. Computational efficiency can be significantly enhanced by using this IMEX LDG method with a large time step.

For cable bridges, the cable tension force plays a crucial role in their construction, assessment and long-term structural health monitoring. Cable tension forces vary in real time with the change of the moving vehicle loads and environmental effects, and this continual variation in tension force may cause fatigue damage of a cable. Traditional vibration-based cable tension force estimation methods can only obtain the time-averaged cable tension force and not the instantaneous force. This paper proposes a new approach to identify the time-varying cable tension forces of bridges based on an adaptive sparse time-frequency analysis method. This is a recently developed method to estimate the instantaneous frequency by looking for the sparsest time-frequency representation of the signal within the largest possible time-frequency dictionary (i.e. set of expansion functions). In the proposed approach, first, the time-varying modal frequencies are identified from acceleration measurements on the cable, then, the time-varying cable tension is obtained from the relation between this force and the identified frequencies. By considering the integer ratios of the different modal frequencies to the fundamental frequency of the cable, the proposed algorithm is further improved to increase its robustness to measurement noise. A cable experiment is implemented to illustrate the validity of the proposed method. For comparison, the Hilbert–Huang transform is also employed to identify the time-varying frequencies, which are then used to calculate the time-varying cable-tension force. The results show that the adaptive sparse time-frequency analysis method produces more accurate estimates of the time-varying cable tension forces than the Hilbert–Huang transform method.

In this follow-up work, we extend the discontinuous Galerkin (DG) methods previously developed on rectangular meshes \cite{previous} to triangular meshes. The DG schemes in \cite{previous} are both optimally convergent and energy conserving. However, as we shall see in the numerical results section, the DG schemes on triangular meshes only have suboptimal convergence rate. We prove the energy conservation and an error estimate for the semi-discrete schemes. The stability of the fully discrete scheme is proved and its error estimate is stated. We present extensive numerical results with convergence consistent of our error estimate, and simulations of wave propagation in Drude metamaterials to demonstrate the flexibility of triangular meshes.