流体力学方程的间断有限元方法

在二维区域三角形网格上应用一阶、二阶和三阶精度间断有限元方法,对流体力学方程和方程组进行了数值模拟.计算结果与差分方法计算结果比较,认为间断有限元方法在求解复杂边界条件和区域问题上有一定的优势.     

求解二维Euler方程有限单元边插值的降维重构算法

数值求解二维Euler方程的有限体积法（如k-exact，WENO重构、紧致重构等），无一例外地要进行耗时的网格单元上的二维重构.然而这些二维重构最后仅用于确定网格单元边界上高斯积分点处的解值，单元上二维重构似乎并非必需的.因此，文章提出用网格边上的一维重构来取代有限体积法中网格单元上的二维重构，分别在一致矩形网格和非结构三角形网格上发展了基于网格边重构的求解二维Euler方程的新方法，称为降维重构算法.数值算例表明该算法可以计算有强激波的无黏流动问题，且有较高的计算效率.      

基于变分原理的二维热传导方程差分格式

研究二维热传导方程的差分数值模拟.用变分原理在不规则结构网格上建立热流通量形式的差分格式.将热流通量作为未知函数求泛函极值,并与温度函数联立求解.克服通常九点格式用插值方法计算网格边界上的热传导系数和网格结点上的温度所引入的误差.     

多流管网格上九点格式的节点插值方法

多流管方法是二维多介质辐射流体力学数值模拟中一类常用的求解方法,它采用Lagrange-Euler混合型四边形网格,称为多流管网格。通常其网格品质高于一般的四边形网格。在这类网格上,可以利用网格特性对九点扩散格式中的节点插值方法进行改进。本文利用调和平均点和梯度离散构造的方法提出几种节点插值方法。并给出数值实验,说明现有应用程序中的节点插值方法损失精度,而新的节点插值方法能够使得九点格式在多流管网格上具有二阶精度。     

一种新的非结构四边形网格生成方法及软件

提出一种间接生成非结构四边形网格的方法.采用Delaunay三角形网格生成的算法生成三角形单元,在三角形单元中定义生成四边形网格的前沿,根据三角形的特点(边、角大小)以及生成四边形前沿的方向,将两个三角形合并为一个四边形单元.此方法的特点是不增加节点,仅改变点与点之间的拓扑连接,使单元基本变为四边形占优的单元.然后研究四边形占优的网格转换为100%的四边形网格的转换模版.并应用于实际物理问题.数值实验显示了该方法生成非结构四边形网格的特点.     

Galerkin时空耦合谱元法求解声波动方程

提出了时空耦合谱元方法,并将其用于带第一类边界条件的非齐次一维、二维、三维波动方程的求解。分别采用四边形、六面体和超六面体作为计算单元,在每个单元内采用Chebyshev多项式的极值点作为Lagrange插值节点,并且探讨了区域剖分方式对计算精度的影响。时空耦合谱元法能够得到精度很高的数值结果,并且其色散随时间推移是稳定的;当总网格节点数相同时,不同的网格剖分方式所得数值误差不同,当空间方向Chebyshev多项式的阶数较高和时间方向Chebyshev多项式的阶数较低时,得到的数值精度较高;在总节点数相同的情况下,与时间全域方式相比,逐时间子区域方式计算所需要的时间更经济,两种方式可以得到相同的精度。结果表明:时空耦合谱元方法使时空方向精度相匹配,可以提高整体精度;空间方向的Chebyshev多项式对数值精度起主要影响作用;时间子区域方式的采用可以扩大问题的计算区域。      

网格单元形状对数值计算的影响

实现了二维非结构化网格上的N-S方程的离散,编制了非结构化同位网格上的SIMPLE程序,针对CFD/NHT 中几个有基准解的问题进行了试算考核。结果表明,在步长基本相同的情况下,四边形网格的收敛速度和解的精度都优于 三角形网格,尤其是收敛速度更明显;使用考虑边界贴体特性的边界层网格,解的精度优于均匀网格。     

基于特征理论的二维可压缩流动的二阶拉氏算法

提出一种求解二维拉氏可压缩流体力学方程的中心型二阶精度有限体积方法.利用特征理论构造网格节点处的局部近似演化算子,算子用来求解网格节点处的速度及压力,利用这些物理量更新节点位置及计算网格界面通量.通过结合一定的重构方案,该方法达到时、空二阶精度,并且形式简单、计算量小,适用于结构网格与非结构网格.典型数值实验表明,本文格式具有良好的收敛性、对称性及鲁棒性,且能自然地求解多物质流动问题.     

基于特征理论的二维无粘Lagrange流体力学有限体积法

提出一个求解二维无粘Lagrange流体力学方程的中心型有限体积方法.采用特征理论求解网格节点处的速度及压力,并利用这些物理量更新节点位置及计算网格界面通量.方法适用于结构网格与非结构网格.典型数值实验的结果表明,格式具有较好的收敛性、对称性、能量守恒性及鲁棒性,且能自然地求解多物质流动问题.     

多流管网格上九点格式的节点插值方法（英文）

多流管方法是二维多介质辐射流体力学数值模拟中一类常用的求解方法,它采用Lagrange-Euler混合型四边形网格,称为多流管网格。通常其网格品质高于一般的四边形网格。在这类网格上,可以利用网格特性对九点扩散格式中的节点插值方法进行改进。本文利用调和平均点和梯度离散构造的方法提出几种节点插值方法。并给出数值实验,说明现有应用程序中的节点插值方法损失精度,而新的节点插值方法能够使得九点格式在多流管网格上具有二阶精度。     

三维不规则区域热传导问题无网格方法的数值模拟

配点型无网格法是纯无网格法,它不需要任何背景网格,效率高。本文用加权最小二乘配点方法(Weighted Least-Squares Collocation Method-WLSCM)计算不规则区域热传导问题,形函数采用径向基函数近似。通过二维具有分析解的实例表明WLSCM方法精度高,稳定性好且具有较高的计算效率。此外,将WLSCM方法应用于工程中常见的三维不规则区域热传导问题,结果表明:WLSCM方法的计算结果与FLUENT的计算结果符合很好。     

Point-wise hierarchical reconstruction for discontinuous Galerkin and finite volume methods for solving conservation laws

We develop a new hierarchical reconstruction (HR) method  and  for limiting solutions of the discontinuous Galerkin and finite volume methods up to fourth order of accuracy without local characteristic decomposition for solving hyperbolic nonlinear conservation laws on triangular meshes. The new HR utilizes a set of point values when evaluating polynomials and remainders on neighboring cells, extending the technique introduced in Hu, Li and Tang [9]. The point-wise HR simplifies the implementation of the previous HR method which requires integration over neighboring cells and makes HR easier to extend to arbitrary meshes. We prove that the new point-wise HR method keeps the order of accuracy of the approximation polynomials. Numerical computations for scalar and system of nonlinear hyperbolic equations are performed on two-dimensional triangular meshes. We demonstrate that the new hierarchical reconstruction generates essentially non-oscillatory solutions for schemes up to fourth order on triangular meshes.     

变系数瞬态热传导问题边界元格式的特征正交分解降阶方法

提出了一种基于边界元法求解变系数瞬态热传导问题的特征正交分解(POD)降阶方法,重组并推导出变系数瞬态热传导问题适合降阶的边界元离散积分方程,建立了变系数瞬态热传导问题边界元格式的POD降阶模型,并用常数边界条件下建立的瞬态热传导问题的POD降阶模态,对光滑时变边界条件瞬态热传导问题进行降阶分析.首先,对一个变系数瞬态热传导问题,建立其边界域积分方程,并将域积分转换成边界积分;其次,离散并重组积分方程,获得可用于降阶分析的矩阵形式的时间微分方程组;最后,用POD模态矩阵对该时间微分方程组进行降阶处理,建立降阶模型并对其求解.数值算例验证了本文方法的正确性和有效性.研究表明:1)常数边界条件下建立的低阶POD模态矩阵,能够用来准确预测复杂光滑时变边界条件下的温度场结果;2)低阶模型的建立,解决了边界元法中采用时间差分推进技术求解大型时间微分方程组时求解速度慢、算法稳定性差的问题.     

A new complex variable meshless method for transient heat conduction problems

<正>In this paper,based on the improved complex variable moving least-square(ICVMLS) approximation,a new complex variable meshless method(CVMM) for two-dimensional(2D) transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations,and the essential boundary conditions are imposed by the penalty method.As the transient heat conduction problems are related to time,the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization.Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained.In order to demonstrate the applicability of the proposed method,numerical examples are given to show the high convergence rate,good accuracy,and high efficiency of the CVMM presented in this paper.     

Three Boundary Meshless Methods for Heat Conduction Analysis in Nonlinear FGMs with Kirchhoff and Laplace Transformation

This paper presents three boundary meshless methods for solving problems of steady-state and transient heat conduction in nonlinear functionally graded materials (FGMs). The three methods are, respectively, the method of fundamental solution (MFS), the boundary knot method (BKM), and the collocation Trefftz method (CTM) in conjunction with Kirchhoff transformation and various variable transformations. In the analysis, Laplace transform technique is employed to handle the time variable in transient heat conduction problem and the Stehfest numerical Laplace inversion is applied to retrieve the corresponding time-dependent solutions. The proposed MFS, BKM and CTM are mathematically simple, easy-to-programming, meshless, highly accurate and integration-free. Three numerical examples of steady state and transient heat conduction in nonlinear FGMs are considered, and the results are compared with those from meshless local boundary integral equation method (LBIEM) and analytical solutions to demonstrate the efficiency of the present schemes.     

Meshless analysis of three-dimensional steady-state heat conduction problems

Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson’s equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method (RKPM). A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.     

Transient heat conduction analysis using the MLPG method and modified precise time step integration method

The meshless local Petrov–Galerkin (MLPG) method in conjunction with the modified precise time step integration method in the time domain is proposed for transient heat conduction analysis in this paper. The MLPG method is often referred to as a truly meshless method because it requires no elements or background cells for either field interpolation or background integration. Local weak forms are developed using weighted residual method locally from the partial differential equation of transient heat conduction. In order to simplify the treatment of essential boundary conditions, the natural neighbour interpolation (NNI) is employed for the construction of trial functions. Moreover, the three-node triangular FEM shape functions are taken as test functions to reduce the order of integrands involved in domain integrals. The semi-discrete heat conduction equation is solved numerically with modified precise time step integration method in the time domain. The availability and accuracy of the present method for transient heat conduction analysis are tested through numerical examples.     

Central discontinuous Galerkin methods for ideal MHD equations with the exactly divergence-free magnetic field

In this paper, central discontinuous Galerkin methods are developed for solving ideal magnetohydrodynamic (MHD) equations. The methods are based on the original central discontinuous Galerkin methods designed for hyperbolic conservation laws on overlapping meshes, and use different discretization for magnetic induction equations. The resulting schemes carry many features of standard central discontinuous Galerkin methods such as high order accuracy and being free of exact or approximate Riemann solvers. And more importantly, the numerical magnetic field is exactly divergence-free. Such property, desired in reliable simulations of MHD equations, is achieved by first approximating the normal component of the magnetic field through discretizing induction equations on the mesh skeleton, namely, the element interfaces. And then it is followed by an element-by-element divergence-free reconstruction with the matching accuracy. Numerical examples are presented to demonstrate the high order accuracy and the robustness of the schemes.