The aim of this paper is to simulate and study the early moments of the reactive ballistics of a large caliber projectile fired from a gun,combining 0D and 2D axisymmetric Computational Fluid Dynamics(CFD)approaches.F...The aim of this paper is to simulate and study the early moments of the reactive ballistics of a large caliber projectile fired from a gun,combining 0D and 2D axisymmetric Computational Fluid Dynamics(CFD)approaches.First,the methodology is introduced with the development of an interior ballistics(IB)lumped parameter code(LPC),integrating an original image processing method for calculating the specific regression of propellant grains that compose the gun propellant.The ONERA CFD code CEDRE,equipped with a Dynamic Mesh Technique(DMT),is then used in conjunction with the developed LPC to build a dedicated methodology to calculate IB.First results obtained on the AGARD gun and 40 mm gun test cases are in a good agreement with the existing literature.CEDRE is also used to calculate inter-mediate ballistics(first milliseconds of free flight of the projectile)with a multispecies and reactive approach either starting from the gun muzzle plane or directly following IB.In the latter case,an inverse problem involving a Latin hypercube sampling method is used to find a gun propellant configuration that allows the projectile to reach a given exit velocity and base pressure when IB ends.The methodology developed in this work makes it possible to study the flame front of the intermediate flash and depressurization that occurs in a base bleed(BB)channel at the gun muzzle.Average pressure variations in the BB channel during depressurization are in good agreement with literature.展开更多
In 2012,the International Agency for Research on Cancer(IARC)classified diesel particulate matter(DPM)as a carcinogen to human.With the increased usage of diesel equipment in underground mines,miners have a high risk ...In 2012,the International Agency for Research on Cancer(IARC)classified diesel particulate matter(DPM)as a carcinogen to human.With the increased usage of diesel equipment in underground mines,miners have a high risk of over-exposure to DPM,which has drawn many concerns from the public.This study used computational fluid dynamics(CFD)to analyse the DPM dispersion and concentration distribution characteristics in an underground development face based on an onsite experiment.The DPM emitted from a moving loader under a forcing auxiliary ventilation system was simulated.The motion of the load-haul-dump(LHD)in the tunnel was represented by a dynamic mesh method.The species transport approach was applied to study the DPM behaviours.High DPM concentration zones were then identified based on the simulation results.The results could provide guidelines for work practices and be helpful to an optimum auxiliary ventilation design to reduce underground miner exposure.展开更多
Autonomous underwater vehicles （AUVs） navigating on the sea surface are usually required to complete the communication tasks in complex sea conditions. The movement forms and flow field characteristics of a multi-mo...Autonomous underwater vehicles （AUVs） navigating on the sea surface are usually required to complete the communication tasks in complex sea conditions. The movement forms and flow field characteristics of a multi-moving state AUV navigating in head sea at high speed were studied. The mathematical model on longitudinal motion of the high-speed AUV in head sea was established with considering the hydrodynamic lift based on strip theory, which was solved to get the heave and pitch of the AUV by Gaussian elimination method. Based on this, computational fluid dynamics （CFD） method was used to establish the mathematical model of the unsteady viscous flow around the AUV with considering free surface effort by using the Reynolds-averaged Navier-Stokes （RANS） equations, shear-stress transport （SST） k-w model and volume of fluid （VOF） model. The three-dimensional numerical wave in the computational field was realized through defining the unsteady inlet boundary condition. The motion forms of the AUV navigating in head sea at high speed were carried out by the program source code of user-defined function （UDF） based on dynamic mesh. The hydrodynamic parameters of the AUV such as drag, lift, pitch torque, velocity, pressure and wave profile were got, which reflect well the real ambient flow field of the AUV navigating in head sea at high speed. The computational wave profile agrees well with the experimental phenomenon of a wave-piercing surface vehicle. The force law of the AUV under the impacts of waves was analyzed qualitatively and quantitatively, which provides an effective theoretical guidance and technical support for the dynamics research and shape design of the AUV in real complex environnaent.展开更多
To improve mesh quality for KIVA-3V a method has been developed for rapid mesh generation and dynamic mesh management with moving valves for internal combustion engines. Two phases are included in rapid mesh generatio...To improve mesh quality for KIVA-3V a method has been developed for rapid mesh generation and dynamic mesh management with moving valves for internal combustion engines. Two phases are included in rapid mesh generation： the initial mesh generation and the mesh pre-treatment. In the second step （pre-treatment）, the connectivity of those cells is generated by a new algorithm added to the KIVA-3V code after the initial mesh generated. In dynamic mesh management phase, a new rezoning algorithm is developed and the basic principle is that the rezoning starts from the moving part. The movement of the adjustment is treated as an ＂earth quake wave＂ propagating to the surrounding vertexes. The amount of coordinate adjustment of the surrounding vertexes is determined by the movement of the epicenter and the distance between the vertexes and the ＂epicenter＂. Finally, a real IC engine mesh is generated and managed aceording to the new method. It gives a new theory and a new method for creating and managing the mesh in IC engine.展开更多
Three-dimensional simulations of diesel particulate matter （DPM） distribution inside a single straight entry for the Load-Haul-Dump loader （LHD）-truck loading and truck hauling operations were conducted by using A...Three-dimensional simulations of diesel particulate matter （DPM） distribution inside a single straight entry for the Load-Haul-Dump loader （LHD）-truck loading and truck hauling operations were conducted by using ANSYS FLUENT computational fluid dynamics software. The loading operation was performed for a fixed period of 3 min. The dynamic mesh technique in FLUENT was used to study the impact of truck motion on DPM distribution. The resultant DPM distributions are presented for the cases when the truck were driving upstream and downstream of the loading face. Interesting phenomena were revealed in the study including the piston effect, layering of DPM in the roof region, and backflow of diesel exhaust against ventilation. The results from the simulation can be used to determine if the areas inside the face area and straight entry exceed the current U.S. regulatory requirement for DPM concentration （〉160 pg/m3）. This research can guide the selection of DPM reduction strategies and improve the working practices for the underground miners.展开更多
Modeling and experiments of the growth of ash deposits during the combustion in a drop tube furnace are presented in this paper.An existing deposition model for ash deposit was used,which involves relationship among t...Modeling and experiments of the growth of ash deposits during the combustion in a drop tube furnace are presented in this paper.An existing deposition model for ash deposit was used,which involves relationship among the force of gravity,elastic rebound and adhesion forces acting at the moment of ash particle impaction.Important parameters of this model were determined by the experimental data.The influence of particle size and velocity on deposit geometry was investigated.The growth of ash deposits involving various inlet velocities,tube diameters and tube arrangements were simulated with dynamic mesh technique in Fluent.The simulation results show that particle normal impacted velocity plays an important role due to elastic bounce force term.Deposited geometry was influenced by side velocity(velocity of inlet⁃2),and large velocity has an obvious effect to change the shape of deposition and postpone the steady time.展开更多
A new method is developed to assess and analyze the dynamic performance of hydrostatic bearing oil film by using an amulets-layer dynamic mesh technique. It is implemented using C Language to compile the UDF program o...A new method is developed to assess and analyze the dynamic performance of hydrostatic bearing oil film by using an amulets-layer dynamic mesh technique. It is implemented using C Language to compile the UDF program of a single oil film of the hydrostatic bearing. The effects of key lubrication parameters of the hydrostatic bearing are evaluated and analyzed under various working conditions,i.e. under no-load,a load of 40 t,a full load of 160 t,and the rotation speed of 1r/min,2r/min,4r/min,8r/min,16r/min,32r/min. The transient data of oil film bearing capacity under different load and rotation speed are acquired for a total of 18 working conditions during the oil film thickness changing. It allows the effective prediction of dynamic performance of large size hydrostatic bearing. Experiments on hydrostatic bearing oil film have been performed and the results were used to define the boundary conditions for the numerical simulations and validate the developed numerical model. The results showed that the oil film thickness became thinner with the increase of the operating time of the hydrostatic bearing,both the oil film rigidity and the oil cavity pressure increased significantly,and the increase of the bearing capacity was inversely proportional to the cube of the change of the film thickness. Meanwhile,the effect of the load condition on carrying capacity of large size static bearing was more important than the speed condition. The error between the simulation value and the experimental value was 4.25%.展开更多
A Natural Ecological Wastewater Treatment System （NEWTS） is usually built on natural terrain with necessary topography modification to improve water flowing route and pattern, and then the topography modified NEWTS ...A Natural Ecological Wastewater Treatment System （NEWTS） is usually built on natural terrain with necessary topography modification to improve water flowing route and pattern, and then the topography modified NEWTS should also have a reasonable water storage volume and hydraulic retention time so as to achieve the anticipated water purification effect. In this study, the dynamic mesh technique based on the finite element method and element storativity coefficients was presented to develop a two-dimensional hydrodynamic and water quality model, which was used to optimize the design of NEWTS under the dynamic land-water boundary due to various water storage volume. The models were employed in the optimized design of NEWTS from a large abandoned coal mine, which purifies the polluted water flowing into a large water storage lake, as part of the East Route South-to-North Water Transfer Project in China. Specifically, the natural topography modification scheme was presented, and further, a reasonable water storage volume and hydraulic residence time was obtained, based on the reasonable estimation of roughness coefficient and pollutant removal rate of the NEWTS with phragmites communis.展开更多
A grate cooler is key equipment in quenching clinker and recovering heat in cement production. A two-dimensional numerical model based on a 5000 t/d cement plant is proposed to for a study on the gas-solid coupled hea...A grate cooler is key equipment in quenching clinker and recovering heat in cement production. A two-dimensional numerical model based on a 5000 t/d cement plant is proposed to for a study on the gas-solid coupled heat transfer process between the cooling air and clinker in a grate cooler. In this study, we use the Fluent dynamic mesh technique and porous media model through which the transient temperature distribution with the clinker motion process and steady temperature and pressure distribution are obtained. We validate the numerical model with the operating data of the cooling air outlet temperature. Then, we discuss the amount of mid-temperature air outlet and average diameter of clinker particles, which affect the heat effective utilization and cooling air pressure drop in clinker layer. We found that after adding one more mid-temperature air outlet, the average temperature of the air flowing into the heat recovery boiler increases by 29.04℃ and the ratio of heat effective utilization increases by 5.3%. This means heat recovery is more effective on adding one more mid-temperature air outlet. Further, the smaller the clinker particles, the more is the pressure drop in clinker layer; thus more power consumption is needed by the cooling fan.展开更多
It is of great significance to develop a high-efficiency and low-noise propeller optimization method for new-generation propeller aircraft design.Coupled with free form deformation method,dynamic mesh interpolation te...It is of great significance to develop a high-efficiency and low-noise propeller optimization method for new-generation propeller aircraft design.Coupled with free form deformation method,dynamic mesh interpolation technology,optimization algorithm,surrogate model,aerodynamic calculation and aeroacoustic prediction model module,the integrated aerodynamic and aeroacoustic design method of propeller is built.The optimization design for the six-blade propeller is carried out.The non-reduction in efficiency,thrust coefficient and the minimum of aerodynamic noise is treated as the optimization design objective.The spatial vorticity distribution of the propeller before and after the design is also analyzed by using unsteady computational fluid dynamics method.The results show that the optimized propeller can effectively reduce the aerodynamic noise level.The maximum total sound pressure level can be reduced by 5 dB without reducing its aerodynamic performance.The developed method has good application potential in low-noise optimization design of propeller and other rotating machinery.展开更多
The flow in the positive displacement blower is very complex.The existing two-dimensional numerical simulation cannot provide the detailed flow information,especially flow characteristics along the axial direction,whi...The flow in the positive displacement blower is very complex.The existing two-dimensional numerical simulation cannot provide the detailed flow information,especially flow characteristics along the axial direction,which is unfavorable to improve the performance of positive displacement blower.To investigate the effects of spiral inlet and outlet on the aerodynamic performance of positive displacement blower,three-dimensional unsteady flow characteristics in a three-lobe positive displacement blower with and without the spiral inlet and outlet are simulated by solving Navier-Stokes equations coupled with RNG k-ε turbulent model.In the numerical simulation,the dynamic mesh technique and overset mesh updating method are used.The computational results are compared with the experimental measurements on the variation of flow rate with the outlet pressure to verify the validity of the numerical method presented.The results show that the mass flow rate with the change of pressure is slightly affected by the application of spiral inlet and outlet,but the internal flow state is largely affected.In the exhaust region,the fluctuations of pressure,velocity and temperature as well as the average values of velocity are significantly reduced.This illustrates that the spiral outlet can effectively suppress the fluctuations of pressure,thus reducing reflux shock and energy dissipation.In the intake area,the average value of pressure,velocity and temperature are slightly declined,but the fluctuations of them are significantly reduced,indicating that the spiral inlet plays the role in making the flow more stable.The numerical results obtained reveal the three-dimensional flow characteristics of the positive displacement blower with spiral inlet and outlet,and provide useful reference to improve performance and empirical correction in the noise-reduction design of the positive displacement blowers.展开更多
The fluid–structure interaction and aerodynamic shape optimization usually involve the moving or deforming boundaries, thus the dynamic mesh techniques are the key techniques to cope with such deformation. A novel dy...The fluid–structure interaction and aerodynamic shape optimization usually involve the moving or deforming boundaries, thus the dynamic mesh techniques are the key techniques to cope with such deformation. A novel dynamic mesh method was developed based on the Delaunay graph in this paper. According to the Delaunay graph, the mesh points were divided into groups. In each group, a factor ranging from 0 to 1 was calculated based on the area/volume ratio. By introducing a proper function for this factor, this method can control the mesh quality with high efficiency. Several test cases were compared with other dynamic mesh methods regarding mesh quality and CPU time, such as radial basis function method and Delaunay graph mapping method.展开更多
Currently, the flow field of annular seals disturbed by the circular whirl motion of rotors is usually solved using computational fluid dynamics(CFD) to evaluate the five rotordynamic coefficients. The simulations are...Currently, the flow field of annular seals disturbed by the circular whirl motion of rotors is usually solved using computational fluid dynamics(CFD) to evaluate the five rotordynamic coefficients. The simulations are based on the traditional quasi-steady method. In this work, an improved quasi-steady method along with the transient method was presented to compute the rotordynamic coefficients of a long seal. By comparisons with experimental data, the shortcomings of quasi-steady methods have been identified. Then, the effects of non-uniform incoming flow on seal dynamic coefficients were studied by transient simulations. Results indicate that the long seal has large cross stiffness k and direct mass M which are not good for rotor stability, while the transient method is more suitable for the long seal for its excellent performance in predicting M. When the incoming flow is non-uniform, the stiffness coefficients vary with the eccentric directions. Based on the rotordynamic coefficients under uniform incoming flow, the linearized fluid force formulas, which can consider the effects of non-uniform incoming flow, have been presented and can well explain the varying-stiffness phenomenon.展开更多
While passing through a lock, a ship usually undergoes a steady forward motion at low speed. Owing to the size restriction of lock chamber, the shallow water and bank effects on the hydrodynamic forces acting on the s...While passing through a lock, a ship usually undergoes a steady forward motion at low speed. Owing to the size restriction of lock chamber, the shallow water and bank effects on the hydrodynamic forces acting on the ship may be remarkable, which may have an adverse effect on navigation safety. However, the complicated hydrodynamics is not yet fully understood. This paper focuses on the hydrodynamic forces acting on a ship passing through a lock. The unsteady viscous flow and hydrodynamic forces are calculated by applying an unsteady RANS code with a RNG k-e turbulence model. User-defined function （UDF） is compiled to define the ship motion. Meanwhile, the grid regeneration is dealt with by using the dynamic mesh method and sliding interface technique. Numerical study is carried out for a bulk carrier ship passing through the Pierre Vandamme Lock in Zeebrugge at the model scale. The proposed method is validated by comparing the numerical results with the data of captive model tests. By analyzing the numerical results obtained at different speeds, water depths and eccentricities, the influences of speed, water depth and eccentricity on the hydrodynamic forces are illustrated. The numerical method proposed in this paper can qualitatively predict the ship-lock hydrodynamic interaction. It can provide certain guidance on the manoeuvring and control of ships passing through a lock.展开更多
This paper focuses on the effect of the phase offset of Leading-Edge(LE)morphing on the aerodynamic characteristics of a pitching NACA0012 airfoil.Assuming an unstretched camber and using polynomial interpolation,an e...This paper focuses on the effect of the phase offset of Leading-Edge(LE)morphing on the aerodynamic characteristics of a pitching NACA0012 airfoil.Assuming an unstretched camber and using polynomial interpolation,an explicit expression for LE nonlinear morphing is proposed and implemented for the large pitching motion of the airfoil.Flow field results and aerodynamic forces are obtained by solving the unsteady Reynolds-averaged Navier-Stokes equations for both the airfoil’s pitching motion and LE morphing.Furthermore,the index of instantaneous aerodynamic power is used to quantify the work done by the airflow in a dynamic process.According to the instantaneous aerodynamic power and energy map,which denotes the energy transfer between the airfoil’s oscillation and flow field,the airfoil is subject to stall flutter.The results show that LE morphing with an optimal phase offset of 315°reduces the energy extraction from the flow field,suppressing the stall flutter instability.This optimal phase offset is effective at different pitching axis positions of the airfoil.The results signify that LE morphing can suppress stall flutter by advancing the occurrence of the first LE vortex and increasing the nose-down moment during the upstroke period.展开更多
The performance of a newly designed tri-lobe industrial lobe pump of high capacity is simulated by using commercial CFD solver Ansys Fluent. A combination of user-defined-functions and meshing strategies is employed t...The performance of a newly designed tri-lobe industrial lobe pump of high capacity is simulated by using commercial CFD solver Ansys Fluent. A combination of user-defined-functions and meshing strategies is employed to capture the rotation of the lobes. The numerical model is validated by comparing the simulated results with the literature values. The processes of suction, displacement, compression and exhaust are accurately captured in the transient simulation. The fluid pressure value remains in the range of inlet pressure value till the processes of suction and displacement are over. The instantaneous process of compression is accurately captured in the simulation. The movement of a particular working chamber is traced along the gradual degree of lobe’s rotation. At five different degrees of lobe’s rotation, pressure contour plots are reported which clearly shows the pressure values inside the working chamber. Each pressure value inside the working chamber conforms to the particular process in which the working chamber is operating. Finally, the power requirement at the shaft of rotation is estimated from the simulated values. The estimated value of power requirement is 3.61 BHP FHP whereas the same calculated theoretically is 3 BHP FHP. The discrepancy is attributed to the assumption of symmetry of blower along the thickness.展开更多
A feasible method was proposed to improve the vibration intensity of screen surface via application of a new type elastic screen surface with multi degree of freedom(NTESSMDF). In the NTESSMDF, the primary robs were c...A feasible method was proposed to improve the vibration intensity of screen surface via application of a new type elastic screen surface with multi degree of freedom(NTESSMDF). In the NTESSMDF, the primary robs were coupled to the main screen structure with ends embedded into the elastomers, and the secondary robs were attached to adjacent two primary robs with elastic bands. The dynamic model of vibrating screen with NTESSMDF was established based on Lagrange's equation and the equivalent stiffnesses of the elastomer and elastic band were calculated. According to numerical simulation using the 4th order Runge-Kutta method, the vibration intensity of screen surface can be enhanced substantially with an averaged acceleration amplitude increasing ratio of 72.36%. The primary robs and secondary robs vibrate inversely in steady state, which would result in the friability of materials and avoid stoppage. The experimental results validate the dynamic characteristics with acceleration amplitude rising by62.93% on average, which demonstrates the feasibility of NTESSMDF.展开更多
Aiming at the high angle of attack pull-up and multi-channel roll pull-up coupling problems of high maneuvering aircraft, this paper establishes the flight attitude control rate by means of unsteady flow numerical sol...Aiming at the high angle of attack pull-up and multi-channel roll pull-up coupling problems of high maneuvering aircraft, this paper establishes the flight attitude control rate by means of unsteady flow numerical solution, dynamic unstructured nested mesh assembly method and numerical solution method of flight mechanics equation. On this basis, a virtual flight simulation platform integrating pneumatics, motion and control is established. Based on this virtual flight simulation platform, F-16 aircraft is simulated by high angle of attack pull-up flight mode and multi-channel roll pull-up coupling flight mode. Finally, the influence of rudder on the yaw control channel is investigated. The results show that the numerical virtual flight simulation platform established in this paper has the ability to simulate maneuvering flight of aircraft.展开更多
A two-dimensional computational fluid dynamics model was established to simulate the friction stir butt-welding of 6061 aluminum alloy. The dynamic mesh method was applied in this model to make the tool move forward a...A two-dimensional computational fluid dynamics model was established to simulate the friction stir butt-welding of 6061 aluminum alloy. The dynamic mesh method was applied in this model to make the tool move forward and rotate in a manner similar to a real tool, and the calculated volumetric source of energy was loaded to establish a similar thermal environment to that used in the experiment. Besides, a small piece of zinc stock was embedded into the workpiece as a trace element. Temperature fields and vector plots were determined using a finite volume method, which was indirectly verified by traditional metallography. The simulation result indicated that the temperature distribution was asymmetric but had a similar tendency on the two sides of the welding line. The maximum temperature on the advancing side was approximately 10 K higher than that on the retreating side. Furthermore, the precise process of material flow behavior in combination with streamtraces was demonstrated by contour maps of the phases. Under the shearing force and forward extrusion pressure, material located in front of the tool tended to move along the tangent direction of the rotating tool. Notably, three whirlpools formed under a special pressure environment around the tool, resulting in a uniform composition distribution.展开更多
基金the French Defense Innovation Agency (AID)the French Procurement Agency for Armament (DGA)ONERA's scientific direction for funding and supporting the present work
文摘The aim of this paper is to simulate and study the early moments of the reactive ballistics of a large caliber projectile fired from a gun,combining 0D and 2D axisymmetric Computational Fluid Dynamics(CFD)approaches.First,the methodology is introduced with the development of an interior ballistics(IB)lumped parameter code(LPC),integrating an original image processing method for calculating the specific regression of propellant grains that compose the gun propellant.The ONERA CFD code CEDRE,equipped with a Dynamic Mesh Technique(DMT),is then used in conjunction with the developed LPC to build a dedicated methodology to calculate IB.First results obtained on the AGARD gun and 40 mm gun test cases are in a good agreement with the existing literature.CEDRE is also used to calculate inter-mediate ballistics(first milliseconds of free flight of the projectile)with a multispecies and reactive approach either starting from the gun muzzle plane or directly following IB.In the latter case,an inverse problem involving a Latin hypercube sampling method is used to find a gun propellant configuration that allows the projectile to reach a given exit velocity and base pressure when IB ends.The methodology developed in this work makes it possible to study the flame front of the intermediate flash and depressurization that occurs in a base bleed(BB)channel at the gun muzzle.Average pressure variations in the BB channel during depressurization are in good agreement with literature.
基金the Minerals Research Institute of Western Australia(M495)the computation resources provided by the Pawsey Supercomputing Centre with funding from the Australian Government and the Government of Western Australia.
文摘In 2012,the International Agency for Research on Cancer(IARC)classified diesel particulate matter(DPM)as a carcinogen to human.With the increased usage of diesel equipment in underground mines,miners have a high risk of over-exposure to DPM,which has drawn many concerns from the public.This study used computational fluid dynamics(CFD)to analyse the DPM dispersion and concentration distribution characteristics in an underground development face based on an onsite experiment.The DPM emitted from a moving loader under a forcing auxiliary ventilation system was simulated.The motion of the load-haul-dump(LHD)in the tunnel was represented by a dynamic mesh method.The species transport approach was applied to study the DPM behaviours.High DPM concentration zones were then identified based on the simulation results.The results could provide guidelines for work practices and be helpful to an optimum auxiliary ventilation design to reduce underground miner exposure.
基金Project(2006AA09Z235)supported by the National High Technology Research and Development Program of ChinaProject(CX2009B003)supported by Hunan Provincial Innovation Foundation For Postgraduate,China
文摘Autonomous underwater vehicles （AUVs） navigating on the sea surface are usually required to complete the communication tasks in complex sea conditions. The movement forms and flow field characteristics of a multi-moving state AUV navigating in head sea at high speed were studied. The mathematical model on longitudinal motion of the high-speed AUV in head sea was established with considering the hydrodynamic lift based on strip theory, which was solved to get the heave and pitch of the AUV by Gaussian elimination method. Based on this, computational fluid dynamics （CFD） method was used to establish the mathematical model of the unsteady viscous flow around the AUV with considering free surface effort by using the Reynolds-averaged Navier-Stokes （RANS） equations, shear-stress transport （SST） k-w model and volume of fluid （VOF） model. The three-dimensional numerical wave in the computational field was realized through defining the unsteady inlet boundary condition. The motion forms of the AUV navigating in head sea at high speed were carried out by the program source code of user-defined function （UDF） based on dynamic mesh. The hydrodynamic parameters of the AUV such as drag, lift, pitch torque, velocity, pressure and wave profile were got, which reflect well the real ambient flow field of the AUV navigating in head sea at high speed. The computational wave profile agrees well with the experimental phenomenon of a wave-piercing surface vehicle. The force law of the AUV under the impacts of waves was analyzed qualitatively and quantitatively, which provides an effective theoretical guidance and technical support for the dynamics research and shape design of the AUV in real complex environnaent.
基金Sponsored by the Ministerial Level Foundation (1003)Beijing Municipal Commission of Education Foundation (KM200710016013)
文摘To improve mesh quality for KIVA-3V a method has been developed for rapid mesh generation and dynamic mesh management with moving valves for internal combustion engines. Two phases are included in rapid mesh generation： the initial mesh generation and the mesh pre-treatment. In the second step （pre-treatment）, the connectivity of those cells is generated by a new algorithm added to the KIVA-3V code after the initial mesh generated. In dynamic mesh management phase, a new rezoning algorithm is developed and the basic principle is that the rezoning starts from the moving part. The movement of the adjustment is treated as an ＂earth quake wave＂ propagating to the surrounding vertexes. The amount of coordinate adjustment of the surrounding vertexes is determined by the movement of the epicenter and the distance between the vertexes and the ＂epicenter＂. Finally, a real IC engine mesh is generated and managed aceording to the new method. It gives a new theory and a new method for creating and managing the mesh in IC engine.
文摘Three-dimensional simulations of diesel particulate matter （DPM） distribution inside a single straight entry for the Load-Haul-Dump loader （LHD）-truck loading and truck hauling operations were conducted by using ANSYS FLUENT computational fluid dynamics software. The loading operation was performed for a fixed period of 3 min. The dynamic mesh technique in FLUENT was used to study the impact of truck motion on DPM distribution. The resultant DPM distributions are presented for the cases when the truck were driving upstream and downstream of the loading face. Interesting phenomena were revealed in the study including the piston effect, layering of DPM in the roof region, and backflow of diesel exhaust against ventilation. The results from the simulation can be used to determine if the areas inside the face area and straight entry exceed the current U.S. regulatory requirement for DPM concentration （〉160 pg/m3）. This research can guide the selection of DPM reduction strategies and improve the working practices for the underground miners.
基金Sponsored by the General Programs of the National Natural Science Foundation of China(Grant No.51676058).
文摘Modeling and experiments of the growth of ash deposits during the combustion in a drop tube furnace are presented in this paper.An existing deposition model for ash deposit was used,which involves relationship among the force of gravity,elastic rebound and adhesion forces acting at the moment of ash particle impaction.Important parameters of this model were determined by the experimental data.The influence of particle size and velocity on deposit geometry was investigated.The growth of ash deposits involving various inlet velocities,tube diameters and tube arrangements were simulated with dynamic mesh technique in Fluent.The simulation results show that particle normal impacted velocity plays an important role due to elastic bounce force term.Deposited geometry was influenced by side velocity(velocity of inlet⁃2),and large velocity has an obvious effect to change the shape of deposition and postpone the steady time.
基金Supported by the National Natural Science Foundation of China(No.51005063,51375123)National Science and Technology Cooperation Projects of China(No.2012DFR70840)
文摘A new method is developed to assess and analyze the dynamic performance of hydrostatic bearing oil film by using an amulets-layer dynamic mesh technique. It is implemented using C Language to compile the UDF program of a single oil film of the hydrostatic bearing. The effects of key lubrication parameters of the hydrostatic bearing are evaluated and analyzed under various working conditions,i.e. under no-load,a load of 40 t,a full load of 160 t,and the rotation speed of 1r/min,2r/min,4r/min,8r/min,16r/min,32r/min. The transient data of oil film bearing capacity under different load and rotation speed are acquired for a total of 18 working conditions during the oil film thickness changing. It allows the effective prediction of dynamic performance of large size hydrostatic bearing. Experiments on hydrostatic bearing oil film have been performed and the results were used to define the boundary conditions for the numerical simulations and validate the developed numerical model. The results showed that the oil film thickness became thinner with the increase of the operating time of the hydrostatic bearing,both the oil film rigidity and the oil cavity pressure increased significantly,and the increase of the bearing capacity was inversely proportional to the cube of the change of the film thickness. Meanwhile,the effect of the load condition on carrying capacity of large size static bearing was more important than the speed condition. The error between the simulation value and the experimental value was 4.25%.
基金supported by the Key Program on the S and T for the Pollution Control and Treatment of Water Bodies (Grant Nos.2009ZX07210-008,2009ZX07316-005)the Science and Technology Commission of Shanghai Municipal People's Government (Grant No.072312050)
文摘A Natural Ecological Wastewater Treatment System （NEWTS） is usually built on natural terrain with necessary topography modification to improve water flowing route and pattern, and then the topography modified NEWTS should also have a reasonable water storage volume and hydraulic retention time so as to achieve the anticipated water purification effect. In this study, the dynamic mesh technique based on the finite element method and element storativity coefficients was presented to develop a two-dimensional hydrodynamic and water quality model, which was used to optimize the design of NEWTS under the dynamic land-water boundary due to various water storage volume. The models were employed in the optimized design of NEWTS from a large abandoned coal mine, which purifies the polluted water flowing into a large water storage lake, as part of the East Route South-to-North Water Transfer Project in China. Specifically, the natural topography modification scheme was presented, and further, a reasonable water storage volume and hydraulic residence time was obtained, based on the reasonable estimation of roughness coefficient and pollutant removal rate of the NEWTS with phragmites communis.
基金supported by the Horizontal Subject(Grant No.11471501)the National Basic Research Program of China("973"Project)(Grant No.2013CB228305)
文摘A grate cooler is key equipment in quenching clinker and recovering heat in cement production. A two-dimensional numerical model based on a 5000 t/d cement plant is proposed to for a study on the gas-solid coupled heat transfer process between the cooling air and clinker in a grate cooler. In this study, we use the Fluent dynamic mesh technique and porous media model through which the transient temperature distribution with the clinker motion process and steady temperature and pressure distribution are obtained. We validate the numerical model with the operating data of the cooling air outlet temperature. Then, we discuss the amount of mid-temperature air outlet and average diameter of clinker particles, which affect the heat effective utilization and cooling air pressure drop in clinker layer. We found that after adding one more mid-temperature air outlet, the average temperature of the air flowing into the heat recovery boiler increases by 29.04℃ and the ratio of heat effective utilization increases by 5.3%. This means heat recovery is more effective on adding one more mid-temperature air outlet. Further, the smaller the clinker particles, the more is the pressure drop in clinker layer; thus more power consumption is needed by the cooling fan.
文摘It is of great significance to develop a high-efficiency and low-noise propeller optimization method for new-generation propeller aircraft design.Coupled with free form deformation method,dynamic mesh interpolation technology,optimization algorithm,surrogate model,aerodynamic calculation and aeroacoustic prediction model module,the integrated aerodynamic and aeroacoustic design method of propeller is built.The optimization design for the six-blade propeller is carried out.The non-reduction in efficiency,thrust coefficient and the minimum of aerodynamic noise is treated as the optimization design objective.The spatial vorticity distribution of the propeller before and after the design is also analyzed by using unsteady computational fluid dynamics method.The results show that the optimized propeller can effectively reduce the aerodynamic noise level.The maximum total sound pressure level can be reduced by 5 dB without reducing its aerodynamic performance.The developed method has good application potential in low-noise optimization design of propeller and other rotating machinery.
基金supported by Fundamental Research Funds for the Central UniversitiesChina(Grant No.xjj20100073)Science and Technology Innovation Project of Shaanxi Province of China(Grant No.2011KTCL01-04)
文摘The flow in the positive displacement blower is very complex.The existing two-dimensional numerical simulation cannot provide the detailed flow information,especially flow characteristics along the axial direction,which is unfavorable to improve the performance of positive displacement blower.To investigate the effects of spiral inlet and outlet on the aerodynamic performance of positive displacement blower,three-dimensional unsteady flow characteristics in a three-lobe positive displacement blower with and without the spiral inlet and outlet are simulated by solving Navier-Stokes equations coupled with RNG k-ε turbulent model.In the numerical simulation,the dynamic mesh technique and overset mesh updating method are used.The computational results are compared with the experimental measurements on the variation of flow rate with the outlet pressure to verify the validity of the numerical method presented.The results show that the mass flow rate with the change of pressure is slightly affected by the application of spiral inlet and outlet,but the internal flow state is largely affected.In the exhaust region,the fluctuations of pressure,velocity and temperature as well as the average values of velocity are significantly reduced.This illustrates that the spiral outlet can effectively suppress the fluctuations of pressure,thus reducing reflux shock and energy dissipation.In the intake area,the average value of pressure,velocity and temperature are slightly declined,but the fluctuations of them are significantly reduced,indicating that the spiral inlet plays the role in making the flow more stable.The numerical results obtained reveal the three-dimensional flow characteristics of the positive displacement blower with spiral inlet and outlet,and provide useful reference to improve performance and empirical correction in the noise-reduction design of the positive displacement blowers.
基金partially funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions of ChinaNational Natural Science Foundation of China (No. 11432007)Natural Science Foundation of Jiangsu Province of China (No. BK20140805)
文摘The fluid–structure interaction and aerodynamic shape optimization usually involve the moving or deforming boundaries, thus the dynamic mesh techniques are the key techniques to cope with such deformation. A novel dynamic mesh method was developed based on the Delaunay graph in this paper. According to the Delaunay graph, the mesh points were divided into groups. In each group, a factor ranging from 0 to 1 was calculated based on the area/volume ratio. By introducing a proper function for this factor, this method can control the mesh quality with high efficiency. Several test cases were compared with other dynamic mesh methods regarding mesh quality and CPU time, such as radial basis function method and Delaunay graph mapping method.
基金Project(51276213)supported by the National Natural Science Foundation of ChinaProject(2013BAF01B00)supported by the National Science and Technology Support Program of China
文摘Currently, the flow field of annular seals disturbed by the circular whirl motion of rotors is usually solved using computational fluid dynamics(CFD) to evaluate the five rotordynamic coefficients. The simulations are based on the traditional quasi-steady method. In this work, an improved quasi-steady method along with the transient method was presented to compute the rotordynamic coefficients of a long seal. By comparisons with experimental data, the shortcomings of quasi-steady methods have been identified. Then, the effects of non-uniform incoming flow on seal dynamic coefficients were studied by transient simulations. Results indicate that the long seal has large cross stiffness k and direct mass M which are not good for rotor stability, while the transient method is more suitable for the long seal for its excellent performance in predicting M. When the incoming flow is non-uniform, the stiffness coefficients vary with the eccentric directions. Based on the rotordynamic coefficients under uniform incoming flow, the linearized fluid force formulas, which can consider the effects of non-uniform incoming flow, have been presented and can well explain the varying-stiffness phenomenon.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51061130548 and 51179019)
文摘While passing through a lock, a ship usually undergoes a steady forward motion at low speed. Owing to the size restriction of lock chamber, the shallow water and bank effects on the hydrodynamic forces acting on the ship may be remarkable, which may have an adverse effect on navigation safety. However, the complicated hydrodynamics is not yet fully understood. This paper focuses on the hydrodynamic forces acting on a ship passing through a lock. The unsteady viscous flow and hydrodynamic forces are calculated by applying an unsteady RANS code with a RNG k-e turbulence model. User-defined function （UDF） is compiled to define the ship motion. Meanwhile, the grid regeneration is dealt with by using the dynamic mesh method and sliding interface technique. Numerical study is carried out for a bulk carrier ship passing through the Pierre Vandamme Lock in Zeebrugge at the model scale. The proposed method is validated by comparing the numerical results with the data of captive model tests. By analyzing the numerical results obtained at different speeds, water depths and eccentricities, the influences of speed, water depth and eccentricity on the hydrodynamic forces are illustrated. The numerical method proposed in this paper can qualitatively predict the ship-lock hydrodynamic interaction. It can provide certain guidance on the manoeuvring and control of ships passing through a lock.
基金co-supported by the National Natural Science Foundation of China(No.11672018)the Fundamental Research Funds for the Central Universities,China(No.YWF-23-SDHK-L-002).
文摘This paper focuses on the effect of the phase offset of Leading-Edge(LE)morphing on the aerodynamic characteristics of a pitching NACA0012 airfoil.Assuming an unstretched camber and using polynomial interpolation,an explicit expression for LE nonlinear morphing is proposed and implemented for the large pitching motion of the airfoil.Flow field results and aerodynamic forces are obtained by solving the unsteady Reynolds-averaged Navier-Stokes equations for both the airfoil’s pitching motion and LE morphing.Furthermore,the index of instantaneous aerodynamic power is used to quantify the work done by the airflow in a dynamic process.According to the instantaneous aerodynamic power and energy map,which denotes the energy transfer between the airfoil’s oscillation and flow field,the airfoil is subject to stall flutter.The results show that LE morphing with an optimal phase offset of 315°reduces the energy extraction from the flow field,suppressing the stall flutter instability.This optimal phase offset is effective at different pitching axis positions of the airfoil.The results signify that LE morphing can suppress stall flutter by advancing the occurrence of the first LE vortex and increasing the nose-down moment during the upstroke period.
文摘The performance of a newly designed tri-lobe industrial lobe pump of high capacity is simulated by using commercial CFD solver Ansys Fluent. A combination of user-defined-functions and meshing strategies is employed to capture the rotation of the lobes. The numerical model is validated by comparing the simulated results with the literature values. The processes of suction, displacement, compression and exhaust are accurately captured in the transient simulation. The fluid pressure value remains in the range of inlet pressure value till the processes of suction and displacement are over. The instantaneous process of compression is accurately captured in the simulation. The movement of a particular working chamber is traced along the gradual degree of lobe’s rotation. At five different degrees of lobe’s rotation, pressure contour plots are reported which clearly shows the pressure values inside the working chamber. Each pressure value inside the working chamber conforms to the particular process in which the working chamber is operating. Finally, the power requirement at the shaft of rotation is estimated from the simulated values. The estimated value of power requirement is 3.61 BHP FHP whereas the same calculated theoretically is 3 BHP FHP. The discrepancy is attributed to the assumption of symmetry of blower along the thickness.
基金Project(51221462)supported by the National Natural Science Foundation of China for Innovative Research GroupProject(20120095110001)supported by the Doctoral Fund of Ministry of Education of China+1 种基金Project supported by the Priority Academic Program Development(PAPD)of Jiangsu Higher Education Institutions,ChinaProject(CXJJ201303)supported by the Innovation Foundation of Xuyi Research and Development Center of Mining Equipment and Materials,China University of Mining and Technology,China
文摘A feasible method was proposed to improve the vibration intensity of screen surface via application of a new type elastic screen surface with multi degree of freedom(NTESSMDF). In the NTESSMDF, the primary robs were coupled to the main screen structure with ends embedded into the elastomers, and the secondary robs were attached to adjacent two primary robs with elastic bands. The dynamic model of vibrating screen with NTESSMDF was established based on Lagrange's equation and the equivalent stiffnesses of the elastomer and elastic band were calculated. According to numerical simulation using the 4th order Runge-Kutta method, the vibration intensity of screen surface can be enhanced substantially with an averaged acceleration amplitude increasing ratio of 72.36%. The primary robs and secondary robs vibrate inversely in steady state, which would result in the friability of materials and avoid stoppage. The experimental results validate the dynamic characteristics with acceleration amplitude rising by62.93% on average, which demonstrates the feasibility of NTESSMDF.
文摘Aiming at the high angle of attack pull-up and multi-channel roll pull-up coupling problems of high maneuvering aircraft, this paper establishes the flight attitude control rate by means of unsteady flow numerical solution, dynamic unstructured nested mesh assembly method and numerical solution method of flight mechanics equation. On this basis, a virtual flight simulation platform integrating pneumatics, motion and control is established. Based on this virtual flight simulation platform, F-16 aircraft is simulated by high angle of attack pull-up flight mode and multi-channel roll pull-up coupling flight mode. Finally, the influence of rudder on the yaw control channel is investigated. The results show that the numerical virtual flight simulation platform established in this paper has the ability to simulate maneuvering flight of aircraft.
基金Project(51475232)supported by the National Natural Science Foundation of China
文摘A two-dimensional computational fluid dynamics model was established to simulate the friction stir butt-welding of 6061 aluminum alloy. The dynamic mesh method was applied in this model to make the tool move forward and rotate in a manner similar to a real tool, and the calculated volumetric source of energy was loaded to establish a similar thermal environment to that used in the experiment. Besides, a small piece of zinc stock was embedded into the workpiece as a trace element. Temperature fields and vector plots were determined using a finite volume method, which was indirectly verified by traditional metallography. The simulation result indicated that the temperature distribution was asymmetric but had a similar tendency on the two sides of the welding line. The maximum temperature on the advancing side was approximately 10 K higher than that on the retreating side. Furthermore, the precise process of material flow behavior in combination with streamtraces was demonstrated by contour maps of the phases. Under the shearing force and forward extrusion pressure, material located in front of the tool tended to move along the tangent direction of the rotating tool. Notably, three whirlpools formed under a special pressure environment around the tool, resulting in a uniform composition distribution.
