Showing 1–50 of 361 results for author: Wang, F

On Green's function of the vorticity formulation for the 3D Navier-Stokes equations

Authors: Igor Kukavica, Fei Wang, Yichun Zhu

Abstract: We give a novel vorticity formulation for the 3D Navier-Stokes equations with Dirichlet boundary conditions. Via a resolvent argument, we obtain Green's function and establish an upper bound, which is the 3D analog of [24]. Moreover, we prove similar results for the corresponding Stokes problem with more general mixed boundary conditions. We give a novel vorticity formulation for the 3D Navier-Stokes equations with Dirichlet boundary conditions. Via a resolvent argument, we obtain Green's function and establish an upper bound, which is the 3D analog of [24]. Moreover, we prove similar results for the corresponding Stokes problem with more general mixed boundary conditions. △ Less

Submitted 15 July, 2024; originally announced July 2024.

The Planar Turán Number of $Θ_6$-graphs

Authors: David Guan, Ervin Győri, Diep Luong-Le, Felicia Wang, Mengyuan Yang

Abstract: There are two particular $Θ_6$-graphs - the 6-cycle graphs with a diagonal. We find the planar Turán number of each of them, i.e. the maximum number of edges in a planar graph $G$ of $n$ vertices not containing the given $Θ_6$ as a subgraph and we find infinitely many extremal constructions showing the sharpness of these results - apart from a small additive constant error in one of the cases. There are two particular $Θ_6$-graphs - the 6-cycle graphs with a diagonal. We find the planar Turán number of each of them, i.e. the maximum number of edges in a planar graph $G$ of $n$ vertices not containing the given $Θ_6$ as a subgraph and we find infinitely many extremal constructions showing the sharpness of these results - apart from a small additive constant error in one of the cases. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Profiled Transfer Learning for High Dimensional Linear Model

Authors: Ziqian Lin, Junlong Zhao, Fang Wang, Hansheng Wang

Abstract: We develop here a novel transfer learning methodology called Profiled Transfer Learning (PTL). The method is based on the \textit{approximate-linear} assumption between the source and target parameters. Compared with the commonly assumed \textit{vanishing-difference} assumption and \textit{low-rank} assumption in the literature, the \textit{approximate-linear} assumption is more flexible and less… ▽ More We develop here a novel transfer learning methodology called Profiled Transfer Learning (PTL). The method is based on the \textit{approximate-linear} assumption between the source and target parameters. Compared with the commonly assumed \textit{vanishing-difference} assumption and \textit{low-rank} assumption in the literature, the \textit{approximate-linear} assumption is more flexible and less stringent. Specifically, the PTL estimator is constructed by two major steps. Firstly, we regress the response on the transferred feature, leading to the profiled responses. Subsequently, we learn the regression relationship between profiled responses and the covariates on the target data. The final estimator is then assembled based on the \textit{approximate-linear} relationship. To theoretically support the PTL estimator, we derive the non-asymptotic upper bound and minimax lower bound. We find that the PTL estimator is minimax optimal under appropriate regularity conditions. Extensive simulation studies are presented to demonstrate the finite sample performance of the new method. A real data example about sentence prediction is also presented with very encouraging results. △ Less

Submitted 5 June, 2024; v1 submitted 2 June, 2024; originally announced June 2024.

Uniform Inviscid Damping and Inviscid Limit of the 2D Navier-Stokes equation with Navier Boundary Conditions

Authors: Jacob Bedrossian, Siming He, Sameer Iyer, Fei Wang

Abstract: We consider the 2D, incompressible Navier-Stokes equations near the Couette flow, $ω^{(NS)} = 1 + εω$, set on the channel $\mathbb{T} \times [-1, 1]$, supplemented with Navier boundary conditions on the perturbation, $ω|_{y = \pm 1} = 0$. We are simultaneously interested in two asymptotic regimes that are classical in hydrodynamic stability: the long time, $t \rightarrow \infty$, stability of back… ▽ More We consider the 2D, incompressible Navier-Stokes equations near the Couette flow, $ω^{(NS)} = 1 + εω$, set on the channel $\mathbb{T} \times [-1, 1]$, supplemented with Navier boundary conditions on the perturbation, $ω|_{y = \pm 1} = 0$. We are simultaneously interested in two asymptotic regimes that are classical in hydrodynamic stability: the long time, $t \rightarrow \infty$, stability of background shear flows, and the inviscid limit, $ν\rightarrow 0$ in the presence of boundaries. Given small ($ε\ll 1$, but independent of $ν$) Gevrey 2- datum, $ω_0^{(ν)}(x, y)$, that is supported away from the boundaries $y = \pm 1$, we prove the following results: \begin{align*} & \|ω^{(ν)}(t) - \frac{1}{2π}\int ω^{(ν)}(t) dx \|_{L^2} \lesssim εe^{-δν^{1/3} t}, & \text{(Enhanced Dissipation)} \\ & \langle t \rangle \|u_1^{(ν)}(t) - \frac{1}{2π} \int u_1^{(ν)}(t) dx\|_{L^2} + \langle t \rangle^2 \|u_2^{(ν)}(t)\|_{L^2} \lesssim εe^{-δν^{1/3} t}, & \text{(Inviscid Damping)} \\ &\| ω^{(ν)} - ω^{(0)} \|_{L^\infty} \lesssim ενt^{3+η}, \quad\quad t \lesssim ν^{-1/(3+η)} & \text{(Long-time Inviscid Limit)} \end{align*} This is the first nonlinear asymptotic stability result of its type, which combines three important physical phenomena at the nonlinear level: inviscid damping, enhanced dissipation, and long-time inviscid limit in the presence of boundaries. The techniques we develop represent a major departure from prior works on nonlinear inviscid damping as physical space techniques necessarily play a central role. In this paper, we focus on the primary nonlinear result, while tools for handling the linearized parabolic and elliptic equations are developed in our separate, companion work. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 157 pages

Pseudo-Gevrey Smoothing for the Passive Scalar Equations near Couette

Authors: Jacob Bedrossian, Siming He, Sameer Iyer, Fei Wang

Abstract: In this article, we study the regularity theory for two linear equations that are important in fluid dynamics: the passive scalar equation for (time-varying) shear flows close to Couette in $\mathbb T \times [-1,1]$ with vanishing diffusivity $ν\to 0$ and the Poisson equation with right-hand side behaving in similar function spaces to such a passive scalar. The primary motivation for this work is… ▽ More In this article, we study the regularity theory for two linear equations that are important in fluid dynamics: the passive scalar equation for (time-varying) shear flows close to Couette in $\mathbb T \times [-1,1]$ with vanishing diffusivity $ν\to 0$ and the Poisson equation with right-hand side behaving in similar function spaces to such a passive scalar. The primary motivation for this work is to develop some of the main technical tools required for our treatment of the (nonlinear) 2D Navier-Stokes equations, carried out in our companion work. Both equations are studied with homogeneous Dirichlet conditions (the analogue of a Navier slip-type boundary condition) and the initial condition is taken to be compactly supported away from the walls. We develop smoothing estimates with the following three features: [1] Uniform-in-$ν$ regularity is with respect to $\partial_x$ and a time-dependent adapted vector-field $Γ$ which approximately commutes with the passive scalar equation (as opposed to `flat' derivatives), and a scaled gradient $\sqrtν \nabla$; [2] $(\partial_x, Γ)$-regularity estimates are performed in Gevrey spaces with regularity that depends on the spatial coordinate, $y$ (what we refer to as `pseudo-Gevrey'); [3] The regularity of these pseudo-Gevrey spaces degenerates to finite regularity near the center of the channel and hence standard Gevrey product rules and other amenable properties do not hold. Nonlinear analysis in such a delicate functional setting is one of the key ingredients to our companion paper, \cite{BHIW24a}, which proves the full nonlinear asymptotic stability of the Couette flow with slip boundary conditions. The present article introduces new estimates for the associated linear problems in these degenerate pseudo-Gevrey spaces, which is of independent interest. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 130 pages

Network shell structure based on hub and non-hub nodes

Authors: Gaogao Dong, Nannan Sun, Fan Wang, Renaud Lambiotte

Abstract: The shell structure holds significant importance in various domains such as information dissemination, supply chain management, and transportation. This study focuses on investigating the shell structure of hub and non-hub nodes, which play important roles in these domains. Our framework explores the topology of Erdös-Rényi (ER) and Scale-Free (SF) networks, considering source node selection strat… ▽ More The shell structure holds significant importance in various domains such as information dissemination, supply chain management, and transportation. This study focuses on investigating the shell structure of hub and non-hub nodes, which play important roles in these domains. Our framework explores the topology of Erdös-Rényi (ER) and Scale-Free (SF) networks, considering source node selection strategies dependent on the nodes' degrees. We define the shell $l$ in a network as the set of nodes at a distance $l$ from a given node and represent $r_l$ as the fraction of nodes outside shell $l$. Statistical properties of the shells are examined for a selected node, taking into account the node's degree. For a network with a given degree distribution, we analytically derive the degree distribution and average degree of nodes outside shell $l$ as functions of $r_l$. Moreover, we discover that $r_l$ follows an iterative functional form $r_l = φ(r_{l-1})$, where $φ$ is expressed in terms of the generating function of the original degree distribution of the network. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Local well-posedness of strong solutions to the 2D nonhomogeneous primitive equations with density-dependent viscosity

Authors: Quansen Jiu, Lin Ma, Fengchao Wang

Abstract: In this paper, we consider the initial-boundary value problem of the nonhomogeneous primitive equations with density-dependent viscosity. Local well-posedness of strong solutions is established for this system with a natural compatibility condition. The initial density does not need to be strictly positive and may contain vacuum. Meanwhile, we also give the corresponding blow-up criterion if the m… ▽ More In this paper, we consider the initial-boundary value problem of the nonhomogeneous primitive equations with density-dependent viscosity. Local well-posedness of strong solutions is established for this system with a natural compatibility condition. The initial density does not need to be strictly positive and may contain vacuum. Meanwhile, we also give the corresponding blow-up criterion if the maximum existence interval with respect to the time is finite. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Convergence Acceleration of Favre-Averaged Non-Linear Harmonic Method

Authors: Feng Wang, Kurt Webber, David Radford, Luca di Mare, Marcus Meyer

Abstract: This paper develops a numerical procedure to accelerate the convergence of the Favre-averaged Non-Linear Harmonic (FNLH) method. The scheme provides a unified mathematical framework for solving the sparse linear systems formed by the mean flow and the time-linearized harmonic flows of FNLH in an explicit or implicit fashion. The approach explores the similarity of the sparse linear systems of FNLH… ▽ More This paper develops a numerical procedure to accelerate the convergence of the Favre-averaged Non-Linear Harmonic (FNLH) method. The scheme provides a unified mathematical framework for solving the sparse linear systems formed by the mean flow and the time-linearized harmonic flows of FNLH in an explicit or implicit fashion. The approach explores the similarity of the sparse linear systems of FNLH and leads to a memory efficient procedure, so that its memory consumption does not depend on the number of harmonics to compute. The proposed method has been implemented in the industrial CFD solver HYDRA. Two test cases are used to conduct a comparative study of explicit and implicit schemes in terms of convergence, computational efficiency, and memory consumption. Comparisons show that the implicit scheme yields better convergence than the explicit scheme and is also roughly 7 to 10 times more computationally efficient than the explicit scheme with 4 levels of multigrid. Furthermore, the implicit scheme consumes only approximately $50\%$ of the explicit scheme with four levels of multigrid. Compared with the full annulus unsteady Reynolds averaged Navier-Stokes (URANS) simulations, the implicit scheme produces comparable results to URANS with computational time and memory consumption that are two orders of magnitude smaller. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Learning-based Multi-continuum Model for Multiscale Flow Problems

Authors: Fan Wang, Yating Wang, Wing Tat Leung, Zongben Xu

Abstract: Multiscale problems can usually be approximated through numerical homogenization by an equation with some effective parameters that can capture the macroscopic behavior of the original system on the coarse grid to speed up the simulation. However, this approach usually assumes scale separation and that the heterogeneity of the solution can be approximated by the solution average in each coarse blo… ▽ More Multiscale problems can usually be approximated through numerical homogenization by an equation with some effective parameters that can capture the macroscopic behavior of the original system on the coarse grid to speed up the simulation. However, this approach usually assumes scale separation and that the heterogeneity of the solution can be approximated by the solution average in each coarse block. For complex multiscale problems, the computed single effective properties/continuum might be inadequate. In this paper, we propose a novel learning-based multi-continuum model to enrich the homogenized equation and improve the accuracy of the single continuum model for multiscale problems with some given data. Without loss of generalization, we consider a two-continuum case. The first flow equation keeps the information of the original homogenized equation with an additional interaction term. The second continuum is newly introduced, and the effective permeability in the second flow equation is determined by a neural network. The interaction term between the two continua aligns with that used in the Dual-porosity model but with a learnable coefficient determined by another neural network. The new model with neural network terms is then optimized using trusted data. We discuss both direct back-propagation and the adjoint method for the PDE-constraint optimization problem. Our proposed learning-based multi-continuum model can resolve multiple interacted media within each coarse grid block and describe the mass transfer among them, and it has been demonstrated to significantly improve the simulation results through numerical experiments involving both linear and nonlinear flow equations. △ Less

Submitted 20 June, 2024; v1 submitted 20 March, 2024; originally announced March 2024.

Comments: Corrected typos

A Selective Review on Statistical Methods for Massive Data Computation: Distributed Computing, Subsampling, and Minibatch Techniques

Authors: Xuetong Li, Yuan Gao, Hong Chang, Danyang Huang, Yingying Ma, Rui Pan, Haobo Qi, Feifei Wang, Shuyuan Wu, Ke Xu, Jing Zhou, Xuening Zhu, Yingqiu Zhu, Hansheng Wang

Abstract: This paper presents a selective review of statistical computation methods for massive data analysis. A huge amount of statistical methods for massive data computation have been rapidly developed in the past decades. In this work, we focus on three categories of statistical computation methods: (1) distributed computing, (2) subsampling methods, and (3) minibatch gradient techniques. The first clas… ▽ More This paper presents a selective review of statistical computation methods for massive data analysis. A huge amount of statistical methods for massive data computation have been rapidly developed in the past decades. In this work, we focus on three categories of statistical computation methods: (1) distributed computing, (2) subsampling methods, and (3) minibatch gradient techniques. The first class of literature is about distributed computing and focuses on the situation, where the dataset size is too huge to be comfortably handled by one single computer. In this case, a distributed computation system with multiple computers has to be utilized. The second class of literature is about subsampling methods and concerns about the situation, where the sample size of dataset is small enough to be placed on one single computer but too large to be easily processed by its memory as a whole. The last class of literature studies those minibatch gradient related optimization techniques, which have been extensively used for optimizing various deep learning models. △ Less

Submitted 17 March, 2024; originally announced March 2024.

Homotopical Minimal Measures for Geodesic flows on Surfaces of Higher Genus

Authors: Fang Wang, Zhihong Xia

Abstract: We study the homotopical minimal measures for positive definite autonomous Lagrangian systems. Homotopical minimal measures are action-minimizers in their homotopy classes, while the classical minimal measures (Mather measures) are action-minimizers in homology classes. Homotopical minimal measures are much more general, they are not necessarily homological action-minimizers. However, some of them… ▽ More We study the homotopical minimal measures for positive definite autonomous Lagrangian systems. Homotopical minimal measures are action-minimizers in their homotopy classes, while the classical minimal measures (Mather measures) are action-minimizers in homology classes. Homotopical minimal measures are much more general, they are not necessarily homological action-minimizers. However, some of them can be obtained from the classical ones by lifting them to finite-fold covering spaces. We apply this idea of finite covering to the geodesic flows on surfaces of higher genus. Let $(M,G)$ be a compact closed surface with genus $g>1$, where $G$ is a complete Riemannian metric on $M$. Consider the positive definite autonomous Lagrangian $L(x,v)=G_x(v,v)$, whose Lagrangian system $φ_t: TM\rightarrow TM$ is exactly the complete geodesic flow on $TM$. We show that for each homotopical minimal ergodic measure $μ$ that is supported on a nontrivial simple closed periodic trajectory, there is a finite-fold covering space $M'$ such that each ergodic preimage of $μ$ on $TM'$ is a minimal measure in the classic Mather theory for the Lagrangian system on $TM'$. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Global existence and uniqueness of strong solutions to the 2D nonhomogeneous primitive equations with density-dependent viscosity

Authors: Quansen Jiu, Lin Ma, Fengchao Wang

Abstract: This paper is concerned with an initial-boundary value problem of the two-dimensional inhomogeneous primitive equations with density-dependent viscosity. The global well-posedness of strong solutions is established, provided the initial horizontal velocity is suitably small, that is, $\|\nabla u_{0}\|_{L^{2}}\leq η_{0}$ for suitably small $η_{0}>0$. The initial data may contain vacuum. The proof i… ▽ More This paper is concerned with an initial-boundary value problem of the two-dimensional inhomogeneous primitive equations with density-dependent viscosity. The global well-posedness of strong solutions is established, provided the initial horizontal velocity is suitably small, that is, $\|\nabla u_{0}\|_{L^{2}}\leq η_{0}$ for suitably small $η_{0}>0$. The initial data may contain vacuum. The proof is based on the local well-posedness and the blow-up criterion proved in \cite{0}, which states that if $T^{*}$ is the maximal existence time of the local strong solutions $(ρ,u,w,P)$ and $T^{*}<\infty$, then \begin{equation*} \sup_{0\leq t<T^{*}}(\left\|\nabla ρ(t)\right\|_{L^{\infty}}+\left\|\nabla^{2}ρ(t)\right\|_{L^{2}}+\left\|\nabla u(t)\right\|_{L^{2}})=\infty. \end{equation*} To complete the proof, it is required to make an estimate on a key term $\|\nabla u_{t}\|_{L_{t}^{1}L_Ω^{2}}$. We prove that it is bounded and could be as small as desired under certain smallness conditions, by making use of the regularity result of hydrostatic Stokes equations and some careful time weighted estimates. △ Less

Submitted 1 March, 2024; originally announced March 2024.

An inexact Bregman proximal point method and its acceleration version for unbalanced optimal transport

Authors: Xiang Chen, Faqiang Wang, Jun Liu, Li Cui

Abstract: The Unbalanced Optimal Transport (UOT) problem plays increasingly important roles in computational biology, computational imaging and deep learning. Scaling algorithm is widely used to solve UOT due to its convenience and good convergence properties. However, this algorithm has lower accuracy for large regularization parameters, and due to stability issues, small regularization parameters can easi… ▽ More The Unbalanced Optimal Transport (UOT) problem plays increasingly important roles in computational biology, computational imaging and deep learning. Scaling algorithm is widely used to solve UOT due to its convenience and good convergence properties. However, this algorithm has lower accuracy for large regularization parameters, and due to stability issues, small regularization parameters can easily lead to numerical overflow. We address this challenge by developing an inexact Bregman proximal point method for solving UOT. This algorithm approximates the proximal operator using the Scaling algorithm at each iteration. The algorithm (1) converges to the true solution of UOT, (2) has theoretical guarantees and robust regularization parameter selection, (3) mitigates numerical stability issues, and (4) can achieve comparable computational complexity to the Scaling algorithm in specific practice. Building upon this, we develop an accelerated version of inexact Bregman proximal point method for solving UOT by using acceleration techniques of Bregman proximal point method and provide theoretical guarantees and experimental validation of convergence and acceleration. △ Less

Submitted 26 February, 2024; originally announced February 2024.

Diffusion Processes on $p$-Wasserstein Space over Banach Space

Authors: Panpan Ren, Feng-Yu Wang, Simon Wittmann

Abstract: To study diffusion processes on the $p$-Wasserstein space $\scr P_p $ for $p\in [1,\infty)$ over a separable, reflexive Banach space $X$, we present a criterion on the quasi-regularity of Dirichlet forms in $L^2(\scr P_p,Λ)$ for a reference probability $Λ$ on $\scr P_p$. It is formulated in terms of an upper bound condition with the uniform norm of the intrinsic derivative. The condition is easy t… ▽ More To study diffusion processes on the $p$-Wasserstein space $\scr P_p $ for $p\in [1,\infty)$ over a separable, reflexive Banach space $X$, we present a criterion on the quasi-regularity of Dirichlet forms in $L^2(\scr P_p,Λ)$ for a reference probability $Λ$ on $\scr P_p$. It is formulated in terms of an upper bound condition with the uniform norm of the intrinsic derivative. The condition is easy to check in relevant applications and allows to construct a type of Ornstein-Uhlenbeck process on $\scr P_p$. We find a versatile class of quasi-regular local Dirichlet forms on $\scr P_p$ by using images of Dirichlet forms on the tangent space $L^p(X\to X,μ_0)$ at a reference point $μ_0\in \scr P_p$. The Ornstein-Uhlenbeck type Dirichlet form is an important example in this class. An $L^2$-estimate for the corresponding heat kernel is derived, based on the eigenvalues of the covariance operator of the underlying Gaussian measure. △ Less

Submitted 23 April, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: Acknowledgments added; Application of Sect. 3.2 revised; Typos corrected

MSC Class: 60J60; 60J25; 60J46; 47D07; 60G57; 60J45

A parallel domain decomposition method for solving elliptic equations on manifolds

Authors: Lizhen Qin, Feng Wang, Yun Wang

Abstract: We propose a new numerical domain decomposition method for solving elliptic equations on compact Riemannian manifolds. One advantage of this method is its ability to bypass the need for global triangulations or grids on the manifolds. Additionally, it features a highly parallel iterative scheme. To verify its efficacy, we conduct numerical experiments on some $4$-dimensional manifolds without and… ▽ More We propose a new numerical domain decomposition method for solving elliptic equations on compact Riemannian manifolds. One advantage of this method is its ability to bypass the need for global triangulations or grids on the manifolds. Additionally, it features a highly parallel iterative scheme. To verify its efficacy, we conduct numerical experiments on some $4$-dimensional manifolds without and with boundary. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: Comments are welcome!

MSC Class: Primary 65N30; Secondary 58J05; 65N55

The rigidity of eigenvalues on Kähler manifolds with positive Ricci lower bound

Authors: Jianchun Chu, Feng Wang, Kewei Zhang

Abstract: In this work, optimal rigidity results for eigenvalues on Kähler manifolds with positive Ricci lower bound are established. More precisely, for those Kähler manifolds whose first eigenvalue agrees with the Ricci lower bound, we show that the complex projective space is the only one with the largest multiplicity of the first eigenvalue. Moreover, there is a specific gap between the largest and the… ▽ More In this work, optimal rigidity results for eigenvalues on Kähler manifolds with positive Ricci lower bound are established. More precisely, for those Kähler manifolds whose first eigenvalue agrees with the Ricci lower bound, we show that the complex projective space is the only one with the largest multiplicity of the first eigenvalue. Moreover, there is a specific gap between the largest and the second largest multiplicity. In the Kähler--Einstein case, almost rigidity results for eigenvalues are also obtained. △ Less

Submitted 28 January, 2024; originally announced January 2024.

Comments: Comments are welcome!

Obtaining the pseudoinverse solution of singular range-symmetric linear systems with GMRES-type methods

Authors: Kui Du, Jia-Jun Fan, Fang Wang

Abstract: It is well known that for singular inconsistent range-symmetric linear systems, the generalized minimal residual (GMRES) method determines a least squares solution without breakdown. The reached least squares solution may be or not be the pseudoinverse solution. We show that a lift strategy can be used to obtain the pseudoinverse solution. In addition, we propose a new iterative method named RSMAR… ▽ More It is well known that for singular inconsistent range-symmetric linear systems, the generalized minimal residual (GMRES) method determines a least squares solution without breakdown. The reached least squares solution may be or not be the pseudoinverse solution. We show that a lift strategy can be used to obtain the pseudoinverse solution. In addition, we propose a new iterative method named RSMAR (minimum $\mathbf A$-residual) for range-symmetric linear systems $\mathbf A\mathbf x=\mathbf b$. At step $k$ RSMAR minimizes $\|\mathbf A\mathbf r_k\|$ in the $k$th Krylov subspace generated with $\{\mathbf A, \mathbf r_0\}$ rather than $\|\mathbf r_k\|$, where $\mathbf r_k$ is the $k$th residual vector and $\|\cdot\|$ denotes the Euclidean vector norm. We show that RSMAR and GMRES terminate with the same least squares solution when applied to range-symmetric linear systems. We provide two implementations for RSMAR. Our numerical experiments show that RSMAR is the most suitable method among GMRES-type methods for singular inconsistent range-symmetric linear systems. △ Less

Submitted 22 January, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: 22 pages, 4 figures

MSC Class: 15A06; 15A09; 65F10; 65F25; 65F50

Energy-Efficient Data Offloading for Earth Observation Satellite Networks

Authors: Lijun He, Ziye Jia, Juncheng Wang, Feng Wang, Erick Lansard, Chau Yuen

Abstract: In Earth Observation Satellite Networks (EOSNs) with a large number of battery-carrying satellites, proper power allocation and task scheduling are crucial to improving the data offloading efficiency. As such, we jointly optimize power allocation and task scheduling to achieve energy-efficient data offloading in EOSNs, aiming to balance the objectives of reducing the total energy consumption and i… ▽ More In Earth Observation Satellite Networks (EOSNs) with a large number of battery-carrying satellites, proper power allocation and task scheduling are crucial to improving the data offloading efficiency. As such, we jointly optimize power allocation and task scheduling to achieve energy-efficient data offloading in EOSNs, aiming to balance the objectives of reducing the total energy consumption and increasing the sum weights of tasks. First, we derive the optimal power allocation solution to the joint optimization problem when the task scheduling policy is given. Second, leveraging the conflict graph model, we transform the original joint optimization problem into a maximum weight independent set problem when the power allocation strategy is given. Finally, we utilize the genetic framework to combine the above special solutions as a two-layer solution for the joint optimization problem. Simulation results demonstrate that our proposed solution can properly balance the sum weights of tasks and the total energy consumption, achieving superior system performance over the current best alternatives. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Spectral integrated neural networks (SINNs) for solving forward and inverse dynamic problems

Authors: Lin Qiu, Fajie Wang, Wenzhen Qu, Yan Gu, Qing-Hua Qin

Abstract: This paper proposes a novel neural network framework, denoted as spectral integrated neural networks (SINNs), for resolving three-dimensional forward and inverse dynamic problems. In the SINNs, the spectral integration method is applied to perform temporal discretization, and then a fully connected neural network is adopted to solve resulting partial differential equations (PDEs) in the spatial do… ▽ More This paper proposes a novel neural network framework, denoted as spectral integrated neural networks (SINNs), for resolving three-dimensional forward and inverse dynamic problems. In the SINNs, the spectral integration method is applied to perform temporal discretization, and then a fully connected neural network is adopted to solve resulting partial differential equations (PDEs) in the spatial domain. Specifically, spatial coordinates are employed as inputs in the network architecture, and the output layer is configured with multiple outputs, each dedicated to approximating solutions at different time instances characterized by Gaussian points used in the spectral method. By leveraging the automatic differentiation technique and spectral integration scheme, the SINNs minimize the loss function, constructed based on the governing PDEs and boundary conditions, to obtain solutions for dynamic problems. Additionally, we utilize polynomial basis functions to expand the unknown function, aiming to enhance the performance of SINNs in addressing inverse problems. The conceived framework is tested on six forward and inverse dynamic problems, involving nonlinear PDEs. Numerical results demonstrate the superior performance of SINNs over the popularly used physics-informed neural networks in terms of convergence speed, computational accuracy and efficiency. It is also noteworthy that the SINNs exhibit the capability to deliver accurate and stable solutions for long-time dynamic problems. △ Less

Submitted 7 January, 2024; originally announced January 2024.

GPBiLQ and GPQMR: Two iterative methods for unsymmetric partitioned linear systems

Authors: Kui Du, Jia-Jun Fan, Fang Wang

Abstract: We introduce two iterative methods, GPBiLQ and GPQMR, for solving unsymmetric partitioned linear systems. The basic mechanism underlying GPBiLQ and GPQMR is a novel simultaneous tridiagonalization via biorthogonality that allows for short-recurrence iterative schemes. Similar to the biconjugate gradient method, it is possible to develop another method, GPBiCG, whose iterate (if it exists) can be o… ▽ More We introduce two iterative methods, GPBiLQ and GPQMR, for solving unsymmetric partitioned linear systems. The basic mechanism underlying GPBiLQ and GPQMR is a novel simultaneous tridiagonalization via biorthogonality that allows for short-recurrence iterative schemes. Similar to the biconjugate gradient method, it is possible to develop another method, GPBiCG, whose iterate (if it exists) can be obtained inexpensively from the GPBiLQ iterate. Whereas the iterate of GPBiCG may not exist, the iterates of GPBiLQ and GPQMR are always well defined as long as the biorthogonal tridiagonal reduction process does not break down. We discuss connections between the proposed methods and some existing methods, and give numerical experiments to illustrate the performance of the proposed methods. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 22 pages, 4 figures

MSC Class: 15A06; 65F10; 65F25; 65F50

Functional Inequalities for Brownian Motion on Riemannian Manifolds with Sticky-Reflecting Boundary Diffusion

Authors: Marie Bormann, Max von Renesse, Feng-Yu Wang

Abstract: We prove geometric upper bounds for the Poincaré and Logarithmic Sobolev constants for Brownian motion on manifolds with sticky reflecting boundary diffusion i.e. extended Wentzell-type boundary condition under general curvature assumptions on the manifold and its boundary. The method is based on an interpolation involving energy interactions between the boundary and the interior of the manifold.… ▽ More We prove geometric upper bounds for the Poincaré and Logarithmic Sobolev constants for Brownian motion on manifolds with sticky reflecting boundary diffusion i.e. extended Wentzell-type boundary condition under general curvature assumptions on the manifold and its boundary. The method is based on an interpolation involving energy interactions between the boundary and the interior of the manifold. As side results we obtain explicit geometric bounds on the first nontrivial Steklov eigenvalue, for the norm of the boundary trace operator on Sobolev functions, and on the boundary trace logarithmic Sobolev constant. The case of Brownian motion with pure sticky reflection is also treated. △ Less

Submitted 3 April, 2024; v1 submitted 30 December, 2023; originally announced January 2024.

MSC Class: 60J65; 35A23; 58C40

Local Randomized Neural Networks with Hybridized Discontinuous Petrov-Galerkin Methods for Stokes-Darcy Flows

Authors: Haoning Dang, Fei Wang

Abstract: This paper introduces a new numerical approach that integrates local randomized neural networks (LRNNs) and the hybridized discontinuous Petrov-Galerkin (HDPG) method for solving coupled fluid flow problems. The proposed method partitions the domain of interest into several subdomains and constructs an LRNN on each subdomain. Then, the HDPG scheme is used to couple the LRNNs to approximate the unk… ▽ More This paper introduces a new numerical approach that integrates local randomized neural networks (LRNNs) and the hybridized discontinuous Petrov-Galerkin (HDPG) method for solving coupled fluid flow problems. The proposed method partitions the domain of interest into several subdomains and constructs an LRNN on each subdomain. Then, the HDPG scheme is used to couple the LRNNs to approximate the unknown functions. We develop LRNN-HDPG methods based on velocity-stress formulation to solve two types of problems: Stokes-Darcy problems and Brinkman equations, which model the flow in porous media and free flow. We devise a simple and effective way to deal with the interface conditions in the Stokes-Darcy problems without adding extra terms to the numerical scheme. We conduct extensive numerical experiments to demonstrate the stability, efficiency, and robustness of the proposed method. The numerical results show that the LRNN-HDPG method can achieve high accuracy with a small number of degrees of freedom. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: 21 pages, 4 figures

MSC Class: 65N30; 41A46

Stability threshold of nearly-Couette shear flows with Navier boundary conditions in 2D

Authors: Jacob Bedrossian, Siming He, Sameer Iyer, Fei Wang

Abstract: In this work, we prove a threshold theorem for the 2D Navier-Stokes equations posed on the periodic channel, $\mathbb{T} \times [-1,1]$, supplemented with Navier boundary conditions $ω|_{y = \pm 1} = 0$. Initial datum is taken to be a perturbation of Couette in the following sense: the shear component of the perturbation is assumed small (in an appropriate Sobolev space) but importantly is indepen… ▽ More In this work, we prove a threshold theorem for the 2D Navier-Stokes equations posed on the periodic channel, $\mathbb{T} \times [-1,1]$, supplemented with Navier boundary conditions $ω|_{y = \pm 1} = 0$. Initial datum is taken to be a perturbation of Couette in the following sense: the shear component of the perturbation is assumed small (in an appropriate Sobolev space) but importantly is independent of $ν$. On the other hand, the nonzero modes are assumed size $O(ν^{\frac12})$ in an anisotropic Sobolev space. For such datum, we prove nonlinear enhanced dissipation and inviscid damping for the resulting solution. The principal innovation is to capture quantitatively the \textit{inviscid damping}, for which we introduce a new Singular Integral Operator which is a physical space analogue of the usual Fourier multipliers which are used to prove damping. We then include this SIO in the context of a nonlinear hypocoercivity framework. △ Less

Submitted 31 October, 2023; originally announced November 2023.

An Unconditionally Stable Iterative Decoupled Algorithm for Multiple-Network Poroelasticity

Authors: Meng Lei, Mingchao Cai, Feng Wang

Abstract: In this work, we introduce an iterative decoupled algorithm designed for addressing the quasi-static multiple-network poroelasticity problem. This problem pertains to the simultaneous modeling of fluid flow and deformations within an elastic porous medium permeated by multiple fluid networks, each with distinct characteristics. Our approach focuses on the total-pressure-based formulation, which tr… ▽ More In this work, we introduce an iterative decoupled algorithm designed for addressing the quasi-static multiple-network poroelasticity problem. This problem pertains to the simultaneous modeling of fluid flow and deformations within an elastic porous medium permeated by multiple fluid networks, each with distinct characteristics. Our approach focuses on the total-pressure-based formulation, which treats the solid displacement, total pressure, and network pressures as primary unknowns. This formulation transforms the original problem into a combination of the generalized Stokes problem and the parabolic problem, offering certain advantages such as mitigating elastic locking effects and streamlining the discretization process. Notably, the algorithm ensures unconditional convergence to the solution of the total-pressure-based coupled algorithm. To validate the accuracy and efficiency of our method, we present numerical experiments. The robustness of the algorithm with respect to the physical parameters and the discretization parameters is carefully investigated. △ Less

Submitted 27 October, 2023; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: to be submitted

Transfer learning for piecewise-constant mean estimation: Optimality, $\ell_1$- and $\ell_0$-penalisation

Authors: Fan Wang, Yi Yu

Abstract: We study transfer learning for estimating piecewise-constant signals when source data, which may be relevant but disparate, are available in addition to the target data. We first investigate transfer learning estimators that respectively employ $\ell_1$- and $\ell_0$-penalties for unisource data scenarios and then generalise these estimators to accommodate multisources. To further reduce estimatio… ▽ More We study transfer learning for estimating piecewise-constant signals when source data, which may be relevant but disparate, are available in addition to the target data. We first investigate transfer learning estimators that respectively employ $\ell_1$- and $\ell_0$-penalties for unisource data scenarios and then generalise these estimators to accommodate multisources. To further reduce estimation errors, especially when some sources significantly differ from the target, we introduce an informative source selection algorithm. We then examine these estimators with multisource selection and establish their minimax optimality. Unlike the common narrative in the transfer learning literature that the performance is enhanced through large source sample sizes, our approaches leverage higher observation frequencies and accommodate diverse frequencies across multiple sources. Our theoretical findings are supported by extensive numerical experiments, with the code available online, see https://github.com/chrisfanwang/transferlearning △ Less

Submitted 29 October, 2023; v1 submitted 9 October, 2023; originally announced October 2023.

Adaptive finite element approximation of bilinear optimal control with fractional Laplacian

Authors: Fangyuan Wang, Qiming Wang, Zhaojie Zhou

Abstract: We investigate the application of a posteriori error estimates to a fractional optimal control problem with pointwise control constraints. Specifically, we address a problem in which the state equation is formulated as an integral form of the fractional Laplacian equation, with the control variable embedded within the state equation as a coefficient. We propose two distinct finite element discreti… ▽ More We investigate the application of a posteriori error estimates to a fractional optimal control problem with pointwise control constraints. Specifically, we address a problem in which the state equation is formulated as an integral form of the fractional Laplacian equation, with the control variable embedded within the state equation as a coefficient. We propose two distinct finite element discretization approaches for an optimal control problem. The first approach employs a fully discrete scheme where the control variable is discretized using piecewise constant functions. The second approach, a semi-discrete scheme, does not discretize the control variable. Using the first-order optimality condition, the second-order optimality condition, and a solution regularity analysis for the optimal control problem, we devise a posteriori error estimates. We subsequently demonstrate the reliability and efficiency of the proposed error estimators. Based on the established error estimates framework, an adaptive refinement strategy is developed to help achieve the optimal convergence rate. The effectiveness of the refinement strategy is verified by numerical experiments. △ Less

Submitted 8 October, 2023; originally announced October 2023.

Unadjusted Langevin Algorithms for SDEs with Hoelder Drift

Authors: Xiang Li, Feng-Yu Wang, Lihu Xu

Abstract: Consider the following stochastic differential equation for $(X_t)_{t\ge 0}$ on $\mathbb R^d$ and its Euler-Maruyama (EM) approximation $(Y_{t_n})_{n\in \mathbb Z^+}$: \begin{align*} &d X_t=b( X_t) d t+σ(X_t) d B_t, \\ & Y_{t_{n+1}}=Y_{t_{n}}+η_{n+1} b(Y_{t_{n}})+σ(Y_{t_{n}})\left(B_{t_{n+1}}-B_{t_{n}}\right), \end{align*} where… ▽ More Consider the following stochastic differential equation for $(X_t)_{t\ge 0}$ on $\mathbb R^d$ and its Euler-Maruyama (EM) approximation $(Y_{t_n})_{n\in \mathbb Z^+}$: \begin{align*} &d X_t=b( X_t) d t+σ(X_t) d B_t, \\ & Y_{t_{n+1}}=Y_{t_{n}}+η_{n+1} b(Y_{t_{n}})+σ(Y_{t_{n}})\left(B_{t_{n+1}}-B_{t_{n}}\right), \end{align*} where $b:\mathbb{R}^d \rightarrow \mathbb{R}^d,\ \ σ: \mathbb R^d \rightarrow \mathbb{R}^{d \times d}$ are measurable, $B_t$ is the $d$-dimensional Brownian motion, $t_0:=0,t_{n}:=\sum_{k=1}^{n} η_{k}$ for constants $η_k>0$ satisfying $\lim_{k \rightarrow \infty} η_k=0$ and $\sum_{k=1}^\inftyη_k =\infty$. Under (partial) dissipation conditions ensuring the ergodicity, we obtain explicit convergence rates of $\mathbb W_p(\mathscr{L}(Y_{t_n}), \mathscr{L}(X_{t_n}))+\mathbb W_p(\mathscr{L}(Y_{t_n}), μ)\rightarrow 0$ as $n\rightarrow \infty$, where $\mathbb W_p$ is the $L^p$-Wasserstein distance for certain $p\in [0,\infty)$, $\mathscr{L}(ξ)$ is the distribution of random variable $ξ$, and $μ$ is the unique invariant probability measure of $(X_t)_{t \ge 0}$. Comparing with the existing results where $b$ is at least $C^2$-smooth, our estimates apply to Hoelder continuous drift and can be sharp in several specific situations. △ Less

Submitted 29 September, 2023; originally announced October 2023.

Adaptive finite element approximation of sparse optimal control with integral fractional Laplacian

Authors: Fangyuan Wang, Qiming Wang, Zhaojie Zhou

Abstract: In this paper we present and analyze a weighted residual a posteriori error estimate for an optimal control problem. The problem involves a nondifferentiable cost functional, a state equation with an integral fractional Laplacian, and control constraints. We employ subdifferentiation in the context of nondifferentiable convex analysis to obtain first-order optimality conditions. Piecewise linear p… ▽ More In this paper we present and analyze a weighted residual a posteriori error estimate for an optimal control problem. The problem involves a nondifferentiable cost functional, a state equation with an integral fractional Laplacian, and control constraints. We employ subdifferentiation in the context of nondifferentiable convex analysis to obtain first-order optimality conditions. Piecewise linear polynomials are utilized to approximate the solutions of the state and adjoint equations. The control variable is discretized using the variational discretization method. Upper and lower bounds for the a posteriori error estimate of the finite element approximation of the optimal control problem are derived. In the region where 3/2 < alpha < 2, the residuals do not satisfy the L2(Omega) regularity. To address this issue, an additional weight is included in the weighted residual estimator, which is based on a power of the distance from the mesh skeleton. Furthermore, we propose an h-adaptive algorithm driven by the posterior view error estimator, utilizing the Dorfler labeling criterion. The convergence analysis results show that the approximation sequence generated by the adaptive algorithm converges at the optimal algebraic rate. Finally, numerical experiments are conducted to validate the theoretical results. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Wasserstein Convergence Rate for Empirical Measures of Markov Processes

Authors: Feng-Yu Wang

Abstract: The convergence rate in Wasserstein distance is estimated for empirical measures of ergodic Markov processes, and the estimate can be sharp in some specific situations. The main result is applied to subordinations of typical models excluded by existing results, which include: stochastic Hamiltonian systems on $\mathbb R^{n}\times \mathbb R^{m}$, spherical velocity Langevin processes on… ▽ More The convergence rate in Wasserstein distance is estimated for empirical measures of ergodic Markov processes, and the estimate can be sharp in some specific situations. The main result is applied to subordinations of typical models excluded by existing results, which include: stochastic Hamiltonian systems on $\mathbb R^{n}\times \mathbb R^{m}$, spherical velocity Langevin processes on $\mathbb R^n\times\mathbb S^{n-1},$ multi-dimensional Wright-Fisher type diffusion processes, and stable type jump processes. △ Less

Submitted 2 October, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 35 pages

Entropy Estimate for Degenerate SDEs with Applications to Nonlinear Kinetic Fokker-Planck Equations

Authors: Zhongmin Qian, Panpan Ren, Feng-Yu Wang

Abstract: The relative entropy for two different degenerate diffusion processes is estimated by using the Wasserstein distance of initial distributions and the difference between coefficients. As applications, the entropy cost inequality and exponential ergodicity in entropy are derived for distribution dependent stochastic Hamiltonian systems associated with nonlinear kinetic Fokker Planck equations. The relative entropy for two different degenerate diffusion processes is estimated by using the Wasserstein distance of initial distributions and the difference between coefficients. As applications, the entropy cost inequality and exponential ergodicity in entropy are derived for distribution dependent stochastic Hamiltonian systems associated with nonlinear kinetic Fokker Planck equations. △ Less

Submitted 30 April, 2024; v1 submitted 30 August, 2023; originally announced August 2023.

Comments: 22 pages

MSC Class: 60J60 60H30

Towards Integrated Traffic Control with Operating Decentralized Autonomous Organization

Authors: Shengyue Yao, Jingru Yu, Yi Yu, Jia Xu, Xingyuan Dai, Honghai Li, Fei-Yue Wang, Yilun Lin

Abstract: With a growing complexity of the intelligent traffic system (ITS), an integrated control of ITS that is capable of considering plentiful heterogeneous intelligent agents is desired. However, existing control methods based on the centralized or the decentralized scheme have not presented their competencies in considering the optimality and the scalability simultaneously. To address this issue, we p… ▽ More With a growing complexity of the intelligent traffic system (ITS), an integrated control of ITS that is capable of considering plentiful heterogeneous intelligent agents is desired. However, existing control methods based on the centralized or the decentralized scheme have not presented their competencies in considering the optimality and the scalability simultaneously. To address this issue, we propose an integrated control method based on the framework of Decentralized Autonomous Organization (DAO). The proposed method achieves a global consensus on energy consumption efficiency (ECE), meanwhile to optimize the local objectives of all involved intelligent agents, through a consensus and incentive mechanism. Furthermore, an operation algorithm is proposed regarding the issue of structural rigidity in DAO. Specifically, the proposed operation approach identifies critical agents to execute the smart contract in DAO, which ultimately extends the capability of DAO-based control. In addition, a numerical experiment is designed to examine the performance of the proposed method. The experiment results indicate that the controlled agents can achieve a consensus faster on the global objective with improved local objectives by the proposed method, compare to existing decentralized control methods. In general, the proposed method shows a great potential in developing an integrated control system in the ITS △ Less

Submitted 25 July, 2023; originally announced August 2023.

Comments: 6 pages, 6 figures. To be published in 2023 IEEE 26th International Conference on Intelligent Transportation Systems (ITSC)

Randomized Neural Networks with Petrov-Galerkin Methods for Solving Linear Elasticity Problems

Authors: Yong Shang, Fei Wang

Abstract: We develop the Randomized Neural Networks with Petrov-Galerkin Methods (RNN-PG methods) to solve linear elasticity problems. RNN-PG methods use Petrov-Galerkin variational framework, where the solution is approximated by randomized neural networks and the test functions are piecewise polynomials. Unlike conventional neural networks, the parameters of the hidden layers of the randomized neural netw… ▽ More We develop the Randomized Neural Networks with Petrov-Galerkin Methods (RNN-PG methods) to solve linear elasticity problems. RNN-PG methods use Petrov-Galerkin variational framework, where the solution is approximated by randomized neural networks and the test functions are piecewise polynomials. Unlike conventional neural networks, the parameters of the hidden layers of the randomized neural networks are fixed randomly, while the parameters of the output layer are determined by the least square method, which can effectively approximate the solution. We also develop mixed RNN-PG methods for linear elasticity problems, which ensure the symmetry of the stress tensor and avoid locking effects. We compare RNN-PG methods with the finite element method, the mixed discontinuous Galerkin method, and the physics-informed neural network on several examples, and the numerical results demonstrate that RNN-PG methods achieve higher accuracy and efficiency. △ Less

Submitted 6 August, 2023; originally announced August 2023.

Comments: 15 pages, 4 figures

MSC Class: 65N30; 41A46

Local Randomized Neural Networks Methods for Interface Problems

Authors: Yunlong Li, Fei Wang

Abstract: Accurate modeling of complex physical problems, such as fluid-structure interaction, requires multiphysics coupling across the interface, which often has intricate geometry and dynamic boundaries. Conventional numerical methods face challenges in handling interface conditions. Deep neural networks offer a mesh-free and flexible alternative, but they suffer from drawbacks such as time-consuming opt… ▽ More Accurate modeling of complex physical problems, such as fluid-structure interaction, requires multiphysics coupling across the interface, which often has intricate geometry and dynamic boundaries. Conventional numerical methods face challenges in handling interface conditions. Deep neural networks offer a mesh-free and flexible alternative, but they suffer from drawbacks such as time-consuming optimization and local optima. In this paper, we propose a mesh-free approach based on Randomized Neural Networks (RNNs), which avoid optimization solvers during training, making them more efficient than traditional deep neural networks. Our approach, called Local Randomized Neural Networks (LRNNs), uses different RNNs to approximate solutions in different subdomains. We discretize the interface problem into a linear system at randomly sampled points across the domain, boundary, and interface using a finite difference scheme, and then solve it by a least-square method. For time-dependent interface problems, we use a space-time approach based on LRNNs. We show the effectiveness and robustness of the LRNNs methods through numerical examples of elliptic and parabolic interface problems. We also demonstrate that our approach can handle high-dimension interface problems. Compared to conventional numerical methods, our approach achieves higher accuracy with fewer degrees of freedom, eliminates the need for complex interface meshing and fitting, and significantly reduces training time, outperforming deep neural networks. △ Less

Submitted 6 August, 2023; originally announced August 2023.

Comments: 22 pages, 15 figures

MSC Class: 65M06; 68T07; 41A46

On Krylov subspace methods for skew-symmetric and shifted skew-symmetric linear systems

Authors: Kui Du, Jia-Jun Fan, Xiao-Hui Sun, Fang Wang, Ya-Lan Zhang

Abstract: Krylov subspace methods for solving linear systems of equations involving skew-symmetric matrices have gained recent attention. Numerical equivalences among Krylov subspace methods for nonsingular skew-symmetric linear systems have been given in Greif et al. [SIAM J. Matrix Anal. Appl., 37 (2016), pp. 1071--1087]. In this work, we extend the results of Greif et al. to singular skew-symmetric linea… ▽ More Krylov subspace methods for solving linear systems of equations involving skew-symmetric matrices have gained recent attention. Numerical equivalences among Krylov subspace methods for nonsingular skew-symmetric linear systems have been given in Greif et al. [SIAM J. Matrix Anal. Appl., 37 (2016), pp. 1071--1087]. In this work, we extend the results of Greif et al. to singular skew-symmetric linear systems. In addition, we systematically study three Krylov subspace methods (called S$^3$CG, S$^3$MR, and S$^3$LQ) for solving shifted skew-symmetric linear systems. They all are based on Lanczos triangularization for skew-symmetric matrices, and correspond to CG, MINRES, and SYMMLQ for solving symmetric linear systems, respectively. To the best of our knowledge, this is the first work that studies S$^3$LQ. We give some new theoretical results on S$^3$CG, S$^3$MR, and S$^3$LQ. We also provide the relationship among the three methods and those based on Golub--Kahan bidiagonalization and Saunders--Simon--Yip tridiagonalization. Numerical examples are given to illustrate our theoretical findings. △ Less

Submitted 31 July, 2023; originally announced July 2023.

Comments: 23 pages, 3 figures

Boundary integrated neural networks (BINNs) for acoustic radiation and scattering

Authors: Wenzhen Qu, Yan Gu, Shengdong Zhao, Fajie wang

Abstract: This paper presents a novel approach called the boundary integrated neural networks (BINNs) for analyzing acoustic radiation and scattering. The method introduces fundamental solutions of the time-harmonic wave equation to encode the boundary integral equations (BIEs) within the neural networks, replacing the conventional use of the governing equation in physics-informed neural networks (PINNs). T… ▽ More This paper presents a novel approach called the boundary integrated neural networks (BINNs) for analyzing acoustic radiation and scattering. The method introduces fundamental solutions of the time-harmonic wave equation to encode the boundary integral equations (BIEs) within the neural networks, replacing the conventional use of the governing equation in physics-informed neural networks (PINNs). This approach offers several advantages. Firstly, the input data for the neural networks in the BINNs only require the coordinates of "boundary" collocation points, making it highly suitable for analyzing acoustic fields in unbounded domains. Secondly, the loss function of the BINNs is not a composite form, and has a fast convergence. Thirdly, the BINNs achieve comparable precision to the PINNs using fewer collocation points and hidden layers/neurons. Finally, the semi-analytic characteristic of the BIEs contributes to the higher precision of the BINNs. Numerical examples are presented to demonstrate the performance of the proposed method. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Surfaces in which every point sounds the same

Authors: Feng Wang, Emmett L. Wyman, Yakun Xi

Abstract: We address a maximally structured case of the question, "Can you hear your location on a manifold," posed in arXiv:2304.04659 for dimension $2$. In short, we show that if a compact surface without boundary sounds the same at every point, then the surface has a transitive action by the isometry group. In the process, we show that you can hear your location on Klein bottles and that you can hear the… ▽ More We address a maximally structured case of the question, "Can you hear your location on a manifold," posed in arXiv:2304.04659 for dimension $2$. In short, we show that if a compact surface without boundary sounds the same at every point, then the surface has a transitive action by the isometry group. In the process, we show that you can hear your location on Klein bottles and that you can hear the lengths and multiplicities of looping geodesics on compact hyperbolic quotients. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 9 pages, 1 figure

Ambarzumyan-type theorem for the Sturm-Liouville operator on the lasso graph

Authors: Feng Wang, Chuan-Fu Yang

Abstract: We consider the Sturm-Liouville operator on the lasso graph with a segment and a loop joined at one point, which has arbitrary length. The Ambarzumyan's theorem for the operator is proved, which says that if the eigenvalues of the operator coincide with those of the zero potential, then the potential is zero. We consider the Sturm-Liouville operator on the lasso graph with a segment and a loop joined at one point, which has arbitrary length. The Ambarzumyan's theorem for the operator is proved, which says that if the eigenvalues of the operator coincide with those of the zero potential, then the potential is zero. △ Less

Submitted 31 May, 2023; originally announced June 2023.

Local Randomized Neural Networks with Discontinuous Galerkin Methods for Diffusive-Viscous Wave Equation

Authors: Jingbo Sun, Fei Wang

Abstract: The diffusive-viscous wave equation is an advancement in wave equation theory, as it accounts for both diffusion and viscosity effects. This has a wide range of applications in geophysics, such as the attenuation of seismic waves in fluid-saturated solids and frequency-dependent phenomena in porous media. Therefore, the development of an efficient numerical method for the equation is of both theor… ▽ More The diffusive-viscous wave equation is an advancement in wave equation theory, as it accounts for both diffusion and viscosity effects. This has a wide range of applications in geophysics, such as the attenuation of seismic waves in fluid-saturated solids and frequency-dependent phenomena in porous media. Therefore, the development of an efficient numerical method for the equation is of both theoretical and practical importance. Recently, local randomized neural networks with discontinuous Galerkin (LRNN-DG) methods have been introduced in \cite{Sun2022lrnndg} to solve elliptic and parabolic equations. Numerical examples suggest that LRNN-DG can achieve high accuracy, and can handle time-dependent problems naturally and efficiently by using a space-time framework. In this paper, we develop LRNN-DG methods for solving the diffusive-viscous wave equation and present numerical experiments with several cases. The numerical results show that the proposed methods can solve the diffusive-viscous wave equation more accurately with less computing costs than traditional methods. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: 20 pages, 5 figures

MSC Class: 65M60; 68T07

Ambarzumyan-type theorem for vectorial Sturm-Liouville operator with impulses

Authors: Feng Wang, Chuan-Fu Yang

Abstract: We consider the vector-impulsive Sturm-Liouville problem with Neumann conditions. The Ambarzumyan$^{\textbf{,}}$s theorem for the problem is proved, which states that if the eigenvalues of the problem coincide with those of the zero potential, then the potential is zero. We consider the vector-impulsive Sturm-Liouville problem with Neumann conditions. The Ambarzumyan$^{\textbf{,}}$s theorem for the problem is proved, which states that if the eigenvalues of the problem coincide with those of the zero potential, then the potential is zero. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Incomplete inverse problem for Dirac operator with constant delay

Authors: Feng Wang, Chuan-Fu Yang

Abstract: In this work, we consider Dirac-type operators with a constant delay less than two-fifths of the interval and not less than one-third of the interval. For our considered Dirac-type operators, an incomplete inverse spectral problem is studied. Specifically, when two complex potentials are known a priori on a certain subinterval, reconstruction of the two potentials on the entire interval is studied… ▽ More In this work, we consider Dirac-type operators with a constant delay less than two-fifths of the interval and not less than one-third of the interval. For our considered Dirac-type operators, an incomplete inverse spectral problem is studied. Specifically, when two complex potentials are known a priori on a certain subinterval, reconstruction of the two potentials on the entire interval is studied from complete spectra of two boundary value problems with one common boundary condition. The uniqueness of the solution of the inverse problem is proved. A constructive method is developed for the solution of the inverse problem. △ Less

Submitted 21 May, 2023; originally announced May 2023.

Comments: arXiv admin note: substantial text overlap with arXiv:2305.10752

Inverse problems for Dirac operators with a constant delay less than half of the interval

Authors: Feng Wang, Chuan-Fu Yang

Abstract: In this work, we consider Dirac-type operators with a constant delay less than half of the interval and not less than two-fifths of the interval. For our considered Dirac-type operators, two inverse spectral problems are studied. Specifically, reconstruction of two complex $L_{2}$-potentials is studied from complete spectra of two boundary value problems with one common Dirichlet boundary conditio… ▽ More In this work, we consider Dirac-type operators with a constant delay less than half of the interval and not less than two-fifths of the interval. For our considered Dirac-type operators, two inverse spectral problems are studied. Specifically, reconstruction of two complex $L_{2}$-potentials is studied from complete spectra of two boundary value problems with one common Dirichlet boundary condition and Neumann boundary condition, respectively. We give answers to the full range of questions usually raised in the inverse spectral theory. That is, we give uniqueness, necessary and sufficient conditions of the solvability, reconstruction algorithm and uniform stability for our considered inverse problems. △ Less

Submitted 30 May, 2023; v1 submitted 18 May, 2023; originally announced May 2023.

Comments: arXiv admin note: text overlap with arXiv:2204.08259 by other authors

Global Existence of Weak Solutions for Compresssible Navier--Stokes--Fourier Equations with the Truncated Virial Pressure Law

Authors: Didier Bresch, Pierre-Emmanuel Jabin, Fei Wang

Abstract: This paper concerns the existence of global weak solutions {\it à la Leray} for compressible Navier--Stokes--Fourier systems with periodic boundary conditions and the truncated virial pressure law which is assumed to be thermodynamically unstable. More precisely, the main novelty is that the pressure law is not assumed to be monotone with respect to the density. This provides the first global weak… ▽ More This paper concerns the existence of global weak solutions {\it à la Leray} for compressible Navier--Stokes--Fourier systems with periodic boundary conditions and the truncated virial pressure law which is assumed to be thermodynamically unstable. More precisely, the main novelty is that the pressure law is not assumed to be monotone with respect to the density. This provides the first global weak solutions result for the compressible Navier-Stokes-Fourier system with such kind of pressure law which is strongly used as a generalization of the perfect gas law. The paper is based on a new construction of approximate solutions through an iterative scheme and fixed point procedure which could be very helpful to design efficient numerical schemes. Note that our method involves the recent paper by the authors published in Nonlinearity (2021) for the compactness of the density when the temperature is given. △ Less

Submitted 11 May, 2023; originally announced May 2023.

Kähler-Ricci flow on $\mathbf G$-spherical Fano manifolds

Authors: Feng Wang, Xiaohua Zhu

Abstract: We prove that the Gromov-Hausdorff limit of Kähler-Ricci flow on a $\mathbf G$-spherical Fano manifold $X$ is a $\mathbf G$-spherical $\mathbb Q$-Fano variety $X_{\infty}$, which admits a (singular) Kähler-Ricci soliton. Moreover, the $\mathbf G$-spherical variety structure of $X_{\infty}$ can be constructed as a center of torus $\mathbb C^*$-degeneration of $X$ induced by an element in the Lie al… ▽ More We prove that the Gromov-Hausdorff limit of Kähler-Ricci flow on a $\mathbf G$-spherical Fano manifold $X$ is a $\mathbf G$-spherical $\mathbb Q$-Fano variety $X_{\infty}$, which admits a (singular) Kähler-Ricci soliton. Moreover, the $\mathbf G$-spherical variety structure of $X_{\infty}$ can be constructed as a center of torus $\mathbb C^*$-degeneration of $X$ induced by an element in the Lie algebra of Cartan torus of $\mathbf G$. △ Less

Submitted 9 May, 2023; originally announced May 2023.

MSC Class: Primary: 53E20; 53C25; Secondary: 32Q20; 53C30; 14L10

The Hopf-Tsuji-Sullivan Dichotomy on Visibility Manifolds Without Conjugate Points

Authors: Fei Liu, Xiaokai Liu, Fang Wang

Abstract: In this article, we establish the Hopf-Tsuji-Sullivan dichotomy for geodesic flows on certain manifolds with no conjugate points: either the geodesic flow is conservative and ergodic, or it is completely dissipative and non-ergodic. We also show several equivalent conditions to the conservativity, such the Poincaré series diverges at the critical exponent, the conical limit set has full Patterson-… ▽ More In this article, we establish the Hopf-Tsuji-Sullivan dichotomy for geodesic flows on certain manifolds with no conjugate points: either the geodesic flow is conservative and ergodic, or it is completely dissipative and non-ergodic. We also show several equivalent conditions to the conservativity, such the Poincaré series diverges at the critical exponent, the conical limit set has full Patterson-Sullivan measure, etc. △ Less

Submitted 7 June, 2023; v1 submitted 29 April, 2023; originally announced May 2023.

Comments: 48 pages, 9 figures

MSC Class: 37D40 (Primary); 37D25 (Secondary)

Singular Degenerate SDEs: Well-Posedness and Exponential Ergodicity

Authors: Panpan Ren, Martin Grothaus, Feng-Yu Wang

Abstract: The well-posedness and exponential ergodicity are proved for stochastic Hamiltonian systems containing a singular drift term which is locally integrable in the component with noise. As an application, the well-posedness and uniform exponential ergodicity are derived for a class of singular degenerated McKean-Vlasov SDEs. The well-posedness and exponential ergodicity are proved for stochastic Hamiltonian systems containing a singular drift term which is locally integrable in the component with noise. As an application, the well-posedness and uniform exponential ergodicity are derived for a class of singular degenerated McKean-Vlasov SDEs. △ Less

Submitted 28 April, 2023; originally announced May 2023.

Comments: 26 pages

ACM Class: G.3

Probability Distance Estimates Between Diffusion Processes and Applications to Singular McKean-Vlasov SDEs

Authors: Xing Huang, Panpan Ren, Feng-Yu Wang

Abstract: The $L^k$-Wasserstein distance $\mathbb{W}_k (k\ge 1)$ and the probability distance $\mathbb{W}_ψ$ induced by a concave function $ψ$, are estimated between different diffusion processes with singular coefficients. As applications, the well-posedness, probability distance estimates and the log-Harnack inequality are derived for McKean-Vlasov SDEs with multiplicative distribution dependent noise, wh… ▽ More The $L^k$-Wasserstein distance $\mathbb{W}_k (k\ge 1)$ and the probability distance $\mathbb{W}_ψ$ induced by a concave function $ψ$, are estimated between different diffusion processes with singular coefficients. As applications, the well-posedness, probability distance estimates and the log-Harnack inequality are derived for McKean-Vlasov SDEs with multiplicative distribution dependent noise, where the coefficients are singular in time-space variables and $(\mathbb{W}_k+\mathbb{W}_ψ)$-Lipschitz continuous in the distribution variable. This improves existing results derived in the literature under the $\mathbb{W}_k$-Lipschitz or derivative conditions in the distribution variable. △ Less

Submitted 5 November, 2023; v1 submitted 15 April, 2023; originally announced April 2023.

Comments: 22 pages

Entropy Estimate Between Diffusion Processes with Application to Nonlinear Fokker-Planck Equations

Authors: Panpan Ren, Feng-Yu Wang

Abstract: By developing a new technique called the bi-coupling argument, we estimate the relative entropy between different diffusion processes in terms of the distances of initial distributions and drift-diffusion coefficients. As an application, the entropy-cost inequality is established for nonlinear Fokker-Planck equations, which has potential applications in optimal transport, information theory and me… ▽ More By developing a new technique called the bi-coupling argument, we estimate the relative entropy between different diffusion processes in terms of the distances of initial distributions and drift-diffusion coefficients. As an application, the entropy-cost inequality is established for nonlinear Fokker-Planck equations, which has potential applications in optimal transport, information theory and mean field particle systems. △ Less

Submitted 26 June, 2023; v1 submitted 26 February, 2023; originally announced February 2023.

Comments: 23 pages

Inducing Schemes with Finite Weighted Complexity

Authors: Jianyu Chen, Fang Wang, Hong-Kun Zhang

Abstract: In this paper, we consider a Borel measurable map of a compact metric space which admits an inducing scheme. Under the finite weighted complexity condition, we establish a thermodynamic formalism for a parameter family of potentials $\varphi+tψ$ in an interval containing $t=0$. Furthermore, if there is a generating partition compatible to the inducing scheme, we show that all ergodic invariant mea… ▽ More In this paper, we consider a Borel measurable map of a compact metric space which admits an inducing scheme. Under the finite weighted complexity condition, we establish a thermodynamic formalism for a parameter family of potentials $\varphi+tψ$ in an interval containing $t=0$. Furthermore, if there is a generating partition compatible to the inducing scheme, we show that all ergodic invariant measures with sufficiently large pressure are liftable. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Discrepancies of subtrees

Authors: Tarun Krishna, Peleg Michaeli, Michail Sarantis, Fenglin Wang, Yiqing Wang

Abstract: We study multicolour, oriented and high-dimensional discrepancies of the set of all subtrees of a tree. As our main result, we show that the $r$-colour discrepancy of the subtrees of any tree is a linear function of the number of leaves $\ell$ of that tree. More concretely, we show that it is bounded by $\lceil(r-1)\ell/r\rceil$ from below and $\lceil(r-1)\ell/2\rceil$ from above, and that these b… ▽ More We study multicolour, oriented and high-dimensional discrepancies of the set of all subtrees of a tree. As our main result, we show that the $r$-colour discrepancy of the subtrees of any tree is a linear function of the number of leaves $\ell$ of that tree. More concretely, we show that it is bounded by $\lceil(r-1)\ell/r\rceil$ from below and $\lceil(r-1)\ell/2\rceil$ from above, and that these bounds are asymptotically sharp. Motivated by this result, we introduce natural notions of oriented and high-dimensional discrepancies and prove bounds for the corresponding discrepancies of the set of all subtrees of a given tree as functions of its number of leaves. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: 12 pages

MSC Class: 11K38; 05C15; 05C05

Low progress math in a high performing system

Authors: A. Jamaludin, A. I. Jabir, F. J. Wang, A. L. Tan

Abstract: Math anxiety negatively relates to math performance. This negative relationship may be exacerbated in low-progress math learners. However, there are limited studies on math anxiety among low progress learners in a paradoxically high performing education system like Singapore. To fill this research gap, this research analysed the anxiety profiles of 151 students who were in the math learning suppor… ▽ More Math anxiety negatively relates to math performance. This negative relationship may be exacerbated in low-progress math learners. However, there are limited studies on math anxiety among low progress learners in a paradoxically high performing education system like Singapore. To fill this research gap, this research analysed the anxiety profiles of 151 students who were in the math learning support intervention program administered by the Ministry of Education, Singapore (MOE). We examined the complex relationship centred in math anxiety with relevant variables such as demographic characteristics, working memory and math performance. Limitations and future directions are discussed. △ Less

Submitted 25 January, 2023; originally announced February 2023.