This paper considers the finite difference(FD)approximations of diffusion operators and the boundary treatments for different boundary conditions.The proposed schemes have the compact form and could achieve arbitrary ...This paper considers the finite difference(FD)approximations of diffusion operators and the boundary treatments for different boundary conditions.The proposed schemes have the compact form and could achieve arbitrary even order of accuracy.The main idea is to make use of the lower order compact schemes recursively,so as to obtain the high order compact schemes formally.Moreover,the schemes can be implemented efficiently by solving a series of tridiagonal systems recursively or the fast Fourier transform(FFT).With mathematical induction,the eigenvalues of the proposed differencing operators are shown to be bounded away from zero,which indicates the positive definiteness of the operators.To obtain numerical boundary conditions for the high order schemes,the simplified inverse Lax-Wendroff(SILW)procedure is adopted and the stability analysis is performed by the Godunov-Ryabenkii method and the eigenvalue spectrum visualization method.Various numerical experiments are provided to demonstrate the effectiveness and robustness of our algorithms.展开更多
In order to control traffic congestion, many mathematical models have been used for several decades. In this paper, we study diffusion-type traffic flow model based on exponential velocity density relation, which prov...In order to control traffic congestion, many mathematical models have been used for several decades. In this paper, we study diffusion-type traffic flow model based on exponential velocity density relation, which provides a non-linear second-order parabolic partial differential equation. The analytical solution of the diffusion-type traffic flow model is very complicated to approximate the initial density of the Cauchy problem as a function of x from given data and it may cause a huge error. For the complexity of the analytical solution, the numerical solution is performed by implementing an explicit upwind, explicitly centered, and second-order Lax-Wendroff scheme for the numerical solution. From the comparison of relative error among these three schemes, it is observed that Lax-Wendroff scheme gives less error than the explicit upwind and explicit centered difference scheme. The numerical, analytical analysis and comparative result discussion bring out the fact that the Lax-Wendroff scheme with exponential velocity-density relation of diffusion type traffic flow model is suitable for the congested area and shows a better fit in traffic-congested regions.展开更多
We give a brief discussion of some of the contributions of Peter Lax to Com- putational Fluid Dynamics. These include the Lax-Friedrichs and Lax-Wendroff numerical schemes. We also mention his collaboration in the 198...We give a brief discussion of some of the contributions of Peter Lax to Com- putational Fluid Dynamics. These include the Lax-Friedrichs and Lax-Wendroff numerical schemes. We also mention his collaboration in the 1983 HLL Riemann solver. We de- velop two-dimensional Lax-Friedrichs and Lax-Wendroff schemes for the Lagrangian form of the Euler equations on triangular grids. We apply a composite scheme that uses a Lax- Friedrichs time step as a dissipative filter after several Lax-Wendroff time steps. Numerical results for Noh＇s infinite strength shock problem, the Sedov blast wave problem, and the Saltzman piston problem are presented.展开更多
This paper studies the geometric boundary representations for Inverse Lax-Wendroff(ILW)method,aiming to develop a practical computer-aided engineering method without body-fitted meshes.We propose the signed distance f...This paper studies the geometric boundary representations for Inverse Lax-Wendroff(ILW)method,aiming to develop a practical computer-aided engineering method without body-fitted meshes.We propose the signed distance function(SDF)representation of the geometric boundary and design an extremely efficient algorithm for foot point calculation,which is particularly in line with the needs of ILW.Theoretical and numerical analyses demonstrate that the SDF representation of geometric boundary can satisfy ILW’s needs better than others.The effectiveness and robustness of our proposed method are verified by simulating initial boundary value computational physical problems of Euler equation for compressible fluids.展开更多
In this paper, we use Hermite weighted essentially non-oscillatory （HWENO） schemes with a Lax-Wendroff time discretization procedure, termed HWENO-LW schemes, to solve Hamilton-Jacobi equations. The idea of the reco...In this paper, we use Hermite weighted essentially non-oscillatory （HWENO） schemes with a Lax-Wendroff time discretization procedure, termed HWENO-LW schemes, to solve Hamilton-Jacobi equations. The idea of the reconstruction in the HWENO schemes comes from the original WENO schemes, however both the function and its first derivative values are evolved in time and are used in the reconstruction. One major advantage of HWENO schemes is its compactness in the reconstruction. We explore the possibility in avoiding the nonlinear weights for part of the procedure, hence reducing the cost but still maintaining non-oscillatory properties for problems with strong discontinuous derivative. As a result, comparing with HWENO with Runge-Kutta time discretizations schemes （HWENO-RK） of Qiu and Shu [19] for Hamilton-Jacobi equations, the major advantages of HWENO-LW schemes are their saving of computational cost and their compactness in the reconstruction. Extensive numerical experiments are performed to illustrate the capability of the method.展开更多
IT is known from Brenner, Thomee and Wahlbin that the well-known second-order Lax-Wendroff scheme is stable in L^2, but unstable in L^p, p≠2. Generally speaking, if the initialdata is smooth enough and if a differenc...IT is known from Brenner, Thomee and Wahlbin that the well-known second-order Lax-Wendroff scheme is stable in L^2, but unstable in L^p, p≠2. Generally speaking, if the initialdata is smooth enough and if a difference scheme, which is stable in L^p for some p, has orderof accuracy μ, then we can expect that the solution of the difference scheme converges to thesolution of the differential equation at the rate of order μ in L^p. But for discontinuous solu-tions, which are essential to hyperbolic equations, the above expectation is not true. Error es-timates for discontinuous solutions not only have theoretical meaning, but also practical value.展开更多
The inverse Lax-Wendroff(ILW)procedure is a numerical boundary treatment technique,which allows finite difference schemes and other schemes to achieve stability and high order accuracy when using cartesian meshes to s...The inverse Lax-Wendroff(ILW)procedure is a numerical boundary treatment technique,which allows finite difference schemes and other schemes to achieve stability and high order accuracy when using cartesian meshes to solve boundary value problems defined on complex computational domain.In this short survey we summarize the main ingredients of the ILW procedure,discuss its applicability and stability properties,and provide possible directions of its future development.展开更多
基金supported by the NSFC grant 11801143J.Lu’s research is partially supported by the NSFC grant 11901213+3 种基金the National Key Research and Development Program of China grant 2021YFA1002900supported by the NSFC grant 11801140,12171177the Young Elite Scientists Sponsorship Program by Henan Association for Science and Technology of China grant 2022HYTP0009the Program for Young Key Teacher of Henan Province of China grant 2021GGJS067.
文摘This paper considers the finite difference(FD)approximations of diffusion operators and the boundary treatments for different boundary conditions.The proposed schemes have the compact form and could achieve arbitrary even order of accuracy.The main idea is to make use of the lower order compact schemes recursively,so as to obtain the high order compact schemes formally.Moreover,the schemes can be implemented efficiently by solving a series of tridiagonal systems recursively or the fast Fourier transform(FFT).With mathematical induction,the eigenvalues of the proposed differencing operators are shown to be bounded away from zero,which indicates the positive definiteness of the operators.To obtain numerical boundary conditions for the high order schemes,the simplified inverse Lax-Wendroff(SILW)procedure is adopted and the stability analysis is performed by the Godunov-Ryabenkii method and the eigenvalue spectrum visualization method.Various numerical experiments are provided to demonstrate the effectiveness and robustness of our algorithms.
文摘In order to control traffic congestion, many mathematical models have been used for several decades. In this paper, we study diffusion-type traffic flow model based on exponential velocity density relation, which provides a non-linear second-order parabolic partial differential equation. The analytical solution of the diffusion-type traffic flow model is very complicated to approximate the initial density of the Cauchy problem as a function of x from given data and it may cause a huge error. For the complexity of the analytical solution, the numerical solution is performed by implementing an explicit upwind, explicitly centered, and second-order Lax-Wendroff scheme for the numerical solution. From the comparison of relative error among these three schemes, it is observed that Lax-Wendroff scheme gives less error than the explicit upwind and explicit centered difference scheme. The numerical, analytical analysis and comparative result discussion bring out the fact that the Lax-Wendroff scheme with exponential velocity-density relation of diffusion type traffic flow model is suitable for the congested area and shows a better fit in traffic-congested regions.
基金performed under the auspices of the National Nuclear Security Administration of the US Department of Energy at Los Alamos National Laboratory under Contract No.DE-AC52-06NA25396supported in part by the Czech Science Foundation GrantP205/10/0814the Czech Ministry of Education grants MSM 6840770022 and LC528
文摘We give a brief discussion of some of the contributions of Peter Lax to Com- putational Fluid Dynamics. These include the Lax-Friedrichs and Lax-Wendroff numerical schemes. We also mention his collaboration in the 1983 HLL Riemann solver. We de- velop two-dimensional Lax-Friedrichs and Lax-Wendroff schemes for the Lagrangian form of the Euler equations on triangular grids. We apply a composite scheme that uses a Lax- Friedrichs time step as a dissipative filter after several Lax-Wendroff time steps. Numerical results for Noh＇s infinite strength shock problem, the Sedov blast wave problem, and the Saltzman piston problem are presented.
文摘This paper studies the geometric boundary representations for Inverse Lax-Wendroff(ILW)method,aiming to develop a practical computer-aided engineering method without body-fitted meshes.We propose the signed distance function(SDF)representation of the geometric boundary and design an extremely efficient algorithm for foot point calculation,which is particularly in line with the needs of ILW.Theoretical and numerical analyses demonstrate that the SDF representation of geometric boundary can satisfy ILW’s needs better than others.The effectiveness and robustness of our proposed method are verified by simulating initial boundary value computational physical problems of Euler equation for compressible fluids.
基金Research partially supported by NNSFC grant 10371118,SRF for ROCS,SEM and Nanjing University Talent Development Foundation.
文摘In this paper, we use Hermite weighted essentially non-oscillatory （HWENO） schemes with a Lax-Wendroff time discretization procedure, termed HWENO-LW schemes, to solve Hamilton-Jacobi equations. The idea of the reconstruction in the HWENO schemes comes from the original WENO schemes, however both the function and its first derivative values are evolved in time and are used in the reconstruction. One major advantage of HWENO schemes is its compactness in the reconstruction. We explore the possibility in avoiding the nonlinear weights for part of the procedure, hence reducing the cost but still maintaining non-oscillatory properties for problems with strong discontinuous derivative. As a result, comparing with HWENO with Runge-Kutta time discretizations schemes （HWENO-RK） of Qiu and Shu [19] for Hamilton-Jacobi equations, the major advantages of HWENO-LW schemes are their saving of computational cost and their compactness in the reconstruction. Extensive numerical experiments are performed to illustrate the capability of the method.
文摘IT is known from Brenner, Thomee and Wahlbin that the well-known second-order Lax-Wendroff scheme is stable in L^2, but unstable in L^p, p≠2. Generally speaking, if the initialdata is smooth enough and if a difference scheme, which is stable in L^p for some p, has orderof accuracy μ, then we can expect that the solution of the difference scheme converges to thesolution of the differential equation at the rate of order μ in L^p. But for discontinuous solu-tions, which are essential to hyperbolic equations, the above expectation is not true. Error es-timates for discontinuous solutions not only have theoretical meaning, but also practical value.
文摘The inverse Lax-Wendroff(ILW)procedure is a numerical boundary treatment technique,which allows finite difference schemes and other schemes to achieve stability and high order accuracy when using cartesian meshes to solve boundary value problems defined on complex computational domain.In this short survey we summarize the main ingredients of the ILW procedure,discuss its applicability and stability properties,and provide possible directions of its future development.
