We develop a new compact scheme for second-order PDE (parabolic and Schrödinger type) with a variable time-independent coefficient. It has a higher order and smaller error than classic implicit scheme. The Dirichlet and Neumann boundary problems are considered. The relative finite-difference operator is almost self-adjoint.

The implicit compact finite-difference scheme was developed for evolutionary partial differential parabolic and Schrödinger-type equations and systems with a weak nonlinearity. To make a temporal step of the compact implicit scheme we need to solve a non-linear system. We use for this step a simple explicit difference scheme and then Newton -- Raphson iterations, which are implemented by the double-sweep method. Numerical experiments confirm the 4-th order of an algorithm. The Richardson extrapolation improves it up to the 6-th order.

Dec 15 2017

math.NA arXiv:1712.05174v1

We present a numerical study to investigate the conditioning of the plane wave discontinuous Galerkin discretization of the Helmholtz problem. We provide empirical evidence that the spectral condition number of the plane wave basis on a single element depends algebraically on the mesh size and the wave number, and exponentially on the number of plane wave directions; we also test its dependence on the element shape. We show that the conditioning of the global system can be improved by orthogonalization of the local basis functions with the modified Gram-Schmidt algorithm, which results in significantly fewer GMRES iterations for solving the discrete problem iteratively.

Dec 15 2017

math.NA arXiv:1712.05145v1

In this paper, we first introduce the reader to the Basic Scheme of Moulinec and Suquet in the setting of quasi-static linear elasticity, which takes advantage of the fast Fourier transform on homogenized microstructures to accelerate otherwise time-consuming computations. By means of an asymptotic expansion, a hierarchy of linear problems is derived, whose solutions are looked at in detail. It is highlighted how these generalized homogenization problems depend on each other. We extend the Basic Scheme to fit this new problem class and give some numerical results for the first two problem orders.

Dec 15 2017

math.NA arXiv:1712.05126v1

Exciton diffusion length plays a vital role in the function of opto-electronic devices. Oftentimes, the domain occupied by an organic semiconductor is subject to surface measurement error. In many experiments, photoluminescence over the domain is measured and used as the observation data to estimate this length parameter in an inverse manner based on the least square method. However, the result is sometimes found to be sensitive to the surface geometry of the domain. In this paper, we employ a random function representation for the uncertain surface of the domain. After non-dimensionalization, the forward model becomes a diffusion-type equation over the domain whose geometric boundary is subject to small random perturbations. We propose an asymptotic-based method as an approximate forward solver whose accuracy is justified both theoretically and numerically. It only requires solving several deterministic problems over a fixed domain. Therefore, for the same accuracy requirements we tested here, the running time of our approach is more than one order of magnitude smaller than that of directly solving the original stochastic boundary-value problem by the stochastic collocation method. In addition, from numerical results, we find that the correlation length of randomness is important to determine whether a 1D reduced model is a good surrogate for the 2D model.

Recently, in the paper "The Schwarz alternating method in solid mechanics" by Alejandro Mota, Irina Tezaur and Coleman Alleman [Comput. Methods Appl. Mech. Engrg. 319 (2017) 1951] the authors extended the well known Schwarz alternating method from linear to finite-deformation solid mechanics. They developed and introduced four variants of the Schwarz alternating method, presented proof of geometric convergence of the method and prepared parallel implementation applied to some examples. Unfortunately, the work contains serious errors, both from the point of view of finite-deformation solid mechanics as well as mathematical elasticity.

The characterization of intermittent, multiscale and transient dynamics using data-driven analysis remains an open challenge. We demonstrate an application of the Dynamic Mode Decomposition (DMD) with sparse sampling for the diagnostic analysis of multiscale physics. The DMD method is an ideal spatiotemporal matrix decomposition that correlates spatial features of computational or experimental data to periodic temporal behavior. DMD can be modified into a multiresolution analysis to separate complex dynamics into a hierarchy of multiresolution timescale components, where each level of the hierarchy divides dynamics into distinct background (slow) and foreground (fast) timescales. The multiresolution DMD is capable of characterizing nonlinear dynamical systems in an equation-free manner by recursively decomposing the state of the system into low-rank spatial modes and their temporal Fourier dynamics. Moreover, these multiresolution DMD modes can be used to determined sparse sampling locations which are nearly optimal for dynamic regime classification and full state reconstruction. Specifically, optimized sensors are efficiently chosen using QR column pivots of the DMD library, thus avoiding an NP-hard selection process. We demonstrate the efficacy of the method on several examples, including global sea-surface temperature data, and show that only a small number of sensors are needed for accurate global reconstructions and classification of El Niño events.

Dec 15 2017

math.NA arXiv:1712.04969v1

The aim of this paper is to analyze the long-time behaviour of extended Runge--Kutta--Nyström (ERKN) integrators for highly oscillatory Hamiltonian systems. It is important to note that the well-known Gautschi-type methods of order two yield examples of ERKN integrators. To this end, we use the modulated Fourier expansion to analyze the long-time numerical energy conservation. We will show that the ERKN integrators with symmetry conditions have two almost-invariants and have a near conservation of the total and oscillatory energy over a long term.

Dec 15 2017

math.NA arXiv:1712.04934v1

We present a novel algorithm for high resolution coherent imaging of sound sources in random scattering media using time resolved measurements of the acoustic pressure at an array of receivers. The sound waves travel a long distance between the sources and receivers so that they are significantly affected by scattering in the random medium. We model the scattering effects by large random wavefront distortions, but the results extend to stronger effects, as long as the waves retain some coherence i.e., before the onset of wave diffusion. It is known that scattering in random media can be mitigated in imaging using coherent interferometry (CINT). This method introduces a statistical stabilization in the image formation, at the cost of image blur. We show how to modify the CINT method in order to image wave sources that are too close to each other to be distinguished by CINT alone. We introduce the algorithm from first principles and demonstrate its performance with numerical simulations.