无网格局部强弱法求解不规则域问题

无网格局部彼得洛夫-伽辽金(meshless local Petrov-Galerkin,MLPG)法是一种具有代表性的无网格方法,在计算力学领域得到广泛应用.然而,这种方法在边界上需执行积分运算,通常很难处理不规则求解域问题.为了克服MLPG法的这种局限性,提出了无网格局部强弱(meshless local strong-weak,MLSW)法.MLSW法采用MLPG法离散内部求解域,采用无网格介点(meshless intervention-point,MIP)法施加自然边界条件,并采用配点法施加本质边界条件,避免执行边界积分运算,可适用于求解各类复杂的不规则域问题.从理论上讲,这种结合式方法,既保持了MLPG法稳定而精确计算的优势,同时兼备配点型方法在处理复杂结构问题时简洁而灵活的优势,实现了弱式法和强式法的优势互补.此外,MLSW法采用移动最小二乘核(moving least squares core,MLSc)近似法来构造形函数,是对传统移动最小二乘(moving least squares,MLS)近似法的一种改进.MLSc使用核基函数代替通常的基函数,有利于数值求解的精确性和稳定性,而且其导数近似计算变得更为简单.数值算例结果初步表明:这种新方法实施简单,求解稳定、精确,表现出适合工程运用的潜力.     

非线性振动系统周期运动及其稳定性的数值研究

§1引言确定型非线性振动系统的运动可分类如下: 1.非定常运动;2.定常运动:(1)周期运动,(2)各态历经运动,(3)浑沌运动。其中非定常运动是一暂态过程,会随着时间的增长逐步衰减乃至实际上消失。定常运动中的各态历经运动,指系统至少有两个互不通约(即其比值为无理数)的振动频率,因此运动虽然局      

热传导问题灵敏度分析的伴随法

在热传导灵敏度分析的直接法的研究基础上，进一步探讨了稳态和瞬态热传导问题灵敏度分析的伴随法.推导了伴随法的计算列式，对于瞬态热传导问题，研究了瞬态约束处理的关键点方法，并提出伴随方程的精细积分解法。算例表明，稳态问题灵敏度计算，伴随法与直接法的结果是一致的；瞬态问题灵敏度计算，两种方法的精度相当。     

结构可靠指标求解的一种新的迭代方法

由国际安全度联合委员会推荐的一次二阶矩可靠度方法是国际工程界公认的一种较好的可靠度分析方法（简称ＪＣ法）。这种方法具有迭代格式简单，收敛快的特点。但许多实际计算表明，这种方法只适用于功能函数非线性程度不高的情况，否则不收敛。本文提出一种新的迭代方法，这种方法迭代格式与ＪＣ法同样简单，一般情况下均能保证收敛．     

An extended unresolved CFD-DEM coupling method for simulation of fluid and non-spherical particles

This study develops an extended unresolved CFD-DEM coupling method for simulation of the fluid–solid flow with non-spherical particles. The limitation of fluid grid size is discussed, by simulating the settling of a cylinder in a Newtonian fluid based on the resolved and unresolved CFD-DEM coupling method. Then, the calculation of porosity and the fluid–particle relative velocity based on the particle shape enlargement method for simulation of non-spherical particles is proposed. The availability of the particle shape enlargement method for the simulation of non-spherical particles with different sphericity is discussed in this work, by comparing it with the results from the equivalent diameter enlargement method. The limitation of the equivalent diameter enlargement method for non-spherical particles is revealed from the simulation results. Several typical cases are employed to elaborate and verify the extended unresolved CFD-DEM method based on particle shape enlargement method, by presenting a good consistency with the experimental results. It proves that the extended unresolved CFD-DEM method is suitable for different CFD grid size ratios, and consolidates that it is a universal calculation method for CFD-DEM coupling simulation.     

Second‐order entropy diminishing scheme for the Euler equations

This paper presents a method of controlling the water levels in a conduit system by employing optimal control theory and the finite element method. A shallow‐water equation is employed for the analysis of flow behaviour. Optimal control theory is utilized to obtain a control value for the target state value. The Sakawa–Shindo method is employed as a minimization technique. For the computational storage requirements, the time domain decomposition method is applied. The Crank–Nicolson method is used for temporal discretization. In addition to a method for optimally controlling water level, a method is presented for determining transversality conditions, the terminal condition of the Lagrange multiplier. Copyright © 2006 John Wiley & Sons, Ltd.     

求解平面和轴对称热冲击问题的热权函数有限元法

根据三维通用权函数法的普遍表达式,给出了轴对称问题和平面问题热权函数法的基本公式,并利用刚度阵导数法,将热权函数法与有限元法直接耦合,给出了二维问题的热权函数有限元法计算格式,实例计算表明,该文给出的计算方法具有极高的计算效率和满意的工程应用精度。     

求解转子系统突加不平衡响应方法的研究

本文提出并研究了用传递函数分析法求转子系统的突加不平衡响应,讨论了求这一响应的最佳数值积分方法。用本文提出的理论解方法及三种常用的数值积分法:尤拉后差法,纽马克法,及呼伯特法,对一模型转子算例进行了分析计算,结果表明;尤拉后差法最适用于求转子系统的突加不平衡响应。本文提出的传递函数分析法,可以和试验模态分析方法相结合,以求解复杂转子系统的突加不平衡响应。     

无网格介点法求解Helmholtz方程

作为一种配点型无网格法,无网格介点MIP法具有数值实施简单、计算精度高、运算高效和适用范围广等优点。Helmholtz方程是科学与工程问题中广泛应用的一类特殊方程,因此对MIP法求解此类方程的适用性进行了验证。利用MIP法的d适应性,给出了MIP法求解该方程的两种计算格式。在数值算例中,分别对平面规则域和不规则域上的一般Helmholtz方程,以及轴对称Helmholtz方程进行了数值分析。结果表明,MIP法完全适用于求解Helmholtz方程。而且,MIP法的计算精度和收敛性都优于普通配点法。此外,MIP法的两种计算格式中,L2C0型通常具有更好的计算效果,故建议将该计算格式作为MIP法求解该类方程的标准形式。     

Prediction of the turbulent boundary layer development by the lag-entrainment method

The lag-entrainment method, which is a well-established integral method for predicting the development of turbulent boundary layers, is used in this study to predict two-dimensional turbulent separated flow. The method is used in an inverse mode, in which the displacement thickness is specified together with other integral parameters of the boundary layer. It is concluded that the prediction of two-dimensional separated flow by an integral method is feasible, but there is a need for accurate data for both equilibrium and general separated flows for making a comparison.     

A Newton type iterative inverse method for heat-conduction problems

An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.     

非定常流函数涡量方程的一种数值解法的研究

对非定常流函数涡量方程的数值求解方法进行了改进,其中流函数一阶导数即速度项采用四阶精度的Hermitian公式,对流项由一般二阶精度的中心差分提高到四阶精度离散差分,包含温度方程在内的离散方程组采用ADI迭代方法求得定常解．以无内热体及有一内热体的封闭方腔内自然对流为例,进行了不同瑞利数（Ra）条件下的数值研究．结果表明,该方法推导简单,求解精度高且计算稳定,适用于封闭腔内高瑞利数复杂混合对流的数值模拟．     

结构可靠度分析FORM迭代算法的混沌控制

利用混沌控制原理对FORM收敛失败进行控制. 理清了全局性和局部性两类混沌反馈 控制各种方法的内在联系,说明稳定转换法和自适应调节法属于全局混沌反馈控制 方法,自适应调节法可视为稳定转换法的特例. 参 数调节混合法不过是松弛牛顿法的另一种表达形式,它们都属于局部混沌反馈控制方法. 阐 明了混沌反馈控制表达式与工程力学收敛控制迭代算法的对应关系. 也揭示了这些迭代算法 收敛控制措施的功效和局限性. 提出了一个以稳定转换法为主联合松弛牛顿法的混 沌反馈控制方法,对可靠度分析FORM迭代算法实现了周期振荡、分岔和混沌控制.     

正交各向异性稳态热传导问题的ICVEFG方法

本文将改进的复变量无单元Galerkin方法（Improved Complex Variable Element-free Galerkin method,ICVEFG）应用于求解正交各向异性介质中的稳态热传导问题,提出了正交各向异性稳态热传导问题的ICVEFG方法。采用罚函数法引入本质边界条件,推导了正交各向异性介质中的稳态热传导问题的Galerkin积分弱形式。采用改进的复变量移动最小二乘近似（Improved Complex Variable Moving least-squares approximation,ICVMLS）建立二维温度场问题的逼近函数,推导了相应的计算公式。编制了计算程序,对三个正交各向异性介质中的热传导问题进行了分析,说明了本文方法的有效性。     

A new method for computing turbulence in compressible laminar-turbulent transition and boundary layers--PSE＋DNS

A new method for computing laminar-turbulent transition and turbulence in compressible boundary layers is proposed. It is especially useful for computation of laminar-turbulent transition and turbulence starting from small-amplitude disturbances. The laminar stage, up to the beginning of the breakdown in laminar-turbulent transition, is computed by parabolized stability equations (PSE). The direct numerical simulation (DNS) method is used to compute the transition process and turbulent flow, for which the inflow condition is provided by using the disturbances obtained by PSE method up to that stage. In the two test cases incfuding a subsonic and a supersonic boundary layer, the transition locations and the turbulent flow obtained with this method agree well with those obtained by using only DNS method for the whole process. The computational cost of the proposed method is much less than using only DNS method.     

不可压黏性流问题的无网格法分析

不可压缩黏性流问题一般采用Navier-Stokes方程来描述，基于加权残值法，推导了问题的无网格伽辽金法（EFGM）离散Navier-Stokes方程，在时间域上采用分步方法计算，速度和压力由相互独立的方程以解耦的形式求解，并采用同阶移动最小二乘近似，在每一时间步中，对压力解和速度解采用了Newton-Raphson迭代法进行修正，最后将所得到的方法应用到剪切驱动空腔流问题中，验证了方法的有效性，且解的精度高、稳定性好。     

谱元法和高阶时间分裂法求解方腔顶盖驱动流

详细推导了谱元方法的具体计算公式和时间分裂法的具体计算过程 ;对一般的时间分裂法进行了改进 ,即对非线性步分别用 3阶 Adams-Bashforth方法和 4阶显式 Runge-Kutta法 ,粘性步采用 3阶隐式 Adams-Moulton形式 ,提高了时间方向的离散精度 ,同时还改进了压力边界条件 ,采用 3阶的压力边界条件 ;利用改进的时间分裂方法分解不可压缩 Navier-Stokes方程 ,并结合谱元法计算了移动顶盖方腔驱动流 ,提高了方法可以计算的 Re数 ,缩短了达到收敛的时间 ,并将结果与基准解进行比较 ;分析了移动顶盖方腔驱动流中 Re数对流场分布的影响。     

A Newton type iterative method for heat-conduction inverse problems

An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.     

计算含大位移动边界非定常流动的无网格算法

研究了可以计算处理包含大位移动边界非定常流动的无网格算法.创建了无网格方法中由于大位移边界运动造成的不合要求点的判断标准,采用高效的填充方法实现点云重构,运用线性插值方法得到新点参数,实现了局部重构处理大位移动点.流场计算方面,在计算域自动布点基础上,采用曲线逼近计算导数及HLLC格式计算数值通量,发展了求解基于无网格的ALE方程组的算法.最后,对活塞问题进行了模拟,结果与解析解相吻合,验证了算法的准确性;另外,计算比较了气流流过静止圆柱以及圆柱在静止流场运动流场,结果表明方法是成功的.     

势问题的重构核粒子边界无单元法

将重构核粒子法和势问题的边界积分方程方法结合,提出了势问题的重构核粒子边界无单元法. 推导了势问题的重构核粒子边界无单元法的公式,研究其数值积分方案,建立了重构核粒子边界无单元法的离散化边界积分方程,并推导了重构核粒子边界无单元法的内点位势的积分公式. 重构核粒子法形成的形函数具有重构核函数的光滑性,且能再现多项式在插值点的精确值,所以该方法具有更高的精度. 最后给出了数值算例,验证了所提方法的有效性和正确性. }