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.展开更多
We study the convergence and asymptotic compatibility of higher order collocation methods for nonlocal operators inspired by peridynamics,a nonlocal formulation of continuum mechanics.We prove that the methods are opt...We study the convergence and asymptotic compatibility of higher order collocation methods for nonlocal operators inspired by peridynamics,a nonlocal formulation of continuum mechanics.We prove that the methods are optimally convergent with respect to the polynomial degree of the approximation.A numerical method is said to be asymptotically compatible if the sequence of approximate solutions of the nonlocal problem converges to the solution of the corresponding local problem as the horizon and the grid sizes simultaneously approach zero.We carry out a calibration process via Taylor series expansions and a scaling of the nonlocal operator via a strain energy density argument to ensure that the resulting collocation methods are asymptotically compatible.We fnd that,for polynomial degrees greater than or equal to two,there exists a calibration constant independent of the horizon size and the grid size such that the resulting collocation methods for the nonlocal difusion are asymptotically compatible.We verify these fndings through extensive numerical experiments.展开更多
In this work,the localized method of fundamental solution(LMFS)is extended to Signorini problem.Unlike the traditional fundamental solution(MFS),the LMFS approximates the field quantity at each node by using the field...In this work,the localized method of fundamental solution(LMFS)is extended to Signorini problem.Unlike the traditional fundamental solution(MFS),the LMFS approximates the field quantity at each node by using the field quantities at the adjacent nodes.The idea of the LMFS is similar to the localized domain type method.The fictitious boundary nodes are proposed to impose the boundary condition and governing equations at each node to formulate a sparse matrix.The inequality boundary condition of Signorini problem is solved indirectly by introducing nonlinear complementarity problem function(NCP-function).Numerical examples are carried out to validate the reliability and effectiveness of the LMFS in solving Signorini problems.展开更多
Under some local superquadratic conditions on <em>W</em> (<em>t</em>, <em>u</em>) with respect to <em>u</em>, the existence of infinitely many solutions is obtained for ...Under some local superquadratic conditions on <em>W</em> (<em>t</em>, <em>u</em>) with respect to <em>u</em>, the existence of infinitely many solutions is obtained for the nonperiodic fractional Hamiltonian systems<img src="Edit_b2a2ac0a-6dde-474f-8c75-e9f5fc7b9918.bmp" alt="" />, where <em>L</em> (<em>t</em>) is unnecessarily coercive.展开更多
Traditional geoacoustic inversions are generally solved by matched-field processing in combination with metaheuristic global searching algorithms which usually need massive computations. This paper proposes a new phys...Traditional geoacoustic inversions are generally solved by matched-field processing in combination with metaheuristic global searching algorithms which usually need massive computations. This paper proposes a new physical framework for geoacoustic retrievals. A parabolic approximation of wave equation with non-local boundary condition is used as the forward propagation model. The expressions of the corresponding tangent linear model and the adjoint operator are derived, respectively, by variational method. The analytical expressions for the gradient of the cost function with respect to the control variables can be formulated by the adjoint operator, which in turn can be used for optimization by the gradient-based method.展开更多
A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements f...A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.展开更多
基金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.
文摘We study the convergence and asymptotic compatibility of higher order collocation methods for nonlocal operators inspired by peridynamics,a nonlocal formulation of continuum mechanics.We prove that the methods are optimally convergent with respect to the polynomial degree of the approximation.A numerical method is said to be asymptotically compatible if the sequence of approximate solutions of the nonlocal problem converges to the solution of the corresponding local problem as the horizon and the grid sizes simultaneously approach zero.We carry out a calibration process via Taylor series expansions and a scaling of the nonlocal operator via a strain energy density argument to ensure that the resulting collocation methods are asymptotically compatible.We fnd that,for polynomial degrees greater than or equal to two,there exists a calibration constant independent of the horizon size and the grid size such that the resulting collocation methods for the nonlocal difusion are asymptotically compatible.We verify these fndings through extensive numerical experiments.
基金supported by the National Science Foundation of China(No.52109089)support of Post Doctor Program(2019M652281)Nature Science Foundation of Jiangxi Province(20192BAB216040).
文摘In this work,the localized method of fundamental solution(LMFS)is extended to Signorini problem.Unlike the traditional fundamental solution(MFS),the LMFS approximates the field quantity at each node by using the field quantities at the adjacent nodes.The idea of the LMFS is similar to the localized domain type method.The fictitious boundary nodes are proposed to impose the boundary condition and governing equations at each node to formulate a sparse matrix.The inequality boundary condition of Signorini problem is solved indirectly by introducing nonlinear complementarity problem function(NCP-function).Numerical examples are carried out to validate the reliability and effectiveness of the LMFS in solving Signorini problems.
文摘Under some local superquadratic conditions on <em>W</em> (<em>t</em>, <em>u</em>) with respect to <em>u</em>, the existence of infinitely many solutions is obtained for the nonperiodic fractional Hamiltonian systems<img src="Edit_b2a2ac0a-6dde-474f-8c75-e9f5fc7b9918.bmp" alt="" />, where <em>L</em> (<em>t</em>) is unnecessarily coercive.
基金Project supported by the Foundation of Key Laboratory of Marine Intelligent Equipment and System of Ministry of Education,China(Grant No.SJTUMIES1908)the National Natural Science Foundation of China(Grant No.41775027)
文摘Traditional geoacoustic inversions are generally solved by matched-field processing in combination with metaheuristic global searching algorithms which usually need massive computations. This paper proposes a new physical framework for geoacoustic retrievals. A parabolic approximation of wave equation with non-local boundary condition is used as the forward propagation model. The expressions of the corresponding tangent linear model and the adjoint operator are derived, respectively, by variational method. The analytical expressions for the gradient of the cost function with respect to the control variables can be formulated by the adjoint operator, which in turn can be used for optimization by the gradient-based method.
基金supported by the National Natural Science Foundation of China(Nos.11271273 and 11271298)
文摘A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.
