A mixed hp FEM for the approximation of fourth‐order singularly perturbed problems on smooth domains

We consider fourth‐order singularly perturbed problems posed on smooth domains and the approximation of their solution by a mixed Finite Element Method on the so‐called Spectral Boundary Layer Mesh. We show that the method converges uniformly, with respect to the singular perturbation parameter, at an exponential rate when the error is measured in the energy norm. Numerical examples illustrate our theoretical findings.     

The local ultraconvergence of high‐order finite element method for second‐order elliptic problems with constant coefficients over a rectangular partition

In this article, we will discuss the local ultraconvergence of high‐degree finite element method based on a rectangular partition for the second‐degree elliptic problem with constant coefficients in Ω ? ?² , u( y ) = 0 on ?Ω . Based on suitable regularity, ultraconvergence of the displacement of the extrapolated kth (k ≥ 3) degree finite element solution has been obtained by an extrapolation technique. Finally, numerical experiments are applied to demonstrate our theoretical findings.     

Evaluation of RBFs for volume data interpolation in CFD

A greedy method for choosing an optimum reduced set of control points is integrated with RBF interpolation and evaluated for the purpose of interpolating large‐volume data sets in CFD. Given a function defined at a set of points, the greedy method selects a small subset of these points that is sufficient to keep the interpolation error at all the remaining points below a chosen bound. This is equivalent to a type of data compression and would have useful storage, post‐processing, and computational applications in CFD. To test the method in terms of both the point selection scheme and the suitability of reduced control point volume interpolation, a trial application of the interpolation to velocity fields in CFD volume meshes is considered. To optimise the point selection process, and attempt to be able to capture multiple length scales, a variable support radius formulation has also been included. Structured and unstructured mesh cases are considered for aerofoils, a wing case and a wing‐body case. For smooth volume functions, the method is shown to work well, producing accurate velocity interpolations using a very small number of the cells in the mesh. For general complex fields including large gradients, the method is still shown to be effective, although large gradients require more interpolation points to be used.Copyright © 2013 John Wiley & Sons, Ltd.     

Search for a Steric Penultimate Effect: Copolymerization of 2,3,4-Trimethyl-3-pentyl Methacrylate with Styrene

The hindered monomer, 2,3,4-trimethyl-3-pentyl methacrylate (I), was synthesized for penultimate effect studies. Since it readily homopoiymerized (km¹¹¹≠ 0) and readily copolymerized with styrene, copolymerizations of I with styrene were carried out at 60°C in benzene with AIBN as initiator. The conversion to copolymer and the copolymer composition were determined by using GLC techniques. Composition-conversion data was analyzed by performing a computerized nonlinear least-squares fitting to the integrated form of the penultimate model equation. The experimental design included the use of optimized M¹°/M²° ratios. The penultimate reactivity ratios calculated from these data were r¹′ = 0.23, r¹′= 0.59, r² = 0.59, r²′ = 1.34. Thus, when I is the penultimate unit, a terminal styryl radical prefers to add styrene, whereas when styrene is the penultimate unit, terminal styryl radicals prefer to add I. These results constitute the best evidence for a steric penultimate effect yet available in the literature from composition-conversion studies. However, the case is not yet proved. Further studies to strengthen this conclusion are proposed.     

How a moving passive observer can perceive its environment ? The Unruh effect revisited

We consider a point-like observer that moves in a medium illuminated by noise sources with Lorentz-invariant spectrum. We show that the autocorrelation function of the signal recorded by the observer allows it to perceive its environment. More precisely, we consider an observer with constant acceleration (along a Rindler trajectory) and we exploit the recent work on the emergence of the Green’s function from the cross correlation of signals transmitted by noise sources. First we recover the result that the signal recorded by the observer has a constant Wigner transform, i.e. a constant local spectrum, when the medium is homogeneous (this is the classical analogue of the Unruh effect). We complete that result by showing that the Rindler trajectory is the only straight-line trajectory that satisfies this property. We also show that, in the presence of an obstacle in the form of an infinite perfect mirror, the Wigner transform is perturbed when the observer comes into the neighborhood of the obstacle. The perturbation makes it possible for the observer to determine its position relative to the obstacle once the entire trajectory has been traversed.     

A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

Evolution of vortices in the wake of an ARJ21 airplane: Application of the lift-drag model

Wake separation is crucial to aircraft landing safety and is an important factor in airport operational efficiency. The near-ground evolution characteristics of wake vortices form the foundation of the wake separation system design. In this study, we analysed the near-ground evolution of vortices in the wake of a domestic aircraft ARJ21 initialised by the lift-drag model using large eddy simulations based on an adaptive mesh. Evolution of wake vortices formed by the main wing, flap and horizontal tail was discussed in detail. The horizontal tail vortices are the weakest and dissipate rapidly, whereas the flap vortices are the strongest and induce the tip vortex to merge with them. The horizontal tail and flap of an ARJ21 do not significantly influence the circulation evolution, height change and movement trajectory of the wake vortices. The far-field evolution of wake vortices can therefore be analysed using the conventional wake vortex model.     

基于移动可变形孔洞方法的超弹性结构拓扑优化

现有隐式拓扑优化方法在进行超弹性结构拓扑优化设计时,具有设计变量多、中间设计有限元分析存在严重的收敛性和设计结果无法直接导入CAD/CAE系统等问题。为解决这些问题,提出了一种基于移动可变形孔洞的显式拓扑优化方法来进行承受大变形的超弹性结构设计,材料本构采用常用的Mooney-Rivlin模型。首先,介绍了移动可变形孔洞方法的基本思想和可变形孔洞的显式描述方法;其次,构造了基于移动可变形孔洞方法的超弹性结构拓扑优化的数学列式,给出了相应的灵敏度结果;最后,通过数值算例验证了本方法的有效性。数值结果表明,该方法可以通过较少的设计变量和非常稳健的优化过程,给出边界由B样条曲线描述且可与CAD/CAE软件无缝连接的超弹性结构设计。     

光学透明频率选择表面的设计研究

以氟化镁为基底材料, 采用基底、金属网栅与频率选择表面一体化设计方法设计了一种雷达波与光学波段双带通的结构. 利用模式匹配法对设计结构的传输特性进行了仿真研究, 并将设计结果与制备样件的测试结果进行了对比分析. 结果发现: 采用一体化设计的方法设计光学透明频率选择表面, 不仅能够快速得到电场基函数而且还能够准确预估其谐振尺寸, 从而在提高计算效率的同时避免了模式互作用零点的出现. 采用一体化设计方法获得了具有稳定滤波特性的光学透明频率选择表面, 为雷达/红外双模制导头罩的电磁屏蔽技术和隐身技术提供了一种有效技术方案. 关键词： 复合制导 频率选择表面 金属网栅 一体化设计     

求解带自由面渗流问题移动网格法的改进及其工程应用

采用有限元固定网格和移动网格相结合的方法,求解带自由面的坝体稳定渗流问题。本方法是针对一般坝体的结构,在有限元网格生成之后,又自动生成了网格移动所需的信息,实现了网格在迭代计算中随自由面的变化自动均匀移动;讨论了网格奇异的不同类型,提出了避免单元形状怪异的移动网格方法;在欧拉空间中分区记录坝体不同的渗透性能,在网格移动到新位置后,自动识别网格所在的物性空间,从而确定单元在新位置上的渗流系数。本文将该方法应用到吉林丰满混凝土重力坝防渗加固的实际工程当中,在考虑大坝原有防渗措施的前提下,利用数值分析,研究了防渗芯墙对大坝渗流性能的影响。计算结果表明,加芯墙后坝体自由面位置明显下移,湿水部分大为减少,坝基扬压力没有明显提升,下游出水点位置没有明显下降。计算过程表明,本文方法及措施具有很好的收敛性与稳定性,且收敛性与自由面初始位置的选择无关。