微重力环境下Marangoni对流的有限元数值模拟

本文以悬浮区为背景研究液桥中气液交界面上由表面张力所驱动的流体对流。我们采用有限元方法对轴对称、定常运动方程,能量方程,自由面切向、法向应力平衡条件迭代求解,首次考虑了边界形状的影响,得到了自洽的流区位形和流场、温度场、表面压强分布。一般说来,流区自由面呈弯月形。结果表明,只要表面张力数 C_a<1,静态平衡的流区位形就是动态情况的良好近似。本文还分析了C_a数及G_r 数对流区位形的影响,得到了不同 M_a 数及散热条件下的温度场和流函数分布。     

用有限元自适网格控制结构应力分析的精度

本文涉及有限元结构应务分析精度控制的三个基本方面：１）自适应的有限元网格自动生成技术；２）有限元应力分析的误差估计：３）用单元尺寸场将误差分布用于形成新网格。     

物理参数变化对短脉冲激光激励温度场的影响

 为研究多物理参数（耦合系数、电子热导率、电子热容、晶格热容）同时随温度变化对短脉冲激光辐照金属材料产生温度场分布的影响,基于双温耦合理论,建立了短脉冲激光辐照金属材料金的加热过程的有限元求解模型。在同时考虑脉冲激光的空间、时间分布和多参数同时随温度变化的情况下,得到短脉冲激光辐照金属材料金激励产生的温度场二维瞬态分布,并进一步比较了多物理参数同时随温度变化和采用室温物理参数两种情况下温度场分布的区别。数值结果表明：多物理参数同时随温度变化使电子温度和晶格温度的上升变快,最大值变大,而且使得材料中激光穿透直接辐照到的区域温度变高。     

Computational performance of Free Mesh Method applied to continuum mechanics problems

The free mesh method (FMM) is a kind of the meshless methods intended for particle-like finite element analysis of problems that are difficult to handle using global mesh generation, or a node-based finite element method that employs a local mesh generation technique and a node-by-node algorithm. The aim of the present paper is to review some unique numerical solutions of fluid and solid mechanics by employing FMM as well as the Enriched Free Mesh Method (EFMM), which is a new version of FMM, including compressible flow and sounding mechanism in air-reed instruments as applications to fluid mechanics, and automatic remeshing for slow crack growth, dynamic behavior of solid as well as large-scale Eigen-frequency of engine block as applications to solid mechanics.     

A symmetric BEM–FEM method for an axisymmetric eddy current problem

This paper deals with the numerical solution of the time-harmonic eddy current model in an axisymmetric unbounded domain. To this end, a new symmetric BEM–FEM formulation is derived and also analyzed. Moreover, error estimates for the corresponding discretization are proven.     

Simulation of oscillatory flow in an aortic bifurcation using FVM and FEM: A comparative study of implementation strategies

Cardiovascular diseases are one of the major causes of long‐term morbidity and mortality in human beings. The nearly epidemic increase in prevalence of such diseases poses a serious threat to public health and calls for efficient methods of diagnosis and treatment. Non‐invasive diagnostic procedures such as MRI are often used in this context; however, these are limited in terms of spatial and temporal resolution and do not provide information on time‐dependent pressures and wall shear stresses—key quantities considered to be partially responsible for the formation and development of related pathologies. The present study is concerned with the numerical simulation of oscillatory flow through the abdominal aortic bifurcation. Computational fluid dynamics simulation of oscillatory flow in a branched geometry at high Reynolds numbers poses considerable challenges. The present study reports a detailed comparison of simulations performed with a finite volume and a finite element method, two approaches with significant differences in their discretization strategy, treatment of boundary conditions and other numerical aspects. Both solvers were parallelized, using loop parallelization of the BiCGStab linear solver for the finite volume and domain decomposition based on the Schur complement method for the finite element technique. The experience gained with these two approaches for the solution of flow in a bifurcation forms the focus of this study. Although similar results were obtained for both methods, the computation time required for convergence was found to be significantly smaller for the finite element approach. Copyright © 2010 John Wiley & Sons, Ltd.     

A least‐square weighted residual method for level set formulation

In this paper, a least‐square weighted residual method (LSWRM) for level set (LS) formulation is introduced to achieve interface capturing in two‐dimensional (2D) and three‐dimensional (3D) problems. An LSWRM was adopted for two semi‐discretized advection and reinitialization equations of the LS formulation. The present LSWRM provided good mathematical properties such as natural numerical diffusion and the symmetry of the resulting algebraic systems for the advection and reinitialization equations. The proposed method was validated by solving some 2D and 3D benchmark problems such as those involving a rotating slotted disk, the rotation of a slotted sphere, and a time‐reversed single‐vortex flow and a deformation problem of a spherical fluid. The numerical results were compared with those obtained from essentially non‐oscillatory type formulations and particle LS methods. Further, the proposed LSWRM for the LS formulation was coupled with a splitting finite element method code to solve the incompressible Navier–Stokes equations, and then, the collapse of a 3D broken dam flow was well simulated; in the simulation, the entrapping of air and the splashing of the surge front of water were reproduced. The mass conservation of the present method was found to be satisfactory during the entire simulation. Copyright © 2011 John Wiley & Sons, Ltd.     

EFFECT OF HEAD SHAPE ON THE ADHESION CAPABILITY OF MUSHROOM-LIKE BIOLOGICAL ADHESIVE STRUCTURES

Theoretical studies have been extensively performed to reveal the mechanism of the biological adhesive structure with a flat head.However,most biological adhesive structures have a big head.Here we investigate the effect of the head shape on the adhesion capability of a mushroomlike adhesion structure based on a finite element model.We show that a mushroom-like head produces significant adhesion enhancement.We investigate also the effect of the elastic modulus of the adhesive structure on the adhesion.Our results are consistent with existing experimental results.     

Convergence on layer-adapted meshes and anisotropic interpolation error estimates of non-standard higher order finite elements

For a general class of finite element spaces based on local polynomial spaces E with Pp⊂E⊂Qp we construct a vertex-edge-cell and point-value oriented interpolation operators that fulfil anisotropic interpolation error estimates.Using these estimates we prove ε-uniform convergence of order p for the Galerkin FEM and the LPSFEM for a singularly perturbed convection-diffusion problem with characteristic boundary layers.