对传统直角坐标网格进行改进,根据计算机图形学对与物面相交的正方形网格进行切割,使单元与物面贴体,并能根据流场的变化对网格区域进行局部加密。结合直角坐标网格与非结构网格的特点,编写了自适应直角坐标网格的生成程序,并发展了基...对传统直角坐标网格进行改进,根据计算机图形学对与物面相交的正方形网格进行切割,使单元与物面贴体,并能根据流场的变化对网格区域进行局部加密。结合直角坐标网格与非结构网格的特点,编写了自适应直角坐标网格的生成程序,并发展了基于切割单元直角坐标网格的直接模拟蒙特卡洛(Direct simulation Monte Carlo,DSMC)方法。同时通过当地模拟分子数、动态时间步长技术优化DSMC方法,提升计算效率。数值结果表明,本方法在计算过程中使搜索效率得到大大提升,同时也保证了物面网格单元的贴体性,提高了物面附近流场的计算精度,还易于实现网格自适应。展开更多
In this paper,we propose an immersed finite volume element method for solving the semi-linear elliptic interface problems with non-homogeneous jump conditions.Furthermore,two-grid techniques are used to improve the co...In this paper,we propose an immersed finite volume element method for solving the semi-linear elliptic interface problems with non-homogeneous jump conditions.Furthermore,two-grid techniques are used to improve the computational efficiency.In this way,we only need to solve a non-linear system on the coarse grid,and a linear system on the fine grid.Numerical results illustrate that the proposed method can solve the semi-linear elliptic interface problems efficiently.Approximate secondorder accuracy for the solution in the L¥norm can be obtained for the considered examples.展开更多
A new mesh free method proposed by the authors was presented, in which the derivatives at each node were constructed using whole derivative formulas through the nodes selected around the node using local Cartesian fra...A new mesh free method proposed by the authors was presented, in which the derivatives at each node were constructed using whole derivative formulas through the nodes selected around the node using local Cartesian frame in an autonomous manner, so that without any element it could be considered as a completely mesh free method. The method was tested with a numerical example, and reliable solution was obtained with high accuracy and efficiency.展开更多
Performance of the LSFD method is compared with conventional FD schemes. Generally, 9-point stencils for 2D cases and 27-point stencils for 3D cases are used for the approximation of the first and second order derivat...Performance of the LSFD method is compared with conventional FD schemes. Generally, 9-point stencils for 2D cases and 27-point stencils for 3D cases are used for the approximation of the first and second order derivatives obtained with conventional central difference schemes. When the same stencils are used, explicit LSFD formulations for approximation of the first and second order derivatives are presented. The LSFD formulations are actually a combination of conventional central difference schemes along relevant mesh lines. It has been found that LSFD formulations need much less iteration steps than the conventional FD schemes to converge, and the ratio of mesh spacing in the x and y directions is an important parameter in the LSFD application, with a great impact on stability of LSFD computation.展开更多
A type of mesh called a body-fi tted Cartesian mesh,very different from the traditional structured body-fi tted mesh,is established.At f irst,the right parallelepiped mesh is generated,then,a feature analysis is done ...A type of mesh called a body-fi tted Cartesian mesh,very different from the traditional structured body-fi tted mesh,is established.At f irst,the right parallelepiped mesh is generated,then,a feature analysis is done on the cross sections.These cross sections are the intersections of the casting shape with the right parallelepiped grids(under the Cartesian coordinate system).On the basis of the feature analysis,two sorts of body-f itted boundary grids,shape-keeping grids and shape-distortion grids,are def ined.Shape-distortion grids can be removed or weaken by increasing the number of grids or moving the coordinates of the mesh generation region,so actually the body-fitted Cartesian mesh generation is to get shape-keeping grids.A shape-keeping grid mainly consists of two sorts of surfaces(I type face and II type face),and each of them is joined by two types of points(I type point and II type point).If only these two types of points were given,the shape-keeping mesh would be constructed.In this paper,the cases of the above two boundary grids being generated were discussed.An algorithm was put forward to get the shape-keeping grids.Several body-fi tted Cartesian meshes generated on castings show the validity of the algorithm.The mesh generation examples show that the body-fi tted Cartesian mesh is more excellent than the right parallelepiped mesh in aspects of decreasing grids number and being closer to the shape of the casting solid.展开更多
In the paper, we study a high order numerical boundary scheme for solving the complex moving boundary problem on a fixed Cartesian mesh, and numerically investigate the moving rigid body with the complex boundary unde...In the paper, we study a high order numerical boundary scheme for solving the complex moving boundary problem on a fixed Cartesian mesh, and numerically investigate the moving rigid body with the complex boundary under the impingement of an inviscid shock wave. Based on the high order inverse Lax-Wendroff(ILW) procedure developed in the previous work(TAN, S. and SHU, C. W. A high order moving boundary treatment for compressible inviscid flows. Journal of Computational Physics, 230(15),6023–6036(2011)), in which the authors only considered the translation of the rigid body,we consider both translation and rotation of the body in this paper. In particular, we reformulate the material derivative on the moving boundary with no-penetration condition, and the newly obtained formula plays a key role in the proposed algorithm. Several numerical examples, including cylinder, elliptic cylinder, and NACA0012 airfoil, are given to indicate the effectiveness and robustness of the present method.展开更多
文摘对传统直角坐标网格进行改进,根据计算机图形学对与物面相交的正方形网格进行切割,使单元与物面贴体,并能根据流场的变化对网格区域进行局部加密。结合直角坐标网格与非结构网格的特点,编写了自适应直角坐标网格的生成程序,并发展了基于切割单元直角坐标网格的直接模拟蒙特卡洛(Direct simulation Monte Carlo,DSMC)方法。同时通过当地模拟分子数、动态时间步长技术优化DSMC方法,提升计算效率。数值结果表明,本方法在计算过程中使搜索效率得到大大提升,同时也保证了物面网格单元的贴体性,提高了物面附近流场的计算精度,还易于实现网格自适应。
基金supported by the National Natural Science Foundation of China Nos.11701283,12171482the Fundamental Research Funds for the Central Universities No.KJQN201839Science Foundation of China University of Petroleum(Beijing)No.2462020XKJS02.
文摘In this paper,we propose an immersed finite volume element method for solving the semi-linear elliptic interface problems with non-homogeneous jump conditions.Furthermore,two-grid techniques are used to improve the computational efficiency.In this way,we only need to solve a non-linear system on the coarse grid,and a linear system on the fine grid.Numerical results illustrate that the proposed method can solve the semi-linear elliptic interface problems efficiently.Approximate secondorder accuracy for the solution in the L¥norm can be obtained for the considered examples.
基金Project supported by the Nation-al Natural Science Foundation of China(No.10372055)the Shanghai Leading Academic Discipline Project(No.Y0103)
文摘A new mesh free method proposed by the authors was presented, in which the derivatives at each node were constructed using whole derivative formulas through the nodes selected around the node using local Cartesian frame in an autonomous manner, so that without any element it could be considered as a completely mesh free method. The method was tested with a numerical example, and reliable solution was obtained with high accuracy and efficiency.
基金supported by the National Natural Science Foundation of China (Nos. 10872005, 10532010)
文摘Performance of the LSFD method is compared with conventional FD schemes. Generally, 9-point stencils for 2D cases and 27-point stencils for 3D cases are used for the approximation of the first and second order derivatives obtained with conventional central difference schemes. When the same stencils are used, explicit LSFD formulations for approximation of the first and second order derivatives are presented. The LSFD formulations are actually a combination of conventional central difference schemes along relevant mesh lines. It has been found that LSFD formulations need much less iteration steps than the conventional FD schemes to converge, and the ratio of mesh spacing in the x and y directions is an important parameter in the LSFD application, with a great impact on stability of LSFD computation.
基金financially supported by the fund of the State Key Laboratory of Solidif ication Processing in NWPU(No:SKLSP201006)the fund of the National Basic Research Program of China(No:2011CB610402)
文摘A type of mesh called a body-fi tted Cartesian mesh,very different from the traditional structured body-fi tted mesh,is established.At f irst,the right parallelepiped mesh is generated,then,a feature analysis is done on the cross sections.These cross sections are the intersections of the casting shape with the right parallelepiped grids(under the Cartesian coordinate system).On the basis of the feature analysis,two sorts of body-f itted boundary grids,shape-keeping grids and shape-distortion grids,are def ined.Shape-distortion grids can be removed or weaken by increasing the number of grids or moving the coordinates of the mesh generation region,so actually the body-fitted Cartesian mesh generation is to get shape-keeping grids.A shape-keeping grid mainly consists of two sorts of surfaces(I type face and II type face),and each of them is joined by two types of points(I type point and II type point).If only these two types of points were given,the shape-keeping mesh would be constructed.In this paper,the cases of the above two boundary grids being generated were discussed.An algorithm was put forward to get the shape-keeping grids.Several body-fi tted Cartesian meshes generated on castings show the validity of the algorithm.The mesh generation examples show that the body-fi tted Cartesian mesh is more excellent than the right parallelepiped mesh in aspects of decreasing grids number and being closer to the shape of the casting solid.
基金Project supported by the National Natural Science Foundation of China (Nos. 11901555, 11901213,11871448, and 11732016)the National Numerical Windtunnel Project (No. NNW2019ZT4-B10)。
文摘In the paper, we study a high order numerical boundary scheme for solving the complex moving boundary problem on a fixed Cartesian mesh, and numerically investigate the moving rigid body with the complex boundary under the impingement of an inviscid shock wave. Based on the high order inverse Lax-Wendroff(ILW) procedure developed in the previous work(TAN, S. and SHU, C. W. A high order moving boundary treatment for compressible inviscid flows. Journal of Computational Physics, 230(15),6023–6036(2011)), in which the authors only considered the translation of the rigid body,we consider both translation and rotation of the body in this paper. In particular, we reformulate the material derivative on the moving boundary with no-penetration condition, and the newly obtained formula plays a key role in the proposed algorithm. Several numerical examples, including cylinder, elliptic cylinder, and NACA0012 airfoil, are given to indicate the effectiveness and robustness of the present method.
