-
Well-posedness and Bilinear Controllability of a Repairable System with Degraded State
Authors:
Daniel Owusu Adu,
Weiwei Hu
Abstract:
In this work, we consider the dynamics of repairable systems characterized by three distinct states: one signifying normal operational states, another representing degraded conditions and a third denoting failed conditions. These systems are characterized by their ability to be repaired when failures and/or degradation occur. Typically described by transport equations, these systems exhibit a coup…
▽ More
In this work, we consider the dynamics of repairable systems characterized by three distinct states: one signifying normal operational states, another representing degraded conditions and a third denoting failed conditions. These systems are characterized by their ability to be repaired when failures and/or degradation occur. Typically described by transport equations, these systems exhibit a coupled nature, interlinked through integro-differential equations and integral boundary conditions that dictate the transitions among all the states. In this paper, we address two less-explored facets: 1) the well-posedness and the asymptotic behavior of such systems with maximum repair time being finite; and 2) the bilinear controllability of the system via repair actions. In particular, we focus on the case where only one degraded and one failed states exist. We first discuss part 1) for given time-independent repair rates and then design the space-time dependent repair strategies that can manipulate system dynamics to achieve the desired level over a finite horizon. Our objective is to enhance the system availability -- the probability of being operational when needed over a fixed period of time. We present rigorous analysis and develop control strategies that leverage the bilinear structure of the system model.
△ Less
Submitted 6 July, 2024;
originally announced July 2024.
-
Bilinear Controllability of a Simple Reparable System
Authors:
Daniel Owusu Adu,
Weiwei Hu
Abstract:
Reparable systems are systems that are characterized by their ability to undergo maintenance actions when failures occur. These systems are often described by transport equations, all coupled through an integro-differential equation. In this paper, we address the understudied aspect of the controllability of reparable systems. In particular, we focus on a two-state reparable system and our goal is…
▽ More
Reparable systems are systems that are characterized by their ability to undergo maintenance actions when failures occur. These systems are often described by transport equations, all coupled through an integro-differential equation. In this paper, we address the understudied aspect of the controllability of reparable systems. In particular, we focus on a two-state reparable system and our goal is to design a control strategy that enhances the system availability -- the probability of being operational when needed. We establish bilinear controllability, demonstrating that appropriate control actions can manipulate system dynamics to achieve desired availability levels. We provide theoretical foundations and develop control strategies that leverage the bilinear structure of the equations.
△ Less
Submitted 6 July, 2024;
originally announced July 2024.
-
Boundedness of variation, oscillation and maximal differential transform on BMO space
Authors:
Wenting Hu,
Kai Wu,
Dongyong Yang,
Chao Zhang
Abstract:
In this paper, we prove that the oscillation operator, variation operator and maximal differential transform associated with the approximate identities are bounded from ${\rm BMO}({\mathbb R}^n)$ to its subspace ${\rm BLO}({\mathbb R}^n)$.
In this paper, we prove that the oscillation operator, variation operator and maximal differential transform associated with the approximate identities are bounded from ${\rm BMO}({\mathbb R}^n)$ to its subspace ${\rm BLO}({\mathbb R}^n)$.
△ Less
Submitted 4 June, 2024;
originally announced June 2024.
-
BO4IO: A Bayesian optimization approach to inverse optimization with uncertainty quantification
Authors:
Yen-An Lu,
Wei-Shou Hu,
Joel A. Paulson,
Qi Zhang
Abstract:
This work addresses data-driven inverse optimization (IO), where the goal is to estimate unknown parameters in an optimization model from observed decisions that can be assumed to be optimal or near-optimal solutions to the optimization problem. The IO problem is commonly formulated as a large-scale bilevel program that is notoriously difficult to solve. Deviating from traditional exact solution m…
▽ More
This work addresses data-driven inverse optimization (IO), where the goal is to estimate unknown parameters in an optimization model from observed decisions that can be assumed to be optimal or near-optimal solutions to the optimization problem. The IO problem is commonly formulated as a large-scale bilevel program that is notoriously difficult to solve. Deviating from traditional exact solution methods, we propose a derivative-free optimization approach based on Bayesian optimization, which we call BO4IO, to solve general IO problems. We treat the IO loss function as a black box and approximate it with a Gaussian process model. Using the predicted posterior function, an acquisition function is minimized at each iteration to query new candidate solutions and sequentially converge to the optimal parameter estimates. The main advantages of using Bayesian optimization for IO are two-fold: (i) it circumvents the need of complex reformulations of the bilevel program or specialized algorithms and can hence enable computational tractability even when the underlying optimization problem is nonconvex or involves discrete variables, and (ii) it allows approximations of the profile likelihood, which provide uncertainty quantification on the IO parameter estimates. We apply the proposed method to three computational case studies, covering different classes of forward optimization problems ranging from convex nonlinear to nonconvex mixed-integer nonlinear programs. Our extensive computational results demonstrate the efficacy and robustness of BO4IO to accurately estimate unknown model parameters from small and noisy datasets. In addition, the proposed profile likelihood analysis has proven to be effective in providing good approximations of the confidence intervals on the parameter estimates and assessing the identifiability of the unknown parameters.
△ Less
Submitted 28 May, 2024;
originally announced May 2024.
-
A fast low-rank inversion algorithm of dielectric matrix in GW approximation
Authors:
Zhengbang Zhou,
Huanhuan Ma,
Wentiao Wu,
Weiguo Gao,
Jinlong Yang,
Meiyue Shao,
Wei Hu
Abstract:
The dielectric response function and its inverse are crucial physical quantities in materials science. We propose an accurate and efficient strategy to invert the dielectric function matrix. The GW approximation, a powerful approach to accurately describe many-body excited states, is taken as an application to demonstrate accuracy and efficiency. We incorporate the interpolative separable density…
▽ More
The dielectric response function and its inverse are crucial physical quantities in materials science. We propose an accurate and efficient strategy to invert the dielectric function matrix. The GW approximation, a powerful approach to accurately describe many-body excited states, is taken as an application to demonstrate accuracy and efficiency. We incorporate the interpolative separable density fitting (ISDF) algorithm with Sherman--Morrison--Woodbury (SMW) formula to accelerate the inversion process by exploiting low-rank properties of dielectric function in plane-wave GW calculations. Our ISDF--SMW strategy produces accurate quasiparticle energies with $O(N_{\mathrm{r}}N_{\mathrm{e}}^2)$ computational cost $(N_{\mathrm{e}}$ is the number of electrons and $N_{\mathrm{r}}=100$--$1000N_{\mathrm{e}}$ is the number of grid points) with negligible small error of $0.03$ eV for both complex molecules and solids. This new strategy for inverting the dielectric matrix can be \(50\times\) faster than the current state-of-the-art implementation in BerkeleyGW, resulting in two orders of magnitude speedup for total GW calculations.
△ Less
Submitted 12 June, 2024; v1 submitted 18 March, 2024;
originally announced March 2024.
-
Entropy of axial product of multiplicative subshifts
Authors:
Jung-Chao Ban,
Wen-Guei Hu,
Guan-Yu Lai,
Lingmin Liao
Abstract:
We obtain the entropy and the surface entropy of the axial products on $\mathbb{N}^d$ and the $d$-tree $T^d$ of two types of systems: the subshift and the multiplicative subshift.
We obtain the entropy and the surface entropy of the axial products on $\mathbb{N}^d$ and the $d$-tree $T^d$ of two types of systems: the subshift and the multiplicative subshift.
△ Less
Submitted 29 February, 2024;
originally announced February 2024.
-
Hausdorff dimensions of affine multiplicative subshifts
Authors:
Jung-Chao Ban,
Wen-Guei Hu,
Guan-Yu Lai,
Lingmin Liao
Abstract:
We calculate the Minkowski and Hausdorff dimensions of affine multiplicative subshifts on $\mathbb{N}$.
We calculate the Minkowski and Hausdorff dimensions of affine multiplicative subshifts on $\mathbb{N}$.
△ Less
Submitted 28 February, 2024;
originally announced February 2024.
-
Hausdorff dimension of random attractors for a stochastic delayed parabolic equation in Banach spaces
Authors:
Wenjie Hu,
TomásCaraballo,
Yueliang Duan
Abstract:
The main purpose of this paper is to give an upper bound of Hausdorff dimension of random attractors for a stochastic delayed parabolic equation in Banach spaces. The estimation of dimensions of random attractors are obtained by combining the squeezing property and a covering lemma of finite subspace of Banach spaces, which generalizes the method established in Hilbert spaces. Unlike the existing…
▽ More
The main purpose of this paper is to give an upper bound of Hausdorff dimension of random attractors for a stochastic delayed parabolic equation in Banach spaces. The estimation of dimensions of random attractors are obtained by combining the squeezing property and a covering lemma of finite subspace of Banach spaces, which generalizes the method established in Hilbert spaces. Unlike the existing works, where orthogonal projectors with finite ranks applied for proving the squeezing property of stochastic partial differential equations in Hilbert spaces, we adopt the state decomposition of phase space based on the exponential dichotomy of the the linear deterministic part of the studied SDPE to obtain similar squeezing property due to the lack of smooth inner product geometry structure. The obtained dimension of the random attractors depend only on the spectrum of the linear part and the random Lipschitz constant of the nonlinear term, while not relating to the compact embedding of the phase space to another Banach space as the existing works did.
△ Less
Submitted 19 February, 2024;
originally announced February 2024.
-
Exponential attractors for a nonlocal delayed reaction-diffusion equation on an unbounded domain
Authors:
Wenjie Hu,
Tomás Caraballo
Abstract:
The main objective of this paper is to investigate exponential attractors for a nonlocal delayed reaction-diffusion equation on an unbounded domain. We first obtain the existence of a globally attractive absorbing set for the dynamical system generated by the equation under the assumption that the nonlinear term is bounded. Then, we construct exponential attractors of the equation directly in its…
▽ More
The main objective of this paper is to investigate exponential attractors for a nonlocal delayed reaction-diffusion equation on an unbounded domain. We first obtain the existence of a globally attractive absorbing set for the dynamical system generated by the equation under the assumption that the nonlinear term is bounded. Then, we construct exponential attractors of the equation directly in its natural phase space, i.e., a Banach space with explicit fractal dimension by combining squeezing properties of the system as well as a covering lemma of finite subspace of Banach spaces. Our result generalizes the methods established in Hilbert spaces and weighted spaces, and the fractal dimension of the obtained exponential attractor does not depend on the entropy number but only depends on some inner characteristic of the studied equation.
△ Less
Submitted 19 February, 2024;
originally announced February 2024.
-
A constructive approach to selective risk control
Authors:
Zijun Gao,
Wenjie Hu,
Qingyuan Zhao
Abstract:
Many modern applications require the use of data to both select the statistical tasks and make valid inference after selection. In this article, we provide a unifying approach to control for a class of selective risks. Our method is motivated by a reformulation of the celebrated Benjamini-Hochberg (BH) procedure for multiple hypothesis testing as the iterative limit of the Benjamini-Yekutieli (BY)…
▽ More
Many modern applications require the use of data to both select the statistical tasks and make valid inference after selection. In this article, we provide a unifying approach to control for a class of selective risks. Our method is motivated by a reformulation of the celebrated Benjamini-Hochberg (BH) procedure for multiple hypothesis testing as the iterative limit of the Benjamini-Yekutieli (BY) procedure for constructing post-selection confidence intervals. Although several earlier authors have made noteworthy observations related to this, our discussion highlights that (1) the BH procedure is precisely the fixed-point iteration of the BY procedure; (2) the fact that the BH procedure controls the false discovery rate is almost an immediate corollary of the fact that the BY procedure controls the false coverage-statement rate. Building on this observation, we propose a constructive approach to control extra-selection risk (selection made after decision) by iterating decision strategies that control the post-selection risk (decision made after selection), and show that many previous methods and results are special cases of this general framework. We further extend this approach to problems with multiple selective risks and demonstrate how new methods can be developed. Our development leads to two surprising results about the BH procedure: (1) in the context of one-sided location testing, the BH procedure not only controls the false discovery rate at the null but also at other locations for free; (2) in the context of permutation tests, the BH procedure with exact permutation p-values can be well approximated by a procedure which only requires a total number of permutations that is almost linear in the total number of hypotheses.
△ Less
Submitted 29 January, 2024;
originally announced January 2024.
-
Shuffle Hopf algebra of Multiple zeta values
Authors:
Wenchuan Hu,
Hongyu Xiang,
Bin Zhang
Abstract:
The shuffle relation among multiple zeta values is algebraically expressed as the shuffle algebra. In this paper, the shuffle algebra structure for multiple zeta values is extended to a Hopf algebra structure, for which the key idea is the lifting of the shuffle multiplication to Chen fractions as the function multiplication. The linear span of Chen fractions can be equipped with a locality Hopf a…
▽ More
The shuffle relation among multiple zeta values is algebraically expressed as the shuffle algebra. In this paper, the shuffle algebra structure for multiple zeta values is extended to a Hopf algebra structure, for which the key idea is the lifting of the shuffle multiplication to Chen fractions as the function multiplication. The linear span of Chen fractions can be equipped with a locality Hopf algebra structure, and the pushforward of the coproduct gives us the desired construction on the shuffle algebra.
△ Less
Submitted 29 December, 2023;
originally announced January 2024.
-
Learning Free Terminal Time Optimal Closed-loop Control of Manipulators
Authors:
Wei Hu,
Yue Zhao,
Weinan E,
Jiequn Han,
Jihao Long
Abstract:
This paper presents a novel approach to learning free terminal time closed-loop control for robotic manipulation tasks, enabling dynamic adjustment of task duration and control inputs to enhance performance. We extend the supervised learning approach, namely solving selected optimal open-loop problems and utilizing them as training data for a policy network, to the free terminal time scenario. Thr…
▽ More
This paper presents a novel approach to learning free terminal time closed-loop control for robotic manipulation tasks, enabling dynamic adjustment of task duration and control inputs to enhance performance. We extend the supervised learning approach, namely solving selected optimal open-loop problems and utilizing them as training data for a policy network, to the free terminal time scenario. Three main challenges are addressed in this extension. First, we introduce a marching scheme that enhances the solution quality and increases the success rate of the open-loop solver by gradually refining time discretization. Second, we extend the QRnet in Nakamura-Zimmerer et al. (2021b) to the free terminal time setting to address discontinuity and improve stability at the terminal state. Third, we present a more automated version of the initial value problem (IVP) enhanced sampling method from previous work (Zhang et al., 2022) to adaptively update the training dataset, significantly improving its quality. By integrating these techniques, we develop a closed-loop policy that operates effectively over a broad domain with varying optimal time durations, achieving near globally optimal total costs.
△ Less
Submitted 29 November, 2023;
originally announced November 2023.
-
Rigidity dimensions of self-injective Nakayama algebras
Authors:
Wei Hu,
Xiaojuan Yin
Abstract:
Rigidity dimension is a new homological dimension which is intended to measure the quality of the best resolution of an algebra. In this paper, we determine the rigidity dimensions of self-injective Nakayama agebras A_{n,m} with n simple modules and the Loewy length m>=n.
Rigidity dimension is a new homological dimension which is intended to measure the quality of the best resolution of an algebra. In this paper, we determine the rigidity dimensions of self-injective Nakayama agebras A_{n,m} with n simple modules and the Loewy length m>=n.
△ Less
Submitted 21 November, 2023;
originally announced November 2023.
-
Existence and dimensions of global attractors for a delayed reaction-diffusion equation on an unbounded domain
Authors:
Wenjie Hu,
Tomás Caraballo,
Alain Miranville
Abstract:
The purpose of this paper is to investigate the existence and Hausdorff dimension as well as fractal dimension of global attractors for a delayed reaction-diffusion equation on an unbounded domain. The noncompactness of the domain causes the Laplace operator has a continuous spectrum, the semigroup generated by the linear part and the Sobolev embeddings are no longer compact, making the problem mo…
▽ More
The purpose of this paper is to investigate the existence and Hausdorff dimension as well as fractal dimension of global attractors for a delayed reaction-diffusion equation on an unbounded domain. The noncompactness of the domain causes the Laplace operator has a continuous spectrum, the semigroup generated by the linear part and the Sobolev embeddings are no longer compact, making the problem more difficult compared with the equations on bounded domains. We first obtain the existence of an absorbing set for the infinite dimensional dynamical system generated by the equation by a priori estimate of the solutions. Then, we show the asymptotic compactness of the solution semiflow by an uniform a priori estimates for far-field values of solutions together with the Arzelà-Ascoli theorem, which facilitates us to show the existence of global attractors. By decomposing the solution into three parts and establishing a squeezing property of each part, we obtain the explicit upper estimation of both Hausdorff and fractal dimension of the global attractors, which only depend on the inner characteristic of the equation, while not related to the entropy number compared with the existing literature.
△ Less
Submitted 16 November, 2023;
originally announced November 2023.
-
A branch-and-price approach for the nurse rostering problem with multiple units
Authors:
Wanzhe Hu,
Xiaozhou He,
Li Luo,
Panos M. Pardalos
Abstract:
In this paper, we study the nurse rostering problem that considers multiple units and many soft time-related constraints. An efficient branch and price solution approach that relies on a fast algorithm to solve the pricing subproblem of the column generation process is presented. For the nurse rostering problem, its pricing subproblem can be formulated as a shortest path problem with resource cons…
▽ More
In this paper, we study the nurse rostering problem that considers multiple units and many soft time-related constraints. An efficient branch and price solution approach that relies on a fast algorithm to solve the pricing subproblem of the column generation process is presented. For the nurse rostering problem, its pricing subproblem can be formulated as a shortest path problem with resource constraints, which has been the backbone of several solutions for several classical problems like vehicle routing problems. However, approaches that perform well on these problems cannot be used since most constraints in the nurse rostering problem are soft. Based on ideas borrowed from global constraints in constraint programming to model rostering problems, an efficient dynamic programming algorithm with novel label definitions and dominating rules specific to soft time-related constraints is proposed. In addition, several acceleration strategies are employed to improve the branch and price algorithm. Computational results on instances of different sizes indicate that the proposed algorithm is a promising solution for the nurse rostering problem with multiple units.
△ Less
Submitted 9 November, 2023;
originally announced November 2023.
-
Understanding Deep Representation Learning via Layerwise Feature Compression and Discrimination
Authors:
Peng Wang,
Xiao Li,
Can Yaras,
Zhihui Zhu,
Laura Balzano,
Wei Hu,
Qing Qu
Abstract:
Over the past decade, deep learning has proven to be a highly effective tool for learning meaningful features from raw data. However, it remains an open question how deep networks perform hierarchical feature learning across layers. In this work, we attempt to unveil this mystery by investigating the structures of intermediate features. Motivated by our empirical findings that linear layers mimic…
▽ More
Over the past decade, deep learning has proven to be a highly effective tool for learning meaningful features from raw data. However, it remains an open question how deep networks perform hierarchical feature learning across layers. In this work, we attempt to unveil this mystery by investigating the structures of intermediate features. Motivated by our empirical findings that linear layers mimic the roles of deep layers in nonlinear networks for feature learning, we explore how deep linear networks transform input data into output by investigating the output (i.e., features) of each layer after training in the context of multi-class classification problems. Toward this goal, we first define metrics to measure within-class compression and between-class discrimination of intermediate features, respectively. Through theoretical analysis of these two metrics, we show that the evolution of features follows a simple and quantitative pattern from shallow to deep layers when the input data is nearly orthogonal and the network weights are minimum-norm, balanced, and approximate low-rank: Each layer of the linear network progressively compresses within-class features at a geometric rate and discriminates between-class features at a linear rate with respect to the number of layers that data have passed through. To the best of our knowledge, this is the first quantitative characterization of feature evolution in hierarchical representations of deep linear networks. Empirically, our extensive experiments not only validate our theoretical results numerically but also reveal a similar pattern in deep nonlinear networks which aligns well with recent empirical studies. Moreover, we demonstrate the practical implications of our results in transfer learning. Our code is available at \url{https://github.com/Heimine/PNC_DLN}.
△ Less
Submitted 9 January, 2024; v1 submitted 6 November, 2023;
originally announced November 2023.
-
A new mixed finite element method for arbitrary element pair for a quasi-static nonlinear permeability thermo-poroelasticity model
Authors:
Zhihao Ge,
Wenshuai Hu
Abstract:
In this paper, we develop a multiphysics finite element method for solving the quasi-static thermo-poroelasticity model with nonlinear permeability. The model involves multiple physical processes such as deformation, pressure, diffusion and heat transfer. To reveal the multi-physical processes of deformation, diffusion and heat transfer, we reformulate the original model into a fluid coupled probl…
▽ More
In this paper, we develop a multiphysics finite element method for solving the quasi-static thermo-poroelasticity model with nonlinear permeability. The model involves multiple physical processes such as deformation, pressure, diffusion and heat transfer. To reveal the multi-physical processes of deformation, diffusion and heat transfer, we reformulate the original model into a fluid coupled problem that is general Stokes equation coupled with two reaction-diffusion equations. Then, we prove the existence and uniqueness of weak solution for the original problem by the $B$-operator technique and by sequence approximation for the reformulated problem. As for the reformulated problem we propose a fully discrete finite element method which can use arbitrary finite element pairs to solve the displacement $\bu$ pressure $τ$ and variable $\varpi,ς$, and the backward Euler method for time discretization. Finally, we give the stability analysis of the above proposed method, also we prove that the fully discrete multiphysics finite element method has an optimal convergence order. Numerical experiments show that the proposed method can achieve good results under different finite element pairs and are consistent with the theoretical analysis.
△ Less
Submitted 10 October, 2023;
originally announced October 2023.
-
Simulating time-harmonic acoustic wave effects induced by periodic holes/inclusions on surfaces
Authors:
Wen Hu,
Zhuojia Fu,
Leevan Ling
Abstract:
This paper introduces the first attempt to employ a localized meshless method to analyze time-harmonic acoustic wave propagation on curved surfaces with periodic holes/inclusions. In particular, the generalized finite difference method is used as a localized meshless technique to discretize the surface gradient and Laplace-Beltrami operators defined extrinsically in the governing equations. An abs…
▽ More
This paper introduces the first attempt to employ a localized meshless method to analyze time-harmonic acoustic wave propagation on curved surfaces with periodic holes/inclusions. In particular, the generalized finite difference method is used as a localized meshless technique to discretize the surface gradient and Laplace-Beltrami operators defined extrinsically in the governing equations. An absorbing boundary condition is introduced to reduce reflections from boundaries and accurately simulate wave propagation on unclosed surfaces with periodic inclusions. Several benchmark examples demonstrate the efficiency and accuracy of the proposed method in simulating acoustic wave propagation on surfaces with diverse geometries, including complex shapes and periodic holes or inclusions.
△ Less
Submitted 19 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 hybrid neural-network and MAC scheme for Stokes interface problems
Authors:
Che-Chia Chang,
Chen-Yang Dai,
Wei-Fan Hu,
Te-Sheng Lin,
Ming-Chih Lai
Abstract:
In this paper, we present a hybrid neural-network and MAC (Marker-And-Cell) scheme for solving Stokes equations with singular forces on an embedded interface in regular domains. As known, the solution variables (the pressure and velocity) exhibit non-smooth behaviors across the interface so extra discretization efforts must be paid near the interface in order to have small order of local truncatio…
▽ More
In this paper, we present a hybrid neural-network and MAC (Marker-And-Cell) scheme for solving Stokes equations with singular forces on an embedded interface in regular domains. As known, the solution variables (the pressure and velocity) exhibit non-smooth behaviors across the interface so extra discretization efforts must be paid near the interface in order to have small order of local truncation errors in finite difference schemes. The present hybrid approach avoids such additional difficulty. It combines the expressive power of neural networks with the convergence of finite difference schemes to ease the code implementation and to achieve good accuracy at the same time. The key idea is to decompose the solution into singular and regular parts. The neural network learning machinery incorporating the given jump conditions finds the singular part solution, while the standard MAC scheme is used to obtain the regular part solution with associated boundary conditions. The two- and three-dimensional numerical results show that the present hybrid method converges with second-order accuracy for the velocity and first-order accuracy for the pressure, and it is comparable with the traditional immersed interface method in literature.
△ Less
Submitted 3 April, 2024; v1 submitted 9 June, 2023;
originally announced June 2023.
-
Large deviation principle of Multiplicative Ising models on Markov-Cayley Trees
Authors:
Jung-Chao Ban,
Wen-Guei Hu,
Zongfan Zhang
Abstract:
In this paper, we study the large deviation principle (LDP) for two types (Type I and Type II) of multiplicative Ising models. For Types I and II, the explicit formulas for the free energy functions and the associated rate functions are derived. Furthermore, we prove that those free energy functions are differentiable, which indicates that both systems are characterized by a lack of phase transiti…
▽ More
In this paper, we study the large deviation principle (LDP) for two types (Type I and Type II) of multiplicative Ising models. For Types I and II, the explicit formulas for the free energy functions and the associated rate functions are derived. Furthermore, we prove that those free energy functions are differentiable, which indicates that both systems are characterized by a lack of phase transition phenomena.
△ Less
Submitted 15 May, 2023;
originally announced May 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.
-
The entropy structures of axial products on $\mathbb{N}^d$ and Trees
Authors:
Jung-Chao Ban,
Wen-Guei Hu,
Guan-Yu Lai
Abstract:
In this paper, we first concentrate on the possible values and dense property of entropies for isotropic and anisotropic axial products of subshifts of finite type (SFTs) on $\mathbb{N}^d$ and $d$-tree $\mathcal{T}_d$. We prove that the entropies of isotropic and anisotropic axial products of SFTs on $\mathbb{N}^d$ are dense in $[0,\infty)$, and the same result also holds for anisotropic axial pro…
▽ More
In this paper, we first concentrate on the possible values and dense property of entropies for isotropic and anisotropic axial products of subshifts of finite type (SFTs) on $\mathbb{N}^d$ and $d$-tree $\mathcal{T}_d$. We prove that the entropies of isotropic and anisotropic axial products of SFTs on $\mathbb{N}^d$ are dense in $[0,\infty)$, and the same result also holds for anisotropic axial products of SFTs on $\mathcal{T}_d$. However, the result is no longer true for isotropic axial products of SFTs on $\mathcal{T}_d$. Next, motivated by the work of Johnson, Kass and Madden [16], and Schraudner [28], we establish the entropy formula and structures for full axial extension shifts on $\mathbb{N}^d$ and $\mathcal{T}_d$. Combining the aforementioned results with the findings on the surface entropy for multiplicative integer systems [8] on $\mathbb{N}^d$ enables us to estimate the surface entropy for the full axial extension shifts on $\mathcal{T}_d$. Finally, we extend the results of full axial extension shifts on $\mathcal{T}_d$ to general trees.
△ Less
Submitted 22 March, 2023;
originally announced March 2023.
-
Exponential attractors with explicit fractal dimensions for retarded functional differential equations
Authors:
Wenjie Hu,
Tomás Caraballo
Abstract:
The aim of this paper is to propose a new method to construct exponential attractors for infinite dimensional dynamical systems in Banach spaces with explicit fractal dimension. The approach is established by combing the squeezing properties and the covering of finite subspace of Banach spaces, which generalize the method established in Hilbert spaces. The constructed exponential attractors posses…
▽ More
The aim of this paper is to propose a new method to construct exponential attractors for infinite dimensional dynamical systems in Banach spaces with explicit fractal dimension. The approach is established by combing the squeezing properties and the covering of finite subspace of Banach spaces, which generalize the method established in Hilbert spaces. The constructed exponential attractors possess explicit fractal dimensions which do not depend on the entropy number but only depend on the spectrum of the linear part and Lipschitz constant of the nonlinear part. The method is especially effective for functional differential equations in Banach spaces for which sate decomposition of the linear part can be adopted to prove squeezing property. The theoretical results are applied to the retarded functional differential equation and retarded reaction-diffusion equations.
△ Less
Submitted 3 May, 2023; v1 submitted 7 March, 2023;
originally announced March 2023.
-
Lyapunov exponents and invariant manifolds for stochastic linear partial functional differential equations
Authors:
Wenjie Hu,
Tomás Caraballo
Abstract:
The main purpose of this work is to characterize the almost sure local structure stability of solutions to a class of linear stochastic partial functional differential equations (SPFDEs) by investigating the Lyapunov exponents and invariant manifolds near the stationary point. It is firstly proved that the trajectory field of the stochastic delayed stochastic partial functional differential equati…
▽ More
The main purpose of this work is to characterize the almost sure local structure stability of solutions to a class of linear stochastic partial functional differential equations (SPFDEs) by investigating the Lyapunov exponents and invariant manifolds near the stationary point. It is firstly proved that the trajectory field of the stochastic delayed stochastic partial functional differential equation admits an almost sure continuous version which is compact for $t>τ$ by a delicate construction based on the random semiflow generated by the diffusion term. Then it is proved that the version generates a random dynamical system(RDS) by the Wong-Zakai approximation of the stochastic partial differential equation constructed by the diffusion term. Subsequently, it is shown that the constructed linear cocycle admits fixed (at most) countable set of Lyapunov exponents and the associate Oseledets random filtration of the Banach space is obtained by adopting the infinite-dimensional multiplicative ergodic theorem in Banach spaces established by Lian and Lu [\textit{Mem Amer Math Soc, 2010, 206: 967}]. As a by product, the stable-manifolds theorem for the linear SPFDE in the hyperbolic case is also established.
△ Less
Submitted 19 October, 2023; v1 submitted 7 March, 2023;
originally announced March 2023.
-
Topological dimensions of attractors for partial functional differential equations in Banach spaces
Authors:
Wenjie Hu,
Tomás Caraballo
Abstract:
The main objective of this paper is to obtain estimations of Hausdorff dimension as well as fractal dimension of global attractors and pullback attractors for both autonomous and nonautonomous functional differential equations (FDEs) in Banach spaces. New criterions for the finite Hausdorff dimension and fractal dimension of attractors in Banach spaces are firslty proposed by combining the squeezi…
▽ More
The main objective of this paper is to obtain estimations of Hausdorff dimension as well as fractal dimension of global attractors and pullback attractors for both autonomous and nonautonomous functional differential equations (FDEs) in Banach spaces. New criterions for the finite Hausdorff dimension and fractal dimension of attractors in Banach spaces are firslty proposed by combining the squeezing property and the covering of finite subspace of Banach spaces, which generalize the method established in Hilbert spaces. In order to surmount the barrier caused by the lack of orthogonal projectors with finite rank, which is the key tool for proving the squeezing property of partial differential equations in Hilbert spaces, we adopt the state decomposition of phase space based on the exponential dichotomy of the studied FDEs to obtain similar squeezing property. The theoretical results are applied to a retarded nonlinear reaction-diffusion equation and a non-autonomous retarded functional differential equation in the natural phase space, for which explicit bounds of dimensions that do not depend on the entropy number but only depend on the spectrum of the linear parts and Lipschitz constants of the nonlinear parts are obtained.
△ Less
Submitted 19 November, 2023; v1 submitted 7 March, 2023;
originally announced March 2023.
-
Topological dimensions of global attractors for a delayed reaction-diffusion equation on an unbounded domain
Authors:
Wenjie Hu,
Tomás Caraballo
Abstract:
The purpose of this paper is to investigate the existence and estimation of Hausdorff and fractal dimension of global attractors for a delayed reaction-diffusion equation on an unbounded domain. The noncompactness of the domain cause the Laplace operator has a continuous spectrum, the semigroup generated by the linear part and the Sobolev embeddings no longer compact, making the problem more diffi…
▽ More
The purpose of this paper is to investigate the existence and estimation of Hausdorff and fractal dimension of global attractors for a delayed reaction-diffusion equation on an unbounded domain. The noncompactness of the domain cause the Laplace operator has a continuous spectrum, the semigroup generated by the linear part and the Sobolev embeddings no longer compact, making the problem more difficult compared with the equation on bounded domain. In order to adopt Hausdorff and fractal dimension dimension estimation tools established for dynamical systems in Hilbert space, we recast the equation in an auxiliary Hilbert space. We first obtain the existence of global solutions on unbounded domain by a perturbation semigroup approach which generates an infinite dimensional dynamical system. Then, we show the existence of global attractors by firstly showing the existence of an absorbing set and then give a an uniform a priori estimates for far-field values of solutions which facilitate us to prove the asymptotic compactness of the generated dynamical system. After establishing the variational equation in the auxiliary Hilbert space and the differentiable properties of the generated dynamical system, the upper estimate of both Hausdorff and fractal dimensions of the global attractors is obtained.
△ Less
Submitted 19 October, 2023; v1 submitted 7 March, 2023;
originally announced March 2023.
-
Random attractors of a stochastic Hopfield neural network model with delays
Authors:
Wenjie Hu,
Quanxin Zhu,
Peter E. Kloeden
Abstract:
The global asymptotic behavior of a stochastic Hopfield neural network model (HNNM) with delays is explored by studying the existence and structure of random attractors. It is first proved that the trajectory field of the stochastic delayed HNNM admits an almost sure continuous version, which is compact for $t>τ$ (where $τ$ is the delay) by a delicate construction based on the random semiflow gene…
▽ More
The global asymptotic behavior of a stochastic Hopfield neural network model (HNNM) with delays is explored by studying the existence and structure of random attractors. It is first proved that the trajectory field of the stochastic delayed HNNM admits an almost sure continuous version, which is compact for $t>τ$ (where $τ$ is the delay) by a delicate construction based on the random semiflow generated by the diffusion term. Then, this version is shown to generate a random dynamical system (RDS) by piece-wise linear approximation, after which the existence of a random absorbing set is obtained by a careful uniform apriori estimate of the solutions. Subsequently, the pullback asymptotic compactness of the RDS generated by the stochastic delayed HNNM is proved and hence the existence of random attractors is obtained. Moreover, sufficient conditions under which the attractors turn out to be an exponential attracting stationary solution are given. Numerical simulations are also conducted at last to illustrate the effectiveness of the established results.
△ Less
Submitted 10 February, 2023;
originally announced February 2023.
-
Invariant manifolds for stochastic delayed partial differential equations of parabolic type
Authors:
Wenjie Hu,
Quanxin Zhu,
Tomás Caraballo
Abstract:
The aim of this paper is to prove the existence and smoothness of stable and unstable invariant manifolds for a stochastic delayed partial differential equation of parabolic type. The stochastic delayed partial differential equation is firstly transformed into a random delayed partial differential equation by a conjugation, which is then recast into a Hilbert space. For the auxiliary equation, the…
▽ More
The aim of this paper is to prove the existence and smoothness of stable and unstable invariant manifolds for a stochastic delayed partial differential equation of parabolic type. The stochastic delayed partial differential equation is firstly transformed into a random delayed partial differential equation by a conjugation, which is then recast into a Hilbert space. For the auxiliary equation, the variation of constants formula holds and we show the existence of Lipschitz continuous stable and unstable manifolds by the Lyapunov-Perron method. Subsequently, we prove the smoothness of these invariant manifolds under appropriate spectral gap condition by carefully investigating the smoothness of auxiliary equation, after which, we obtain the invariant manifolds of the original equation by projection and inverse transformation. Eventually, we illustrate the obtained theoretical results by their application to a stochastic single-species population model.
△ Less
Submitted 12 June, 2023; v1 submitted 10 February, 2023;
originally announced February 2023.
-
Topological dimensions of random attractors for stochastic partial differential equations with delay
Authors:
Wenjie Hu,
Tomás Caraballo
Abstract:
The aim of this paper is to obtain an estimation of Hausdorff as well as fractal dimensions of random attractors for a class of stochastic partial differential equations with delay. The stochastic equation is first transformed into a delayed random partial differential equation by means of a random conjugation, which is then recast into an auxiliary Hilbert space. For the obtained equation, it is…
▽ More
The aim of this paper is to obtain an estimation of Hausdorff as well as fractal dimensions of random attractors for a class of stochastic partial differential equations with delay. The stochastic equation is first transformed into a delayed random partial differential equation by means of a random conjugation, which is then recast into an auxiliary Hilbert space. For the obtained equation, it is firstly proved that it generates a random dynamical system (RDS) in the auxiliary Hilbert space. Then it is shown that the equation possesses random attractors by a uniform estimate of the solution and the asymptotic compactness of the generated RDS. After establishing the variational equation in the auxiliary Hilbert space and the $\mathbb{P}$ almost surely differentiable properties of the RDS, an upper estimate of both Hausdorff and fractal dimensions of the random attractors are obtained.
△ Less
Submitted 10 February, 2023;
originally announced February 2023.
-
Random attractors for a stochastic nonlocal delayed reaction-diffusion equation on a semi-infinite interval
Authors:
Wenjie Hu,
Quanxin Zhu,
Tomás Caraballo
Abstract:
The aim of this paper is to prove the existence and qualitative property of random attractors for a stochastic nonlocal delayed reaction-diffusion equation (SNDRDE) on a semi-infinite interval with a Dirichlet boundary condition on the finite end. This equation models the spatial-temporal evolution of the mature individuals for a two-stage species whose juvenile and adults both diffuse that lives…
▽ More
The aim of this paper is to prove the existence and qualitative property of random attractors for a stochastic nonlocal delayed reaction-diffusion equation (SNDRDE) on a semi-infinite interval with a Dirichlet boundary condition on the finite end. This equation models the spatial-temporal evolution of the mature individuals for a two-stage species whose juvenile and adults both diffuse that lives on a semi-infinite domain and subject to random perturbations. By transforming the SNDRDE into a random evolution equation with delay, by means of a stationary conjugate transformation, we first establish the global existence and uniqueness of solutions to the equation, after which we show the solutions generate a random dynamical system. Then, we deduce uniform a priori estimates of the solutions and show the existence of bounded random absorbing sets. Subsequently, we prove the pullback asymptotic compactness of the random dynamical system generated by the SNDRDE with respect to the compact open topology, and hence obtain the existence of random attractors. At last, it is proved that the random attractor is an exponentially attracting stationary solution under appropriate conditions.
△ Less
Submitted 10 February, 2023;
originally announced February 2023.
-
Typed topology and its application to data set
Authors:
Wanjun Hu
Abstract:
The concept of $typed$ $topology$ is introduced. In a typed topological space, some open sets are assigned "types", and topological concepts such as closure, connectedness can be defined using types. A finite data set in $R^2$ is a typically typed topological space. Clusters calculated by the DBSCAN algorithm for data clustering can be well represent in a finite typed topological space. Other conc…
▽ More
The concept of $typed$ $topology$ is introduced. In a typed topological space, some open sets are assigned "types", and topological concepts such as closure, connectedness can be defined using types. A finite data set in $R^2$ is a typically typed topological space. Clusters calculated by the DBSCAN algorithm for data clustering can be well represent in a finite typed topological space. Other concepts such as tracks, port (starting points), type-p-connectedness, p-closure-connectedness, indexing, branches are also introduced for a finite typed topological space. Finally, $left-r$ and $up-left-r$ type open sets are introduced for data sets in $R^2$, so that tracks, port, branches can be calculated.
△ Less
Submitted 10 February, 2024; v1 submitted 10 February, 2023;
originally announced February 2023.
-
LS-DYNA Machine Learning-based Multiscale Method for Nonlinear Modeling of Short Fiber-Reinforced Composites
Authors:
Haoyan Wei,
C. T. Wu,
Wei Hu,
Tung-Huan Su,
Hitoshi Oura,
Masato Nishi,
Tadashi Naito,
Stan Chung,
Leo Shen
Abstract:
Short-fiber-reinforced composites (SFRC) are high-performance engineering materials for lightweight structural applications in the automotive and electronics industries. Typically, SFRC structures are manufactured by injection molding, which induces heterogeneous microstructures, and the resulting nonlinear anisotropic behaviors are challenging to predict by conventional micromechanical analyses.…
▽ More
Short-fiber-reinforced composites (SFRC) are high-performance engineering materials for lightweight structural applications in the automotive and electronics industries. Typically, SFRC structures are manufactured by injection molding, which induces heterogeneous microstructures, and the resulting nonlinear anisotropic behaviors are challenging to predict by conventional micromechanical analyses. In this work, we present a machine learning-based multiscale method by integrating injection molding-induced microstructures, material homogenization, and Deep Material Network (DMN) in the finite element simulation software LS-DYNA for structural analysis of SFRC. DMN is a physics-embedded machine learning model that learns the microscale material morphologies hidden in representative volume elements of composites through offline training. By coupling DMN with finite elements, we have developed a highly accurate and efficient data-driven approach, which predicts nonlinear behaviors of composite materials and structures at a computational speed orders-of-magnitude faster than the high-fidelity direct numerical simulation. To model industrial-scale SFRC products, transfer learning is utilized to generate a unified DMN database, which effectively captures the effects of injection molding-induced fiber orientations and volume fractions on the overall composite properties. Numerical examples are presented to demonstrate the promising performance of this LS-DYNA machine learning-based multiscale method for SFRC modeling.
△ Less
Submitted 6 January, 2023;
originally announced January 2023.
-
Liouville--type Theorems for Steady MHD and Hall--MHD Equations in $\R^2 \times \T$
Authors:
Wentao Hu,
Zhengce Zhang
Abstract:
In this paper, we study the Liouville--type theorems for three--dimensional stationary incompressible MHD and Hall--MHD systems in a slab with periodic boundary condition. We show that, under the assumptions that $(u^θ,b^θ)$ or $(u^r,b^r)$ is axisymmetric, or $(ru^r,rb^r)$ is bounded, any smooth bounded solution to the MHD or Hall--MHD system with local Dirichlet integral growing as an arbitrary p…
▽ More
In this paper, we study the Liouville--type theorems for three--dimensional stationary incompressible MHD and Hall--MHD systems in a slab with periodic boundary condition. We show that, under the assumptions that $(u^θ,b^θ)$ or $(u^r,b^r)$ is axisymmetric, or $(ru^r,rb^r)$ is bounded, any smooth bounded solution to the MHD or Hall--MHD system with local Dirichlet integral growing as an arbitrary power function must be constant. This hugely improves the result of \cite[Theorem 1.2]{pan2021Liouville}, where the Dirichlet integral of $\mathbf{u}$ is assumed to be finite. Motivated by \cite[Bang--Gui--Wang--Xie, 2022, {\it arXiv:2205.13259}]{bang2022Liouvilletype}, our proof relies on establishing Saint--Venant's estimates associated with our problem, and the result in the current paper extends that for stationary Navier--Stokes equations shown by \cite{bang2022Liouvilletype} to MHD and Hall--MHD equations. To achieve this, more intricate estimates are needed to handle the terms involving $\mathbf{b}$ properly.
△ Less
Submitted 12 December, 2022;
originally announced December 2022.
-
RVE Analysis in LS-DYNA for High-fidelity Multiscale Material Modeling
Authors:
Haoyan Wei,
Dandan Lyu,
Wei Hu,
C. T. Wu
Abstract:
In modern engineering designs, advanced materials (e.g., fiber/particle-reinforced polymers, metallic alloys, laminar composites, etc.) are widely used, where microscale heterogeneities such as grains, inclusions, voids, micro-cracks, and interfaces significantly affect the macroscopic constitutive behaviors. Obviously, an accurate description of the multiscale material behaviors is of great impor…
▽ More
In modern engineering designs, advanced materials (e.g., fiber/particle-reinforced polymers, metallic alloys, laminar composites, etc.) are widely used, where microscale heterogeneities such as grains, inclusions, voids, micro-cracks, and interfaces significantly affect the macroscopic constitutive behaviors. Obviously, an accurate description of the multiscale material behaviors is of great importance to the success of material design and structural analysis. The Representative Volume Element (RVE) analysis method provides a rigorous means to obtain homogenized macroscopic material properties at the upper length scale from the properties of the material constituents and structures at a lower length scale. Recently, we have developed an RVE module (keyword: *RVE_ANALYSIS_FEM) in the multiphysics simulation software LS-DYNA to enable high-fidelity virtual testing of numerically re-constructed material samples at user-specified characteristic length scales. In this article, a brief introduction to this new feature will be given.
△ Less
Submitted 21 October, 2022;
originally announced October 2022.
-
Boundary complexity and surface entropy of 2-multiplicative integer systems on $\mathbb{N}^d$
Authors:
Jung-Chao Ban,
Wen-Guei Hu,
Guan-Yu Lai
Abstract:
In this article, we introduce the concept of the boundary complexity and prove that for a 2-multiplicative integer system (2-MIS) $X^{p}_Ω$ on $\mathbb{N}$ (or $X^{\bf p}_Ω$ on $\mathbb{N}^d,d\geq 2$), every point in $[h(X^p_Ω), \log r]$ can be realized as a boundary complexity of a 2-MIS with a specific speed, where r stands for the number of the alphabets. The result is new and quite different f…
▽ More
In this article, we introduce the concept of the boundary complexity and prove that for a 2-multiplicative integer system (2-MIS) $X^{p}_Ω$ on $\mathbb{N}$ (or $X^{\bf p}_Ω$ on $\mathbb{N}^d,d\geq 2$), every point in $[h(X^p_Ω), \log r]$ can be realized as a boundary complexity of a 2-MIS with a specific speed, where r stands for the number of the alphabets. The result is new and quite different from $\mathbb{N}^d$ subshifts of finite type (SFT) for $d\geq 1$. Furthermore, the rigorous formula of surface entropy for a $\mathbb{N}^d$ 2-MIS is also presented. This provides an efficient method to calculate the topological entropy for $\mathbb{N}^d$ 2-MIS and also provides an intrinsic differences between $\mathbb{N}^d$ $k$-MIS and SFTs for $d\geq 1$ and $k\geq 2$.
△ Less
Submitted 17 October, 2022;
originally announced October 2022.
-
A cusp-capturing PINN for elliptic interface problems
Authors:
Yu-Hau Tseng,
Te-Sheng Lin,
Wei-Fan Hu,
Ming-Chih Lai
Abstract:
In this paper, we propose a cusp-capturing physics-informed neural network (PINN) to solve discontinuous-coefficient elliptic interface problems whose solution is continuous but has discontinuous first derivatives on the interface. To find such a solution using neural network representation, we introduce a cusp-enforced level set function as an additional feature input to the network to retain the…
▽ More
In this paper, we propose a cusp-capturing physics-informed neural network (PINN) to solve discontinuous-coefficient elliptic interface problems whose solution is continuous but has discontinuous first derivatives on the interface. To find such a solution using neural network representation, we introduce a cusp-enforced level set function as an additional feature input to the network to retain the inherent solution properties; that is, capturing the solution cusps (where the derivatives are discontinuous) sharply. In addition, the proposed neural network has the advantage of being mesh-free, so it can easily handle problems in irregular domains. We train the network using the physics-informed framework in which the loss function comprises the residual of the differential equation together with certain interface and boundary conditions. We conduct a series of numerical experiments to demonstrate the effectiveness of the cusp-capturing technique and the accuracy of the present network model. Numerical results show that even using a one-hidden-layer (shallow) network with a moderate number of neurons and sufficient training data points, the present network model can achieve prediction accuracy comparable with traditional methods. Besides, if the solution is discontinuous across the interface, we can simply incorporate an additional supervised learning task for solution jump approximation into the present network without much difficulty.
△ Less
Submitted 16 April, 2023; v1 submitted 15 October, 2022;
originally announced October 2022.
-
An efficient neural-network and finite-difference hybrid method for elliptic interface problems with applications
Authors:
Wei-Fan Hu,
Te-Sheng Lin,
Yu-Hau Tseng,
Ming-Chih Lai
Abstract:
A new and efficient neural-network and finite-difference hybrid method is developed for solving Poisson equation in a regular domain with jump discontinuities on embedded irregular interfaces. Since the solution has low regularity across the interface, when applying finite difference discretization to this problem, an additional treatment accounting for the jump discontinuities must be employed. H…
▽ More
A new and efficient neural-network and finite-difference hybrid method is developed for solving Poisson equation in a regular domain with jump discontinuities on embedded irregular interfaces. Since the solution has low regularity across the interface, when applying finite difference discretization to this problem, an additional treatment accounting for the jump discontinuities must be employed. Here, we aim to elevate such an extra effort to ease our implementation by machine learning methodology. The key idea is to decompose the solution into singular and regular parts. The neural network learning machinery incorporating the given jump conditions finds the singular solution, while the standard five-point Laplacian discretization is used to obtain the regular solution with associated boundary conditions. Regardless of the interface geometry, these two tasks only require supervised learning for function approximation and a fast direct solver for Poisson equation, making the hybrid method easy to implement and efficient. The two- and three-dimensional numerical results show that the present hybrid method preserves second-order accuracy for the solution and its derivatives, and it is comparable with the traditional immersed interface method in the literature. As an application, we solve the Stokes equations with singular forces to demonstrate the robustness of the present method.
△ Less
Submitted 2 March, 2023; v1 submitted 11 October, 2022;
originally announced October 2022.
-
Paradoxes and resolutions for semiparametric fusion of individual and summary data
Authors:
Wenjie Hu,
Ruoyu Wang,
Wei Li,
Wang Miao
Abstract:
Suppose we have available individual data from an internal study and various types of summary statistics from relevant external studies. External summary statistics have been used as constraints on the internal data distribution, which promised to improve the statistical inference in the internal data; however, the additional use of external summary data may lead to paradoxical results: efficiency…
▽ More
Suppose we have available individual data from an internal study and various types of summary statistics from relevant external studies. External summary statistics have been used as constraints on the internal data distribution, which promised to improve the statistical inference in the internal data; however, the additional use of external summary data may lead to paradoxical results: efficiency loss may occur if the uncertainty of summary statistics is not negligible and large estimation bias can emerge even if the bias of external summary statistics is small. We investigate these paradoxical results in a semiparametric framework. We establish the semiparametric efficiency bound for estimating a general functional of the internal data distribution, which is shown to be no larger than that using only internal data. We propose a data-fused efficient estimator that achieves this bound so that the efficiency paradox is resolved. Besides, a debiased estimator is further proposed which has selection consistency property by employing adaptive lasso penalty so that the resultant estimator can achieve the same asymptotic distribution as the oracle one that uses only unbiased summary statistics, which resolves the bias paradox. Simulations and application to a Helicobacter pylori infection dataset are used to illustrate the proposed methods.
△ Less
Submitted 17 July, 2023; v1 submitted 1 October, 2022;
originally announced October 2022.
-
Initial Value Problem Enhanced Sampling for Closed-Loop Optimal Control Design with Deep Neural Networks
Authors:
Xuanxi Zhang,
Jihao Long,
Wei Hu,
Weinan E,
Jiequn Han
Abstract:
Closed-loop optimal control design for high-dimensional nonlinear systems has been a long-standing challenge. Traditional methods, such as solving the associated Hamilton-Jacobi-Bellman equation, suffer from the curse of dimensionality. Recent literature proposed a new promising approach based on supervised learning, by leveraging powerful open-loop optimal control solvers to generate training dat…
▽ More
Closed-loop optimal control design for high-dimensional nonlinear systems has been a long-standing challenge. Traditional methods, such as solving the associated Hamilton-Jacobi-Bellman equation, suffer from the curse of dimensionality. Recent literature proposed a new promising approach based on supervised learning, by leveraging powerful open-loop optimal control solvers to generate training data and neural networks as efficient high-dimensional function approximators to fit the closed-loop optimal control. This approach successfully handles certain high-dimensional optimal control problems but still performs poorly on more challenging problems. One of the crucial reasons for the failure is the so-called distribution mismatch phenomenon brought by the controlled dynamics. In this paper, we investigate this phenomenon and propose the initial value problem enhanced sampling method to mitigate this problem. We theoretically prove that this sampling strategy improves over the vanilla strategy on the classical linear-quadratic regulator by a factor proportional to the total time duration. We further numerically demonstrate that the proposed sampling strategy significantly improves the performance on tested control problems, including the optimal landing problem of a quadrotor and the optimal reaching problem of a 7 DoF manipulator.
△ Less
Submitted 9 July, 2023; v1 submitted 8 September, 2022;
originally announced September 2022.
-
Rigidity degrees of indecomposable modules over representation-finite self-injective algebras
Authors:
Wei Hu,
Xiaojuan Yin
Abstract:
The rigidity degree of a generator-cogenerator determines the dominant dimension of its endomorphism algebra, and is closely related to a recently introduced homological dimension -- rigidity dimension. In this paper, we give explicit formulae for the rigidity degrees of all indecomposable modules over representation-finite self-injective algebras by developing combinatorial methods from the Eucli…
▽ More
The rigidity degree of a generator-cogenerator determines the dominant dimension of its endomorphism algebra, and is closely related to a recently introduced homological dimension -- rigidity dimension. In this paper, we give explicit formulae for the rigidity degrees of all indecomposable modules over representation-finite self-injective algebras by developing combinatorial methods from the Euclidean algorithm. As an application, the rigidity dimensions of some algebras of types $A$ and $E$ are given.
△ Less
Submitted 6 August, 2022;
originally announced August 2022.
-
On spatial entropy and periodic entropies of Two-dimensional Shifts of Finite Type
Authors:
Wen-Guei Hu,
Guan-Yu Lai,
Song-Sun Lin
Abstract:
Topological entropy or spatial entropy is a way to measure the complexity of shift spaces. This study investigates the relationships between the spatial entropy and the various periodic entropies which are computed by skew-coordinated systems $γ\in GL_2(\mathbb{Z})$ on two dimensional shifts of finite type.
Topological entropy or spatial entropy is a way to measure the complexity of shift spaces. This study investigates the relationships between the spatial entropy and the various periodic entropies which are computed by skew-coordinated systems $γ\in GL_2(\mathbb{Z})$ on two dimensional shifts of finite type.
△ Less
Submitted 22 July, 2022;
originally announced July 2022.
-
Solving optimal control of rigid-body dynamics with collisions using the hybrid minimum principle
Authors:
Wei Hu,
Jihao Long,
Yaohua Zang,
Weinan E,
Jiequn Han
Abstract:
Collisions are common in many dynamical systems with real applications. They can be formulated as hybrid dynamical systems with discontinuities automatically triggered when states transverse certain manifolds. We present an algorithm for the optimal control problem of such hybrid dynamical systems based on solving the equations derived from the hybrid minimum principle (HMP). The algorithm is an i…
▽ More
Collisions are common in many dynamical systems with real applications. They can be formulated as hybrid dynamical systems with discontinuities automatically triggered when states transverse certain manifolds. We present an algorithm for the optimal control problem of such hybrid dynamical systems based on solving the equations derived from the hybrid minimum principle (HMP). The algorithm is an iterative scheme following the spirit of the method of successive approximations (MSA), and it is robust to undesired collisions observed in the initial guesses. We analyze the discontinuities in the system and propose a stable collision condition, which is crucial for the convergence of iterative algorithms in systems experiencing collisions. Subsequently, we establish a convergence theorem demonstrating linear convergence for the MSA algorithm when collisions are present. We also address several numerical challenges introduced by the discontinuities. The algorithm is tested on disc collision problems whose optimal solutions exhibit one or multiple collisions. Linear convergence in terms of iteration steps and asymptotic first-order accuracy in terms of time discretization are observed when the algorithm is implemented with the forward-Euler scheme. The numerical results demonstrate that the proposed algorithm has better accuracy and convergence than direct methods based on gradient descent. Furthermore, the algorithm is also simpler, more accurate, and more stable than a deep reinforcement learning method.
△ Less
Submitted 10 May, 2023; v1 submitted 17 May, 2022;
originally announced May 2022.
-
Thermodynamic formalism and large deviation principle of multiplicative Ising models
Authors:
Jung-Chao Ban,
Wen-Guei Hu,
Guan-Yu Lai
Abstract:
The aim of this study is tree-fold. First, we investigate the thermodynamics of the Ising models with respect to 2-multiple Hamiltonians. This extends the previous results of [Chazotte and Redig, Electron. J. Probably., 2014] to $\mathbb{N}^d$. Second, we establish the large deviation principle (LDP) of the average $\frac{1}{N} S_N^G$, where $S_N^G$ is a 2-multiple sum along a semigroup generated…
▽ More
The aim of this study is tree-fold. First, we investigate the thermodynamics of the Ising models with respect to 2-multiple Hamiltonians. This extends the previous results of [Chazotte and Redig, Electron. J. Probably., 2014] to $\mathbb{N}^d$. Second, we establish the large deviation principle (LDP) of the average $\frac{1}{N} S_N^G$, where $S_N^G$ is a 2-multiple sum along a semigroup generated by k numbers which are k co-primes. This extends the previous results [Ban et al. Indag. Math., 2021] to a board class of the long-range interactions. Finally, the results described above are generalized to the multidimensional lattice $\mathbb{N}^d, d\geq1$.
△ Less
Submitted 16 March, 2022;
originally announced March 2022.
-
A Machine Learning Enhanced Algorithm for the Optimal Landing Problem
Authors:
Yaohua Zang,
Jihao Long,
Xuanxi Zhang,
Wei Hu,
Weinan E,
Jiequn Han
Abstract:
We propose a machine learning enhanced algorithm for solving the optimal landing problem. Using Pontryagin's minimum principle, we derive a two-point boundary value problem for the landing problem. The proposed algorithm uses deep learning to predict the optimal landing time and a space-marching technique to provide good initial guesses for the boundary value problem solver. The performance of the…
▽ More
We propose a machine learning enhanced algorithm for solving the optimal landing problem. Using Pontryagin's minimum principle, we derive a two-point boundary value problem for the landing problem. The proposed algorithm uses deep learning to predict the optimal landing time and a space-marching technique to provide good initial guesses for the boundary value problem solver. The performance of the proposed method is studied using the quadrotor example, a reasonably high dimensional and strongly nonlinear system. Drastic improvement in reliability and efficiency is observed.
△ Less
Submitted 13 March, 2022;
originally announced March 2022.
-
A shallow physics-informed neural network for solving partial differential equations on surfaces
Authors:
Wei-Fan Hu,
Yi-Jun Shih,
Te-Sheng Lin,
Ming-Chih Lai
Abstract:
In this paper, we introduce a shallow (one-hidden-layer) physics-informed neural network for solving partial differential equations on static and evolving surfaces. For the static surface case, with the aid of level set function, the surface normal and mean curvature used in the surface differential expressions can be computed easily. So instead of imposing the normal extension constraints used in…
▽ More
In this paper, we introduce a shallow (one-hidden-layer) physics-informed neural network for solving partial differential equations on static and evolving surfaces. For the static surface case, with the aid of level set function, the surface normal and mean curvature used in the surface differential expressions can be computed easily. So instead of imposing the normal extension constraints used in literature, we write the surface differential operators in the form of traditional Cartesian differential operators and use them in the loss function directly. We perform a series of performance study for the present methodology by solving Laplace-Beltrami equation and surface diffusion equation on complex static surfaces. With just a moderate number of neurons used in the hidden layer, we are able to attain satisfactory prediction results. Then we extend the present methodology to solve the advection-diffusion equation on an evolving surface with given velocity. To track the surface, we additionally introduce a prescribed hidden layer to enforce the topological structure of the surface and use the network to learn the homeomorphism between the surface and the prescribed topology. The proposed network structure is designed to track the surface and solve the equation simultaneously. Again, the numerical results show comparable accuracy as the static cases. As an application, we simulate the surfactant transport on the droplet surface under shear flow and obtain some physically plausible results.
△ Less
Submitted 20 January, 2023; v1 submitted 3 March, 2022;
originally announced March 2022.
-
On the Posterior Distribution of a Random Process Conditioned on Empirical Frequencies of a Finite Path: the i.i.d and finite Markov chain case
Authors:
Wenqing Hu,
Hong Qian
Abstract:
We obtain the posterior distribution of a random process conditioned on observing the empirical frequencies of a finite sample path. We find under a rather broad assumption on the "dependence structure" of the process, {\em c.f.} independence or Markovian, the posterior marginal distribution of the process at a given time index can be identified as certain empirical distribution computed from the…
▽ More
We obtain the posterior distribution of a random process conditioned on observing the empirical frequencies of a finite sample path. We find under a rather broad assumption on the "dependence structure" of the process, {\em c.f.} independence or Markovian, the posterior marginal distribution of the process at a given time index can be identified as certain empirical distribution computed from the observed empirical frequencies of the sample path. We show that in both cases of discrete-valued i.i.d. sequence and finite Markov chain, a certain "conditional symmetry" given by the observation of the empirical frequencies leads to the desired result on the posterior distribution. Results for both finite-time observations and its asymptotic infinite-time limit are connected via the idea of Gibbs conditioning. Finally, since our results demonstrate a central role of the empirical frequency in understanding the information content of data, we use the Large Deviations Principle (LDP) to construct a general notion of "data-driven entropy", from which one can apply a formalism from the recent study of statistical thermodynamics to data.
△ Less
Submitted 1 March, 2022; v1 submitted 23 February, 2022;
originally announced February 2022.
-
The equivalence of Friedlander-Mazur and standard conjectures for threefolds
Authors:
Jin Cao,
Wenchuan Hu
Abstract:
We show that the Friedlander-Mazur conjecture holds for a complex smooth projective variety X of dimension three implies the standard conjectures hold for X. This together with a result of Friedlander yields the equivalence of the two conjectures in dimension three. From this we provide some new examples whose standard conjectures hold.
We show that the Friedlander-Mazur conjecture holds for a complex smooth projective variety X of dimension three implies the standard conjectures hold for X. This together with a result of Friedlander yields the equivalence of the two conjectures in dimension three. From this we provide some new examples whose standard conjectures hold.
△ Less
Submitted 4 November, 2021;
originally announced November 2021.
-
Characterization and Topological Behavior of Homomorphism Tree-Shifts
Authors:
Jung-Chao Ban,
Chih-Hung Chang,
Wen-Guei Hu,
Guan-Yu Lai,
Yu-Liang Wu
Abstract:
The purpose of this article is twofold. On one hand, we reveal the equivalence of shift of finite type between a one-sided shift $X$ and its associated hom tree-shift $\mathcal{T}_{X}$, as well as the equivalence in the sofic shift. On the other hand, we investigate the interrelationship among the comparable mixing properties on tree-shifts as those on multidimensional shift spaces. They include i…
▽ More
The purpose of this article is twofold. On one hand, we reveal the equivalence of shift of finite type between a one-sided shift $X$ and its associated hom tree-shift $\mathcal{T}_{X}$, as well as the equivalence in the sofic shift. On the other hand, we investigate the interrelationship among the comparable mixing properties on tree-shifts as those on multidimensional shift spaces. They include irreducibility, topologically mixing, block gluing, and strong irreducibility, all of which are defined in the spirit of classical multidimensional shift, complete prefix code (CPC), and uniform CPC. In summary, the mixing properties defined in all three manners coincide for $\mathcal{T}_{X}$. Furthermore, an equivalence between irreducibility on $\mathcal{T}_{A}$ and irreducibility on $X_A$ are seen, and so is one between topologically mixing on $\mathcal{T}_{A}$ and mixing property on $X_A$, where $X_A$ is the one-sided shift space induced by the matrix $A$ and $T_A$ is the associated tree-shift. These equivalences are consistent with the mixing properties on $X$ or $X_A$ when viewed as a degenerate tree-shift.
△ Less
Submitted 30 August, 2021;
originally announced August 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.