The difference discrete system of Euler-beam with arbitrary supports was constructed by using the two order central difference formulas. This system is equivalent to the spring-mass-rigidrod model. By using the theory...The difference discrete system of Euler-beam with arbitrary supports was constructed by using the two order central difference formulas. This system is equivalent to the spring-mass-rigidrod model. By using the theory of oscillatory matrix, the signoscillatory property of stiffness matrices of this system was proved, and the necessary and sufficient condition for the system to be positive was obtained completely.展开更多
The discrete scheme called discrete operator difference for differential equations was given. Several difference elements for plate bending problems and plane problems were given. By investigating these elements, the ...The discrete scheme called discrete operator difference for differential equations was given. Several difference elements for plate bending problems and plane problems were given. By investigating these elements, the ability of the discrete forms expressing to the element functions was talked about. In discrete operator difference method, the displacements of the elements can be reproduced exactly in the discrete forms whether the displacements are conforming or not. According to this point, discrete operator difference method is a method with good performance.展开更多
We introduce the Euler-Lagrange cohomology to study the symplectic and multisymplectic structures and their preserving properties in finite and infinite dimensional Lagrangian systems respectively. We also explore the...We introduce the Euler-Lagrange cohomology to study the symplectic and multisymplectic structures and their preserving properties in finite and infinite dimensional Lagrangian systems respectively. We also explore their certain difference discrete counterparts in the relevant regularly discretized finite and infinite dimensional Lagrangian systems by means of the difference discrete variational principle with the difference being regarded as an entire geometric object and the noncommutative differential calculus on regular lattice. In order to show that in all these cases the symplectic and multisymplectic preserving properties do not necessarily depend on the relevant Euler-Lagrange equations, the Euler-Lagrange cohomological concepts and content in the configuration space are employed.展开更多
In this paper,a new finite element and finite difference(FE-FD)method has been developed for anisotropic parabolic interface problems with a known moving interface using Cartesian meshes.In the spatial discretization,...In this paper,a new finite element and finite difference(FE-FD)method has been developed for anisotropic parabolic interface problems with a known moving interface using Cartesian meshes.In the spatial discretization,the standard P,FE discretization is applied so that the part of the coefficient matrix is symmetric positive definite,while near the interface,the maximum principle preserving immersed interface discretization is applied.In the time discretization,a modified Crank-Nicolson discretization is employed so that the hybrid FE-FD is stable and second order accurate.Correction terms are needed when the interface crosses grid lines.The moving interface is represented by the zero level set of a Lipschitz continuous function.Numerical experiments presented in this paper confirm second orderconvergence.展开更多
In this paper,we focus on inferring graph Laplacian matrix from the spatiotemporal signal which is defined as“time-vertex signal”.To realize this,we first represent the signals on a joint graph which is the Cartesia...In this paper,we focus on inferring graph Laplacian matrix from the spatiotemporal signal which is defined as“time-vertex signal”.To realize this,we first represent the signals on a joint graph which is the Cartesian product graph of the time-and vertex-graphs.By assuming the signals follow a Gaussian prior distribution on the joint graph,a meaningful representation that promotes the smoothness property of the joint graph signal is derived.Furthermore,by decoupling the joint graph,the graph learning framework is formulated as a joint optimization problem which includes signal denoising,timeand vertex-graphs learning together.Specifically,two algorithms are proposed to solve the optimization problem,where the discrete second-order difference operator with reversed sign(DSODO)in the time domain is used as the time-graph Laplacian operator to recover the signal and infer a vertex-graph in the first algorithm,and the time-graph,as well as the vertex-graph,is estimated by the other algorithm.Experiments on both synthetic and real-world datasets demonstrate that the proposed algorithms can effectively infer meaningful time-and vertex-graphs from noisy and incomplete data.展开更多
For a sparse non-singular matrix A, generally A- 1 is a dense matrix. However, for a class of matrices, A-1 can be a matrix with off-diagonal decay properties, i.e., |Aij^-1| decays fast to 0 with respect to the inc...For a sparse non-singular matrix A, generally A- 1 is a dense matrix. However, for a class of matrices, A-1 can be a matrix with off-diagonal decay properties, i.e., |Aij^-1| decays fast to 0 with respect to the increase of a properly defined distance between i and j. Here we consider the off-diagonal decay properties of discretized Green's functions for SchrSdinger type operators. We provide decay estimates for discretized Green's functions obtained from the finite difference discretization, and from a variant of the pseudo-spectral discretization. The asymptotic decay rate in our estimate is independent of the domain size and of the discretization parameter. We verify the decay estimate with numerical results for one-dimensional Schr6dinger type operators.展开更多
During the past years,the recovery of unconventional gas formation has attracted lots of attention and achieved huge success.To produce gas from the low-permeability unconventional formations,hydraulic fracturing tech...During the past years,the recovery of unconventional gas formation has attracted lots of attention and achieved huge success.To produce gas from the low-permeability unconventional formations,hydraulic fracturing technology is essential and critical.In this paper,we present the development of a three-dimensional thermalhydraulic-mechanical numerical simulator for the simulation of hydraulic fracturing operations in tight sandstone reservoirs.Our simulator is based on integrated finite difference(IFD)method.In this method,the simulation domain is subdivided into sub domains and the governing equations are integrated over a sub domain with flux terms expressed as an integral over the sub domain boundary using the divergence theorem.Our simulator conducts coupled thermal-hydraulic-mechanical simulation of the initiation and extension of hydraulic fractures.It also calculates the mass/heat transport of injected hydraulic fluids as well as proppants.Our simulator is able to handle anisotropic formations with multiple layers.Our simulator has been validated by comparing with an analytical solution as well as Ribeiro and Sharma model.Our model can simulate fracture spacing effect on fracture profile when combining IFD with Discontinuous Displacement Method(DDM).展开更多
基金Project supported by the National Natural Science Foundation of China (No.60034010)
文摘The difference discrete system of Euler-beam with arbitrary supports was constructed by using the two order central difference formulas. This system is equivalent to the spring-mass-rigidrod model. By using the theory of oscillatory matrix, the signoscillatory property of stiffness matrices of this system was proved, and the necessary and sufficient condition for the system to be positive was obtained completely.
文摘The discrete scheme called discrete operator difference for differential equations was given. Several difference elements for plate bending problems and plane problems were given. By investigating these elements, the ability of the discrete forms expressing to the element functions was talked about. In discrete operator difference method, the displacements of the elements can be reproduced exactly in the discrete forms whether the displacements are conforming or not. According to this point, discrete operator difference method is a method with good performance.
文摘We introduce the Euler-Lagrange cohomology to study the symplectic and multisymplectic structures and their preserving properties in finite and infinite dimensional Lagrangian systems respectively. We also explore their certain difference discrete counterparts in the relevant regularly discretized finite and infinite dimensional Lagrangian systems by means of the difference discrete variational principle with the difference being regarded as an entire geometric object and the noncommutative differential calculus on regular lattice. In order to show that in all these cases the symplectic and multisymplectic preserving properties do not necessarily depend on the relevant Euler-Lagrange equations, the Euler-Lagrange cohomological concepts and content in the configuration space are employed.
基金partially supported by the National Natural Science Foundation of China(Grant No.12261070)the Ningxia Key Research and Development Project of China(Grant No.2022BSB03048)+2 种基金partially supported by the Simons(Grant No.633724)and by Fundacion Seneca grant 21760/IV/22partially supported by the Spanish national research project PID2019-108336GB-I00by Fundacion Séneca grant 21728/EE/22.Este trabajo es resultado de las estancias(21760/IV/22)y(21728/EE/22)financiadas por la Fundacion Séneca-Agencia de Ciencia y Tecnologia de la Region de Murcia con cargo al Programa Regional de Movilidad,Colaboracion Internacional e Intercambio de Conocimiento"Jimenez de la Espada".(Plan de Actuacion 2022).
文摘In this paper,a new finite element and finite difference(FE-FD)method has been developed for anisotropic parabolic interface problems with a known moving interface using Cartesian meshes.In the spatial discretization,the standard P,FE discretization is applied so that the part of the coefficient matrix is symmetric positive definite,while near the interface,the maximum principle preserving immersed interface discretization is applied.In the time discretization,a modified Crank-Nicolson discretization is employed so that the hybrid FE-FD is stable and second order accurate.Correction terms are needed when the interface crosses grid lines.The moving interface is represented by the zero level set of a Lipschitz continuous function.Numerical experiments presented in this paper confirm second orderconvergence.
基金supported by the National Natural Science Foundation of China(Grant No.61966007)Key Laboratory of Cognitive Radio and Information Processing,Ministry of Education(No.CRKL180106,No.CRKL180201)+1 种基金Guangxi Key Laboratory of Wireless Wideband Communication and Signal Processing,Guilin University of Electronic Technology(No.GXKL06180107,No.GXKL06190117)Guangxi Colleges and Universities Key Laboratory of Satellite Navigation and Position Sensing.
文摘In this paper,we focus on inferring graph Laplacian matrix from the spatiotemporal signal which is defined as“time-vertex signal”.To realize this,we first represent the signals on a joint graph which is the Cartesian product graph of the time-and vertex-graphs.By assuming the signals follow a Gaussian prior distribution on the joint graph,a meaningful representation that promotes the smoothness property of the joint graph signal is derived.Furthermore,by decoupling the joint graph,the graph learning framework is formulated as a joint optimization problem which includes signal denoising,timeand vertex-graphs learning together.Specifically,two algorithms are proposed to solve the optimization problem,where the discrete second-order difference operator with reversed sign(DSODO)in the time domain is used as the time-graph Laplacian operator to recover the signal and infer a vertex-graph in the first algorithm,and the time-graph,as well as the vertex-graph,is estimated by the other algorithm.Experiments on both synthetic and real-world datasets demonstrate that the proposed algorithms can effectively infer meaningful time-and vertex-graphs from noisy and incomplete data.
基金supported by Laboratory Directed Research and Development Funding from Berkeley Labprovided by the Director,Office of Science,of the US Department of Energy(Grant No.DE-AC02-05CH11231)+3 种基金the Alfred P Sloan Foundationthe DOE Scientific Discovery through the Advanced Computing Programthe DOE Center for Applied Mathematics for Energy Research Applications Programthe National Science Foundation of USA(Grant Nos.DMS-1312659 and DMS-1454939)
文摘For a sparse non-singular matrix A, generally A- 1 is a dense matrix. However, for a class of matrices, A-1 can be a matrix with off-diagonal decay properties, i.e., |Aij^-1| decays fast to 0 with respect to the increase of a properly defined distance between i and j. Here we consider the off-diagonal decay properties of discretized Green's functions for SchrSdinger type operators. We provide decay estimates for discretized Green's functions obtained from the finite difference discretization, and from a variant of the pseudo-spectral discretization. The asymptotic decay rate in our estimate is independent of the domain size and of the discretization parameter. We verify the decay estimate with numerical results for one-dimensional Schr6dinger type operators.
文摘During the past years,the recovery of unconventional gas formation has attracted lots of attention and achieved huge success.To produce gas from the low-permeability unconventional formations,hydraulic fracturing technology is essential and critical.In this paper,we present the development of a three-dimensional thermalhydraulic-mechanical numerical simulator for the simulation of hydraulic fracturing operations in tight sandstone reservoirs.Our simulator is based on integrated finite difference(IFD)method.In this method,the simulation domain is subdivided into sub domains and the governing equations are integrated over a sub domain with flux terms expressed as an integral over the sub domain boundary using the divergence theorem.Our simulator conducts coupled thermal-hydraulic-mechanical simulation of the initiation and extension of hydraulic fractures.It also calculates the mass/heat transport of injected hydraulic fluids as well as proppants.Our simulator is able to handle anisotropic formations with multiple layers.Our simulator has been validated by comparing with an analytical solution as well as Ribeiro and Sharma model.Our model can simulate fracture spacing effect on fracture profile when combining IFD with Discontinuous Displacement Method(DDM).