Lax-Wendroff时间离散的自适应间断有限元方法求解三维可压缩欧拉方程

应用自适应LWDG方法求解三维双曲守恒律方程组,与传统的二阶RKDG方法相比,该方法具有计算量小和精度高的特点.给出一种自适应策略,其中均衡折中策略适用于非相容四面体网格.将二维情形下的后验误差指示子推广到三维双曲守恒律方程组中,数值实验证明了方法的有效性.     

含运动物体流场计算的自适应非结构二维网格生成方法

采用自动插点的Delaunay方法和局部网格重新生成方法,对含运动边界自适应非结构二维网格生成方法进行了系统研究,实现了带有非定常运动物体流场的数值计算,并进行了数值实验.     

强爆炸火球辐射流体自适应网格高精度数值模拟

提出一种自适应网格方法,应用于基于Euler方法的强爆炸辐射流体高精度数值求解.通过与Zinn数值结果对比,验证该方法的正确性.研究自适应网格对冲击波和光辐射输出模拟精度的影响,对比不同网格尺度下的计算耗时.在相同的条件下,使用自适应网格与均匀网格加密3倍得到的冲击波超压分布、光辐射输出演化接近,计算效率提升约8.5倍...     

自适应间断有限元方法求解双曲守恒律方程

研究自适应Runge-Kutta间断Galerkin (RKDG)方法求解双曲守恒律方程组,并提出两种生成相容三角形网格的自适应算法.第一种算法适用于规则网格,实现简单、计算速度快.第二种算法基于非结构网格,设计一类基于间断界面的自适应网格加密策略,方法灵活高效.两种方法都具有令人满意的计算效果,而且降低了RKDG的计算量.     

求解双曲型守恒律方程的一类自适应多分辨方法

针对双曲型守恒律方程问题,发展一种有效的自适应多分辨分析方法.通过对嵌套网格上的数值解构造离散多分辨分析,建立小波系数与多层嵌套网格点之间的对应关系.对于小波系数较大的网格点采用高精度WENO格式计算,其余区域则直接采用多项式插值.数值试验表明,该方法在保持原规则网格方法的精度和分辨率的同时,显著地减少计算的CPU时间.     

基于任意拉格朗日—欧拉型移动网格的反应流广义黎曼问题方法

提出一种在自由重映移动网格下的广义黎曼问题方法模拟反应流.该方法基于显式的自由重映移动网格广义黎曼问题的解.为保证在时间和空间上的高精度，应用广义黎曼问题方法构造数值通量.为保证反应区的高分辨率，采用变分法生成自适应移动网格.该方法不仅能够保证网格质量，而且能有效地避免任意拉格朗日—欧拉方法中由于显式重映过程而带来的数值误差.包括CJ爆轰及不稳定爆轰的数值实验说明该格式的精确性和鲁棒性，证明这种移动网格下的二阶广义黎曼问题方法可以较好地捕捉反应流的间断与光滑结构.     

基于WENO重构的高阶移动网格动理学格式

提出一种基于WENO重构的高阶(至少三阶)移动网格动理学格式.利用流体力学方程的积分形式得到移动网格上离散格式,再利用自适应移动网格方法移动网格,进而得到网格速度,利用WENO重构得到高阶插值多项式,最后使用时间方向上精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高精度、高分辨率的特点.     

反应流体的移动网格动理学格式

在移动网格上构造一种反应流的动理学格式.首先利用BGK模型推导含化学反应的流体力学方程组,并利用其积分形式构造移动网格上离散格式,再利用自适应移动网格方法得到网格速度,最后利用时间精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.一维与二维的数值实验表明这种格式同时具有高精度、高分辨率的特点.     

振动NACA0012翼型的移动网格数值模拟

利用流体力学方程的积分形式给出非结构移动网格上离散格式,利用自适应移动网格方法移动网格,进而得到网格速度.对振动Naca0012翼型问题,分三种类型确定网格速度,再结合Riemann问题的解法器构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高效、高分辨率的特点.     

一种基于预估-修正策略的非定常流网格自适应加密/稀疏方法

给出一种非定常流动数值模拟的网格自适应处理方法.在"求解流动方程-自适应调整网格"的流程中,引入预估-修正步.根据自适应周期内每个时间步上的流场预估解,计算单元上的事后误差估算值.建立考虑解演变的网格自适应指示器,并进行多层次单元加密-稀疏的动态网格自适应处理.在自适应网格上重新计算流场.每个自适应周期中,流动演变区域的网格获得加密;而前一个周期中的特征现象已离开区域的网格被稀疏.应用边界非协调的当地DFD(Domain-Free Discretization)方法求解流动方程.为验证网格自适应处理方法,针对静止圆柱和自推进游鱼的流动进行了数值实验.     

Novel energy dissipative method on the adaptive spatial discretization for the Allen–Cahn equation

We propose a novel energy dissipative method for the Allen–Cahn equation on nonuniform grids. For spatial discretization, the classical central difference method is utilized, while the average vector field method is applied for time discretization. Compared with the average vector field method on the uniform mesh, the proposed method can involve fewer grid points and achieve better numerical performance over long time simulation. This is due to the moving mesh method, which can concentrate the grid points more densely where the solution changes drastically. Numerical experiments are provided to illustrate the advantages of the proposed concrete adaptive energy dissipative scheme under large time and space steps over a long time.     

ENO守恒插值(重映)方法及其在流体计算中的应用

在ENO(Essentially Non-oscillatory)守恒插值方法的基础上,分析和研究现今流体力学计算中涉及的几类网格技术:重叠网格技术、自适应加密技术和运动网格技术.基于ENO插值多项式构造的重映方法具有良好的守恒性,可以有效保证数据传递中物理量的总体守恒.提出该类守恒插值方法在以上几种网格技术中的一些应用前景,并给出一些数值算例.     

Adaptive mesh refinement for stochastic reaction–diffusion processes

We present an algorithm for adaptive mesh refinement applied to mesoscopic stochastic simulations of spatially evolving reaction–diffusion processes. The transition rates for the diffusion process are derived on adaptive, locally refined structured meshes. Convergence of the diffusion process is presented and the fluctuations of the stochastic process are verified. Furthermore, a refinement criterion is proposed for the evolution of the adaptive mesh. The method is validated in simulations of reaction–diffusion processes as described by the Fisher–Kolmogorov and Gray–Scott equations.     

块结构自适应网格上任意区域和任意界面的数值积分

研究块结构自适应网格上计算任意界面上和任意区域内的数值积分方法,其中任意界面和任意区域通过一个水平集函数表示.首先介绍在一致网格上任意界面上和任意区域内的数值积分方法.然后,将该方法推广到块结构的自适应网格上.数值算例表明,自适应网格方法有二阶精度.同一致网格方法相比,自适应网格方法显著地减少了计算机存储量的需求.     

L-stability of vertex-based MUSCL finite volume schemes on unstructured grids: Simulation of incompressible flows with high density ratios

This work is devoted to the design of multi-dimensional finite volume schemes for solving transport equations on unstructured grids. In the framework of MUSCL vertex-based methods we construct numerical fluxes such that the local maximum property is guaranteed under an explicit Courant–Friedrichs–Levy condition. The method can be naturally completed by adaptive local mesh refinements and it turns out that the mesh generation is less constrained than when using the competitive cell-centered methods. We illustrate the effectiveness of the scheme by simulating variable density incompressible viscous flows. Numerical simulations underline the theoretical predictions and succeed in the computation of high density ratio phenomena such as a water bubble falling in air.     

Adaptive variational curve smoothing based on level set method

This paper presents an adaptive method for variational curve smoothing based on level set implementation. A suitable cost functional is minimized via solving the derived Euler–Lagrangian equation, of which the discretization is conducted on unstructured triangular meshes by employing a simple and effective finite volume scheme. Through adaptive refinement of the mesh, the geometry features of the given curve can be well resolved in a cost-effective way. Various numerical experiments demonstrate the effectiveness and efficiency of the proposed approach.     

非结构网格自适应的新策略

对于用推进波阵面法生成二维非结构网格,提出了一种网格自适应的新方法。在网格需要加密处,布置一些小的三角形单元,并以这些单元的边作为初始推进波阵面的一部分。在推进过程中,小单元附近自然生长出尺寸逐渐变大的单元。对有激波的流动,在网格图形上,可形成十分清楚的激波图象。