No abstract uploaded!

We present a uniformly accurate finite difference method and establish rigorously its uniform error bounds for the Zakharov system (ZS) with a dimensionless parameter 0 < ε ≤ 1, which is inversely proportional to the speed of sound. In the subsonic limit regime, i.e., 0 < ε << 1, the solution propagates highly oscillatory waves and/or rapid outgoing initial layers due to the perturbation of the wave operator in the ZS and/or the incompatibility of the initial data, which is characterized by two parameters α ≥ 0 and β ≥ −1. Specifically, the solution propagates waves with wavelength of O(ε) and O(1) in time and space, respectively, and amplitude at O(ε^{min{2,α,1+β}}). This high oscillation of the solution in time brings significant difficulties in designing numerical methods and establishing their error bounds, especially in the subsonic limit regime. A uniformly accurate finite difference method is proposed by reformulating the ZS into an asymptotic consistent formulation and adopting an integral approximation of the oscillatory term. By applying the energy method and using the limiting equation via a nonlinear Schr ̈odinger equation with an oscillatory potential, we rigorously establish two independent error bounds at O(h^2 + τ^2/ε) and O(h^2 + τ^2 + τε^{α^∗} + ε^{1+α^∗}), respectively, with h the mesh size, τ the time step and α^∗ = min{1, α, 1 + β}. Thus we obtain error bounds at O(h^2 + τ^{4/3}) and O(h^2 + τ^{1+α^∗/(2+α^∗)}) for well-prepared (α^∗ = 1) and ill-prepared (0 ≤ α^∗ < 1) initial data, respectively, which are uniform in both space and time for 0 < ε ≤ 1 and optimal at the second order in space. Other techniques in the analysis include the cut-off technique for treating the nonlinearity and inverse estimates to bound the numerical solution. Numerical results are reported to demonstrate our error bounds.

Two fundamental difficulties are encountered in the numerical evaluation of time- dependent layer potentials. One is the quadratic cost of history dependence, which has been successfully addressed by splitting the potentials into two parts—a local part that contains the most recent contributions and a history part that contains the con- tributions from all earlier times. The history part is smooth, easily discretized using high-order quadratures, and straightforward to compute using a variety of fast algo- rithms. The local part, however, involves complicated singularities in the underlying Green’s function. Existing methods, based on exchanging the order of integration in space and time, are able to achieve high-order accuracy, but are limited to the case of stationary boundaries. Here, we present a new quadrature method that leaves the order of integration unchanged, making use of a change of variables that converts the singular integrals with respect to time into smooth ones. We have also derived asymptotic formulas for the local part that lead to fast and accurate hybrid schemes, extending earlier work for scalar heat potentials and applicable to moving boundaries. The performance of the overall scheme is demonstrated via numerical examples.

In the paper [4], the authors proposed a derivative-free descent algorithm for the nonlinear complementarity problems (NCPs) by the generalized Fischer-Burmeister merit function: p (a, b)= 1

Recent technical advances lead to the coupling of PET and MRI scanners, enabling one to acquire functional and anatomical data simultaneously. In this paper, we propose a tight frame based PET-MRI joint reconstruction model via the joint sparsity of tight frame coefficients. In addition, a nonconvex balanced approach is adopted to take the different regularities of PET and MRI images into account. To solve the nonconvex and nonsmooth model, a proximal alternating minimization algorithm is proposed, and the global convergence is present based on the Kurdyka--Łojasiewicz property. Finally, the numerical experiments show that our proposed models achieve better performance over the existing PET-MRI joint reconstruction models.