目前的气体动理学统一算法(unified gas kinetic scheme,简称UGKS)在求解高速流动问题时的计算效率,难以满足求解复杂工程问题的需求.为了提高该算法的计算效率,本文对模型方程的对流项和碰撞项进行了隐式处理,并针对UGKS界面通量与演...目前的气体动理学统一算法(unified gas kinetic scheme,简称UGKS)在求解高速流动问题时的计算效率,难以满足求解复杂工程问题的需求.为了提高该算法的计算效率,本文对模型方程的对流项和碰撞项进行了隐式处理,并针对UGKS界面通量与演化时间相关的特点,引入了演化时间平均界面通量,通过对控制方程矩阵进行近似LU分解(lower-upper decomposition),实现了隐式UGKS.不同来流马赫数的圆柱绕流算例测试表明,只要演化时间选取得当,隐式方法可以得到与显式方法完全相同的结果,且计算效率可以提高1～2个量级.展开更多
辐射输运方程的数值模拟在天体物理、武器物理和惯性约束与磁约束聚变等研究中都起着非常重要的作用.在实际问题中,背景介质的不透明度系数决定了辐射光子在其中的传输行为.光性薄(不透明度系数小)的介质对辐射光子是透明的,光子与背景...辐射输运方程的数值模拟在天体物理、武器物理和惯性约束与磁约束聚变等研究中都起着非常重要的作用.在实际问题中,背景介质的不透明度系数决定了辐射光子在其中的传输行为.光性薄(不透明度系数小)的介质对辐射光子是透明的,光子与背景介质的相互作用弱,光子传输具有输运传播性质;而光性厚(不透明度系数大)的介质对辐射光子是不透明的,光子与背景介质的相互作用强,光子传输具有扩散性质.因此在辐射输运方程的计算中,如何设计既能得到光子输运传播性质又能捕捉光子扩散传播性质的渐近保持离散格式是目前一个非常活跃和前沿的研究方向.本文简要介绍近几年在辐射输运方程的渐近保持统一气体动理学格式(unified gas kinetic scheme,UGKS)研究方面的进展.本文主要以灰体辐射输运方程为例,详细介绍UGKS的构造方法并给出其渐近分析.同时,结合角度有限元方法和球谐函数展开的方法,介绍如何减弱/去除基于离散纵标法的UGKS具有射线效应的问题,以及相应的改进渐近保持格式.此外,也介绍了将渐近保持的UGKS应用拓展到考虑流体运动的完全辐射流体力学方程组.最后,用一些数值例子验证了格式的渐近保持性和保正性等性质.展开更多
In this paper,we consider the multi-dimensional asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations on distorted quadrilateral meshes.Different from the former scheme [J.Comput.Phys....In this paper,we consider the multi-dimensional asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations on distorted quadrilateral meshes.Different from the former scheme [J.Comput.Phys.285（2015）,265-279] on uniform meshes,in this paper,in order to obtain the boundary fluxes based on the framework of unified gas kinetic scheme（UGKS）,we use the real multi-dimensional reconstruction for the initial data and the macro-terms in the equation of the gray transfer equations.We can prove that the scheme is asymptotic preserving,and especially for the distorted quadrilateral meshes,a nine-point scheme [SIAM J.SCI.COMPUT.30（2008）,1341-1361] for the diffusion limit equations is obtained,which is naturally reduced to standard five-point scheme for the orthogonal meshes.The numerical examples on distorted meshes are included to validate the current approach.展开更多
工程应用中的介质热辐射问题是典型的多尺度问题.基于Boltzmann输运方程建立的各类气体动理学格式,在多尺度瞬态问题中得到了广泛应用.为了克服显式求解方案中CFL条件等的限制,文章通过气体动理学格式实现稳态辐射输运方程的直接求解.SD...工程应用中的介质热辐射问题是典型的多尺度问题.基于Boltzmann输运方程建立的各类气体动理学格式,在多尺度瞬态问题中得到了广泛应用.为了克服显式求解方案中CFL条件等的限制,文章通过气体动理学格式实现稳态辐射输运方程的直接求解.SDUGKS格式由离散统一气体动理学格式(discrete unified gas kinetic scheme,DUGKS)的核心思想发展而来,应用于稳态问题计算.将SDUGKS格式进一步拓展到多尺度的稳态热辐射输运计算.SDUGKS格式继承了DUGKS格式沿特征线离散实现的界面重构,并通过隐式增量格式的单元更新实现对辐射强度的较正,采用逐次迭代法将辐射强度渐近收敛到稳定值.选用多组一维和二维不同尺度的辐射传热算例,通过与特定的解析解以及其他数值方法结果对比,检验了SDUGKS的计算精度和计算效率,并论证了它在多尺度问题中的渐进保持性质.展开更多
Sound wave propagation in rarefied monatomic gases is simulated using a newly developed unified gaskinetic scheme （UGKS）. The numerical calculations are carfled out for a wide range of wave oscillating frequencies. ...Sound wave propagation in rarefied monatomic gases is simulated using a newly developed unified gaskinetic scheme （UGKS）. The numerical calculations are carfled out for a wide range of wave oscillating frequencies. The corresponding rarefaction parameter is defined as the ratio of sound wave frequency to the intermolecular particle collision frequency. The simulation covers the flow regime from the continuum to free molecule one. The treatment of the os- cillating wall boundary condition and the methods for eval- uating the absorption coefficient and sound wave speed are presented in detail. The simulation results from the UGKS are compared to the Navier-Stokes solutions, the direct sim- ulation Monte Carlo （DSMC） simulation, and experimental measurements. Good agreement with the experimental data has been obtained in the whole flow regimes for the corresponding Knudsen number from 0.08 to 32. The cur- rent study clearly demonstrates the capability of the UGKS method in capturing the sound wave propagation and its usefulness for the rarefied flow study.展开更多
Compressible flows exhibit a diverse set of behaviors, where individual particle transports and their collective dynamics play different roles at different scales. At the same time, the atmosphere is composed of diffe...Compressible flows exhibit a diverse set of behaviors, where individual particle transports and their collective dynamics play different roles at different scales. At the same time, the atmosphere is composed of different components that require additional degrees of freedom for representation in computational fluid dynamics. It is challenging to construct an accurate and efficient numerical algorithm to faithfully represent multiscale flow physics across different regimes. In this paper, a unified gas-kinetic scheme(UGKS) is developed to study non-equilibrium multicomponent gaseous flows. Based on the Boltzmann kinetic equation, an analytical space-time evolving solution is used to construct the discretized equations of gas dynamics directly according to cell size and scales of time steps, i.e., the so-called direct modeling method. With the variation in the ratio of the numerical time step to the local particle collision time(or the cell size to the local particle mean free path), the UGKS automatically recovers all scale-dependent flows over the given domain and provides a continuous spectrum of the gas dynamics. The performance of the proposed unified scheme is fully validated through numerical experiments.The UGKS can be a valuable tool to study multiscale and multicomponent flow physics.展开更多
The unified gas-kinetic scheme （UGKS） is presented and used in this letter to study the slider air bearing problem. The UGKS solutions are first val- idated by comparison with direct simulation Monte Carlo results. ...The unified gas-kinetic scheme （UGKS） is presented and used in this letter to study the slider air bearing problem. The UGKS solutions are first val- idated by comparison with direct simulation Monte Carlo results. After valida- tion, the UGKS is used to study the air-bearing problem under different non- equilibrium conditions. On the surface of the slider, the dependency of the gas pressure and normal force on the Mach and Knudsen numbers are fully evaluated. The non-equilibrium effect on the force loading in the whole transition regime up to the free molecular limit is also studied.展开更多
文摘目前的气体动理学统一算法(unified gas kinetic scheme,简称UGKS)在求解高速流动问题时的计算效率,难以满足求解复杂工程问题的需求.为了提高该算法的计算效率,本文对模型方程的对流项和碰撞项进行了隐式处理,并针对UGKS界面通量与演化时间相关的特点,引入了演化时间平均界面通量,通过对控制方程矩阵进行近似LU分解(lower-upper decomposition),实现了隐式UGKS.不同来流马赫数的圆柱绕流算例测试表明,只要演化时间选取得当,隐式方法可以得到与显式方法完全相同的结果,且计算效率可以提高1～2个量级.
文摘辐射输运方程的数值模拟在天体物理、武器物理和惯性约束与磁约束聚变等研究中都起着非常重要的作用.在实际问题中,背景介质的不透明度系数决定了辐射光子在其中的传输行为.光性薄(不透明度系数小)的介质对辐射光子是透明的,光子与背景介质的相互作用弱,光子传输具有输运传播性质;而光性厚(不透明度系数大)的介质对辐射光子是不透明的,光子与背景介质的相互作用强,光子传输具有扩散性质.因此在辐射输运方程的计算中,如何设计既能得到光子输运传播性质又能捕捉光子扩散传播性质的渐近保持离散格式是目前一个非常活跃和前沿的研究方向.本文简要介绍近几年在辐射输运方程的渐近保持统一气体动理学格式(unified gas kinetic scheme,UGKS)研究方面的进展.本文主要以灰体辐射输运方程为例,详细介绍UGKS的构造方法并给出其渐近分析.同时,结合角度有限元方法和球谐函数展开的方法,介绍如何减弱/去除基于离散纵标法的UGKS具有射线效应的问题,以及相应的改进渐近保持格式.此外,也介绍了将渐近保持的UGKS应用拓展到考虑流体运动的完全辐射流体力学方程组.最后,用一些数值例子验证了格式的渐近保持性和保正性等性质.
基金supported by the Science and Technology Development foundation of China Academy of Engineering Physics(Grant Nos.2015B0202041,2015B0202040)the Science and Technology Development foundation of China Academy of Engineering Physics(Grant 2015B0202040)+2 种基金the Science and Technology Development foundation of China Academy of Engineering Physics(Grant No.2015B0202033)for LiNSFC(Grant No.11371068)for SunNSFC(Grant No.11371068)for Zeng
文摘In this paper,we consider the multi-dimensional asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations on distorted quadrilateral meshes.Different from the former scheme [J.Comput.Phys.285（2015）,265-279] on uniform meshes,in this paper,in order to obtain the boundary fluxes based on the framework of unified gas kinetic scheme（UGKS）,we use the real multi-dimensional reconstruction for the initial data and the macro-terms in the equation of the gray transfer equations.We can prove that the scheme is asymptotic preserving,and especially for the distorted quadrilateral meshes,a nine-point scheme [SIAM J.SCI.COMPUT.30（2008）,1341-1361] for the diffusion limit equations is obtained,which is naturally reduced to standard five-point scheme for the orthogonal meshes.The numerical examples on distorted meshes are included to validate the current approach.
文摘工程应用中的介质热辐射问题是典型的多尺度问题.基于Boltzmann输运方程建立的各类气体动理学格式,在多尺度瞬态问题中得到了广泛应用.为了克服显式求解方案中CFL条件等的限制,文章通过气体动理学格式实现稳态辐射输运方程的直接求解.SDUGKS格式由离散统一气体动理学格式(discrete unified gas kinetic scheme,DUGKS)的核心思想发展而来,应用于稳态问题计算.将SDUGKS格式进一步拓展到多尺度的稳态热辐射输运计算.SDUGKS格式继承了DUGKS格式沿特征线离散实现的界面重构,并通过隐式增量格式的单元更新实现对辐射强度的较正,采用逐次迭代法将辐射强度渐近收敛到稳定值.选用多组一维和二维不同尺度的辐射传热算例,通过与特定的解析解以及其他数值方法结果对比,检验了SDUGKS的计算精度和计算效率,并论证了它在多尺度问题中的渐进保持性质.
基金supported by Hong Kong Research Grant Council(621709,621011)HKUST grants SRFI11SC05 and RPC10SC11the Nanoscience and Nanotechnology Program at HKUST
文摘Sound wave propagation in rarefied monatomic gases is simulated using a newly developed unified gaskinetic scheme （UGKS）. The numerical calculations are carfled out for a wide range of wave oscillating frequencies. The corresponding rarefaction parameter is defined as the ratio of sound wave frequency to the intermolecular particle collision frequency. The simulation covers the flow regime from the continuum to free molecule one. The treatment of the os- cillating wall boundary condition and the methods for eval- uating the absorption coefficient and sound wave speed are presented in detail. The simulation results from the UGKS are compared to the Navier-Stokes solutions, the direct sim- ulation Monte Carlo （DSMC） simulation, and experimental measurements. Good agreement with the experimental data has been obtained in the whole flow regimes for the corresponding Knudsen number from 0.08 to 32. The cur- rent study clearly demonstrates the capability of the UGKS method in capturing the sound wave propagation and its usefulness for the rarefied flow study.
基金Project supported by the National Natural Science Foundation of China(Nos.11772281,91530319,and 11521091)the Hong Kong Research Grant Council(Nos.16207715 and 16206617)
文摘Compressible flows exhibit a diverse set of behaviors, where individual particle transports and their collective dynamics play different roles at different scales. At the same time, the atmosphere is composed of different components that require additional degrees of freedom for representation in computational fluid dynamics. It is challenging to construct an accurate and efficient numerical algorithm to faithfully represent multiscale flow physics across different regimes. In this paper, a unified gas-kinetic scheme(UGKS) is developed to study non-equilibrium multicomponent gaseous flows. Based on the Boltzmann kinetic equation, an analytical space-time evolving solution is used to construct the discretized equations of gas dynamics directly according to cell size and scales of time steps, i.e., the so-called direct modeling method. With the variation in the ratio of the numerical time step to the local particle collision time(or the cell size to the local particle mean free path), the UGKS automatically recovers all scale-dependent flows over the given domain and provides a continuous spectrum of the gas dynamics. The performance of the proposed unified scheme is fully validated through numerical experiments.The UGKS can be a valuable tool to study multiscale and multicomponent flow physics.
基金supported by Hong Kong Research Grant Council(621011 and 620813)HKUST(SRFI11SC05 and FSGRF13SC21)
文摘The unified gas-kinetic scheme （UGKS） is presented and used in this letter to study the slider air bearing problem. The UGKS solutions are first val- idated by comparison with direct simulation Monte Carlo results. After valida- tion, the UGKS is used to study the air-bearing problem under different non- equilibrium conditions. On the surface of the slider, the dependency of the gas pressure and normal force on the Mach and Knudsen numbers are fully evaluated. The non-equilibrium effect on the force loading in the whole transition regime up to the free molecular limit is also studied.
