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.

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.

Equation governing the probability density evolution of multi-dimensional linear fractional differential systems subject to Gaussian white noise

Stochastic fractional differential systems are important and useful in the mathematics,physics,and engineering fields.However,the determination of their probabilistic responses is difficult due to their non-Markovian property.The recently developed globally-evolving-based generalized density evolution equation(GE-GDEE),which is a unified partial differential equation(PDE)governing the transient probability density function(PDF)of a generic path-continuous process,including non-Markovian ones,provides a feasible tool to solve this problem.In the paper,the GE-GDEE for multi-dimensional linear fractional differential systems subject to Gaussian white noise is established.In particular,it is proved that in the GE-GDEE corresponding to the state-quantities of interest,the intrinsic drift coefficient is a time-varying linear function,and can be analytically determined.In this sense,an alternative low-dimensional equivalent linear integer-order differential system with exact closed-form coefficients for the original highdimensional linear fractional differential system can be constructed such that their transient PDFs are identical.Specifically,for a multi-dimensional linear fractional differential system,if only one or two quantities are of interest,GE-GDEE is only in one or two dimensions,and the surrogate system would be a one-or two-dimensional linear integer-order system.Several examples are studied to assess the merit of the proposed method.Though presently the closed-form intrinsic drift coefficient is only available for linear stochastic fractional differential systems,the findings in the present paper provide a remarkable demonstration on the existence and eligibility of GE-GDEE for the case that the original high-dimensional system itself is non-Markovian,and provide insights for the physical-mechanism-informed determination of intrinsic drift and diffusion coefficients of GE-GDEE of more generic complex nonlinear systems.

The Existence of Meromorphic Solutions to Non-Linear Delay Differential Equations

In this paper, we study the existence of the transcendental meromorphic solution of the delay differential equations , where a(z) is a rational function, and are polynomials in w(z) with rational coefficients, k is a positive integer. Under the assumption when above equations own transcendental meromorphic solutions with minimal hyper-type, we derive the concrete conditions on the degree of the right side of them. Specially, when w(z)=0 is a root of , its multiplicity is at most k. Some examples are given here to illustrate that our results are accurate.

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.

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.

Novel method based on ant colony opti mization for solving ill-conditioned linear systems of equations

A novel method based on ant colony optimization （ACO）, algorithm for solving the ill-conditioned linear systems of equations is proposed. ACO is a parallelized bionic optimization algorithm which is inspired from the behavior of real ants. ACO algorithm is first introduced, a kind of positive feedback mechanism is adopted in ACO. Then, the solu- tion problem of linear systems of equations was reformulated as an unconstrained optimization problem for solution by an ACID algorithm. Finally, the ACID with other traditional methods is applied to solve a kind of multi-dimensional Hilbert ill-conditioned linear equations. The numerical results demonstrate that ACO is effective, robust and recommendable in solving ill-conditioned linear systems of equations.

Balancing Chemical Equations by Systems of Linear Equations

In this paper, a formal and systematic method for balancing chemical reaction equations was presented. The results satisfy the law of conservation of matter, and confirm that there is no contradiction to the existing way(s) of balancing chemical equations. A chemical reaction which possesses atoms with fractional oxidation numbers that have unique coefficients was studied. In this paper, the chemical equations were balanced by representing the chemical equation into systems of linear equations. Particularly, the Gauss elimination method was used to solve the mathematical problem with this method, it was possible to handle any chemical reaction with given reactants and products.

Analytical Solutions of System of Non-Linear Differential Equations in the Single-Enzyme, Single-Substrate Reaction with Non-Mechanism-Based Enzyme Inactivation

A closed form of an analytical expression of concentration in the single-enzyme, single-substrate system for the full range of enzyme activities has been derived. The time dependent analytical solution for substrate, enzyme-substrate complex and product concentrations are presented by solving system of non-linear differential equation. We employ He’s Homotopy perturbation method to solve the coupled non-linear differential equations containing a non-linear term related to basic enzymatic reaction. The time dependent simple analytical expressions for substrate, enzyme-substrate and free enzyme concentrations have been derived in terms of dimensionless reaction diffusion parameters ε, λ1, λ2 and λ3 using perturbation method. The numerical solution of the problem is also reported using SCILAB software program. The analytical results are compared with our numerical results. An excellent agreement with simulation data is noted. The obtained results are valid for the whole solution domain.

The <i>LA = U</i>Decomposition Method for Solving Systems of Linear Equations

A method for solving systems of linear equations is presented based on direct decomposition of the coefficient matrix using the form LAX = LB = B’ . Elements of the reducing lower triangular matrix L can be determined using either row wise or column wise operations and are demonstrated to be sums of permutation products of the Gauss pivot row multipliers. These sums of permutation products can be constructed using a tree structure that can be easily memorized or alternatively computed using matrix products. The method requires only storage of the L matrix which is half in size compared to storage of the elements in the LU decomposition. Equivalence of the proposed method with both the Gauss elimination and LU decomposition is also shown in this paper.

Criteria for System of Three Second-Order Ordinary Differential Equations to Be Reduced to a Linear System via Restricted Class of Point Transformation

This paper is devoted to the study of the linearization problem of system of three second-order ordinary differential equations and . The necessary conditions for linearization by general point transformation and are found. The sufficient conditions for linearization by restricted class of point transformation and are obtained. Moreover, the procedure for obtaining the linearizing transformation is provided in explicit forms. Examples demonstrating the procedure of using the linearization theorems are presented.

The Continuous Analogy of Newton’s Method for Solving a System of Linear Algebraic Equations

We propose a continuous analogy of Newton’s method with inner iteration for solving a system of linear algebraic equations. Implementation of inner iterations is carried out in two ways. The former is to fix the number of inner iterations in advance. The latter is to use the inexact Newton method for solution of the linear system of equations that arises at each stage of outer iterations. We give some new choices of iteration parameter and of forcing term, that ensure the convergence of iterations. The performance and efficiency of the proposed iteration is illustrated by numerical examples that represent a wide range of typical systems.

A Mathod to Solve Systems of Fuzzy Linear Equations

Many systems of fuzzy linear equations do not have solutions when the solution concept is based on α cuts and interval arithmetic. In this paper,we establish the relations between the systems of fuzzy linear equations and the possibilistic linear programming problems and present an alternative method of solving the systems of fuzzy linear equations.

Certain Algebraic Test for Analyzing Aperiodic Stability of Two-Dimensional Linear Discrete Systems

This paper addresses the new algebraic test to check the aperiodic stability of two dimensional linear time invariant discrete systems. Initially, the two dimensional characteristics equations are converted into equivalent one-dimensional equation. Further Fuller’s idea is applied on the equivalent one-dimensional characteristics equation. Then using the co-efficient of the characteristics equation, the routh table is formed to ascertain the aperiodic stability of the given two-dimensional linear discrete system. The illustrations were presented to show the applicability of the proposed technique.

A Modified Stability Analysis of Two-Dimensional Linear Time Invariant Discrete Systems within the Unity-Shifted Unit Circle

This paper proposes a method to ascertain the stability of two dimensional linear time invariant discrete system within the shifted unit circle which is represented by the form of characteristic equation. Further an equivalent single dimensional characteristic equation is formed from the two dimensional characteristic equation then the stability formulation in the left half of Z-plane, where the roots of characteristic equation f(Z) = 0 should lie within the shifted unit circle. The coefficient of the unit shifted characteristic equation is suitably arranged in the form of matrix and the inner determinants are evaluated using proposed Jury’s concept. The proposed stability technique is simple and direct. It reduces the computational cost. An illustrative example shows the applicability of the proposed scheme.

SOLUTIONS FOR A SYSTEM OF NONLINEAR RANDOM INTEGRAL AND DIFFERENTIAL EQUATIONS

In this paper we first prove a Darbao type fixed point theorem for a system of continuous random operators with random domains. Thenb, by using the theorem. wegive the existence criteria of solutions for a systems of nonlinear random Volterraintegral equations and for the Cauchy problem of a system of nonlinear random differential equations. The existence of extremal random solutions and random comparison results for these systems of random equations are also obtained Our theorems improve and generalize the corresponding results of Vaughn Lakshmikantham Lakshmidantham-Leela De blasi-Myjak and Ding

THE ESTIMATION OF SOLUTION OF THE BOUNDARY VALUE PROBLEM OF THE SYSTEMS FOR QUASI-LINEAR ORDINARY DIFFERENTIAL EQUATIONS

This paper deals with the singular perturbation of the boundary value problem of the systems for quasi-linear ordinary differential equationswhere x,f, y , h, A, B and C all belong to Rn , and g is an n×n matrix function. Under suitable conditions we prove the existence of the solutions by diagonalization and the fixed point theorem and also estimate the remainder.

Domain Decomposition of an Optimal Control Problem for Semi-Linear Elliptic Equations on Metric Graphs with Application to Gas Networks

We consider optimal control problems for the flow of gas in a pipe network. The equations of motions are taken to be represented by a semi-linear model derived from the fully nonlinear isothermal Euler gas equations. We formulate an optimal control problem on a given network and introduce a time discretization thereof. We then study the well-posedness of the corresponding time-discrete optimal control problem. In order to further reduce the complexity, we consider an instantaneous control strategy. The main part of the paper is concerned with a non-overlapping domain decomposition of the semi-linear elliptic optimal control problem on the graph into local problems on a small part of the network, ultimately on a single edge.

THE SOLUTION FOR THE GENERALIZED RICCATIALGEBRAIC EQUATIONS OF LINEAR EQUALITY CONSTRAINT SYSTEM

Based on the dynamic equation, the performance functional and the system constraint equation of time-invariant discrete LQ control problem, the generalized Riccati equations of linear equality constraint system are obtained according to the minimum principle, then a deep discussion about the above equations is given, and finally numerical example is shown in this paper.