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.

Fractional cyclic integrals and Routh equations of fractional Lagrange system with combined Caputo derivatives

In this paper, we develop a fractional cyclic integral and a Routh equation for fractional Lagrange system defined in terms of fractional Caputo derivatives. The fractional Hamilton principle and the fractional Lagrange equations of the system are obtained under a combined Caputo derivative. Furthermore, the fractional cyclic integrals based on the Lagrange equations are studied and the associated Routh equations of the system are presented. Finally, two examples are given to show the applications of the results.

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.

ROUTH’S EQUATIONS FOR GENERAL NONHOLONOMIC MECHANICAL SYSTEMS OF VARIABLE MASS

In this paper, Routh 's 'equations for the mechanical systems of the variable mass with nonlinear nonholonomic constraints of arbitrary orders in a noninertial reference system have been deduced not from any variational principles, but from the dynamical equations of Newtonian mechanics. And then again the other forms of equations for nonholonomic systems of variable mass are obtained from Routh's equations.

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.

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.

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.

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.

Noether Symmetry of Three-Order Lagrangian Equations

Based on the three-order Lagrangian equation, pseudo-Hamilton actoon I^＊ is defined and the three-order Hamilton＇s principle and the conditions are obtained in the paper. Then, the Noether symmetry about three-order Lagrangian equations is deduced. Finally, an example is given to illustrate the application of the result.

Fractional Pfaff-Birkhoff Principle and Birkhoff′s Equations in Terms of Riesz Fractional Derivatives

The dynamical and physical behavior of a complex system can be more accurately described by using the fractional model.With the successful use of fractional calculus in many areas of science and engineering,it is necessary to extend the classical theories and methods of analytical mechanics to the fractional dynamic system.Birkhoffian mechanics is a natural generalization of Hamiltonian mechanics,and its core is the Pfaff-Birkhoff principle and Birkhoff′s equations.The study on the Birkhoffian mechanics is an important developmental direction of modern analytical mechanics.Here,the fractional Pfaff-Birkhoff variational problem is presented and studied.The definitions of fractional derivatives,the formulae for integration by parts and some other preliminaries are firstly given.Secondly,the fractional Pfaff-Birkhoff principle and the fractional Birkhoff′s equations in terms of RieszRiemann-Liouville fractional derivatives and Riesz-Caputo fractional derivatives are presented respectively.Finally,an example is given to illustrate the application of the results.

GRADIENT ESTIMATES FOR SOLUTIONS TO QUASILINEAR ELLIPTIC EQUATIONS WITH CRITICAL SOBOLEV GROWTH AND HARDY POTENTIAL

This note is a continuation of the work[17].We study the following quasilinear elliptic equations- △pu-μ/｜x｜p ｜u｜p-2 u=Q（x）｜u｜Np/N-p -2u,x∈R N,where 1 〈 p 〈 N,0 ≤ μ 〈（（N-p）/p）p and Q ∈ L∞（RN）.Optimal asymptotic estimates on the gradient of solutions are obtained both at the origin and at the infinity.

Adaptive Maxwell’s Equations Derived Optimization and Its Application in Antenna Array Synthesis

In this paper,a self-adaptive method for the Maxwell’s Equations Derived Optimization(MEDO)is proposed.It is implemented by applying the Sequential Model-Based Optimization(SMBO)algorithm to the iterations of the MEDO,and achieves the automatic adjustment of the parameters.The proposed method is named as adaptive Maxwell’s equations derived optimization(AMEDO).In order to evaluate the performance of AMEDO,eight benchmarks are used and the results are compared with the original MEDO method.The results show that AMEDO can greatly reduce the workload of manual adjustment of parameters,and at the same time can keep the accuracy and stability.Moreover,the convergence of the optimization can be accelerated due to the dynamical adjustment of the parameters.In the end,the proposed AMEDO is applied to the side lobe level suppression and array failure correction of a linear antenna array,and shows great potential in antenna array synthesis.

Global Well-posedness of the Generalized Long-short Wave Equations

In the present paper, we investigate the well-posedness of the global solutionfor the Cauchy problem of generalized long-short wave equations. Applying Kato＇s methodfor abstract quasi-linear evolution equations and a priori estimates of solution,we get theexistence of globally smooth solution.

The Inertial Manifolds for a Class of Higher-Order Coupled Kirchhoff-Type Equations

In this paper, we mainly deal with a class of higher-order coupled Kirch-hoff-type equations. At first, we take advantage of Hadamard’s graph to get the equivalent form of the original equations. Then, the inertial manifolds are proved by using spectral gap condition. The main result we gained is that the inertial manifolds are established under the proper assumptions of M(s) and gi(u,v), i=1, 2.

A GPU-based general numerical framework for plasma simulations in terms of microscopic kinetic equations with full collision terms

We have proposed a general numerical framework for plasma simulations on graphics processing unit clusters based on microscopic kinetic equations with full collision terms.Our numerical algorithm consistently deals with both long-range(classical forces in the Vlasov term)and short-range(quantum processes in the collision term)interactions.Providing the relevant particle masses,charges and types(classical,fermionic or bosonic),as well as the external forces and the matrix elements(in the collisional integral),the algorithm consistently solves the coupled multi-particle kinetic equations.Currently,the framework is being tested and applied in the field of relativistic heavy-ion collisions;extensions to other plasma systems are straightforward.Our framework is a potential and competitive numerical platform for consistent plasma simulations.