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 g...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.展开更多
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...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.展开更多
In this paper,we investigate the growth of meromorphic solutions of some kind of non-homogeneous linear difference equations with special meromorphic coefficients.When there are more than one coefficient having the sa...In this paper,we investigate the growth of meromorphic solutions of some kind of non-homogeneous linear difference equations with special meromorphic coefficients.When there are more than one coefficient having the same maximal order and the same maximal type,the estimates on the lower bound of the order of meromorphic solutions of the involved equations are obtained.Meanwhile,the above estimates are sharpened by combining the relative results of the corresponding homogeneous linear difference equations.展开更多
In this paper, we study the threshold result for the initial boundary value problem of non-homogeneous semilinear parabolic equations {μt-△μ=g(μ)+λf(x),(x,t)∈Ω×(0,T),μ=0,(x,t)∈ Ω×[0,T)...In this paper, we study the threshold result for the initial boundary value problem of non-homogeneous semilinear parabolic equations {μt-△μ=g(μ)+λf(x),(x,t)∈Ω×(0,T),μ=0,(x,t)∈ Ω×[0,T),μ(x,0)=μ0(x)≥0,x∈Ω.By combining a priori estimate of global solution with property of stationary solution set of problem (P), we prove that the minimal stationary solution Uλ(x) of problem (P) is stable, whereas, any other stationary solution is an initial datum threshold for the existence and nonexistence of global solution to problem (P).展开更多
This note contains three main results.Firstly,a complete solution of the Linear Non-Homogeneous Matrix Differential Equations(LNHMDEs)is presented that takes into account both the non-zero initial conditions of the ps...This note contains three main results.Firstly,a complete solution of the Linear Non-Homogeneous Matrix Differential Equations(LNHMDEs)is presented that takes into account both the non-zero initial conditions of the pseudo state and the nonzero initial conditions of the input.Secondly,in order to characterise the dynamics of the LNHMDEs correctly,some important concepts such as the state,slow state(smooth state)and fast state(impulsive state)are generalized to the LNHMDE case and the solution of the LNHMDEs is separated into the smooth(slow)response and the fast(implusive)response.As a third result,a new characterization of the impulsive free initial conditions of the LNHMDEs is given.展开更多
In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled b...In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.展开更多
Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of ...Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of linear second-order partial differential equations. This paper discusses the Hessian and its applications in optimization, and then proceeds to introduce and derive the notion of the Jaffa Transform, a new linear operator that directly maps a Hessian square matrix space to the initial corresponding scalar field in nth dimensional Euclidean space. The Jaffa Transform is examined, including the properties of the operator, the transform of notable matrices, and the existence of an inverse Jaffa Transform, which is, by definition, the Hessian matrix operator. The Laplace equation is then noted and investigated, particularly, the relation of the Laplace equation to Poisson’s equation, and the theoretical applications and correlations of harmonic functions to Hessian matrices. The paper concludes by introducing and explicating the Jaffa Theorem, a principle that declares the existence of harmonic Jaffa Transforms, which are, essentially, Jaffa Transform solutions to the Laplace partial differential equation.展开更多
In this work,we study the linearized Landau equation with soft potentials and show that the smooth solution to the Cauchy problem with initial datum in L^(2)(ℝ^(3))enjoys an analytic regularization effect,and that the...In this work,we study the linearized Landau equation with soft potentials and show that the smooth solution to the Cauchy problem with initial datum in L^(2)(ℝ^(3))enjoys an analytic regularization effect,and that the evolution of the analytic radius is the same as the heat equations.展开更多
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 ...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.展开更多
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 c...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.展开更多
Let M be an n-dimensional complete noncompact Riemannian manifold. In this paper, we will give the elliptic gradient estimate for positive smooth solutions to the non-homogeneous heat equation(?_t-△)u(x, t) = A(x, t)...Let M be an n-dimensional complete noncompact Riemannian manifold. In this paper, we will give the elliptic gradient estimate for positive smooth solutions to the non-homogeneous heat equation(?_t-△)u(x, t) = A(x, t)when the metric evolves under the Ricci flow. As applications, we get Harnack inequalities to compare solutions at the same time.展开更多
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 tremend...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.展开更多
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 o...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.展开更多
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...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.展开更多
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-m...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.展开更多
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 smoot...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.展开更多
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 ...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.展开更多
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 C<sup>m</sup>) and th...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 C<sup>m</sup>) and the generalized in A<sub>T,S</sub> For the special case when b ∈AT and dim(T)=dim(AT), this splitting iterative methverse A<sub>T,S</sub> hod converges to A<sub>T,S</sub>b (the unique solution of the general restricted system Ax=bx ∈T).展开更多
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)=...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.展开更多
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 iterativ...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.展开更多
基金supported by the Simons Foundation:Collaboration Grantssupported by the AFOSR grant FA9550-18-1-0383.
文摘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.
基金supported by the NSF under Grant DMS-2208391sponsored by the NSF under Grant DMS-1753581.
文摘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.
基金Supported by the National Natural Science Foundation of China(No.11761035)the Natural Science Foundation of Jiangxi Province in China(No.20171BAB201002)
文摘In this paper,we investigate the growth of meromorphic solutions of some kind of non-homogeneous linear difference equations with special meromorphic coefficients.When there are more than one coefficient having the same maximal order and the same maximal type,the estimates on the lower bound of the order of meromorphic solutions of the involved equations are obtained.Meanwhile,the above estimates are sharpened by combining the relative results of the corresponding homogeneous linear difference equations.
基金supported by Natural Science Foundation of China(10971061)Hunan Provincial Innovation Foundation For Postgraduate(CX2010B209)
文摘In this paper, we study the threshold result for the initial boundary value problem of non-homogeneous semilinear parabolic equations {μt-△μ=g(μ)+λf(x),(x,t)∈Ω×(0,T),μ=0,(x,t)∈ Ω×[0,T),μ(x,0)=μ0(x)≥0,x∈Ω.By combining a priori estimate of global solution with property of stationary solution set of problem (P), we prove that the minimal stationary solution Uλ(x) of problem (P) is stable, whereas, any other stationary solution is an initial datum threshold for the existence and nonexistence of global solution to problem (P).
文摘This note contains three main results.Firstly,a complete solution of the Linear Non-Homogeneous Matrix Differential Equations(LNHMDEs)is presented that takes into account both the non-zero initial conditions of the pseudo state and the nonzero initial conditions of the input.Secondly,in order to characterise the dynamics of the LNHMDEs correctly,some important concepts such as the state,slow state(smooth state)and fast state(impulsive state)are generalized to the LNHMDE case and the solution of the LNHMDEs is separated into the smooth(slow)response and the fast(implusive)response.As a third result,a new characterization of the impulsive free initial conditions of the LNHMDEs is given.
文摘In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.
文摘Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of linear second-order partial differential equations. This paper discusses the Hessian and its applications in optimization, and then proceeds to introduce and derive the notion of the Jaffa Transform, a new linear operator that directly maps a Hessian square matrix space to the initial corresponding scalar field in nth dimensional Euclidean space. The Jaffa Transform is examined, including the properties of the operator, the transform of notable matrices, and the existence of an inverse Jaffa Transform, which is, by definition, the Hessian matrix operator. The Laplace equation is then noted and investigated, particularly, the relation of the Laplace equation to Poisson’s equation, and the theoretical applications and correlations of harmonic functions to Hessian matrices. The paper concludes by introducing and explicating the Jaffa Theorem, a principle that declares the existence of harmonic Jaffa Transforms, which are, essentially, Jaffa Transform solutions to the Laplace partial differential equation.
基金supported by the Natural Science Foundation of Hubei Province,China (2022CFB444)the Key Laboratory of Mathematical Modelling and High Performance Computing of Air Vehicles (NUAA)+1 种基金supported by the NSFC (12031006)the Fundamental Research Funds for the Central Universities of China.
文摘In this work,we study the linearized Landau equation with soft potentials and show that the smooth solution to the Cauchy problem with initial datum in L^(2)(ℝ^(3))enjoys an analytic regularization effect,and that the evolution of the analytic radius is the same as the heat equations.
基金The supports of the National Natural Science Foundation of China(Grant Nos.51725804 and U1711264)the Research Fund for State Key Laboratories of Ministry of Science and Technology of China(SLDRCE19-B-23)the Shanghai Post-Doctoral Excellence Program(2022558)。
文摘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.
文摘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.
文摘Let M be an n-dimensional complete noncompact Riemannian manifold. In this paper, we will give the elliptic gradient estimate for positive smooth solutions to the non-homogeneous heat equation(?_t-△)u(x, t) = A(x, t)when the metric evolves under the Ricci flow. As applications, we get Harnack inequalities to compare solutions at the same time.
文摘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.
基金supported by the National Natural Science Foundation of China (11171119 and 10871076)
文摘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.
文摘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.
文摘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.
基金Project supported by the National Science Foundation of China (10271097)
文摘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.
基金the National Natural Science Foundation of China(10161006,10571044)the Natural Science Foundation of Guangdong Prov(06025059)
文摘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.
基金This project is supported by Science and Technology Foundation of Shanghai Higher Eduction,Doctoral Program Foundation of Higher Education in China.National Nature Science Foundation of China and Youth Science Foundation of Universities in Shanghai.
文摘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 C<sup>m</sup>) and the generalized in A<sub>T,S</sub> For the special case when b ∈AT and dim(T)=dim(AT), this splitting iterative methverse A<sub>T,S</sub> hod converges to A<sub>T,S</sub>b (the unique solution of the general restricted system Ax=bx ∈T).
基金Research supported by China National Science Foundation
文摘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.
基金Supported by the National Natural Science Foundation of China(61272300)
文摘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.