完全耦合层边界条件在圆柱绕流DNS中的应用

数值模拟圆柱绕流问题是计算流体力学和计算气动声学中的一个经典问题,要准确模拟这一类问题,远场无反射边界条件是一个关键环节。本文采用适用于主控方程为Navier-Stokes方程的完全耦合层无反射边界条件对低雷诺数条件下圆柱绕流进行直接数值模拟,研究了计算域大小的选取对数值结果的影响。数值模拟研究表明:(1)由完全耦合层边界条件产生的相对反射误差很小;(2)如果仅是为了获得圆柱附近的流场和声场信息,采用完全耦合层边界条件可以选用较小的计算域;(3)适当增加计算域的上下边界宽度,可以很大程度地减小右边界的宽度,从而在保证精度的情况下减小计算量。     

80K分置式斯特林制冷机连接软管动力特性的模拟与验证

建立了 80 K分置式斯特林制冷机连接软管的理论模型 ,通过对控制方程进行处理及合理简化 ,对边界条件的处理引入特征线法 ,通过数值计算对控制方程进行了求解 ,通过实验对理论模拟结果进行了验证 ,证明了本模拟方法具有良好的准确性。     

考虑空气柱的水力旋流器内流场的数值模拟

对水力旋流器内流场进行了数值模拟.旋流器内液相和气相的湍流模拟采用了雷诺应力模型(RSM),边界条件采用速度进口、压力出口.在同一工况下选用三维模型及二维简化模型、分别用VOF模型及刚盖法对空气柱进行处理,就空气柱形成、切向及轴向速度与相同条件下试验测量的结果进行比较.结果表明,三维模型计算结果较二维简化模型与实验值吻合得更好,但二维简化模型在降低精度不大的前提下大大减少了计算时间.对空气柱的处理,选用VOF模型得到的结果更准确,而刚盖法减少了计算用时.     

地铁火灾CFD模拟时边界条件的研究

本文对地铁车站及隧道内火灾通风排烟过程中的空气流动特性进行了数值模拟.对模拟边界提出了简化模型,即当某一个车站发生火灾时,计算区域可以只考虑其两端各一个车站和各两个区间;当某一区间发生火灾时,计算区域只需考虑该区间两端各两个车站和各一个区间.上面两种情况下,所有出入口均可假定为自由边界条件.模拟结果与现场实验结果基本一致.     

流体力学方程的间断有限元方法

在二维区域三角形网格上应用一阶、二阶和三阶精度间断有限元方法,对流体力学方程和方程组进行了数值模拟.计算结果与差分方法计算结果比较,认为间断有限元方法在求解复杂边界条件和区域问题上有一定的优势.     

方柱绕流大涡模拟

采用有限体/有限元混合格式、非结构网格和大涡模拟方法求解可压缩的N-S方程,对Re=22 000的方柱绕流进行数值模拟,并对不同的边界条件进行详细的分析比较.通过对以往研究经验的总结和利用精细的边界条件,使得采用二阶精度的数值格式和较稀疏的网格仍然得到了令人满意的计算结果,甚至优于以往采用密网格的模拟结果.     

时域有限差分法中电磁场的吸收边界条件

吸收边界条件在时域有限差分法计算电磁场中有着广泛的应用,从吸收边界所要满足的条件出发,将Lindman吸收边界条件推广,得到了便于计算的简化形式.两维数值模拟结果表明:以各种角度向左或右传播的p极化和s极化电磁波在边界面上的反射率均低于0.5%.     

三维空腔辐射场计算的一种简化数值方法

给出三维空腔辐射场计算的一种简化数值方法。在一些物理假定条件下,把问题简化为解反照边界条件的 Boltzmann方程。对方法作了详细描述,讨论了迭代过程的收敛性,最后给出了数值例子。     

风机三维气动性能数值分析研究

为了提高风机性能数值计算的准确性,本文应用Fluent软件对低压轴流风机整机性能进行数值模拟,分析不同湍流模型条件下,风机气动性能的差别。同时,通过合理建立数值模型、设置边界条件、选取与实际工况一致的风机静压来提高数值模拟精度和可靠性。最后,针对不同工况下数值模拟的收敛判定提出建议。本文研究结果将对聚能型风力机的数值模拟研究起指导性作用。     

Rijke管自激热声振荡的数值模拟

为了对热声不稳定的发生及控制机理进行研究,对Rijke管内的自激热声振荡现象进行了数值模拟。采用具有低频散低耗散特点的计算气动声学方法,对带有非线性热源项的声波方程进行数值求解,并比较了不同的热源模型及边界条件对非线性效应的影响。结果表明,计算气动声学方法可以成功捕捉到Rijke管内压力的起振过程,而且在速度扰动达到平均流速度的1/3时,振荡会由线性增长转为非线性增长,最终达到有限幅值极限循环。相比热源项,考虑管口辐射耗散的非线性边界条件在振荡幅值和频谱方面对结果的影响都比较小。数值模拟得到的结果与实验符合较好,表明计算气动声学方法适合于热声振荡问题的研究。     

N=2 boundary supersymmetry in integrable models and perturbed boundary conformal field theory

Boundary integrable models with N=2 supersymmetry are considered. For the simplest boundary N=2 superconformal minimal model with a Chebyshev bulk perturbation we show explicitly how fermionic boundary degrees of freedom arise naturally in the boundary perturbation in order to maintain integrability and N=2 supersymmetry. A new boundary reflection matrix is obtained for this model and N=2 boundary superalgebra is studied. A factorized scattering theory is proposed for a N=2 supersymmetric extension of the boundary sine-Gordon model with either (i) fermionic or (ii) bosonic and fermionic boundary degrees of freedom. Exact results are obtained for some quantum impurity problems: the boundary scaling Lee–Yang model, a massive deformation of the anisotropic Kondo model at the filling values g=2/(2n+3) and the boundary Ashkin–Teller model.     

粘接界面的非线性弹簧模型及实验验证

研究了纵波及横波垂直入射时在粘接界面处的非线性边界条件。在二阶微扰近似及粘接界面层厚度远小于入射声波长的条件下,给出了严格的界面非线性弹簧模型的表达式。数值计算结果表明,利用该模型得到的透射二次谐波与利用连续性边界条件得到的精确解吻合良好。针对钢(铝)的楔形粘接结构,在纵波垂直入射的条件下对该非线性弹簧模型进行了实验验证。实际测量结果显示该模型具有较高的精度;且被粘物与胶层的阻抗比越大,该模型的误差越小。      

Theoretical and experimental research on nonlinear spring models of a bonding interface

<正>The nonlinear boundary conditions on a bonding interface between two solids of a longitudinal or shear horizontal(SH) wave under normal incidence were explored.By applying the second-order perturbation method,the nonlinear spring models are rigorously developed in the limit of small thickness to wavelength ratio.Numerical results agree well with the exact solutions obtained by continuous boundary conditions.The nonlinear spring model for longitudinal wave is verified by measuring the nonlinearity of a wedge-shaped bonding structure of steel or aluminum substrates.The experimental results reveal that the above model is rather accurate and as the impedance ratio of adherend to adhesive increases,the model becomes more accurate.     

Three-dimensional finite-element simulation of Zener pinning dynamics

The Zener pinning dynamics of a moving boundary interacting with one or more particles is described by a three-dimensional (3D) finite-element model. The model, based upon a variational formulation for boundary motion by viscous drag, is solved by a finite-element method to obtain the velocity at each node of triangular linear elements on the grain boundary. It is first applied to relatively simple and validated cases, for which analytical and numerical results are available. These cases correspond to an axisymmetrical geometry, in which the grain boundary interacts with a centred particle. A simple analytical pinning criterion is derived from these simulations. The model is then applied to general 3D cases, in which the grain boundary interacts with arbitrarily localized and sized particles. The aim of these 3D simulations is to quantify the influence of the position and the number of particles on the average grain-boundary velocity. It is shown, for example, that the drag effect is enhanced when the particle, or the cluster of particles, is off-centre and that pinning is less efficient with several particles than with a single particle producing the same Zener force.     

ZA27合金晶界处铁、稀土元素的有序化与交互作用

为从理论上揭示铁、稀土元素在锌铝合金晶界处的行为本质，建立了ZA27合金中α相大角度重位点阵晶界模型，利用递归法(Recursion)计算了晶界的电子结构(状态密度、费米能级、结构能).用晶界结构能定义合金的团簇能(有序能)，并计算了偏聚铁及稀土晶界的团簇能.计算结果表明：铁、稀土元素在锌铝合金晶界处团簇能为正值，不能形成团簇，具有有序化倾向，趋于形成稳定的金属间化合物.铁与稀土元素在晶界形成负电中心，降低晶界的费米能级. 关键词： 稀土 晶界 电子结构 有序化     

A simple enthalpy-based lattice Boltzmann scheme for complicated thermal systems

To extend the lattice Boltzmann (LB) method to describe the applicable energy systems, the first key step is to build a suitable thermal LB model and corresponding boundary treatments. There are two main shortcomings in the existing related works: either some additional energy source terms are inconvenient to be naturally incorporated or the implementation of non-Dirichlet-type thermal boundary conditions is extremely difficult and sometimes impossible in them for complicated thermal systems, which restrict their applicability to only a few special classes of problems. In order to overcome these drawbacks by a simple way, in this paper a thermal LB model and corresponding boundary treatments are constructed based on the total enthalpy. The specific benefits due to the introduction of the total enthalpy are analyzed and it is found that the numerical results obtained by the present scheme agree well with the analytical solutions and/or the data reported in previous studies.     

基于Rudin-Osher-Fatemi模型的图像除模糊和除噪音新模型

介绍一种依赖时间的新模型来解决图像除噪音和除模糊问题.分别使用逆反射、中值两种边界条件的数值试验比较本文新模型和Rudin-Osher-Fatemi模型.试验结果表明中值边界条件的误差比逆反射边界条件误差小.     

A unified framework for the Kondo problem and for an impurity in a Luttinger liquid

We develop a unified theoretical framework for the anisotropic Kondo model and the boundary sine-Gordon model. They are both boundary integrable quantum field theories with a quantum-group spin at the boundary which takes values, respectively, in standard or cyclic representations of the quantum groupSU(2) _q. This unification is powerful, and allows us to find new results for both models. For the anisotropic Kondo problem, we find exact expressions (in the presence of a magnetic field) for all the coefficients in the Anderson-Yuval perturbative expansion. Our expressions hold initially in the very anisotropic regime, but we show how to continue them beyond the Toulouse point all the way to the isotropic point using an analog of dimensional regularization. The analytic structure is transparent, involving only simple poles which we determine exactly, together with their residues. For the boundary sine-Gordon model, which describes an impurity in a Luttinger liquid, we find the nonequilibrium conductance for all values of the Luttinger coupling. This is an intricate computation because the voltage operator and the boundary scattering do not commute with each other.     

Grain boundary pinning in two-phase materials

The calculation of the Zener stress involves two parts (1) the determination of the force between a particle and a grain boundary, and the shape of the resulting dimple on the grain boundary and (2) the statistical problem of summing the forces over the boundary. Calculations of the forces from spherical, needle-shaped and plate shaped particles are used to estimate the Zener stress by means of a point obstacle model and statistical dimple model. The results are compared with experimental determinations of stable grain sizes, with a computer simulation growth and with experimental measurements of retarded primary recrystallisation. When the particles are non-spherical or inhomogeneously distributed the Zener stress is anisotropic; this effect has been observed experimentally.     

晶粒生长演变相场法模拟界面表达的物理模型

讨论了当前固态组织演变过程的相场法模拟模型,论证了相场法中界面的概念以及模型中界面的各种处理方法,以AZ31镁合金再结晶系统为例,研究了模型参数取值对界面特征的影响,提出了晶界作用域的概念,阐述了晶界作用域就是相场模型中界面处有序化变量的变化范围,其物理意义是界面能量的分布范围,并对应于成分界面偏析的范围.模拟得出,晶界作用域宽度主要由梯度项系数决定,晶界能则由梯度项系数和耦合项系数共同决定.对于AZ31镁合金,模拟研究了晶界作用域宽度取值的合理性和对显微组织影响的关系,得出取值为1.18 μm时,模拟符 关键词： 相场法 界面 计算机模拟 显微组织