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.展开更多
Fixed-point fast sweeping WENO methods are a class of efficient high-order numerical methods to solve steady-state solutions of hyperbolic partial differential equations(PDEs).The Gauss-Seidel iterations and alternati...Fixed-point fast sweeping WENO methods are a class of efficient high-order numerical methods to solve steady-state solutions of hyperbolic partial differential equations(PDEs).The Gauss-Seidel iterations and alternating sweeping strategy are used to cover characteristics of hyperbolic PDEs in each sweeping order to achieve fast convergence rate to steady-state solutions.A nice property of fixed-point fast sweeping WENO methods which distinguishes them from other fast sweeping methods is that they are explicit and do not require inverse operation of nonlinear local systems.Hence,they are easy to be applied to a general hyperbolic system.To deal with the difficulties associated with numerical boundary treatment when high-order finite difference methods on a Cartesian mesh are used to solve hyperbolic PDEs on complex domains,inverse Lax-Wendroff(ILW)procedures were developed as a very effective approach in the literature.In this paper,we combine a fifthorder fixed-point fast sweeping WENO method with an ILW procedure to solve steadystate solution of hyperbolic conservation laws on complex computing regions.Numerical experiments are performed to test the method in solving various problems including the cases with the physical boundary not aligned with the grids.Numerical results show highorder accuracy and good performance of the method.Furthermore,the method is compared with the popular third-order total variation diminishing Runge-Kutta(TVD-RK3)time-marching method for steady-state computations.Numerical examples show that for most of examples,the fixed-point fast sweeping method saves more than half CPU time costs than TVD-RK3 to converge to steady-state solutions.展开更多
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.展开更多
本文首先概要地介绍了国际与日共存(International Living with a Star, ILWS)计划的目的、组织机构及ILWS第一次和第二次工作会议的情况。中国国家航天局于2003 年4 月正式参加了ILWS计划,并参加了第一次和第二次工作会议。中方代表在...本文首先概要地介绍了国际与日共存(International Living with a Star, ILWS)计划的目的、组织机构及ILWS第一次和第二次工作会议的情况。中国国家航天局于2003 年4 月正式参加了ILWS计划,并参加了第一次和第二次工作会议。中方代表在会上提出了中国参加ILWS的建议,报告了中国新提出的“空间风暴探测计划”。现双星计划已被列为2006年以前的ILWS项目,新提出的空间风暴探测计划已被列为ILWS第一阶段的重要项目。其次介绍了空间风暴探测计划项目的需求,主要目标和研究内容,有效载荷、轨道方案以及空间风暴计划的进度计划等。展开更多
基金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.
基金Research was supported by the NSFC Grant 11872210Research was supported by the NSFC Grant 11872210 and Grant No.MCMS-I-0120G01+1 种基金Research supported in part by the AFOSR Grant FA9550-20-1-0055NSF Grant DMS-2010107.
文摘Fixed-point fast sweeping WENO methods are a class of efficient high-order numerical methods to solve steady-state solutions of hyperbolic partial differential equations(PDEs).The Gauss-Seidel iterations and alternating sweeping strategy are used to cover characteristics of hyperbolic PDEs in each sweeping order to achieve fast convergence rate to steady-state solutions.A nice property of fixed-point fast sweeping WENO methods which distinguishes them from other fast sweeping methods is that they are explicit and do not require inverse operation of nonlinear local systems.Hence,they are easy to be applied to a general hyperbolic system.To deal with the difficulties associated with numerical boundary treatment when high-order finite difference methods on a Cartesian mesh are used to solve hyperbolic PDEs on complex domains,inverse Lax-Wendroff(ILW)procedures were developed as a very effective approach in the literature.In this paper,we combine a fifthorder fixed-point fast sweeping WENO method with an ILW procedure to solve steadystate solution of hyperbolic conservation laws on complex computing regions.Numerical experiments are performed to test the method in solving various problems including the cases with the physical boundary not aligned with the grids.Numerical results show highorder accuracy and good performance of the method.Furthermore,the method is compared with the popular third-order total variation diminishing Runge-Kutta(TVD-RK3)time-marching method for steady-state computations.Numerical examples show that for most of examples,the fixed-point fast sweeping method saves more than half CPU time costs than TVD-RK3 to converge to steady-state solutions.
文摘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.
文摘本文首先概要地介绍了国际与日共存(International Living with a Star, ILWS)计划的目的、组织机构及ILWS第一次和第二次工作会议的情况。中国国家航天局于2003 年4 月正式参加了ILWS计划,并参加了第一次和第二次工作会议。中方代表在会上提出了中国参加ILWS的建议,报告了中国新提出的“空间风暴探测计划”。现双星计划已被列为2006年以前的ILWS项目,新提出的空间风暴探测计划已被列为ILWS第一阶段的重要项目。其次介绍了空间风暴探测计划项目的需求,主要目标和研究内容,有效载荷、轨道方案以及空间风暴计划的进度计划等。