适用欧拉-拉格朗日方法模拟气液泡状流的气泡破碎模型

欧拉-拉格朗日方法已被广泛应用于模拟鼓泡塔等气-液反应器内的流型、气泡尺寸(或气含率)及其分布。文献中该方法主要基于临界Weber数观点来描述气泡破碎行为,且破碎后的子气泡尺寸由随机数确定。但现有实验和理论研究表明,临界Weber数约束不能体现气体密度等物性参数和泡内气体重分布对气泡破碎行为的影响。针对这些不足,提出了适用欧拉-拉格朗日框架且考虑泡内气体重分布贡献的气泡破碎机理模型,并利用开源软件OpenFOAM开发了基于新破碎模型的求解器。新模型预测结果能较好地吻合实验测量的时均轴向液速、气泡尺寸及其分布等实验数据。特别地,考虑泡内气体重分布现象的破碎机理模型成功预测了实验观测的气泡尺寸双峰分布特征。

基于任意拉格朗日欧拉方法的射流破土仿真分析

基于任意拉格朗日-欧拉方法对射流破土过程进行仿真分析。以琼州海峡黏土质砂为研究对象,选用MAT-FHWA-SOIL模型建立土体材料本构。运用LS_DYNA软件建立包括射流源、喷嘴、水域及土体的仿真模型,探究射流破土机理。并讨论射流关键参数对破土效果的影响,对射流压力、喷嘴直径、射流靶距进行单因素研究。研究发现在其他射流参数不变的情况下,射流压力增大1倍,最大冲坑深度也增大1倍;土体的冲蚀体积会随着射流压力的增大而增加,与射流时间近似呈线性关系;喷嘴直径从5 cm增大到9 cm时,最大冲坑深度增幅为43%。在20~60 cm之间的射流靶距对清理效果的影响较小。该研究可为海底电缆抢修的射流破土作业过程提供依据。

钻井船插桩对邻近桩影响的耦合欧拉-拉格朗日有限元方法研究

依据模拟钻井船在黏土层中插桩对邻近桩影响的离心模型试验结果,研究了通过耦合欧拉-拉格朗日（CEL）有限元计算并结合非线性地基梁有限元计算,分析钻井船插桩对邻近桩影响问题的可行性。CEL有限元方法将产生大变形的土体设为欧拉体,采用欧拉有限元方法计算该区域的变形响应,计算过程中,欧拉体的空间网格形状、大小位置保持不变,物质可在网格之间运动;其他土体设为拉格朗日体,采用拉格朗日有限元方法计算变形响应,计算过程中,物质的运动和网格的变形保持一致。运用罚函数方法实现欧拉体与拉格朗日体的耦合。通过CEL有限元计算,可以确定钻井船插桩导致的邻近桩桩身水平位移。进一步通过非线性地基梁有限元模型计算确定桩身弯矩。计算结果表明,利用CEL有限元方法并结合非线性地基梁有限元方法计算出的桩身位移和弯矩沿桩长的变化与离心模型试验结果基本一致。说明采用CEL有限元方法并结合非线性地基梁有限元方法分析黏土层中插桩对邻近桩的影响问题是可行的。CEL有限元模型中欧拉土体范围的设置对计算结果有明显影响。研究表明,若插桩深度小于0.75倍桩靴直径,可将欧拉土体范围设置成1.00~1.25倍桩靴直径;若插桩深度大于0.75倍桩靴直径,将欧拉土体范围取为插桩深度以下0.5倍桩靴直径是恰当的。

基于耦合欧拉-拉格朗日方法的射弹高速入水尾拍数值分析

针对超空泡射弹火炮武器在航空反水雷作战中的应用,基于有限元分析软件Abaqus耦合欧拉-拉格朗日方法,开展射弹高速入水过程的数值模拟研究。通过对射弹垂直入水过程的数值模拟,获得高速入水载荷的一般形式;分析射弹倾斜入水过程中入水空泡的包覆状态与入水载荷、入水姿态变化之间的耦合关系,揭示射弹入水尾拍这一动态稳定运动的形成机理。不同入水攻角下射弹高速入水尾拍的仿真结果表明,在相同发射状态下,入水攻角是影响射弹密集度的重要因素。

基于代理拉格朗日松弛的电-氢耦合网络优化调度

针对电-氢耦合网络的调度问题,将电网与氢供应网协调优化,决定氢运输拖车的路径,并考虑氢运输的时效性,提出基于代理拉格朗日松弛的分布式求解方法。首先,考虑氢运输中拖车的行程和到达加氢站的时间等因素,以电网和氢供应网的运行成本最小为目标函数,建立电-氢耦合网络的混合整数线性规划调度模型;其次,根据代理拉格朗日松弛法将该调度问题分解成电制氢与氢供应子问题,并进行同步协调;最后,通过算例验证了所提模型及求解方法的有效性。

桶形基础液压下沉过程的耦合欧拉-拉格朗日有限元法分析

由于结构物贯入时网格变形过大而产生扭曲畸变等问题,常会造成岩土工程下沉贯入领域的数值分析收敛困难甚至计算结果失真。采用合适的数值方法分析此类问题颇具挑战性。传统的有限元模拟方法往往会出现收敛困难、作出不合理的假设以及需要依赖用户的专业网格重划分和插值程序等问题。耦合的欧拉-拉格朗日(CEL)分析方法结合了拉格朗日网格与欧拉网格的优点,可以有效地解决有关大变形和材料破坏等诸多问题。通过位移控制法和力控制法两种下沉方式,进行桶形基础室内液压下沉模型试验,得出不同强度的黏土中桶形基础下沉阻力和下沉深度的关系及土塞高度。应用CEL有限元法进行模拟,计算结果与试验结果较为符合。采用的CEL有限元模拟方法不仅可对桶形基础自重下沉和液压下沉进行预测,也可为其他海洋基础结构的贯入模拟提供有益参考。

基于欧拉-拉格朗日方法的水下气体泄漏扩散行为研究

针对水下气体泄漏扩散问题,基于计算流体动力学(CFD)理论采用流体体积模型(VOF)与离散相模型(DPM)耦合的方法,对气体在水中的扩散过程进行模拟与分析,水和空气作为连续相,泄漏气体作为离散相,离散相粒子与水下气泡具有相同的物理性质,其密度变化服从理想气体状态方程。基于建立的数值模型,研究水下气体羽流的形成和发展过程以及在水面形成的涌流效应,评估气体上浮时间、水面气池尺寸和涌流高度等参数。研究表明:泄漏气体以喷射状涌入水中,上升过程中体积逐渐膨胀增大,运动至水面时形成倒立的锥形羽流结构;气体带动表层水运动,引起羽流两侧表层水回流,在水面产生涌流效应和圆形的气池;涌流高度逐渐增大后呈小幅波动状发展,气池半径逐渐增大后稳定。仿真结果与小尺度实验数据对比验证了数值模型的可行性。

基于修正剑桥模型的软黏土中沉桩过程欧拉-拉格朗日耦合模拟分析

利用CEL(Coupled Eulerian-Lagrangian)有限元大变形方法模拟了不排水条件下软黏土中沉桩施工对桩周土体的影响.利用有效应力的方法实现了修正剑桥(MCC)模型Vumat子程序中超孔压的计算,通过与室内土工试验(一维压缩试验和三轴试验)结果对比证明了子程序的可靠性.将修正剑桥模型的Vumat子程序嵌入到Abaqus/Explicit的CEL方法中用于软黏土中沉桩效应的模拟分析.对计算得到的桩周土体流动、水平位移、超孔压、水平应力的分析表明,所采用的考虑孔压的CEL数值分析方法可以较好地模拟软黏土非排水条件下沉桩效应问题.

海冰动力学的混合拉格朗日-欧拉数值方法

综合考虑欧拉坐标下有限差分法(FDM)在海冰动力学计算中的效率,以及拉格朗日坐标下光滑质点流体动力学方法(SPH)对海冰流变行为的精确模拟,本文发展了一种海冰动力学的混合拉格朗日-欧拉(HLE)数值方法。该方法首先在拉格朗日坐标下将海冰离散为若干个具有厚度、密集度的海冰质点,并由这些海冰质点通过Gauss函数对欧拉网格上的海冰参量进行积分插值;然后,在欧拉坐标下对海冰动量方程进行差分计算以确定各网格节点的海冰速度,并由此采用Gauss函数积分插值出拉格朗日坐标下各海冰质点的速度分布;最后,通过对海冰质点运动和分布的计算,确定出各海冰质点的位置、厚度和密集度等参量。采用该HLE方法对规则区域内的海冰堆积过程和涡动风场作用下的海冰动力演化趋势进行了数值试验;最后,采用该HLE方法对渤海海冰的动力过程进行了72 h数值模拟,其计算结果与卫星遥感图像和现场观测资料吻合较好。以上计算结果均表明该HLE方法在海冰动力学数值模拟中具有较高的计算效率和模拟精度,可用于海冰动力过程的数值模拟。

基于耦合欧拉- 拉格朗日算法的航行体缓冲头帽冲击性能

航行体在入水过程中会受到强烈的冲击载荷,过高的冲击载荷极易造成结构破坏和内部设备失灵,研究结构高速入水降载技术具有十分重要的意义。设计一种以预制裂纹碳纤维材料作为外罩、泡沫铝作为内部缓冲材料的缓冲头帽。基于有限元方法,采用耦合欧拉- 拉格朗日求解器研究缓冲头帽的降载效果,通过对比数值结果与小球入水实验图像验证数值方法的有效性。分析在不同入水速度和入水角度条件下,缓冲头帽的抗冲击性能、损伤模式以及空泡的发展过程。结果表明:缓冲头帽撞击水面后沿着预制裂纹破碎,泡沫铝压缩吸能;随着自由液面以下的空泡不断向外扩张,破碎的头帽贴着空泡内壁向四周散开;缓冲头帽对于不同入水速度和入水角度的航行体有着良好的降载效果,在最恶劣的垂直入水条件下降载率可达75%以上。

基于欧拉-拉格朗日方法的携病毒飞沫扩散过程的数值模拟

新型冠状病毒肺炎短时间内在全球范围迅速传播,已经成为受国际关注的大流行病。携病毒飞沫的传播和扩散过程是新型冠状病毒传播的重要途径。文章通过多相流仿真手段,基于欧拉-拉格朗日方法,对房间、电梯、高铁车厢和飞机座舱等相对封闭环境中携病毒飞沫扩散过程进行了数据模拟研究。结果表明,相对封闭空间中采用对侧面进风和出风的通风形式,尤其是上下通风的方式,可有效减少飞沫在封闭空间中的扰动和停留时间,抑制携病毒飞沫的扩散,增强飞沫的排出。

基于耦合欧拉-拉格朗日方法的浮筒着水数值仿真

基于耦合欧拉-拉格朗日(CEL)方法对浮筒着水进行数值仿真。建立浮筒着水三维模型,在保持浮筒质量、水平速度不变的情况下,利用三维模型对不同初始姿态角和下沉速度时浮筒着水工况进行数值仿真,最后对比各工况下垂直加速度、垂直速度、姿态角时历曲线。结果表明:下沉速度越大,浮筒所受载荷越大,初始姿态角越小,浮筒发生第二次弹跳的时间就越早;垂直加速度的仿真结果与水池着水试验结果比较吻合。

一种基于任意拉格朗日-欧拉方法的降落伞充气展开数值模型

降落伞充气展开过程是一个典型时变非线性流构耦合问题,大多研究所建立的数值计算模型未考虑织物透气性和初始投放速度的影响。为修正模型,基于任意拉格朗日-欧拉(arbitrary Lagrange-Euler,ALE)方法和动网格技术,考虑织物透气性并预设初始投放速度,建立了降落伞充气展开的有限元模型。计算了充气过程中,有限质量情况下C-9伞的开伞力,并与实验数据对比,验证了该模型的正确性与合理性。随后分析了其开伞过程中伞衣应力分布特征与气动减速特性,为降落伞数值建模和结构设计提供一定参考。

下沉过程耦合欧拉-拉格朗日有限元法分析

在进行基础或结构物下沉贯入的数值分析时,常会遇到接触的处理困难、网格大变形造成的扭曲畸变等问题,引起程序收敛困难、计算失真。应用何种数值分析方法,准确的模拟下沉过程结构物和土的相互作用,已成为此领域的一个重要课题。耦合的欧拉-拉格朗日方法结合了欧拉计算方法与拉格朗日计算方法的特点,通过固定的欧拉网格及可以在网格中运动的材料,可有效的处理网格畸变等大变形分析中存在的解法困难,其接触算法可得到基础或结构物准确的应力应变状态。通过砂土中桩的静压下沉、吸力锚的自重下沉和海底管线自沉分析三个应用实例,表明CEL有限元法可有效模拟大变形问题,所得的计算结果有较好的适用性。其分析方法和成果可为其它类型基础结构的下沉分析提供有益参考。

基于分布鲁棒优化的分布式配-微协同日前调度方法

微网是消纳新能源和提高配电网可靠性的重要形态,通常有独立的能量管理系统。含多微网的配电网需要协同各个微网,满足安全经济运行的约束。在进行配-微协同调度的过程中,需要采用考虑新能源不确定性、离散调节设备特性的分布式协同方法。针对以上需求,采用分布鲁棒优化建模的思想对新能源的不确定性进行建模,并基于单纯分解方法-非线性分块高斯-赛德尔方法-增广拉格朗日方法(SDM-GS-ALM)对离散调节设备进行处理。该算法通过凸包对整数变量进行逼近,理论上保证了每次计算问题都是凸问题。在凸问题的基础上,将列与约束生成(C&CG)作为外层算法、SDM-GS-ALM作为内层算法内嵌在C&CG的各阶段,实现配-微协同分布式计算框架的构建。最后,基于北美240节点实际工程系统,采用线路运行状态的聚合等值方法验证所提方法的有效性。

利用欧拉-拉格朗日方法预测层流泡状流的含气率分布

利用欧拉拉格朗日方法,提出了用于预测竖直管道内绝热层流泡状流中含气率分布的三维模型。该模型能够跟踪单独的气泡轨迹,从而获取更多的界面力信息;同时气泡尺寸可以作为参数之一引入模型,使含气率分布的计算更为方便。模型中分析了绝热层流泡状流中气泡的各项受力表达式,建立了两种描述方法下的气液两相间的耦合关系。利用现有实验数据对模型进行的检验表明,该模型能够预测一定尺寸范围内气泡的分布;气泡径向分布主要取决于气泡所受侧向提升力。对于更大尺寸的气泡,气泡变形和气泡尾迹与当地流场间相互作用将对侧向提升力产生很大影响。

基于欧拉—拉格朗日模型的单相有源电力滤波器无源性控制新方法

有源电力滤波器(active power filter,APF)的控制方法是决定其补偿性能的关键因素之一。该文基于单相APF的Euler-Lagrange(EL)系统状态平均模型,提出了一种无源性控制新方法。该方法首先基于替代定理建立了考虑源阻抗和非线性负荷影响的单相APFEL系统状态平均模型,在此基础上设计了无源性间接控制律,确保对单相APF控制目标的渐近跟踪。由于在线计算直流电容指令电压波动量十分困难,在控制方法的执行中忽略波动量的影响,简化了控制方法的实现;此外,该方法还研究了所注入的阻尼大小对控制效果的影响以及APF内层控制过调制产生的原因,针对无源性控制律跟踪精度的要求与内层控制过调制限幅条件间的矛盾,通过构造模糊逻辑推理环节实现阻尼系数的在线调整,在满足内层控制限幅条件下确保跟踪精度和补偿效果。仿真结果验证了所提出方法的正确性和有效性。

基于耦合欧拉-拉格朗日方法的砂土中静压桩挤土效应数值模拟

为探讨砂土地基中静压桩沉桩过程的挤土效应,引入Simsand砂土临界状态本构模型,采用耦合欧拉-拉格朗日(coupled Eulerian-Lagrangian, CEL)方法,模拟分析了砂土中静压沉桩过程的挤土效应。基于相关试验数据校定本构模型参数,利用模拟离心机沉桩试验,验证数值模型的准确性。结果表明,挤土效应的影响区域直径随着桩尖角度和砂土相对密度的增大而增大,压入过程中,平底桩的挤土效应倾向于竖向向前,增大了土体阻力。桩基压入引起近地表土层发生隆起,密砂地层的隆起较松砂地层显著。随着深度增加,桩周土体的竖向位移由隆起变为向下移动,且趋于稳定,密砂中发生的侧向位移大于松砂地层。挤土效应的影响区域直径随桩端竖向向下方向迅速衰减,随桩端竖向向上方向发生小幅度减小,在密砂中沿径向方向呈现先增大后减小的趋势。

基于任意拉格朗日-欧拉法的盾构刀盘土体切削仿真

针对盾构刀盘切削土体的问题,提出了基于有限元和高性能计算方法的大规模数值模拟方法.建立了盾构和土体的三维分析模型,并且充分考虑到土体的材料非线性,利用动态显示算法进行了基于国产“曙光4000A”高性能计算机的大规模流固耦合分析.仿真采用任意拉格朗日-欧拉法(ALE),解决了土体大变形造成的有限元单元畸变问题,预测土体-刀具相互作用的过程;将修正的Drucker-Prager塑性土体模型作为土体的非线性有限元模型,其仿真方法和结果可以为未来盾构机的刀盘设计提供参考.

阵性泥石流运动与堆积的欧拉-拉格朗日模型──Ⅱ.应用

采用文献［１］中提出的阵性泥石流运动与堆积的欧拉－拉格朗日模型，模拟了阵性泥石流的运 动过程和堆积形态，得出了与地学研究中观测结果较为一致的认识．文中针对影响阵性泥石流的关 键条件，重点通过改变密度、坡度和各阵泥石流的时间间隔等参数，分析了这些参数变化对泥石流 运动及堆积规律的影响，提出了简化分析整个阵性泥石流的条件。研究表明，应用经试验资料验证 的数学模型不但可以方便、快捷地定量描述阵性泥石流的一般特性，而且能够提供关于阵性泥石流 运动及堆积的更多细节。