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 ASYMPTOTIC BEHAVIOR AND OSCILLATION FOR A CLASS OF THIRD-ORDER NONLINEAR DELAY DYNAMIC EQUATIONS

In this paper,we consider a class of third-order nonlinear delay dynamic equations.First,we establish a Kiguradze-type lemma and some useful estimates.Second,we give a sufficient and necessary condition for the existence of eventually positive solutions having upper bounds and tending to zero.Third,we obtain new oscillation criteria by employing the Potzsche chain rule.Then,using the generalized Riccati transformation technique and averaging method,we establish the Philos-type oscillation criteria.Surprisingly,the integral value of the Philos-type oscillation criteria,which guarantees that all unbounded solutions oscillate,is greater than θ_(4)(t_(1),T).The results of Theorem 3.5 and Remark 3.6 are novel.Finally,we offer four examples to illustrate our results.

Dynamics of fundamental and double-pole breathers and solitons for a nonlinear Schrodinger equation with sextic operator under non-zero boundary conditions

We study the dynamics of fundamental and double-pole breathers and solitons for the focusing and defocusing nonlinear Schrodinger equation with the sextic operator under non-zero boundary conditions. Our analysis mainly focuses onthe dynamical properties of simple- and double-pole solutions. Firstly, through verification, we find that solutions undernon-zero boundary conditions can be transformed into solutions under zero boundary conditions, whether in simple-pole ordouble-pole cases. For the focusing case, in the investigation of simple-pole solutions, temporal periodic breather and thespatial-temporal periodic breather are obtained by modulating parameters. Additionally, in the case of multi-pole solitons,we analyze parallel-state solitons, bound-state solitons, and intersecting solitons, providing a brief analysis of their interactions.In the double-pole case, we observe that the two solitons undergo two interactions, resulting in a distinctive “triangle”crest. Furthermore, for the defocusing case, we briefly consider two situations of simple-pole solutions, obtaining one andtwo dark solitons.

Bifurcations, Analytical and Non-Analytical Traveling Wave Solutions of (2 + 1)-Dimensional Nonlinear Dispersive Boussinesq Equation

For the (2 + 1)-dimensional nonlinear dispersive Boussinesq equation, by using the bifurcation theory of planar dynamical systems to study its corresponding traveling wave system, the bifurcations and phase portraits of the regular system are obtained. Under different parametric conditions, various sufficient conditions to guarantee the existence of analytical and non-analytical solutions of the singular system are given by using singular traveling wave theory. For certain special cases, some explicit and exact parametric representations of traveling wave solutions are derived such as analytical periodic waves and non-analytical periodic cusp waves. Further, two-dimensional wave plots of analytical periodic solutions and non-analytical periodic cusp wave solutions are drawn to visualize the dynamics of the equation.

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.

Energy Stable Nodal DG Methods for Maxwell’s Equations of Mixed-Order Form in Nonlinear Optical Media

In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.

A Comparative Study of Adomian Decomposition Method with Variational Iteration Method for Solving Linear and Nonlinear Differential Equations

This study compares the Adomian Decomposition Method (ADM) and the Variational Iteration Method (VIM) for solving nonlinear differential equations in engineering. Differential equations are essential for modeling dynamic systems in various disciplines, including biological processes, heat transfer, and control systems. This study addresses first, second, and third-order nonlinear differential equations using Mathematica for data generation and graphing. The ADM, developed by George Adomian, uses Adomian polynomials to handle nonlinear terms, which can be computationally intensive. In contrast, VIM, developed by He, directly iterates the correction functional, providing a more straightforward and efficient approach. This study highlights VIM’s rapid convergence and effectiveness of VIM, particularly for nonlinear problems, where it simplifies calculations and offers direct solutions without polynomial derivation. The results demonstrate VIM’s superior efficiency and rapid convergence of VIM compared with ADM. The VIM’s minimal computational requirements make it practical for real-time applications and complex system modeling. Our findings align with those of previous research, confirming VIM’s efficiency of VIM in various engineering applications. This study emphasizes the importance of selecting appropriate methods based on specific problem requirements. While ADM is valuable for certain nonlinearities, VIM’s approach is ideal for many engineering scenarios. Future research should explore broader applications and hybrid methods to enhance the solution’s accuracy and efficiency. This comprehensive comparison provides valuable guidance for selecting effective numerical methods for differential equations in engineering.

Crank-Nicolson Quasi-Compact Scheme for the Nonlinear Two-Sided Spatial Fractional Advection-Diffusion Equations

The higher-order numerical scheme of nonlinear advection-diffusion equations is studied in this article, where the space fractional derivatives are evaluated by using weighted and shifted Grünwald difference operators and combining the compact technique, in the time direction is discretized by the Crank-Nicolson method. Through the energy method, the stability and convergence of the numerical scheme in the sense of L2-norm are proved, and the convergence order is . Some examples are given to show that our numerical scheme is effective.

Existence and Stability of Standing Waves for the Nonlinear Schrödinger Equation with Combined Nonlinearities and a Partial Harmonic Potential

In this paper, we study the existence of standing waves for the nonlinear Schrödinger equation with combined power-type nonlinearities and a partial harmonic potential. In the L2-supercritical case, we obtain the existence and stability of standing waves. Our results are complements to the results of Carles and Il’yasov’s artical, where orbital stability of standing waves have been studied for the 2D Schrödinger equation with combined nonlinearities and harmonic potential.

Stability of Standing Waves for the Nonlinear Schrödinger Equation with Mixed Power-Type and Hartree-Type Nonlinearities

This paper studies the existence of stable standing waves for the nonlinear Schrödinger equation with Hartree-type nonlinearity i∂tψ+Δψ+| ψ |pψ+(| x |−γ∗| ψ |2)ψ=0, (t,x)∈[ 0,T )×ℝN.Where ψ=ψ(t,x)is a complex valued function of (t,x)∈ℝ+×ℝN. The parameters N≥3, 0p4Nand 0γmin{ 4,N }. By using the variational methods and concentration compactness principle, we prove the orbital stability of standing waves.

Thermomechanical Dynamics (TMD) and Bifurcation-Integration Solutions in Nonlinear Differential Equations with Time-Dependent Coefficients

The new independent solutions of the nonlinear differential equation with time-dependent coefficients (NDE-TC) are discussed, for the first time, by employing experimental device called a drinking bird whose simple back-and-forth motion develops into water drinking motion. The solution to a drinking bird equation of motion manifests itself the transition from thermodynamic equilibrium to nonequilibrium irreversible states. The independent solution signifying a nonequilibrium thermal state seems to be constructed as if two independent bifurcation solutions are synthesized, and so, the solution is tentatively termed as the bifurcation-integration solution. The bifurcation-integration solution expresses the transition from mechanical and thermodynamic equilibrium to a nonequilibrium irreversible state, which is explicitly shown by the nonlinear differential equation with time-dependent coefficients (NDE-TC). The analysis established a new theoretical approach to nonequilibrium irreversible states, thermomechanical dynamics (TMD). The TMD method enables one to obtain thermodynamically consistent and time-dependent progresses of thermodynamic quantities, by employing the bifurcation-integration solutions of NDE-TC. We hope that the basic properties of bifurcation-integration solutions will be studied and investigated further in mathematics, physics, chemistry and nonlinear sciences in general.

GLOBAL SOLUTIONS IN THE CRITICAL SOBOLEV SPACE FOR THE LANDAU EQUATION

The Landau equation is studied for hard potential with-2≤γ≤1.Under a perturbation setting,a unique global solution of the Cauchy problem to the Landau equation is established in a critical Sobolev space H_(x)^(d)L_(v)^(2)(d>3/2),which extends the results of[11]in the torus domain to the whole space R_(x)^(3).Here we utilize the pseudo-differential calculus to derive our desired result.

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))∈R×H^(1)(R^(N))to the general Kirchhoff problem-M■,satisfying the normalization constraint f_(R)^N u^2dx=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.

Analysis of Extended Fisher-Kolmogorov Equation in 2D Utilizing the Generalized Finite Difference Method with Supplementary Nodes

In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

Implementation of the Integrated Green’s Function Method for 3D Poisson’s Equation in a Large Aspect Ratio Computational Domain

The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.

Derivation of a Revised Tsiolkovsky Rocket Equation That Predicts Combustion Oscillations

Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx1, and slow-jet, vx2. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.

High-Order Spatial FDTD Solver of Maxwell’s Equations for Terahertz Radiation Production

We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.

MULTIPLE POSITIVE SOLUTIONS TO A CLASS OF MODIFIED NONLINEAR SCHR?DINGER EQUATION IN A HIGH DIMENSION

The existence and multiplicity of positive solutions for equation(1.1)with the new critical exponent 4<p<2·2*shall be investigated in a high dimension.The conclusions extend the relative results recently attained in[1]for the one-dimensional case.More precisely,as the coefficient a(x)in the nonlinearity is sign-changing,the modified term 2(Δ(|u|2))u is still helpful for obtaining multiple positive solutions in a high dimension,even if a sign condition like∫_(R)N a(x)e_(1)^(p)dx<0(also named“a necessary condition”see[2,3])does not hold.

Coupled-generalized nonlinear Schr¨odinger equations solved by adaptive step-size methods in interaction picture

We extend two adaptive step-size methods for solving two-dimensional or multi-dimensional generalized nonlinear Schr ¨odinger equation(GNLSE): one is the conservation quantity error adaptive step-control method(RK4IP-CQE), and the other is the local error adaptive step-control method(RK4IP-LEM). The methods are developed in the vector form of fourthorder Runge–Kutta iterative scheme in the interaction picture by converting a vector equation in frequency domain. By simulating the supercontinuum generated from the high birefringence photonic crystal fiber, the calculation accuracies and the efficiencies of the two adaptive step-size methods are discussed. The simulation results show that the two methods have the same global average error, while RK4IP-LEM spends more time than RK4IP-CQE. The decrease of huge calculation time is due to the differences in the convergences of the relative photon number error and the approximated local error between these two adaptive step-size algorithms.