BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jingwei Hu

Jingwei Hu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.APastro-ph.HEmath-phmath.MPphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2026arXiv

A Separable and Asymptotic-Preserving Dynamical Low-Rank Method for the Vlasov--Poisson--Fokker--Planck System

We present a dynamical low-rank (DLR) method for the Vlasov--Poisson--Fokker--Planck (VPFP) system. Our main contributions are two-fold: (i) a conservative spatial discretization of the Fokker--Planck operator that factors into velocity-only and space-only components, enabling efficient low-rank projection, and (ii) a time discretization within the DLR framework that properly handles stiff collisions. We propose both first-order and second-order low-rank IMEX schemes. For the first-order scheme, we prove an asymptotic-preserving (AP) property when the field fluctuation is small. Numerical experiments demonstrate accuracy, robustness, and AP property at modest ranks.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

A Structure-Preserving Penalization Method for the Single-species Rosenbluth-Fokker-Planck Equation

The Rosenbluth-Fokker-Planck (RFP) equation describes Coulomb collisional dynamics within and across species in plasmas. It belongs to the broader class of anisotropic-diffusion-advection equations, whose numerical approximation is highly-nontrivial due to its nonlinearity, stiffness, and structural properties such as conservation and entropy dissipation (hence with the Maxwellian distribution as the equilibrium state). In this paper, we propose a structure-preserving penalization scheme for the stiff, single-species RFP equation. The scheme features three novel components: 1) a novel generalization of the well-known Chang-Cooper discretization for the RFP equation that is equilibrium-preserving and enables positivity while preserving mass, momentum, and energy; 2) an easy-to-invert isotropic variable-coefficient penalization operator to deal with the temporal stiffness without resorting to a fully implicit scheme, borrowing ideas from explicit-implicit-null (EIN) methods, and 3) an adaptive timestepping strategy that preserves the positivity of the full penalized scheme. The resulting scheme conserves mass, momentum, and energy strictly, is unconditionally stable, and robustly positivity preserving. The scheme is demonstrated with linear and nonlinear anisotropic diffusion examples of increasing complexity, including several single-species RFP examples.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A structure preserving numerical scheme for Fokker-Planck equations of structured neural networks with learning rules

In this work, we are concerned with a Fokker-Planck equation related to the nonlinear noisy leaky integrate-and-fire model for biological neural networks which are structured by the synaptic weights and equipped with the Hebbian learning rule. The equation contains a small parameter $\varepsilon$ separating the time scales of learning and reacting behavior of the neural system, and an asymptotic limit model can be derived by letting $\varepsilon\to 0$, where the microscopic quasi-static states and the macroscopic evolution equation are coupled through the total firing rate. To handle the endowed flux-shift structure and the multi-scale dynamics in a unified framework, we propose a numerical scheme for this equation that is mass conservative, unconditionally positivity preserving, and asymptotic preserving. We provide extensive numerical tests to verify the schemes' properties and carry out a set of numerical experiments to investigate the model's learning ability, and explore the solution's behavior when the neural network is excitatory.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Systematic Study of the Short-Term X-ray Variability of Seyfert Galaxies I. Diversity of the X-ray Rms Spectra

The X-ray variability of active galactic nuclei (AGN) carries crucial information about the X-ray radiation mechanism. We performed a systematic study of the X-ray short-term (1-100 ks timescale) variability for a large sample of 78 Seyferts with 426 deep XMM-Newton observations. In this paper, we present the time-averaged spectra and rms spectra for the entire sample, which show a variety of properties. Based on the spectral shape, we divide the rms spectra into five subtypes and the time-averaged spectra into four subtypes. The most common shape of the rms spectra is concave-down where the rms peaks at $\sim$ 1 keV. We find that different sources can show similar time-averaged spectra and rms spectra. However, there is no one-to-one mapping between the subtypes of the time-averaged spectra and rms spectra, as similar time-averaged spectra can be accompanied by different rms spectra, and vice versa. This is likely because different physical mechanisms can produce similar rms spectra. For every subtype of the time-averaged spectra, we report its preferred subtypes of the rms spectra in both low- and high-frequency bands. We also compare the statistical properties for different subtypes, such as the black hole mass and Eddington ratio. Finally, we investigate the rms in the Fe K$α$ line regime and find that those with a broad and extended red-wing profile tend to show stronger variability than those showing a narrow or relatively symmetric profile. Our results demonstrate the necessity of performing joint spectral and variability modeling in order to understand the mechanism of the X-ray emission in AGN. All of the rms spectra have been made publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Efficient dynamical low-rank approximation for the Vlasov-Ampère-Fokker-Planck system

Kinetic equations are difficult to solve numerically due to their high dimensionality. A promising approach for reducing computational cost is the dynamical low-rank algorithm, which decouples the dimensions of the phase space by proposing an ansatz as the sum of separable (rank-1) functions in position and velocity respectively. The fluid asymptotic limit of collisional kinetic equations, obtained in the small-Knudsen number limit, admits a low-rank representation when written as $f = Mg$, where $M$ is the local Maxwellian, and $g$ is low-rank. We apply this decomposition to the Vlasov-Ampère-Fokker-Planck equation of plasma dynamics, considering the asymptotic limit of strong collisions and electric field. We implement our proposed algorithm and demonstrate the expected improvement in computation time by comparison to an implementation that evolves the full solution tensor $f$. We also demonstrate that our algorithm can capture dynamics in both the kinetic regime, and in the fluid regime with relatively lower computational effort, thereby efficiently capturing the asymptotic fluid limit.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Implicit Regularization Effects of the Sobolev Norms in Image Processing

In this paper, we propose to use the general $L^2$-based Sobolev norms, i.e., $H^s$ norms where $s\in \mathbb{R}$, to measure the data discrepancy due to noise in image processing tasks that are formulated as optimization problems. As opposed to a popular trend of developing regularization methods, we emphasize that an implicit regularization effect can be achieved through the class of Sobolev norms as the data-fitting term. Specifically, we analyze that the implicit regularization comes from the weights that the $H^s$ norm imposes on different frequency contents of an underlying image. We further analyze the underlying noise assumption of using the Sobolev norm as the data-fitting term from a Bayesian perspective, build the connections with the Sobolev gradient-based methods and discuss the preconditioning effects on the convergence rate of the gradient descent algorithm, leading to a better understanding of functional spaces/metrics and the optimization process involved in image processing. Numerical results in full waveform inversion, image denoising and deblurring demonstrate the implicit regularization effects.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Luminosity function and event rate density of XMM-Newton-selected supernova shock-breakout candidates

A dozen X-ray supernova shock breakout (SN SBO) candidates were reported recently based on XMM-Newton archival data, which increased the X-ray selected SN SBO sample by an order of magnitude. Assuming they are genuine SN SBOs, we study the luminosity function (LF) by improving upon the method used in our previous work. The light curves and the spectra of the candidates were used to derive the maximum volume within which these objects could be detected with XMM-Newton by simulation. The results show that the SN SBO LF can be described by either a broken power law (BPL) with indices (at the 68$\%$ confidence level) of $0.48 \pm 0.28$ and $2.11 \pm 1.27$ before and after the break luminosity at $\log (L_b/\rm erg\,s^{-1})=$ $45.32 \pm 0.55$ or a single power law (SPL) with index of $0.80 \pm 0.16$. The local event rate densities of SN SBOs above $5\times 10^{42}$ $\rm erg\,s^{-1}$ are consistent for two models, i.e., $4.6 ^{+1.7}_{-1.3} \times 10^4$ and $4.9 ^{+1.9}_{-1.4} \times 10^4$ $\rm Gpc^{-3}\,yr^{-1}$ for BPL and SPL models, respectively. The number of fast X-ray transients of SN SBO origin can be significantly increased by the wide-field X-ray telescopes such as the Einstein Probe.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Positivity-preserving and energy-dissipative finite difference schemes for the Fokker-Planck and Keller-Segel equations

In this work, we introduce semi-implicit or implicit finite difference schemes for the continuity equation with a gradient flow structure. Examples of such equations include the linear Fokker-Planck equation and the Keller-Segel equations. The two proposed schemes are first order accurate in time, explicitly solvable, and second order and fourth order accurate in space, which are obtained via finite difference implementation of the classical continuous finite element method. The fully discrete schemes are proven positivity-preserving and energy-dissipative: the second order scheme can achieve so unconditionally; the fourth order scheme only requires a mild time step and mesh size constraint. Furthermore, the fourth order scheme is the first high order spatial discretization that can achieve both positivity and energy decay properties, which is suitable for long time simulation and to obtain accurate steady state solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A new stability and convergence proof of the Fourier-Galerkin spectral method for the spatially homogeneous Boltzmann equation

Numerical approximation of the Boltzmann equation is a challenging problem due to its high-dimensional, nonlocal, and nonlinear collision integral. Over the past decade, the Fourier-Galerkin spectral method has become a popular deterministic method for solving the Boltzmann equation, manifested by its high accuracy and potential of being further accelerated by the fast Fourier transform. Albeit its practical success, the stability of the method is only recently proved by Filbet, F. & Mouhot, C. in [$ Trans. Amer. Math. Soc.$ 363, no. 4 (2011): 1947-1980.] by utilizing the "spreading" property of the collision operator. In this work, we provide a new proof based on a careful $L^2$ estimate of the negative part of the solution. We also discuss the applicability of the result to various initial data, including both continuous and discontinuous functions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A particle method for the homogeneous Landau equation

We propose a novel deterministic particle method to numerically approximate the Landau equation for plasmas. Based on a new variational formulation in terms of gradient flows of the Landau equation, we regularize the collision operator to make sense of the particle solutions. These particle solutions solve a large coupled ODE system that retains all the important properties of the Landau operator, namely the conservation of mass, momentum and energy, and the decay of entropy. We illustrate our new method by showing its performance in several test cases including the physically relevant case of the Coulomb interaction. The comparison to the exact solution and the spectral method is strikingly good maintaining 2nd order accuracy. Moreover, an efficient implementation of the method via the treecode is explored. This gives a proof of concept for the practical use of our method when coupled with the classical PIC method for the Vlasov equation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the uniform accuracy of implicit-explicit backward differentiation formulas (IMEX-BDF) for stiff hyperbolic relaxation systems and kinetic equations

Many hyperbolic and kinetic equations contain a non-stiff convection/transport part and a stiff relaxation/collision part (characterized by the relaxation or mean free time $\varepsilon$). To solve this type of problems, implicit-explicit (IMEX) multistep methods have been widely used and their performance is understood well in the non-stiff regime ($\varepsilon=O(1)$) and limiting regime ($\varepsilon\rightarrow 0$). However, in the intermediate regime (say, $\varepsilon=O(Δt)$), uniform accuracy has been reported numerically without a complete theoretical justification (except some asymptotic or stability analysis). In this work, we prove the uniform accuracy -- an optimal {\it a priori} error bound -- of a class of IMEX multistep methods, IMEX backward differentiation formulas (IMEX-BDF), for linear hyperbolic systems with stiff relaxation. The proof is based on the energy estimate with a new multiplier technique. For nonlinear hyperbolic and kinetic equations, we numerically verify the same property using a series of examples.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Recent development in kinetic theory of granular materials: analysis and numerical methods

Over the past decades, kinetic description of granular materials has received a lot of attention in mathematical community and applied fields such as physics and engineering. This article aims to review recent mathematical results in kinetic granular materials, especially for those which arose since the last review by Villani on the same subject. We will discuss both theoretical and numerical developments. We will finally showcase some important open problems and conjectures by means of numerical experiments based on spectral methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Uncertainty quantification for the BGK model of the Boltzmann equation using multilevel variance reduced Monte Carlo methods

We propose a control variate multilevel Monte Carlo method for the kinetic BGK model of the Boltzmann equation subject to random inputs. The method combines a multilevel Monte Carlo technique with the computation of the optimal control variate multipliers derived from local or global variance minimization problems. Consistency and convergence analysis for the method equipped with a second-order positivity-preserving and asymptotic-preserving scheme in space and time is also performed. Various numerical examples confirm that the optimized multilevel Monte Carlo method outperforms the classical multilevel Monte Carlo method especially for problems with discontinuities.

0 citations0 reviewsNo code linkedOpen access