一种模拟动边界绕流的锐利界面浸入边界法

When the structural wall moves over a fixed grid, the structure coverage will change, resulting in many dead and emerging elements. To avoid the influence of malformation and reconstruction of body-fitted grids on the calculation efficiency and accuracy of the fluid-structure interaction problems with coupled boundary movement on the fixed grid, an improved numerical method for describing the interaction between an immersed rigid body and fluid based on a sharp-interface is proposed. In this method, both the fluid and solid are regarded as pure fluid domains in the whole computational domain, and the solid boundary is divided into several Lagrangian grid points. The flow parameter or velocity is reconstructed by interpolation at the interface element, which is then directly used as the boundary condition of the flow field, thus reflecting the influence of the wall boundary conditions. The method constructs the calculation structure of “virtual point, force point and vertical foot point”, and the velocity of the virtual point is obtained by bilinear interpolation. Then, the velocity of the force point is calculated by forcing the solid boundary to meet the no-slip condition, and the equations of the coupling system based on the immersion boundary method are finally solved to realize the numerical simulation of the flow with a complex moving boundary. The numerical program for this immersed boundary method is established using C++, then the accuracy and reliability of the proposed method are validated by comparison with the literature and experimental results of the basic numerical example of flow around a cylinder. Furthermore, the effects of the structural shape and the angle of attack on the trailing vortex structure, the vortex shedding frequency, and the lift/ coefficient characteristics of the flow around the elliptical cylinder have been analyzed. The anti-symmetric S-type, “P+S” Ⅰ-type and “P+S” Ⅱ-type trailing vortex shedding modes, as well as the variation laws of the vortex structure size, vortex shedding frequency and lift-drag coefficients ratio with axis ratio and angle of attack, are captured. The critical angle of attack (25°) corresponding to the maximum lift-drag ratio is determined as 25°.     

模拟相变问题的精细积分边界元法

将精细积分边界元法和界面追踪法相结合求解相变问题。因为边界元法只需要将待求解空间域的边界离散,方便连续追踪移动界面位置和重构网格,所以边界元法适合应用于移动边界问题的模拟。首先,利用精细积分边界元法在固相区域和液相区域分别求解相应的瞬态热传导控制方程,从而求得温度场和边界热流密度。然后,根据固-液相变界面上的能量平衡方程,利用热流密度求得相变界面的移动速度,再采用界面追踪法预测移动相变界面的位置变化。最后,给出了几个数值算例,并通过与参考解的对比验证本文方法的准确性。     

维生素 E型人工膝关节在ISO标准步态与上楼梯运动下的磨损情况对比

为研究人工膝关节假体在ISO标准步态与上楼梯运动负载下的摩擦磨损情况,通过对日常活动频次及动力学和运动学的研究与分析,基于人工膝关节磨损试验机,以新型Vitamin E型人工膝关节为对象进行500万次标准步态磨损测试加200万次上楼梯运动磨损测试. 结果表明正常步态阶段的磨损率为4.34 mg/mc,上楼梯运动阶段的磨损率为5.95 mg/mc. 上楼梯运动相比标准步态对半月板衬垫的磨损量造成更大的影响,半月板衬垫表面形貌相较标准步态磨损后表面形貌也发生了显著变化,出现了条纹图案和局部凹陷,这种表面形貌结果更接近于临床半月板衬垫取出物的表面形貌,目前对于人工膝关节假体的标准磨损测试方式不够全面,需要进一步改进整个测试的方案才能更好地体现体外磨损测试的意义.      

集中载荷作用下不同边界条件的弯曲厚矩形板的受迫振动

运用边界积分法研究了四边简支、两对边固定另两对边简支、四边固定三种复杂边界条件下厚矩形板的受迫振动问题,求解过程清晰,从而给出了受迫振动控制方程和挠曲面方程。通过在Matlab平台上进行数值计算,得出了图表形式的计算结果,并与有限元模拟值进行对照。研究表明,边界积分法用于求解厚矩形板的受迫振动问题的准确性,本文推导的控制方程和挠曲面方程的正确性,进而对工程实际中的各种相关问题具有一定的现实意义,也为求解此类问题提供了一种新途径,可以直接运用到工程实际中。     

边界振动的微腔中薛定谔猫态的产生

采用全量子理论,考虑了边界振动的微腔(振动边界视作频率为ωm的量子谐振子)与单模辐射场构成的系统,在薛定谔绘景中,给出了系统的时间演化算符,得到了系统的态函数随时间的演化关系.结果表明,当振动边界回到初态时的腔场态能给出许多与已有文献不同的薛定谔猫态,且通过调整单模辐射场的频率可以选择想要的薛定谔猫态.     

New bifurcations of basin boundaries involving Wada and a smooth Wada basin boundary

This paper demonstrates and analyses double heteroclinic tangency in a three-well potential model, which can produce three new types of bifurcations of basin boundaries including from smooth to Wada basin boundaries, from fractal to Wada basin boundaries in which no changes of accessible periodic orbits happen, and from Wada to Wada basin boundaries. In a model of mechanical oscillator, it shows that a Wada basin boundary can be smooth.     

Numeral eddy current sensor modelling based on genetic neural network

This paper presents a method used to the numeral eddy current sensor modelling based on the genetic neural network to settle its nonlinear problem. The principle and algorithms of genetic neural network are introduced. In this method, the nonlinear model parameters of the numeral eddy current sensor are optimized by genetic neural network （GNN） according to measurement data. So the method remains both the global searching ability of genetic algorithm and the good local searching ability of neural network. The nonlinear model has the advantages of strong robustness, on-line modelling and high precision. The maximum nonlinearity error can be reduced to 0.037% by using GNN. However, the maximum nonlinearity error is 0.075% using the least square method.     

Bound entangled states invariant under Ux

This paper obtains an entangled condition for isotropic-like states by using an atomic map. It constructs a class of bound entangled states from the entangled condition and shows that the partial transposition of the state from the constructed bound entangled class is an edge bound entangled state by using range criterion.