The Modified Ghost Fluid Method Applied to Fluid-Elastic Structure Interaction

In this work, the modified ghost fluid method is developed to deal with 2D compressible fluid interacting with elastic solid in an Euler-Lagrange coupled system. In applying the modified Ghost Fluid Method to treat the fluid-elastic solid coupling, the Navier equations for elastic solid are cast into a system similar to the Euler equations but in Lagrangian coordinates. Furthermore, to take into account the influence of material deformation and nonlinear wave interaction at the interface, an Euler-Lagrange Riemann problem is constructed and solved approximately along the normal direction of the interface to predict the interfacial status and then define the ghost fluid and ghost solid states. Numerical tests are presented to verify the resultant method.     

模拟多介质界面问题的虚拟流体方法综述

虚拟流体方法为模拟具有清晰物质界面的多介质流动问题提供了一种简便途径.尤其基于多介质Riemann问题解的修正虚拟流体方法及其变体，能够真实考虑到界面附近非线性波的相互作用和物质性质的影响，可以有效解决各种界面强间断等挑战性难题，具有巨大的工程应用潜力.文章重点回顾了虚拟流体方法的发展历史，总结和对比了各种代表性版本在模拟可压缩多介质流时的界面条件定义方式和多维推广方式，并介绍了该方法的设计原则和精度分析方面的研究进展.文章还回顾了该方法在其他更广泛和更具挑战性典型科学问题中的最新应用进展，并对方法的优势和特点进行了总结.      

Accuracies and conservation errors of various ghost fluid methods for multi-medium Riemann problem

Since the (original) ghost fluid method (OGFM) was proposed by Fedkiw et al. in 1999 [5], a series of other GFM-based methods such as the gas–water version GFM (GWGFM), the modified GFM (MGFM) and the real GFM (RGFM) have been developed subsequently. Systematic analysis, however, has yet to be carried out for the various GFMs on their accuracies and conservation errors. In this paper, we develop a technique to rigorously analyze the accuracies and conservation errors of these different GFMs when applied to the multi-medium Riemann problem with a general equation of state (EOS). By analyzing and comparing the interfacial state provided by each GFM to the exact one of the original multi-medium Riemann problem, we show that the accuracy of interfacial treatment can achieve “third-order accuracy” in the sense of comparing to the exact solution of the original mutli-medium Riemann problem for the MGFM and the RGFM, while it is of at most “first-order accuracy” for the OGFM and the GWGFM when the interface approach is actually near in balance. Similar conclusions are also obtained in association with the local conservation errors. A special test method is exploited to validate these theoretical conclusions from the numerical viewpoint.     

虚拟流体方法的设计原则

研究模拟可压缩多介质流的虚拟流体方法,建立一般状态方程下定义虚拟流体状态的基本原则.根据波系结构和使用的自由变量分别推导虚拟流体状态的定义方式.结果表明在这些方式下求解多介质Riemann问题理论上完全精确.进一步总结几种简单有效的虚拟流体方法,这些定义方式不依赖于虚拟流体区域可能产生的波系结构.其中一种类似于反射边界条件,只是界面速度需要首先精确预测出来.数值算例验证了研究结果的合理性.     

Scattering of sound by an elastic plate with flow

An elastic plate, set in an infinite baffle and immersed in a fluid moving with a uniform subsonic velocity, is excited by an acoustic source. The scattered sound field is analyzed when fluid-plate coupling is large, and a solution is found by the use of matched asymptotic expansions. The far field is found to approximate to the solution obtained when the elastic plate is absent. At a plate resonance, however, the outer field must include eigensolutions with singularities at the plate edges, and close to the plate the dominant terms are travelling plate waves. These plate waves are found to have a wavelength independent of the frequency of the source. It is also shown that a plate resonance corresponds to a divergence instability of aerodynamic flutter theory and that the stability results found in this paper are in agreement with those obtained by using modal expansions. The limit as the Mach number goes to zero is found to be singular, suggesting an analysis of the model for small flow velocity. This calculation is performed and the results match smoothly to the respective solutions for a stationary fluid and for a large subsonic flow.     

Numerical simulation of underwater explosion bubble with a refined interface treatment

With the intermediate flow states predicted by local two phase Riemann problem,the modified ghost fluid method(MGFM)and its variant(r GFM)have been widely employed to resolve the interface condition in the simulation of compressible multi-medium flows.In this work,the drawback of the construction procedure of local two phase Riemann problem in r GFM was investigated in detail,and a refined version of the construction procedure was specially developed to make the simulation of underwater explosion bubbles more accurate and robust.Beside the refined r GFM,the fast and accurate particle level set method was also adopted to achieve a more effective and computationally efficient capture of the evolving multi-medium interfaces during the simulation.To demonstrate the improvement brought by current refinement,several typical numerical examples of underwater explosion bubbles were performed with original r GFM and refined r GFM,respectively.The results indicate that,when compared with original r GFM,numerical oscillations were effectively removed with the proposed refinement.Accordingly,with present refined treatment of interface condition,a more accurate and robust simulation of underwater explosion bubbles was accomplished in this work.     

运动激波和气泡串相互作用的初步数值模拟

通过对激波和流体界面相互作用诱导的大变形界面演化的数值模拟,验证Level set方法精确模拟多个流体界面的有效性.采用2阶迎风TVD求解欧拉方程得到流场解,采用5阶WENO求解Level set方程追踪多流体界面,采用GFM方法处理流体内界面.利用文[1]的计算结果校核本文程序.在此基础上,对运动激波和气泡串相互作用过程进行了初步数值模拟,得到了不同时刻运动激波和圆管内的两个气泡作用后的演化图象,包括压力和密度等值线分布.计算结果表明:针对推广后的多界面Level set方程,该方法仍可高质量地捕捉多个流体界面.     

Curvilinear Immersed Boundary Method for Simulating Fluid Structure Interaction with Complex 3D Rigid Bodies

The sharp-interface CURVIB approach of Ge and Sotiropoulos [L. Ge, F. Sotiropoulos, A Numerical Method for Solving the 3D Unsteady Incompressible Navier-Stokes Equations in Curvilinear Domains with Complex Immersed Boundaries, Journal of Computational Physics 225 (2007) 1782-1809] is extended to simulate fluid structure interaction (FSI) problems involving complex 3D rigid bodies undergoing large structural displacements. The FSI solver adopts the partitioned FSI solution approach and both loose and strong coupling strategies are implemented. The interfaces between immersed bodies and the fluid are discretized with a Lagrangian grid and tracked with an explicit front-tracking approach. An efficient ray-tracing algorithm is developed to quickly identify the relationship between the background grid and the moving bodies. Numerical experiments are carried out for two FSI problems: vortex induced vibration of elastically mounted cylinders and flow through a bileaflet mechanical heart valve at physiologic conditions. For both cases the computed results are in excellent agreement with benchmark simulations and experimental measurements. The numerical experiments suggest that both the properties of the structure (mass, geometry) and the local flow conditions can play an important role in determining the stability of the FSI algorithm. Under certain conditions unconditionally unstable iteration schemes result even when strong coupling FSI is employed. For such cases, however, combining the strong-coupling iteration with under-relaxation in conjunction with the Aitken's acceleration technique is shown to effectively resolve the stability problems. A theoretical analysis is presented to explain the findings of the numerical experiments. It is shown that the ratio of the added mass to the mass of the structure as well as the sign of the local time rate of change of the force or moment imparted on the structure by the fluid determine the stability and convergence of the FSI algorithm. The stabilizing role of under-relaxation is also clarified and an upper bound of the required for stability under-relaxation coefficient is derived.     

Smoothed particle hydrodynamics(SPH) for modeling fluid-structure interactions

Fluid-structure interaction(FSI) is a class of mechanics-related problems with mutual dependence between the fluid and structure parts and it is observable nearly everywhere, in natural phenomena to many engineering systems. The primary challenges in developing numerical models with conventional grid-based methods are the inherent nonlinearity and timedependent nature of FSI, together with possible large deformations and moving interfaces. Smoothed particle hydrodynamics(SPH) method is a truly Lagrangian and meshfree particle method that conveniently treats large deformations and naturally captures rapidly moving interfaces and free surfaces. Since its invention, the SPH method has been widely applied to study different problems in engineering and sciences, including FSI problems. This article presents a review of the recent developments in SPH based modeling techniques for solving FSI-related problems. The basic concepts of SPH along with conventional and higher order particle approximation schemes are first introduced. Then, the implementation of FSI in a pure SPH framework and the hybrid approaches of SPH with other grid-based or particle-based methods are discussed. The SPH models of FSI problems with rigid, elastic and flexible structures, with granular materials, and with extremely intensive loadings are demonstrated. Some discussions on several key techniques in SPH including the balance of accuracy, stability and efficiency, the treatment of material interface, the coupling of SPH with other methods, and the particle regularization and adaptive particle resolution are provided as concluding marks.     

Viscoelastic theory for nematic interfaces

A complete macroscopic theory for compressible nematic-viscous fluid interfaces is developed and used to characterize the interfacial elastic, viscous, and viscoelastic material properties. The derived expression for the interfacial stress tensor includes elastic and viscous components. Surface gradients of the interfacial elastic stress tensor generates tangential Marangoni forces as well as normal forces. The latter may be present even in planar surfaces, implying that in principle static planar interfaces may accommodate pressure jumps. The asymmetric interfacial viscous stress tensor takes into account the surface nematic ordering and is given in terms of the interfacial rate of deformation and interfacial Jaumann derivative. The material function that describes the anisotropic viscoelasticity is the dynamic interfacial tension, which includes the interfacial tension and dilational viscosities. Viscous dissipation due to interfacial compressibility is described by the anisotropic dilational viscosity, and it is shown to describe the Boussinesq surface fluid appropriate for Newtonian interfaces when the director is homeotropic. Three characteristic interfacial shear viscosities are defined according to whether the surface orientation is along the velocity direction, the velocity gradient, or the unit normal. In the last case the expression reduces to the interfacial shear viscosity of the Boussinesq surface fluid. The theory provides a theoretical framework to study interfacial stability, thin liquid film stability and hydrodynamics, and any other interfacial rheology phenomena.     

三角形网格上多介质流动问题界面处理方法

多介质流动问题的求解一般是在结构网格上实现,而三角形网格对于复杂计算区域具有更好的适应性,本文结合rGFM方法,给出三角形网格上多介质流动问题界面处理方法.利用level-set方法跟踪界面,在界面处构造Riemann问题,得到界面处流体准确的流动状态.通过定义界面边界条件,将多介质流动问题转化为单介质流动问题,利用高精度RKDG方法求解.采用多个算例验证该方法的稳健性和有效性,结果表明该方法能准确捕捉界面和激波的位置,保持界面清晰.     

基于求解Riemann问题的界面处理方法

通过在界面处构造Riemann问题,根据流体的法向速度和压力在界面(接触间断)处连续的特性,利用Riemann问题的解不仅定义了ghost流体的值,而且对真实流体中邻近界面的点值进行了更新,使得在界面处的流体的状态满足接触间断的性质,给出了更加精确的界面边界条件,守恒误差分析表明该方法在界面计算过程中引入较小的误差.数值试验表明该方法能准确地捕捉界面和激波的位置.     

多介质流动数值计算中的界面处理方法

 求解Riemann问题得到界面接触间断的流动状态,并以此构造带状区域的虚拟流体状态,对于多维问题设计了一种方便有效的算法。同时求解耦合的守恒形式欧拉方程组和非守恒界面捕捉方程,并用Level-Set函数捕捉界面,数值计算采用高分辨率MWENO格式。最后对可压缩多介质流动问题进行了数值模拟。     

弹性板水下爆炸冲击载荷的修正虚拟流体方法分析

研究非接触水下爆炸中板结构厚度和密度对冲击载荷和空化演变的影响.将近期发展的修正的虚拟流体方法推广应用于处理可压缩流体与弹性板结构的非线性相互作用.研究发现,假设结构不发生断裂破坏,较小的板厚度或密度可以减弱爆炸波对结构的冲击,冲击载荷随着厚度或密度的增大非线性地趋向于固壁边界下的结果.而且较小的板厚度或密度可以促使空化较早形成,产生较大的空化区,并推迟和减弱空化破裂对结构的二次冲击.     

Vibrations of rectangular plates reinforced by any number of beams of arbitrary lengths and placement angles

This paper presents an analytical method for the vibration analysis of plates reinforced by any number of beams of arbitrary lengths and placement angles. Both the plate and stiffening beams are generally modeled as three-dimensional (3-D) structures having six displacement components at a point, and the coupling at an interface is generically described by a set of distributed elastic springs. Each of the displacement functions is here invariably expressed as a modified Fourier series, which consists of a standard Fourier cosine series plus several supplementary series/functions used to ensure and improve uniform convergence of the series representation. Unlike most existing techniques, the current method offers a unified solution to the vibration problems for a wide spectrum of stiffened plates, regardless of their boundary conditions, coupling conditions, and reinforcement configurations. Several numerical examples are presented to validate the methodology and demonstrate the effect on modal parameters for a stiffened plate with various boundary conditions, coupling conditions, and reinforcement configurations.     

基于数值仿真船体结构毁伤测点布设方法研究

为有效测量水下爆炸作用下船体结构的毁伤程度,优化测点布设,建立了典型舰船有限元模型,运用大型有限元分析软件ABAQUS提供的声固耦合方法,对船体结构在爆炸载荷作用下的冲击响应进行仿真。依据数值仿真结果,分析船体结构薄弱环节,并将这些薄弱环节作为重点测量部位,结合工程实践和海上试验经验,在准确测量结构毁伤的基础上,合理布设传感器测点,研究设计了水下爆炸作用下船体结构毁伤测点布设方法。经海上爆炸试验验证,数值仿真结果与试验测试结果吻合良好,该毁伤测点布设方法设计合理。     

SPH方法气-液界面边界条件研究

针对界面附近粒子光滑函数截断和非物理穿透问题,提出一种气-液界面边界条件的处理方法.当界面附近支持域出现不同材料粒子,每步计算可在支持域设置虚粒子,按照密度分配方法给虚粒子物理量赋值,并对界面附近粒子引入气-液两相阻力.采用SPH方法和Level-Set方法,计算运动激波对气-液界面作用问题,两者计算结果一致,初步验证了气-液界面边界条件处理的适用性.用SPH方法分别计算超声速气流中的圆截面液柱绕流和下落问题,界面两侧粒子压力和法向速度连续,给出弓形激波、回流区和下游回流区等定性合理结果.表明本文方法可适度避免界面附近流体粒子光滑截断和粒子非物理穿透现象、界面附近流场数值振荡.     

Deforming composite grids for solving fluid structure problems

We describe a mixed Eulerian–Lagrangian approach for solving fluid–structure interaction (FSI) problems. The technique, which uses deforming composite grids (DCG), is applied to FSI problems that couple high speed compressible flow with elastic solids. The fluid and solid domains are discretized with composite overlapping grids. Curvilinear grids are aligned with each interface and these grids deform as the interface evolves. The majority of grid points in the fluid domain generally belong to background Cartesian grids which do not move during a simulation. The FSI-DCG approach allows large displacements of the interfaces while retaining high quality grids. Efficiency is obtained through the use of structured grids and Cartesian grids. The governing equations in the fluid and solid domains are evolved in a partitioned approach. We solve the compressible Euler equations in the fluid domains using a high-order Godunov finite-volume scheme. We solve the linear elastodynamic equations in the solid domains using a second-order upwind scheme. We develop interface approximations based on the solution of a fluid–solid Riemann problem that results in a stable scheme even for the difficult case of light solids coupled to heavy fluids. The FSI-DCG approach is verified for three problems with known solutions, an elastic-piston problem, the superseismic shock problem and a deforming diffuser. In addition, a self convergence study is performed for an elastic shock hitting a fluid filled cavity. The overall FSI-DCG scheme is shown to be second-order accurate in the max-norm for smooth solutions, and robust and stable for problems with discontinuous solutions for a wide range of constitutive parameters.     

A matrix-free,implicit, incompressible fractional-step algorithm for fluid–structure interaction applications

In this paper we detail a fast, fully-coupled, partitioned fluid–structure interaction (FSI) scheme. For the incompressible fluid, new fractional-step algorithms are proposed which make possible the fully implicit, but matrix-free, parallel solution of the entire coupled fluid–solid system. These algorithms include artificial compressibility pressure-poisson solution in conjunction with upwind velocity stabilisation, as well as simplified pressure stabilisation for improved computational efficiency. A dual-timestepping approach is proposed where a Jacobi method is employed for the momentum equations while the pressures are concurrently solved via a matrix-free preconditioned GMRES methodology. This enables efficient sub-iteration level coupling between the fluid and solid domains. Parallelisation is effected for distributed-memory systems. The accuracy and efficiency of the developed technology is evaluated by application to benchmark problems from the literature. The new schemes are shown to be efficient and robust, with the developed preconditioned GMRES solver furnishing speed-ups ranging between 50 and 80.     

Stress-Assisted Reaction at a Solid-Fluid Interface

On the interface between a solid and a fluid, a reaction can occur in which atoms either leave the solid to join the fluid, or leave the fluid to join the solid. If the solid is in addition subject to a mechanical load, two outcomes may be expected. The reaction may proceed uniformly, so that the interface remains flat as the solid recedes or extends. Alternatively, the reaction may cause the interface to roughen and develop sharp cracks, leading to fracture. This paper reviews the current understanding of the subject. The solid-fluid is a thermodynamic system: the solid is in elastic equilibrium with the mechanical load, but not in chemical equilibrium with the fluid. Thermodynamic forces that drive the interfacial reaction include chemical energy difference between the solid and the fluid, elastic energy stored in the solid, and interfacial energy. The reaction is taken to be thermally activated. A kinetic law is adopted in which the stress affects both the activation energy and the driving force of the interface reaction. A linear perturbation analysis identifies the stability condition, which differs substantially from the well known stability condition based on the driving force alone. Large perturbations are examined by assuming that the interface varies as a family of cycloids, from slight waviness to sharp cracks. An analytic elasticity solution is used to compute the stress field in the solid, and a variational method to evolve the shape of the interface.