Gas-dynamic signs of explosive eruptions of volcanoes. 2. Model of homogeneous-heterogeneous nucleation. Specific features of destruction of the cavitating magma

The dynamics of state of the crystallite-containing magma is studied within the framework of the gas-dynamic model of bubble cavitation. The effect of crystallites on flow evolution is considered for two cases: where the crystallites are cavitation nuclei (homogeneous-heterogeneous nucleation model) and where large clusters of crystallites are formed in the magma in the period between eruptions. In the first case, decompression jumps are demonstrated to arise as early as in the wave precursor; the intensity of these jumps turns out to be sufficient to form a series of discrete zones of nucleation ahead of the front of the main decompression wave. Results of experimental modeling of an explosive eruption with ejection of crystallite clusters (magmatic “bombs”) suggest that a cocurrent flow of the cavitating magma with dynamically varying properties (mean density and viscosity) transforms to an independent unsteady flow whose velocity is greater than the magma flow velocity. Experimental results on modeling the flow structure during the eruption show that coalescence of bubbles in the flow leads to the formation of spatial “slugs” consisting of the gas and particles. This process is analyzed within a combined nucleation model including the two-phase Iordansky-Kogarko-van Wijngaarden model and the model of the “frozen” field of mass velocities in the cavitation zone. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 2, pp. 167–177, March–April, 2009.     

Limit analysis of multi-layered plates. Part II: Shear effects

In the first part of this work [Dallot, J., Sab, K., 2007. Limit analysis of multi-layered plates. Part I: the homogenized Love-Kirchhoff model. J. Mech. Phys. Solids, in press, doi:10.1016/j.jmps.2007.05.005], the limit analysis of a multi-layered plastic plate submitted to out-of-plane loads was studied. The authors have shown that a homogeneous equivalent Love-Kirchhoff plate can be substituted for the heterogeneous multi-layered plate, as the slenderness (length-to-thickness) ratio goes to infinity. In fact, the out-of-plane shear stresses are shown to become asymptotically negligible when compared to in-plane stresses, as the slenderness ratio goes to infinity. Actually, failure of thick multi-layered structures often occurs by shearing in the core layers and sliding at the interfaces between the layers. Both shearing and sliding are caused by the out-of-plane shear stresses. The purpose of the present paper is to build an enhanced Multi-particular Model for Multi-layered Material (M4) taking into account shear stress effects. In this model, each layer is seen as a Reissner-Mindlin plate interacting with its neighboring layers through interfaces. The proposed model is asymptotically consistent with the homogenized Love-Kirchhoff model described in the first part of the work, as the slenderness ratio goes to infinity. Kinematic and static methods for the determination of the limit load of a thick multi-layered plate which is submitted to out-of-plane distributed forces are described. The special case of multi-layered plates under cylindrical bending conditions is studied. These conditions lead to simplifications which often allow for the analytical resolution of the Love-Kirchhoff and the M4 limit analysis problems. The benefit of the proposed M4 model is demonstrated on an example. A comparison between the heterogeneous 3D model, the Love-Kirchhoff model and the M4 model is performed on a three-layer sandwich plate under cylindrical bending conditions. Finite element calculations are used to solve the 3D problem, while both the Love-Kirchhoff and the M4 problems are analytically solved. It is shown that, when the contrast between the core and the skins strengths is high, the Love-Kirchhoff model fails to capture the plastic collapse modes that cause the ruin of the sandwich plate. These modes are well captured by the M4 model which predicts limit loads that are very consistent with the limit loads predicted by the heterogeneous 3D model (the relative error is found to be smaller than 1%).     

关于二相流、多相流、多流体模型和非牛顿流等概念的探讨

本文分析了单相流、二相流和多相流等概念上的差异,也分析了单流体模型、双流体模型和多流体模型等概念上的差异,指出前面三种概念是按流动介质的客观物理构成划分的,而后者是按主观采用的研究方法划分的．目前这些概念在使用中存在一些混乱,如二相流与多相流,多相流与多流体模型等．本文还研究了扩散模型、非牛顿流模型和颗粒流模型等,指出前两种模型在分类上属于单流体模型,分析了非牛顿流模型、扩散模型和双（多）流体模型的特点和应用范围,最后,以泥石流为例讨论了以上概念的应用．      

RESEARCH ON AMSAA-BISE MODEL—AN ANSWER TO MEI MEN-HUA

The AMSAA-BISE model is derived from another approach. This certainly shows the correctness of the AMSAA-BISE model, and indicates the incorrectness of the approximate model given in this paper. The engineering example illustrating these conclusions is given. Merits and demerits of AMSAA and AMSAA-BISE model are discussed.     

RESEARCH ON AMSAA-BISE MODEL—AN ANSWER TO MEI MEN-HUA

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎＭＥＩＷｅｎ＿ｈｕａ’ｓｐａｐｅｒ [1 ] ,ａｃｏｍｍｅｎｔｏｎＡＭＳＡＡ＿ＢＩＳＥｍｏｄｅｌｐｒｏｐｏｓｅｄｂｙａｕｔｈｏｒｅｔａｌ.(ｃｆ[2 ] )ｉｓｇｉｖｅｎ ,ａｎｄａｎａｐｐｒｏｘｉｍａｔｅｍｏｄｅｌａｔｔｅｍｐｔｉｎｇｔｏｓｏｌｖｅｔｈｉｓｐｒｏｂｌｅｍｉｓｐｒｅｓｅｎｔｅｄ .ＴｈｅＡＭＳＡＡ＿ＢＩＳＥｍｏｄｅｌｉｓｄｅｒｉｖｅｄｄｉｒｅｃｔｌｙｆｒｏｍａｎｏｔｈｅｒａｐｐｒｏａｃｈｔｈａｔｉｓｄｉｆｆｅｒｅｎｔｆｒｏｍｔｈｅａｐｐｒｏａｃｈｅｓｏｆ[2 ]ａｎｄ [3 ]…     

RESEARCH ON AMSAA-BISE MODEL—AN ANSWER TO MEI WEN-HUA

The AMSAA-BISE model is derived from another approach. This certainly shows the correctness of the AMSAA- BISE model, and indicates the incorrectness of the approximate model given in this paper. The engineering example illustrating these conclusions is given.Merits and demerits of AMSAA and AMSAA-BISE model are discussed.     

无壁面参数低雷诺数非线性涡黏性模式研究

建立了一个低雷诺数的非线性涡黏性湍流模式,该模式的一个显著特性是它不包含壁面参数（如y^ ,n等）,因而特别适用于复杂几何流场的计算,本模式在几种包括回流、分离、激波等典型流动中进行了验证,结果令人满意。     

非完整系统动力学中的Vakonomic模型和Четаев模型

综述了非完整系统动力学中Чｅтａｅв模型和Ｖａｋｏｎｏｍｉｃ模型．论证了Ｈｅｒｔｚ,Ｃａｐｏｎ和Ｌｉｎｄｅｌｏｆ的工作与Ｖａｋｏｎｏｍｉｃ模型相吻合,而Ｈｏｌｄｅｒ,Ｐａｒｓ和Чаплыгин的工作与Четаев模型相吻合;Ｖａｋｏｎｏｍｉｃ模型和четаев模型两类模型并存的局面可以追溯到非完整力学发展的初期．说明了伴随两类模型并存的现象而出现的两个问题:ｄ－δ交换性问题和∫－δ交换性问题．综述了我国学者为统一两类模型作出的贡献．最后,展望了该领域的研究方向．      

Multi‐scale simulation with a hybrid Boussinesq‐RANS hydrodynamic model

A hybrid wave model is developed for simulation of water wave propagation from deep water to shoreline. The constituent wave models are the irrotational, 1‐D horizontal Boussinesq and 2‐D vertical Reynolds‐averaged Navier–Stokes (RANS). The models are two‐way coupled, and the interface is placed at a location where turbulence is relatively small. Boundary conditions on the interfacing side of each model are provided by its counterpart model through data exchange. Prior to the exchange, a data transformation step is carried out due to the differences in physical variables and approximations employed in both models. The hybrid model is tested for both accuracy and speedup performance. Tests consisting of idealized solitary and standing wave motions and wave overtopping of nearshore structures show that: (1) the simulation results of the current hybrid model compare well with the idealized data, experimental data, and pure RANS model results and (2) the hybrid model saves computational time by a factor proportional to the reduction in the size of the RANS model domain. Finally, a large‐scale tsunami simulation is provided for a numerical setup that is practically unapproachable using RANS model alone; not only does the hybrid model offer more rapid simulation of relatively small‐scale problems, it provides an opportunity to examine very large total domains with the fine resolution typical of RANS simulations. Copyright © 2009 John Wiley & Sons, Ltd.     

航天器动力学模型试验验证技术研究进展

耦合载荷分析是航天器研制过程中的一个重要环节, 能够为航天器结构设计, 地面验证试验条件制定以及批准型号发射提供重要依据, 而精确地得到试验验证的航天器动力学分析模型是开展耦合载荷分析的基础. 对于大型复杂航天器结构系统, 动力学模型的试验验证需要统筹安排初始建模、动力学试验、相关分析和模型修正等工作, 这是一项极具挑战的任务. 本文首先给出了结构动力学模型验证的基本流程, 然后重点讨论动力学试验, 相关分析与模型修正等关键技术, 最后结合工程实际的需求, 提出了今后航天器动力学模型试验验证技术研究的重点.      

椭球凸模型非概率可靠性度量和区间安全系数的关系

研究了椭球凸模型非概率可靠性度量和区间模型安全系数的关系。根据基于椭球凸模型的非概率可靠性指标和非概率可靠度的定义,建立了两者之间的函数关系;按照区间模型安全系数的定义,给出了由椭球参数确定的3种区间模型安全系数,分析了它们的意义;建立了非概率可靠性指标和区间模型安全系数之间的解析关系,讨论了它们在评估结构可靠性或安全程度上的意义;通过数值算例验证了分析结果。     

Immiscible Displacement in the Interacting Capillary Bundle Model Part I. Development of Interacting Capillary Bundle Model

An interacting capillary bundle model is developed for analysing immiscible displacement processes in porous media. In this model, pressure equilibration among the capillaries is stipulated and capillary forces are included. This feature makes the model entirely different from the traditional tube bundle models in which fluids in different capillaries are independent of each other. In this work, displacements of a non-wetting phase by a wetting phase at different injection rates were analysed using the interacting capillary bundle model. The predicted evolutions of saturation profiles were consistent with both numerical simulation and experimental results for porous media reported in literature which cannot be re-produced with the non-interacting tube bundle models.     

再关于非完整力学——答《争议》

本文进一步讨论了非完整力学的Chetaev模型和Vacco模型,得到了这两模型并不是合理模型的结论,同时对文[9]的争议作了回答。     

无磨损摩擦的原子理论

对描述原子级光滑表面之间的无磨损界面摩擦的微观过程的2个理论模型进行了综述,探讨了Frenkel-Kontorova(FK)模型的非线性静态和动态性质的复杂性,根据FK模型,静态特性在很大程度上依赖于表面晶格常数间的公度性,而动态特性则可分为周期解、准周期解和混沌解加以讨论。     

Multi-grain analysis versus self-consistent estimates of ferroelectric polycrystals

Ferroelectrics are polycrystalline materials consisting of intragranular regions with different polarization directions, called domains. The domains can be switched into different states by the application of an electric field or mechanical stress. We study the influence of grain-to-grain interactions on the overall and local switching behavior. The behavior inside each grain is represented by the micromechanics model of Huber et al. [1999. A constitutive model for ferroelectric polycrystals. J. Mech. Phys. Solids 47 (8), 1663-1697]. The predictions of a self-consistent model of the polycrystal response are compared with those of a multi-grain model in which grains are represented individually. In one flavor of the multi-grain model, each grain is represented by a single finite element, while in the other the fields inside each grain are captured in more detail through a fine discretization. Different electrical and mechanical loading situations are investigated. It is found that the overall response is only mildly dependent on the accuracy with which grain-to-grain interactions are captured, while the distribution of grain-average stresses is quite sensitive to the resolution of the intragranular fields.     

基于J-C本构模型的Comp.B炸药落锤冲击数值模拟

为研究Comp.B炸药在惯性冲击下的力学响应特性, 基于Johnson-Cook本构模型, 对 点火前的Comp.B炸药大落锤(400kg)冲击实验进行了数值模拟. 在不同应力加载条件下模 拟计算得到了用状态方程和不用状态方程的$\sigma$-$t$曲线, 对比应力峰值、应 力上升时间等主要力学参数, 结果表明考虑了状态方程的数值模 拟结果与实验吻合最好, 可以很好地模拟炸药惯性冲击下的力学响应.     

A mesoscopic tube model of polymer/layered silicate nanocomposites

Rheology of isothermal suspensions of completely exfoliated silicate lamellae in polymer melts is investigated. In order to express more faithfully the physics involved in low shear rates and low frequencies, we model the polymer molecules composing the melt as chains whose motion is confined to a tube formed by surrounding chains and lamellae. In the absence of lamellae, the model reduces to the mesoscopic model of reptating chains developed in Eslami and Grmela (Rheol Acta, 2008). If the chains are seen only as FENE-P dumbbells, the model reduces to the model developed in Eslami et al. (J Rheol 51:1189–1222, 2007). Responses to oscillatory, transient, and steady shear flows are calculated and compared with available experimental data. Particular attention is payed to the region of low shear rate and low frequency.     

优化方法用于模型校核

在工程实际中,常常会遇到需要选择本构模型和校准参数的模型校核问题。文中采用优化方法实现了复杂模型校核的自动化。文中,首先将各种本构模型映射为整数编号,并根据这些编号来激活相应的本构模型和本构模型参数,解决了本构模型的参数化问题;另外,根据模型校核问题的特点选择了试验结果和计算结果相对误差的平方和作为目标函数。通过这些步骤,用实现了优化方法模型校核自动化。     

Turbulence and cavitation models for time-dependent turbulent cavitating flows

Cavitation typically occurs when the fluid pressure is lower than the vapor pressure at a local thermodynamic state,and the flow is frequently unsteady and turbulent.To assess the state-of-the-art of computational capabilities for unsteady cavitating flows,different cavitation and turbulence model combinations are conducted.The selected cavitation models include several widely-used models including one based on phenomenological argument and the other utilizing interface dynamics.The kε turbulence model with additional implementation of the filter function and density correction function are considered to reduce the eddy viscosity according to the computed turbulence length scale and local fluid density respectively.We have also blended these alternative cavitation and turbulence treatments,to illustrate that the eddy viscosity near the closure region can significantly influence the capture of detached cavity.From the experimental validations regarding the force analysis,frequency,and the cavity visualization,no single model combination performs best in all aspects.Furthermore,the implications of parameters contained in different cavitation models are investigated.The phase change process is more pronounced around the detached cavity,which is better illus-trated by the interfacial dynamics model.Our study provides insight to aid further modeling development.     

Recent progress on finite element model updating: From linearity to nonlinearity

有限元分析在实际工程中得到了广泛应用.然而有限元模型由于受到网格划分、边界条件和材料物理参数不确定性等的影响,与真实结构有差异. 因此须通过试验数据加以修正,使其尽可能接近实际结构,以保证之后的结构动力模拟分析和监测等具有实际意义. 经过多年发展,有限元模型修正技术已经能够成功应用于一些实际工程,但现代工程技术的进步对有限元模型修正提出了更高要求,修正后的有限元模型不仅要有较高的精确度,还需要为后续应用给出具有指导意义的置信度.而现有的有限元模型修正、确认方法多基于结构线性的假设,而未能考虑实际结构中广泛存在的非线性.因此本文以土木工程结构模型修正的一些研究成果为例,通过对传统有限元模型修正的发展历程进行全面回顾;总结评述传统有限元修正技术的主要方法,以及包括有限元模型确认在内的最新研究进展;重点探讨有限元模型修正技术向非线性发展的技术路线和目前主要研究成果,展望其未来发展方向, 并提出值得研究的问题.