In this paper,we apply local discontinuous Galerkin(LDG)methods for pattern formation dynamical model in polymerizing actin focks.There are two main dificulties in designing effective numerical solvers.First of all,th...In this paper,we apply local discontinuous Galerkin(LDG)methods for pattern formation dynamical model in polymerizing actin focks.There are two main dificulties in designing effective numerical solvers.First of all,the density function is non-negative,and zero is an unstable equilibrium solution.Therefore,negative density values may yield blow-up solutions.To obtain positive numerical approximations,we apply the positivitypreserving(PP)techniques.Secondly,the model may contain stif source.The most commonly used time integration for the PP technique is the strong-stability-preserving Runge-Kutta method.However,for problems with stiff source,such time discretizations may require strictly limited time step sizes,leading to large computational cost.Moreover,the stiff source any trigger spurious filament polarization,leading to wrong numerical approximations on coarse meshes.In this paper,we combine the PP LDG methods with the semi-implicit Runge-Kutta methods.Numerical experiments demonstrate that the proposed method can yield accurate numerical approximations with relatively large time steps.展开更多
辐射输运方程的数值模拟在天体物理、武器物理和惯性约束与磁约束聚变等研究中都起着非常重要的作用.在实际问题中,背景介质的不透明度系数决定了辐射光子在其中的传输行为.光性薄(不透明度系数小)的介质对辐射光子是透明的,光子与背景...辐射输运方程的数值模拟在天体物理、武器物理和惯性约束与磁约束聚变等研究中都起着非常重要的作用.在实际问题中,背景介质的不透明度系数决定了辐射光子在其中的传输行为.光性薄(不透明度系数小)的介质对辐射光子是透明的,光子与背景介质的相互作用弱,光子传输具有输运传播性质;而光性厚(不透明度系数大)的介质对辐射光子是不透明的,光子与背景介质的相互作用强,光子传输具有扩散性质.因此在辐射输运方程的计算中,如何设计既能得到光子输运传播性质又能捕捉光子扩散传播性质的渐近保持离散格式是目前一个非常活跃和前沿的研究方向.本文简要介绍近几年在辐射输运方程的渐近保持统一气体动理学格式(unified gas kinetic scheme,UGKS)研究方面的进展.本文主要以灰体辐射输运方程为例,详细介绍UGKS的构造方法并给出其渐近分析.同时,结合角度有限元方法和球谐函数展开的方法,介绍如何减弱/去除基于离散纵标法的UGKS具有射线效应的问题,以及相应的改进渐近保持格式.此外,也介绍了将渐近保持的UGKS应用拓展到考虑流体运动的完全辐射流体力学方程组.最后,用一些数值例子验证了格式的渐近保持性和保正性等性质.展开更多
基金supported by the Natural Science Foundation of Shandong Province(ZR2021MA001)the Fundamental Research Funds for the Central Universities(20CX05011A)+1 种基金supported by National Natural Science Foundation of China Grant 11801569supported by NSF grant DMS-1818467 and Simons Foundation 961585.
文摘In this paper,we apply local discontinuous Galerkin(LDG)methods for pattern formation dynamical model in polymerizing actin focks.There are two main dificulties in designing effective numerical solvers.First of all,the density function is non-negative,and zero is an unstable equilibrium solution.Therefore,negative density values may yield blow-up solutions.To obtain positive numerical approximations,we apply the positivitypreserving(PP)techniques.Secondly,the model may contain stif source.The most commonly used time integration for the PP technique is the strong-stability-preserving Runge-Kutta method.However,for problems with stiff source,such time discretizations may require strictly limited time step sizes,leading to large computational cost.Moreover,the stiff source any trigger spurious filament polarization,leading to wrong numerical approximations on coarse meshes.In this paper,we combine the PP LDG methods with the semi-implicit Runge-Kutta methods.Numerical experiments demonstrate that the proposed method can yield accurate numerical approximations with relatively large time steps.
文摘辐射输运方程的数值模拟在天体物理、武器物理和惯性约束与磁约束聚变等研究中都起着非常重要的作用.在实际问题中,背景介质的不透明度系数决定了辐射光子在其中的传输行为.光性薄(不透明度系数小)的介质对辐射光子是透明的,光子与背景介质的相互作用弱,光子传输具有输运传播性质;而光性厚(不透明度系数大)的介质对辐射光子是不透明的,光子与背景介质的相互作用强,光子传输具有扩散性质.因此在辐射输运方程的计算中,如何设计既能得到光子输运传播性质又能捕捉光子扩散传播性质的渐近保持离散格式是目前一个非常活跃和前沿的研究方向.本文简要介绍近几年在辐射输运方程的渐近保持统一气体动理学格式(unified gas kinetic scheme,UGKS)研究方面的进展.本文主要以灰体辐射输运方程为例,详细介绍UGKS的构造方法并给出其渐近分析.同时,结合角度有限元方法和球谐函数展开的方法,介绍如何减弱/去除基于离散纵标法的UGKS具有射线效应的问题,以及相应的改进渐近保持格式.此外,也介绍了将渐近保持的UGKS应用拓展到考虑流体运动的完全辐射流体力学方程组.最后,用一些数值例子验证了格式的渐近保持性和保正性等性质.
