Source author record

XiaoFei Zhao

XiaoFei Zhao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2020arXiv

Embedded exponential-type low-regularity integrators for KdV equation under rough data

In this paper, we introduce a novel class of embedded exponential-type low-regularity integrators (ELRIs) for solving the KdV equation and establish their optimal convergence results under rough initial data. The schemes are explicit and efficient to implement. By rigorous error analysis, we first show that the ELRI scheme provides the first order accuracy in $H^γ$ for initial data in $H^{γ+1}$ for $γ>\frac12$. Moreover, by adding two more correction terms to the first order scheme, we show a second order ELRI that provides the second order accuracy in $H^γ$ for initial data in $H^{γ+3}$ for $γ\ge0$. The proposed ELRIs further reduce the regularity requirement of existing methods so far for optimal convergence. The theoretical results are confirmed by numerical experiments, and comparisons with existing methods illustrate the efficiency of the new methods.

preprint2020arXiv

Error estimates of some splitting schemes for charged-particle dynamics under strong magnetic field

In this work, we consider the error estimates of some splitting schemes for the charged-particle dynamics under a strong magnetic field. We first propose a novel energy-preserving splitting scheme with computational cost per step independent from the strength of the magnetic field. Then under the maximal ordering scaling case, we establish for the scheme and in fact for a class of Lie-Trotter type splitting schemes, a uniform (in the strength of the magnetic field) and optimal error bound in the position and in the velocity parallel to the magnetic field. For the general strong magnetic field case, the modulated Fourier expansions of the exact and the numerical solutions are constructed to obtain a convergence result. Numerical experiments are presented to illustrate the error and energy behaviour of the splitting schemes.

preprint2020arXiv

Numerical integrators for continuous disordered nonlinear Schrödinger equation

In this paper, we consider the numerical solution of the continuous disordered nonlinear Schrödinger equation, which contains a spatial random potential. We address the finite time accuracy order reduction issue of the usual numerical integrators on this problem, which is due to the presence of the random/rough potential. By using the recently proposed low-regularity integrator (LRI) from (33, SIAM J. Numer. Anal., 2019), we show how to integrate the potential term by losing two spatial derivatives. Convergence analysis is done to show that LRI has the second order accuracy in $L^2$-norm for potentials in $H^2$. Numerical experiments are done to verify this theoretical result. More numerical results are presented to investigate the accuracy of LRI compared with classical methods under rougher random potentials from applications.

preprint2020arXiv

Optimal convergence of a second order low-regularity integrator for the KdV equation

In this paper, we establish the optimal convergence result of a second order exponential-type integrator from (136, Numer. Math., 2017) for solving the KdV equation under rough initial data. The scheme is explicit and efficient to implement. By rigorous error analysis, we show that the scheme provides the second order accuracy in $H^γ$ for initial data in $H^{γ+4}$ for any $γ\geq0$, where the regularity requirement is lower than the classical methods. The result is confirmed by numerical experiments and comparisons are made with the Strang splitting scheme.

preprint2019arXiv

Estimating the sample mean and standard deviation from commonly reported quantiles in meta-analysis

Researchers increasingly use meta-analysis to synthesize the results of several studies in order to estimate a common effect. When the outcome variable is continuous, standard meta-analytic approaches assume that the primary studies report the sample mean and standard deviation of the outcome. However, when the outcome is skewed, authors sometimes summarize the data by reporting the sample median and one or both of (i) the minimum and maximum values and (ii) the first and third quartiles, but do not report the mean or standard deviation. To include these studies in meta-analysis, several methods have been developed to estimate the sample mean and standard deviation from the reported summary data. A major limitation of these widely used methods is that they assume that the outcome distribution is normal, which is unlikely to be tenable for studies reporting medians. We propose two novel approaches to estimate the sample mean and standard deviation when data are suspected to be non-normal. Our simulation results and empirical assessments show that the proposed methods often perform better than the existing methods when applied to non-normal data.

preprint2018arXiv

Uniformly accurate methods for Vlasov equations with non-homogeneous strong magnetic field

In this paper, we consider the numerical solution of highly-oscillatory Vlasov and Vlasov-Poisson equations with non-homogeneous magnetic field. Designed in the spirit of recent uniformly accurate methods, our schemes remain insensitive to the stiffness of the problem, in terms of both accuracy and computational cost. The specific difficulty (and the resulting novelty of our approach) stems from the presence of a non-periodic oscillation, which necessitates a careful ad-hoc reformulation of the equations. Our results are illustrated numerically on several examples.

preprint2017arXiv

One-sample aggregate data meta-analysis of medians

An aggregate data meta-analysis is a statistical method that pools the summary statistics of several selected studies to estimate the outcome of interest. When considering a continuous outcome, typically each study must report the same measure of the outcome variable and its spread (e.g., the sample mean and its standard error). However, some studies may instead report the median along with various measures of spread. Recently, the task of incorporating medians in meta-analysis has been achieved by estimating the sample mean and its standard error from each study that reports a median in order to meta-analyze the means. In this paper, we propose two alternative approaches to meta-analyze data that instead rely on medians. We systematically compare these approaches via simulation study to each other and to methods that transform the study-specific medians and spread into sample means and their standard errors. We demonstrate that the proposed median-based approaches perform better than the transformation-based approaches, especially when applied to skewed data and data with high inter-study variance. In addition, when meta-analyzing data that consists of medians, we show that the median-based approaches perform considerably better than or comparably to the best-case scenario for a transformation approach: conducting a meta-analysis using the actual sample mean and standard error of the mean of each study. Finally, we illustrate these approaches in a meta-analysis of patient delay in tuberculosis diagnosis.

preprint2016arXiv

A modulation equations approach for numerically solving the moving soliton and radiation solutions of NLS

Based on our previous work for solving the nonlinear Schrodinger equation with multichannel dynamics that is given by a localized standing wave and radiation, in this work we deal with the multichannel solution which consists of a moving soliton and radiation. We apply the modulation theory to give a system of ODEs coupled to the radiation term for describing the solution, which is valid for all times. The modulation equations are solved accurately by the proposed numerical method. The soliton and radiation are captured separately in the computation, and they are solved on the translated domain that is moving with them. Thus for a fixed finite physical domain in the lab frame, the multichannel solution can pass through the boundary naturally, which can not be done by imposing any existing boundary conditions. We comment on the differences of this method from the collective coordinates.

preprint2016arXiv

Modulation equations approach for solving vortex and radiation in nonlinear Schrodinger equation

We apply the modulation theory to study the vortex and radiation solution in the two-dimensional nonlinear Schrödinger equation. The full modulation equations which describe the dynamics of the vortex and radiation separately are derived. A general algorithm is proposed to efficiently and accurately find vortices with different values of energy and spin index. The modulation equations are solved by accurate numerical method. Numerical tests and simulations of scattering are given.

preprint2016arXiv

Symmetric high order Gautschi-type exponential wave integrators pseudospectral method for the nonlinear Klein-Gordon equation in the nonrelativistic limit regime

A group of high order Gautschi-type exponential wave integrators (EWIs) Fourier pseudospectral method are proposed and analyzed for solving the nonlinear Klein-Gordon equation (KGE) in the nonrelativistic limit regime, where a parameter $0<\varepsilon\ll1$ which is inversely proportional to the speed of light, makes the solution propagate waves with wavelength $O(\varepsilon^2)$ in time and $O(1)$ in space. With the Fourier pseudospectral method to discretize the KGE in space, we propose a group of EWIs with designed Gautschi's type quadratures for the temporal integrations, which can offer any intended even order of accuracy provided that the solution is smooth enough, while all the current existing EWIs offer at most second order accuracy. The scheme is explicit, time symmetric and rigorous error estimates show the meshing strategy of the proposed method is time step $τ=O(\varepsilon^2)$ and mesh size $h=O(1)$ as $0<\varepsilon\ll1$, which is `optimal' among all classical numerical methods towards solving the KGE directly in the limit regime, and which also distinguish our methods from other high order approaches such as Runge-Kutta methods which require $τ=O(\varepsilon^3)$. Numerical experiments with comparisons are done to confirm the error bound and show the superiority of the proposed methods over existing classical numerical methods.

preprint2016arXiv

Uniformly accurate numerical schemes for the nonlinear Dirac equation in the nonrelativistic limit regime

We apply the two-scale formulation approach to propose uniformly accurate (UA) schemes for solving the nonlinear Dirac equation in the nonrelativistic limit regime. The nonlinear Dirac equation involves two small scales $\varepsilon$ and $\varepsilon^2$ with $\varepsilon\to0$ in the nonrelativistic limit regime. The small parameter causes high oscillations in time which brings severe numerical burden for classical numerical methods. We transform our original problem as a two-scale formulation and present a general strategy to tackle a class of highly oscillatory problems involving the two small scales $\varepsilon$ and $\varepsilon^2$. Suitable initial data for the two-scale formulation is derived to bound the time derivatives of the augmented solution. Numerical schemes with uniform (with respect to $\varepsilon\in (0,1]$) spectral accuracy in space and uniform first order or second order accuracy in time are proposed. Numerical experiments are done to confirm the UA property.

preprint2015arXiv

On multichannel solutions of nonlinear Schrödinger equations: algorithm, analysis and numerical explorations

We apply the method of modulation equations to numerically solve the NLS with multichannel dynamics, given by a trapped localized state and radiation. This approach employs the modulation theory of Soffer-Weinstein, which gives a system of ODE's coupled to the radiation term, which is valid for all times. We comment on the differences of this method from the well-known method of collective coordinates.

preprint2014arXiv

A uniformly accurate multiscale time integrator pseudospectral method for the Klein-Gordon equation in the nonrelativistic limit regime

We propose and analyze a multiscale time integrator Fourier pseudospectral (MTI-FP) method for solving the Klein-Gordon (KG) equation with a dimensionless parameter $0<\varepsilon\leq1$ which is inversely proportional to the speed of light. In the nonrelativistic limit regime, i.e. $0<\varepsilon\ll1$, the solution to the KG equation propagates waves with amplitude at $O(1)$ and wavelength at $O(\varepsilon^2)$ in time and $O(1)$ in space, which causes significantly numerical burdens due to the high oscillation in time. The MTI-FP method is designed by adapting a multiscale decomposition by frequency (MDF) to the solution at each time step and applying an exponential wave integrator to the nonlinear Schrödinger equation with wave operator under well-prepared initial data for $\varepsilon^2$-frequency and $O(1)$-amplitude waves and a KG-type equation with small initial data for the reminder waves in the MDF. We rigorously establish two independent error bounds in $H^2$-norm to the MTI-FP method at $O(h^{m_0}+τ^2+\varepsilon^2)$ and $O(h^{m_0}+τ^2/\varepsilon^2)$ with $h$ mesh size, $τ$ time step and $m_0\ge2$ an integer depending on the regularity of the solution, which immediately imply that the MTI-FP converges uniformly and optimally in space with exponential convergence rate if the solution is smooth, and uniformly in time with linear convergence rate at $O(τ)$ for all $\varepsilon\in(0,1]$ and optimally with quadratic convergence rate at $O(τ^2)$ in the regimes when either $\varepsilon=O(1)$ or $0<\varepsilon\le τ$. Numerical results are reported to confirm the error bounds and demonstrate the efficiency and accuracy of the MTI-FP method for the KG equation, especially in the nonrelativistic limit regime.

preprint2014arXiv

Scalar-field theory of dark matter

We develop a theory of dark matter based on a previously proposed picture, in which a complex vacuum scalar field makes the universe a superfluid, with the energy density of the superfluid giving rise to dark energy, and variations from vacuum density giving rise to dark matter. We formulate a nonlinear Klein-Gordon equation to describe the superfluid, treating galaxies as external sources. We study the response of the superfluid to the galaxies, in particular, the emergence of the dark-matter galactic halo, contortions during galaxy collisions, and the creation of vortices due to galactic rotation.