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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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%.展开更多
The reduced weight and improved efficiency of modern aeronautical structures result in a decreasing separation of frequency ranges of rigid and elastic modes.Particularly,a high-aspect-ratio flexible flying wing is pr...The reduced weight and improved efficiency of modern aeronautical structures result in a decreasing separation of frequency ranges of rigid and elastic modes.Particularly,a high-aspect-ratio flexible flying wing is prone to body freedomflutter(BFF),which is a result of coupling of the rigid body short-periodmodewith 1st wing bendingmode.Accurate prediction of the BFF characteristics is helpful to reflect the attitude changes of the vehicle intuitively and design the active flutter suppression control law.Instead of using the rigid body mode,this work simulates the rigid bodymotion of the model by using the six-degree-of-freedom(6DOF)equation.A dynamicmesh generation strategy particularly suitable for BFF simulation of free flying aircraft is developed.An accurate Computational Fluid Dynamics/Computational Structural Dynamics/six-degree-of-freedom equation(CFD/CSD/6DOF)-based BFF prediction method is proposed.Firstly,the time-domain CFD/CSD method is used to calculate the static equilibrium state of the model.Based on this state,the CFD/CSD/6DOF equation is solved in time domain to evaluate the structural response of themodel.Then combinedwith the variable stiffnessmethod,the critical flutter point of the model is obtained.This method is applied to the BFF calculation of a flyingwing model.The calculation results of the BFF characteristics of the model agree well with those fromthe modalmethod andNastran software.Finally,the method is used to analyze the influence factors of BFF.The analysis results show that the flutter speed can be improved by either releasing plunge constraint or moving the center ofmass forward or increasing the pitch inertia.展开更多
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.展开更多
This paper investigates the hydrodynamic performance of a cylindrical-dual or rectangular-single pontoon floating breakwater using the numerical method and experimental study. The numerical simulation work is based on...This paper investigates the hydrodynamic performance of a cylindrical-dual or rectangular-single pontoon floating breakwater using the numerical method and experimental study. The numerical simulation work is based on the multi-physics computational fluid dynamics(CFD) code and an innovative full-structured dynamic grid method applied to update the three-degree-of-freedom(3-DOF) rigid structure motions. As a time-marching scheme, the trapezoid analogue integral method is used to update the time integration combined with remeshing at each time step.The application of full-structured mesh elements can prevent grids distortion or deformation caused by large-scale movement and improve the stability of calculation. In movable regions, each moving zone is specified with particular motion modes(sway, heave and roll). A series of experimental studies are carried out to validate the performance of the floating body and verify the accuracy of the proposed numerical model. The results are systematically assessed in terms of wave coefficients, mooring line forces, velocity streamlines and the 3-DOF motions of the floating breakwater. When compared with the wave coefficient solutions, excellent agreements are achieved between the computed and experimental data, except in the vicinity of resonant frequency. The velocity streamlines and wave profile movement in the fluid field can also be reproduced using this numerical model.展开更多
Generally the underwater bio-robots take the tail fin as propulsor, and combined with pectoral fin they can manoeuvre agilely and control their position and movement at will. In nature, a lot of fishes realize to susp...Generally the underwater bio-robots take the tail fin as propulsor, and combined with pectoral fin they can manoeuvre agilely and control their position and movement at will. In nature, a lot of fishes realize to suspend itself in water to go forward and to move back up by the pectoral fin moving complexly. So that it is significant theoretically and valuable for practical application to investigate the propulsive principle and hydrodynamic performance of pectoral fin, and find the method utilizing the pectoral fin to manoeuvre the underwater bio-robot agilely. In this paper, a two degree of freedom （DoF） motion model is established for a rigid pectoral fin, and the hydrodynamic performances of the pectoral fin are studied by use of the pectoral fin propulsive experimental platform developed by Harbin Engineering University, simultaneously the hydrodynamic performance of the pectoral fin is analyzed when some parameters change. Then, through the secondary development of FLUENT （CFD code） software, the hydrodynamic performances of rigid pectoral fin in viscous flows are calculated and the results are compared with the latest experimental results. The research in this paper will provide the theoretical reference for the design of the manoeuvring system imitating pectoral fin, at the same time will become the foundation for the development of the small underwater bio-robot.展开更多
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.展开更多
基金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.
基金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.
文摘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.
文摘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.
文摘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.
基金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.
基金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.
基金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.
基金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%.
基金This work was supported by the National Natural Science Foundation of China(No.11872212)and a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions.
文摘The reduced weight and improved efficiency of modern aeronautical structures result in a decreasing separation of frequency ranges of rigid and elastic modes.Particularly,a high-aspect-ratio flexible flying wing is prone to body freedomflutter(BFF),which is a result of coupling of the rigid body short-periodmodewith 1st wing bendingmode.Accurate prediction of the BFF characteristics is helpful to reflect the attitude changes of the vehicle intuitively and design the active flutter suppression control law.Instead of using the rigid body mode,this work simulates the rigid bodymotion of the model by using the six-degree-of-freedom(6DOF)equation.A dynamicmesh generation strategy particularly suitable for BFF simulation of free flying aircraft is developed.An accurate Computational Fluid Dynamics/Computational Structural Dynamics/six-degree-of-freedom equation(CFD/CSD/6DOF)-based BFF prediction method is proposed.Firstly,the time-domain CFD/CSD method is used to calculate the static equilibrium state of the model.Based on this state,the CFD/CSD/6DOF equation is solved in time domain to evaluate the structural response of themodel.Then combinedwith the variable stiffnessmethod,the critical flutter point of the model is obtained.This method is applied to the BFF calculation of a flyingwing model.The calculation results of the BFF characteristics of the model agree well with those fromthe modalmethod andNastran software.Finally,the method is used to analyze the influence factors of BFF.The analysis results show that the flutter speed can be improved by either releasing plunge constraint or moving the center ofmass forward or increasing the pitch inertia.
基金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.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51579122,51609109,and 51622902)the Natural Science Found of Jiangsu Province(Grant No.BK20160556)+1 种基金the University Natural Science Research Project of Jiangsu Province(Grant No.16kjb70003)the Key Lab Foundation for Advanced Manufacturing Technology of Jiangsu Province(Grant No.CJ1506)
文摘This paper investigates the hydrodynamic performance of a cylindrical-dual or rectangular-single pontoon floating breakwater using the numerical method and experimental study. The numerical simulation work is based on the multi-physics computational fluid dynamics(CFD) code and an innovative full-structured dynamic grid method applied to update the three-degree-of-freedom(3-DOF) rigid structure motions. As a time-marching scheme, the trapezoid analogue integral method is used to update the time integration combined with remeshing at each time step.The application of full-structured mesh elements can prevent grids distortion or deformation caused by large-scale movement and improve the stability of calculation. In movable regions, each moving zone is specified with particular motion modes(sway, heave and roll). A series of experimental studies are carried out to validate the performance of the floating body and verify the accuracy of the proposed numerical model. The results are systematically assessed in terms of wave coefficients, mooring line forces, velocity streamlines and the 3-DOF motions of the floating breakwater. When compared with the wave coefficient solutions, excellent agreements are achieved between the computed and experimental data, except in the vicinity of resonant frequency. The velocity streamlines and wave profile movement in the fluid field can also be reproduced using this numerical model.
基金supported by the National Natural Science Foundation of China (Grant Nos .50579007 and 50879014)the Specialized Research Fund for the Doctoral Program of Higher Education (Grant No.200802170010)
文摘Generally the underwater bio-robots take the tail fin as propulsor, and combined with pectoral fin they can manoeuvre agilely and control their position and movement at will. In nature, a lot of fishes realize to suspend itself in water to go forward and to move back up by the pectoral fin moving complexly. So that it is significant theoretically and valuable for practical application to investigate the propulsive principle and hydrodynamic performance of pectoral fin, and find the method utilizing the pectoral fin to manoeuvre the underwater bio-robot agilely. In this paper, a two degree of freedom （DoF） motion model is established for a rigid pectoral fin, and the hydrodynamic performances of the pectoral fin are studied by use of the pectoral fin propulsive experimental platform developed by Harbin Engineering University, simultaneously the hydrodynamic performance of the pectoral fin is analyzed when some parameters change. Then, through the secondary development of FLUENT （CFD code） software, the hydrodynamic performances of rigid pectoral fin in viscous flows are calculated and the results are compared with the latest experimental results. The research in this paper will provide the theoretical reference for the design of the manoeuvring system imitating pectoral fin, at the same time will become the foundation for the development of the small underwater bio-robot.
基金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.
