三维PTT黏弹性液滴撞击固壁面问题的改进SPH模拟

基于光滑粒子动力学(smoothed particle hydrodynamics, SPH)方法,对三维Phan-Thien Tanner(PTT)黏弹性液滴撞击固壁面问题进行了数值模拟.为了有效地防止粒子穿透固壁,且缩减三维数值模拟所消耗的计算时间,提出了一种适合三维数值模拟的改进固壁边界处理方法.为了消除张力不稳定性问题,采用一种简化的人工应力技术.应用改进SPH方法对三维PTT黏弹性液滴撞击固壁面问题进行了数值模拟,精细地捕捉了液滴在不同时刻的自由面,讨论了PTT黏弹性液滴不同于Newton(牛顿)液滴的流动特征,分析了PTT拉伸参数对液滴宽度、高度和弹性收缩比等的影响.模拟结果表明,改进SPH方法能够有效而准确地描述三维PTT黏弹性液滴撞击固壁面问题的复杂流变特性和自由面变化特征.     

尖劈吸波体的研究和微波暗室的模拟

对尖劈形状吸波体的吸波性能进行了研究,并对导弹导引仿真实验用的微波暗室的性能进行分析和仿真.首先对尖劈体的二维反射性能进行了研究,从单条波线反射的原理出发,得到波束平面反射的统计模型.单条波线的反射通过数值模拟得到;波束反射模型则通过对数值模拟的结果进行统计和拟合得到,最终用多项式表示.对于一些简单或特殊的情况,也给出了解析解.通过分析发现,三维反射和二维反射之间有明确的关系.这种关系可以由三维入射角和反射次数决定,而反射次数可以通过二维模型得到.据此将平面反射模型扩展为三维反射模型,从而得到尖劈形状吸波体的三维反射模型.无回波暗室用于模拟没有背景微波辐射的环境,其关键在于选择合适的吸波材料.基于微波反射通量平衡原理,建立了考虑暗室墙面各点之间的相互影响的耦合模型,从而可以求解出在指定的发射源照射之下墙面各点的辐射强度分布.对模型的求解精度和收敛性进行了验证.基于此模型,对一个导弹引导试验进行了数值模拟,推算出了使用两种不同吸波材料时静区接收到的微波信号的信噪比.     

应用纵向面积分割多次曝光彩虹全息术合成多幅断层图像的三维单色像

采用纵向面积分割多次曝光彩虹全息术对多幅模拟层析图进行全息记录及白光再现合成三维单色像的实验和理论分析结果. 实验结果及其理论分析表明, 采用纵向面积分割法将主全息图分割为单元主全息图的记录方法可将多幅断层图像合成为无视觉灰度及几何畸变的三维单色像. 首次给出了采用该方法合成的一组医学断层图像的三维单色像.     

基于生态模型的企业竞争稳定性研究

从生态学视角来研究企业竞争的动态演化.针对现有研究文献中只讨论企业自身线性制约的三维系统的不足,构造了企业自身非线性制约的三维竞争模型,运用微分方程稳定性理论分析其稳定性,并由此揭示出企业间竞争的动力学机制.数值仿真结果表明,该三维模型能有效地模拟企业间竞争的动态演化规律.     

基于流-固耦合的水平缝压裂滤失规律研究

使用ABAQUS软件,建立了包含隔层和储层的层状多孔介质三维有限元模型,采用三维流固耦合实体单元和cohesive粘结单元描述岩石和裂缝的损伤行为.模拟直井水平缝的起裂和扩展以及液体滤失过程,分析得到了地质因素和工程因素对水平裂缝扩展参数的影响以及对滤失距离的影响,并进行了参数显著性分析,参数包括弹性模量、施工排量、渗透率、压力液粘度、地应力、孔隙压力和施工液量,最后利用模拟的规律认识指导实际井的压裂施工参数设计,压后效果显著.     

三维变系数椭圆型方程数值求解的交替方向法

利用交替方向隐格式研究了一类三维变系数椭圆方程的边值问题,给出了交替方向法的推导过程,建立了相应的误差分析,并进行了数值模拟,结果表明,该格式具有易于计算、求解精确度高等优点.     

基于元胞自动机理论的孔隙介质渗流的模拟应用

基于元胞自动机理论,建立了四方形网格元胞空间和摩尔型元胞邻居为类型的时空动态演化模型,并设计了模拟流程.利用两点分布表征岩石骨架和孔隙空间的非均匀分布,设定了可表征流体在多孔介质中渗流的演化规则,对二维和三维岩石多孔介质流体渗流的元胞矩阵进行了模拟,分析得到了两点分布不同概率对渗流数学参数的影响,分析了不同注入方式对流体最终占据孔隙空间的分布形态的影响,注入方式包括点注入、线注入和面注入,最终对真实CT扫描的数字二维岩心进行了流体渗流过程模拟.基于以上的研究,为流体在多孔介质中渗流模拟提供了有效的方法.     

时滞反馈控制下三维节能减排系统的稳定性分析

考虑三维节能减排系统,施加时滞反馈控制,建立反馈控制下的节能减排系统,求解平衡点,分析系统稳定性,得到出现Hopf分支的条件.利用中心流形定理和Poincare规范型方法讨论Hopf分支的性质,给定参数进行数值模拟,验证理论结果的有效性.     

Gause型食物链模型Hopf分支的存在性分析

考虑了一类三维时滞Gause型食物链模型.首先分析了共存平衡点稳定的条件,然后利用多项式理论分析了特征方程特征根的分布,得到了Hopf分支存在的条件,最后给出了几组数值模拟验证文中得到的结论,进一步预测了Hopf分支的全局存在性.     

核事故应急决策中水体放射性核素扩散仿真分析

放射性核素浓度的确定是放射性事故应急救援和辐射防护等工作的基础和前提,是放射性事故应急救援的重要组成部分.在一维扩散模型的基础上,建立了放射性核素瞬时污染点源三维扩散模型、连续污染点源三维扩散模型、连续污染点源三维扩散稳态模型,并编写程序,实现了对放射性核素污染扩散的快速仿真分析,较好模拟了水体中放射性核素的扩散趋势.研究结果在核事故应急救援过程中,可为相关部门制定救援方案及应急决策提供科学依据.     

Optimal rectangular MITC finite elements for Reissner–Mindlin plates

In recent years a family of finite elements named mixed interpolated tensorial components (MITC) has been introduced for the numerical approximation of Reissner–Mindlin plates. The elements have been proved to be locking free. In this article, we consider the MITC rectangular finite elements and show that it is possible to reduce the number of internal degrees of freedom in the approximation of the rotation field without losing order of convergence. Our mathematical analysis is carried out combining some results for the Stokes problem with the special features of the MITC finite elements. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 575–585, 1997     

平面弹性方程的非协调有限元分析

针对纯位移平面弹性问题,构造了两个无闭锁非协调有限元,单元对于Lamé常数λ一致收敛,证明了能量模和L2模误差分别为O(h2)和O(h3)．最后给出了数值试验验证了理论分析的正确性．     

A BE based finite element for crack propagation processes

A boundary element based finite macro element for the simulation of 3D crack propagation in the framework of linear elastic fracture mechanics is presented. While the major part of the numerical model is discretized with finite elements, a small domain containing the crack is meshed with boundary elements. By means of the Symmetric Galerkin BEM a stiffness formulation for the cracked BE domain is obtained which enables a direct FEM/BEM coupling. All necessary operations for the crack propagation are carried out within this boundary element based finite macro element and exploit the potential of the boundary integral formulation. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

New anisotropic a priori error estimates

Summary. We prove a priori anisotropic estimates for the and interpolation error on linear finite elements. The full information about the mapping from a reference element is employed to separate the contribution to the elemental error coming from different directions. This new error estimate does not require the “maximal angle condition”. The analysis has been carried out for the 2D case, but may be extended to three dimensions. Numerical experiments have been carried out to test our theoretical results. Received March 3, 2000 / Revised version received June 27, 2000 / Published online April 5, 2001     

Reconstructing initial data using observers: error analysis of the semi-discrete and fully discrete approximations

A new iterative algorithm for solving initial data inverse problems from partial observations has been recently proposed in Ramdani et al. (Automatica 46(10), 1616–1625, 2010). Based on the concept of observers (also called Luenberger observers), this algorithm covers a large class of abstract evolution PDE’s. In this paper, we are concerned with the convergence analysis of this algorithm. More precisely, we provide a complete numerical analysis for semi-discrete (in space) and fully discrete approximations derived using finite elements in space and an implicit Euler method in time. The analysis is carried out for abstract Schr?dinger and wave conservative systems with bounded observation (locally distributed).     

Numerical solution of the three‐dimensional parabolic equation with an integral condition

Developement of numerical methods for obtaining approximate solutions to the three dimensional diffusion equation with an integral condition will be carried out. The numerical techniques discussed are based on the fully explicit (1,7) finite difference technique and the fully implicit (7,1) finite difference method and the (7,7) Crank‐Nicolson type finite difference formula. The new developed methods are tested on a problem. Truncation error analysis and numerical examples are used to illustrate the accuracy of the new algorithms. The results of numerical testing show that the numerical methods based on the finite difference techniques discussed in the present article produce good results. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 193–202, 2002; DOI 10.1002/num.1040     

Stabilized spectral element approximation for the Navier–Stokes equations

The conforming spectral element methods are applied to solve the linearized Navier–Stokes equations by the help of stabilization techniques like those applied for finite elements. The stability and convergence analysis is carried out and essential numerical results are presented demonstrating the high accuracy of the method as well as its robustness. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 115–141, 1998     

Hybrid finite element methods for the Signorini problem

We study three mixed linear finite element methods for the numerical simulation of the two-dimensional Signorini problem. Applying Falk's Lemma and saddle point theory to the resulting discrete mixed variational inequality allows us to state the convergence rate of each of them. Two of these finite elements provide optimal results under reasonable regularity assumptions on the Signorini solution, and the numerical investigation shows that the third method also provides optimal accuracy.

     

Error analysis of discontinuous Galerkin finite element method for optimal control problem governed by the transport equation

This article discusses a priori and a posteriori error estimates of discontinuous Galerkin finite element method for optimal control problem governed by the transport equation. We use variational discretization concept to discretize the control variable and discontinuous piecewise linear finite elements to approximate the state and costate variable. Based on the error estimates of discontinuous Galerkin finite element method for the transport equation, we get a priori and a posteriori error estimates for the transport equation optimal control problem. Finally, two numerical experiments are carried out to confirm the theoretical analysis.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1493–1512, 2017     

A uniformly convergent B-spline collocation method on a nonuniform mesh for singularly perturbed one-dimensional time-dependent linear convection–diffusion problem

A numerical method is proposed for solving singularly perturbed one-dimensional parabolic convection–diffusion problems. The method comprises a standard implicit finite difference scheme to discretize in temporal direction on a uniform mesh by means of Rothe's method and B-spline collocation method in spatial direction on a piecewise uniform mesh of Shishkin type. The method is shown to be unconditionally stable and accurate of order O((Δx)²+Δt). An extensive amount of analysis has been carried out to prove the uniform convergence with respect to the singular perturbation parameter. Several numerical experiments have been carried out in support of the theoretical results. Comparisons of the numerical solutions are performed with an upwind finite difference scheme on a piecewise uniform mesh and exponentially fitted method on a uniform mesh to demonstrate the efficiency of the method.