An axisymmetric hypersingular boundary integral formulation for simulating acoustic wave propagation in supercavitating flows

Flow supercavitation begins when fluid is accelerated over a sharp edge, usually at the nose of an underwater vehicle, where a phase change occurs and causes a low density gaseous cavity to gradually envelop the whole object (supercavity) thereby allowing for higher speeds of underwater vehicles. The supercavity may be maintained through ventilated cavitation caused by injection of gases into the cavity, which causes fluctuations at the vapor–water interface. A major issue that concerns the efficient operation of an underwater object’s guidance system (which is achieved by high frequency acoustic sensors mounted within the nose region), is the hydrodynamic noise produced due to the fluctuating vapor–water interface. It is important to carry out a detailed study on the effect of self-noise at the vehicle’s nose that is generated by the ventilating gas jet impingement on the supercavity wall. For this purpose, the present study uses a boundary element method which is more versatile compared to other numerical techniques such as the finite element/finite difference methods. The variation of acoustic pressure at the vehicle nose for various shapes of cavitators, boundary conditions and jet impact diameters are presented. Comparisons are made with the semi-analytical procedure of Howe et al. (Howe et al., On self-noise at the nose of a supercavitating vehicle. Journal of Sound and Vibration, 322 (2009a), 772–784) and finite element based COMSOL commercial package. Several issues pertaining to the behaviour of analytical and numerical results are highlighted. Finally, the proposed boundary element technique is used to study arbitrary shapes of supercavities which may encountered at various stages of supercavity development.     

Finite element analysis of nano-scale Timoshenko beams using the integral model of nonlocal elasticity

Stress-strain relation in Eringen's nonlocal elasticity theory was originally formulated within the framework of an integral model. Due to difficulty of working with that integral model, the differential model of nonlocal constitutive equation is widely used for nanostructures. However, paradoxical results may be obtained by the differential model for some boundary and loading conditions. Presented in this article is a finite element analysis of Timoshenko nano-beams based on the integral model of nonlocal continuum theory without employing any simplification in the model. The entire procedure of deriving equations of motion is carried out in the matrix form of representation, and hence, they can be easily used in the finite element analysis. For comparison purpose, the differential counterparts of equations are also derived. To study the outcome of analysis based on the integral and differential models, some case studies are presented in which the influences of boundary conditions, nonlocal length scale parameter and loading factor are analyzed. It is concluded that, in contrast to the differential model, there is no paradox in the numerical results of developed integral model of nonlocal continuum theory for different situations of problem characteristics. So, resolving the mentioned paradoxes by means of a purely numerical approach based on the original integral form of nonlocal elasticity theory is the major contribution of present study.     

Modelling of spall damage in ductile materials and its application to the simulation of plate impact on copper

A statistical model of dynamic spall damage due to void nucleation and growth is proposed for ductile materials under intense loading, which takes into account inertia, elastic-plastic effect, and initial void size. To some extent, void interaction could be accounted for in this approach. Based on this model, the simulation of spall experiments for copper is performed with the Lagrangian finite element method. The simulation results are in good agreement with experimental data for the free surface velocity profile, stress record behind copper target, final porosity, and void concentrations across the target. The influence of elastic-plastic effect upon the damage evolution is explored. The correlation between the damage evolution and the history of the stress near the spall plane is also analyzed.     

Modelling of spall damage in ductile materials and its application to the simulation of the plate impact on copper

A statistical model of dynamic spall damage due to void nucleation and growth is proposed for ductile materials under intense loading, which takes into account inertia, the elastic-plastic effect, and initial void size. To some extent, void interaction could be accounted for in this approach. Based on this model, the simulation of spall experiments for copper is performed by using the Lagrangian finite element method. The simulation results are in good agreement with experimental data for the free surface velocity profile, stress record behind copper target, final porosity, and void concentrations across the target. The influence of elastic-plastic effect upon the damage evolution is explored. The correlation between the damage evolution and the history of the stress near the spall plane is also analyzed.     

The surface cracking behavior in air plasma sprayed thermal barrier coating system incorporating interface roughness effect

The objective of this work is to understand the effect of interface roughness on the strain energy release rate and surface cracking behavior in air plasma sprayed thermal barrier coating system. This is achieved by a parameter investigation of the interfacial shapes, in which the extended finite element method (XFEM) and periodic boundary condition are used. Predictions for the stress field and driving force of multiple surface cracks in the film/substrate system are presented. It is seen that the interface roughness has significant effects on the strain energy release rate, the interfacial stress distribution, and the crack propagation patterns. One can see the completely different distributions of stress and strain energy release rate in the regions of convex and concave asperities of the substrate. Variation of the interface asperity is responsible for the oscillatory characteristics of strain energy release rate, which can cause the local arrest of surface cracks. It is concluded that artificially created rough interface can enhance the durability of film/substrate system with multiple cracks.     

多晶材料晶粒生长粗化过程的相场方法模拟

基于采用晶体有序化程度参量ψ和晶体学取向θ来表示多晶粒结构的相场模型,利用自适应有限元方法模拟了多晶材料等温过程中的晶粒粗化现象.模拟结果显示,在曲率作用下,通过晶界迁移弯曲晶界逐渐平直化,小晶粒逐渐被大晶粒吞并,当晶界之间的取向差较小时,满足一定能量和几何条件的两晶粒在界面能作用下会发生转动,合并为单个晶粒.模拟结果与实验结果符合较好.因此,该相场模型可以很好地用来模拟固态相变中多晶材料的生长粗化等现象. 关键词： 相场 晶界迁移 晶粒转动 粗化     

石墨烯-纳米探针相互作用有限元准静态计算

纳米尺度探针是研究纳米薄膜材料的重要工具.针对纳米探针和石墨烯相互作用有限元模型静态计算中难以收敛的困难,应用动态显式算法通过间歇式探针进给方式进行能量耗散,得出静态计算结果.模型中界面作用力由界面黏结能和原子间作用势导出并植入Abaqus软件中界面作用子程序,实现对石墨烯、探针,基体系统内相互作用的仿真计算.通过对比计算结果和实验数据,对实验结果给出了一致性解释.     

三维粗糙面电磁双站散射的直接型区域分解计算

提出三维粗糙面双站电磁散射的直接型有限元-区域分解方法.首先建立含有迭代Robin边界条件(IRBC)的区域分解法耦合模型,再用内视法导出高度稀疏分块的分区耦合矩阵,之后给出缩减耦合矩阵带宽的子区域排序方法和IRBC的FFT加速算法.用有限元-完全匹配层和未分区的有限元-IRBC方法验证数值结果.     

A Method of Lines Based on Immersed Finite Elements for Parabolic Moving Interface Problems

This article extends the finite element method of lines to a parabolic initial boundary value problem whose diffusion coefficient is discontinuous across an interface that changes with respect to time. The method presented here uses immersed finite element (IFE) functions for the discretization in spatial variables that can be carried out over a fixed mesh (such as a Cartesian mesh if desired), and this feature makes it possible to reduce the parabolic equation to a system of ordinary differential equations (ODE) through the usual semi-discretization procedure. Therefore, with a suitable choice of the ODE solver, this method can reliably and efficiently solve a parabolic moving interface problem over a fixed structured (Cartesian) mesh. Numerical examples are presented to demonstrate features of this new method.     

超声脉冲窄带激励的水泥环第二界面探测方法*

如何评价水泥环的第二界面固井质量一直是声波测井领域的世界性难题。本文利用有限元方法研究了基于超声脉冲反射法的水泥环第二界面固井质量探测方法。数值模拟了不同声源脉冲宽度及不同声阻抗地层对第二界面超声反射回波的影响。研究结果表明采用窄带的激励信号源可有效的提高水泥环第二界面反射波回波幅度,有望解决水泥环第二界面固井质量评价问题。制作了水泥环第二界面具有胶结缺陷的实验室样品及大型固井质量刻度井,在此基础上开展了基于窄带超声脉冲激励的实验室及刻度井群实验研究。实验结果表明满足水泥环第一界固井质量良好条件,采用窄带的超声脉冲回波技术可有效评价水泥环第二界面固井质量。     

Numerical investigation of the effect of insoluble surfactant on drop formation in microfluidic device

Surfactants are commonly used in droplet-based microfluidics to stabilize the droplet interface. In this study, we investigate the effect of insoluble surfactant on drop formation in a capillary microfluidic device. We use a diffuse-interface method to describe the evolution of interface involving insoluble surfactant. The Navier-Stokes/Cahn-Hilliard equations and the surfactant conservation equation are solved by a finite element method along with a grid deformation method. As the surfactant has a non-uniform distribution during the drop formation in general, the surface tension has a gradient on the interface, which affects the flow field and interface evolution. The surfactant effect is discussed for dripping and jetting regimes.     

Stress influence on high temperature oxide scale growth: Modeling and investigation on a thermal barrier coating system

In thermal barrier coating (TBC) systems, an oxide layer develops at high temperature below the ceramic coating, leading in the long term to the mechanical failure of the structure upon cooling. The mechanism of stress-affected oxidation likely to induce the growth of a non-uniform oxide scale detrimental to the TBC lifetime was investigated. A continuum thermodynamics formulation is derived accounting for the influence of the stress and strain situation at the sharp metal/oxide phase boundary on the local oxidation kinetics. It specially includes the contributions of the large volumetric strain and the mass consumption associated with metal oxidation. A continuum mechanics/mass diffusion framework is used along with the developed formulation for the interface evolution to study the growth of an oxide layer coupled with local stress development. The implementation of the model has required the development of a specific simulation tool, based on a finite element method completed with an external routine for the phase boundary propagation. Results on an electron-beam physical vapor deposited (EB-PVD) TBC case are presented. The processes resulting in a non-uniform oxide scale growth are analyzed and the main influences are discussed.     

强脉冲离子束辐照薄金属靶的热力学过程研究

在回顾和总结现有强脉冲离子束诊断技术和能量沉积模型的基础上, 结合红外成像诊断分析, 基于能量平衡, 提出了强脉冲离子束在固体靶中功率密度分布模型, 并采用蒙特卡罗方法对其进行计算. 以该功率密度模型作为源项, 使用有限元分析方法模拟强脉冲离子束入射100 μm不锈钢靶后内部温度场在毫秒时间范围内的分布和演化. 结果显示, 在微秒时间范围内, 热场以存在于近表面区域数倍于离子射程范围内的冲击热场为主要特征; 而在毫秒时间范围内, 靶的前后表面(纵向)已达到温度平衡, 且靶后表面温度场和入射前表面的离子束横截面能量密度具有空间分布的相似性. 这证明了, 在采用具有毫秒响应速度的红外拍摄系统的情况下, 背面红外诊断技术可以实现以较高的精度对强脉冲离子束横截面的能量分布进行诊断和分析.     

A coupled solid–fluid method for modelling subduction

We present a novel dynamic approach for solid–fluid coupling by joining two different numerical methods: the boundary-element method (BEM) and the finite element method (FEM). The FEM results describe the thermomechanical evolution of the solid while the fluid is solved with the BEM. The bidirectional feedback between the two domains evolves along a Lagrangian interface where the FEM domain is embedded inside the BEM domain. The feedback between the two codes is based on the calculation of a specific drag tensor for each boundary on finite element. The approach is presented here to solve the complex problem of the descent of a cold subducting oceanic plate into a hot fluid-like mantle. The coupling technique is shown to maintain the proper energy dissipation caused by the important secondary induced mantle flow induced by the lateral migrating of the subducting plate. We show how the method can be successfully applied for modelling the feedback between deformation of the oceanic plate and the induced mantle flow. We find that the mantle flow drag is singular at the edge of the retreating plate causing a distinct hook shape. In nature, such hooks can be observed at the northern end of the Tonga trench and at the southern perimeter, of the South American trench.     

基于孔洞体胞模型铸造镁合金ML308的本构方程

铸造镁合金不可避免地包含许多微孔洞,这些微孔洞在材料的后续加工及服役过程中将发生演化,并对材料的力学行为产生重要影响.基于球形孔洞体胞模型,提出微孔洞长大及形核方程,它们构成微孔洞的演化方程.根据孔洞演化将造成材料性质弱化的物理机制,将微孔洞演化以弱化函数的形式引入到非经典弹塑性本构方程,得到考虑孔洞演化的铸造镁合金弹塑性本构方程.发展与本构方程相应的有限元数值分析程序,用其模拟了铸造镁合金ML308的微孔洞演化及力学行为,计算结果与实验结果符合较好. 关键词： 铸造镁合金 孔洞体胞模型 孔洞演化方程 本构方程     

Coarsening and self-organization in dilute diblock copolymer melts and mixtures

This paper explores the evolution of a sharp interface model for phase separation of copolymers in the limit of low volume fraction. Particles both exchange material as in usual Ostwald ripening, and migrate because of an effectively repulsive nonlocal energetic term. Coarsening via mass diffusion only occurs while particle radii are small, and they eventually approach a finite equilibrium size. Migration, on the other hand, is responsible for producing self-organized patterns.We construct approximations based upon an ansatz of spherical particles similar to the classical LSW theory to derive finite dimensional dynamics for particle positions and radii. For large systems, kinetic-type equations which describe the evolution of a probability density are constructed. For systems larger than the screening length, we obtain an analog of the homogenization result of Niethammer & Otto [B. Niethammer, F. Otto, Ostwald ripening: The screening length revisited, Calc. Var. Partial Differential Equations 13-1 (2001) 33-68]. A separation of timescales between particle growth and migration allows for a variational characterization of spatially inhomogeneous quasi-equilibrium states.     

3维Mortar谱元法模拟Kelvin-Helmhtz不稳定性混合层

采用Mortar谱元法和多处理器并行计算技术模拟了Kelvin-Helmhtz界面不稳定性湍流的混合发展过程,通过对混合层动量厚度、能谱和总动能的计算,评估了Kelvin-Helmhtz混合层的演化机理。计算结果表明:3维Mortar谱元法具有高计算精度和光滑区域的指数收敛特性,可以有效模拟混合层流动的湍流混合和演化,能够捕捉到涡的合并现象和大涡到小涡的级联过程;初期的混合层层流运动发展成具有连续谱结构的湍流运动过程,实现了Kelvin-Helmhtz界面不稳定性混合层流动从2维发展到3维的转捩特征,总湍流统计动能的变化反映了粘性耗散过程的作用。通过对Kelvin-Helmhtz 3维界面不稳定性混合层流动和3维层流向湍流转捩过程的数值模拟,程序的有效性得到了验算,表明谱元法应用于湍流混合模拟是可行的。     

Finite element prediction of debonding onset in short fiber reinforced polymers incorporating a hybrid interphase region

A robust finite element procedure for investigating damage evolution in short fiber reinforced polymeric composites under external loads is developed. This procedure is based on an axisymmetric unit cell composed of a fiber, surrounding interphase and bulk matrix. The hybrid interphase concept involves a degraded material phase, the extent of which is material and property dependent. One of the most significant features of the model relies on establishment of variable adhesion conditions between the primary material phases. The unit cell is discretized into linearly elastic elements for the fiber and the matrix and interface elements which allow debonding in the fiber–matrix interface. The interface elements fail according to critical stress and critical energy release rate criteria. The tension and shear aspects of failure are uncoupled, although the resulting nonlinear problem is solved implicitly utilizing quasi-static incremental loading conditions. Final failure resulting from saturation and breakage is modeled by the vanishing interface element technique. Details of the propagation of interface cracks and the initiation of debonds are also observed and discussed for various shapes of fiber end. Numerical results reveal an intense effect of the fiber-end geometry on the initial fiber–matrix de-cohesion. The present finite element procedures can generate meaningful results in the analysis of fiber-reinforced composites.