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 m...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 5munder 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.展开更多
An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （...An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （AUSM） scheme, and a 5-stage explicit Runge-Kutta scheme is adopted in the time integration. To avoid the influence of the magnetic field divergence created during the simulation, the hyperbolic divergence cleaning method is introduced. The shock-capturing properties of the method are verified by solving the MHD shock-tube problem. Then the 2-D nozzle flow with the magnetic field is numerically simulated on the unstructured mesh. Computational results demonstrate the effects of the magnetic field and agree well with those from references.展开更多
Computational fluid dynamics(CFD) has recently emerged as an effective tool for the investigation of the hydraulic parameters and efficiency of tray towers.The computation domain was established for two types of orien...Computational fluid dynamics(CFD) has recently emerged as an effective tool for the investigation of the hydraulic parameters and efficiency of tray towers.The computation domain was established for two types of oriented valves within a tray and meshed into two parts with different grid types and sizes.The volume fraction correlation concerning inter-phase momentum transfer source was fitted based on experimental data,and built in UDF for simulation.The flow pattern of oriented valve tray under different operating conditions was simulated under Eulerian-Eulerian framework with realizable k-ε model.The predicted liquid height from CFD simulation was in good agreement with the results of pressure drop and volume fraction correlations.Meanwhile,the velocity distribution and volume fraction of the two phases were demonstrated and analyzed,which are useful in design and analysis of the column trays.展开更多
In the smoothed particle hydrodynamics (SPH) method, a meshless interpolation scheme is needed for the unknown function in order to discretize the governing equation.A particle approximation method has so far been use...In the smoothed particle hydrodynamics (SPH) method, a meshless interpolation scheme is needed for the unknown function in order to discretize the governing equation.A particle approximation method has so far been used for this purpose.Traditional particle interpolation (TPI) is simple and easy to do, but its low accuracy has become an obstacle to its wider application.This can be seen in the cases of particle disorder arrangements and derivative calculations.There are many different methods to improve accuracy, with the moving least square (MLS) method one of the most important meshless interpolation methods.Unfortunately, it requires complex matrix computing and so is quite time-consuming.The authors developed a simpler scheme, called higher-order particle interpolation (HPI).This scheme can get more accurate derivatives than the MLS method, and its function value and derivatives can be obtained simultaneously.Although this scheme was developed for the SPH method, it has been found useful for other meshless methods.展开更多
For the numerical simulation of compressible flows,normally different mesh sizes are expected in different regions.For example,smaller mesh sizes are required to improve the local numerical resolution in the regions w...For the numerical simulation of compressible flows,normally different mesh sizes are expected in different regions.For example,smaller mesh sizes are required to improve the local numerical resolution in the regions where the physical variables vary violently(for example,near the shock waves or in the boundary layers)and larger elements are expected for the regions where the solution is smooth.h-adaptive mesh has been widely used for complex flows.However,there are two difficulties when employing h-adaptivity for high-order discontinuous Galerkin(DG)methods.First,locally curved elements are required to precisely match the solid boundary,which significantly increases the difficulty to conduct the"refining"and"coarsening"operations since the curved information has to be maintained.Second,h-adaptivity could break the partition balancing,which would significantly affect the efficiency of parallel computing.In this paper,a robust and automatic h-adaptive method is developed for high-order DG methods on locally curved tetrahedral mesh,for which the curved geometries are maintained during the h-adaptivity.Furthermore,the reallocating and rebalancing of the computational loads on parallel clusters are conducted to maintain the parallel efficiency.Numerical results indicate that the introduced h-adaptive method is able to generate more reasonable mesh according to the structure of flow-fields.展开更多
Local scour, a non-negligible factor in hydraulic engineering, endangers the safety of hydraulic structures. In this work, a numerical model for simulating local scour was constructed, based on the open source code co...Local scour, a non-negligible factor in hydraulic engineering, endangers the safety of hydraulic structures. In this work, a numerical model for simulating local scour was constructed, based on the open source code computational fluid dynamics model Open FOAM. We consider both the bedload and suspended load sediment transport in the scour model and adopt the dynamic mesh method to simulate the evolution of the bed elevation. We use the finite area method to project data between the three-dimensional flow model and the two-dimensional(2D) scour model. We also improved the 2D sand slide method and added it to the scour model to correct the bed bathymetry when the bed slope angle exceeds the angle of repose. Moreover, to validate our scour model, we conducted and compared the results of three experiments with those of the developed model. The validation results show that our developed model can reliably simulate local scour.展开更多
Numerical simulation on conjugate heat transfer of an internal cooled turbine vane was carried out. Numerical techniques employed included the third-order accuracy TVD scheme, multi-block structured grids and the tech...Numerical simulation on conjugate heat transfer of an internal cooled turbine vane was carried out. Numerical techniques employed included the third-order accuracy TVD scheme, multi-block structured grids and the technique of arbitrary curved mesh. Comparison between results of commercial CFD codes with several turbulence models and those of this code shows that it is incorrect of commercial CFD codes to predict the thermal boundary layer with traditional turbulence models, and that turbulence models considering transition lead to more accurate heat transfer in thermal boundary layer with some reliability and deficiency yet. The results of this code are close to those of CFX with transition model.展开更多
基金National Natural Science Foundation of China(No.51175481)
文摘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 5munder 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.
文摘An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （AUSM） scheme, and a 5-stage explicit Runge-Kutta scheme is adopted in the time integration. To avoid the influence of the magnetic field divergence created during the simulation, the hyperbolic divergence cleaning method is introduced. The shock-capturing properties of the method are verified by solving the MHD shock-tube problem. Then the 2-D nozzle flow with the magnetic field is numerically simulated on the unstructured mesh. Computational results demonstrate the effects of the magnetic field and agree well with those from references.
文摘Computational fluid dynamics(CFD) has recently emerged as an effective tool for the investigation of the hydraulic parameters and efficiency of tray towers.The computation domain was established for two types of oriented valves within a tray and meshed into two parts with different grid types and sizes.The volume fraction correlation concerning inter-phase momentum transfer source was fitted based on experimental data,and built in UDF for simulation.The flow pattern of oriented valve tray under different operating conditions was simulated under Eulerian-Eulerian framework with realizable k-ε model.The predicted liquid height from CFD simulation was in good agreement with the results of pressure drop and volume fraction correlations.Meanwhile,the velocity distribution and volume fraction of the two phases were demonstrated and analyzed,which are useful in design and analysis of the column trays.
基金Supported by the National Natural Science Foundation of China under Grant No.10572041,50779008Doctoral Fund of Ministry of Education of China under Grant No.20060217009
文摘In the smoothed particle hydrodynamics (SPH) method, a meshless interpolation scheme is needed for the unknown function in order to discretize the governing equation.A particle approximation method has so far been used for this purpose.Traditional particle interpolation (TPI) is simple and easy to do, but its low accuracy has become an obstacle to its wider application.This can be seen in the cases of particle disorder arrangements and derivative calculations.There are many different methods to improve accuracy, with the moving least square (MLS) method one of the most important meshless interpolation methods.Unfortunately, it requires complex matrix computing and so is quite time-consuming.The authors developed a simpler scheme, called higher-order particle interpolation (HPI).This scheme can get more accurate derivatives than the MLS method, and its function value and derivatives can be obtained simultaneously.Although this scheme was developed for the SPH method, it has been found useful for other meshless methods.
基金supported by the funding of the Key Laboratory of Aerodynamic Noise Control(No.ANCL20190103)the State Key Laboratory of Aerodynamics(No.SKLA20180102)+1 种基金the Aeronautical Science Foundation of China(Nos.2018ZA52002,2019ZA052011)the National Natural Science Foundation of China(Nos.61672281,61732006)。
文摘For the numerical simulation of compressible flows,normally different mesh sizes are expected in different regions.For example,smaller mesh sizes are required to improve the local numerical resolution in the regions where the physical variables vary violently(for example,near the shock waves or in the boundary layers)and larger elements are expected for the regions where the solution is smooth.h-adaptive mesh has been widely used for complex flows.However,there are two difficulties when employing h-adaptivity for high-order discontinuous Galerkin(DG)methods.First,locally curved elements are required to precisely match the solid boundary,which significantly increases the difficulty to conduct the"refining"and"coarsening"operations since the curved information has to be maintained.Second,h-adaptivity could break the partition balancing,which would significantly affect the efficiency of parallel computing.In this paper,a robust and automatic h-adaptive method is developed for high-order DG methods on locally curved tetrahedral mesh,for which the curved geometries are maintained during the h-adaptivity.Furthermore,the reallocating and rebalancing of the computational loads on parallel clusters are conducted to maintain the parallel efficiency.Numerical results indicate that the introduced h-adaptive method is able to generate more reasonable mesh according to the structure of flow-fields.
基金the State Key Laboratory of Hydraulic Engineering Simulation and Safety Foundation (No. HESS-1412)the National Science Fund (No. 51179178)the 111 Project (No. B14028)
文摘Local scour, a non-negligible factor in hydraulic engineering, endangers the safety of hydraulic structures. In this work, a numerical model for simulating local scour was constructed, based on the open source code computational fluid dynamics model Open FOAM. We consider both the bedload and suspended load sediment transport in the scour model and adopt the dynamic mesh method to simulate the evolution of the bed elevation. We use the finite area method to project data between the three-dimensional flow model and the two-dimensional(2D) scour model. We also improved the 2D sand slide method and added it to the scour model to correct the bed bathymetry when the bed slope angle exceeds the angle of repose. Moreover, to validate our scour model, we conducted and compared the results of three experiments with those of the developed model. The validation results show that our developed model can reliably simulate local scour.
基金Sponsored by the National Natural Science Foundation of China (Grant No.5047028 and 50476017)
文摘Numerical simulation on conjugate heat transfer of an internal cooled turbine vane was carried out. Numerical techniques employed included the third-order accuracy TVD scheme, multi-block structured grids and the technique of arbitrary curved mesh. Comparison between results of commercial CFD codes with several turbulence models and those of this code shows that it is incorrect of commercial CFD codes to predict the thermal boundary layer with traditional turbulence models, and that turbulence models considering transition lead to more accurate heat transfer in thermal boundary layer with some reliability and deficiency yet. The results of this code are close to those of CFX with transition model.
