Interface‐fitted moving mesh method for axisymmetric two‐phase flow in microchannels

A boundary‐fitted moving mesh scheme is presented for the simulation of two‐phase flow in two‐dimensional and axisymmetric geometries. The incompressible Navier‐Stokes equations are solved using the finite element method, and the mini element is used to satisfy the inf‐sup condition. The interface between the phases is represented explicitly by an interface adapted mesh, thus allowing a sharp transition of the fluid properties. Surface tension is modelled as a volume force and is discretized in a consistent manner, thus allowing to obtain exact equilibrium (up to rounding errors) with the pressure gradient. This is demonstrated for a spherical droplet moving in a constant flow field. The curvature of the interface, required for the surface tension term, is efficiently computed with simple but very accurate geometric formulas. An adaptive moving mesh technique, where smoothing mesh velocities and remeshing are used to preserve the mesh quality, is developed and presented. Mesh refinement strategies, allowing tailoring of the refinement of the computational mesh, are also discussed. Accuracy and robustness of the present method are demonstrated on several validation test cases. The method is developed with the prospect of being applied to microfluidic flows and the simulation of microchannel evaporators used for electronics cooling. Therefore, the simulation results for the flow of a bubble in a microchannel are presented and compared to experimental data.     

变角度纤维复合材料层合斜板的颤振分析

假设纤维方向角沿层合板的长度方向线性变化，研究了变角度纤维复合材料层合斜板的颤振．通过坐标变换将斜板变换为正方形板，采用层合板表面连续变化的速度环量来模拟空气对其的作用，速度环量分布利用Cauchy积分公式计算．建立了系统的Lagrange方程并采用Ritz法得到了层合板的自振频率和颤振/不稳定性分离临界速度．通过数值算例验证了本文模型和方法的正确性和收敛性，分析了各个铺层内纤维方向角的变化对自振频率和颤振/不稳定性分离临界速度的影响．研究结果表明，通过纤维的变角度铺设，可有效地提高层合板的基频和颤振/不稳定性分离临界速度．经合理设计的变角度复合材料层合板具有抑制颤振的作用．     

具任意次非线性项的Camassa-Holm方程的双孤子新解

给出辅助方程、函数变换与变量分离解相结合的方法,构造了具任意次非线性项的Camassa-Holm方程的双孤子和双周期新解.首先,通过两个辅助方程、函数变换与变量分离解,将具任意次非线性项的Camassa-Holm方程的求解问题转化为非线性代数方程的求解问题.然后,借助符号计算系统Mathematica求出该方程组的解,并用辅助方程的相关结论,构造了双周期解和双孤子新解.     

基于离散网格的流形曲面构造综述

基于离散网格的流形曲面构造技术不仅能够生成具有高阶光滑性的曲面,并且该曲面可以是任意拓扑结构的.此外,在构造流形曲面时,无需进行额外的拼接操作,克服了传统曲面造型技术在进行面片之间的拼接时,计算量增大以及曲面光滑性难以保证的难题.本文介绍了流形曲面构造的流程以及构造过程中的难点,然后将目前已有的流形曲面构造技术分为三大类:传统意义上的流形构造方法;基于规范区域的流形构造方法;基于样条曲面推广的流形构造方法.并对每一类都进行详细地分类介绍.最后,对其作一个总结以及对未来的展望.     

On shapes of body waves in periodically inhomogeneous, magnetostrictive, dielectric materials

The shapes of shear body waves in periodically inhomogeneous, magnetostrictive, dielectric media are studied with emphasis on the partial (elastic and magnetostrictive) wave motions coupled to produce magnetoelastic waves __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 7, pp. 57–63, July 2006.     

A fast method for solving fluid–structure interaction problems numerically

The paper presents a semi‐implicit algorithm for solving an unsteady fluid–structure interaction problem. The algorithm for solving numerically the fluid–structure interaction problems was obtained by combining the backward Euler scheme with a semi‐implicit treatment of the convection term for the Navier–Stokes equations and an implicit centered scheme for the structure equations. The structure is governed either by the linear elasticity or by the non‐linear St Venant–Kirchhoff elasticity models. At each time step, the position of the interface is predicted in an explicit way. Then, an optimization problem must be solved, such that the continuity of the velocity as well as the continuity of the stress hold at the interface. During the Broyden, Fletcher, Goldforb, Shano (BFGS) iterations for solving the optimization problem, the fluid mesh does not move, which reduces the computational effort. The term ‘semi‐implicit’ used for the fully algorithm means that the interface position is computed explicitly, while the displacement of the structure, velocity and the pressure of the fluid are computed implicitly. Numerical results are presented. Copyright © 2008 John Wiley & Sons, Ltd.     

An engineering approach to blade designs for low to medium pressure rise rotor-only axial fans

Low to medium pressure rise axial fan equipment of the arbitrary vortex flow rotor-only type is widely used in industrial and commercial applications, with many of the installations and rotor designs being far from optimum. Complex computational methods exist for analyzing flows in, for example, high-speed axial flow compressors with multistage blade rows; however, the designers and manufacturers of low-speed, general-purpose axial flow fan equipment have been reluctant to embrace this technology. A simpler yet reliable design technique is presented that allows this category of ducted axial fan rotors, in the presence of swirl-free inlet flow, to be designed to achieve a specified duty with sufficient accuracy for engineering purposes. Practical blade design recommendations and limits, similar to those that exist for free vortex flow axial rotors, have been established for the arbitrary vortex flow rotor-only case.

The technique employs a straightforward engineering approach to arbitrary vortex flow axial fan rotor design, and the equation set can be solved by using relatively simple numerical methods. Estimates of pressure rise and shaft power characteristics for a proposed fan/rotor design can be computed and the design loop iterated until an acceptable set of blade parameters is identified. It is also possible to analyze the performance of an existing axial fan installation as a prelude to the design of a more efficient and effective replacement rotor.

Experimental data used in validating the design and analysis techniques are also presented. These data include comprehensive Cobra pressure probe surveys of local flow parameters downstream of three different low boss ratio, low solidity, arbitrary vortex flow rotors (all with circular arc camber line type blades) as well as fan performance characteristics for one of the experimental rotors configured as a direct-exhaust fan unit. Installation-dependent factors such as direct-exhaust losses and tip clearance effects are also examined. The analytical technique is shown to provide acceptable estimates of fan/rotor pressure rise performance and shaft power characteristics over a moderately wide range of blade angles and operating conditions.     

Static equilibrium of an elastic orthotropic medium with an arbitrarily oriented triaxial ellipsoidal inclusion

The stress state of an elastic orthotropic medium with an arbitrarily oriented triaxial ellipsoidal inclusion is analyzed. A solution is obtained using the triple Fourier transform and the Fourier-transformed Green’s function for an infinite anisotropic medium. The high efficiency of the approach is demonstrated by solving the problem for a transversely isotropic material with a spheroidal cavity for which the exact solution is known. A numerical analysis is conducted to study the stress distribution over the surface of the inclusion with different orientations in the orthotropic space. It is revealed that in some cases the orientation of the inclusion has a strong effect on the stress concentration __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 4, pp. 55–61, April 2007.     

Bending of Composite Plates under Static and Dynamic Loading

The stationary bending of a two-dimensional elastic system of joined rectangular plates with different mechanical properties and boundary conditions is studied. A technique for solving the corresponding problem is proposed. It is based on approximate approaches in combination with a generalization of the force method. It is established that the static strain state of the system is essentially dependent on the difference between the elastic moduli of the constituent plates. It is also shown that with certain boundary conditions, the dependence of the first resonant frequency on the relative position of plates along the contact line is nonmonotonic__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 1, pp. 77–84, January 2005.     

Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD

Based on the successive iteration in the Taylor series expansion method,a three-point explicit compact difference scheme with arbitrary order of accuracy is derived in this paper.Numerical characteristics of the scheme are studied by the Fourier analysis. Unlike the conventional compact difference schemes which need to solve the equation to obtain the unknown derivatives in each node,the proposed scheme is explicit and can achieve arbitrary order of accuracy in space.Application examples for the convection- diffusion problem with a sharp front gradient and the typical lid-driven cavity flow are given.It is found that the proposed compact scheme is not only simple to implement and economical to use,but also is effective to simulate the convection-dominated problem and obtain high-order accurate solution in coarse grid systems.