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.

High Order IMEX Stochastic Galerkin Schemes for Linear Transport Equation with Random Inputs and Diffusive Scalings

In this paper,we consider the high order method for solving the linear transport equations under diffusive scaling and with random inputs.To tackle the randomness in the problem,the stochastic Galerkin method of the generalized polynomial chaos approach has been employed.Besides,the high order implicit-explicit scheme under the micro-macro decomposition framework and the discontinuous Galerkin method have been employed.We provide several numerical experiments to validate the accuracy and the stochastic asymptotic-preserving property.

Stability and Time-Step Constraints of Implicit-Explicit Runge-Kutta Methods for the Linearized Korteweg-de Vries Equation

This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.

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.

SOLVERS FOR SYSTEMS OF LARGE SPARSE LINEAR AND NONLINEAR EQUATIONS BASED ON MULTI-GPUS

Numerical treatment of engineering application problems often eventually results in a solution of systems of linear or nonlinear equations.The solution process using digital computational devices usually takes tremendous time due to the extremely large size encountered in most real-world engineering applications.So,practical solvers for systems of linear and nonlinear equations based on multi graphic process units（GPUs）are proposed in order to accelerate the solving process.In the linear and nonlinear solvers,the preconditioned bi-conjugate gradient stable（PBi-CGstab）method and the Inexact Newton method are used to achieve the fast and stable convergence behavior.Multi-GPUs are utilized to obtain more data storage that large size problems need.

ZEROS OF ENTIRE SOLUTIONS TO COMPLEX LINEAR DIFFERENCE EQUATIONS

In this article, we study the complex oscillation problems of entire solutions to homogeneous and nonhomogeneous linear difference equations, and obtain some relations of the exponent of convergence of zeros and the order of growth of entire solutions to complex linear difference equations.

THE STABILITY OF LINEAR MULTISTEP METHODS FOR SYSTEMS OF DELAY DIFFERENTIAL EQUATIONS

This paper deals with the numerical solution of initial value problems for systems of differential equations with a delay argument. The numerical stability of a linear multistep method is investigated by analysing the solution of the lest equation y’(t)=Ay(t) + By(1-t),where A,B denote constant complex N×N-matrices,and t】0.We investigate carefully the characterization of the stability region.

ASYNCHRONOUS RELAXED ITERATIVE METHODS FOR SOLVING LINEAR SYSTEMS OF EQUATIONS

In this paper, the asynchronous versions of classical iterative methods for solving linear systems of equations are considered. Sufficient conditions for convergence of asynchronous relaxed processes are given for H-matrix by which nor only the requirements of [3] on coefficient matrix are lowered, but also a larger region of convergence than that in [3] is obtained.

LINEAR SINGULAR INTEGRAL EQUATION ON DOMAINS COMPOSED BY BALLS

For domains composed by balls in C^n, this paper studies the boundary behaviour of Cauchy type integrals with discrete holomorphic kernels and the corresponding linear singular integral equation on each piece of smooth lower dimensional edges on the boundary of the domain.

ESTIMATES FOR THE HYPER-ORDER OF SOLUTIONS OF CERTAIN LINEAR DIFFERENTIAL EQUATIONS

In this article, the authors study the growth of certain second order linear differential equation f″＋A（z）f′＋B（z）f=0 and give precise estimates for the hyperorder of solutions of infinite order. Under similar conditions, higher order differential equations will be considered.

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

Roughness of Exponential Estimates for Linear Functional Differential Equations

Assume that the fundamental solution matrix U (t, s ) of x’(t)=L(t, x,) satisfies |U(t,s)|≤Ke-e(t-s) for t≥s.If|(t,φ)|≤δ|φ(0)|with δ【a/K, then the fundamental solution matrix of the perturbed equation x’(t)=L(t,x,)+(t ,x,) also possesses similar exponential estimate. For α=0, a similar result is given.

Geometric interpretation of several classical iterative methods for linear system of equations and diverse relaxation parameter of the SOR method

Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterative methods and provide a reference for further study and design. Finally, a new iterative method is designed named as the diverse relaxation parameter of the SOR method which, in particular, demonstrates the geometric characteristics. Many examples prove that the method is quite effective.

FINITE TIME BLOW UP OF THE SOLUTIONS TO BOUSSINESQ EQUATION WITH LINEAR RESTORING FORCE AND ARBITRARY POSITIVE ENERGY

Finite time blow up of the solutions to Boussinesq equation with linear restoring force and combined power nonlinearities is studied. Sufficient conditions on the initial data for nonexistence of global solutions are derived. The results are valid for initial data with arbitrary high positive energy. The proofs are based on the concave method and new sign preserving functionals.

GRADIENT ESTIMATES AND LIOUVILLE THEOREMS FOR LINEAR AND NONLINEAR PARABOLIC EQUATIONS ON RIEMANNIAN MANIFOLDS

In this article, we will derive local elliptic type gradient estimates for positive solutions of linear parabolic equations （△-e/et）u（x,t）＋q（x,t）u^p（x,t）=0 and nonlinear parabolic equations （△-e/et）u（x,t）＋h（x,t）u^p（x,t）=0（p 〉 1） on Riemannian manifolds.As applications, we obtain some theorems of Liouville type for positive ancient solutions of such equations. Our results generalize that of Souplet-Zhang （[1], Bull. London Math. Soc. 38（2006）, 1045-1053） and the author （[2], Nonlinear Anal. 74 （2011）, 5141-5146）.

Exact solution of the one-dimensional Klein-Gordon equation with scalar and vector linear potentials in the presence of a minimal length

Using the momentum space representation, we solve the Klein--Gordon equation in one spatial dimension for the case of mixed scalar and vector linear potentials in the context of deformed quantum mechanics characterized by a finite minimal uncertainty in position. The expressions of bound state energies and the associated wave functions are exactly obtained.

Equivalent method for accurate solution to linear interval equations

Based on linear interval equations, an accurate interval finite element method for solving structural static problems with uncertain parameters in terms of optimization is discussed. On the premise of ensuring the consistency of solution sets, the original interval equations are equivalently transformed into some deterministic inequations. On this basis, calculating the structural displacement response with interval parameters is predigested to a number of deterministic linear optimization problems. The results are proved to be accurate to the interval governing equations. Finally, a numerical example is given to demonstrate the feasibility and efficiency of the proposed method.