局部人工边界稳定性的一种分析方法

采用时域局部人工边界同有限元或有限差分相结合的方法分析无界域中的波动问题时可能出现由边界引起的数值失稳现象．本文提出一种分析局部人工边界稳定性的实用方法．该方法从时域逐步积分格式出发，考虑了边界节点运动和内部节点运动的耦合以及离散化的影响，可以给出局部人工边界在多维波动数值模拟中的稳定性近似准则．文中对多次透射边界在一维及二维出平面波动模拟中稳定性的分析结果，表明了该方法的合理性     

双相各向异性介质弹性波场有限差分正演模拟

从双相各向异性介质模型出发,以Boit理论为基础,推导了斜方晶系各向异性介质-阶弹性波动方程,引入固、流体密度比和孔隙几何参数,将Biot方程系数简化为测量简单、物理意义明确的物理量,采用交错网格技术建立了各向异性孔隙介质波动方程的高精度差分格式,并首次对这类差分格式的频散特性和稳定性作了详细分析讨论,解决了计算稳定性和边界反射问题,与解析解的对比以及理论模型的数值模拟都表明,该方法不仅大大降低了计算量,提高了正演速度,并且具有良好的稳定性和精确性。     

高精度边界格式的研究

利用有精确解的Ｒｉｎｇｌｅｂ流动，构造了对流场数值解精度进行检验的“Ｒｉｎｇｌｅｂ机器”．重点讨论了边界格式对流场数值解的影响及高精度边界格式的建立．计算表明，在场内应用二阶精度格式情况下，采用二阶精度的边界格式所得到的流场解精度将大大高于采用一阶边界格式所得到的精度．为提高流场解的精度，不仅需要高精度的边界格式，还必须注意边界格式与场内格式的匹配．计算还表明，采用特征线修正的方法能有效地提高边界处理的精度     

基于MUSTA格式的全非线性Boussinesq波浪传播数值模型

建立了求解二维全非线性布氏（Boussinesq）水波方程的有限差分/有限体积混合数值格式. 针对守恒形式的控制方程,采用有限体积方法并结合 MUSTA格式计算数值通量, 剩余项则采用有限差分方法求解, 采用具有总变差减小（totalvariation diminishing, TVD）性质的三阶龙格-库塔法进行时间积分.该格式具备间断捕捉、程序实现简单、数值稳定性强、海岸动边界以及波浪破碎处理方便和可调参数少等优点.利用典型算例对数值模型进行了验证,计算结果与实验数据吻合较好.      

A FULLY NONLINEAR BOUSSINESQ MODELFOR WAVE PROPAGATION BASED ON MUSTA SCHEME

建立了求解二维全非线性布氏（Boussinesq）水波方程的有限差分/有限体积混合数值格式. 针对守恒形式的控制方程,采用有限体积方法并结合 MUSTA格式计算数值通量, 剩余项则采用有限差分方法求解, 采用具有总变差减小（totalvariation diminishing, TVD）性质的三阶龙格-库塔法进行时间积分.该格式具备间断捕捉、程序实现简单、数值稳定性强、海岸动边界以及波浪破碎处理方便和可调参数少等优点.利用典型算例对数值模型进行了验证,计算结果与实验数据吻合较好.     

奇异边界法求解多连通域的瞬态热传导问题

提出一种基于奇异边界法结合双重互易法的数值模型来求解瞬态热传导问题。奇异边界法属于配点型边界无网格方法,相对于网格方法,其具有无需划分网格,只需边界配点的优势。运用差分格式来处理热传导方程中的时间变量,将原热传导方程化为非齐次修正Helmholtz方程。修正Helmholtz方程的解由齐次解和特解两部分组成,齐次解通过奇异边界法求出,特解由双重互易法求出,源项由径向基函数近似。通过数值算例检验了本文数值模型的精度及有效性;算例结果表明,该数值模型计算精度较高,误差基本都在1%以内,具有很好的稳定性,能有效地应用于求解多连通域的瞬态热传导问题。     

基于时间依赖基本解的奇异边界法模拟二维狄利克雷边界标量波方程

奇异边界法是一个半解析边界配点强格式方法,具有无数值积分和无网格、编程容易以及数学简单等优点。本文首次将时间依赖基本解运用于奇异边界法,计算模拟二维标量波方程;结合确定源点强度因子的反插值技术,提出了二维狄利克雷边界标量波方程源点强度因子的一个经验公式;引进了解决波方程基本解G奇异性的一种无奇异积分处理方法。数值实验证明,基于时间依赖基本解的奇异边界法可精确高效地模拟二维狄利克雷边界标量波方程,在计算效率、精度、稳定性和适应性等方面有明显优势。     

网格与高精度差分计算问题

研究NS方程差分求解时来流雷诺数、计算格式精度和计算网格之间的关系．给出了判定空间三个方向上的粘性贡献在给定雷诺数、格式精度和网格下是否能够正确计入的估计方法．指出在NS方程的二阶差分方法的数值模拟中,由于物面法向采用了压缩网格技术,物面附近的网格间距很小,该方向上的粘性贡献可被计入．但是如果流向和周向的网格较粗,相应的差分方程中的粘性贡献可能落入截断误差相同的量级,因此在精度上等于仍是求解略去流向和周向粘性项的薄层近似方程．指出,高阶精度的差分计算格式,可以避免对网格要求苛刻的困难．并进一步讨论了建立高阶精度格式的问题,提出了建立高阶精度格式应该满足的原则：耗散控制原则以及色散控制原则．为了避免激波附近可能出现的微小非物理振荡,建议发展混合高阶精度格式,即在激波区,采用网格自适应的NND格式,在激波以外的区域,采用按上述原则发展的高阶格式．     

基于隐式扩散的直接力格式浸没边界格子Boltzmann方法

采用浸没边界格子Boltzmann (immersed boundary-lattice Boltzmann, IB-LB)模型执行动边界绕流数值模拟时,信息交互界面和边界力计算格式直接影响流动求解器的数值精度和计算效率.基于隐式扩散界面,一种改进的直接力格式IB-LB模型被提出.边界力表达式基于欧拉/拉格朗日变量同一性准则推导,转换矩阵描述的信息交互界面耦合了拉格朗日节点间的非同步运动.采用Richardson迭代数值求解关联边界力与无滑移速度约束的线性方程组,不仅克服了传统速度修正格式中矩阵求逆引起的计算效率问题,而且摆脱了算法稳定性对拉格朗日点分布的依赖.根据解析解已知的Taylor-Green涡流评估本文模型的数值模拟精度,结果表明改进的IB模型能够完整保留背景LB模型的二阶数值精度.静止圆柱和振荡圆柱绕流数值实验结果表明,当前模型在涉及复杂外形和运动界面的流动模拟中能够提供可靠的数值预测,满足力同一性的IB-LB模型能够有效抑制非定常流体力的伪物理震荡.波动翼型绕流模拟验证了当前模型的实用性,可在大变形柔性体流固耦合动力学问题中进一步推广.     

一种求解反应扩散方程（组）的初边值问题的快速数值方法

本文提出了一种数值求解反应扩散方程（组）的快速差分方法。差分近似的总截断误差为Ｏ（△ｘ＾２，△ｔ＾２，（△ｔ／△ｘ）＾２）。相容性、稳定性分析表明此差分格式能够有效地抑制非线性项带来的刚性效应和短波不稳定性。对１维和２维的Ｆｉｔｚｈｕｇｈ－Ｎａｇｕｍｏ（ＦＨＮ）方程的数值模拟表明此方法具有计算速度快，稳定性好（△ｘ，△ｔ）的稳定区域较宽，晚皇在高维空间中推广等优点，是数值求解反应扩散方程较好的数值     

The three-point fifth-order accurate generalized compact scheme and its applications

A three-point fifth-order accurate generalized compact scheme (GC scheme) with a spectral-like resolution is constructed in a general way. The scheme satisfies the principle of stability and the principle about suppression of the oscillations, therefore numerical errors can decay automatically and no spurious oscillations are generated around shocks. The third-order TVD type Runge-Kutta method is employed for the time integration, thus making the GC scheme best suited for unsteady problems. Numerical results show that the GC scheme is shock-capturing. The time-dependent boundary conditions proposed by Thompson are well employed when the algorithm is applied to the Euler equations of gas dynamics. The project supported by the National Natural Science Foundation of China (19972038) and Foundation of the National CFD Laboratory of China     

RESTUDY OF COUPLED FIELD THEORIES FOR MICROPOLAR CONTINUA (Ⅰ)──MICROPOLAR THERMOELASTICITY

ＩｎｔｒｏｄｕｃｔｉｏｎＳｏｍｅａｕｔｈｏｒｓｓｔｕｄｉｅｄｔｈｅｃｏｕｐｌｅｄｆｉｅｌｄｐｒｏｂｌｅｍｓｆｏｒｍｉｃｒｏｐｏｌａｒｃｏｎｔｉｎｕａ .Ｅｓｐｅｃｉａｌｌｙ ,Ｗ .Ｎｏｗａｃｋｉｐｕｂｌｉｓｈｅｄａｓｅｒｉｅｓｏｆａｂｏｕｔ 4 0ｓｃｉｅｎｔｉｆｉｃｐａｐｅｒｓｄｅａｌｉｎｇｗｉｔｈｔｈｅｍｉｃｒｏｐｏｌａｒｔｈｅｒｍｏｅｌａｓｔｉｃｉｔｙａｓｗｅｌｌａｓｔｈｅｐｒｏｂｌｅｍｓｏｆｄｉｓｔｏｒｔｉｏｎ ,ｔｈｅｒｍｏｄｉｆｆｕｓｉｏｎ ,ｔｈｅｒｍｏｐｉｅｚｏｅｌｅｃｔｒｉｃｉｔｙａｎｄｍ…     

A NEW IHIERARCHICAL BOUNDARY ELEMENT METHOD AND ITS ADAPTIVE PROCESSES FOR PLATE BENDING PROBLEM

In this paper a new hierarchical boundary element method is introducedfor solving the problem of plate bending.Exact solutions of the governing equationsare used inside the domain together with independent deflections and rotations on theboundary.A generalized variational principle is employed to achieve the boundaryelement formulation.Further,the adaptive processes of the method are also discussed.By virtue of the error estimate technique proposed by Zienkiewicz et al.areasonable error indicator and adaptive scheme are suggested.Numerical examplesillustrate the high accuracy of the new elements and show the excellent efficiency ofthe adaptive computation described in this paper.     

FE analysis of stress and strain fields in finite random composite bodies

In this paper the mechanical behaviour of finite random heterogeneous bodies is considered. The analysis of non-local interactions between heterogeneities in microscopically heterogeneous materials is necessary when the spatial variation of the load or the dimensions of the body, relative to the scale of the microstructure, cannot be ignored. Microstructures can be periodic but generically they are random. In the first case, an exact calculation can be performed but in the second case recourse has to be made either to simulation or to some scheme of approximation. One such scheme is based on a stochastic variational principle. The novelty of the present work is that a stochastic variational principle is projected directly onto a finite-element basis so that all subsequent analysis is performed within a finite-element framework. The proposed formulation provides expressions for the local stress and strain fields in any realization of the medium, from which expressions for statistically-averaged quantities can be derived. Then an approximation of Hashin-Shtrikman type is developed, which generates a FE-based numerical procedure able to take account of interactions between random inclusions and boundary layer effects in finite composite structures. Finally, two examples are presented, namely a cylinder with square cross-section subjected to mixed boundary conditions of different types on different faces and a rectangular body containing a centre crack. The results show that in the vicinity of the boundary or close to the crack tip, the strain and the stress in the matrix and in the inclusions differ considerably from those obtained by the formal application of conventional homogenization.     

Restudy of coupled field theories for micropolar continua (I)—Micropolar thermoelasticity

Problems of micropolar thermoelasticity have been presented and discussed by some authors in the traditional framework of micropolar continuum field theory. In this paper the theory of micropolar thermoelasticity is restudied. The reason why it was restricted to a linear one is analyzed. The rather general principle of virtual work and the new formulation for the virtual work of internal forces as well as the rather complete Hamilton principle in micropolar thermoelasticity are established. From this new Hamilton principle not only the equations of motion, the balance equation of entropy, the boundary conditions of stress, couple stress and heat, but also the boundary conditions of displacement, microrotation and temperature are simultaneously derived. Contributed by DAI Tian-min Foundation item: the National Natural Science Foundation of China (10072024); the International Cooperation Project of the NSFC (10011130235) and the DFG (51520001); the Research Foundation of Liaoning Education Committee (990111001) Biography: DAI Tian-min (1931-)     

On deforming of a linear viscoelastic orthotropic half-space under the turning rigid cylindrical roller

In this paper we consider the problem of rigid cylinder turning on a linear viscoelastic orthotropic half-space with Coulomb's friction acting along the contact area. Results for extents of contact area and pressure under the cylinder are obtained using Volterra's principle. The obtained functions of viscoelastic operators are interpreted by a method based on expansion of such functions in operator continued fractions. A solution is given for the general type of resolvent viscoelastic operators expressing rheological properties of half-space material. Algebra of resolvent Volterrian operators is used to facilitate the calculations. An example is given to illustrate the results for real viscoelastic material with the rheological properties expressed by the operators of Yu.N. Rabotnov.     

A finite volume scheme preserving extremum principle for convection–diffusion equations on polygonal meshes

We propose a nonlinear finite volume scheme for convection–diffusion equation on polygonal meshes and prove that the discrete solution of the scheme satisfies the discrete extremum principle. The approximation of diffusive flux is based on an adaptive approach of choosing stencil in the construction of discrete normal flux, and the approximation of convection flux is based on the second‐order upwind method with proper slope limiter. Our scheme is locally conservative and has only cell‐centered unknowns. Numerical results show that our scheme can preserve discrete extremum principle and has almost second‐order accuracy. Copyright © 2017 John Wiley & Sons, Ltd.     

Numerical analysis of mixed convection in three‐dimensional rectangular channel with flush‐mounted heat sources based on field synergy principle

In this paper the fluid flow and heat transfer characteristics of mixed convection in three‐dimensional rectangular channel with four heat sources are investigated numerically. The SIMPLEC algorithm is applied to deal with the coupling between pressure and velocity, and a new high‐order stability‐guaranteed second‐order difference (SGSD) scheme is adopted to discretize the convection term. The influence of four parameters is studied: Richardson number, heat source distribution, channel height and inclination angle. The numerical results are analysed from the viewpoint of the field synergy principle, which says that the enhanced convective heat transfer is related not only to the velocity field and temperature field, but also to the synergy between them. It is found that the effects of the four parameters on the thermal performance can all be explained with the field synergy principle. To obtain better electronic cooling, the synergy between the velocity and temperature gradient should be increased when other conditions are unchanged. Copyright © 2006 John Wiley & Sons, Ltd.     

Generalized hydrodynamical equations obeying the entropy principle

Three generalizations of classical hydrodynamic theories that are compatible with equilibrium thermodynamics and that are suitable for an appropriate macroscopic dynamical theory of polymeric liquids are considered. The strain tensor, the stress tensor and the chain segment distribution function (introduced in the network theory of polymeric liquids) are accepted as new state variables. We find that the generalized hydrodynamic equations are compatible with equilibrium thermodynamics provided certain conditions restricting the freedom of choice of constitutive relations are satisfied. In some particular cases the conditions are known from other considerations. We say that a dynamical theory is compatible with equilibrium thermodynamics, or equivalently, that it obeys the entropy principle if the properties listed in section 2.1 are satisfied.     

Further study of the equivalent theorem of Hellinger-Reissner and Hu-Washizu variational principles

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎ１９５４,ｂｙｔｈｅｓｏ＿ｃａｌｅｄｔｒｉａｌ＿ａｎｄ＿ｅｒｏｒｍｅｔｈｏｄＨｕ［１］ｄｅｄｕｃｅｄｔｈｅｗｅｌ＿ｋｎｏｗｎＨｕ＿Ｗａｓｈｉｚｕｐｒｉｎｃｉｐｌｅ,ｗｈｉｃｈｐｌａｙｓａｎｉｍｐｏｒｔａｎｔｒｏｌｅｉｎｔｈ...