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.展开更多
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.展开更多
On the basis of experimental observations on animals, applications to clinical data on patients and theoretical statistical reasoning, the author developed a com-puter-assisted general mathematical model of the ‘prob...On the basis of experimental observations on animals, applications to clinical data on patients and theoretical statistical reasoning, the author developed a com-puter-assisted general mathematical model of the ‘probacent’-probability equation, Equation (1) and death rate (mortality probability) equation, Equation (2) derivable from Equation (1) that may be applica-ble as a general approximation method to make use-ful predictions of probable outcomes in a variety of biomedical phenomena [1-4]. Equations (1) and (2) contain a constant, γ and c, respectively. In the pre-vious studies, the author used the least maximum- difference principle to determine these constants that were expected to best fit reported data, minimizing the deviation. In this study, the author uses the method of computer-assisted least sum of squares to determine the constants, γ and c in constructing the ‘probacent’-related formulas best fitting the NCHS- reported data on survival probabilities and death rates in the US total adult population for 2001. The results of this study reveal that the method of com-puter-assisted mathematical analysis with the least sum of squares seems to be simple, more accurate, convenient and preferable than the previously used least maximum-difference principle, and better fit-ting the NCHS-reported data on survival probabili-ties and death rates in the US total adult population. The computer program of curved regression for the ‘probacent’-probability and death rate equations may be helpful in research in biomedicine.展开更多
Using Jacobi elliptic function linear superposition approach for the(1+1)-dimensional Caudrey-Dodd-Gibbon-Sawada-Kotera(CDGSK)equation and the(2+ 1)-dimensional Nizhnik-Novikov-Veselov(NNV)equation,many new periodic t...Using Jacobi elliptic function linear superposition approach for the(1+1)-dimensional Caudrey-Dodd-Gibbon-Sawada-Kotera(CDGSK)equation and the(2+ 1)-dimensional Nizhnik-Novikov-Veselov(NNV)equation,many new periodic travelling wave solutions with different periods and velocities are obtained based on the known periodic solutions.This procedure is crucially dependent on a sequence of cyclic identities involving Jacobi elliptic functions sn(ζ,m),cn(ζ,m),and dn(ζ,m).展开更多
In this article, we establish the exponential time decay of smooth solutions around a global Maxwellian to the non-linear Vlasov–Poisson–Fokker–Planck equations in the whole space by uniform-in-time energy estimate...In this article, we establish the exponential time decay of smooth solutions around a global Maxwellian to the non-linear Vlasov–Poisson–Fokker–Planck equations in the whole space by uniform-in-time energy estimates. The non-linear coupling of macroscopic part and Fokker–Planck operator in the model brings new difficulties for the energy estimates, which is resolved by adding tailored weighted-in-v energy estimates suitable for the Fokker–Planck operator.展开更多
The following article has been retracted due to the investigation of complaints received against it. Mr. Mohammadali Ghorbani (corresponding author and also the last author) cheated the authors’ name: Alireza Heidari...The following article has been retracted due to the investigation of complaints received against it. Mr. Mohammadali Ghorbani (corresponding author and also the last author) cheated the authors’ name: Alireza Heidari and Seyedali Vedad. The scientific community takes a very strong view on this matter and we treat all unethical behavior such as plagiarism seriously. This paper published in Vol.3 No.5 420-429, 2012, has been removed from this site.展开更多
The non-linear Fokker-Planck equation arises in describing the evolution of stochastic system, which is a variant of the Boltzmann equation modeling the evolution of the random system with Brownian motion, where the c...The non-linear Fokker-Planck equation arises in describing the evolution of stochastic system, which is a variant of the Boltzmann equation modeling the evolution of the random system with Brownian motion, where the collision term is replaced by a drift-diffusion operator. This model conserves mass, momentum and energy;the dissipation is much weaker than that in a simplified model we considered before which conserved only mass, thus more difficult to analyze. The macro-micro decomposition of the solution around the local Maxwellian introduced by T.-P. Liu, T. Yang and S.-H. Yu for Boltzmann equation is used, to reformulate the model into a fluid-type system incorporate viscosity and heat diffusion terms, coupled with an equation of the microscopic part. The viscosity and heat diffusion terms can give dissipative mechanism for the analysis of the model.展开更多
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.展开更多
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 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.展开更多
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.展开更多
基金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.
文摘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.
文摘On the basis of experimental observations on animals, applications to clinical data on patients and theoretical statistical reasoning, the author developed a com-puter-assisted general mathematical model of the ‘probacent’-probability equation, Equation (1) and death rate (mortality probability) equation, Equation (2) derivable from Equation (1) that may be applica-ble as a general approximation method to make use-ful predictions of probable outcomes in a variety of biomedical phenomena [1-4]. Equations (1) and (2) contain a constant, γ and c, respectively. In the pre-vious studies, the author used the least maximum- difference principle to determine these constants that were expected to best fit reported data, minimizing the deviation. In this study, the author uses the method of computer-assisted least sum of squares to determine the constants, γ and c in constructing the ‘probacent’-related formulas best fitting the NCHS- reported data on survival probabilities and death rates in the US total adult population for 2001. The results of this study reveal that the method of com-puter-assisted mathematical analysis with the least sum of squares seems to be simple, more accurate, convenient and preferable than the previously used least maximum-difference principle, and better fit-ting the NCHS-reported data on survival probabili-ties and death rates in the US total adult population. The computer program of curved regression for the ‘probacent’-probability and death rate equations may be helpful in research in biomedicine.
文摘Using Jacobi elliptic function linear superposition approach for the(1+1)-dimensional Caudrey-Dodd-Gibbon-Sawada-Kotera(CDGSK)equation and the(2+ 1)-dimensional Nizhnik-Novikov-Veselov(NNV)equation,many new periodic travelling wave solutions with different periods and velocities are obtained based on the known periodic solutions.This procedure is crucially dependent on a sequence of cyclic identities involving Jacobi elliptic functions sn(ζ,m),cn(ζ,m),and dn(ζ,m).
基金partially supported by Fundamental Research Funds for the Central Universities,NSFC(11871335)by the SJTU’s SMC Projection
文摘In this article, we establish the exponential time decay of smooth solutions around a global Maxwellian to the non-linear Vlasov–Poisson–Fokker–Planck equations in the whole space by uniform-in-time energy estimates. The non-linear coupling of macroscopic part and Fokker–Planck operator in the model brings new difficulties for the energy estimates, which is resolved by adding tailored weighted-in-v energy estimates suitable for the Fokker–Planck operator.
文摘The following article has been retracted due to the investigation of complaints received against it. Mr. Mohammadali Ghorbani (corresponding author and also the last author) cheated the authors’ name: Alireza Heidari and Seyedali Vedad. The scientific community takes a very strong view on this matter and we treat all unethical behavior such as plagiarism seriously. This paper published in Vol.3 No.5 420-429, 2012, has been removed from this site.
文摘The non-linear Fokker-Planck equation arises in describing the evolution of stochastic system, which is a variant of the Boltzmann equation modeling the evolution of the random system with Brownian motion, where the collision term is replaced by a drift-diffusion operator. This model conserves mass, momentum and energy;the dissipation is much weaker than that in a simplified model we considered before which conserved only mass, thus more difficult to analyze. The macro-micro decomposition of the solution around the local Maxwellian introduced by T.-P. Liu, T. Yang and S.-H. Yu for Boltzmann equation is used, to reformulate the model into a fluid-type system incorporate viscosity and heat diffusion terms, coupled with an equation of the microscopic part. The viscosity and heat diffusion terms can give dissipative mechanism for the analysis of the model.
基金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.
基金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.
文摘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.
基金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.