We introduce adaptive moving mesh central-upwind schemes for one-and two-dimensional hyperbolic systems of conservation and balance laws.The proposed methods consist of three steps.First,the solution is evolved by sol...We introduce adaptive moving mesh central-upwind schemes for one-and two-dimensional hyperbolic systems of conservation and balance laws.The proposed methods consist of three steps.First,the solution is evolved by solving the studied system by the second-order semi-discrete central-upwind scheme on either the one-dimensional nonuniform grid or the two-dimensional structured quadrilateral mesh.When the evolution step is complete,the grid points are redistributed according to the moving mesh differential equation.Finally,the evolved solution is projected onto the new mesh in a conservative manner.The resulting adaptive moving mesh methods are applied to the one-and two-dimensional Euler equations of gas dynamics and granular hydrodynamics systems.Our numerical results demonstrate that in both cases,the adaptive moving mesh central-upwind schemes outperform their uniform mesh counterparts.展开更多
Linear Least Squares(LLS) problems are particularly difficult to solve because they are frequently ill-conditioned, and involve large quantities of data. Ill-conditioned LLS problems are commonly seen in mathematics...Linear Least Squares(LLS) problems are particularly difficult to solve because they are frequently ill-conditioned, and involve large quantities of data. Ill-conditioned LLS problems are commonly seen in mathematics and geosciences, where regularization algorithms are employed to seek optimal solutions. For many problems, even with the use of regularization algorithms it may be impossible to obtain an accurate solution. Riley and Golub suggested an iterative scheme for solving LLS problems. For the early iteration algorithm, it is difficult to improve the well-conditioned perturbed matrix and accelerate the convergence at the same time. Aiming at this problem, self-adaptive iteration algorithm(SAIA) is proposed in this paper for solving severe ill-conditioned LLS problems. The algorithm is different from other popular algorithms proposed in recent references. It avoids matrix inverse by using Cholesky decomposition, and tunes the perturbation parameter according to the rate of residual error decline in the iterative process. Example shows that the algorithm can greatly reduce iteration times, accelerate the convergence,and also greatly enhance the computation accuracy.展开更多
To improve the efficiency of the discrete unified gas kinetic scheme(DUGKS)in capturing cross-scale flow physics,an adaptive partitioning-based discrete unified gas kinetic scheme(ADUGKS)is developed in this work.The ...To improve the efficiency of the discrete unified gas kinetic scheme(DUGKS)in capturing cross-scale flow physics,an adaptive partitioning-based discrete unified gas kinetic scheme(ADUGKS)is developed in this work.The ADUGKS is designed from the discrete characteristic solution to the Boltzmann-BGK equation,which contains the initial distribution function and the local equilibrium state.The initial distribution function contributes to the calculation of free streaming fluxes and the local equilibrium state contributes to the calculation of equilibrium fluxes.When the contribution of the initial distribution function is negative,the local flow field can be regarded as the continuous flow and the Navier-Stokes(N-S)equations can be used to obtain the solution directly.Otherwise,the discrete distribution functions should be updated by the Boltzmann equation to capture the rarefaction effect.Given this,in the ADUGKS,the computational domain is divided into the DUGKS cell and the N-S cell based on the contribu-tion of the initial distribution function to the calculation of free streaming fluxes.In the N-S cell,the local flow field is evolved by solving the N-S equations,while in the DUGKS cell,both the discrete velocity Boltzmann equation and the correspond-ing macroscopic governing equations are solved by a modified DUGKS.Since more and more cells turn into the N-S cell with the decrease of the Knudsen number,a significant acceleration can be achieved for the ADUGKS in the continuum flow regime as compared with the DUGKS.展开更多
基金The work of A.Kurganov was supported in part by the National Natural Science Foundation of China grant 11771201by the fund of the Guangdong Provincial Key Laboratory of Computational Science and Material Design(No.2019B030301001).
文摘We introduce adaptive moving mesh central-upwind schemes for one-and two-dimensional hyperbolic systems of conservation and balance laws.The proposed methods consist of three steps.First,the solution is evolved by solving the studied system by the second-order semi-discrete central-upwind scheme on either the one-dimensional nonuniform grid or the two-dimensional structured quadrilateral mesh.When the evolution step is complete,the grid points are redistributed according to the moving mesh differential equation.Finally,the evolved solution is projected onto the new mesh in a conservative manner.The resulting adaptive moving mesh methods are applied to the one-and two-dimensional Euler equations of gas dynamics and granular hydrodynamics systems.Our numerical results demonstrate that in both cases,the adaptive moving mesh central-upwind schemes outperform their uniform mesh counterparts.
基金supported by Open Fund of Engineering Laboratory of Spatial Information Technology of Highway Geological Disaster Early Warning in Hunan Province(Changsha University of Science&Technology,kfj150602)Hunan Province Science and Technology Program Funded Projects,China(2015NK3035)+1 种基金the Land and Resources Department Scientific Research Project of Hunan Province,China(2013-27)the Education Department Scientific Research Project of Hunan Province,China(13C1011)
文摘Linear Least Squares(LLS) problems are particularly difficult to solve because they are frequently ill-conditioned, and involve large quantities of data. Ill-conditioned LLS problems are commonly seen in mathematics and geosciences, where regularization algorithms are employed to seek optimal solutions. For many problems, even with the use of regularization algorithms it may be impossible to obtain an accurate solution. Riley and Golub suggested an iterative scheme for solving LLS problems. For the early iteration algorithm, it is difficult to improve the well-conditioned perturbed matrix and accelerate the convergence at the same time. Aiming at this problem, self-adaptive iteration algorithm(SAIA) is proposed in this paper for solving severe ill-conditioned LLS problems. The algorithm is different from other popular algorithms proposed in recent references. It avoids matrix inverse by using Cholesky decomposition, and tunes the perturbation parameter according to the rate of residual error decline in the iterative process. Example shows that the algorithm can greatly reduce iteration times, accelerate the convergence,and also greatly enhance the computation accuracy.
基金the National Natural Science Foundation of China(12202191,92271103)Natural Science Foundation of Jiangsu Province(BK20210273)+1 种基金Fund of Prospective Layout of Scientific Research for NUAA(Nanjing University of Aeronautics and Astronautics)Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD).
文摘To improve the efficiency of the discrete unified gas kinetic scheme(DUGKS)in capturing cross-scale flow physics,an adaptive partitioning-based discrete unified gas kinetic scheme(ADUGKS)is developed in this work.The ADUGKS is designed from the discrete characteristic solution to the Boltzmann-BGK equation,which contains the initial distribution function and the local equilibrium state.The initial distribution function contributes to the calculation of free streaming fluxes and the local equilibrium state contributes to the calculation of equilibrium fluxes.When the contribution of the initial distribution function is negative,the local flow field can be regarded as the continuous flow and the Navier-Stokes(N-S)equations can be used to obtain the solution directly.Otherwise,the discrete distribution functions should be updated by the Boltzmann equation to capture the rarefaction effect.Given this,in the ADUGKS,the computational domain is divided into the DUGKS cell and the N-S cell based on the contribu-tion of the initial distribution function to the calculation of free streaming fluxes.In the N-S cell,the local flow field is evolved by solving the N-S equations,while in the DUGKS cell,both the discrete velocity Boltzmann equation and the correspond-ing macroscopic governing equations are solved by a modified DUGKS.Since more and more cells turn into the N-S cell with the decrease of the Knudsen number,a significant acceleration can be achieved for the ADUGKS in the continuum flow regime as compared with the DUGKS.