一维含化学反应流体力学方程组的ALE间断有限元方法

研究一维含化学反应流体力学方程组的数值模拟方法.结合理想气体状态方程并利用HLLC解法器在各个单元边界处的数值通量，给出ALE间断有限元方法.高阶计算时，使用TVD斜率限制器对数值解可能产生的非物理振荡进行抑制.结果表明：该算法能够保持物理量的守恒性和高精度，并能够清晰地捕捉爆轰波的结构特征.     

摄动有限体积法重构近似高精度的意义

研讨有限体积(FV)方法重构近似高精度的作用问题．FV方法中积分近似采用中点规则为二阶精度时，重构近似高精度(精度高于二阶)的意义和作用是一个有争议的问题．利用数值摄动技术构造了标量输运方程的积分近似为二阶精度、重构近似为任意阶精度的迎风型和中心型摄动有限体积(PFV)格式．迎风PFV格式无条件满足对流有界准则(CBC)，中心型PFV格式为正型格式，两者均不会产生数值振荡解．利用PFV格式求解模型方程的数值结果表明：与一阶迎风和二阶中心格式相比，PFV格式精度高、对解的间断分辨率高、稳定性好、雷诺数的适用范围大，数值地“证实”重构近似高精度和PFV格式的实际意义和好处．     

球几何中子输运保正线性间断有限元格式

针对球几何中子输运方程线性间断有限元方法计算的负中子通量问题，构造了保正线性间断有限元格式，该格式保持中子角通量0阶矩和1阶矩。现有方法计算中子角通量非负时，采用传统的线性间断有限元方法，求解线性方程组；原方法计算出现负通量，则采用构造的保正格式，求解非线性方程组。编制了球几何中子输运问题保正格式程序模块，并集成到应用程序。数值算例表明构造的保正格式计算的中子通量非负，有效降低数值误差，提高数值计算的精度。     

局部间断Petrov-Galerkin方法在大气污染模型中的应用

构造数值模拟两类大气污染模型的局部间断Petrov-Galerkin方法。首先通过变量代换将大气污染模型方程转化为与之等价的一阶微分方程组,再利用间断Petrov-Galerkin方法求解微分方程组。该方法既可以选取不同的检验函数和试探函数空间,又可以保持间断Petrov-Galerkin方法的优势。同局部间断有限元方法相比,该方法的计算公式较为简便,数值算例表明该方法具有三阶精度。与有限体积元方法相比,该方法具有较小的误差。本算法可为大气污染模型的数值模拟提供实用工具。     

间断解问题的有限体积法

有限体积法（FVM）是集有限差分法和有限元方法之优点而发展起来的一种新的数值方法,首先利用守恒问题给出FVM的基本概念和基本过程,接下来就间断解问题的求解进一步介绍了FVM的一些比较活跃的构造方法,并给出了用其中几个方法模拟二维溃坝问题的数值结果,然后按问题分类介绍了FVM数值方法和实际应用的进展,最后给出了FVM理论分析的一些进展。     

粒子输运方程的线性间断有限元方法

将空间线性间断有限元方法应用于动态粒子输运方程的求解,数值算例表明,空间线性间断有限元方法在网格边界的数值精度方面明显高于指数格式和菱形格式,并且通量在时间上的微分曲线相对光滑,避免了指数格式、菱形格式数值解的非物理振荡现象.     

龙格库塔间断有限元方法在计算爆轰问题中的应用

构造求解带源项守恒律方程组的龙格库塔间断有限元(RKDG)方法,并分别结合源项的Strang分裂法和无分裂法数值求解模型守恒律方程和反应欧拉方程.为了和有限体积型WENO方法进行比较,设计计算源项的WENO重构格式.对一维带源项守恒律的计算表明,对于非刚性问题,RKDG方法比有限体积型WENO方法的误差更小;对于刚性问题,RKDG方法对于间断面位置的捕捉更为精确.对于一二维爆轰波问题的计算结果表明,RKDG方法对爆轰波结构的分辨和爆轰波位置的捕捉能力更强.     

定常不可压Navier—Stokes型方程的非线性Galerkin有限元算法的收敛性

给出数值求解二维定常不可压Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ型方程的非线性Ｇａｌｅｒｋｉｎ有限元算法，并分析了数值解的正则性和收敛性，当粗网格参数Ｈ和细网格参数ｈ满足关系式Ｈ＝Ｏ（ｈ＾１／２）时，该算法具有和Ｇａｌｅｒｋｉｎ有限元算法同阶的收敛精度，然而在计算上比Ｇａｌｅｒｋｉｎ有限元算法更为简单，可以节省可观的计算量。最后给出了数值试验，验证了上述结果。     

间断有限元方法求解一维非平衡辐射扩散方程

研究一维非平衡辐射扩散方程的数值方法.通过求解间断系数热传导方程的广义黎曼问题,得到一种带加权数值流量,基于该数值流量构造了一类新型的间断有限元方法.在时间离散上采用向后Euler方法,形成的非线性方程组采用Picard迭代求解.数值试验表明该方法具有捕捉大梯度的能力,而且能适应扩散系数间断的情形.     

一种高精度的修正Hermite-ENO格式

设计一种基于三单元具有六阶精度的修正Hermite-ENO格式（CHENO），求解一维双曲守恒律问题.CHENO格式利用有限体积法进行空间离散，在空间层上，使用ENO格式中的Newton差商法自适应选择模板.在重构半节点处的函数值及其一阶导数值时，利用Taylor展开给出修正Hermite插值使其提高到六阶精度，并设计了间断识别法与相应的处理方法以抑制间断处的虚假振荡；在时间层上采用三阶TVD Runge-Kutta法进行函数值及一阶导数值的推进.其主要优点是在达到高阶精度的同时具有紧致性.数值实验表明对一维双曲守恒律问题的求解达到了理论分析结果，是有效可行的.     

SOLUTION OF DYNAMIC RESPONSE OF FRAMED STRUCTURE USING PIECEWISE BIRKHOFF POLYNOMIAL

In this paper, the piecewise Birkhoff interpolation polynomials and the modal superposition method were employed for the solution of dynamic response of m.d.o.f. system. The related formulae are derived. Because an exact result can be obtained when each loading can be represented by a piecewise polynomial, the proposed method not only can considerably reduce computational effort compared to the traditional step-by-step integration solution technique, but also can thoroughly avoid the problems of accuracy, convergence and stability encountered in many other numerical procedures.     

Nonaxisymmetric radiative transfer in inhomogeneous cylindrical media with anisotropic scattering

In this paper, the control volume finite element method (CVFEM) is applied for the first time to solve nonaxisymmetric radiative transfer in inhomogeneous, emitting, absorbing and anisotropic scattering cylindrical media. Mathematical formulations as well as numerical implementation are given and the final discretized equations are based on similar meshes used for convective and conductive heat transfer in computational fluid dynamic analysis. In order to test the efficiency of the developed method, four nonaxisymmetric problems have been examined. Also, the grid dependence and the false scattering of the CVFEM are investigated and compared with the finite volume method and the discrete ordinates interpolation method.     

Time domain finite volume method for the transient response and natural characteristics of structural-acoustic coupling in an enclosed cavity

A time domain finite volume method（TDFVM）based on wave theory is developed to analyze the transient response and natural characteristics of structural-acoustic coupling problems in an enclosed cavity.In the present method,the elastic dynamic equations and acoustic equation in heterogeneous medium are solved in solid domains and fluid domains respectively.The structural-acoustic coupling is implemented according to the continuity condition of the particle velocity along the normal direction and the normal traction equilibrium condition on the interface.Several numerical examples are presented to validate the effectiveness and accuracy of the present TDFVM.Then the effects of water depth on the acoustic and vibration characteristics and the natural characteristics of a structural-acoustic coupling system are analyzed.The numerical results show that the increase of water depth leads to a stronger coupling between the water and structure and the decrease of natural frequencies of coupling system,The computational cost and memory of this method are small and it can be applicable to structural-acoustic coupling problems in the heterogeneous fluid.     

Padé–Legendre approximants for uncertainty analysis with discontinuous response surfaces

A novel uncertainty propagation method for problems characterized by highly non-linear or discontinuous system responses is presented. The approach is based on a Padé–Legendre (PL) formalism which does not require modifications to existing computational tools (non-intrusive approach) and it is a global method. The paper presents a novel PL method for problems in multiple dimensions, which is non-trivial in the Padé literature. In addition, a filtering procedure is developed in order to minimize the errors introduced in the approximation close to the discontinuities. The numerical examples include fluid dynamic problems characterized by shock waves: a simple dual throat nozzle problem with uncertain initial state, and the turbulent transonic flow over a transonic airfoil where the flight conditions are assumed to be uncertain. Results are presented in terms of statistics of both shock position and strength and are compared to Monte Carlo simulations.     

Boundary element-free method for elastodynamics

The moving least-square approximation is discussed first. Sometimes the method can form an ill-conditioned equation system, and thus the solution cannot be obtained correctly. A Hilbert space is presented on which an orthogonal function system mixed a weight function is defined. Next the improved moving least-square approximation is discussed in detail. The improved method has higher computational efficiency and precision than the old method, and cannot form an ill-conditioned equation system. A boundary element-free method (BEFM) for elastodynamics problems is presented by combining the boundary integral equation method for elastodynamics and the improved moving least-square approximation. The boundary element-free method is a meshless method of boundary integral equation and is a direct numerical method compared with others, in which the basic unknowns are the real solutions of the nodal variables and the boundary conditions can be applied easily. The boundary element-free method has a higher computational efficiency and precision. In addition, the numerical procedure of the boundary element-free method for elastodynamics problems is presented in this paper. Finally, some numerical examples are given.     

Boundary element-free method for elastodynamics

1 Introduction In recent years, more and more attention has been paid to researches on the meshless (or meshfree) method, which makes it a hot direction of computational mechanics[1,2]. The meshless method is the approximation based on nodes, then the large deformation and crack growth problems can be simulated with the method without the re-meshing technique. And the meshless method has some advantages over the traditional computa- tional methods, such as finite element method (FEM) and boun…     

Compact Accurately Boundary-Adjusting high-REsolution Technique for fluid dynamics

A novel high-resolution numerical method is presented for one-dimensional hyperbolic problems based on the extension of the original Upwind Leapfrog scheme to quasi-linear conservation laws. The method is second-order accurate on non-uniform grids in space and time, has a very small dispersion error and computational stencil defined within one space–time cell. For shock-capturing, the scheme is equipped with a conservative non-linear correction procedure which is directly based on the maximum principle. Plentiful numerical examples are provided for linear advection, quasi-linear scalar hyperbolic conservation laws and gas dynamics and comparisons with other computational methods in the literature are discussed.     

液氦温区两级4K脉冲管制冷机数值模拟

为模拟和动态显示工作在液氦温区的两级4K脉冲管制冷机内部工作过程和参数变化规律,发展了一种新的欧拉法-拉格朗日法数值计算模型。采用拉格朗日方法,直接跟踪脉冲管中气体微元随周期性压力波动的具体运行轨迹;采用欧拉法,直接模拟蓄冷器内部的动态参数变化。本文简单介绍该模型,并模拟了一典型两级4 K脉冲管制冷机各参数的变化情况,分析了多层磁性蓄冷材料对制冷机性能的影响情况。     

非规则形状介质内辐射-导热耦合传热的间断有限元求解

采用间断有限元法(discontinuous finite element method,DFEM)求解非规则形状介质内的辐射导热耦合传热问题,得到了典型非规则形状介质内辐射导热耦合传热问题的高精度数值结果.和传统连续型有限元方法不同,DFEM将计算区域划分成相互独立的离散单元,形函数的构造、未知量的加权近似以及控制方程的求解均在每一个离散单元上进行.通过在单元之间施加迎风格式的数值通量,DFEM保证了整个计算区域的连续性,因此这种方法兼具良好的几何灵活性和局部守恒性.推导了辐射传输方程和能量扩散方程的射导热耦合传热问题,得到了典型非规则形状介质内辐射导热耦合传热的高精度数值结果.     

Iterative method for dynamic condensation combined with substructuring scheme

An iterated improved reduced system (IIRS) procedure combined with substructuring scheme for both undamped and nonclassically damped structures is presented. Iterated IIRS method is an efficient reduction technique because the highly accurate eigenproperties from the repeatedly updated condensed matrices can be obtained without consuming expensive computational cost. However, single domain direct approach of this method to large structures requires much computational resources and even makes analysis intractable in the case only limited computer storage is available. These problems can be overcome by combining the substructuring scheme with IIRS procedure. The newly developed IIRS method combined with a substructuring scheme can provide an efficient methodology for large-scale eigenvalue problems. The validation of the present method and the evaluation of computational efficiency are demonstrated through the numerical examples.