狭窄流域内气泡破裂现象数值模拟

基于边界积分法建立狭窄流域内气泡破裂数值模型,开发相应的计算程序,分别模拟对称破裂与非对称破裂两类典型工况并与相应实验结果进行对比,计算值与实验值吻合很好,表明三维数值模型的有效性．从狭窄流域内气泡运动的基本现象入手,基于开发的程序系统地研究气泡的对称破裂与非对称破裂,在已有数值研究成果和实验数据基础上,提出气泡破裂的可行性准则,研究分裂后子气泡的动力学特性,并分析距离参数对气泡破裂特性及子气泡动力学行为的影响,总结相应规律,旨在为相关气泡破裂特性研究提供参考．     

水下爆炸气泡与复杂弹塑性结构的相互作用研究

计及结构的弹塑性,将边界元法(BEM)与有限元法(FEM)耦合提出了气泡与弹塑性结构耦合动力学计算方法,并开发了全套的三维水下气泡分析程序（UBA）,计算值与实验值之间误差在10%以内．以水面舰船为例,将三维计算程序工程化．并分析了水下爆炸气泡载荷作用下船体的弹塑性响应,从船体结构典型单元上的应力时历曲线可以看出,在气泡坍塌时出现应力峰值,证实了气泡坍塌压力及射流引起的压力对舰船等结构造成严重毁伤．从气泡与舰船的相互作用中可以看出,舰船低阶垂向振型被激起,在气泡作用下呈鞭状运动,同时舰船随着气泡的膨胀和收缩作升沉运动,通过本文的分析得到了适合于工程应用的规律及结论．     

泡沫覆盖层对水下爆炸气泡射流防护机理缩比试验研究

水下爆炸会对水中结构物造成严重威胁。柔性覆盖层或夹层板能够降低水中结构物水下爆炸冲击响应，因此成为研究的热点。以往的研究多集中于覆盖层对冲击波的防护机理，适用于较远距离的水下爆炸情况。近距离水下爆炸除了冲击波外，爆炸气泡溃灭时产生的朝向结构物的高速水射流更为致命。该文针对这种情况，基于量纲原理，推导缩比相似关系，通过缩比模型水下爆炸试验发现了覆盖层表面空化微气泡群对爆炸气泡形成高速水射流过程产生干扰，提出了泡沫覆盖层钢板水下爆炸气泡射流防护机理。     

水下爆炸异相气泡动力学特性的Euler有限元数值模拟研究北大核心CSCD

该研究基于Euler有限元方法建立了水下爆炸异相气泡动力学模型,将计算结果和气泡统一理论以及异相爆炸试验进行对比验证计算模型的有效性.通过和自由场单气泡进行对比分析,发现异相爆炸冲击波对气泡做功是异相爆炸气泡总能量增加的原因,相位差绝对值越接近π、距离参数越小气泡总能量损失越少,后产生的气泡会使先产生的气泡提前坍缩.气泡的射流方向受相位差影响,相位差为零时产生对向射流,其他相位差产生反向射流.     

计入气穴的水下爆炸作用下结构流固耦合三维数值模拟

水下爆炸在结构物面附近产生的气穴现象,严重影响水下爆炸作用下的流固耦合动响应,是舰船水下爆炸领域的难点,传统的边界元方法、有限元方法(FEM)难以解决水下爆炸气穴现象这类强非线性问题．针对此问题,计及流体中的气穴现象,考虑流体的可压缩型,忽略流体粘性,建立水下爆炸瞬态强非线性流固耦合三维数值模型,采用流体谱单元方法(SEM)和结构有限元方法求解该模型．计算结果表明:相对有限元法,谱单元法具有更高的计算精度,且谱单元解与解析解、试验值吻合良好．在此基础上,结合ABAQUS软件,分别探讨三维球壳、船体板架在水下爆炸作用下的瞬态流固耦合机理,给出气穴区域及其对水中结构物动响应的影响特征,旨在为舰船水下爆炸瞬态流固耦合问题的相关研究提供参考．     

梢涡流场中气泡轨迹及形态数值模拟

针对气泡在舰船尾迹涡流场运动特性,根据其是否为尾涡所捕获,将数值模拟过程分为两个阶段:准球状运动阶段和非球状运动阶段．分别应用单向耦合质点粒子追踪法（PTM）和边界元法（BEM）模拟这两个阶段,将第1阶段结束的物理量作为第2阶段的初始条件,从而完成整个数值模拟过程．在已有数值研究结果和实验数据基础上,探讨空化发生条件,追踪尾迹空泡运动轨迹,模拟尾迹气泡的运动、变形、溃灭等,以及被尾涡捕获后的撕裂等运动特性,旨在为优化设计尾流场提供参考．     

强激波与流-固界面相互作用的数值模拟

1 引言强激波与流-固界面相互作用的数值模拟在实际问题中具有重要的应用,比如水下炸弹爆炸时产生的强大的冲击波对附近的舰船或水下潜艇会造成非常大的破坏,造成舰体严重变形甚至断裂．在此类问题中,界面两边的流体具有完全不同的特性,流体的密度、     

钢管混凝土柱抗爆性能数值模拟与实验验证

利用LS-DYNA非线性有限元程序,基于多物质流固耦合方法,建立了爆炸荷载作用下钢管混凝土柱的动态响应数值模型.对比分析了模拟结果与足尺构件的爆炸破坏实验结果,验证了数值模型和计算方法的有效性,并运用参数化分析方法,研究了截面形式、比例距离、混凝土强度及钢材等级、截面形状特性等关键参数对钢管混凝土柱抗爆性能的影响.研究结果表明:钢管混凝土柱具有优越的抗爆性能,所建立的数值模型能够有效地分析钢管混凝土柱在爆炸荷载作用下的动态影响及破坏形态;圆形截面钢管混凝土柱的抗爆性能优于方形截面;提高材料等级、减小圆形截面钢管混凝土柱的径厚比、增大矩形截面钢管混凝土柱的长宽比,均有利于提升钢管混凝土柱的抗爆特性.     

应用激光全息干涉术及有限元法对2050热连轧机机座自振特性的三维分析

根据相似理论,进行了宝山钢铁总厂热轧厂2050热连轧精轧机机座自振特性的三维模拟研究,并应用全息干涉术进行了前三阶主振型分析.本文还进行了三维有限元数值分析与模拟值进行了映照.二者结果是吻合的.     

基于生态模型的企业竞争稳定性研究

从生态学视角来研究企业竞争的动态演化.针对现有研究文献中只讨论企业自身线性制约的三维系统的不足,构造了企业自身非线性制约的三维竞争模型,运用微分方程稳定性理论分析其稳定性,并由此揭示出企业间竞争的动力学机制.数值仿真结果表明,该三维模型能有效地模拟企业间竞争的动态演化规律.     

Vortex shedding behind a rising bubble and two-bubble coalescence: A numerical approach

In this work, we present the computational results on the wake instability in wobbling bubble regime as well as on the coalescence of two bubbles in different shape regimes. This is a continuation of our previous studies on the dynamics of a single gas bubble rising in a viscous liquid (see [A. Smolianski, H. Haario, P. Luukka, Computational Study of Bubble Dynamics, Research Report 86, Lappeenranta University of Technology, Finland]), and we use the same, finite-element/level-set/operator-splitting method that was proposed in [A. Smolianski, Numerical Modeling of Two-Fluid Interfacial Flows, Ph.D. Thesis, University of Jyväskylä, 2001]. The numerical method allows to simulate a wide range of flow regimes, accurately capturing the shape of the deforming interface of the bubble and the surface tension effect, while maintaining a good mass conservation. Due to the highly unstable and small-scale nature of the considered problems there are very few experimental investigations, but the comparison with available experimental data confirms a good accuracy of our numerical predictions. Our studies show that plausible results can be obtained with two-dimensional numerical simulations, when a single buoyant bubble or a coalescence of two bubbles is considered.     

Numerical Modeling of Laser-Induced Cavitation Bubbles with a Finite Volume Method

We report on the numerical modeling of laser-induced cavitation bubbles with the finite volume method using the open source software package OpenFOAM. The numerical model is validated by comparison to experimental data for the two cases of a bubble collapsing under normal ambient conditions in an unbounded liquid as well as close to a solid wall. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of single gas bubbles under shear flow conditions

Disperse gas bubbles play an important role in many industrial applications. Knowing the rising velocity, the interfacial area, or the critical size for break-up or coalescence in different systems can be crucial for the process design. Usually the flow experienced by bubbles is not uniform but sheared. Under shear-flow conditions bubbles develop a lift force perpendicular to the flow direction. In the present work direct numerical simulations are applied to examine the dependency of the lift force on the shear rate for bubbles in pure liquids. A level-set based volume-tracking method is implemented into the CFD-code OpenFOAM, to follow the free interface of the gas bubble. Results show good agreement with available experimental results from single bubbles in a rotating chamber. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Large eddy simulation and PIV experiments of air–water mixing tanks

The simulations and experiments of a turbulent bubbly flow are carried out in a cylindrical mixing vessel. Dynamics of the turbulent bubbly flow is visualized using a novel two-phase particle image velocimetry (PIV) with a combination of back lighting, digital masking and fluorescent tracer particles. Using an advanced technique, Mie’s scattering at surfaces of bubbles is totally filtered out and, henceforth, images of tracer particles and of bubbles are obtained with high quality. In parallel to the comprehensive experimental studies, numerical results are obtained from large eddy simulations (LES) of the two-phase air–water mixer. The impeller-induced flow at the blade tip radius is modeled by using sliding mesh method. The results demonstrate the existence of large structures such as tip-vortex tips, and also some finer details. In addition, the stability of the jet is found to be connected with the fluctuations of the tip vortices whose dynamics are affected by the presence of bubbles. Numerical results are used to interpret the measurement data and to guide the refinement of consistent theoretical analyses. Such information is invaluable in the development of advanced theories capable of describing bubbly flows in the presence of complex liquid flow. This detailed information is of real significance in facilitating the design and scale-up of practical stirred tanks.     

Numerical study of dynamics of single bubbles and bubble swarms

A systematic computational study of the dynamics of gas bubbles rising in a viscous liquid is presented. Two-dimensional simulations are carried out. Both the dynamics of single bubbles and small groups of bubbles (bubble swarms) are considered. This is a continuation of our previous studies on the two-bubble coalescence and vortex shedding [A. Smolianski, H. Haario, P. Luukka, Vortex shedding behind a rising bubble and two-bubble coalescence: a numerical approach, Appl. Math. Model. 29 (2005) 615–632]. The proposed numerical method allows us to simulate a wide range of flow regimes, accurately capturing the shape of the deforming interface of the bubble and the surface tension effect, while maintaining the mass conservation. The computed time-evolution of bubble’s position and rise velocity shows a good agreement with the available experimental data. At the same time, the results on the dynamics of bubble interface area, which are, up to our knowledge, presented for the first time, show how much the overall mass transfer would be affected by the interface deformation in the case of the bubble dissolution. Another set of experiments that are of interest for chemical engineers modelling bubbly flows concerns the bubble swarms and their behavior in different bubble-shape regimes. The ellipsoidal and spherical shape regimes are considered to represent, respectively, the coalescing and non-coalescing bubble swarms. The average rise velocities of the bubble swarms are computed and analyzed for both regimes.     

Numerical investigations of bubbles rising in water

The hydrodynamics of rising air bubbles in water filled bubble columns are examined using a transient numerical model following the Euler-Euler method. In this particular case the bubble-size distribution is assumed to be mono-disperse and the influence of this assumption on the specific mass-transfer area is discussed. Additionally, the influence of the bubble size is discussed with regard to the mass-transfer area, the air volume-fraction, and the axial air velocity profile. To validate the numerical model, experimental results from tomographic measurements at the Helmholtz-Zentrum Dresden-Rossendorf (HZDR) are used. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

基于时间序列的混合神经网络数据融合算法

针对传统的数据融合算法对高噪声、大规模、数据结构复杂的时间序列数据融合性能较差的问题,该文提出了一种混合神经网络的数据融合算法(即SCLG算法).SCLG算法的思想是首先利用奇异谱分析算法对数据分解重构以达到去噪的目的;其次,通过深层卷积神经网络提取数据的空间特征和短期时间特征;然后,利用长短期记忆(LSTM)网络和门控循环单元(GRU)网络双层网络,进一步深度提取数据时间维度上的特征;最后,利用全连接网络综合主要信息输出最终的决策.通过SP&500和AQI数据集上的实验结果表明,该算法在融合性能及稳定性方面均优于DCNN、CNN LSTM、FDL数据融合算法.     

Application of a multi-component Eulerian approach to capture interface evolution and vorticity generation in compressible multiphase flow

The hyperbolic Eularian model is used as a mathematical framework for compressible multiphase flows. The formulation was obtained after an averaging process of the single phase Navier-Stokes equations. The closure of multi-component system leads to the volume fraction equation containing a non-conservative term and a pressure equilibrium condition. As a result the model equations cannot be written in a conservative form. To solve the equations a finite volume Godunov type computational approach is developed which uses an approximate Riemann solver combined with a numerical scheme to tackle the non-conservative terms. The approach accounts for pressure non-equilibrium. It enables resolving interfaces separating compressible fluids and captures the baroclinic source of vorticity generation. The computations are performed for various initial conditions and compared with theoretical and experimental data for a shock-bubble interaction problem. The investigated cases include acoustic wave transmission through isolated bubbles of helium and krypton. The numerical results illustrate the characteristic features of the evolving interfaces. The impulsively generated flow perturbations are dominated by the reflection and refraction of the shock and by the vorticity generation within the media. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Measurement of surface and interfacial tension from the maximum pressure in sessile and pendant bubbles and drops

The measurement of surface and interfacial tension from the maximum pressure in sessile and pendant bubble and drops is, described. New theoretical data are presented which enable the surface and interfacial tension to be calculated without iteration. Although the maximum pressure in pendant bubbles is less than that in sessile bubbles, they retain their symmetry during detachment from the capillary tip and so havt experimental advantages. Reproducible results in good agreement with literature values are obtained for the surface and interfacial tension of several liquids and immiscible liquid pairs.