Maximum Interval of Stability and Convergence of Solution of a Forced Mathieu’s Equation

This paper investigates the maximum interval of stability and convergence of solution of a forced Mathieu’s equation, using a combination of Frobenius method and Eigenvalue approach. The results indicated that the equilibrium point was found to be unstable and maximum bounds were found on the derivative of the restoring force showing sharp condition for the existence of periodic solution. Furthermore, the solution to Mathieu’s equation converges which extends and improves some results in literature.

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.

Dynamic characteristics of resonant gyroscopes study based on the Mathieu equation approximate solution

Dynamic characteristics of the resonant gyroscope are studied based on the Mathieu equation approximate solution in this paper.The Mathieu equation is used to analyze the parametric resonant characteristics and the approximate output of the resonant gyroscope.The method of small parameter perturbation is used to analyze the approximate solution of the Mathieu equation.The theoretical analysis and the numerical simulations show that the approximate solution of the Mathieu equation is close to the dynamic output characteristics of the resonant gyroscope.The experimental analysis shows that the theoretical curve and the experimental data processing results coincide perfectly,which means that the approximate solution of the Mathieu equation can present the dynamic output characteristic of the resonant gyroscope.The theoretical approach and the experimental results of the Mathieu equation approximate solution are obtained,which provides a reference for the robust design of the resonant gyroscope.

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.

SOME EXTENDED RESULTS OF“SUBHARMONIC RESONANCE BIFURCATION THEORY OF NONLINEAR MATHIEU EQUATION”

The authors of [1] discussed the subharmonic resonance bifurcation theory of nonlinear Mathieu equation and obtained six bifurcation diagrams in -plane. In this paper, we extended the results of[1] and pointed out that there may exist as many as fourteen bifurcation diagrams which are not topologically equivalent to each other.

Review of Kelvin's Equation and Its Modification in Characterization of Mesoporous Materials

Physical and mathematical models as well as calculation methods of nitrogen bed on porous media have been introduced to evaluate the structural parameters of mesoporous materials. Kelvin＇s equation is a link between the relative adsorbate pressure, the mean pore radius, and pore capillarity on the basis of macroscopic capillary condensation. However, Kelvin＇s equation has been identified that it underestimates the calculated pore size of a material especially in the boundary of pore size which is between 2 and 4 nm.Various modifications on Kelvin＇s equation were mentioned in order to develop a new model to improve the accuracy of pore size calculation. The problems on conventional mathematical models were analyzed and discussed. A number of calculation methods on physisorption and pore size, especially fundamental theories of physisorption, basis of models and their deficiencies are reviewed. It can provide guidance on developing a modified Kelvin＇s equation for pore size calculation.

Efficient Cell-Centered Multigrid Scheme for the Three-Dimensional Navier-Stokes Equations

A cell centered scheme for three dimensional Navier Stokes equations, which is based on central difference approximations and Runge Kutta time stepping, is described. By using local time stepping, implicit residual smoothing, a multigrid method, and carefully controlled artificial dissipative terms, good convergence rates are obtained for two and three dimensional flows. The emphases are on the implicit smoothing and artificial dissipative terms with locally variable coefficients which depend on cel...

PROPER SPLITTINGS FOR RESTRICTED LINEAR EQUATIONS AND THE GENERALIZED INVERSE A_(T,S)~2

This paper presents a proper splitting iterative method for comparing the general restricted linear euqations Ax=b, x ∈T （where, b ∈AT, and T is an arbitrary but fixed subspace of Cm） and the generalized in AT,S For the special case when b ∈AT and dim（T）=dim（AT）, this splitting iterative methverse AT,S hod converges to AT,Sb （the unique solution of the general restricted system Ax=bx ∈T）.

Exact invariants and adiabatic invariants of Raitzin＇s canonical equations of motion for a nonlinear nonholonomic mechanical system

The exact invariants and the adiabatic invariants of Raitzin＇s canonical equations of motion for a nonlinear nonholonomic mechanical system are studied. The relations between the invariants and the symmetries of the system are established. Based on the concept of higher-order adiabatic invariant of a mechanical system under the action of a small perturbation, the forms of the exact invariants and adiabatic invariants and the conditions for their existence are proved. Finally, the inverse problem of the perturbation to symmetries of the system is studied and an example is also given to illustrate the application of the results.

EXISTENCE OF SOLUTION FOR BOUNDARY VALUE PROBLEM OF NONLINEAR FRACTIONAL DIFFERENTIAL EQUATION

This paper is concerned with the boundary value problem of a nonlinear fractional differential equation. By means of Schauder fixed-point theorem, an existence result of solution is obtained.

Generalized Extended tanh-function Method for Traveling Wave Solutions of Nonlinear Physical Equations

In this paper, the generalized extended tanh-function method is used for constructing the traveling wave solutions of nonlinear evolution equations. We choose Fisher＇s equation, the nonlinear schr＆#168;odinger equation to illustrate the validity and ad-vantages of the method. Many new and more general traveling wave solutions are obtained. Furthermore, this method can also be applied to other nonlinear equations in physics.

SOLVING COLORED YANG-BAXTER EQUATION BY WU’S METHOD

In this article, we discuss nonsymmetric solutions of the colored Yang-Baxter equation dependent on spectral as well as colored parameters and give all seven-vertex solutions by Wu＇s method. It is also proved that the solutions are composed of six groups of basic solutions up to five solution transformations. Moreover, al l solutions can be classified into two categories called Baxter type and free-fermion type.

Some Kind of Equations Involving Euler＇s Function

For any given positive integer n ≥ 1, the Euler function φ（n） is defined to be the number of positive integers not exceeding n which are relatively prime to n. w（n） is defined to be the number of different prime divisors of n. Some kind of equations involving Euler＇s function is studied in the paper.

Numerical Solution of Constrained Mechanical System Motions Equations and Inverse Problems of Dynamics

In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.

Computing Bifurcation Diagrams of Steady State Kuramoto Sivashinsky Equation by Difference Method

Utilizing difference formulae, we obtained the discrete systems of steady state Kuramoto Sivashinsky (K S) equation. Applied Newton's method and continuation technology to the systems, the bifurcated solutions are derived, and the bifurcation diagrams are constructed. All the results are successful and satisfactory.

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR NONLINEAR FRACTIONAL DIFFERENTIAL EQUATIONS WITH NON-SEPARATED TYPE INTEGRAL BOUNDARY CONDITIONS

In this paper, we study a boundary value problem of nonlinear fractional dif- ferential equations of order q （1 〈 q 〈 2） with non-separated integral boundary conditions. Some new existence and uniqueness results are obtained by using some standard fixed point theorems and Leray-Schauder degree theory. Some illustrative examples are also presented. We extend previous results even in the integer case q = 2.

A Family of Fifth-order Iterative Methods for Solving Nonlinear Equations

In this paper, we present and analyze a family of fifth-order iterative methods free from second derivative for solving nonlinear equations. It is established that the family of iterative methods has convergence order five. Numerical examples show that the new methods are comparable with the well known existing methods and give better results in many aspects.