We establish uniform error bounds of time-splitting Fourier pseudospectral (TSFP) methods for the nonlinear Klein–Gordon equation (NKGE) with weak power-type nonlinearity and O(1) initial data, while the nonlinearity strength is characterized by εp with a constant p ∈ N+ and a dimensionless parameter ε ∈ (0, 1], for the long-time dynamics up to the time at O(ε^{−β}) with 0 ≤ β ≤ p. In fact, when 0 < ε ≪ 1, the problem is equivalent to the long-time dynamics of NKGE with small initial data and O(1) nonlinearity strength, while the amplitude of the initial data (and the solution) is at O(ε). By reformulating the NKGE into a relativistic nonlinear Schr ̈odinger equation, we adapt the TSFP method to discretize it numerically. By using the method of mathematical induction to bound the numerical solution, we prove uniform error bounds at O(h^m + ε^{p−β} τ^2) of the TSFP method with h mesh size, τ time step and m ≥ 2 depending on the regularity of the solution. The error bounds are uniformly accurate for the long-time simulation up to the time at O(ε^{−β}) and uniformly valid for ε ∈ (0, 1]. Especially, the error bounds are uniformly at the second order rate for the large time step τ = O(ε^{−(p−β)/2}) in the parameter regime 0 ≤ β < p. Numerical results are reported to confirm our error bounds in the long-time regime. Finally, the TSFP method and its error bounds are extended to a highly oscillatory complex NKGE which propagates waves with wavelength at O(1) in space and O(ε^β) in time and wave velocity at O(ε^{−β}).

In this paper, we propose an unfitted Nitsche's method to compute the band structures of phononic crystal with impurities of general geometry. The proposed method does not require the background mesh to fit the interfaces of impurities, and thus avoids the expensive cost of generating body-fitted meshes and simplifies the inclusion of interface conditions in the formulation. The quasi-periodic boundary conditions are handled by the Floquet-Bloch transform, which converts the computation of band structures into an eigenvalue problem with periodic boundary conditions. More importantly, we show the well-posedness of the proposed method using a delicate argument based on the trace inequality, and further prove the convergence by the Babuška-Osborn theory. We achieve the optimal convergence rate at the presence of the impurities of general geometry. We confirm the theoretical results by two numerical examples, and show the capability of the proposed methods for computing the band structures without fitting the interfaces of impurities.

In this paper, we present two adaptive methods for the basis enrichment of the mixed Generalized Multiscale Finite Element Method (GMsFEM) for solving the flow problem in heterogeneous media. We develop an a-posteriori error indicator which depends on the norm of a local residual operator. Based on this indicator, we construct an offline adaptive method to increase the number of basis functions locally in coarse regions with large local residuals. We also develop an online adaptive method which iteratively enriches the function space by adding new functions computed based on the residual of the previous solution and special minimum energy snapshots. We show theoretically and numerically the convergence of the two methods. The online method is, in general, better than the offline method as the online method is able to capture distant effects (at a cost of online computations), and both

In this paper, we establish a connection between the recently developed data-driven timefrequency analysis [10, 11] and the classical second order differential equations. The main idea of the data-driven time-frequency analysis is to decompose a multiscale signal into the sparsest collection of Intrinsic Mode Functions (IMFs) over the largest possible dictionary via nonlinear optimization. These IMFs are of the form a(t) cos(θ(t)) where the amplitude a(t) is positive and slowly varying. The non-decreasing phase function θ(t) is determined by the data and in general depends on the signal in a nonlinear fashion. One of the main results of this paper is that we show that each IMF can be associated with a solution of a second order ordinary differential equation of the form x¨ + p(x, t)x˙ + q(x, t) = 0. Further, we propose a localized variational formulation for this problem and develop an effective l1-based optimization method to recover p(x, t) and q(x, t) by looking for a sparse representation of p and q in terms of the polynomial basis. Depending on the form of nonlinearity in p(x, t) and q(x, t), we can define the order of nonlinearity for the associated IMF. This generalizes a concept recently introduced by Prof. N. E. Huang et al. [15]. Numerical examples will be provided to illustrate the robustness and stability of the proposed method for data with or without noise. This manuscript should be considered as a proof of concept.

The equilibrium metric for minimizing a continuous congested traffic model is the solution of a variational problem involving geodesic distances. The continuous equilibrium metric and its associated variational problem are closely related to the classical discrete Wardrops equilibrium. We propose an adjoint state method to numerically approximate continuous traffic congestion equilibria through the continuous formulation. The method formally derives an adjoint state equation to compute the gradient descent direction so as to minimize a nonlinear functional involving the equilibrium metric and the resulting geodesic distances. The geodesic distance needed for the state equation is computed by solving a factored eikonal equation, and the adjoint state equation is solved by a fast sweeping method. Numerical examples demonstrate that the proposed adjoint state method produces desired equilibrium metrics and