A Structure-Preserving JKO Scheme for the SizeModified Poisson-Nernst-Planck-Cahn-Hilliard Equations

In this paper,we propose a structure-preserving numerical scheme for the size-modified Poisson-Nernst-Planck-Cahn-Hilliard(SPNPCH)equations derived from the free energy including electrostatic energies,entropies,steric energies,and Cahn-Hilliard mixtures.Based on the Jordan-Kinderlehrer-Otto(JKO)framework and the Benamou-Brenier formula of quadratic Wasserstein distance,the SPNPCH equations are transformed into a constrained optimization problem.By exploiting the convexity of the objective function,we can prove the existence and uniqueness of the numerical solution to the optimization problem.Mass conservation and unconditional energy-dissipation are preserved automatically by this scheme.Furthermore,by making use of the singularity of the entropy term which keeps the concentration from approaching zero,we can ensure the positivity of concentration.To solve the optimization problem,we apply the quasi-Newton method,which can ensure the positivity of concentration in the iterative process.Numerical tests are performed to confirm the anticipated accuracy and the desired physical properties of the developed scheme.Finally,the proposed scheme can also be applied to study the influence of ionic sizes and gradient energy coefficients on ion distribution.

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.

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.

Benchmark Computations of the Phase Field Crystal and Functionalized Cahn-Hilliard Equations via Fully Implicit,Nesterov Accelerated Schemes

We introduce a fast solver for the phase field crystal(PFC)and functionalized Cahn-Hilliard(FCH)equations with periodic boundary conditions on a rectangular domain that features the preconditioned Nesterov’s accelerated gradient descent(PAGD)method.We discretize these problems with a Fourier collocation method in space,and employ various second-order schemes in time.We observe a significant speedup with this solver when compared to the preconditioned gradient descent(PGD)method.With the PAGD solver,fully implicit,second-order-in-time schemes are not only feasible to solve the PFC and FCH equations,but also do so more efficiently than some semi-implicit schemes in some cases where accuracy issues are taken into account.Benchmark computations of four different schemes for the PFC and FCH equations are conducted and the results indicate that,for the FCH experiments,the fully implicit schemes(midpoint rule and BDF2 equipped with the PAGD as a nonlinear time marching solver)perform better than their IMEX versions in terms of computational cost needed to achieve a certain precision.For the PFC,the results are not as conclusive as in the FCH experiments,which,we believe,is due to the fact that the nonlinearity in the PFC is milder nature compared to the FCH equation.We also discuss some practical matters in applying the PAGD.We introduce an averaged Newton preconditioner and a sweeping-friction strategy as heuristic ways to choose good preconditioner parameters.The sweeping-friction strategy exhibits almost as good a performance as the case of the best manually tuned parameters.

Analytical solutions fractional order partial differential equations arising in fluid dynamics

This article describes the solution procedure of the fractional Pade-Ⅱ equation and generalized Zakharov equation(GSEs)using the sine-cosine method.Pade-Ⅱ is an important nonlinear wave equation modeling unidirectional propagation of long-wave in dispersive media and GSEs are used to model the interaction between one-dimensional high,and low-frequency waves.Classes of trigonometric and hyperbolic function solutions in fractional calculus are discussed.Graphical simulations of the numerical solutions are flaunted by MATLAB.

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.

Besov Estimates for Sub-Elliptic Equations in the Heisenberg Group

In this article, we deal with weak solutions to non-degenerate sub-elliptic equations in the Heisenberg group, and study the regularities of solutions. We establish horizontal Calderón-Zygmund type estimate in Besov spaces with more general assumptions on coefficients for both homogeneous equations and non-homogeneous equations. This study of regularity estimates expands the Calderón-Zygmund theory in the Heisenberg group.

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].

ELLIPTIC EQUATIONS IN DIVERGENCE FORM WITH DISCONTINUOUS COEFFICIENTS IN DOMAINS WITH CORNERS

We study equations in divergence form with piecewise Cαcoefficients.The domains contain corners and the discontinuity surfaces are attached to the edges of the corners.We obtain piecewise C1,αestimates across the discontinuity surfaces and provide an example to illustrate the issue regarding the regularity at the corners.

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.

On Two Types of Stability of Solutions to a Pair of Damped Coupled Nonlinear Evolution Equations

The stability of a set of spatially constant plane wave solutions to a pair of damped coupled nonlinear Schrödinger evolution equations is considered. The equations could model physical phenomena arising in fluid dynamics, fibre optics or electron plasmas. The main result is that any small perturbation to the solution remains small for all time. Here small is interpreted as being both in the supremum sense and the square integrable sense.

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.

THE GLOBAL EXISTENCE OF STRONG SOLUTIONS TO THERMOMECHANICAL CUCKER-SMALE-STOKES EQUATIONS IN THE WHOLE DOMAIN

We study the global existence and uniqueness of a strong solution to the kinetic thermomechanical Cucker-Smale(for short,TCS) model coupled with Stokes equations in the whole space.The coupled system consists of the kinetic TCS equation for a particle ensemble and the Stokes equations for a fluid via a drag force.In this paper,we present a complete analysis of the existence of global-in-time strong solutions to the coupled model without any smallness restrictions on the initial data.

GLOBAL UNIQUE SOLUTIONS FOR THE INCOMPRESSIBLE MHD EQUATIONS WITH VARIABLE DENSITY AND ELECTRICAL CONDUCTIVITY

We study the global unique solutions to the 2-D inhomogeneous incompressible MHD equations,with the initial data(u0,B0)being located in the critical Besov space■and the initial densityρ0 being close to a positive constant.By using weighted global estimates,maximal regularity estimates in the Lorentz space for the Stokes system,and the Lagrangian approach,we show that the 2-D MHD equations have a unique global solution.

NORMALIZED SOLUTIONS FOR THE GENERAL KIRCHHOFF TYPE EQUATIONS

In the present paper,we prove the existence,non-existence and multiplicity of positive normalized solutions(λ_(c),u_(c))∈ℝ×H^(1)(ℝ^(N))to the general Kirchhoff problem-M\left(\int_{\mathbb{R}^N}\vert\nabla u\vert^2{\rm d}x\right)\Delta u+\lambda u=g(u)~\hbox{in}~\mathbb{R}^N,u\in H^1(\mathbb{R}^N),N\geq 1,satisfying the normalization constraint\int_{\mathbb{R}^N}u^2{\rm d}x=c,where M∈C([0,∞))is a given function satisfying some suitable assumptions.Our argument is not by the classical variational method,but by a global branch approach developed by Jeanjean et al.[J Math Pures Appl,2024,183:44–75]and a direct correspondence,so we can handle in a unified way the nonlinearities g(s),which are either mass subcritical,mass critical or mass supercritical.

Multi-soliton solutions of coupled Lakshmanan-Porsezian-Daniel equations with variable coefficients under nonzero boundary conditions

This paper aims to investigate the multi-soliton solutions of the coupled Lakshmanan–Porsezian–Daniel equations with variable coefficients under nonzero boundary conditions.These equations are utilized to model the phenomenon of nonlinear waves propagating simultaneously in non-uniform optical fibers.By analyzing the Lax pair and the Riemann–Hilbert problem,we aim to provide a comprehensive understanding of the dynamics and interactions of solitons of this system.Furthermore,we study the impacts of group velocity dispersion or the fourth-order dispersion on soliton behaviors.Through appropriate parameter selections,we observe various nonlinear phenomena,including the disappearance of solitons after interaction and their transformation into breather-like solitons,as well as the propagation of breathers with variable periodicity and interactions between solitons with variable periodicities.

THE EXACT MEROMORPHIC SOLUTIONS OF SOME NONLINEAR DIFFERENTIAL EQUATIONS

We find the exact forms of meromorphic solutions of the nonlinear differential equations■,n≥3,k≥1,where q,Q are nonzero polynomials,Q■Const.,and p_(1),p_(2),α_(1),α_(2)are nonzero constants withα_(1)≠α_(2).Compared with previous results on the equation p(z)f^(3)+q(z)f"=-sinα(z)with polynomial coefficients,our results show that the coefficient of the term f^((k))perturbed by multiplying an exponential function will affect the structure of its solutions.

THE GLOBAL EXISTENCE AND ANALYTICITY OF A MILD SOLUTION TO THE 3D REGULARIZED MHD EQUATIONS

In this paper,we study the three-dimensional regularized MHD equations with fractional Laplacians in the dissipative and diffusive terms.We establish the global existence of mild solutions to this system with small initial data.In addition,we also obtain the Gevrey class regularity and the temporal decay rate of the solution.