-
On the equivalence problem of Smith forms for multivariate polynomial matrices
Authors:
Dong Lu,
Dingkang Wang,
Fanghui Xiao,
Xiaopeng Zheng
Abstract:
This paper delves into the equivalence problem of Smith forms for multivariate polynomial matrices. Generally speaking, multivariate ($n \geq 2$) polynomial matrices and their Smith forms may not be equivalent. However, under certain specific condition, we derive the necessary and sufficient condition for their equivalence. Let $F\in K[x_1,\ldots,x_n]^{l\times m}$ be of rank $r$,…
▽ More
This paper delves into the equivalence problem of Smith forms for multivariate polynomial matrices. Generally speaking, multivariate ($n \geq 2$) polynomial matrices and their Smith forms may not be equivalent. However, under certain specific condition, we derive the necessary and sufficient condition for their equivalence. Let $F\in K[x_1,\ldots,x_n]^{l\times m}$ be of rank $r$, $d_r(F)\in K[x_1]$ be the greatest common divisor of all the $r\times r$ minors of $F$, where $K$ is a field, $x_1,\ldots,x_n$ are variables and $1 \leq r \leq \min\{l,m\}$. Our key findings reveal the result: $F$ is equivalent to its Smith form if and only if all the $i\times i$ reduced minors of $F$ generate $K[x_1,\ldots,x_n]$ for $i=1,\ldots,r$.
△ Less
Submitted 9 July, 2024;
originally announced July 2024.
-
Local convergence analysis of L1/finite element scheme for a constant delay reaction-subdiffusion equation with uniform time mesh
Authors:
Weiping Bu,
Xin Zheng
Abstract:
The aim of this paper is to develop a refined error estimate of L1/finite element scheme for a reaction-subdiffusion equation with constant delay $τ$ and uniform time mesh. Under the non-uniform multi-singularity assumption of exact solution in time, the local truncation errors of the L1 scheme with uniform mesh is investigated. Then we introduce a fully discrete finite element scheme of the consi…
▽ More
The aim of this paper is to develop a refined error estimate of L1/finite element scheme for a reaction-subdiffusion equation with constant delay $τ$ and uniform time mesh. Under the non-uniform multi-singularity assumption of exact solution in time, the local truncation errors of the L1 scheme with uniform mesh is investigated. Then we introduce a fully discrete finite element scheme of the considered problem. Next, a novel discrete fractional Grönwall inequality with constant delay term is proposed, which does not include the increasing Mittag-Leffler function comparing with some popular other cases. By applying this Grönwall inequality, we obtain the pointwise-in-time and piecewise-in-time error estimates of the finite element scheme without the Mittag-Leffler function. In particular, the latter shows that, for the considered interval $((i-1)τ,iτ]$, although the convergence in time is low for $i=1$, it will be improved as the increasing $i$, which is consistent with the factual assumption that the smoothness of the solution will be improved as the increasing $i$. Finally, we present some numerical tests to verify the developed theory.
△ Less
Submitted 6 July, 2024;
originally announced July 2024.
-
New construction methods for uninorms via functions with q and uninorms on bounded lattices
Authors:
Zhen-Yu Xiu,
Zheng-Yuan Si,
Xu Zheng
Abstract:
In this paper, we study the construction methods for uninorms on bounded lattices via functions with the given uninorms and $q\in \mathbb{L_{B}}$ (or $p\in \mathbb{L_{B}}$). Specifically, we investigate the conditions under which these functions can be uninorms on bounded lattices when $q\in (0,\mathfrak{e})\cup I_{\mathfrak{e}}^{\varrho}$ and $q\in I_{\varrho}^{\mathfrak{e}}$ (or…
▽ More
In this paper, we study the construction methods for uninorms on bounded lattices via functions with the given uninorms and $q\in \mathbb{L_{B}}$ (or $p\in \mathbb{L_{B}}$). Specifically, we investigate the conditions under which these functions can be uninorms on bounded lattices when $q\in (0,\mathfrak{e})\cup I_{\mathfrak{e}}^{\varrho}$ and $q\in I_{\varrho}^{\mathfrak{e}}$ (or $p\in (\mathfrak{e},1)\cup I_{\mathfrak{e}}^σ$ and $p\in I_σ^{\mathfrak{e}}$), respectively. Moreover, some illustrative examples and figures are provided.
△ Less
Submitted 2 July, 2024;
originally announced July 2024.
-
Sub-Gaussian High-Dimensional Covariance Matrix Estimation under Elliptical Factor Model with 2 + εth Moment
Authors:
Yi Ding,
Xinghua Zheng
Abstract:
We study the estimation of high-dimensional covariance matrices under elliptical factor models with 2 + εth moment. For such heavy-tailed data, robust estimators like the Huber-type estimator in Fan, Liu and Wang (2018) can not achieve sub-Gaussian convergence rate. In this paper, we develop an idiosyncratic-projected self-normalization (IPSN) method to remove the effect of heavy-tailed scalar par…
▽ More
We study the estimation of high-dimensional covariance matrices under elliptical factor models with 2 + εth moment. For such heavy-tailed data, robust estimators like the Huber-type estimator in Fan, Liu and Wang (2018) can not achieve sub-Gaussian convergence rate. In this paper, we develop an idiosyncratic-projected self-normalization (IPSN) method to remove the effect of heavy-tailed scalar parameter, and propose a robust pilot estimator for the scatter matrix that achieves the sub-Gaussian rate. We further develop an estimator of the covariance matrix and show that it achieves a faster convergence rate than the generic POET estimator in Fan, Liu and Wang (2018).
△ Less
Submitted 26 June, 2024;
originally announced June 2024.
-
Numerical approximation for variable-exponent fractional diffusion-wave equation
Authors:
Xiangcheng Zheng,
Hong Wang,
Wenlin Qiu
Abstract:
This work considers the variable-exponent fractional diffusion-wave equation, which describes, e.g. the propagation of mechanical diffusive waves in viscoelastic media with varying material properties. Rigorous mathematical and numerical analysis for this model is not available in the literature, partly because the variable-exponent Abel kernel may not be positive definite or monotonic. We overcom…
▽ More
This work considers the variable-exponent fractional diffusion-wave equation, which describes, e.g. the propagation of mechanical diffusive waves in viscoelastic media with varying material properties. Rigorous mathematical and numerical analysis for this model is not available in the literature, partly because the variable-exponent Abel kernel may not be positive definite or monotonic. We overcome these difficulties to design two numerical schemes and derive their stability and error estimate based on the proved solution regularity, with $α(0)$-order and second-order accuracy in time, respectively. Numerical experiments are presented to substantiate the theoretical findings.
△ Less
Submitted 2 July, 2024; v1 submitted 5 June, 2024;
originally announced June 2024.
-
Weighted Sonine conditions and application
Authors:
Xiangcheng Zheng,
Shangqin Zhu,
Yiqun Li
Abstract:
The Sonine kernel described by the classical Sonine condition of convolution form is an important class of kernels used in integral equations and nonlocal differential equations. This work extends this idea to introduce weighted Sonine conditions where the non-convolutional weight functions accommodate the inhomogeneity in practical applications. We characterize tight relations between classical S…
▽ More
The Sonine kernel described by the classical Sonine condition of convolution form is an important class of kernels used in integral equations and nonlocal differential equations. This work extends this idea to introduce weighted Sonine conditions where the non-convolutional weight functions accommodate the inhomogeneity in practical applications. We characterize tight relations between classical Sonine condition and its weighted versions, which indicates that the non-degenerate weight functions may not introduce significant changes on the set of Sonine kernels. To demonstrate the application of weighted Sonine conditions, we employ them to derive equivalent but more feasible formulations of weighted integral equations and nonlocal differential equations to prove their well-posedness, and discuss possible application to corresponding partial differential equation models.
△ Less
Submitted 29 May, 2024;
originally announced May 2024.
-
Tests for principal eigenvalues and eigenvectors
Authors:
Jianqing Fan,
Yingying Li,
Ningning Xia,
Xinghua Zheng
Abstract:
We establish central limit theorems for principal eigenvalues and eigenvectors under a large factor model setting, and develop two-sample tests of both principal eigenvalues and principal eigenvectors. One important application is to detect structural breaks in large factor models. Compared with existing methods for detecting structural breaks, our tests provide unique insights into the source of…
▽ More
We establish central limit theorems for principal eigenvalues and eigenvectors under a large factor model setting, and develop two-sample tests of both principal eigenvalues and principal eigenvectors. One important application is to detect structural breaks in large factor models. Compared with existing methods for detecting structural breaks, our tests provide unique insights into the source of structural breaks because they can distinguish between individual principal eigenvalues and/or eigenvectors. We demonstrate the application by comparing the principal eigenvalues and principal eigenvectors of S\&P500 Index constituents' daily returns over different years.
△ Less
Submitted 11 May, 2024;
originally announced May 2024.
-
Two methods addressing variable-exponent fractional initial and boundary value problems and Abel integral equation
Authors:
Xiangcheng Zheng
Abstract:
Variable-exponent fractional models attract increasing attentions in various applications, while the rigorous analysis is far from well developed. This work provides general tools to address these models. Specifically, we first develop a convolution method to study the well-posedness, regularity, an inverse problem and numerical approximation for the sundiffusion of variable exponent. For models s…
▽ More
Variable-exponent fractional models attract increasing attentions in various applications, while the rigorous analysis is far from well developed. This work provides general tools to address these models. Specifically, we first develop a convolution method to study the well-posedness, regularity, an inverse problem and numerical approximation for the sundiffusion of variable exponent. For models such as the variable-exponent two-sided space-fractional boundary value problem (including the variable-exponent fractional Laplacian equation as a special case) and the distributed variable-exponent model, for which the convolution method does not apply, we develop a perturbation method to prove their well-posedness. The relation between the convolution method and the perturbation method is discussed, and we further apply the latter to prove the well-posedness of the variable-exponent Abel integral equation and discuss the constraint on the data under different initial values of variable exponent.
△ Less
Submitted 11 May, 2024; v1 submitted 14 April, 2024;
originally announced April 2024.
-
Further remarks on absorbing Markov decision processes
Authors:
Yi Zhang,
Xinran Zheng
Abstract:
In this note, based on the recent remarkable results of Dufour and Prieto-Rumeau, we deduce that for an absorbing MDP with a given initial state, under a standard compactness-continuity condition, the space of occupation measures has the same convergent sequences, when it is endowed with the weak topology and with the weak-strong topology. We provided two examples demonstrating that imposed condit…
▽ More
In this note, based on the recent remarkable results of Dufour and Prieto-Rumeau, we deduce that for an absorbing MDP with a given initial state, under a standard compactness-continuity condition, the space of occupation measures has the same convergent sequences, when it is endowed with the weak topology and with the weak-strong topology. We provided two examples demonstrating that imposed condition cannot be replaced with its popular alternative, and the above assertion does not hold for the space of marginals of occupation measures on the state space. Moreover, the examples also clarify some results in the previous literature.
△ Less
Submitted 29 March, 2024;
originally announced March 2024.
-
Addressing complex boundary conditions of miscible flow and transport in two and three dimensions with application to optimal control
Authors:
Yiqun Li,
Hong Wang,
Xiangcheng Zheng
Abstract:
We investigate complex boundary conditions of the miscible displacement system in two and three space dimensions with the commonly-used Bear-Scheidegger diffusion-dispersion tensor, which describes, e.g., the porous medium flow processes in petroleum reservoir simulation or groundwater contaminant transport. Specifically, we incorporate the no-flux boundary condition for the Darcy velocity to prov…
▽ More
We investigate complex boundary conditions of the miscible displacement system in two and three space dimensions with the commonly-used Bear-Scheidegger diffusion-dispersion tensor, which describes, e.g., the porous medium flow processes in petroleum reservoir simulation or groundwater contaminant transport. Specifically, we incorporate the no-flux boundary condition for the Darcy velocity to prove that the general no-flux boundary condition for the transport equation is equivalent to the normal derivative boundary condition of the concentration, based on which we further prove several complex boundary conditions by the Bear-Scheidegger tensor and its derivative. The derived boundary conditions not only provide new insights and distinct properties of the Bear-Scheidegger diffusion-dispersion tensor, but accommodate the coupling and the nonlinearity of the miscible displacement system and the Bear-Scheidegger tensor in deriving the first-order optimality condition of the corresponding optimal control problem for practical application.
△ Less
Submitted 20 March, 2024;
originally announced March 2024.
-
Understanding high-index saddle dynamics via numerical analysis
Authors:
Lei Zhang,
Pingwen Zhang,
Xiangcheng Zheng
Abstract:
High-index saddle dynamics (HiSD) serves as a competitive instrument in searching the any-index saddle points and constructing the solution landscape of complex systems. The Lagrangian multiplier terms in HiSD ensure the Stiefel manifold constraint, which, however, are dropped in the commonly-used discrete HiSD scheme and are replaced by an additional Gram-Schmidt orthonormalization. Though this s…
▽ More
High-index saddle dynamics (HiSD) serves as a competitive instrument in searching the any-index saddle points and constructing the solution landscape of complex systems. The Lagrangian multiplier terms in HiSD ensure the Stiefel manifold constraint, which, however, are dropped in the commonly-used discrete HiSD scheme and are replaced by an additional Gram-Schmidt orthonormalization. Though this scheme has been successfully applied in various fields, it is still unclear why the above modification does not affect its effectiveness. We recover the same form as HiSD from this scheme, which not only leads to error estimates naturally, but indicates that the mechanism of Stiefel manifold preservation by Lagrangian multiplier terms in HiSD is nearly a Gram-Schmidt process (such that the above modification is appropriate). The developed methods are further extended to analyze the more complicated constrained HiSD on high-dimensional sphere, which reveals more mechanisms of the constrained HiSD in preserving several manifold properties.
△ Less
Submitted 17 February, 2024;
originally announced February 2024.
-
A Priori Error Estimation of Physics-Informed Neural Networks Solving Allen--Cahn and Cahn--Hilliard Equations
Authors:
Guangtao Zhang,
Jiani Lin,
Qijia Zhai,
Huiyu Yang,
Xujun Chen,
Xiaoning Zheng,
Ieng Tak Leong
Abstract:
This paper aims to analyze errors in the implementation of the Physics-Informed Neural Network (PINN) for solving the Allen--Cahn (AC) and Cahn--Hilliard (CH) partial differential equations (PDEs). The accuracy of PINN is still challenged when dealing with strongly non-linear and higher-order time-varying PDEs. To address this issue, we introduce a stable and bounded self-adaptive weighting scheme…
▽ More
This paper aims to analyze errors in the implementation of the Physics-Informed Neural Network (PINN) for solving the Allen--Cahn (AC) and Cahn--Hilliard (CH) partial differential equations (PDEs). The accuracy of PINN is still challenged when dealing with strongly non-linear and higher-order time-varying PDEs. To address this issue, we introduce a stable and bounded self-adaptive weighting scheme known as Residuals-RAE, which ensures fair training and effectively captures the solution. By incorporating this new training loss function, we conduct numerical experiments on 1D and 2D AC and CH systems to validate our theoretical findings. Our theoretical analysis demonstrates that feedforward neural networks with two hidden layers and tanh activation function effectively bound the PINN approximation errors for the solution field, temporal derivative, and nonlinear term of the AC and CH equations by the training loss and number of collocation points.
△ Less
Submitted 4 February, 2024;
originally announced February 2024.
-
An equivalent formulation of Sonine condition
Authors:
Xiangcheng Zheng
Abstract:
Sonine kernel is characterized by the Sonine condition (denoted by SC) and is an important class of kernels in nonlocal differential equations and integral equations. This work proposes a SC with a more general form (denoted by gSC), which is more convenient than SC to accommodate complex kernels and equations. A typical kernel is given, and the first-kind Volterra integral equation under gSC is a…
▽ More
Sonine kernel is characterized by the Sonine condition (denoted by SC) and is an important class of kernels in nonlocal differential equations and integral equations. This work proposes a SC with a more general form (denoted by gSC), which is more convenient than SC to accommodate complex kernels and equations. A typical kernel is given, and the first-kind Volterra integral equation under gSC is accordingly transformed and then analyzed. Based on these results, it is finally proved that the gSC is indeed equivalent to the original SC, which indicates that the Sonine kernel may be essentially characterized by the behavior of its convolution with the associated kernel at the starting point.
△ Less
Submitted 13 March, 2024; v1 submitted 2 February, 2024;
originally announced February 2024.
-
Local modification of subdiffusion by initial Fickian diffusion: Multiscale modeling, analysis and computation
Authors:
Xiangcheng Zheng,
Yiqun Li,
Wenlin Qiu
Abstract:
We propose a local modification of the standard subdiffusion model by introducing the initial Fickian diffusion, which results in a multiscale diffusion model. The developed model resolves the incompatibility between the nonlocal operators in subdiffusion and the local initial conditions and thus eliminates the initial singularity of the solutions of the subdiffusion, while retaining its heavy tai…
▽ More
We propose a local modification of the standard subdiffusion model by introducing the initial Fickian diffusion, which results in a multiscale diffusion model. The developed model resolves the incompatibility between the nonlocal operators in subdiffusion and the local initial conditions and thus eliminates the initial singularity of the solutions of the subdiffusion, while retaining its heavy tail behavior away from the initial time. The well-posedness of the model and high-order regularity estimates of its solutions are analyzed by resolvent estimates, based on which the numerical discretization and analysis are performed. Numerical experiments are carried out to substantiate the theoretical findings.
△ Less
Submitted 30 January, 2024;
originally announced January 2024.
-
Zeroth-Order Non-Convex Optimization for Cooperative Multi-Agent Systems with Diminishing Step Size and Smoothing Radius
Authors:
Xinran Zheng,
Tara Javidi,
Behrouz Touri
Abstract:
We study a class of zeroth-order distributed optimization problems, where each agent can control a partial vector and observe a local cost that depends on the joint vector of all agents, and the agents can communicate with each other with time delay. We propose and study a gradient descent-based algorithm using two-point gradient estimators with diminishing smoothing parameters and diminishing ste…
▽ More
We study a class of zeroth-order distributed optimization problems, where each agent can control a partial vector and observe a local cost that depends on the joint vector of all agents, and the agents can communicate with each other with time delay. We propose and study a gradient descent-based algorithm using two-point gradient estimators with diminishing smoothing parameters and diminishing step-size and we establish the convergence rate to a first-order stationary point for general nonconvex problems. A byproduct of our proposed method with diminishing step size and smoothing parameters, as opposed to the fixed-parameter scheme, is that our proposed algorithm does not require any information regarding the local cost functions. This makes the solution appealing in practice as it allows for optimizing an unknown (black-box) global function without prior knowledge of its smoothness parameters. At the same time, the performance will adaptively match the problem instance parameters.
△ Less
Submitted 8 January, 2024;
originally announced January 2024.
-
Some further construction methods for uninorms via uninorms
Authors:
Zhenyu Xiu,
Xu Zheng
Abstract:
In this paper, we further investigate new construction methods for uninorms on bounded lattices via given uninorms. More specifically, we first construct new uninorms on arbitrary bounded lattices by extending a given uninorm on a subinterval of the lattices under necessary and sufficient conditions on the given uninorm. Moreover, based on the new uninorms, we can give another sufficient and neces…
▽ More
In this paper, we further investigate new construction methods for uninorms on bounded lattices via given uninorms. More specifically, we first construct new uninorms on arbitrary bounded lattices by extending a given uninorm on a subinterval of the lattices under necessary and sufficient conditions on the given uninorm. Moreover, based on the new uninorms, we can give another sufficient and necessary condition under which $S_{1}^{*}$ is a $t$-conorm ($T_{1}^{*}$ is a $t$-norm) in \cite{SS06}. Furthermore, using closure operators (interior operators), we also provide new construction methods for uninorms by extending the given uninorm on a subinterval of a bounded lattice under some additional constraints and simultaneously investigate the additional constraints carefully and systematically. Meanwhile, some illustrative examples for the above construction methods of uninorms on bounded lattices are provided.
△ Less
Submitted 3 December, 2023;
originally announced December 2023.
-
Accelerated high-index saddle dynamics method for searching high-index saddle points
Authors:
Yue Luo,
Xiangcheng Zheng,
Lei Zhang
Abstract:
The high-index saddle dynamics (HiSD) method [J. Yin, L. Zhang, and P. Zhang, {\it SIAM J. Sci. Comput., }41 (2019), pp.A3576-A3595] serves as an efficient tool for computing index-$k$ saddle points and constructing solution landscapes. Nevertheless, the conventional HiSD method often encounters slow convergence rates on ill-conditioned problems. To address this challenge, we propose an accelerate…
▽ More
The high-index saddle dynamics (HiSD) method [J. Yin, L. Zhang, and P. Zhang, {\it SIAM J. Sci. Comput., }41 (2019), pp.A3576-A3595] serves as an efficient tool for computing index-$k$ saddle points and constructing solution landscapes. Nevertheless, the conventional HiSD method often encounters slow convergence rates on ill-conditioned problems. To address this challenge, we propose an accelerated high-index saddle dynamics (A-HiSD) by incorporating the heavy ball method. We prove the linear stability theory of the continuous A-HiSD, and subsequently estimate the local convergence rate for the discrete A-HiSD. Our analysis demonstrates that the A-HiSD method exhibits a faster convergence rate compared to the conventional HiSD method, especially when dealing with ill-conditioned problems. We also perform various numerical experiments including the loss function of neural network to substantiate the effectiveness and acceleration of the A-HiSD method.
△ Less
Submitted 22 November, 2023; v1 submitted 6 October, 2023;
originally announced October 2023.
-
Semi-implicit method of high-index saddle dynamics and application to construct solution landscape
Authors:
Yue Luo,
Lei Zhang,
Pingwen Zhang,
Zhiyi Zhang,
Xiangcheng Zheng
Abstract:
We analyze the semi-implicit scheme of high-index saddle dynamics, which provides a powerful numerical method for finding the any-index saddle points and constructing the solution landscape. Compared with the explicit schemes of saddle dynamics, the semi-implicit discretization relaxes the step size and accelerates the convergence, but the corresponding numerical analysis encounters new difficulti…
▽ More
We analyze the semi-implicit scheme of high-index saddle dynamics, which provides a powerful numerical method for finding the any-index saddle points and constructing the solution landscape. Compared with the explicit schemes of saddle dynamics, the semi-implicit discretization relaxes the step size and accelerates the convergence, but the corresponding numerical analysis encounters new difficulties compared to the explicit scheme. Specifically, the orthonormal property of the eigenvectors at each time step could not be fully employed due to the semi-implicit treatment, and computations of the eigenvectors are coupled with the orthonormalization procedure, which further complicates the numerical analysis. We address these issues to prove error estimates of the semi-implicit scheme via, e.g. technical splittings and multi-variable circulating induction procedure. We further analyze the convergence rate of the generalized minimum residual solver for solving the semi-implicit system. Extensive numerical experiments are carried out to substantiate the efficiency and accuracy of the semi-implicit scheme in constructing solution landscapes of complex systems.
△ Less
Submitted 6 October, 2023;
originally announced October 2023.
-
Drinfeld Modular Curves Subordinate to Conjugacy Classes of Nilpotent Upper-Triangular Matrices
Authors:
Zhuo Chen,
Chuangqiang Hu,
Tao Zhang,
Xiaopeng Zheng
Abstract:
We introduce normalized Drinfeld modular curves that parameterize rank $m$ Drinfeld modules compatible with a $T$-torsion structure arising from a given conjugacy class of nilpotent upper-triangular $n\times n$ matrices with rank $\geqslant n-m$ over a finite field $\mathbb{F}_q$. This creates a deep link connecting the classification of nilpotent upper-triangular matrices and the decomposition of…
▽ More
We introduce normalized Drinfeld modular curves that parameterize rank $m$ Drinfeld modules compatible with a $T$-torsion structure arising from a given conjugacy class of nilpotent upper-triangular $n\times n$ matrices with rank $\geqslant n-m$ over a finite field $\mathbb{F}_q$. This creates a deep link connecting the classification of nilpotent upper-triangular matrices and the decomposition of Drinfeld modular curves. The conjugacy classes of nilpotent upper-triangular matrices one-to-one corresponds to certain $T$-torsion flags, and form a tree structure. As a result, the associated Drinfeld modular curves are organized in the same tree. This generalizes the tower structure introduced by Bassa, Beelen, Garcia, Stichtenoth, and others. Additionally,we prove the geometric irreducibility of $(3,2)$-type normalized Drinfeld modular curves, and characterize their associated function fields.
△ Less
Submitted 1 September, 2023;
originally announced September 2023.
-
Error estimate for semi-implicit method of sphere-constrained high-index saddle dynamics
Authors:
Lei Zhang,
Pingwen Zhang,
Xiangcheng Zheng
Abstract:
We prove error estimates for the semi-implicit numerical scheme of sphere-constrained high-index saddle dynamics, which serves as a powerful instrument in finding saddle points and constructing the solution landscapes of constrained systems on the high-dimensional sphere. Due to the semi-implicit treatment and the novel computational procedure, the orthonormality of numerical solutions at each tim…
▽ More
We prove error estimates for the semi-implicit numerical scheme of sphere-constrained high-index saddle dynamics, which serves as a powerful instrument in finding saddle points and constructing the solution landscapes of constrained systems on the high-dimensional sphere. Due to the semi-implicit treatment and the novel computational procedure, the orthonormality of numerical solutions at each time step could not be fully employed to simplify the derivations, and the computations of the state variable and directional vectors are coupled with the retraction, the vector transport and the orthonormalization procedure, which significantly complicates the analysis. We address these issues to prove error estimates for the proposed semi-implicit scheme and then carry out numerical experiments to substantiate the theoretical findings.
△ Less
Submitted 3 July, 2023;
originally announced July 2023.
-
Numerical algorithms and simulations of boundary dynamic control for optimal mixing in unsteady Stokes flows
Authors:
Xiaoming Zheng,
Weiwei Hu,
Jiahong Wu
Abstract:
This work develops an efficient and accurate optimization algorithm to study the optimal mixing problem driven by boundary control of unsteady Stokes flows, based on the theoretical foundation laid by Hu and Wu in a series of work. The scalar being mixed is purely advected by the flow and the control is a force exerted tangentially on the domain boundary through the Navier slip conditions. The con…
▽ More
This work develops an efficient and accurate optimization algorithm to study the optimal mixing problem driven by boundary control of unsteady Stokes flows, based on the theoretical foundation laid by Hu and Wu in a series of work. The scalar being mixed is purely advected by the flow and the control is a force exerted tangentially on the domain boundary through the Navier slip conditions. The control design has potential applications in many industrial processes such as rotating wall driven mixing, micromixers with acoustic waves, and artificial cilia mixing.
The numerical algorithms have high complexity, high accuracy demand, and high computing expense, due to the multiscale nature of the mixing problem and the optimization requirements. A crucial problem is the computation of the G$â$teaux derivative of the cost functional. To this end, a hybrid approach based on variational formula and finite difference is built with high accuracy and efficiency to treat various types of control input functions. We have experimented with various optimization schemes including the steepest descent algorithm, the conjugate gradient method and two line search options (backtracking and exact line search). We are able to identify and implement the best combinations.
The numerical simulations show that the mixing efficacy is limited when only one single type of control is applied, but can be enhanced when more diverse control types and more time segmentation are utilized. The mix-norm in the optimal mixings decays exponentially. The numerical study in this work demonstrates that boundary control alone could be an effective strategy for mixing in incompressible flows.
△ Less
Submitted 19 June, 2023;
originally announced June 2023.
-
A generalized fractional Laplacian
Authors:
Xiangcheng Zheng,
V. J. Ervin,
Hong Wang
Abstract:
In this article we show that the fractional Laplacian in $R^{2}$ can be factored into a product of the divergence operator, a Riesz potential operator, and the gradient operator. Using this factored form we introduce a generalization of the fractional Laplacian, involving a matrix $K(x)$, suitable when the fractional Laplacian is applied in a non homogeneous medium. For the case of $K(x)$ a consta…
▽ More
In this article we show that the fractional Laplacian in $R^{2}$ can be factored into a product of the divergence operator, a Riesz potential operator, and the gradient operator. Using this factored form we introduce a generalization of the fractional Laplacian, involving a matrix $K(x)$, suitable when the fractional Laplacian is applied in a non homogeneous medium. For the case of $K(x)$ a constant, symmetric positive definite matrix we show that the fractional Poisson equation is well posed, and determine the regularity of the solution in terms of the regularity of the right hand side function.
△ Less
Submitted 2 May, 2024; v1 submitted 24 April, 2023;
originally announced April 2023.
-
Iterative projection method for unsteady Navier-Stokes equations with high Reynolds numbers
Authors:
Xiaoming Zheng,
Kun Zhao,
Jiahong Wu,
Weiwei Hu,
Dapeng Du
Abstract:
A new approach, iteration projection method, is proposed to solve the saddle point problem obtained after the full discretization of the unsteady Navier-Stokes equations. The proposed method iterates projections in each time step with a proper convection form. We prove that the projection iterations converge with a certain parameter regime. Optimal iteration convergence can be achieved with the mo…
▽ More
A new approach, iteration projection method, is proposed to solve the saddle point problem obtained after the full discretization of the unsteady Navier-Stokes equations. The proposed method iterates projections in each time step with a proper convection form. We prove that the projection iterations converge with a certain parameter regime. Optimal iteration convergence can be achieved with the modulation of parameter values.
This new method has several significant improvements over the Uzawa method and the projection method both theoretically and practically. First, when the iterative projections are fully convergent in each time step, the numerical velocity is weakly divergence free (pointwise divergence free in divergence free finite element methods), and the stability and error estimate are rigorously proven. With proper parameters, this method converges much faster than the Uzawa method. Second, numerical simulations show that with rather relaxed stopping criteria which require only a few iterations each time step, the numerical solution preserves stability and accuracy for high Reynolds numbers, where the convectional projection method would fail. Furthermore, this method retains the efficiency of the traditional projection method by decoupling the velocity and pressure fields, which splits the saddle point system into small elliptic problems. Three dimensional simulations with Taylor-Hood P2/P1 finite elements are presented to demonstrate the performance and efficiency of this method.
More importantly, this method is a generic approach and thus there are many potential improvements and extensions of the iterative projection method, including utilization of various convection forms, association with stabilization techniques for high Reynolds numbers, applications to other saddle point problems.
△ Less
Submitted 16 April, 2023;
originally announced April 2023.
-
Asymptotic stability of homogeneous solutions to Navier-Stokes equations under $L^{p}$-perturbations
Authors:
Zhiwen Zhao,
Xiaoxin Zheng
Abstract:
It is known that there has been classified for all $(-1)$-homogeneous axisymmetric no-swirl solutions of the three-dimensional Navier-Stokes equations with a possible singular ray. The main purpose of this paper is to show that the least singular solutions among such solutions other than Landau solutions to the Navier-Stokes equations are asymptotically stable under $L^{3}$-perturbations. Moreover…
▽ More
It is known that there has been classified for all $(-1)$-homogeneous axisymmetric no-swirl solutions of the three-dimensional Navier-Stokes equations with a possible singular ray. The main purpose of this paper is to show that the least singular solutions among such solutions other than Landau solutions to the Navier-Stokes equations are asymptotically stable under $L^{3}$-perturbations. Moreover, we establish the $L^{q}$ decay estimate with an explicit decay rate and a sharp constant for any $q>3$. For that purpose, we first study the global well-posedness of solutions to the perturbed equations under small initial data in $L_σ^{3}$ space and the local well-posedness with any initial data in $L_σ^{p}$ spaces for $p\geq3$.
△ Less
Submitted 4 April, 2023; v1 submitted 3 April, 2023;
originally announced April 2023.
-
BINN: A deep learning approach for computational mechanics problems based on boundary integral equations
Authors:
Jia Sun,
Yinghua Liu,
Yizheng Wang,
Zhenhan Yao,
Xiaoping Zheng
Abstract:
We proposed the boundary-integral type neural networks (BINN) for the boundary value problems in computational mechanics. The boundary integral equations are employed to transfer all the unknowns to the boundary, then the unknowns are approximated using neural networks and solved through a training process. The loss function is chosen as the residuals of the boundary integral equations. Regulariza…
▽ More
We proposed the boundary-integral type neural networks (BINN) for the boundary value problems in computational mechanics. The boundary integral equations are employed to transfer all the unknowns to the boundary, then the unknowns are approximated using neural networks and solved through a training process. The loss function is chosen as the residuals of the boundary integral equations. Regularization techniques are adopted to efficiently evaluate the weakly singular and Cauchy principle integrals in boundary integral equations. Potential problems and elastostatic problems are mainly concerned in this article as a demonstration. The proposed method has several outstanding advantages: First, the dimensions of the original problem are reduced by one, thus the freedoms are greatly reduced. Second, the proposed method does not require any extra treatment to introduce the boundary conditions, since they are naturally considered through the boundary integral equations. Therefore, the method is suitable for complex geometries. Third, BINN is suitable for problems on the infinite or semi-infinite domains. Moreover, BINN can easily handle heterogeneous problems with a single neural network without domain decomposition.
△ Less
Submitted 11 January, 2023;
originally announced January 2023.
-
A model-free shrinking-dimer saddle dynamics for finding saddle point and solution landscape
Authors:
Lei Zhang,
Pingwen Zhang,
Xiangcheng Zheng
Abstract:
We propose a model-free shrinking-dimer saddle dynamics for finding any-index saddle points and constructing the solution landscapes, in which the force in the standard saddle dynamics is replaced by a surrogate model trained by the Gassian process learning. By this means, the exact form of the model is no longer necessary such that the saddle dynamics could be implemented based only on some obser…
▽ More
We propose a model-free shrinking-dimer saddle dynamics for finding any-index saddle points and constructing the solution landscapes, in which the force in the standard saddle dynamics is replaced by a surrogate model trained by the Gassian process learning. By this means, the exact form of the model is no longer necessary such that the saddle dynamics could be implemented based only on some observations of the force. This data-driven approach not only avoids the modeling procedure that could be difficult or inaccurate, but also significantly reduces the number of queries of the force that may be expensive or time-consuming. We accordingly develop a sequential learning saddle dynamics algorithm to perform a sequence of local saddle dynamics, in which the queries of the training samples and the update or retraining of the surrogate force are performed online and around the latent trajectory in order to improve the accuracy of the surrogate model and the value of each sampling. Numerical experiments are performed to demonstrate the effectiveness and efficiency of the proposed algorithm.
△ Less
Submitted 29 December, 2022;
originally announced December 2022.
-
High-dimensional covariance matrices under dynamic volatility models: asymptotics and shrinkage estimation
Authors:
Yi Ding,
Xinghua Zheng
Abstract:
We study the estimation of the high-dimensional covariance matrix andits eigenvalues under dynamic volatility models. Data under such modelshave nonlinear dependency both cross-sectionally and temporally. We firstinvestigate the empirical spectral distribution (ESD) of the sample covariancematrix under scalar BEKK models and establish conditions under which thelimiting spectral distribution (LSD)…
▽ More
We study the estimation of the high-dimensional covariance matrix andits eigenvalues under dynamic volatility models. Data under such modelshave nonlinear dependency both cross-sectionally and temporally. We firstinvestigate the empirical spectral distribution (ESD) of the sample covariancematrix under scalar BEKK models and establish conditions under which thelimiting spectral distribution (LSD) is either the same as or different fromthe i.i.d. case. We then propose a time-variation adjusted (TV-adj) sample co-variance matrix and prove that its LSD follows the same Marcenko-Pasturlaw as the i.i.d. case. Based on the asymptotics of the TV-adj sample co-variance matrix, we develop a consistent population spectrum estimator and an asymptotically optimal nonlinear shrinkage estimator of the unconditionalcovariance matrix
△ Less
Submitted 20 November, 2022; v1 submitted 18 November, 2022;
originally announced November 2022.
-
Disk potential functions for polygon spaces
Authors:
Yoosik Kim,
Siu-Cheong Lau,
Xiao Zheng
Abstract:
We derive a Floer theoretical SYZ mirror for an equilateral and generic polygon space. The disk potential function of the monotone torus fiber of the caterpillar bending system is calculated by computing non-trivial open Gromov--Witten invariants from the structural result of the monotone Fukaya category, the topology of fibers of completely integrable systems, and toric degenerations. Then, combi…
▽ More
We derive a Floer theoretical SYZ mirror for an equilateral and generic polygon space. The disk potential function of the monotone torus fiber of the caterpillar bending system is calculated by computing non-trivial open Gromov--Witten invariants from the structural result of the monotone Fukaya category, the topology of fibers of completely integrable systems, and toric degenerations. Then, combining the result with the work of Nohara--Ueda [NU20] and Marsh--Rietsch [MR20], we obtain the disk potential functions of bending systems and produce a mirror cluster variety of type A without frozen variables via Lagrangian Floer theory.
△ Less
Submitted 2 February, 2024; v1 submitted 7 November, 2022;
originally announced November 2022.
-
Approximate Pricing of Derivatives Under Fractional Stochastic Volatility Model
Authors:
Yuecai Han,
Xudong Zheng
Abstract:
We investigate the problem of pricing derivatives under a fractional stochastic volatility model. We obtain an approximate expression of the derivative price where the stochastic volatility can be composed of deterministic functions of time and fractional Ornstein-Uhlenbeck process. Numerical simulations are given to illustrate the feasibility and operability of the approximation, and also demonst…
▽ More
We investigate the problem of pricing derivatives under a fractional stochastic volatility model. We obtain an approximate expression of the derivative price where the stochastic volatility can be composed of deterministic functions of time and fractional Ornstein-Uhlenbeck process. Numerical simulations are given to illustrate the feasibility and operability of the approximation, and also demonstrate the effect of long-range on derivative prices.
△ Less
Submitted 20 October, 2022;
originally announced October 2022.
-
Identifiability Analysis of Noise Covariances for LTI Stochastic Systems with Unknown Inputs
Authors:
He Kong,
Salah Sukkarieh,
Travis J. Arnold,
Tianshi Chen,
Biqiang Mu,
Wei Xing Zheng
Abstract:
Most existing works on optimal filtering of linear time-invariant (LTI) stochastic systems with arbitrary unknown inputs assume perfect knowledge of the covariances of the noises in the filter design. This is impractical and raises the question of whether and under what conditions one can identify the process and measurement noise covariances (denoted as $Q$ and $R$, respectively) of systems with…
▽ More
Most existing works on optimal filtering of linear time-invariant (LTI) stochastic systems with arbitrary unknown inputs assume perfect knowledge of the covariances of the noises in the filter design. This is impractical and raises the question of whether and under what conditions one can identify the process and measurement noise covariances (denoted as $Q$ and $R$, respectively) of systems with unknown inputs. This paper considers the identifiability of $Q$/$R$ using the correlation-based measurement difference approach. More specifically, we establish (i) necessary conditions under which $Q$ and $R$ can be uniquely jointly identified; (ii) necessary and sufficient conditions under which $Q$ can be uniquely identified, when $R$ is known; (iii) necessary conditions under which $R$ can be uniquely identified, when $Q$ is known. It will also be shown that for achieving the results mentioned above, the measurement difference approach requires some decoupling conditions for constructing a stationary time series, which are proved to be sufficient for the well-known strong detectability requirements established by Hautus.
△ Less
Submitted 15 September, 2022;
originally announced September 2022.
-
Global regularity and decay behavior for Leray equations with critical-dissipation and Its Application to Self-similar Solutions
Authors:
Changxing Miao,
Xiaoxin Zheng
Abstract:
In this paper, we show the global regularity and the optimal decay of weak solutions to the generalized Leray problem with critical dissipation. Our method is based on the maximal smoothing effect, $L^{p}$-type elliptic regularity of linearization, and the action of the heat semigroup generated by the fractional powers of Laplace operator on distributions with Fourier transforms supported in an an…
▽ More
In this paper, we show the global regularity and the optimal decay of weak solutions to the generalized Leray problem with critical dissipation. Our method is based on the maximal smoothing effect, $L^{p}$-type elliptic regularity of linearization, and the action of the heat semigroup generated by the fractional powers of Laplace operator on distributions with Fourier transforms supported in an annulus. As a by-product, we shall construct a self-similar solution to the three-dimensional incompressible Navier-Stokes equations, and more importantly, prove the global regularity and the optimal decay without additional requirement of existing literatures.
△ Less
Submitted 31 August, 2022;
originally announced September 2022.
-
Distributed Algorithm Over Time-Varying Unbalanced Topologies for Optimization Problem Subject to Multiple Local Constraints
Authors:
Hongzhe Liu,
Wenwu Yu,
Guanghui Wen,
Wei Xing Zheng
Abstract:
This paper studies the distributed optimization problem with possibly nonidentical local constraints, where its global objective function is composed of $N$ convex functions. The aim is to solve the considered optimization problem in a distributed manner over time-varying unbalanced directed topologies by using only local information and performing only local computations. Towards this end, a new…
▽ More
This paper studies the distributed optimization problem with possibly nonidentical local constraints, where its global objective function is composed of $N$ convex functions. The aim is to solve the considered optimization problem in a distributed manner over time-varying unbalanced directed topologies by using only local information and performing only local computations. Towards this end, a new distributed discrete-time algorithm is developed by synthesizing the row stochastic matrices sequence and column stochastic matrices sequence analysis technique. Furthermore, for the developed distributed discrete-time algorithm, its convergence property to the optimal solution as well as its convergence rate are established under some mild assumptions. Numerical simulations are finally presented to verify the theoretical results.
△ Less
Submitted 25 August, 2022;
originally announced August 2022.
-
Mathematical and numerical analysis to shrinking-dimer saddle dynamics with local Lipschitz conditions
Authors:
Lei Zhang,
Pingwen Zhang,
Xiangcheng Zheng
Abstract:
We present a mathematical and numerical investigation to the shrinkingdimer saddle dynamics for finding any-index saddle points in the solution landscape. Due to the dimer approximation of Hessian in saddle dynamics, the local Lipschitz assumptions and the strong nonlinearity for the saddle dynamics, it remains challenges for delicate analysis, such as the the boundedness of the solutions and the…
▽ More
We present a mathematical and numerical investigation to the shrinkingdimer saddle dynamics for finding any-index saddle points in the solution landscape. Due to the dimer approximation of Hessian in saddle dynamics, the local Lipschitz assumptions and the strong nonlinearity for the saddle dynamics, it remains challenges for delicate analysis, such as the the boundedness of the solutions and the dimer error. We address these issues to bound the solutions under proper relaxation parameters, based on which we prove the error estimates for numerical discretization to the shrinking-dimer saddle dynamics by matching the dimer length and the time step size. Furthermore, the Richardson extrapolation is employed to obtain a high-order approximation.
The inherent reason of requiring the matching of the dimer length and the time step size lies in that the former serves a different mesh size from the later, and thus the proposed numerical method is close to a fully-discrete numerical scheme of some spacetime PDE model with the Hessian in the saddle dynamics and its dimer approximation serving as a "spatial operator" and its discretization, respectively, which in turn indicates the PDE nature of the saddle dynamics.
△ Less
Submitted 20 July, 2022;
originally announced July 2022.
-
Discretization and index-robust error analysis for constrained high-index saddle dynamics on high-dimensional sphere
Authors:
Lei Zhang,
Pingwen Zhang,
Xiangcheng Zheng
Abstract:
We develop and analyze numerical discretization to the constrained high-index saddle dynamics, the dynamics searching for the high-index saddle points confined on the high-dimensional unit sphere. Compared with the saddle dynamics without constraints, the constrained high-index saddle dynamics has more complex dynamical forms, and additional operations such as the retraction and vector transport a…
▽ More
We develop and analyze numerical discretization to the constrained high-index saddle dynamics, the dynamics searching for the high-index saddle points confined on the high-dimensional unit sphere. Compared with the saddle dynamics without constraints, the constrained high-index saddle dynamics has more complex dynamical forms, and additional operations such as the retraction and vector transport are required due to the constraint, which significantly complicate the numerical scheme and the corresponding numerical analysis. Furthermore, as the existing numerical analysis results usually depend on the index of the saddle points implicitly, the proved numerical accuracy may be reduced if the index is high in many applications, which indicates the lack of robustness with respect to the index. To address these issues, we derive the error estimates for numerical discretization of the constrained high-index saddle dynamics on high-dimensional sphere, and then improve it by providing an index-robust error analysis in an averaged norm by adjusting the relaxation parameters. The developed results provide mathematical supports for the accuracy of numerical computations.
△ Less
Submitted 21 June, 2022;
originally announced June 2022.
-
Elliptic Hypertoric Varieties
Authors:
Naichung Conan Leung,
Xiao Zheng
Abstract:
We introduce elliptic hypertoric varieties, which is an elliptic analogue of hypertoric varieties and multiplicative hypertoric varieties. We also prove an elliptic version of Hikita conjecture, which relates elliptic (resp. additive and mulitplicative) hypertoric varieties to the equivariant elliptic cohomology (resp. ordinary cohomology and K-theory) of their 3d mirror hypertoric varieties.
We introduce elliptic hypertoric varieties, which is an elliptic analogue of hypertoric varieties and multiplicative hypertoric varieties. We also prove an elliptic version of Hikita conjecture, which relates elliptic (resp. additive and mulitplicative) hypertoric varieties to the equivariant elliptic cohomology (resp. ordinary cohomology and K-theory) of their 3d mirror hypertoric varieties.
△ Less
Submitted 26 April, 2022;
originally announced April 2022.
-
Convergence analysis of discrete high-index saddle dynamics
Authors:
Yue Luo,
Xiangcheng Zheng,
Xiangle Cheng,
Lei Zhang
Abstract:
Saddle dynamics is a time continuous dynamics to efficiently compute the any-index saddle points and construct the solution landscape. In practice, the saddle dynamics needs to be discretized for numerical computations, while the corresponding numerical analysis are rarely studied in the literature, especially for the high-index cases. In this paper we propose the convergence analysis of discrete…
▽ More
Saddle dynamics is a time continuous dynamics to efficiently compute the any-index saddle points and construct the solution landscape. In practice, the saddle dynamics needs to be discretized for numerical computations, while the corresponding numerical analysis are rarely studied in the literature, especially for the high-index cases. In this paper we propose the convergence analysis of discrete high-index saddle dynamics. To be specific, we prove the local linear convergence rates of numerical schemes of high-index saddle dynamics, which indicates that the local curvature in the neighborhood of the saddle point and the accuracy of computing the eigenfunctions are main factors that affect the convergence of discrete saddle dynamics. The proved results serve as compensations for the convergence analysis of high-index saddle dynamics and are substantiated by numerical experiments.
△ Less
Submitted 31 August, 2022; v1 submitted 31 March, 2022;
originally announced April 2022.
-
Analysis and Petrov-Galerkin numerical approximation for variable coefficient two-sided fractional diffusion, advection, reaction equations
Authors:
Xiangcheng Zheng,
V. J. Ervin,
Hong Wang
Abstract:
In this paper we investigate the variable coefficient two-sided fractional diffusion, advection, reaction equations on a bounded interval. It is known that the fractional diffusion operator may lose coercivity due to the variable coefficient, which makes both the mathematical and numerical analysis challenging. To resolve this issue, we design appropriate test and trial functions to prove the inf-…
▽ More
In this paper we investigate the variable coefficient two-sided fractional diffusion, advection, reaction equations on a bounded interval. It is known that the fractional diffusion operator may lose coercivity due to the variable coefficient, which makes both the mathematical and numerical analysis challenging. To resolve this issue, we design appropriate test and trial functions to prove the inf-sup condition of the variable coefficient fractional diffusion, advection, reaction operators in suitable function spaces. Based on this property, we prove the well-posedness and regularity of the solutions, as well as analyze the Petrov-Galerkin approximation scheme for the proposed model. Numerical experiments are presented to substantiate the theoretical findings and to compare the behaviors of different models.
△ Less
Submitted 22 March, 2022;
originally announced March 2022.
-
Moment properties for two-type continuous-state branching processes in random environments
Authors:
Shukai Chen,
Xiangqi Zheng
Abstract:
We first derive the recurisions for integer moments of two-type continuous-state branching processes in Lévy random environments. Result shows that the $n$th moment of the process is a polynomial of the initial value of the process with at most $n$ degree. Under some natural condition, the criteria for the existence of $f$-moment of the process are also proved.
We first derive the recurisions for integer moments of two-type continuous-state branching processes in Lévy random environments. Result shows that the $n$th moment of the process is a polynomial of the initial value of the process with at most $n$ degree. Under some natural condition, the criteria for the existence of $f$-moment of the process are also proved.
△ Less
Submitted 4 August, 2022; v1 submitted 11 March, 2022;
originally announced March 2022.
-
PERCEPT: a new online change-point detection method using topological data analysis
Authors:
Xiaojun Zheng,
Simon Mak,
Liyan Xie,
Yao Xie
Abstract:
Topological data analysis (TDA) provides a set of data analysis tools for extracting embedded topological structures from complex high-dimensional datasets. In recent years, TDA has been a rapidly growing field which has found success in a wide range of applications, including signal processing, neuroscience and network analysis. In these applications, the online detection of changes is of crucial…
▽ More
Topological data analysis (TDA) provides a set of data analysis tools for extracting embedded topological structures from complex high-dimensional datasets. In recent years, TDA has been a rapidly growing field which has found success in a wide range of applications, including signal processing, neuroscience and network analysis. In these applications, the online detection of changes is of crucial importance, but this can be highly challenging since such changes often occur in a low-dimensional embedding within high-dimensional data streams. We thus propose a new method, called PERsistence diagram-based ChangE-PoinT detection (PERCEPT), which leverages the learned topological structure from TDA to sequentially detect changes. PERCEPT follows two key steps: it first learns the embedded topology as a point cloud via persistence diagrams, then applies a non-parametric monitoring approach for detecting changes in the resulting point cloud distributions. This yields a non-parametric, topology-aware framework which can efficiently detect online changes from high-dimensional data streams. We investigate the effectiveness of PERCEPT over existing methods in a suite of numerical experiments where the data streams have an embedded topological structure. We then demonstrate the usefulness of PERCEPT in two applications in solar flare monitoring and human gesture detection.
△ Less
Submitted 8 March, 2022;
originally announced March 2022.
-
Wasserstein-type distances of two-type continuous-state branching processes in Lévy random environments
Authors:
Shukai Chen,
Rongjuan Fang,
Xiangqi Zheng
Abstract:
Under natural conditions, we proved the exponential ergodicity in Wasserstein distance of two-type continuous-state branching processes in Lévy random environments with immigration. Furthermore, we expressed accurately the parameters of the exponent. The coupling method and the conditioned branching property play an important role in the approach. Using the tool of superprocesses, the ergodicity i…
▽ More
Under natural conditions, we proved the exponential ergodicity in Wasserstein distance of two-type continuous-state branching processes in Lévy random environments with immigration. Furthermore, we expressed accurately the parameters of the exponent. The coupling method and the conditioned branching property play an important role in the approach. Using the tool of superprocesses, the ergodicity in total variance distance is also proved.
△ Less
Submitted 11 February, 2022;
originally announced February 2022.
-
The Noise Covariances of Linear Gaussian Systems with Unknown Inputs Are Not Uniquely Identifiable Using Autocovariance Least-squares
Authors:
He Kong,
Salah Sukkarieh,
Travis J. Arnold,
Tianshi Chen,
Wei Xing Zheng
Abstract:
Existing works in optimal filtering for linear Gaussian systems with arbitrary unknown inputs assume perfect knowledge of the noise covariances in the filter design. This is impractical and raises the question of whether and under what conditions one can identify the noise covariances of linear Gaussian systems with arbitrary unknown inputs. This paper considers the above identifiability question…
▽ More
Existing works in optimal filtering for linear Gaussian systems with arbitrary unknown inputs assume perfect knowledge of the noise covariances in the filter design. This is impractical and raises the question of whether and under what conditions one can identify the noise covariances of linear Gaussian systems with arbitrary unknown inputs. This paper considers the above identifiability question using the correlation-based autocovariance least-squares (ALS) approach. In particular, for the ALS framework, we prove that (i) the process noise covariance Q and the measurement noise covariance R cannot be uniquely jointly identified; (ii) neither Q nor R is uniquely identifiable, when the other is known. This not only helps us to have a better understanding of the applicability of existing filtering frameworks under unknown inputs (since almost all of them require perfect knowledge of the noise covariances) but also calls for further investigation of alternative and more viable noise covariance methods under unknown inputs. Especially, it remains to be explored whether the noise covariances are uniquely identifiable using other correlation-based methods. We are also interested to use regularization for noise covariance estimation under unknown inputs, and investigate the relevant property guarantees for the covariance estimates. The above topics are the main subject of our current and future work.
△ Less
Submitted 10 February, 2022;
originally announced February 2022.
-
Quasi-Framelets: Another Improvement to GraphNeural Networks
Authors:
Mengxi Yang,
Xuebin Zheng,
Jie Yin,
Junbin Gao
Abstract:
This paper aims to provide a novel design of a multiscale framelets convolution for spectral graph neural networks. In the spectral paradigm, spectral GNNs improve graph learning task performance via proposing various spectral filters in spectral domain to capture both global and local graph structure information. Although the existing spectral approaches show superior performance in some graphs,…
▽ More
This paper aims to provide a novel design of a multiscale framelets convolution for spectral graph neural networks. In the spectral paradigm, spectral GNNs improve graph learning task performance via proposing various spectral filters in spectral domain to capture both global and local graph structure information. Although the existing spectral approaches show superior performance in some graphs, they suffer from lack of flexibility and being fragile when graph information are incomplete or perturbated. Our new framelets convolution incorporates the filtering func-tions directly designed in the spectral domain to overcome these limitations. The proposed convolution shows a great flexibility in cutting-off spectral information and effectively mitigate the negative effect of noisy graph signals. Besides, to exploit the heterogeneity in real-world graph data, the heterogeneous graph neural network with our new framelet convolution provides a solution for embedding the intrinsic topological information of meta-path with a multi-level graph analysis.Extensive experiments have been conducted on real-world heterogeneous graphs and homogeneous graphs under settings with noisy node features and superior performance results are achieved.
△ Less
Submitted 10 January, 2022;
originally announced January 2022.
-
Global-in-time Boundedness of solution for Cauchy problem to the Parabolic-Parabolic Keller-Segel system with logistic growth
Authors:
Yao Nie,
Xiaoxin Zheng
Abstract:
We study global-in-time well-posedness and the behaviour and of the solution to Cauchy problem in the classical Keller-Segel system with logistic term \begin{equation*} \left. \aligned
\partial_tn-Δn=&-χ\nabla\cdot(n\nabla c)+\la n-μn^2
τ\partial_tc-Δc=&-c+n \endaligned \right\}\quad\text{in}\,\,\,\RR^d\times\RR^+, \end{equation*} where $d\ge 1$, $τ,\, χ,\, μ>0$ and $λ\ge 0$. It's inspired by…
▽ More
We study global-in-time well-posedness and the behaviour and of the solution to Cauchy problem in the classical Keller-Segel system with logistic term \begin{equation*} \left. \aligned
\partial_tn-Δn=&-χ\nabla\cdot(n\nabla c)+\la n-μn^2
τ\partial_tc-Δc=&-c+n \endaligned \right\}\quad\text{in}\,\,\,\RR^d\times\RR^+, \end{equation*} where $d\ge 1$, $τ,\, χ,\, μ>0$ and $λ\ge 0$. It's inspired by a previous result \cite[M. Winkler, Commun. Part. Diff. Eq., 35 (2010), 1516-1537]{Win10}, where the global-in-time boundedness of the above Keller-Segel system in smooth \emph{bounded }convex domains is established for large $μ$. However, his approach in bounded domain ceases to directly apply in the entire space $\RR^d$, and then they raised an interesting question whether a similar global-in-time boundedness statement remains true of Cauchy problem. In this paper, we answer this open problem by developing local-in-space estimates. More precisely, we prove that the above Keller-Segel system possesses a uniquely global-in-time bounded solution for any $τ>0$ under the assumption that $μ$ is large. The key point of our proof heavily relies on localization in space of solution caused by "local effect" of $L^\infty(\RR^d)$-norm.
△ Less
Submitted 5 January, 2022;
originally announced January 2022.
-
Supercritical Spatial SIR Epidemics: Spreading Speed and Herd Immunity
Authors:
Xinghua Zheng,
Qingsan Zhu
Abstract:
We study supercritical spatial SIR epidemics on $\mathbb{Z}^2\times \{1,2,\ldots, N\}$, where each site in $\mathbb{Z}^2$ represents a village and $N$ stands for the village size. We establish several key asymptotic results as $N\to\infty$. In particular, we derive the probability that the epidemic will last forever if the epidemic is started by one infected individual. Moreover, conditional on th…
▽ More
We study supercritical spatial SIR epidemics on $\mathbb{Z}^2\times \{1,2,\ldots, N\}$, where each site in $\mathbb{Z}^2$ represents a village and $N$ stands for the village size. We establish several key asymptotic results as $N\to\infty$. In particular, we derive the probability that the epidemic will last forever if the epidemic is started by one infected individual. Moreover, conditional on that the epidemic lasts forever, we show that the epidemic spreads out linearly in all directions and derive an explicit formula for the spreading speed. Furthermore, we prove that the ultimate proportion of infection converges to a number that is constant over space and find its explicit value. An important message is that if there is no vaccination, then the ultimate proportion of population who will be infected can be \emph{much higher} than the vaccination proportion that is needed in order to prevent sustained spread of the infection.
△ Less
Submitted 29 November, 2021;
originally announced November 2021.
-
Error estimates for Euler discretization of high-index saddle dynamics
Authors:
Lei Zhang,
Pingwen Zhang,
Xiangcheng Zheng
Abstract:
High-index saddle dynamics provides an effective means to compute the any-index saddle points and construct the solution landscape. In this paper we prove error estimates for Euler discretization of high-index saddle dynamics with respect to the time step size, which remains untreated in the literature. We overcome the main difficulties that lie in the strong nonlinearity of the saddle dynamics an…
▽ More
High-index saddle dynamics provides an effective means to compute the any-index saddle points and construct the solution landscape. In this paper we prove error estimates for Euler discretization of high-index saddle dynamics with respect to the time step size, which remains untreated in the literature. We overcome the main difficulties that lie in the strong nonlinearity of the saddle dynamics and the orthonormalization procedure in the numerical scheme that is uncommon in standard discretization of differential equations. The derived methods are further extended to study the generalized high-index saddle dynamics for non-gradient systems and provide theoretical support for the accuracy of numerical implementations.
△ Less
Submitted 4 August, 2022; v1 submitted 9 November, 2021;
originally announced November 2021.
-
Numerical approximation for a nonlinear variable-order fractional differential equation via an integral equation method
Authors:
Xiangcheng Zheng
Abstract:
We study a numerical approximation for a nonlinear variable-order fractional differential equation via an integral equation method. Due to the lack of the monotonicity of the discretization coefficients of the variable-order fractional derivative in standard approximation schemes, existing numerical analysis techniques do not apply directly. By an approximate inversion technique, the proposed mode…
▽ More
We study a numerical approximation for a nonlinear variable-order fractional differential equation via an integral equation method. Due to the lack of the monotonicity of the discretization coefficients of the variable-order fractional derivative in standard approximation schemes, existing numerical analysis techniques do not apply directly. By an approximate inversion technique, the proposed model is transformed as a second kind Volterra integral equation, based on which a collocation method under uniform or graded mesh is developed and analyzed. In particular, the error estimates improve the existing results by proving a consistent and sharper mesh grading parameter and characterizing the convergence rates in terms of the initial value of the variable order, which demonstrates its critical role in determining the smoothness of the solutions and thus the numerical accuracy.
△ Less
Submitted 10 October, 2021;
originally announced October 2021.
-
Approximate inversion for Abel integral operators of variable exponent and applications to fractional Cauchy problems
Authors:
Xiangcheng Zheng
Abstract:
We investigate the variable-exponent Abel integral equations and corresponding fractional Cauchy problems. The main contributions of the work are enumerated as follows: (i) We develop an approximate inversion technique for variable-exponent Abel integral operators, based on which we analyze the corresponding integral and differential equations; (ii) We prove that the sensitive dependence of the we…
▽ More
We investigate the variable-exponent Abel integral equations and corresponding fractional Cauchy problems. The main contributions of the work are enumerated as follows: (i) We develop an approximate inversion technique for variable-exponent Abel integral operators, based on which we analyze the corresponding integral and differential equations; (ii) We prove that the sensitive dependence of the well-posedness of classical Riemann-Liouville fractional differential equations on the initial value could be resolved by adjusting the initial value of the variable exponent; (iii) We prove that the singularity of the solutions to the Riemann-Liouville fractional differential equations could also be eliminated by adjusting the variable exponent and its derivatives at the initial time, which, together with (ii), demonstrates the advantages of introducing the variable exponent. The proposed approximate inversion technique provides a potential means to convert the intricate variable-exponent integral and differential problems to feasible forms, and the above findings suggest that the variable-exponent fractional problems may serve as a connection between integer-order and fractional models by adjusting the variable exponent at the initial time.
△ Less
Submitted 10 October, 2021; v1 submitted 2 October, 2021;
originally announced October 2021.
-
Evaluation of mechanical and energy properties for the phase field modeling of failure
Authors:
Yuanfeng Yu,
Xiaoya Zheng,
Peng Li,
Zhongzhou Zhang,
Jinyou Xiao
Abstract:
In recent years, various phase field models have been developed in variational methods to simulate the failure of brittle solids. However, there is a lack of objective evaluation of the existing results, and in particular, there are few studies on model nonhomogeneous resolution, stress-strain linear elastic properties, and failure stress estimation. To compensate for the above gaps, the commonly…
▽ More
In recent years, various phase field models have been developed in variational methods to simulate the failure of brittle solids. However, there is a lack of objective evaluation of the existing results, and in particular, there are few studies on model nonhomogeneous resolution, stress-strain linear elastic properties, and failure stress estimation. To compensate for the above gaps, the commonly used variational phase field model is systematically analyzed to solve the problem of evaluating the mechanics and energy properties of the model in this paper. The unified expression of the analytical solution and the nonhomogeneous solution under specific boundary conditions is analyzed and verified. Additionally, we theoretically analyze the energy properties of the phase field model and study the influence of the critical strain energy on the damage field, stress and strain of different models. Finally, the effect of different effective material parameters on the material failure stress is analyzed, the stress estimation formula for the failure of mode I of various phase field models is given, and the accuracy of the theory is verified by some examples.
△ Less
Submitted 22 September, 2021;
originally announced September 2021.
-
Analysis of viscoelastic flow with a generalized memory and its exponential convergence to steady state
Authors:
Yingwen Guo,
Xiangcheng Zheng
Abstract:
We investigate a viscoelastic flow model with a generalized memory, in which a weak-singular component is introduced in the exponential convolution kernel of classical viscoelastic flow equations that remains untreated in the literature. We prove the well-posedness and regularity of the solutions, based on which we prove the exponential convergence of the solutions to the steady state. The propose…
▽ More
We investigate a viscoelastic flow model with a generalized memory, in which a weak-singular component is introduced in the exponential convolution kernel of classical viscoelastic flow equations that remains untreated in the literature. We prove the well-posedness and regularity of the solutions, based on which we prove the exponential convergence of the solutions to the steady state. The proposed model serves as an extension of classical viscoelastic flow equations by adding a dimension characterized by the power of the weak-singular kernel, and the derived results provide theoretical supports for designing numerical methods for both the considered equation and its steady state.
△ Less
Submitted 1 October, 2021; v1 submitted 1 September, 2021;
originally announced September 2021.
-
Boundary control for optimal mixing via Stokes flows and numerical implementation
Authors:
Weiwei Hu,
Xiaoming Zheng
Abstract:
This work develops scientific computing techniques to further the exploration of using boundary control alone to optimize mixing in Stokes flows. The theoretical foundation including mathematical model and the optimality conditions for solving the optimal control has been established by Hu and Wu in a series of work. The scalar being mixed is purely advected by the flow and the control is exerted…
▽ More
This work develops scientific computing techniques to further the exploration of using boundary control alone to optimize mixing in Stokes flows. The theoretical foundation including mathematical model and the optimality conditions for solving the optimal control has been established by Hu and Wu in a series of work. The scalar being mixed is purely advected by the flow and the control is exerted tangentially on the domain boundary through the Navier slip conditions. The control design is motivated by the physical observations that the moving or rotating walls accelerate mixing.
A gradient descent-based optimization algorithm is designed. A critical problem is the computation of the Gateaux derivative or the gradient of the cost functional. Two methods are proposed: one is based on the Variational Formula (VF) and one utilizes Algorithmic Differentiation (AD). The convergence of the algorithm is studied and various designs of boundary control using cosine and sine functions with time segmentation are computed.
The algorithm has a first order convergence rate and the VF method is more efficient by taking only one third of the time as the AD method when the dimension of control basis is large. The numerical implementations show that the boundary control produces similar mixing results as internal mixings in the existing literature. The mixing effect becomes better when more diverse basis control functions and more time segmentation are utilized. It is shown that the mixing decay rate in time follows power rules, approximately. The numerical study in this work suggests that boundary control alone could be an effective strategy for mixing in incompressible fluid flows.
△ Less
Submitted 21 February, 2024; v1 submitted 21 August, 2021;
originally announced August 2021.