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.

Finite Element Approach for the Solution of First-Order Differential Equations

The finite element method has established itself as an efficient numerical procedure for the solution of arbitrary-shaped field problems in space. Basically, the finite element method transforms the underlying differential equation into a system of algebraic equations by application of the method of weighted residuals in conjunction with a finite element ansatz. However, this procedure is restricted to even-ordered differential equations and leads to symmetric system matrices as a key property of the finite element method. This paper aims in a generalization of the finite element method towards the solution of first-order differential equations. This is achieved by an approach which replaces the first-order derivative by fractional powers of operators making use of the square root of a Sturm-Liouville operator. The resulting procedure incorporates a finite element formulation and leads to a symmetric but dense system matrix. Finally, the scheme is applied to the barometric equation where the results are compared with the analytical solution and other numerical approaches. It turns out that the resulting numerical scheme shows excellent convergence properties.

Comparative Studies between Picard’s and Taylor’s Methods of Numerical Solutions of First Ordinary Order Differential Equations Arising from Real-Life Problems

To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.

An Instability Result to a Certain Vector Differential Equation of the Sixth Order

The nonlinear vector differential equation of the sixth order with constant delay is considered in this article. New criteria for instability of the zero solution are established using the Lyapunov-Krasovskii functional approach and the differential inequality techniques. The result of this article improves previously known results.

Extraction of P-and S-wave angle-domain common-image gathers based on first-order velocity-dilatation-rotation equations

Accuracy of angle-domain common-image gathers(ADCIGs)is the key to multiwave AVA inversion and migration velocity analysis,and of which Poynting vectors of pure P-and S-wave are the decisive factors in obtaining multi-component seismic data ADCIGs.A Poynting vector can be obtained from conventional velocity-stress elastic wave equations,but it focused on the propagation direction of mixed P-and S-wave fields,and neither on the propagation direction of the P-wave nor the direction of the S-wave.The Poynting vectors of pure P-or pure S-wave can be calculated from first-order velocity-dilatation-rotation equations.This study presents a method of extracting ADCIGs based on first order velocitydilatation-rotation elastic wave equations reverse-time migration algorithm.The method is as follows:calculating the pure P-wave Poynting vector of source and receiver wavefields by multiplication of P-wave particle-velocity vector and dilatation scalar,calculating the pure S-wave Poynting vector by vector multiplying S-wave particle-velocity vector and rotation vector,selecting the Poynting vector at the time of maximum P-wave energy of source wavefield as the propagation direction of incident P-wave,and obtaining the reflected P-wave(or converted S-wave)propagation direction of the receiver wavefield by the Poynting vector at the time of maximum P-(S-)wave energy in each grid point.Then,the P-wave incident angle is computed by the two propagation directions.Thus,the P-and S-wave ADGICs can obtained Numerical tests show that the proposed method can accurately compute the propagation direction and incident angle of the source and receiver wavefields,thereby achieving high-precision extraction of P-and S-wave ADGICs.

The Exact Solutions of Such Coupled Linear Matrix Fractional Differential Equations of Diagonal Unknown Matrices by Using Hadamard Product

In this paper, we present the general exact solutions of such coupled system of matrix fractional differential equations for diagonal unknown matrices in Caputo sense by using vector extraction operators and Hadamard product. Some illustrated examples are also given to show our new approach.

Linear Quadratic Optimal Control for Systems Governed by First-Order Hyperbolic Partial Differential Equations

This paper focuses on linear-quadratic(LQ)optimal control for a class of systems governed by first-order hyperbolic partial differential equations(PDEs).Different from most of the previous works,an approach of discretization-then-continuousization is proposed in this paper to cope with the infinite-dimensional nature of PDE systems.The contributions of this paper consist of the following aspects:(1)The differential Riccati equations and the solvability condition of the LQ optimal control problems are obtained via the discretization-then-continuousization method.(2)A numerical calculation way of the differential Riccati equations and a practical design way of the optimal controller are proposed.Meanwhile,the relationship between the optimal costate and the optimal state is established by solving a set of forward and backward partial difference equations(FBPDEs).(3)The correctness of the method used in this paper is verified by a complementary continuous method and the comparative analysis with the existing operator results is presented.It is shown that the proposed results not only contain the classic results of the standard LQ control problem of systems governed by ordinary differential equations as a special case,but also support the existing operator results and give a more convenient form of computation.

Multi-sources distinguishing of vector transducer via differential evolution

Vector transducer can simultaneously measure components of particle velocity as well as pressure at some point in sound field. In this paper, a series of equations are obtained from the correlation of particle velocity and pressure of the incident wave field, the error of each equation with white noise is studied, and Differential Evolution is used in solving the equations to distinguish multi-sources. Results of computer simulation show that Differential Evolution has more superiority than Genetic Algorithms on the rate and precision of convergence under the same condition.

An Alternative Method for Solving Lagrange＇s First-Order Partial Differential Equation with Linear Function Coefficients

An alternative method of solving Lagrange＇s first-order partial differential equation of the form（a1x ＋b1y＋C1z）p＋ （a2x ＋b2y＋c2z）q =a3x ＋b3y＋c3z,where p = Эz/Эx, q = Эz/Эy and ai, bi, ci （i = 1,2,3） are all real numbers has been presented here.

EXPONENTIAL FOURIER COLLOCATION METHODS FOR SOLVING FIRST-ORDER DIFFERENTIAL EQUATIONS

In this paper, a novel class of exponential Fourier collocation methods （EFCMs） is presented for solving systems of first-order ordinary differential equations. These so-called exponential Fourier collocation methods are based on the variation-of-constants formula, incorporating a local Fourier expansion of the underlying problem with collocation meth- ods. We discuss in detail the connections of EFCMs with trigonometric Fourier colloca- tion methods （TFCMs）, the well-known Hamiltonian Boundary Value Methods （HBVMs）, Gauss methods and Radau IIA methods. It turns out that the novel EFCMs are an es- sential extension of these existing methods. We also analyse the accuracy in preserving the quadratic invariants and the Hamiltonian energy when the underlying system is a Hamiltonian system. Other properties of EFCMs including the order of approximations and the convergence of fixed-point iterations are investigated as well. The analysis given in this paper proves further that EFCMs can achieve arbitrarily high order in a routine manner which allows us to construct higher-order methods for solving systems of first- order ordinary differential equations conveniently. We also derive a practical fourth-order EFCM denoted by EFCM（2,2） as an illustrative example. The numerical experiments using EFCM（2,2） are implemented in comparison with an existing fourth-order HBVM, an energy-preserving collocation method and a fourth-order exponential integrator in the literature. The numerical results demonstrate the remarkable efficiency and robustness of the novel EFCM（2,2）.

Contribution to the Analytical Equation Resolution Using Charts for Analysis and Design of Cylindrical and Conical Open Surge Tanks

In the event of an instantaneous valve closure, the pressure transmitted to a surge tank induces the mass fluctuations that can cause high amplitude of water-level fluctuation in the surge tank for a reasonable cross-sectional area. The height of the surge tank is then designed using this high water level mark generated by the completely closed penstock valve. Using a conical surge tank with a non-constant cross-sectional area can resolve the problems of space and height. When addressing issues in designing open surge tanks, key parameters are usually calculated by using complex equations, which may become cumbersome when multiple iterations are required. A more effective alternative in obtaining these values is the use of simple charts. Firstly, this paper presents and describes the equations used to design open conical surge tanks. Secondly, it introduces user-friendly charts that can be used in the design of cylindrical and conical open surge tanks. The contribution can be a benefit for practicing engineers in this field. A case study is also presented to illustrate the use of these design charts. The case study’s results show that key parameters obtained via successive approximation method required 26 iterations or complex calculations, whereas these values can be obtained by simple reading of the proposed chart. The use of charts to help surge tanks designing, in the case of preliminary designs, can save time and increase design efficiency, while reducing calculation errors.

Numeric Solution of the Fokker-Planck-Kolmogorov Equation

The solution of an n-dimensional stochastic differential equation driven by Gaussian white noises is a Markov vector. In this way, the transition joint probability density function (JPDF) of this vector is given by a deterministic parabolic partial differential equation, the so-called Fokker-Planck-Kolmogorov (FPK) equation. There exist few exact solutions of this equation so that the analyst must resort to approximate or numerical procedures. The finite element method (FE) is among the latter, and is reviewed in this paper. Suitable computer codes are written for the two fundamental versions of the FE method, the Bubnov-Galerkin and the Petrov-Galerkin method. In order to reduce the computational effort, which is to reduce the number of nodal points, the following refinements to the method are proposed: 1) exponential (Gaussian) weighting functions different from the shape functions are tested;2) quadratic and cubic splines are used to interpolate the nodal values that are known in a limited number of points. In the applications, the transient state is studied for first order systems only, while for second order systems, the steady-state JPDF is determined, and it is compared with exact solutions or with simulative solutions: a very good agreement is found.

Differential Equations and Singular Vectors in Verma Modules over sl(n,C)

Xu introduced a system of partial differential equations to investigate singular vectors in the Verma module of highest weight λ, over sl（n,C）. He gave a differential-operator representation of the symmetric group Sn on the corresponding space of truncated power series and proved that the solution space of the system is spanned by {σ（1）｜ σ ∈ Sn}. It is known that Sn is also the Weyl group of sl（n, C） and generated by all reflections sα with positive roots α. We present an explicit formula of the solution sα（1） for every positive root α and show directly that sα（1） is a polynomial if and only if （λ＋p, α） is a nonnegative integer. From this, we can recover a formula of singular vectors given by Malikov et al..

UNSTABLE SOLUTIONS TO A CLASS OF VECTOR DIFFERENTIAL EQUATIONS OF SIXTH ORDER

In this paper, we consider a class of nonlinear vector differential equations of sixth order. By constructing appropriate Lyapunov functions, the non-existence of periodic solutions is established. Moreover, we provide an example to show the feasibility of our results. Our results extend and improve two related results in the previous literature from scalar cases to vectorial cases.

ON THE INSTABILITY OF SOLUTIONS TO A NONLINEAR VECTOR DIFFERENTIAL EQUATION OF FOURTH ORDER

This paper presents a new result related to the instability of the zero solution to a nonlinear vector differential equation of fourth order.Our result includes and improves an instability result in the previous literature,which is related to the instability of the zero solution to a nonlinear scalar differential equation of fourth order.

ON STABILITY OF SOLUTION OF A CERTAIN FOURTH ORDER VECTOR DIFFERENTIAL EQUATIONS

This paper gives sufficient conditions for the global asmptotic stability of the zero solution of the differential equation (1. 1). The result improves and generalizes the wellknown results.

ON THE STABILITY OF SOLUTIONS TO A CERTAIN FOURTH-ORDER VECTOR DELAY DIFFERENTIAL EQUATION

In this paper, by constructing a Lyapunov functional, sufficient conditions for the uniform stability of the zero solution to a fourth-order vector delay differential equation are given.

Existence and Uniqueness of Fixed Point in Partially Ordered Sets and Applications to Ordinary Differential Equations

We prove some fixed point theorems in partially ordered sets, providing an extension of the Banach contractive mapping theorem. Having studied previously the nondecreasing case, we consider in this paper nonincreasing mappings as well as non monotone mappings. We also present some applications to first-order ordinary differential equations with periodic boundary conditions, proving the existence of a unique solution admitting the existence of a lower solution.

Euler-type schemes for weakly coupled forward-backward stochastic differential equations and optimal convergence analysis

We introduce a new Euler-type scheme and its iterative algorithm for solving weakly coupled forward-backward stochastic differential equations （FBSDEs）. Although the schemes share some common features with the ones proposed by C. Bender and J. Zhang [Ann. Appl. Probab., 2008, 18： 143-177], less computational work is needed for our method. For both our schemes and the ones proposed by Bender and Zhang, we rigorously obtain first-order error estimates, which improve the half-order error estimates of Bender and Zhang. Moreover, numerical tests are given to demonstrate the first-order accuracy of the schemes.

EXISTENCE AND UNIQUENESS FOR SECOND-ORDER VECTOR BOUNDARY VALUE PROBLEM OF NONLINEAR SYSTEMS

This paper is concerned with the following second-order vector boundary value problem ：x^R=f（t,Sx,x,x＇）,0〈t〈1,x（0）=A,g（x（1）,x＇（1））=B,where x,f,g,A and B are n-vectors. Under appropriate assumptions,existence and uniqueness of solutions are obtained by using upper and lower solutions method.