HAMILTON-JACOBI EQUATIONS FOR A REGULAR CONTROLLED HAMILTONIAN SYSTEM AND ITS REDUCED SYSTEMS

In this paper,we give the geometric constraint conditions of a canonical symplectic form and regular reduced symplectic forms for the dynamical vector fields of a regular controlled Hamiltonian(RCH)system and its regular reduced systems,which are called the Type I and Type II Hamilton-Jacobi equations.First,we prove two types of Hamilton-Jacobi theorems for an RCH system on the cotangent bundle of a configuration manifold by using the canonical symplectic form and its dynamical vector field.Second,we generalize the above results for a regular reducible RCH system with symmetry and a momentum map,and derive precisely two types of Hamilton-Jacobi equations for the regular point reduced RCH system and the regular orbit reduced RCH system.Third,we prove that the RCH-equivalence for the RCH system,and the RpCH-equivalence and RoCH-equivalence for the regular reducible RCH systems with symmetries,leave the solutions of corresponding Hamilton-Jacobi equations invariant.Finally,as an application of the theoretical results,we show the Type I and Type II Hamilton-Jacobi equations for the Rp-reduced controlled rigid body-rotor system and the Rp-reduced controlled heavy top-rotor system on the generalizations of the rotation group SO(3)and the Euclidean group SE(3),respectively.This work reveals the deeply internal relationships of the geometrical structures of phase spaces,the dynamical vector fields and the controls of the RCH system.

THE LARGE TIME GENERIC FORM OF THE SOLUTION TO HAMILTON-JACOBI EQUATIONS

We use Hopf-Lax formula to study local regularity of solution to Hamilton- Jacobi （HJ） equations of multi-dimensional space variables with convex Hamiltonian. Then we give the large time generic form of the solution to HJ equation, i.e. for most initial data there exists a constant T 〉 0, which depends only on the Hamiltonian and initial datum, for t 〉 T the solution of the IVP （1.1） is smooth except for ~ smooth n-dimensional hypersurface, across which Du（x, t） is discontinuous. And we show that the hypersurface 1 tends asymptotically to a given hypersurface with rate t-1/4.

A NOTE ON GRADIENT BLOWUP RATE OF THE INHOMOGENEOUS HAMILTON-JACOBI EQUATIONS

The gradient blowup of the equation ut = △u ＋ a（x）｜△u｜p ＋ h（x）, where p 〉 2, is studied. It is shown that the gradient blowup rate will never match that of the self-similar variables. The exact blowup rate for radial solutions is established under the assumptions on the initial data so that the solution is monotonically increasing in time.

Arc Length-Based WENO Scheme for Hamilton-Jacobi Equations

In this article,novel smoothness indicators are presented for calculating the nonlinear weights of the weighted essentially non-oscillatory scheme to approximate the viscosity numerical solutions of Hamilton-Jacobi equations.These novel smoothness indicators are constructed from the derivatives of reconstructed polynomials over each sub-stencil.The constructed smoothness indicators measure the arc-length of the reconstructed polynomials so that the new nonlinear weights could get less absolute truncation error and give a high-resolution numerical solution.Extensive numerical tests are conducted and presented to show the performance capability and the numerical accuracy of the proposed scheme with the comparison to the classical WENO scheme.

UNIQUENESS OF VISCOSITY SOLUTIONS OF STOCHASTIC HAMILTON-JACOBI EQUATIONS

This article is devoted to the study of fully nonlinear stochastic Hamilton-Jacobi (HJ) equations for the optimal stochastic control problem of ordinary differential equations with random coefficients. Under the standard Lipschitz continuity assumptions on the coefficients, the value function is proved to be the unique viscosity solution of the associated stochastic HJ equation.

HIGH ORDER FINITE DIFFERENCE HERMITE WENO FAST SWEEPING METHODS FOR STATIC HAMILTON-JACOBI EQUATIONS

In this paper,we propose a novel Hermite weighted essentially non-oscillatory(HWENO)fast sweeping method to solve the static Hamilton-Jacobi equations efficiently.During the HWENO reconstruction procedure,the proposed method is built upon a new finite difference fifth order HWENO scheme involving one big stencil and two small stencils.However,one major novelty and difference from the traditional HWENO framework lies in the fact that,we do not need to introduce and solve any additional equations to update the derivatives of the unknown functionϕ.Instead,we use the currentϕand the old spatial derivative ofϕto update them.The traditional HWENO fast sweeping method is also introduced in this paper for comparison,where additional equations governing the spatial derivatives ofϕare introduced.The novel HWENO fast sweeping methods are shown to yield great savings in computational time,which improves the computational efficiency of the traditional HWENO scheme.In addition,a hybrid strategy is also introduced to further reduce computational costs.Extensive numerical experiments are provided to validate the accuracy and efficiency of the proposed approaches.

Symmetric Reduction and Hamilton-Jacobi Equations of the Controlled Underwater Vehicle-Rotor System

.As an application of the theoretical results,in this paper,we study the symmetric reduction and Hamilton-Jacobi theory for the underwater ve-hicle with two internal rotors as a regular point reducible RCH system,in the cases of coincident and non-coincident centers of the buoyancy and the gravity.At first,we give the regular point reduction and the two types of Hamilton-Jacobi equations for a regular controlled Hamiltonian(RCH)system with sym-metry and a momentum map on the generalization of a semidirect product Lie group.Next,we derive precisely the geometric constraint conditions of the reduced symplectic forms for the dynamical vector fields of the regular point reducible controlled underwater vehicle-rotor system,that is,the two types of Hamilton-Jacobi equations for the reduced controlled underwater vehicle-rotor system,by calculations in detail.These work reveal the deeply internal relationships of the geometrical structures of the phase spaces,the dynamical vector fields and the controls of the system.

Sparse-Grid Implementation of Fixed-Point Fast Sweeping WENO Schemes for Eikonal Equations

Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of fast sweeping schemes,fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order.The resulting iterative schemes have a fast convergence rate to steady-state solutions.Moreover,an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve the inverse operation of any nonlinear local system.Hence,they are robust and flexible,and have been combined with high-order accurate weighted essentially non-oscillatory(WENO)schemes to solve various hyperbolic PDEs in the literature.For multidimensional nonlinear problems,high-order fixed-point fast sweeping WENO methods still require quite a large amount of computational costs.In this technical note,we apply sparse-grid techniques,an effective approximation tool for multidimensional problems,to fixed-point fast sweeping WENO methods for reducing their computational costs.Here,we focus on fixed-point fast sweeping WENO schemes with third-order accuracy(Zhang et al.2006[41]),for solving Eikonal equations,an important class of static Hamilton-Jacobi(H-J)equations.Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse-grid computations of the fixed-point fast sweeping WENO schemes achieve large savings of CPU times on refined meshes,and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.

Some Modified Equations of the Sine-Hilbert Type

Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based on these bilinear equations, some exact solutions to the three modified equations are derived.

Data-Driven Ai-and Bi-Soliton of the Cylindrical Korteweg-de Vries Equation via Prior-Information Physics-Informed Neural Networks

By the modifying loss function MSE and training area of physics-informed neural networks(PINNs),we propose a neural networks model,namely prior-information PINNs(PIPINNs).We demonstrate the advantages of PIPINNs by simulating Ai-and Bi-soliton solutions of the cylindrical Korteweg-de Vries(cKdV)equation.

Matrix Riccati Equations in Optimal Control

In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied thoroughly, matrix Riccati equation of which scalar Riccati equations is a particular case, is much less investigated. This article proposes a change of variable that allows to find explicit solution of the Matrix Riccati equation. We then apply this solution to Optimal Control.

THE STABILITY OF BOUSSINESQ EQUATIONS WITH PARTIAL DISSIPATION AROUND THE HYDROSTATIC BALANCE

This paper is devoted to understanding the stability of perturbations around the hydrostatic equilibrium of the Boussinesq system in order to gain insight into certain atmospheric and oceanographic phenomena.The Boussinesq system focused on here is anisotropic,and involves only horizontal dissipation and thermal damping.In the 2D case R^(2),due to the lack of vertical dissipation,the stability and large-time behavior problems have remained open in a Sobolev setting.For the spatial domain T×R,this paper solves the stability problem and gives the precise large-time behavior of the perturbation.By decomposing the velocity u and temperatureθinto the horizontal average(ū,θ)and the corresponding oscillation(ū,θ),we can derive the global stability in H~2 and the exponential decay of(ū,θ)to zero in H^(1).Moreover,we also obtain that(ū_(2),θ)decays exponentially to zero in H^(1),and thatū_(1)decays exponentially toū_(1)(∞)in H^(1)as well;this reflects a strongly stratified phenomenon of buoyancy-driven fluids.In addition,we establish the global stability in H^(3)for the 3D case R^(3).

The Maxwell-Heaviside Equations Explained by the Theory of Informatons

In the articles “Newtons Law of Universal Gravitation Explained by the Theory of Informatons” and “The Gravitational Interaction between Moving Mass Particles Explained by the Theory of Informatons” the gravitational interaction has been explained by the hypothesis that information carried by informatons is the substance of gravitational fields, i.e. the medium that the interaction in question makes possible. From the idea that “information carried by informatons” is its substance, it has been deduced that—on the macroscopic level—a gravitational field manifests itself as a dual entity, always having a field- and an induction component (Egand Bg) simultaneously created by their common sources. In this article we will mathematically deduce the Maxwell-Heaviside equations from the kinematics of the informatons. These relations describe on the macroscopic level how a gravitational field (Eg, Bg) is generated by whether or not moving masses and how spatial and temporal changes of Egand Bgare related. We show that there is no causal link between Egand Bg.

THE SMOOTHING EFFECT IN SHARP GEVREY SPACE FOR THE SPATIALLY HOMOGENEOUS NON-CUTOFF BOLTZMANN EQUATIONS WITH A HARDPOTENTIAL

In this article, we study the smoothing effect of the Cauchy problem for the spatially homogeneous non-cutoff Boltzmann equation for hard potentials. It has long been suspected that the non-cutoff Boltzmann equation enjoys similar regularity properties as to whose of the fractional heat equation. We prove that any solution with mild regularity will become smooth in Gevrey class at positive time, with a sharp Gevrey index, depending on the angular singularity. Our proof relies on the elementary L^(2) weighted estimates.

On entire solutions of some Fermat type differential-difference equations

On one hand,we study the existence of transcendental entire solutions with finite order of the Fermat type difference equations.On the other hand,we also investigate the existence and growth of solutions of nonlinear differential-difference equations.These results extend and improve some previous in[5,14].

Equivalence between the internal observability and equation

This paper is concerned with a third order in time linear Moore-Gibson-Thompson equation which describes the acoustic velocity potential in ultrasound wave program.Influenced by the work of Kaltenbacher,Lasiecka and Marchand(Control Cybernet.2011,40:971-988),we establish an observability inequality of the conservative problem,and then discuss the equivalence between the exponential stabilization of a dissipative system and the internal observational inequality of the corresponding conservative system.

Multi-soliton solutions, breather-like and bound-state solitons for complex modified Korteweg–de Vries equation in optical fibers

Under investigation in this paper is a complex modified Korteweg–de Vries(KdV) equation, which describes the propagation of short pulses in optical fibers. Bilinear forms and multi-soliton solutions are obtained through the Hirota method and symbolic computation. Breather-like and bound-state solitons are constructed in which the signs of the imaginary parts of the complex wave numbers and the initial separations of the two parallel solitons are important factors for the interaction patterns. The periodic structures and position-induced phase shift of some solutions are introduced.

ON THE STABILITY OF PERIODIC SOLUTIONS OF PIECEWISE SMOOTH PERIODIC DIFFERENTIAL EQUATIONS

In this paper,we address the stability of periodic solutions of piecewise smooth periodic differential equations.By studying the Poincarémap,we give a sufficient condition to judge the stability of a periodic solution.We also present examples of some applications.

Exact and heuristic formulae to compute the geodetic height from the ellipse equation

The conversion of the cartesian coordinates of a point to its geodetic equivalent coordinates in reference to the geodetic ellipsoid is one of the main challenges in geodesy.The ellipse equation in the meridian plane significantly influences the value of the geodetic coordinates.This research analyzes this influence and how it can contribute to their solutions.The study investigates the mathematical relation between them and presents an exact formula relating to the geodetic height and the ellipse equation.In addition,a heuristic formula for the relation between the geodetic height and the ellipse equation is proposed,which is independent of the geodetic latitude and has a relative accuracy better than 99.9 %.The calculation is stable,and the cost is low.

A deep learning method based on prior knowledge with dual training for solving FPK equation

The evolution of the probability density function of a stochastic dynamical system over time can be described by a Fokker–Planck–Kolmogorov(FPK) equation, the solution of which determines the distribution of macroscopic variables in the stochastic dynamic system. Traditional methods for solving these equations often struggle with computational efficiency and scalability, particularly in high-dimensional contexts. To address these challenges, this paper proposes a novel deep learning method based on prior knowledge with dual training to solve the stationary FPK equations. Initially, the neural network is pre-trained through the prior knowledge obtained by Monte Carlo simulation(MCS). Subsequently, the second training phase incorporates the FPK differential operator into the loss function, while a supervisory term consisting of local maximum points is specifically included to mitigate the generation of zero solutions. This dual-training strategy not only expedites convergence but also enhances computational efficiency, making the method well-suited for high-dimensional systems. Numerical examples, including two different two-dimensional(2D), six-dimensional(6D), and eight-dimensional(8D) systems, are conducted to assess the efficacy of the proposed method. The results demonstrate robust performance in terms of both computational speed and accuracy for solving FPK equations in the first three systems. While the method is also applicable to high-dimensional systems, such as 8D, it should be noted that computational efficiency may be marginally compromised due to data volume constraints.