航行体水下发射流固耦合效应分析

对于水下发射过程来说,掌握水动力载荷形成机理与结构响应特征是一个亟待解决的问题.研究该问题需要考虑含相变的复杂多相流动,变约束的结构运动以及这二者之间的耦合效应.本文采用松耦合的方法,以流体求解器为主体,将自编的固体结构程序接入流体求解器中,在每个时间步长内分别对流体动力学方程和固体结构动力学方程进行求解,通过流固界面之间的数据交换实现耦合计算.其中,流体求解器基于雷诺平均纳维斯托克斯方程,采用单流体模型处理多相流问题,引入空化模型描述空化相变,采用修正的湍流模型模拟混合物的湍流效应,并采用动网格技术处理移动边界问题.航行体的刚体运动和结构振动分开求解.结构求解器采用等效梁模型描述结构的振动,通过坐标变换给出了随体坐标系下的结构振动方程,求解方法采用时域积分法.所建立的流固耦合方法不仅能够捕捉到自然空化的演化情况,还可获得航行体所受水动力、结构振动响应以及截面的弯矩,获得了实验的验证.基于该方法研究了结构刚度、发射速度对空泡溃灭与结构振动耦合效应的影响规律.结果表明,同步溃灭是影响结构载荷的主要因素,包括溃灭压力幅值,溃灭压力作用位置,以及溃灭压力与结构振动的相位关系.     

Plectoneme formation in twisted fluctuating rods

The mechanics of DNA supercoiling is a subject of crucial importance to uncover the mechanism and kinetics of several enzymes. It is therefore being investigated using several biochemical and biophysical methods including single molecule experimental techniques. An interesting problem within this realm is that of torsional buckling and plectoneme formation in DNA as it is simultaneously put under tensile and torsional stress. Analytical solutions to this problem are difficult to find since it involves nonlinear kinematics and thermal fluctuations. In this paper we use ideas from the Kirchhoff theory of filaments to find semi-analytical solutions for the average shape of the fluctuating DNA under the assumption that there is no self-contact. The basic step in our method consists of combining a helical solution of the rod with a non-planar localizing solution in such a way that the force, moment, position and slope remain continuous everywhere along the rod. Our solutions allow us to predict the extension vs. linking number behavior of long pieces of DNA for various values of the tension and temperature. An interesting outcome of our calculations is the prediction of a sudden change in extension at buckling which does not seem to have been emphasized in earlier theoretical models or experiments. Our predictions are amenable to falsification by recently developed single molecule techniques which can simultaneously track the force-extension as well as the torque-rotation behavior of DNA.     

Interface motion by interface diffusion driven by bulk energy: justification of a diffusive interface model

We construct an asymptotic solution of a system consisting of the partial differential equations of linear elasticity theory coupled with a degenerate parabolic equation, and show that when a regularity parameter tends to zero, this solution converges to a solution of a sharp interface model, which describes the phase interface in an elastically deformable solid moving by interface diffusion. Therefore, the coupled system can be used as diffusive interface model. Differently from diffusive interface models based on the Cahn–Hilliard equation, the interface diffusion is solely driven by the bulk energy, hence the Laplacian of the curvature is not part of the driving force. Also, no rescaling of the parabolic equation is necessary. Since the asymptotic solution does not solve the system exactly, the proof is formal.     

Finite element methods for a class of continuum models for immiscible flows with moving contact lines

In this paper, we present a finite element method for two‐phase incompressible flows with moving contact lines. We use a sharp interface Navier–Stokes model for the bulk phase fluid dynamics. Surface tension forces, including Marangoni forces and viscous interfacial effects, are modeled. For describing the moving contact lines, we consider a class of continuum models that contains several special cases known from the literature. For the whole model, describing bulk fluid dynamics, surface tension forces, and contact line forces, we derive a variational formulation and a corresponding energy estimate. For handling the evolving interface numerically, the level‐set technique is applied. The discontinuous pressure is accurately approximated by using a stabilized extended finite element space. We apply a Nitsche technique to weakly impose the Navier slip conditions on the solid wall. A unified approach for discretization of the (different types of) surface tension forces and contact line forces is introduced. Results of numerical experiments are presented, which illustrate the performance of the solver. Copyright © 2016 John Wiley & Sons, Ltd.     

Balance equations and structural models for phase interfaces

The general balance equations are developed for an interface represented by a dividing surface and for a moving common line represented as an intersection of dividing surfaces. The surface excess variables associated with a dividing surface are expressed both in terms of those variables describing the three-dimensional interfacial region of finite thickness and in terms of those variables describing bulk phases that extend up to the dividing surface.A structural model for the interface is suggested in which a suspension of solid bodies representing surfactant molecules is distributed about a singular surface separating two adjacent bulk solvent phases. The suspension is required to have the same average behavior as the interfacial region. This is interpreted as meaning that the general jump balance for a continuum dividing surface represented by an interfacial suspension is a local area average. Specific results are derived for two structural models, each in the same simple shear field. One consists of a dilute suspension of neutrally buoyant spheres floating with their centers restricted to the dividing surface. The other is a dilute suspension of chains of neutrally buoyant spheres with the sphere at one end of the chain floating in the dividing surface.     

Dynamics of Phase Transitions and Hysteresis in a Viscoelastic Ericksen's Bar on an Elastic Foundation

This work is a follow-up on the study [32] of interface dynamics and hysteresis in materials undergoing solid-solid phase transitions. We consider the dynamics of a viscoelastic bar with a nonmonotone stress-strain relation and viscous stress linearly proportional to the strain rate. The bar is placed on an elastic foundation with stiffness β mimicking the interaction of phases in higher dimensions. Time-dependent displacement-controlled loading of the bar results in a tilted and serrated hysteresis loop, in qualitative agreement with some experimental observations in shape-memory alloys. The model exhibits three phase transition processes: phase nucleation, interface slip and phase annihilation. Between these dynamic processes the system gets stuck in local minimizers of the potential energy. As β increases from zero, a slip-dominated solution behavior transforms to the one where slip and annihilation events are preceded by a step-by-step nucleation process. We show that this transition is caused by an interplay between the slip-favoring inertia term and the nucleation-favoring elastic foundation terms. This revised version was published online in August 2006 with corrections to the Cover Date.     

Asymmetric Potentials and Motor Effect: A Large Deviation Approach

We provide a mathematical analysis for the appearance of concentrations (as Dirac masses) in the solutions to Fokker–Planck systems with asymmetric potentials. This problem has been proposed as a model to describe motor proteins moving along molecular filaments. The components of the system describe the densities of the different conformations of the proteins. Our results are based on the study of a Hamilton–Jacobi equation arising at the zero diffusion limit after an exponential transformation change of the phase function that yields a viscous Hamilton–Jacobi equation. We consider different classes of conformation transitions coefficients (bounded, unbounded and locally vanishing).     

An ALE finite element method for a coupled Stefan problem and Navier–Stokes equations with free capillary surface

In this article, an ALE finite element method to simulate the partial melting of a workpiece of metal is presented. The model includes the heat transport in both the solid and liquid part, fluid flow in the liquid phase by the Navier–Stokes equations, tracking of the melt interface solid/liquid by the Stefan condition, treatment of the capillary boundary accounting for surface tension effects and a radiative boundary condition. We show that an accurate treatment of the moving boundaries is crucial to resolve their respective influences on the flow field and thus on the overall energy transport correctly. This is achieved by a mesh‐moving method, which explicitly tracks the phase boundary and makes it possible to use a sharp interface model without singularities in the boundary conditions at the triple junction. A numerical example describing the welding of a thin‐steel wire end by a laser, where all aforementioned effects have to be taken into account, proves the effectiveness of the approach.Copyright © 2012 John Wiley & Sons, Ltd.     

Understanding the influence of structural hierarchy and its coupling with chemical environment on the strength of idealized tropocollagen–hydroxyapatite biomaterials

Hard biomaterials such as bone, dentin, and nacre have primarily an organic phase (e.g. tropocollagen (TC)) and a mineral phase (e.g. hydroxyapatite (HAP) or aragonite) arranged in a staggered arrangement at the nanoscopic length scale. Interfacial interactions between the organic phase and the mineral phase as well as the structural effects arising due to the staggered arrangement significantly affect the strength of such biomaterials. The effect of such factors is intricately intertwined with the chemical environment of such materials. In the present investigation, an idealized TC–HAP composite system under tensile loading is analyzed using explicit three-dimensional (3-D) molecular dynamics (MD) simulations to develop an understanding of these factors. The material system is analyzed in three different environments: (1) in the absence of water molecules (non-hydrated), (2) in the presence of water molecules (hydrated), and (3) in the presence of water molecules with calcium ions (ionized water). The analyses focus on understanding the correlations among factors such as the structural arrangement, the peak stress during deformation, Young's modulus, the peak interfacial strength, and the length scale of the localization of peak stress during deformation. Analyses show that maximizing the contact area between the TC and HAP phases results in higher interfacial strength as well as higher fracture strength. Due to the staggered arrangement, the orientation of HAP crystals has insignificant effect on the biomaterial strength. Analyses based on strength scaling as a function of structural hierarchy level reveal that while peak strength follows a multiscaling relation, the fracture strength does not. The peak strain for failure was found to be independent of the changes in levels of structural hierarchy. Overall, the analyses, being limited in size due to the computational time constraint, point out important correlations between the mechanical strength and chemically influenced structural hierarchy of biomaterials.     

Drying of porous oil shales

Rate processes including change of phase are modelled analytically for a half-space porous substance exposed to a jump in external temperature. The model predicts, in a closed form, the pressure build-up and the rate of evaporation of volatiles from a porous matrix. It assumes two distinct regions separated by a moving interface where the change of phase takes place. One region maintains its initial concentration of volatiles while the second is devoid of volatiles. Different thermophysical properties are considered for the two regions. The model was applied for the evaporation of moisture from oil-shale. Results are given in a parametric form     

一种考虑界面相物性影响的界面单元法

界面是由复杂的界面相简化而成的,界面破坏实际是界面相材料的破坏。数值计算为了方便,如经典模型和内聚力模型等,都把很薄的界面相作无厚度化处理。导致只能考虑界面的面力,而无法考虑界面相内的应力(平行于界面方向的应力)。使界面失效准则先天性地排除了界面相内部应力的影响,从界面相材料失效机理的角度来看这是不够严谨的。本文将界面相材料等效为一种弹性连续体,由界面本构关系推导得到了一种新的界面单元。该单元具有界面参数易确定、对界面相物性可以进行等效描述等优点。通过商用有限元软件ABAQUS和用户子程序UEL实现了数值分析,并与直接物理模型的数值模拟结果进行对比,证明了本方法的简便及准确性。通过对不同界面相厚度结构的进一步分析,探讨了本文方法的可行范围。     

模拟相变问题的精细积分边界元法

将精细积分边界元法和界面追踪法相结合求解相变问题。因为边界元法只需要将待求解空间域的边界离散,方便连续追踪移动界面位置和重构网格,所以边界元法适合应用于移动边界问题的模拟。首先,利用精细积分边界元法在固相区域和液相区域分别求解相应的瞬态热传导控制方程,从而求得温度场和边界热流密度。然后,根据固-液相变界面上的能量平衡方程,利用热流密度求得相变界面的移动速度,再采用界面追踪法预测移动相变界面的位置变化。最后,给出了几个数值算例,并通过与参考解的对比验证本文方法的准确性。     

水的汽液界面系统中势能与力的EMD模拟研究

分子间势能作用是研究分子界面行为的一个重点所在. 采用平衡态分子动力学 模拟(equilibrium molecular dynamics simulation, EMDS)方法,对由水分子 构成的汽液界面系统进行了模拟和研究. 分析统计结 果符合势能分布在液相区和气相区内存在明显落差的已知结论,并发现不同种力对分子穿越 两相区时所起的作用不同,Lennard-Jones(简写L-J)力阻碍分子凝结,而静电力则推动分子凝结并且在合力中起 主要作用. 同时,着重对发生相变行为的典型分子进行了追踪和分析,从能量的角度显 示了凝结(蒸发)相变过程对应着一个气态(液态)分子由高(低)势能位落入势阱(翻越 势垒)的能量降落(抬升)过程.     

Stable two-layer flows at all Re; visco-plastic lubrication of shear-thinning and viscoelastic fluids

Multi-fluid flows are frequently thought of as being less stable than single phase flows. Consideration of different non-Newtonian models can give rise to different types of hydrodynamic instability. Here we show that with careful choice of fluid rheologies and flow paradigm, one can achieve multi-layer flows that are linearly stable for Re = ∞. The basic methodology consists of two steps. First we eliminate interfacial instabilities by using a yield stress fluid in one fluid layer and ensuring that for the base flow configurations studied we maintain an unyielded plug region at the interface. Secondly we eliminate linear shear instabilities by ensuring a strong enough Couette component in the second fluid layer, imposed via the moving interface. We show that this technique can be applied to both shear-thinning and visco-elastic fluids.     

Numerical solution of solidification in a square prism using an algebraic grid generation technique

The solidification of an infinitely long square prism was analyzed numerically. A front fixing technique along with an algebraic grid generation scheme was used, where the finite difference form of the energy equation is solved for the temperature distribution in the solid phase and the solid–liquid interface energy balance is integrated for the new position of the moving solidification front. Results are given for the moving solidification boundary with a circular phase change interface. An algebraic grid generation scheme was developed for two-dimensional domains, which generates grid points separated by equal distances in the physical domain. The current scheme also allows the implementation of a finer grid structure at desired locations in the domain. The method is based on fitting a constant arc length mesh in the two computational directions in the physical domain. The resulting simultaneous, nonlinear algebraic equations for the grid locations are solved using the Newton-Raphson method for a system of equations. The approach is used in a two-dimensional solidification problem, in which the liquid phase is initially at the melting temperature, solved by using a front-fixing approach. The difference of the current study lies in the fact that front fixing is applied to problems, where the solid–liquid interface is curved such that the position of the interface, when expressed in terms of one of the coordinates is a double valued function. This requires a coordinate transformation in both coordinate directions to transform the complex physical solidification domain to a Cartesian, square computational domain. Due to the motion of the solid–liquid interface in time, the computational grid structure is regenerated at every time step.     

Two-dimensional Green's functions in anisotropic multiferroic bimaterials with a viscous interface

We derive, by virtue of the unified Stroh formalism, the extremely concise and elegant solutions for two-dimensional and (quasi-static) time-dependent Green's functions in anisotropic magnetoelectroelastic multiferroic bimaterials with a viscous interface subjected to an extended line force and an extended line dislocation located in the upper half-plane. It is found for the first time that, in the multiferroic bimaterial Green's functions, there are 25 static image singularities and 50 moving image singularities in the form of the extended line force and extended line dislocation in the upper or lower half-plane. It is further observed that, as time evolves, the moving image singularities, which originate from the locations of the static image singularities, will move further away from the viscous interface with explicit time-dependent locations. Moreover, explicit expression of the time-dependent image force on the extended line dislocation due to its interaction with the viscous interface is derived, which is also valid for mathematically degenerate materials. Several special cases are discussed in detail for the image force expression to illustrate the influence of the viscous interface on the mobility of the extended line dislocation, and various interesting features are observed. These Green's functions can not only be directly applied to the study of dislocation mobility in the novel multiferroic bimaterials, they can also be utilized as kernel functions in a boundary integral formulation to investigate more complicated boundary value problems where multiferroic materials/composites are involved.     

Freezing over a cylinder with finite heat transfer at the wall

A simple and straightforward implicit numerical technique for phase change problems in cylindrical geometry is proposed. When implicit schemes are used the moving interface along with the convective or radiative boundary condition pose a problem because of the requirement to calculate the interface location and boundary temperature implicitly. Due to this difficulty some of the methods available in literature used approximations near the wall and the rest used iterative methods. The present technique isolates the nonlinearity associated with the moving interface as well as that of any nonlinear wall boundary condition and permits the simultaneous evaluation of the unknown interface location and the wall temperature at the new time level. Thereafter the temperatures at the other nodes are obtained without nodal iterations. Numerical results are obtained for outward solidification with both finite and infinite heat transfer rates at the wall. Good agreement is obtained between the present numerical results and the results available in literature for the limiting cases.     

A New Model and its Application to Investigate Transpiration Cooling with Liquid Coolant Phase Change

This paper presents a new mathematical model, semi-mixing model (SMM), to describe transpiration cooling with coolant phase change from liquid into vapor through two-phase process. The local heat exchange of fluid-solid within pores is considered in this model, and therefore it is closer to real transpiration cooling condition. The differences from the separated phase model and two-phase mixture model are that SMM can overcome the trouble of tracking phase change interface and avoid the inveracious numerical phenomenon, i.e., a thermal insulating layer occurs within the porous matrix. Using SMM, the corresponding numerical method is realized to simulate the processes of coolant moving, absorbing heat and evaporating within porous matrix. To validate SMM and the numerical strategy, an experiment is conducted. Using the validated SMM and numerical strategy, the effects of two-dimensional coolant injection rate and two-dimensional heat flux on transpiration cooling characteristics are simulated and analyzed. The simulations and analysis discover several interesting and valuable phenomena.     

On surface tension modelling using the level set method

The paper describes and compares the performance of two options for numerically representing the surface tension force in combination with the level set interface‐tracking method. In both models, the surface tension is represented as a body force, concentrated near the interface, but the technical implementation is different: the first model is based on a traditional level set approach in which the force is distributed in a band around the interface using a regularized delta function, whereas in the second, the force is partly distributed in a band around the interface and partly localized to the actual computational cells containing the interface. A comparative study, involving analysis of several two‐phase flows with moving interfaces, shows that in general the two surface tension models produce results of similar accuracy. However, in the particular case of merging and pinching‐off of interfaces, the traditional level set model of surface tension produces an error that results in non‐converging solutions for film‐like interfaces (i.e. ones involving large contact areas). In contrast, the second model, based on the localized representation of the surface tension force, displays consistent first‐order convergence. Copyright © 2008 John Wiley & Sons, Ltd.     

Cahn–Hilliard/Navier–Stokes Model for the Simulation of Three-Phase Flows

In this article, we describe some aspects of the diffuse interface modelling of incompressible flows, composed of three immiscible components, without phase change. In the diffuse interface methods, system evolution is driven by the minimisation of a free energy. The originality of our approach, derived from the Cahn–Hilliard model, comes from the particular form of energy we proposed in Boyer and Lapuerta (M2AN Math Model Numer Anal, 40:653–987,2006), which, among other interesting properties, ensures consistency with the two-phase model. The modelling of three-phase flows is further completed by coupling the Cahn–Hilliard system and the Navier–Stokes equations where surface tensions are taken into account through volume capillary forces. These equations are discretized in time and space paying attention to the fact that most of the main properties of the original model (volume conservation and energy estimate) have to be maintained at the discrete level. An adaptive refinement method is finally used to obtain an accurate resolution of very thin moving internal layers, while limiting the total number of cells in the grids all along the simulation. Different numerical results are given, from the validation case of the lens spreading between two phases (contact angles and pressure jumps), to the study of mass transfer through a liquid/liquid interface crossed by a single rising gas bubble. The numerical applications are performed with large ratio between densities and viscosities and three different surface tensions.