My research interest includes numerical methods (especially spectral method), compressive sensing, elastic problem in solid mechanics, age of information, and other general applied math problems. As a versatile applied math researcher, I am excited to learn and implement various tools in applied math to conduct research. Below are descriptions of my projects.
My research focuses on using physics-informed machine learning (PIML) to solve complex fluid dynamics problems in astronomy. Specifically, we tackle the time-dependent 2D compressible Navier-Stokes equations to model the evolution of rotating accretion flows within protoplanetary disks, influenced by planetary gravitational forces. These forces create sharp, time-dependent features that are challenging to solve with traditional methods.
We developed an enhanced PIML approach that trains neural networks using physical laws, without relying on large sets of labeled data. This method includes innovations such as a self-scalable activation function, adaptive time-marching strategy, biased loss function, and spatial-temporal adaptive sampling. It effectively captured sharp spatial features and provided reliable solutions with less computational effort. This work offers a promising avenue for solving similar complex fluid dynamics problems in various scientific fields.
The animation presents the gas density distribution and the time evolution within a sample protoplanetary disk problem.
Architecture of the neural network used in our work.
On the forefront of scientific computing, Deep Learning (DL), specifically using Deep Neural Networks (DNNs), has emerged as a powerful tool for solving Partial Differential Equations (PDEs). We leverage recent advancements in function approximation using sparsity-based techniques and random sampling to develop an efficient high-dimensional PDE solver based on deep learning (DL). Our approach demonstrates both theoretically and numerically that it competes with a novel, stable, and accurate compressive spectral collocation method. We introduce a new practical existence theorem, establishing the existence of a class of trainable deep neural networks (DNNs) with suitable bounds on the network architecture and a sufficient condition on sample complexity, with logarithmic scaling in dimension. This ensures that the resulting networks stably and accurately approximate a diffusion-reaction PDE with high probability.
Work in progress: Physics-informed deep learning and compressive collocation for high-dimensional diffusion-reaction equations: \\practical existence theory and numerics.
We propose a new method for the numerical solution of high-dimensional PDEs with periodic boundary conditions called compressive Fourier collocation. Combining techniques from sparse high-dimensional approximation and Monte Carlo sampling, our method requires a number of collocation points that scales mildly (i.e., logarithmically) with respect to the dimension of the PDE domain. Focusing on the case of high-dimensional diffusion equations, we provide recovery guarantees based on the recently proposed framework of sparse recovery in subsampled bounded Riesz systems and carry out numerical experiments up to dimension 20. We currently focus on implementing our technique on even higher dimensional domains and comparing it to recently proposed methods based on deep learning.
Publication: Weiqi Wang, Simone Brugiapaglia, Compressive Fourier collocation methods for high-dimensional diffusion equations with periodic boundary conditions, IMA J. Numer. Anal. (2024), pp. drad102. https://doi.org/10.1093/imanum/drad102.
A two-dimensional compressible PDE solution example.
An example of a patrolling route.
Age of Information (AoI) in graph patrolling problem
Our study focuses on optimizing the patrolling routes of unmanned aerial vehicles (UAVs) in grid systems, especially in hazardous environments, to keep information as current as possible. We aim to minimize the time-average Age of Information (AoI) on network edges. For networks that have Eulerian cycles (routes that visit every edge exactly once), we discovered an optimal patrol strategy. However, for more complex networks without such cycles, finding the best strategy is challenging due to the many possible patrol routes.
To simplify this, we developed two approximation methods. First, we construct multigraphs, which are simplified versions of the original network. Then, we embed Eulerian cycles within these multigraphs. These methods ensure that the patrol routes are nearly optimal, achieving an approximation ratio of 2. Our research shows that distributing visits to network edges evenly over time is more effective in reducing the AoI than simply minimizing the total travel distance of the UAVs. Based on this insight, we propose a heuristic approach that maintains efficient performance across different network complexities.
Work in progress: Weiqi Wang, Jin Xu, Balancing Edge Traverses: Minimizing Age of Information in UAV Patrol Routing.
Does elastic stress modify the equilibrium corner angle?
We consider the influence of elasticity and anisotropic surface energy on the corner angle of the energy-minimizing shape of a two-dimensional void under biaxial loading. We develop a numerical method for determining the stress for a class of arbitrary void shapes and corner angles. Our results show that the precise corner angles that minimize the total energy are not affected by the presence of the singular elastic fields. However, the stress singularity on the void surface at the corner must be balanced by a singularity in the curvature at the corner that effectively changes the macroscopic geometry of the shape and effectively changes the apparent corner angle.
Publications: Weiqi Wang, Brian J. Spencer, Does elastic stress modify the equilibrium corner angle?, J Mech Phys Solids 167 (2022), pp. 105003. https://doi.org/10.1016/j.jmps.2022.105003.
Weiqi Wang, Brian J. Spencer, Numerical solution for the stress near a hole with corners in an infinite plate under biaxial loading, J Eng Math 127, 13 (2021). https://doi.org/10.1007/s10665-021-10104-8
Distribution of stress (xx-component) for a overlapping circles hole
I list my past research below.
An example of transition matrix recovered using Gaussian-kernel method.
We propose a non-homogeneous Markov chain model to predict the delay over stations. We build our prediction model by assuming the delay evolution over stations follows a non-homogeneous Markov chain. We propose a chi-square Markov property test to verify our assumption that the delay evolution over stations has a first-order Markov property. Each transition matrix in our model can characterize the unique pattern of delay evolution between two adjacent stations that a train travels through. We propose a Gaussian-kernel-based method to recover the transition matrices for the Markov chain when the training data are limited. This recovery method captures the delay transition probabilities from the existing data. We show that this recovery method achieves a higher prediction accuracy than the other naive approaches.
Publication: Xu, J., Wang, W., Gao, Z., Luo, H., & Wu, Q. (2023). A novel Markov model for near-term railway delay prediction. Computers & Industrial Engineering, 109302. https://doi.org/10.48550/arXiv.2205.10682.
Original power trace (top) and place of leakage (bottom).
I participated in this research project during my internship at Shanghai Fudan Microelectronics Group Company Limited. In computer security, a side-channel attack is any attack based on extra information that can be gathered because of the fundamental way a computer protocol or algorithm is implemented rather than flaws in the design of the protocol or algorithm itself. We developed methods analyzing the power consumption of FPGA provided by the DPA (Differential Power analysis) contest to decode the key bits of the cryptographic algorithms using power analysis. We first use PCA(Principal Component Analysis) to denoise and search for the leakages of the power traces. Then we decode the key by applying clustering algorithms (e.g., maximum likelihood estimation).
The results can be found on DPA contest webpage: https://www.dpacontest.org/v4/42_hall_of_fame.php