亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        Numerical simulation of two-phase flow field in underwater sealing device based on dynamic mesh

        2015-03-03 08:37:16ZHANGXueweiLIQiangLVMengrou
        關(guān)鍵詞:中北大學(xué)激波工程學(xué)院

        ZHANG Xue-wei, LI Qiang, LV Meng-rou

        (College of Mechatronic Engineering, North University of China, Taiyuan 030051, China)

        張學(xué)偉, 李 強(qiáng), 呂夢柔

        (中北大學(xué) 機(jī)電工程學(xué)院, 山西 太原 030051)

        ?

        Numerical simulation of two-phase flow field in underwater sealing device based on dynamic mesh

        ZHANG Xue-wei, LI Qiang, LV Meng-rou

        (CollegeofMechatronicEngineering,NorthUniversityofChina,Taiyuan030051,China)

        In order to speed underwater launch of minor-caliber weapons, a sealing device can be set in front of underwater muzzle to separate water, preventing the muzzle from water immersion. By establishing and simplifying the model of underwater weapon sealing device and unstructured mesh computing domain model based on computational fluid dynamics (CFD), dynamic mesh and user defined function (UDF), the N-S equation is solved and the numerical analysis and calculation of the complex two-phase flow inside the sealing device are carried out. The results show that the gas discharged from the sealing device is conducive to the formation of the projectile supercavity. When the projectile is launched at 5 m under water, the shock wave before and after the projectile has impact on the box body up to 100 MPa, therefore the sealing device must be strong enough. The research results have the vital significance to the design of underwater weapon sealing device and the formation of the projectile supercavitation.

        two-phase flow; supercavitation; sealing device; computational fluid dynamics (CFD); dynamic mech

        As a result of the change of international situation, the water weapon has got rapid development. The density of water is about 850 times as dense as that of air, thus water resistance is too large for projectile weapons to hit the target at a distance. Therefore, how to increase the initial velocity and decrease the resistance of underwater weapon inevitably becomes the key of the research.

        In recent years, many experts and scholars have done a lot of work in the numerical simulation of underwater jet and supercavity based on computational fluid dynamics (CFD). TANG, et al.[1]used volume of fluid (VOF) model to simulate the evolution of underwater gas jet. XU Jia-wei, et al.[2-3]used mixture model based on FLUENT software to study the unsteady evolution process of nozzle air-water two-phase flow field at different underwater depths in rocket ignition experiment. CHEN Huan-long, et al.[4-5]calculated the initial flow of underwater gas jet using axisymmetric, inviscid and compressible Euler equations, revealed the gas-water interaction and shock wave shape as well as the formation and evolution process of jet air bags at jet initial stage, and confirmed the similarity of near flow fields of the gas jet in the water and gas. ZHU Lin, et al. used the dynamic mesh technique and mixed multiple flow model to get the effect of thrusting force on production, variation and development processes of water cavity and natural cavity hydrodynamics of projectile bodies[6]. YI Wen-jun, et al. used mixed multi-phase flow model of Fluet6.2 to conduct the simulation of resistance characteristics of underwater projectile[7]. Uhlman and Kinnas[8-10]and Fine used the velocity potential modified boundary element method to carry out nonlinear simulation on many cases of partial cavity flow. But the analysis of muzzle gas injecting sealed container is rare. This paper, using mixed multiple flow model and dynamic mesh technology, by solving N-S equations, carries numerical analysis and calculation for complex water-gas two-phase flow inside it and gets the distribution of the flow field inside the sealing device. Furthermore, the strength of the sealing device and bubble formation of supercavitation are analyzed.

        1 Control equation and numerical method

        1.1 Control equation

        The basic laws of conservation in the element include the following equations.

        Continuity equation

        (1)

        Momentum conservation theorem

        (2)

        Energy conservation equation

        (3)

        whereρstands for the density;t, the time;V, the velocity vector;Sm, the mass added to the continuous phase;p, the pressure born on upon fluid micelle;gandFrepresent gravity and other external volume forces exerted on the infinitesimal;τ, the viscous stress tensor on infinitesimal surface due to molecular viscosity;E, the fluid micelle gross energy;T, the temperature;Jj, the diffusion flux ofj;keff, effective thermal conductivity; andSh, the other customized heat sources.

        1.2 Cavitation model

        Regardless of the influence of the latent heat of evaporation, it is done in the isothermal process. Considering the pressure and bubble volume, Rayleigh-Plesset equation is written as

        (4)

        wherepBis the pressure inside the cavity, which is the sum of steam pressurepvand noncondensable gas partial pressurep;σis coefficient of surface tension. Simulation calculation functions for growth and rupture process of cavitation are given by

        (5)

        2 Dynamic meshing theory and computation model

        2.1 Dynamic mesh updating method

        FLUENT provides three methods of dynamic meshing movement to update the mesh after deformation, including spring smoothing method (smoothing), dynamic layer method (layering) and local mesh reconstruction method (remeshing). Spring smoothing method is to idealize the meshes between the nodes as a spring system. The movement of boundary nodes produces spring force between the nodes, and the force spreads along the downstream nodes in turn and eventually produces a new spring system. Dynamic layer method is that whether to increase or to decrease the mesh layer number is determined by the height of mesh layer close to the moving object surface. In FLUENT, when the mesh cell layer near the border increases or reduces to a certain extent, the mesh automatically splits or merges. Local mesh reconstruction method refers to the interpolation reconstruction method for distorted mesh because the movement of boundary may result in serious quality decline, even negative volume, which increases the difficultly for next solution. The methods adopted in this article to control mesh deformation are the spring smoothing method and the local mesh reconstruction method.

        2.2 Establishment of computation model

        A 3D model profile of a sealing device is shown in Fig.1. The whole sealing device is put under water at a depth of 5 m whose inside pressure and water pressure are equal before firing. The valve (unshown in Fig.1) of the sealing device under the water pressure acts as a seal. Mesh partition method is used for the calculation region, and mesh refinement is processed within the tube and sealing device. The farther from the tank, the sparser the mesh is. Boundary conditions include pressure export, pressure inlet and the solid wall boundary. The related parameters of pressure export and the depth of water are the same, and the pressure inlet parametes are defined based on field function. At the same time, FLUENT UDF is used to control projectile motion parameters. The initial velocity of the projectile is 890 m/s.

        Fig.1 Sealing device with a tube

        Fig.2 exhibits the computation mesh in flow field of the sealing device with a tube.

        Fig.2 Computation mesh in flow field of sealing device with a tube

        In the process of the projectile entering into the water, air, vapor and liquid water are allowed to interpenetrate each other and move at different speeds. VOF model is not applicable, mixture model should be shosen.

        FLUENT simulation settings: ① Calculation is based on three-phase flow, whose major phase is water and minor phases are vapor and compressible gas; ② Realizablek-εmodel is selected as turbulence model, PRESTO is as the pressure and velocity fields coupled mode; ③ Gravity effect should be considered, and the change curve of water pressure with depth is defined based on field function, Patch for repairing it in the computational domain; ④ Using the cavitation model is opened.

        Making unsteady calculation according to the above settings, we can get better convergence by selecting the appropriate iteration step length and the relaxion factor.

        3 Simulation and analysis

        Fig.3 is the nephogram of gas-liquid phase distribution at different times after projectile is out of the chamber. It can be seen that the moment the projectile firing, shock wave in front of projectile is generated by squeezing the air in front of projectile when the projectile moves at high speed in the tube. The gas expands rapidly after entering the sealing device, which results in pressure rising and part of the gas is pushed through the valve into the water forming bubbles. With the projectile velocity increasing, the shock waves in front of the projectile keep swelling and the gas inside the sealing device is constantly discharged into the water, and bubbles increase quickly. After 3.3 ms, gunpowder gas was totally released into the case and enough pressure is created to form a big bubble to wrap the projectile outside the sealing device. At 3.5 ms, the projectile moved throughout the sealing device into the big bubble in the water. At 3.8 ms, the projectile pierced the big bubble into the water, but was still wrapped by cavitation bubbles which are natural supercavitation. The bubbles in the valve in the seal device is helpful to provide the projectile with ventilated supercavitation. The device combines the ventilated supercavitation with natural supercavitation perfectly.

        Fig.3 Nephogram of gas-liquid phase distribution at different times after firing

        Fig.4 is the pressure distribution curve in axial direction of gun when the projectile is fired out of muzzle and goes into the sealing device at different times.

        It can be seen that at 3.3 ms, warhead is just fired out of the tube. The pressure distribution behind the projectile is linear, but the pressure in front of the projectile is higher than atmospheric pressure. This is because the moving projectile in the tube pushes the air in the tube leaked behind the projectile as well as gunpowder gas into the sealing device.

        At 3.5 ms, the projectile goes through the tube into the sealing device, the pressure behind the projectile spreads instantly in the sealing device. But unlike in the air, the pressure diffusion is restricted by the sealing device, resulting in high pressure in the sealing device. The projectile under high pressure continues to move forward and warhead appears in high pressure area. For high-speed projectile, it moves in the sealing device in a very short time, therefore, the effect on the projectile velocity can be neglected. At the time of 3.8 ms, the projectile goes into the water and the end of the projectile appears in low pressure area, which is conducive to the formation of cavitation bubble.

        At 4.0 ms, the pressure of gunpowder gas is still greater than the water pressure. It may be impossible for the gunpowder gas to be instantly released into the sealing device because of the small mouth of the valve, which means the water not going back to the sealing device under super-high firing frequency can be realized absolutely.

        During the whole process, the gunpowder gas is kept in the sealing device and the maxium shock wave can reach up to 100 MPa generated by gunpowder gas at high temperature and high pressure, which requires strong strength capacity of the sealing device.

        Fig.4 Pressure distribution curves in axial direction of gun at different times after projectile is out of muzzle

        4 Conclusion

        The following conclusion can be drawn based on numerical simulation on the complex two-phase flow in the sealing device of underwater weapon by using FLUENT. When the underwater weapon with the sealing device is fired, gunpowder gas pressure distribution in the sealing device is extremely complex. The high pressure also requires strong strength of the sealing device. With high-speed movement, the gas in front of the projectile is extruded and shock wave is formed. The gas in front of the projectile and gunpowder gas leaked from the clearance between projectile and tube expand rapidly after entering sealing device, which results in pressure rising in the sealing device. With the increasing of pressure, the gas inside the sealing device is discharged into the water. As the projectile entering the water, the pressure inside the sealing device is strong enough to prevent the reverse flow of water. Under ultra-high frequency shooting, the pressure of the sealing device can ensure that the water will not flow back at the depth of 25 m under water; at the same time, the actuating pressure of the sealing device caused by the projectile shock wave reach up to 100 MPa. The research results may have some significance to the design of underwater weapon sealing device and the formation of the projectile supercavitation.

        [1] Tang J N, Wang N F, Wei Shyy. Flow structure of gaseous jets injected into water for underwater propulsion. Acta Mechanica Sinica, 2011, 27(4): 461-472.

        [2] CHEN Huan-long, WANG Ning, LIU Hua-ping, et al. Investigation of nozzle gas jet characteristics with different launch depth underwater. Chinese Journal of Hydrodynamics, 2012, 27(6): 659-664.

        [3] XU Jia-wei. Numerical simulation of gas jet due to underwater missle of different ocean depths. Harbin: Harbin Institute of Technology, 2009: 33-89.

        [4] HE Xiao-yan, MA Han-dong, JI Chu-qun. Numerical simulation of the underwater gas jets formed. Journal of Hydrodynamics A, 2004, 19(2): 207-212.

        [5] HE Xiao-yan, MA Han-dong, WANG Qiang. Temporal instability of compressible 3D gas injected in liquids. Acta Aerodynamica Sinica, 2002, 20(1):78-83.

        [6] ZHU lin, YUAN Xu-long, WANG Ya-dong. The influence of the thrusting force on the characters of the high speed flow field. Journal of Air Force Engineering University(Natural Science Edition), 2014, 15(1): 10-14.

        [7] YI Wen-jun, WANG Zhong-yuan, XIONG Tian-hong, et al. Research on drag reduction characteristics of a underwater high-speed supercavitation projectile. Journal of Ballistics, 2008, 20(4): 1-4.

        [8] Uhalman J S. The surface singularity method applied to partially cavitating hydofoils. Journal of Ship Research, 1987: 23-36.

        [9] Uhalman J S. The surface singulaity or boundary integral method applied to supercavitating hydrofoils. Journal of ship Research, 1989, 31(1): 51-65.

        [10] Kinnas S A, Fine N E. A numerical nonlinear analysis of the flow around two-and three-dimensional partially cavitating hydrofoils. Journal of Fluid Mechanics, 1993(254): 151-181.

        基于動網(wǎng)格的水下發(fā)射裝置兩相流研究

        為防止水下武器發(fā)射時水流浸沒身管, 提高小口徑武器在水下發(fā)射的初速, 水下武器膛口前需設(shè)計密封裝置。 建立水下武器密封裝置模型的非結(jié)構(gòu)網(wǎng)格計算域模型, 利用計算流體動力學(xué)(CFD)、 動網(wǎng)格和UDF等相關(guān)知識, 通過求解N-S方程, 對密封裝置內(nèi)復(fù)雜的兩相流進(jìn)行了數(shù)值分析計算。 仿真結(jié)果表明, 密封裝置內(nèi)排出的氣體有利于射彈超空炮的形成; 此外, 在水下5 m發(fā)射時彈前激波和彈后激波對箱體沖擊作用力峰值可達(dá)100 MPa, 因此要求密封裝置要有足夠的強(qiáng)度。

        兩相流; 超空泡; 密封裝置; 計算流體動力學(xué); 動網(wǎng)格

        ZHANG Xue-wei, LI Qiang, LV Meng-rou. Numerical simulation of two-phase flow field in underwater sealing device based on dynamic mesh. Journal of Measurement Science and Instrumentation, 2015, 6(3): 253-257.[

        張學(xué)偉, 李 強(qiáng), 呂夢柔

        (中北大學(xué) 機(jī)電工程學(xué)院, 山西 太原 030051)

        10.3969/j.issn.1674-8042.2015.03.009]

        Received date: 2015-05-06 Foundation item: National Natural Science Foundation of China (No.51175481)

        ZHANG Xue-wei (1192579296@qq.com)

        1674-8042(2015)03-0253-05 doi: 10.3969/j.issn.1674-8042.2015.03.009

        CLD number: TJ65; TP391.9 Document code: A

        猜你喜歡
        中北大學(xué)激波工程學(xué)院
        《中北大學(xué)學(xué)報(社會科學(xué)版)》征稿啟事
        福建工程學(xué)院
        福建工程學(xué)院
        中北大學(xué)信創(chuàng)產(chǎn)業(yè)學(xué)院入選首批現(xiàn)代產(chǎn)業(yè)學(xué)院
        一種基于聚類分析的二維激波模式識別算法
        基于HIFiRE-2超燃發(fā)動機(jī)內(nèi)流道的激波邊界層干擾分析
        《中北大學(xué)學(xué)報(自然科學(xué)版)》征稿簡則
        福建工程學(xué)院
        有機(jī)相化學(xué)鍍鋁法制備Al/石墨烯復(fù)合材料粉末
        斜激波入射V形鈍前緣溢流口激波干擾研究
        亚洲—本道中文字幕久久66| 人妻少妇精品视频三区二区一区| 48沈阳熟女高潮嗷嗷叫| 亚洲欧美日韩中文天堂| 久久久久久久性潮| 岛国熟女一区二区三区| 中文字幕中文字幕三区| 精品国产品香蕉在线| 内射人妻少妇无码一本一道| 大胆欧美熟妇xxbbwwbw高潮了| 99在线国产视频| 中文乱码字幕在线亚洲av | 无码国内精品久久人妻| 永久免费av无码入口国语片| 一级呦女专区毛片| 亚洲人妻av在线播放| 精品天堂色吊丝一区二区| 亚洲色欲色欲综合网站| 色吧综合网| 91在线视频视频在线| 丝袜美腿av在线观看| 国产亚洲午夜高清国产拍精品| 久久久亚洲经典视频| 天堂影院久久精品国产午夜18禁| 亚洲成av人片不卡无码| 国产在线无码制服丝袜无码| 亚欧同人精品天堂| 日韩女优图播一区二区| 国产精品无码无片在线观看3d| 中文字幕第1页中文字幕在| 国产自拍精品视频免费观看| 日本中文字幕一区二区有码在线| 国产在线精品一区二区| 亚洲AV秘 无套一区二区三区 | 91精品国产无码在线观看| 美女视频一区二区三区在线| 午夜福利试看120秒体验区| 日韩欧美亚洲中字幕在线播放| 性感的小蜜桃在线观看| 国产乱人偷精品人妻a片| 久久99精品久久久久久久清纯|