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.

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.

Numerical Analysis of Upwind Difference Schemes for Two-Dimensional First-Order Hyperbolic Equations with Variable Coefficients

In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.

A POSTERIORI ERROR ESTIMATE OF THE DSD METHOD FOR FIRST-ORDER HYPERBOLIC EQUATIONS

A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illustrate the accuracy and feasibility of this method.

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.

P-and S-wavefield simulations using both the firstand second-order separated wave equations through a high-order staggered grid finite-difference method

In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.

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.

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.

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）.

Numerical modeling of seismic wavefields in transversely isotropic media with a compact staggered-grid finite difference scheme

To deal with the numerical dispersion problem, by combining the staggeredgrid technology with the compact finite difference scheme, we derive a compact staggered- grid finite difference scheme from the first-order velocity-stress wave equations for the transversely isotropic media. Comparing the principal truncation error terms of the compact staggered-grid finite difference scheme, the staggered-grid finite difference scheme, and the compact finite difference scheme, we analyze the approximation accuracy of these three schemes using Fourier analysis. Finally, seismic wave numerical simulation in transversely isotropic （VTI） media is performed using the three schemes. The results indicate that the compact staggered-grid finite difference scheme has the smallest truncation error, the highest accuracy, and the weakest numerical dispersion among the three schemes. In summary, the numerical modeling shows the validity of the compact staggered-grid finite difference scheme.

A study of perfectly matched layers for joint multicomponent reverse-time migration

Reverse-time migration in finite space requires effective boundary processing technology to eliminate the artificial truncation boundary effect in the migration result.On the basis of the elastic velocity-stress equations in vertical transversely isotropic media and the idea of the conventional split perfectly matched layer（PML）,the PML wave equations in reverse-time migration are derived in this paper and then the high order staggered grid discrete schemes are subsequently given.Aiming at the＂reflections＂from the boundary to the computational domain,as well as the effect of seismic event＇s abrupt changes at the two ends of the seismic array,the PML arrangement in reverse-time migration is given.The synthetic and real elastic,prestack,multi-component,reverse-time depth migration results demonstrate that this method has much better absorbing effects than other methods and the joint migration produces good imaging results.

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.

K-Dimensional Optimal Parallel Algorithm for the Solution of a General Class of Recurrence Equations

This paper proposes a parallel algorithm, called KDOP (K-DimensionalOptimal Parallel algorithm), to solve a general class of recurrence equations efficiently. The KDOP algorithm partitions the computation into a series of sub-computations, each of which is executed in the fashion that all the processors work simultaneously with each one executing an optimal sequential algorithm to solve a subcomputation task. The algorithm solves the equations in O(N/p)steps in EREW PRAM model (Exclusive Read Exclusive Write Parallel Ran-dom Access Machine model) using p<N1-e processors, where N is the size of the problem, and e is a given constant. This is an optimal algorithm (itsspeedup is O(p)) in the case of p<N1-e. Such an optimal speedup for this problem was previously achieved only in the case of p<N0.5. The algorithm can be implemented on machines with multiple processing elements or pipelined vector machines with parallel memory systems.

Precise integration method for solving singular perturbation problems

This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.

A Self-Stabilized Field Theory of Neutrinos

In “A Self-linking Field Formalism” I establish a self-dual field structure with higher order self-induced symmetries that reinforce the first-order dynamics. The structure was derived from Gauss-linking integrals in R3 based on the Biot-Savart law and Ampere’s law applied to Heaviside’s equations, derived in strength-independent fashion in “Primordial Principle of Self-Interaction”. The derivation involves Geometric Calculus, topology, and field equations. My goal in this paper is to derive the simplest solution of a self-stabilized solitonic structure and discuss this model of a neutrino.

A New Modification of the Method of Lines for First Order Hyperbolic PDEs

A new modification of the Method of Lines is proposed for the solution of first order partial differential equations. The accuracy of the method is shown with the matrix analysis. The method is applied to a number of test problems, on uniform grids, to compare the accuracy and computational efficiency with the standard method.

VPVnet:A Velocity-Pressure-Vorticity Neural Network Method for the Stokes’Equations under Reduced Regularity

We present VPVnet,a deep neural network method for the Stokes’equa-tions under reduced regularity.Different with recently proposed deep learning meth-ods[40,51]which are based on the original form of PDEs,VPVnet uses the least square functional of thefirst-order velocity-pressure-vorticity(VPV)formulation([30])as loss functions.As such,onlyfirst-order derivative is required in the loss functions,hence the method is applicable to a much larger class of problems,e.g.problems with non-smooth solutions.Despite that several methods have been proposed recently to reduce the regularity requirement by transforming the original problem into a corresponding variational form,while for the Stokes’equations,the choice of approximating spaces for the velocity and the pressure has to satisfy the LBB condition additionally.Here by making use of the VPV formulation,lower regularity requirement is achieved with no need for considering the LBB condition.Convergence and error estimates have been established for the proposed method.It is worth emphasizing that the VPVnet method is divergence-free and pressure-robust,while classical inf-sup stable mixedfinite elements for the Stokes’equations are not pressure-robust.Various numerical experiments including 2D and 3D lid-driven cavity test cases are conducted to demon-strate its efficiency and accuracy.