多级柱面炸药内爆磁通量压缩技术研究

柱面内爆磁通量压缩发生器是利用炸药内爆压缩其内部磁通量至轴线附近小体积内从而实现超高磁场,传统的单级装置因受到金属套筒内爆失稳等影响性能指标受限。开展了多级内爆磁压缩技术研究,突破多项关键技术,包括研制特殊结构的密绕螺线管、脉冲功率源及大电流放电开关等,具备在直径135 mm套筒空间内实现20 T以上初始磁场产生能力,并建立了动态磁光测量系统。利用磁流体力学编码SSS-MHD开展多级装置设计,计算显示,设计的多级装置能够将约42%的初始磁通量压缩至轴线附近直径7 mm的空间内。最终研制成功多级内爆磁压缩装置CJ-150,在亚立方厘米以上空间实现轴向峰值磁场强度906 T,数据不确定度5.35%。10余发动态考核实验显示,CJ-150装置工作稳定,能够满足物理实验需要。利用经实验验证的磁流体模型计算显示,CJ-150具备1 000 T以上超强磁场产生能力,能够对大尺寸样品实现500 GPa以上的准等熵加载。     

内爆加载金属界面不稳定性的数值分析

采用自行研制的多介质弹塑性流体力学欧拉程序,对柱形内爆加载金属界面不稳定性进行了数值研究,数值模拟结果与文献实验数据吻合较好。数值结果表明:材料强度对界面不稳定性发展有不可忽略的抑制作用;材料屈服强度对较高模数不稳定性增长的抑制较强,而剪切模量对不稳定性发展的影响相似但敏感性相对较弱;金属界面不稳定性增长存在最不稳定模数,最不稳定模数随屈服强度增加而减小,并近似与屈服强度的对数呈线性关系;随着壳的厚度减小,扰动增长加快。     

柱形和球形壳体内爆界面不稳定性数值分析

为了更好地研究柱形和球形构型下果冻界面不稳定性发展,避免内爆聚心反弹前后直角坐标网格计算导致的误差,提高对流场和界面位置的计算精度,通过应用考虑了MVFT程序的网格适应性,使其能够适用于柱形网格和球形网格下的界面不稳定性数值模拟,特别是能够保证内爆聚心反弹前后流场和界面计算的稳定性。应用改进的计算程序对两种构型下的界面不稳定性进行了数值模拟,并对二者界面演化规律进行了详细讨论和归纳。结果表明:对于内外半径相同的柱形和球形果冻,聚心反弹时前者半径较小,而后者反弹时刻早于前者,其向内聚心和向外运动的速度最大值大于前者,对内部气体的压缩强度强于前者。对于外边界带有正弦扰动情况,除遵循上述规律外,计算还给出了峰谷转换现象。该项研究结果为进一步深入进行复杂构型下界面不稳定性高精度数值模拟研究提供了一种分析工具。     

界面不稳定性引起混合过程的二维数值计算

在可压缩多介质流体动力学高精度欧拉计算方法多介质流体分段抛物方法(multi-fluid piecewise parabolic method, MFPPM)基础上,运用算子分裂技术,增加二阶空间中心差方法和两步Rung-Kutta时间推进方法计算动力学黏性以及热流部分对流场的影响,发展适用于NS(Navier-Stokes)方程的可压缩多介质黏性流体计算方法多介质黏性流体分段抛物方法(multi-viscousity-fluid piecewise parabolic method, MVPPM). 文中采用MVPPM对英国AWE(atomic weapons establishment)激波管实验进行二维计算,给出了与实验图像基本一致的计算结果;应用MFPPM和MVPPM分别对二维柱对称内爆动力学界面不稳定性及其后期混合过程进行数值模拟,给出内外界面演化、速度历史以及后期中心气穴不同半径内因RT(Rayleigh-Taylor)界面不稳定性引起的混合量分布情况,从计算结果比较可见黏性对物质界面处混合量的分布影响明显.      

不同初始磁场对激波冲击R22重气柱过程影响的数值模拟

基于磁流体动力(magneto-hydrodynamic,MHD)方程,采用CTU+CT方法,对在不同初始磁场作用下的平面入射激波与磁化R22重质气柱作用过程进行了数值研究。数值结果清晰地描述了不同初始磁场条件时激波诱导R22气柱界面不稳定性的过程,揭示了磁场控制界面不稳定性的机理。另外,还分析了磁感应强度对界面不稳定性的影响,发现在磁场较小时,涡层附着于界面,但随着磁感应强度的增大,平均涡量随之增大,涡层与界面逐渐分离,最终更好地抑制了界面不稳定性。同时,还发现平均涡度拟能随着磁感应强度的增大而减小,而垂直磁场比平行磁场更能降低平均涡度拟能,因而平均涡度拟能可较好地反映磁场对不稳定性的影响效果。     

内爆炸载荷下泡沫铝夹芯圆管塑性动力响应及能量耗散机理

通过实验、理论和数值模拟研究泡沫铝夹芯空心圆管在内爆载荷作用下的动态变形模式及吸能机制.采用不同质量的球形乳化炸药进行爆炸试验.结构轴向变形分为3个区:大塑性变形区、绕塑性铰的刚性旋转区和无变形区.在考虑环向膜力和轴向弯矩的情况下,提出内爆作用下夹芯圆管动态响应的显式计算方法.通过建立基于三维Voronoi算法的泡沫芯有限元模型,探索结构能量耗散机制.通过实验观测到的泡沫铝夹芯空心圆管在内爆载荷作用下的变形机制,结合内外管的弯曲变形及芯层压缩,给出了结构在响应过程中能量吸收的理论解.以结构比吸能和外管中心挠度为控制参量求得夹芯圆管的最优解集.进一步研究炸药质量、内外管直径、壁厚及芯层轴向梯度排列方式对结构动态变形模式和吸能机制的影响.结果表明:内管壁厚对外管中心线的挠度影响较大,而芯层厚度和外管壁厚的影响较小;如果夹芯圆管的轴向梯度结构从管轴向对称面到两端边缘呈对称递减分布时,具有较好的抗爆性;数值计算和实验测试结果均与理论预测吻合.     

柱形汇聚几何中内爆驱动金属界面不稳定性

采用自研的高保真度爆轰与冲击动力学程序,对柱形汇聚几何中内爆驱动金属材料界面不稳定性的动力学行为,进行了数值模拟研究。结果表明,首次冲击后至约12 μs,界面发展以RM（Richtmyer-Meshkov）不稳定性为主;12 μs后至冲击波聚心反弹加载前,界面聚心运动处于加速减速状态,界面发展由RT (Rayleigh-Taylor)不稳定性主导;冲击波聚心反弹加载后,界面发展又由RM不稳定性主导。另外,还研究了初始条件（初始振幅、初始波长、钢壳初始厚度和几何构型）对柱形内爆驱动金属材料界面不稳定性的影响。结果显示:初始振幅较大时振幅增长也较大;初始波长较小（模数较大）时振幅增长较小,而且存在一个截止波长;钢壳厚度会抑制扰动增长,也存在一个截止厚度;几何汇聚效应会使扰动增长速度更快。     

可压缩多介质粘性流体的数值计算

将考虑热传导和粘性情况下的Navier Stokes方程描述的物理过程分解成3个子过程进行数值计算,即把整个流量计算分解成无粘性流量、粘性流量和热流量3部分,采用多介质流体高精度parabolic piecewise method（PPM）方法、二阶空间中心差方法和两步Rung-Kutta时间推进方法相结合进行数值计算。给出了激波管中Riemann问题和二维、三维Richtmyer-Meshkov界面不稳定性的Navier Stokes方程和Euler方程对比计算结果,显示了粘性对界面不稳定性的影响。     

电磁加载下的高能量密度物理问题研究（续1）

3电磁加载高能量密度物理研究进展3.1Z箍缩内爆ICF研究3.1.1Z箍缩内爆等离子体负载的构形—真空黑腔和动力黑腔利用Z箍缩内爆等离子体进行惯性约束聚变(ICF)研究,最初来自1973年Turchi和Bak-er[9]提出的电磁内爆概念,这里包含着两条不同的途径:一是类似于武器的概念,利用高速内爆的固体套筒压缩其内部的热核材料达到高温、高密度状态,发生热核点火。例如俄罗斯实验物理院利用爆炸磁压缩强磁场发生器(MC-1)进行的工作,以及美国Los Alamos实验室的Athena计划等。这条途径的困难在于内爆后期固体套筒流体动力学不稳定性发展严重,很难把…     

果冻内爆实验研究

利用最近研制的果冻内爆装置开展了果冻环在内爆加载条件下Ray—leigh—Taylor不稳定性实验研究。实验结果表明,高速摄影测试图像的效果有明显提高,能够初步观测果冻环反弹阶段内界面的扰动发展情况,获得了果冻环在加速压缩、减速压缩以及反弹过程中内外界面扰动的变化过程以及混合区宽度的发展历程。     

ATOMIZATION OF A LIQUID DROP BY PULSATION

ＡＴＯＭＩＺＡＴＩＯＮＯＦＡＬＩＱＵＩＤＤＲＯＰＢＹＰＵＬＳＡＴＩＯＮ￥(林松飘，周哲玮)Ｓ．Ｐ．Ｌｉｎ；（ＤｅｐａｒｔｍｅｎｔｏｆＭｅｃｈａｎｉｃａｌａｎｄＡｅｒｏｎａｕｔｉｃａｌＥｎｇｉｎｅｅｒｉｎｇＣｌａｒｋｓｏｎＵｎｉｖｅｒｓｉｔｙ，Ｐｏｔｓｄ...     

Computation of transition to turbulence in the compression stage of a reciprocating engine

The transition to turbulent flow in the compression stage of a reciprocating engine is studied by obtaining the finite-difference numerical solutions to the governing Navier-Stokes equations without using explicit turbulence models. A computational method is developed under the assumption that the flow is in a low-subsonic regime with strong compression. The numerical method is a simple extension of the well known MAC method. Computations were performed for three different chamber geometries at the engine speed of 1400 rpm. The results of the computations clearly demonstrate the transient process in which large tumbling vortices break down into smaller ones near the end of the compression process. The transition process is also caught experimentally by using Mach-Zehnder interferometry.     

STABILITY OF ULTRATHIN LIQUID FILM EVOLUTION WITH SURFACTANT

针对固体基底上厚度小于100 nm的含活性剂超薄液膜演化过程, 基于润滑理论推导出包含分离压影响的液膜厚度和活性剂浓度的演化方程, 采用正则模态法导出了描述液膜线性稳定性的特征方程, 分析了多个特征参数对线性稳定性的影响, 数值模拟了液膜厚度和活性剂浓度演化历程, 对比了模拟所得非线性结果与线性分析预测结果的一致性.结果表明:范德华力具有促进扰动增长的作用, 较强的玻恩斥力促使扰动衰减, 使液膜趋于稳定;较小的毛细力数易使液膜凹陷处发生二次失稳, 并最终导致去润湿现象发生;液膜厚度和溶于液膜内部的活性剂浓度初值越大, 液膜稳定性越强, 液膜表面活性剂浓度影响则相反;增大吸附系数不利于液膜稳定性.     

纵向冲击压缩下LY12铝合金圆环的塑性失稳

通过对铝合金圆环的纵向冲击压缩研究发现,一定条件下在试件的宏观塑性硬化阶段会出现明显的应力降过程。为揭示此应力降的发生机制,对润滑、细磨、粗磨3种端面粗糙条件下,外径、内径和高度比值为6:3:2的LY12铝合金圆环进行系统的Hopkinson压杆纵向冲击实验。结果表明:应力降主要发生在较大的应变和较高的应变率条件。进一步对实验样品的金相观察发现:应力降产生的内在机制为绝热剪切带的形成和发展,此现象是一种动态塑性失稳的过程。以上结果为金属材料在冲击条件下绝热剪切带产生的研究提供了参考。     

不同温度下2种钨结构增强Zr基非晶合金的动态压缩性能

利用SHPB装置对钨丝增强Zr基非晶复合材料和钨骨架增强Zr基非晶复合材料进行了3种环境温度下多种应变率的动态压缩性能测试。比较了2种材料的动态力学性能,发现二者均具有应变率敏感性和较强的塑性变形能力。但二者承载机制存在较大差异。钨丝增强结构变形主要表现在钨丝的失稳,由数值模拟初步分析了这种局部结构失稳控制的变形以及热失稳现象;钨骨架增强结构变形前期钨骨架起主要承载作用,而不是各成分的共同作用,这导致材料的屈服强度比纯非晶和纯钨的低。     

Theoretical analysis for the mechanical behavior caused by an electromagnetic cycle in ITER Nb_3Sn cable-in-conduit conductors

The central solenoid(CS)is one of the key components of the International Thermonuclear Experimental Reactor(ITER)tokamak and which is often considered as the heart of this fusion reactor.This solenoid will be built by using Nb3Sn cablein-conduit conductors(CICC),capable of generating a 13 T magnetic field.In order to assess the performance of the Nb3Sn CICC in nearly the ITER condition,many short samples have been evaluated at the SULTAN test facility(the background magnetic field is of 10.85 T with the uniform length of 400 mm at 1%homogeneity)in Centre de Recherches en Physique des Plasma(CRPP).It is found that the samples with pseudo-long twist pitch(including baseline specimens)show a significant degradation in the current-sharing temperature(Tcs),while the qualification tests of all short twist pitch(STP)samples,which show no degradation versus electromagnetic cycling,even exhibits an increase of Tcs.This behavior was perfectly reproduced in the coil experiments at the central solenoid model coil(CSMC)facility last year.In this paper,the complex structure of the Nb3Sn CICC would be simplified into a wire rope consisting of six petals and a cooling spiral.An analytical formula for the Tcs behavior as a function of the axial strain of the cable is presented.Based on this,the effects of twist pitch,axial and transverse stiffness,thermal mismatch,cycling number,magnetic distribution,etc.,on the axial strain are discussed systematically.The calculated Tcs behavior with cycle number show consistency with the previous experimental results qualitatively and quantitatively.Lastly,we focus on the relationship between Tcs and axial strain of the cable,and we conclude that the Tcs behavior caused by electromagnetic cycles is determined by the cable axial strain.Once the cable is in a compression situation,this compression strain and its accumulation would lead to the Tcs degradation.The experimental observation of the Tcs enhancement in the CS STP samples should be considered as a contribution of the shorter length of the high field zone in SULTAN and CSMC devices,as well as the tight cable structure.     

端部切向载荷作用下的压杆失稳分析及其应用

细直杆件在压应力作用下会产生横向屈曲即失稳．直杆撞击刚性平面或拉断卸载后将形成压缩波，因承载压缩载荷的长度增加可以引起失稳．冲击速度转换的压应力沿着杆件切线方向，该处弯矩和剪力为零；而众多文献设定的失稳段固支边界条件并不准确．基于精确的杆件变形曲率方程得到端部载荷指向杆件中固定点时的受压失稳条件，得到其极限状态即载荷沿杆端切向作用时失稳长度相当于两端简支的1.5 倍．对于钢丝绳拉断形成的冲击失稳，载荷恒定而长度增加，可以产生高阶屈曲即在侧向出现多次曲折，并基于尼龙-橡胶带的模拟试验给出了定性说明．     

Suppression or enhancement of interfacial instability in two-layer plane Couette flow of FENE-P fluids past a deformable solid layer

The linear stability of two-layer plane Couette flow of FENE-P fluids past a deformable solid layer is analyzed in order to examine the effect of solid deformability on the interfacial instability due to elasticity and viscosity stratification at the two-fluid interface. The solid layer is modeled using both linear viscoelastic and neo-Hookean constitutive equations. The limiting case of two-layer flow of upper-convected Maxwell (UCM) fluids is used as a starting point, and results for the FENE-P case are obtained by numerically continuing the UCM results for the interfacial mode to finite values of the chain extensibility parameter. For the case of two-layer plane Couette flow past a rigid solid surface, our results show that the finite extensibility of the polymer chain significantly alters the neutral stability boundaries of the interfacial instability. In particular, the two-layer Couette flow of FENE-P fluids is found to be unstable in a larger range of nondimensional parameters when compared to two-layer flow of UCM fluids. The presence of the deformable solid layer is shown to completely suppress the interfacial instability in most of the parameter regimes where the interfacial mode is unstable, while it could have a completely destabilizing effect in other parameter regimes even when the interfacial mode is stable in rigid channels. When compared with two-layer UCM flow, the two-layer FENE-P case is found in general to require solid layers with relatively lower shear modulii in order to suppress the interfacial instability. The results from the linear elastic solid model are compared with those obtained using the (more rigorous) neo-Hookean model for the solid, and good agreement is found between the two models for neutral stability curves pertaining to the two-fluid interfacial mode. The present study thus provides an important extension of the earlier analysis of two-layer UCM flow [V. Shankar, Stability of two-layer viscoelastic plane Couette flow past a deformable solid layer: implications of fluid viscosity stratification, J. Non-Newtonian Fluid Mech. 125 (2005) 143–158] to more accurate constitutive models for the fluid and solid layers, and reaffirms the central conclusion of instability suppression in two-layer flows of viscoelastic fluids by soft elastomeric coatings in more realistic settings.     

A strain localization analysis using a viscoplastic softening model for clay

Strain localization has become an attractive subject in geomechanics during the past decade. Shear bands are well known to develop in clay specimens during the straining process. Strain localization is closely related to plastic instability. In the present paper, a non-linear instability condition for the viscoplastic strain softening model during the creep process is firstly obtained. It is found that the proposed viscoplastic model is capable of describing plastic instability. Secondly, a two-dimensional linear instability analysis is performed and the preferred orientation for the growth of fluctuation and the instability condition are derived. It is worth noting that the two instability conditions are equivalent. Finally, the behavior of the clay is numerically analyzed in undrained plane-strain compression tests by the finite element method, considering a transport of pore water in the material at a quasi-static strain rate. The numerical results show that the model can predict strain localization phenomena, such as shear banding. From the numerical calculations, the effects of strain rate and permeability are discussed.     

Investigation of the effect of SSP in stabilizing the structure of condensation polymer blends via rheological measurements

We investigated the influence of solid-state polymerization (SSP) process on the reactions that could be taken place at the interphase of polycondensation polymer blends to stabilize the structure obtained after melt mixing. Polyethylene terephthalate (PET) and polyamide 66 (PA66) were melt blended in a mixer, and subsequent SSP process was performed for each sample. FTIR spectra indicated reactions between two polymers. Viscosity behavior and interfacial slip were investigated by measuring shear viscosity of components and blends before and after SSP and then compared with the viscosity calculated from the log-additivity model. The results showed that after SSP, there was no sign of interfacial slip, the slope of viscosity reduction with increasing shear rate became smaller, and the viscosity of blends showed positive deviation at all examined shear rate from the log-additivity model while this deviation was negative at higher shear rate before SSP. SEM micrographs, which were taken after shear stress was imposed on the samples, also indicated the morphological stability after SSP. Furthermore, we studied the effect of functional groups concentration on the reactions at the interphase by using hydrolyzed PET as a precursor for blends. The results showed that slip at the interface would decrease with increasing functional groups of the precursors. These results are particularly valuable for using recycled polymers.