In this article, we develop a posteriori error analysis of a nonconforming finite element method for a linear quadratic elliptic distributed optimal control problem with two different sets of constraints, namely (i) integral state constraint and integral control constraint; (ii) integral state constraint and pointwise control constraints. In the analysis, we have taken the approach of reducing the state-control constrained minimization problem into a state minimization problem obtained by eliminating the control variable. The reliability and efficiency of a posteriori error estimator are discussed. Numerical results are reported to illustrate the behavior of the error estimator.
This paper proposes a finite element method for solving the periodic steady-state problem for the scalar-valued and vector-valued Poisson equations, a simple reduction model of the Maxwell equations under the Coulomb gauge. Introducing a new potential variable, we reformulate two systems composed of the scalar-valued and vector-valued Poisson problems to a single Hodge-Laplace problem for the 1-form in using the standard de Rham complex. Consequently, we can directly apply the finite-element exterior calculus (FEEC) theory in to deduce the well-posedness, stability, and convergence. Numerical examples using the cubical element are reported to validate the theoretical results.
Anomalous and non-ergodic diffusion is ubiquitous in the natural world. Fractional Feynman-Kac equations are used to characterize the functional distribution of the trajectories of the sub-diffusion particles which evolve slower than the Brownian motion. In this paper, by introducing the tempered fractional Jacobi functions (TFJFs), an efficient spectral collocation method is presented for solving the spatio-dependent temporal tempered fractional Feynman-Kac (STTFFK) equation. The tempered fractional differentiation matrix in the collocation scheme is generated by a recurrence relation which is fast and stable. The error estimate for the scheme is derived which shows that the proposed method is “spectral accuracy” if the solution is smooth. Several numerical experiments with different tempered factors are performed to support the theoretical results.
We establish a double logarithmic stability inequality for the problem of determining the initial data in an IBVP for the wave equation outside a non-trapping obstacle from two localized measurements.
In this paper, we consider solving a least-squares problem to the generalized Sylvester quaternion tensor equation. From the properties of quaternions and the Kronecker product, a tensor form of the LSQR algorithm is proposed for solving this problem, the convergence analysis of which is then established. Numerical results are reported to illustrate the feasibility and validity of the proposed algorithm compared with the tensor form of the CGLS method, including when the algorithm is tested with some randomly generated data and color video restoration problems.
Recent advances have been made in a wide range of imaging and pattern recognition applications, including picture categorization and object identification systems. These systems necessitate a robust feature extraction method. This study proposes a new class of orthogonal functions known as orthogonal mountain functions (OMFs). Using these functions, a novel set of orthogonal moments and associated scaling, rotation, and translation (SRT) invariants are presented for building a color image’s feature vector components. These orthogonal moments are presented as quaternion orthogonal mountain Fourier moments (QOMFMs). To demonstrate the validity of our theoretically recommended technique, we conduct a number of image analysis and pattern recognition experiments, including a comparison of the performance of the feature vectors proposed above to preexisting orthogonal invariant moments. The result of this study experimentally proves the effectiveness and quality of our QOMFMs.
A broad framework for time-fractional advection-diffusion equations is considered by incorporating a general class of Caputo-type fractional derivatives that includes many well-known fractional derivatives as particular cases. Within this general context, in this paper, we propose a numerical approach to provide approximate solutions to time-fractional advection-diffusion models that involve fractional derivatives of the generalized class. The developed approach is based on discretizing the studied models with respect to spatio-temporal domains using non-uniform meshes. Moreover, we discuss the stability of the proposed approach, which does not require solving large systems of linear equations. Numerical results and 3D graphics are presented for some time-fractional advection-diffusion models using several fractional derivatives. The feasibility and the reliability of the proposed approach are clearly demonstrated by comparing the numerical solutions of the studied models with the exact solutions in cases with a known solution.
This paper presents a novel approach to preserving invariants in the time-implicit numerical discretization of high-order nonlinear wave equations. Many high-order nonlinear wave equations have an infinite number of conserved quantities. Designing time-implicit numerical schemes that preserve many conserved quantities simultaneously is challenging. The proposed method utilizes Lagrange multipliers to reformulate the local discontinuous Galerkin (LDG) discretization as a conservative discretization. Combining an implicit spectral deferred correction method for time discretization can achieve a high-order scheme in both temporal and spatial dimensions. We use the Korteweg-de Vries equations as an example to illustrate the implementation of the algorithm. The algorithm can be easily generalized to other nonlinear wave problems with conserved quantities. Numerical examples for various high-order nonlinear wave equations demonstrate the effectiveness of the proposed methods in preserving invariants while maintaining the high accuracy.
We investigate the inverse scattering transform for the focusing nonlinear Schrödinger (NLS) equation with a particular class of nonvanishing boundary conditions (NVBCs), especially in the case of reflectionless potentials that give rise to a transmission coefficient with an arbitrary finite number N pairs of higher-order poles. The inverse problem is characterized in terms of a matrix Riemann-Hilbert (RH) problem equipped with several residue conditions at N pairs of higher-order poles. In the reflectionless case, we point out that the N-multipole soliton solutions including higher-order Kuznetsov-Ma breathers and Akhmediev breathers can be reconstructed by a linear algebraic system. Furthermore, we verify these special solutions by numerical simulations and display their density structures, also derive some Peregrine solitons by choosing appropriate parameters and taking limits.
Reaction-diffusion equations model various biological, physical, sociological, and environmental phenomena. Often, numerical simulations are used to understand and discover the dynamics of such systems. Following the extension of the nonlinear Lyapunov theory applied to some class of reaction-diffusion partial differential equations (PDEs), we develop the first fully discrete Lyapunov discretizations that are consistent with the stability properties of the continuous parabolic reaction-diffusion models. The proposed framework provides a systematic procedure to develop fully discrete schemes of arbitrary order in space and time for solving a broad class of equations equipped with a Lyapunov functional. The new schemes are applied to solve systems of PDEs, which arise in epidemiology and oncolytic M1 virotherapy. The new computational framework provides physically consistent and accurate results without exhibiting scheme-dependent instabilities and converging to unphysical solutions. The proposed approach represents a capstone for developing efficient, robust, and predictive technologies for simulating complex phenomena.
Very recently, Qi and Cui extended the Perron-Frobenius theory to dual number matrices with primitive and irreducible nonnegative standard parts and proved that they have a Perron eigenpair and a Perron-Frobenius eigenpair. The Collatz method was also extended to find the Perron eigenpair. Qi and Cui proposed two conjectures. One is the k-order power of a dual number matrix, which tends to zero if and only if the spectral radius of its standard part is less than one, and another is the linear convergence of the Collatz method. In this paper, we confirm these conjectures and provide the theoretical proof. The main contribution is to show that the Collatz method R-linearly converges with an explicit rate.
A stable and high-order accurate solver for linear and nonlinear parabolic equations is presented. An additive Runge-Kutta method is used for the time stepping, which integrates the linear stiff terms by an explicit singly diagonally implicit Runge-Kutta (ESDIRK) method and the nonlinear terms by an explicit Runge-Kutta (ERK) method. In each time step, the implicit solution is performed by the recently developed Hierarchical Poincaré-Steklov (HPS) method. This is a fast direct solver for elliptic equations that decomposes the space domain into a hierarchical tree of subdomains and builds spectral collocation solvers locally on the subdomains. These ideas are naturally combined in the presented method since the singly diagonal coefficient in ESDIRK and a fixed time step ensures that the coefficient matrix in the implicit solution of HPS remains the same for all time stages. This means that the precomputed inverse can be efficiently reused, leading to a scheme with complexity (in two dimensions) for the precomputation where the solution operator to the elliptic problems is built, and then for the solution in each time step. The stability of the method is proved for first order in time and any order in space, and numerical evidence substantiates a claim of the stability for a much broader class of time discretization methods. Numerical experiments supporting the accuracy of the efficiency of the method in one and two dimensions are presented.
This paper examines an optimal control problem for mean-field systems under partial observation. The state system is described by a controlled mean-field forward-backward stochastic differential equation that features correlated noises between the system and the observation. Furthermore, the observation coefficients are allowed to depend not only on the control process but also on its probability distribution. Assuming a convex control domain and allowing all coefficients of the systems to be random, we derive directly a variational formula for the cost function in a given control process direction in terms of the Hamiltonian and the associated adjoint system without relying on variational systems under standard assumptions on the coefficients. As an application, we present the necessary and sufficient conditions for the optimality of our control problem using Pontryagin’s maximum principle in a unified manner. The research provides insights into the optimization of mean-field systems under partial observation, which has practical implications for various applications.
Solutions of convection-dominated convection-diffusion problems usually possess layers, which are regions where the solution has a steep gradient. It is well known that many classical numerical discretization techniques face difficulties when approximating the solution to these problems. In recent years, physics-informed neural networks (PINNs) for approximating the solution to (initial-)boundary value problems ((I)BVPs) received a lot of interest. This paper studies various loss functionals for PINNs that are especially designed for convection-dominated convection-diffusion problems and that are novel in the context of PINNs. They are numerically compared to the vanilla and an hp-variational loss functional from the literature based on two steady-state benchmark problems whose solutions possess different types of layers. We observe that the best novel loss functionals reduce the error by 17.3% for the first and 5.5% for the second problem compared to the methods from the literature.
We consider a one-dimensional fluid-solid interaction model governed by the Burgers equation with a time varying interface. We discuss the inverse problem of determining the shape of the interface from Dirichlet and Neumann data at one endpoint of the spatial interval. In particular, we establish unique results and some conditional stability estimates. For the proofs, we use and adapt some lateral estimates that, in turn, rely on appropriate Carleman and interpolation inequalities.
It is important to compute the Hilbert transform of a given function defined on a finite interval. In 2013, Micchelli and his collaborators proposed a fast algorithm, which is called the Hilbert spline transform, to calculate the Hilbert transform of a given function on a finite interval with the computational complexity , where the spline knots were chosen to be the midpoints of sampling points. A natural question is that, whether or not the spline knots can be chosen to be the same as the sampling points. This paper gives a positive answer to this question. Besides, the analytic expression of the Hilbert transform of B-splines of any order is also established. Furthermore, the problem of how to choose spline coefficients, using the quasi-interpolation method or interpolation method, is also considered, although both make sure an optimal approximation order. Several interesting numerical examples are implemented and compared with most of the existing methods. Numerical results show that the proposed algorithm has a relatively high computational accuracy as well as a relatively low computational complexity.
In this paper, we propose a numerical approach for the fractional reaction sub-diffusion equation with a Caputo-Hadamard derivative of fractional order , whose solutions typically exhibit singular behavior at the initial time. The approach involves the use of an L2-type discrete operator of order to approximate the fractional derivative on nonuniform temporal meshes, while a standard second-order difference method is employed for uniform spatial meshes. Detailed discussions are presented on the truncation errors and coefficients of the proposed discrete fractional derivative operator. The stability as well as the accuracy of the resulting numerical scheme is rigorously analyzed on a unique nonuniform mesh. Several numerical examples are provided to validate the theoretical results and to demonstrate the efficiency of the proposed method.
We aim to propose a Riemann-problem-solver-free nonstaggered central difference scheme (NCDS), which is well-balanced and positivity-preserving for two-dimensional shallow water flows over nonflat bottom topography with wet-dry fronts. The well-balanced property of the NCDS for the two-dimensional shallow water equations at wet-dry fronts is still unclear. The NCDS is “generously multidimensional” and consists of three steps: a forward step, a corrector step, and a backward step. Each step needs to guarantee the water depth to be nonnegative and retain the lake at rest. The key ingredient of the NCDS is the backward step. It is very nontrivial and technical to obtain both the well-balanced and positivity-preserving properties when the computational domain contains wet-dry fronts for the NCDS. The traditional NCDS fails to retain the lake at rest when the computational domain contains wet-dry fronts. We introduce a novel backward step, which is well-balanced and positivity-preserving, to recover numerical solutions on nonstaggered cells. Finally, we show several numerical results of two-dimensional shallow water equations with wet-dry fronts to verify these properties.