Compact difference approximation with consistent boundary condition

For simulating multi-scale complex flow fields it should be noted that all the physical quantities we are interested in must be simulated well. With limitation of the computer resources it is preferred to use high order accurate difference schemes. Because of their high accuracy and small stencil of grid points computational fluid dynamics (CFD) workers pay more attention to compact schemes recently. For simulating the complex flow fields the treatment of boundary conditions at the far field boundary points and near far field boundary points is very important. According to authors' experience and published results some aspects of boundary condition treatment for far field boundary are presented, and the emphasis is on treatment of boundary conditions for the upwind compact schemes. The consistent treatment of boundary conditions at the near boundary points is also discussed. At the end of the paper are given some numerical examples. The computed results with presented method are satisfactory.     

边界单元法三维应力分析计算程序

以三维弹性力学问题的基本解 (Kelvin解 )为基础开发了边界单元法三维弹性应力分析计算程序 ,并对其进行了验证。结果表明 ,该程序可用来求解三维弹性应力问题 ,尤其适用于三维应力集中问题。程序中采用动态分配内存 ,自动选点积分 ,并利用分块解法求解方程组 ,有效地节省了计算机资源 ,扩大了求解问题的规模。     

Shakedown Analysis of 3-D Structures Using the Boundary Element Method Based on the Static Theorem

The static shakedown theorem was reformulated for the boundary element method (BEM) rather than the finite element method with Melan‘s theorem, then used to develop a numerical solution procedure for shakedown analysis. The self-equilibrium stress field was constructed by a linear combination of several basis self-equilibrium stress fields with undetermined parameters. These basis self-equilibrium stress fields were expressed as elastic responses of the body to imposed permanent strains obtained using a 3-D BEM elastic-plastic incremental analysis. The lower bound for the shakedown load was obtained from a series of nonlinear mathematical programming problems solved using the Complex method. Numerical examples verified the precision of the present method.     

测度链上非线性微分方程的特征值问题

讨论了非线性特征值问题 u△△(t) λa(t)f(u(δ(t) ) ) =0 ,t∈ [0 ,1]u(0 ) =0 =u(δ(1) ) 正解的存在性 .这里 [0 ,1]是一可测链 ,a与f取正值 ,且limx→ 0 f(x)x 与limx→∞f(x)x 不一定存在     

具有Neumann边界及临界Sobolev指数的半线性抛物方程

主要运用能量方法及稳定集和不稳定集的观点,研究一类半线性抛物方程的整体解和局部解的存在性及爆破问题.这里Ω是RN(N≥3)上的光滑有界区域,2 是临界Sobolev指数.     

拟线性双曲方程组带有一般形式的边值问题的整体经典解

证明了拟线性双曲方程组带有一般形式的边值问题的整体经典解的存在惟一性。     

拟线性双曲方程组带有一般形式的自由边值问题的整体经典解

证明了拟线性双曲方程组带有一般形式的自由边值问题的整体经典解的存在惟一性。     

应用于微波电磁场FDTD数值模拟的吸收边界条件

在应用FDTD方法模拟复杂空间中微波电磁场分布时,不仅要考虑入射电磁波的贡献,还要考虑反射电磁波的影响。对于后者,一般常采用吸收边界条件的方法予以处理。针对TEM电磁波的传播特点,推导出了一种新的吸收边界条件的表达式,并利用其对微波法金刚石膜沉积系统沉积室中电磁场的分布进行了数值模拟。     

三维接触边界元法的一种误差直接估计

将作者所在研究组提出的二维弹性力学问题边界元解误差的直接估计推广到三维问题,给出了确定与域内解连续的边界位移的一种精确有效的方法。在此基础上提出将接触体接触单元间与域内解连续的边界位移之差的某种度量作为三维弹性接触问题边界元法的一种误差直接估计,并且提出了三维弹性接触问题的一种自适应边界元法计算方案。这种方案为确定没有解析解可作比较的复杂接触问题的边界元解精度提供了可能。文中对于三维弹性接触问题,给出了一个计算误差直接估计及自适应边界元法的算例。     

平整轧制压力预测的边界润滑弹塑性模型

平整工艺是使冷轧薄钢板获得一定表面形貌和粗糙度的关键工序,准确预测平整轧制压力的大小是进一步研究其形貌复映关系的理论基础。针对平整工艺压下率s一般小于10%的特点及其润滑特征建立了轧制压力预测的边界润滑弹塑性模型。通过理论计算表明:在s<10%时,轧辊的弹性变形和钢板的弹性恢复对轧制力预测有显著影响,是不能忽略的。使用 130型两辊轧制实验机,压下率为4.6%～24.4%的实验证明:该模型的理论计算结果与实验结果相对误差小于20%,比边界润滑刚性模型(其相对误差小于30%)具有更高的准确性。