Researcher profile

Xiaozhe Hu

Xiaozhe Hu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
8topics
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

14 published item(s)

preprint2025arXiv

A unified spatiotemporal formulation with physics-preserving structure for time-dependent convection-diffusion problems

We propose a unified four-dimensional (4D) spatiotemporal formulation for time-dependent convection-diffusion problems that preserves underlying physical structures. By treating time as an additional space-like coordinate, the evolution problem is reformulated as a stationary convection-diffusion equation on a 4D space-time domain. Using exterior calculus, we extend this framework to the full family of convection-diffusion problems posed on $H(\textbf{grad})$, $H(\textbf{curl})$, and $H(\text{div})$. The resulting formulation is based on a 4D Hodge-Laplacian operator with a spatiotemporal diffusion tensor and convection field, augmented by a small temporal perturbation to ensure nondegeneracy. This formulation naturally incorporates fundamental physical constraints, including divergence-free and curl-free conditions. We further introduce an exponentially-fitted 4D spatiotemporal flux operator that symmetrizes the convection-diffusion operator and enables a well-posed variational formulation. Finally, we prove that the temporally-perturbed formulation converges to the original time-dependent convection-diffusion model as the perturbation parameter tends to zero.

preprint2022arXiv

A Deep Neural Network/Meshfree Method for Solving Dynamic Two-phase Interface Problems

In this paper, a meshfree method using the deep neural network (DNN) approach is developed for solving two kinds of dynamic two-phase interface problems governed by different dynamic partial differential equations on either side of the stationary interface with the jump and high-contrast coefficients. The first type of two-phase interface problem to be studied is the fluid-fluid (two-phase flow) interface problem modeled by Navier-Stokes equations with high-contrast physical parameters across the interface. The second one belongs to fluid-structure interaction (FSI) problems modeled by Navier-Stokes equations on one side of the interface and the structural equation on the other side of the interface, both the fluid and the structure interact with each other via the kinematic- and the dynamic interface conditions across the interface. The DNN/meshfree method is respectively developed for the above two-phase interface problems by representing solutions of PDEs using the DNNs' structure and reformulating the dynamic interface problem as a least-squares minimization problem based upon a space-time sampling point set. Approximation error analyses are also carried out for each kind of interface problem, which reveals an intrinsic strategy about how to efficiently build a sampling-point training dataset to obtain a more accurate DNNs' approximation. In addition, compared with traditional discretization approaches, the proposed DNN/meshfree method and its error analysis technique can be smoothly extended to many other dynamic interface problems with fixed interfaces. Numerical experiments are conducted to illustrate the accuracies of the proposed DNN/meshfree method for the presented two-phase interface problems. Theoretical results are validated to some extent through three numerical examples.

preprint2022arXiv

A Multigrid Preconditioner for Spatially Adaptive High-order Meshless Method on Fluid-solid Interaction Problems

We present a monolithic geometric multigrid preconditioner for solving fluid-solid interaction problems in Stokes limit. The problems are discretized by a spatially adaptive high-order meshless method, the generalized moving least squares (GMLS) with adaptive $h$-refinement. In Stokes limit, solid kinematics can be dominated by the singularities governing the lubrication effects. Resolving those singularities with adaptive $h$-refinement can lead to an ill-conditioned linear system of equations. For constructing the interpolation and restriction operators - the key ingredients of the multigrid preconditioner, we utilize the geometric information of hierarchical sets of GMLS nodes generated in adaptive $h$-refinement. We build decoupled smoothers through physics-based splitting and then combine them via a multiplicative overlapping Schwarz approach. Through numerical examples with the inclusion of different numbers and shapes of solid bodies, we demonstrate the performance and assess the scalability of the designed preconditioner. As the total degrees of freedom and the number of solid bodies $N_s$ increase, the proposed monolithic geometric multigrid preconditioner can ensure convergence and good scalability when using the Krylov iterative method for solving the linear systems of equations generated from the spatially adaptive GMLS discretization. More specifically, for a fixed number of solid bodies, as the discretization resolution is incrementally refined, the number of iterations of the linear solver can be maintained at the same level, indicating nearly linear scalability of our preconditioner with respect to the total degrees of freedom. When $N_s$ increases, the number of iterations is nearly proportional to $\sqrt{N_s}$, implying the sublinear optimality with respect to the number of solid bodies.

preprint2022arXiv

An Enriched Galerkin Method for the Stokes Equations

We present a new enriched Galerkin (EG) scheme for the Stokes equations based on piecewise linear elements for the velocity unknowns and piecewise constant elements for the pressure. The proposed EG method augments the conforming piecewise linear space for velocity by adding an additional degree of freedom which corresponds to one discontinuous linear basis function per element. Thus, the total number of degrees of freedom is significantly reduced in comparison with standard conforming, non-conforming, and discontinuous Galerkin schemes for the Stokes equation. We show the well-posedness of the new EG approach and prove that the scheme converges optimally. For the solution of the resulting large-scale indefinite linear systems we propose robust block preconditioners, yielding scalable results independent of the discretization and physical parameters. Numerical results confirm the convergence rates of the discretization and also the robustness of the linear solvers for a variety of test problems.

preprint2022arXiv

Monolithic multigrid for a reduced-quadrature discretization of poroelasticity

Advanced finite-element discretizations and preconditioners for models of poroelasticity have attracted significant attention in recent years. The equations of poroelasticity offer significant challenges in both areas, due to the potentially strong coupling between unknowns in the system, saddle-point structure, and the need to account for wide ranges of parameter values, including limiting behavior such as incompressible elasticity. This paper was motivated by an attempt to develop monolithic multigrid preconditioners for the discretization developed in [48]; we show here why this is a difficult task and, as a result, we modify the discretization in [48] through the use of a reduced quadrature approximation, yielding a more "solver-friendly" discretization. Local Fourier analysis is used to optimize parameters in the resulting monolithic multigrid method, allowing a fair comparison between the performance and costs of methods based on Vanka and Braess-Sarazin relaxation. Numerical results are presented to validate the LFA predictions and demonstrate efficiency of the algorithms. Finally, a comparison to existing block-factorization preconditioners is also given.

preprint2022arXiv

Solving Graph Laplacians via Multilevel Sparsifiers

We consider effective preconditioners for solving Laplacians of general weighted graphs. Theoretically, spectral sparsifiers (SSs) provide preconditioners of optimal computational complexity. However, they are not easy to use for real-world applications due to the implementation complications. Multigrid (MG) methods, on the contrary, are computationally efficient but lack of theoretical justifications. To bridge the gap between theory and practice, we adopt ideas of MG and SS methods and proposed preconditioners that can be used in practice with theoretical guarantees. We expand the original graph based on a multilevel structure to obtain an equivalent expanded graph. Although the expanded graph has a low diameter, a favorable property for constructing SSs, it has negatively weighted edges, which is an unfavorable property for the SSs. We design an algorithm to properly eliminate the negatively weighted edges and prove that the resulting expanded graph with positively weighted edges is spectrally equivalent to the expanded graph, thus, the original graph. Due to the low-diameter property of the positively-weighted expanded graph preconditioner (PEGP), existing algorithms for finding SSs can be easily applied. To demonstrate the advantage of working with the PEGP, we propose a type of SS, multilevel sparsifier preconditioner (MSP), that can be constructed in an easy and deterministic manner. We provide some preliminary numerical experiments to verify our theoretical findings and illustrate the practical effectiveness of PEGP and MSP in real-world applications.

preprint2021arXiv

An Efficient High-order Numerical Solver for Diffusion Equations with Strong Anisotropy

In this paper, we present an interior penalty discontinuous Galerkin finite element scheme for solving diffusion problems with strong anisotropy arising in magnetized plasmas for fusion applications. We demonstrate the accuracy produced by the high-order scheme and develop an efficient preconditioning technique to solve the corresponding linear system, which is robust to the mesh size and anisotropy of the problem. Several numerical tests are provided to validate the accuracy and efficiency of the proposed algorithm.

preprint2020arXiv

A Mesh-free Method Using Piecewise Deep Neural Network for Elliptic Interface Problems

In this paper, we propose a novel mesh-free numerical method for solving the elliptic interface problems based on deep learning. We approximate the solution by the neural networks and, since the solution may change dramatically across the interface, we employ different neural networks in different sub-domains. By reformulating the interface problem as a least-squares problem, we discretize the objective function using mean squared error via sampling and solve the proposed deep least-squares method by standard training algorithms such as stochastic gradient descent. The discretized objective function utilizes only the point-wise information on the sampling points and thus no underlying mesh is required. Doing this circumvents the challenging meshing procedure as well as the numerical integration on the complex interface. To improve the computational efficiency for more challenging problems, we further design an adaptive sampling strategy based on the residual of the least-squares function and propose an adaptive algorithm. Finally, we present several numerical experiments in both 2D and 3D to show the flexibility, effectiveness, and accuracy of the proposed deep least-square method for solving interface problems.

preprint2020arXiv

A Stabilized Hybrid Mixed Finite Element Method for Poroelasticity

In this work, we consider a hybrid mixed finite element method for Biot's model. The hybrid P1-RT0-P0 discretization of the displacement-pressure-Darcy's velocity system of Biot's model presented in \cite{C. Niu} is not uniformly stable with respect to the physical parameters, resulting in some issues in numerical simulations. To alleviate such problems, following \cite{V. Girault}, we stabilize the hybrid scheme with face bubble functions and show the well-posedness with respect to physical and discretization parameters, which provide optimal error estimates of the stabilized method. We introduce a perturbation of the bilinear form of the displacement which allows for the elimination of the bubble functions. Together with eliminating Darcy's velocity by hybridization, we obtain an eliminated system whose size is the same as the classical P1-RT0-P0 discretization. Based on the well-posedness of the eliminated system, we design block preconditioners that are parameter-robust. Numerical experiments are presented to confirm the theoretical results of the stabilized scheme as well as the block preconditioners.

preprint2020arXiv

Diffusion State Distances: Multitemporal Analysis, Fast Algorithms, and Applications to Biological Networks

Data-dependent metrics are powerful tools for learning the underlying structure of high-dimensional data. This article develops and analyzes a data-dependent metric known as diffusion state distance (DSD), which compares points using a data-driven diffusion process. Unlike related diffusion methods, DSDs incorporate information across time scales, which allows for the intrinsic data structure to be inferred in a parameter-free manner. This article develops a theory for DSD based on the multitemporal emergence of mesoscopic equilibria in the underlying diffusion process. New algorithms for denoising and dimension reduction with DSD are also proposed and analyzed. These approaches are based on a weighted spectral decomposition of the underlying diffusion process, and experiments on synthetic datasets and real biological networks illustrate the efficacy of the proposed algorithms in terms of both speed and accuracy. Throughout, comparisons with related methods are made, in order to illustrate the distinct advantages of DSD for datasets exhibiting multiscale structure.

preprint2020arXiv

Modeling, well-posedness and discretization for a class of models for mixed-dimensional problems with high dimensional gap

In this work, we illustrate the underlying mathematical structure of mixed-dimensional models arising from the composition of graphs and continuous domains. Such models are becoming popular in applications, in particular, to model the human vasculature. We first discuss the model equations in the strong form which describes the conservation of mass and Darcy's law in the continuum and network as well as the coupling between them. By introducing proper scaling, we propose a weak form that avoids degeneracy. Well-posedness of the weak form is shown through standard Babuška-Brezzi theory. We also develop the mixed formulation finite-element method and prove its well-posedness. A mass-lumping technique is introduced to derive the two-point flux approximation type discretization as well, due to its importance in applications. Based on the Babuška-Brezzi theory, error estimates can be obtained for both the finite-element scheme and the TPFA scheme. We also discuss efficient linear solvers for discrete problems. Finally, we present some numerical examples to verify the theoretical results and demonstrate the robustness of our proposed discretization schemes.

preprint2020arXiv

Randomized Fast Subspace Descent Methods

Randomized Fast Subspace Descent (RFASD) Methods are developed and analyzed for smooth and non-constraint convex optimization problems. The efficiency of the method relies on a space decomposition which is stable in $A$-norm, and meanwhile, the condition number $κ_A$ measured in $A$-norm is small. At each iteration, the subspace is chosen randomly either uniformly or by a probability proportional to the local Lipschitz constants. Then in each chosen subspace, a preconditioned gradient descent method is applied. RFASD converges sublinearly for convex functions and linearly for strongly convex functions. Comparing with the randomized block coordinate descent methods, the convergence of RFASD is faster provided $κ_A$ is small and the subspace decomposition is $A$-stable. This improvement is supported by considering a multilevel space decomposition for Nesterov's `worst' problem.

preprint2020arXiv

Robust preconditioners for a new stabilized discretization of the poroelastic equations

In this paper, we present block preconditioners for a stabilized discretization of the poroelastic equations developed in [45]. The discretization is proved to be well-posed with respect to the physical and discretization parameters, and thus provides a framework to develop preconditioners that are robust with respect to such parameters as well. We construct both norm-equivalent (diagonal) and field-of-value-equivalent (triangular) preconditioners for both the stabilized discretization and a perturbation of the stabilized discretization that leads to a smaller overall problem after static condensation. Numerical tests for both two- and three-dimensional problems confirm the robustness of the block preconditioners with respect to the physical and discretization parameters.

preprint2020arXiv

Using hierarchical matrices in the solution of the time-fractional heat equation by multigrid waveform relaxation

This work deals with the efficient numerical solution of the time-fractional heat equation discretized on non-uniform temporal meshes. Non-uniform grids are essential to capture the singularities of "typical" solutions of time-fractional problems. We propose an efficient space-time multigrid method based on the waveform relaxation technique, which accounts for the nonlocal character of the fractional differential operator. To maintain an optimal complexity, which can be obtained for the case of uniform grids, we approximate the coefficient matrix corresponding to the temporal discretization by its hierarchical matrix (${\cal H}$-matrix) representation. In particular, the proposed method has a computational cost of ${\cal O}(k N M \log(M))$, where $M$ is the number of time steps, $N$ is the number of spatial grid points, and $k$ is a parameter which controls the accuracy of the ${\cal H}$-matrix approximation. The efficiency and the good convergence of the algorithm, which can be theoretically justified by a semi-algebraic mode analysis, are demonstrated through numerical experiments in both one- and two-dimensional spaces.