Stability Analysis of Periodic Solutions of Some Duffing’s Equations

In this paper, some stability results were reviewed. A suitable and complete Lyapunov function for the hard spring model was constructed using the Cartwright method. This approach was compared with the existing results which confirmed a superior global stability result. Our contribution relies on its application to high damping door constructions. (2010 Mathematics Subject Classification: 34B15, 34C15, 34C25, 34K13.)

A Study of Periodic Solution of a Duffing’s Equation Using Implicit Function Theorem

In this paper, the well known implicit function theorem was applied to study existence and uniqueness of periodic solution of Duffing-type equation. Un-der appropriate conditions around the origin, a unique periodic solution was obtained.

The Extended Non-Elementary Amplitude Functions as Solutions to the Damped Pendulum Equation, the Van der Pol Equation, the Damped Duffing Equation, the Lienard Equation and the Lorenz Equations

In this paper, we define some non-elementary amplitude functions that are giving solutions to some well-known second-order nonlinear ODEs and the Lorenz equations, but not the chaos case. We are giving the solutions a name, a symbol and putting them into a group of functions and into the context of other functions. These solutions are equal to the amplitude, or upper limit of integration in a non-elementary integral that can be arbitrary. In order to define solutions to some short second-order nonlinear ODEs, we will make an extension to the general amplitude function. The only disadvantage is that the first derivative to these solutions contains an integral that disappear at the second derivation. We will also do a second extension: the two-integral amplitude function. With this extension we have the solution to a system of ODEs having a very strange behavior. Using the extended amplitude functions, we can define solutions to many short second-order nonlinear ODEs.

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.

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.

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.

基于Duffing振子的铁路轨道移频信号测量系统设计

铁路轨道作为移频信号的传输载体,其对于信号对象传输轨迹疏密度的辨别能力,决定了轨道体系对移频信号的测量准确性;为实现对信号对象传输频率的准确辨别,设计基于Duffing振子的铁路轨道移频信号测量系统;将电源模块输出的电量信号,按需分配至STM32F103微处理器、移频信号辨别模块与DSP测量单元之中,完善各级应用部件之间的连接关系,实现铁路轨道移频信号测量系统的硬件方案设计;确定移频信号的Duffing振子描述条件,针对其定义形式,构造移频信号的相空间,完成对铁路轨道移频信号的提取;求解铁路轨道的Lyapunov指数,并构建移频信号的Duffing方程,确定混沌测量参数的取值范围,实现对铁路轨道移频信号的测量,联合各级应用部件,完成基于Duffing振子的铁路轨道移频信号测量系统设计;实验结果表明,上述系统测量所得信号传输轨迹疏密度与真实轨迹疏密度相似度较高,并且测量结果准确率最高值达到了84%左右,说明其能够实现对铁路轨道移频信号的准确测量。

2π-periodic Solution for a Kind High Order Delay Duffing Equation

In this paper, the high order delay Duffing equation ax^（2n） ＋bx＋g（x（t-r）） = p（t） are considered,using the theory of coincidence degree, the sufficient condition for its there being at least a 2π-periodic solution is obtained.

Approximate Solutions of Primary Resonance for Forced Duffing Equation by Means of the Homotopy Analysis Method

Nonlinear dynamic equation is a common engineering model.There is not precise analytical solution for most of nonlinear differential equations.These nonlinear differential equations should be solved by using approximate methods.Classical perturbation methods such as LP method,KBM method,multi-scale method and the averaging method on weakly nonlinear vibration system is effective,while the strongly nonlinear system is difficult to apply.Approximate solutions of primary resonance for forced Duffing equation is investigated by means of homotopy analysis method （HAM）.Different from other approximate computational method,the HAM is totally independent of small physical parameters,and thus is suitable for most nonlinear problems.The HAM provides a great freedom to choose base functions of solution series,so that a nonlinear problem may be approximated more effectively.The HAM provides us a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter and the auxiliary function.Therefore,HAM not only may solve the weakly non-linear problems but also may be suitable for the strong non-linear problem.Through the approximate solution of forced Duffing equation with cubic non-linearity,the HAM and fourth order Runge-Kutta method of numerical solution were compared,the results show that the HAM not only can solve the steady state solution,but also can calculate the unsteady state solution,and has the good computational accuracy.

Exact Solutions of the Cubic Duffing Equation by Leaf Functions under Free Vibration

Exact solutions of the cubic Duffing equation with the initial conditions are presented.These exact solutions are expressed in terms of leaf functions and trigonometric functions.The leaf function r=sleafn(t)or r=cleafn(t)satisfies the ordinary differential equation dx2/dt2=-nr2n-1.The second-order differential of the leaf function is equal to-n times the function raised to the(2n-1)power of the leaf function.By using the leaf functions,the exact solutions of the cubic Duffing equation can be derived under several conditions.These solutions are constructed using the integral functions of leaf functions sleaf2(t)and cleaf2(t)for the phase of a trigonometric function.Since the leaf function and the trigonometric function are used in combination,a highly accurate solution of the Duffing equation can be easily obtained based on the data of leaf functions.In this study,seven types of the exact solutions are derived from leaf functions;the derivation of the seven exact solutions is detailed in the paper.Finally,waves obtained by the exact solutions are graphically visualized with the numerical results.

Damped and Divergence Exact Solutions for the Duffing Equation Using Leaf Functions and Hyperbolic Leaf Functions

According to the wave power rule,the second derivative of a functionχ(t)with respect to the variable t is equal to negative n times the functionχ(t)raised to the power of 2n?1.Solving the ordinary differential equations numerically results in waves appearing in the figures.The ordinary differential equation is very simple;however,waves,including the regular amplitude and period,are drawn in the figure.In this study,the function for obtaining the wave is called the leaf function.Based on the leaf function,the exact solutions for the undamped and unforced Duffing equations are presented.In the ordinary differential equation,in the positive region of the variableχ(t),the second derivative d^2χ(t)/dt^2 becomes negative.Therefore,in the case that the curves vary with the time under the conditionχ(t)>0,the gradient dχ(t)/d constantly decreases as time increases.That is,the tangential vector on the curve of the graph(with the abscissa and the ordinate χ(t)changes from the upper right direction to the lower right direction as time increases.On the other hand,in the negative region of the variableχ(t),the second derivative d^2χ(t)/dt^2 becomes positive.The gradient d χ(t)/d constantly increases as time decreases.That is,the tangent vector on the curve changes from the lower right direction to the upper right direction as time increases.Since the behavior occurring in the positive region of the variable χ(t)and the behavior occurring in the negative region of the variableχ(t)alternately occur in regular intervals,waves appear by these interactions.In this paper,I present seven types of damped and divergence exact solutions by combining trigonometric functions,hyperbolic functions,hyperbolic leaf functions,leaf functions,and exponential functions.In each type,I show the derivation method and numerical examples,as well as describe the features of the waveform.

A New Method for the Exact Solution of Duffing Equation

A lot of methods, such as Jacobian elliptic function analysis, are used to look for the explicit exact solution of Duffing differential equation. The key of the analysis is to construct quotient trigonometric function, and then nonlinear algebraic equation set theory and method are used for the solution of some kinds of nonlinear Duffing differential equation. In this paper, the exact solution of Duffing equation is obtained by using constant variation method, making use of the formula to solve cubic equations and general solution of the homogeneous equation of Duffing equation with appropriate Constant m and function f(t) .

Existence and Uniqueness of 2π-Periodic Solution About Duffing Equation and Its Numerical Method

A new provement of the existence and uniqueness about periodic boundary value Duffing equation is established by using global inverse function theorem. An algorithm for solving differential equation that has a large convergence domain is given. Finally, a numerical example is given.

1/3 PURE SUB-HARMONIC SOLUTION AND FRACTAL CHARACTERISTIC OF TRANSIENT PROCESS FOR DUFFING'S EQUATION

The 1/3 sub-harmonic solution for the Duffing＇s with damping equation was investigated by using the methods of harmonic balance and numerical integration. The assumed solution is introduced, and the domain of sub-harmonic frequencies was found. The asymptotical stability of the subharmonic resonances and the sensitivity of the amplitude responses to the variation of damping coefficient were examined. Then, the subharmonic resonances were analyzed by using the techniques from the general fractal theory. The analysis indicates that the sensitive dimensions of the system time-field responses show sensitivity to the conditions of changed initial perturbation, changed damping coefficient or the amplitude of excitation, thus the sensitive dimension can clearly describe the characteristic of the transient process of the subharmonic resonances.

A STUDY OF THE STATIC AND GLOBAL BIFURCATIONS FOR DUFFING EQUATION

In this paper, the static and global bifurcations of the forced Duffing equation have been studied by means of the averaged system. Bifurcation condition has been obtained in the whole parametric space. The change of the phase plane structure has been investigated.

A SIMPLIFIED PROOF TO A THEOREM BY DING TONG-REN ON PERIODIC SOLUTIONS OF DUFFING EQUATIONS

In 1982, DING Tong-ren gave a basic theorem about existence of periodic solutions of Duffing equations with double resonance. A simplified proof will be given by making use of the Leray-Schauder principle.

Finding Periodic Solutions of High Order Duffing Equations via Homotopy Method

This paper presents a detailed analysis of finding the periodic solutions for the high order Duffing equation x^（2n） ＋ g（x） = e（t） （n ≥ 1）. Firstly, we give a constructive proof for the existence of periodic solutions via the homotopy method. Then we establish an efficient and global convergence method to find periodic solutions numerically.

Harmonic Solutions of Duffing Equation with Singularity via Time Map

This paper is devoted to the study of second-order Duffing equation with singularity at the origin, where? tends to positive infinity as , and the primitive function as . By applying the phase-plane analysis methods and Poincaré-Bohl theorem, we obtain the existence of harmonic solutions of the given equation under a kind of nonresonance condition for the time map.