二维GEL耦合方法的研究及对爆炸容器的数值模拟

 将以Euler方法为基础的MF PPM（Piecewise-Parabolic Method）程序和以Lagrange方法为基础的DEFEL（2-D Finite Elements Code,二维流体弹塑性动力有限元）程序,根据压力和法向速度连续准则进行耦合,发展了基于Level Set的GEL（Ghost-Fluid Euler-Lagrange）方法。该方法在处理大变形流场与小变形结构以及复杂流动与多物体相互作用等问题具有优越性。通过二维算例的计算结果与文献比较,检验了GEL方法和耦合程序的正确性,并对球形和椭球封头的爆炸容器进行了数值模拟,通过与实验结果的比较分析,表明本研究程序可以比较好地处理内爆引起的壳体流固耦合问题。     

基于IB-LBFS和绝对节点坐标法的降落伞柔性结构流固耦合数值模拟

在空气动力学、水动力学和生物流体力学领域中,大变形柔性结构的流固耦合现象是一个重要的非线性力学问题.对该系统的数值模拟是分析这一问题的有效手段.将近年提出的一种Descartes流场求解器,即浸润边界-格子Boltzmann通量求解器(immersed boundary-lattice Boltzmann flux solver,IB-LBFS)作为流场求解方法,并引入绝对节点坐标法(absolute nodal coordinate formulation,ANCF)作为大变形结构分析手段,构建了流固耦合求解器以模拟三维流场中的大变形柔性体运动.使用三维来流中的旗帜摆动算例对该流固耦合求解器进行了验证计算.基于该流固耦合求解器对三维不可压流场中的矩形降落伞和十字形降落伞的展开过程进行了非定常流固耦合数值模拟.     

基于气体动理学方法的高超声速支杆降热模拟

基于气体动理学方法(gas kinetic scheme, GKS)分别从单独流场和流固耦合两方面开展了针对支杆结构热防护系统的数值模拟研究。首先利用二维标准圆管模型验证算法在单独流场计算中的可靠性,然后通过比较不带支杆和带支杆情形下圆管外壁面的气动热特性及相应流场特征,分析得到支杆降热的有效性及其降热机理。接着在将GKS应用到流场/结构温度场耦合模拟并进行验证后,进一步分析流固耦合作用下的支杆降热性能,结果发现支杆的降热效果随着流固耦合时间的推进不断增强,从而有利于高超声速飞行器长时间巡航飞行。     

抗爆容器内爆炸流场数值模拟

 采用计算流体动力学中的二阶精度TVD差分格式和特殊算子分裂法,按轴对称问题,对半球顶圆柱筒密闭式抗爆容器内部爆炸流场进行了数值模拟。计算得到的容器壁面载荷分布与实验结果基本一致。不同时刻爆炸流场压力分布图像清晰地描述了容器壁面的冲击波加载过程。     

高超声速飞行器前缘流固耦合计算方法研究

对高超声速流场和结构温度场进行了耦合计算分析, 同时基于准静态假设对结构应力进行了分析. 流场部分采用基于非定常Navier-Stokes (N-S)方程的有限体积法, 湍流模型采用SST k-ω 模型, 固体部分采用基于非稳态热传导方程的有限元法, 同时基于准静态假设对固体结构的应力应变进行了分析. 在流固交界面处, 高速流体从固体结构得到温度边界条件, 固体结构从高速流体得到热流边界条件, 从而实现了流场和固体温度场的紧耦合计算.通过与超声速无限长圆管绕流试验结果进行对比, 验证了该方法的可靠性. 同时对二维圆管结构在气动加热过程中的温度、应力等的变化进行了比较详细的分析. 研究结果表明: 随着气动加热时间的推进, 由于圆管结构的高温区在不断扩大, 导致了结构的热变形在不断地增大; 圆管最小变形区出现在θ为60°处; 同时研究发现在计算时间内圆管热变形对外部流场的影响可以忽略不计.     

超声波振动辅助下填充填缝工艺的数值研究

本文基于数值方法研究了超声波振动辅助倒装芯片成型下填充填缝工艺。采用了一种基于紧耦合的流固耦合算法,对结构和流场之间的耦合运动进行了数值模拟,其中流场部分采用有限体积法对任意拉格朗日-欧拉方法的不可压缩N-S方程进行离散计算,而结构部分采用有限元法对拉格朗日坐标下的弹性动力学方程进行离散计算。在一个时间步内,结构与流场的计算区域的交界面上进行了多次的数值传递和插值,以保证满足耦合面边界条件。计算了超声波振动作为边界条件作用于结构场,而引起的流场中流体流速的变化规律,并通过实验验证了本方法的正确性。同时,监测了流场入口体积流量的变化规律,探究了超声波振幅、频率以及流体粘度等因素对流体流速的影响,为超声波振动辅助倒装芯片成型下填充工艺的应用提供了理论依据。     

内部爆炸载荷作用下混凝土动力响应的数值模拟

 目前,采用实验方法对爆炸载荷破坏效应进行全面、准确描述仍是一个难题。运用ANSYS/LS-DYNA非线性显式动力学有限元程序,采用流固耦合算法,对混凝土在内部爆炸载荷作用下的空腔形成和发展规律等动力响应问题进行了数值模拟。通过实际算例,比较研究了混凝土敷设钢板与否对成坑大小、损伤区、压力云图以及应力时程曲线等的影响,进一步研究了有钢板加固时混凝土在爆炸荷载作用下的破坏机理。研究表明：采用流固耦合算法可以避免由于物质的大变形引起的单元畸变和网格消蚀等问题,提高了计算效率;很好地模拟了混凝土内部爆炸时,中心爆炸空腔压缩面和自由面的拉伸层裂破坏过程,且与理论计算值吻合较好;通过表面敷设钢板,可以减小混凝土的拉应力峰值,抑制混凝土中裂纹的发展,从而提高其防震塌能力以及结构的完整性。     

流固耦合声子晶体管路冲击振动特性研究

流固耦合管路系统广泛应用于各种装备中,通常用来传递物质和能量或者动量.由于流固耦合效应,管壁在流体作用下易产生强烈的振动与噪声,对装备安全性、隐蔽性产生严重影响,甚至造成严重破坏.流固耦合管路振动抑制需求迫切,意义重大.声子晶体可以利用其带隙特性抑制特定频率范围内弹性波的传播,在减振降噪领域具有广泛的应用前景.本文基于声子晶体理论,研究了流固耦合条件下的布拉格声子晶体管路冲击振动传递特性.将传递矩阵法和有限元法相结合,计算了能带结构与带隙特性,重点考虑了流固耦合效应下,不同冲击激励条件下声子晶体管路振动特性,分析了流固耦合对声子晶体管路振动传递特性的影响.研究结果为流固耦合条件下管路系统的振动控制提供了技术参考.     

基于浸入边界-多松弛时间格子玻尔兹曼通量求解法的流固耦合算法研究

针对流固耦合问题,发展了基于浸入边界-多松弛时间格子玻尔兹曼通量求解法(immersed boundary method multi-relaxation-time lattice Boltzmann flux solver,IB-MRT-LBFS)的弱耦合算法.依据多尺度Chapman-Enskog展开,建立不可压宏观方程状态变量和通量与格子玻尔兹曼方程中粒子密度分布函数之间的关系;采用强制浸入边界法处理流固界面使固壁表面满足无滑移边界条件,根据修正的速度求解动量方程力源项;结构运动方程采用四阶龙格-库塔法求解.格子模型与浸入边界法的引入使流固耦合计算可以在笛卡尔网格下进行,无需生成贴体网格及运用动网格技术,简化了计算过程.数值模拟了单圆柱横向涡激振动、单圆柱及串列双圆柱双自由度涡激振动问题.结果表明,IB-MRT-LBFS能够准确预测圆柱涡激振动的锁定区间、振动响应、受力情况以及捕捉尾流场结构形态,验证了该算法在求解流固耦合问题的有效性和可行性.     

轴流压气机转子叶片的流固耦合分析

本文针对压气机转子叶片流固耦合问题,采用耦合平台MpCCI连接商业CFD软件FINE/Turbo与CSM软件ABAQUS模拟了设计转速下不同流量工况点的流固耦合问题,对比了流场定常计算与流固耦合非定常计算的结果,得到了叶片振动位移频谱振幅图以及效率和压比变化曲线。研究结果表明在失速点附近叶片振动对流场有显著影响,由此可以对压气机的颤振问题进行预测。     

电爆炸箔驱动绝缘飞片的一维数值模拟

 在一维有限差分反应流体动力学程序SSS中引入电爆炸模型,在此基础上结合新引入的物态方程库,开展了电爆炸加载下飞片运动过程的数值模拟研究。计算结果与实验结果基本吻合,表明SSS对电爆炸驱动飞片的过程的模拟是可靠的,为电炮加载技术的优化设计和电爆炸驱动飞片物理过程的进一步研究提供了一套有效的计算方法和程序。     

Numerical Simulation of 2-D Radiation-Drive Ignition Implosion Process

A Lagrangian compatible radiation hydrodynamic algorithm and the nuclear dynamics computing module are developed and implemented in the LARED Integration code, which is a radiation hydrodynamic code based on the 2-D cylindrical coordinates for the numerical simulation of the indirect-drive Inertial Confined Fusion. A number of 1-D and 2-D ignition implosion numerical simulations by using the improved LARED Integration code (ILARED) are presented which show that the 1-D numerical results are consistent with those computed by the 1-D radiation hydrodynamic code RDMG, while the simulation results of the 2-D low-mode radiative asymmetry and hydrodynamic instability growth,according to the physical analysis and anticipation, are satisfactory. The capsules driven by the sources from SGII experiments are also simulated by ILARED, and the fuel shapes agree well with the experimental results. The numerical simulations demonstrate that ILARED can be used in the simulation of the 1-D and 2-D ignition capsule implosion using the multi-group diffusion model for radiation.     

The numerical study of shock-induced hydrodynamic instability and mixing

Based on multi-fluid volume fraction and piecewise parabolic method (PPM), a multi-viscosity-fluid hydrodynamic code MVPPM (Multi-Viscosity-Fluid Piecewise Parabolic Method) is developed and applied to the problems of shock-induced hydrodynamic interfacial instability and mixing. Simulations of gas/liquid interface instability show that the influences of initial perturbations on the fluid mixing zone (FMZ) growth are significant, especially at the late stages, while grids have only a slight effect on the FMZ width, when the interface is impulsively accelerated by a shock wave passing through it. A numerical study of the hydrodynamic interfacial instability and mixing of gaseous flows impacted by re-shocks is presented. It reveals that the numerical results are in good agreement with the experimental results and the mixing growth rate strongly depends on initial conditions. Ultimately, the jelly layer experiment relevant to the instability impacted by exploding is simulated. The shape of jelly interface, position of front face of jelly layer, crest and trough of perturbation versus time are given; their simulated results are in good agreement with experimental results.     

Fully implicit 1D radiation hydrodynamics: Validation and verification

A fully implicit finite difference scheme has been developed to solve the hydrodynamic equations coupled with radiation transport. Solution of the time-dependent radiation transport equation is obtained using the discrete ordinates method and the energy flow into the Lagrangian meshes as a result of radiation interaction is fully accounted for. A tridiagonal matrix system is solved at each time step to determine the hydrodynamic variables implicitly. The results obtained from this fully implicit radiation hydrodynamics code in the planar geometry agrees well with the scaling law for radiation driven strong shock propagation in aluminium. For the point explosion problem the self similar solutions are compared with results for pure hydrodynamic case in spherical geometry. Results obtained when radiation interaction is also accounted agree with those of point explosion with heat conduction for lower input energies. Having, thus, benchmarked the code, self convergence of the method w.r.t. time step is studied in detail for both the planar and spherical problems. Spatial as well as temporal convergence rates are ?1 as expected from the difference forms of mass, momentum and energy conservation equations. This shows that the asymptotic convergence rate of the code is realized properly.     

Numerical modeling of the exterior-to-interior transmission of impulsive sound through three-dimensional,thin-walled elastic structures

Exterior propagation of impulsive sound and its transmission through three-dimensional, thin-walled elastic structures, into enclosed cavities, are investigated numerically in the framework of linear dynamics. A model was developed in the time domain by combining two numerical tools: (i) exterior sound propagation and induced structural loading are computed using the image-source method for the reflected field (specular reflections) combined with an extension of the Biot–Tolstoy–Medwin method for the diffracted field, (ii) the fully coupled vibro-acoustic response of the interior fluid–structure system is computed using a truncated modal-decomposition approach. In the model for exterior sound propagation, it is assumed that all surfaces are acoustically rigid. Since coupling between the structure and the exterior fluid is not enforced, the model is applicable to the case of a light exterior fluid and arbitrary interior fluid(s). The structural modes are computed with the finite-element method using shell elements. Acoustic modes are computed analytically assuming acoustically rigid boundaries and rectangular geometries of the enclosed cavities. This model is verified against finite-element solutions for the cases of rectangular structures containing one and two cavities, respectively.     

Influence of wall motion on particle sedimentation using hybrid LB-IBM scheme

We integrate the lattice Boltzmann method (LBM) and immersed boundary method (IBM) to capture the coupling between a rigid boundary surface and the hydrodynamic response of an enclosed particle laden fluid. We focus on a rigid box filled with a Newtonian fluid where the drag force based on the slip velocity at the wall and settling particles induces the interaction. We impose an external harmonic oscillation on the system boundary and found interesting results in the sedimentation behavior. Our results reveal that the sedimentation and particle locations are sensitive to the boundary walls oscillation amplitude and the subsequent changes on the enclosed flow field. Two different particle distribution analyses were performed and showed the presence of an agglomerate structure of particles. Despite the increase in the amplitude of wall motion, the turbulence level of the flow field and distribution of particles are found to be less in quantity compared to the stationary walls. The integrated LBM-IBM methodology promised the prospect of an efficient and accurate dynamic coupling between a non-compliant bounding surface and flow field in a wide-range of systems. Understanding the dynamics of the fluid-filled box can be particularly important in a simulation of particle deposition within biological systems and other engineering applications.     

Order and fluctuations in nonequilibrium molecular dynamics simulations of two-dimensional fluids

Finite systems of hard disks placed in a temperature gradient and in an external constant field have been studied, simulating a fluid heated from below. We used the methods of nonequilibrium molecular dynamics. The goal was to observe the onset of convection in the fluid. Systems of more than 5000 particles have been considered and the choice of parameters has been made in order to have a Rayleigh number larger than the critical one calculated from the hydrodynamic equations. The appearance of rolls and the large fluctuations in the velocity field are the main features of these simulations.     

Existence of a Hartmann layer in the peristalsis of Sisko fluid

Analytical solutions for the peristaltic flow of a magneto hydrodynamic(MHD) Sisko fluid in a channel, under the effects of strong and weak magnetic fields, are presented. The governing nonlinear problem, for the strong magnetic field,is solved using the matched asymptotic expansion. The solution for the weak magnetic field is obtained using a regular perturbation method. The main observation is the existence of a Hartman boundary layer for the strong magnetic field at the location of the two plates of the channel. The thickness of the Hartmann boundary layer is determined analytically. The effects of a strong magnetic field and the shear thinning parameter of the Sisko fluid on the velocity profile are presented graphically.     

Diffusion of colloidal particles in swimming suspensions

Previously, we have proposed to analyse the hydrodynamic interactions in a suspension of swimmers with respect to an effective hydrodynamic diffusion coefficient, which only considers the fluctuating motion caused by the stirring of the fluid. In this work, we study the diffusion of colloidal particles immersed in a bath of swimmers. To accurately resolve the many-body hydrodynamic interactions responsible for this diffusion, we use a direct numerical simulation scheme based on the smooth profile method. We consider a squirmer model for the self-propelled swimmers, as it accurately reproduces the flow field generated by real microorganisms, such as bacteria or spermatozoa. We show that the diffusion coefficients of the colloids are comparable with the effective diffusion coefficients of the swimmers, provided that the concentration of swimmers is high enough. At low concentrations, the difference in the way colloids and swimmers react to the flow leads to a reduction in the diffusion coefficient of the colloids. This is clearly seen in the appearance of a negative-correlation region for the velocity-correlation function of the colloids, which does not exist for the swimmers.