The forming defects, including thinning, rupture, wrinkling and springback, usually arising in producing a side-door impact beam, were investigated by trial and numerical simulation. A temperature-related constitutive...The forming defects, including thinning, rupture, wrinkling and springback, usually arising in producing a side-door impact beam, were investigated by trial and numerical simulation. A temperature-related constitutive model specific to the temperature range from 350 °C to 500 °C was established and used for the numerical simulation. The trial and numerical simulation were conducted to clarify the quantitative characteristics of forming defects and to analyze the effects of process parameters on the forming defects. Results show that the rupture situation is ameliorated and the springback is eliminated in the aluminum alloy hot stamping. The wrinkling severity decreases with increasing blank holder force (BHF), but the BHF greater than 15 kN causes the rupture at the deepest drawing position of workpiece. The forming defects are avoided with lubricant in the feasible ranges of process parameters: the BHF of 3 to 5 kN and the stamping speed of 50 to 200 mm/s.展开更多
A VBHF(Variable Blank Holder Force) optimization strategy was employed to determine the optimal time-variable and spatial-variable BHF trajectories,aiming at improving the formability of automobile panels with aluminu...A VBHF(Variable Blank Holder Force) optimization strategy was employed to determine the optimal time-variable and spatial-variable BHF trajectories,aiming at improving the formability of automobile panels with aluminum alloy sheet.The strategy was implemented based on adaptive simulation to calculate the critical wrinkling BHF for each segmented binder of the Numisheet' 05 deck lid in a single round of simulation.The thickness comparison of the stamped part under optimal VBHF and constant BHF shows that the variance of the four sections is decreased by 70%,44%,64% and 61%,respectively,which indicates significant improvement in thickness distribution and variation control.The investigation through strain path comparison reveals the fundamental reason of formability improvement.The study proves the applicability of the new VBHF optimization strategy to complex parts with aluminum alloy sheet.展开更多
The idea that the collapse proceeds from the outer boundary of the cavity cloud towards its center for the ultrasonic cavitation proposed by Hasson and Morch in 1980s is further developed for calculating the collapse ...The idea that the collapse proceeds from the outer boundary of the cavity cloud towards its center for the ultrasonic cavitation proposed by Hasson and Morch in 1980s is further developed for calculating the collapse pressure and boundaries of cavity cloud at the collapse stage of bubbles for hydraulic cavitation flow in Venturi in present research. The numerical simulation is carried out based on Gilmore's eouations of bubble dynamics, which take account of the compressibility of fluid besides the viscosity and interfacial tension. The collapse of the cavity cloud is considered to proceed layer by layer from the outer cloud towards its inner part. The simulation results indicate that thepredicted boundaries of the cavity cloudat the collapse stage agree.well with the exPerimental ones.It is also found that the maximum collapse pressure of the cavity cloud is several times as high as the collapse pressure of outside boundary, and it is located at a point in the axis, where the cavity cloud disappears completely. This means that a cavity cloud has higher collapse pressure or strength than that of a single bubble due to the interactions of the bubbles. The effects of operation and structural parameters on the collapse pressure are also analyzed in detail.展开更多
In this paper, a numerical study of flow in the turbulence boundary layer with adverse and pressure gradients (APGs) is conducted by using Reynolds-averaged Navier-Stokes (RANS) equations. This research chooses si...In this paper, a numerical study of flow in the turbulence boundary layer with adverse and pressure gradients (APGs) is conducted by using Reynolds-averaged Navier-Stokes (RANS) equations. This research chooses six typical turbulence models, which are critical to the computing precision, and to evaluating the issue of APGs. Local frictional resistance coefficient is compared between numerical and experimental results. The same comparisons of dimensionless averaged velocity profiles are also performed. It is found that results generated by Wilcox (2006) k-co are most close to the experimental data. Meanwhile, turbulent quantities such as turbulent kinetic energy and Reynolds-stress are also studied.展开更多
Competition of multiple Gortler modes in hypersonic boundary layer flows are investigated with the local and marching methods. The wall-layer mode (mode W) and the trapped-layer mode (mode T) both occur in the com...Competition of multiple Gortler modes in hypersonic boundary layer flows are investigated with the local and marching methods. The wall-layer mode (mode W) and the trapped-layer mode (mode T) both occur in the compressible boundary layer where there exists a temperature adjustment layer near the upper edge. The mode T has the largest growth rate at a lower Gortler number while the mode W dominates at larger G/Srtler numbers. These two modes are both responsible for the flow transition in the hypersonic flows especially when Gortler number is in the high value range in which the crossover of these two modes takes place. Such high Gortler numbers are virtually far beyond the neutral regime. The nonparallel base flows, therefore, cease to influence the stability behavior of the Gortler modes. The effects of the Mach number on the multiple Gortler modes are studied within a chosen Mach number of 0.95, 2, 4 and 6. When the flow Mach number is sufficiently large, e.g., Ma ≥4, the growth rate crossover of the mode T and mode W occurs both in the conventional G-β map as well as on the route downstream for a fixed wavelength disturbance. Four particular regions (Region T, T-W, W-T and W) around the crossover point are highlighted with the marching analysis and the result matches that of the local analysis. The initial disturbance of a normal mode maintains the shape in its corresponding dominating region while a shape-transformation occurs outside this region.展开更多
In order to develop a wall function boundary condition for high-speed flows so as to reduce the grid-dependence of the simula- tion for the skin friction and heat flux, a research was performed to improve the compress...In order to develop a wall function boundary condition for high-speed flows so as to reduce the grid-dependence of the simula- tion for the skin friction and heat flux, a research was performed to improve the compressible wall function boundary condition proposed by Nichols. Values of parameters in the velocity law-of-the-wall were revised according to numerical experiments and the expression of temperature law-of-the-wall was modified based on theoretical analysis and numerical simulation. Be- sides, the formula of the heat conduction term in near-wall region was derived so that the coupling between the wall function boundary condition and CFD code was realized more accurately. Whereafter, the application study of the modified wall func- tion was carried out. The numerical case of supersonic turbulent boundary layer on a flat plate illustrated that the modified wall function produces reasonable results of skin friction and heat flux, and profiles of velocity, temperature and turbulent eddy viscosity for coarse grids with the initial wall spacing of y+〈400, and that the modifications to the original wall function can obviously improve the simulation precision. As for the application of separation flows, it was found from the numerical cases of supersonic cavity flow and hypersonic axisymmetric compression comer that the compressible velocity law-of-the-wall originally established based on the fully-developed attached turbulent boundary layer approximately holds in the near-wall re- gion inside the separation flows, which ensures that reliable skin friction and heat flux can be given by the wall function inside the separation flows, while for the region near separation and reattachment points, the wall function gives results with a rela- tively large error, because the velocity law-of-the-wall used in the wall function takes on obvious deviation from the real ve- locity profiles near the separation and reattachment points.展开更多
Gas diffusion layer(GDL) plays a great important role in proton exchange membrane fuel cell(PEMFC).Water transport mechanism in GDL is still not clear.In the present study,an ex-situ transparent setup is built to visu...Gas diffusion layer(GDL) plays a great important role in proton exchange membrane fuel cell(PEMFC).Water transport mechanism in GDL is still not clear.In the present study,an ex-situ transparent setup is built to visualize the transport phenomena and to measure the threshold pressure of water in GDL at different temperatures.It is found that the relationship between the breakthrough pressure and the temperature is nearly linear(i.e.the pressure decreases linearly with the increase of temperature).To avoid the problems faced by the continuum models,the pore network model is developed to simulate the liquid water transport through the carbon paper.A uniform pressure boundary condition is used in simulation and the results are similar to the ones obtained in the experiment.The reason is that the contact angle and surface tension coefficient of water in GDLs change accordingly with the change of temperature.展开更多
The nonlinear evolution of a finite-amplitude disturbance in a 3-D supersonic boundary layer over a cone was investigated recently by Liu et al. using direct numerical simulation (DNS). It was found that certain sma...The nonlinear evolution of a finite-amplitude disturbance in a 3-D supersonic boundary layer over a cone was investigated recently by Liu et al. using direct numerical simulation (DNS). It was found that certain small-scale 3-D disturbances amplified rapidly. These disturbances exhibit the characteristics of second modes, and the most amplified components have a well- defined spanwise wavelength, indicating a clear selectivity of the amplification. In the case of a cone, the three-dimensionality of the base flow and the disturbances themselves may be responsible for the rapid amplification. In order to ascertain which of these two effects are essential, in this study we carried out DNS of the nonlinear evolution of a spanwise localized disturbance (wavepacket) in a flat-plate boundary layer. A similar amplification of small-scale disturbances was observed, suggesting that the direct reason for the rapid amplification is the three-dimensionality of the disturbances rather than the three-dimensional nature of the base flow, even though the latter does alter the spanwise distribution of the disturbance. The rapid growth of 3-D waves may be attributed to the secondary instability mechanism. Further simulations were performed for a wavepacket of first modes in a supersonic boundary layer and of Tollmien-Schlichting (T-S) waves in an incompressible boundary layer. The re- suits show that the amplifying components are in the band centered at zero spanwise wavenumber rather than at a finite spanwise wavenumber. It is therefore concluded that the rapid growth of 3-D disturbances in a band centered at a preferred large spanwise wavenumber is the main characteristic of nonlinear evolution of second mode disturbances in supersonic boundary layers.展开更多
In this paper, we have discussed the food movement in stomach with thermal bound- ary conditions. Eyring-Prandtl fluid model is considered. Formulation of the considered phenomena have been developed for both fixed an...In this paper, we have discussed the food movement in stomach with thermal bound- ary conditions. Eyring-Prandtl fluid model is considered. Formulation of the considered phenomena have been developed for both fixed and moving frame of references. Regular perturbation is used to find the solution of stream function, temperature profile and pressure gradient. Analysis has been carried out for velocity, "stream function, temper- ature, pressure gradient and heat transfer". Appearance of pressure gradient is quite complicated so to get the expression for pressure rise we have used numerical integra- tion. It is perceived that the velocity close to the channel walls is not same in outlook of the Eyrin^Prandtl fluid parameter taken as fl and Hartman number M. The velocity decreases by increasing β and M.展开更多
The first part of this study is focused on the numerical modelling and experimental investigation of transonic flow through a 2D model of the male rotor-housing gap in a dry screw compressor.Numerical simulations of t...The first part of this study is focused on the numerical modelling and experimental investigation of transonic flow through a 2D model of the male rotor-housing gap in a dry screw compressor.Numerical simulations of the clearance flow are performed with the help of the in-house compressible Navier-Stokes solver.Experimental measurements based on the Schlieren method in Toepler configuration are carried out.The objective of the second part of the study is to derive the analytical solution of gas microflow development in a gap between two parallel plates.The microflow is assumed to be laminar,incompressible and the velocity slip boundary conditions are considered at the walls.The constant velocity profile is prescribed at the inlet.For the mathematical description of the problem,the Oseen equation is used.The analytical results are compared with the numerical ones obtained using the developed incompressible Navier-Stokes solver including the slip flow boundary conditions.展开更多
The effects of Reynolds number on the performance of a high pressure-ratio turbocharger compressor were investigated by both experiments and numerical simulation. The experimental results show that the pressure ratio ...The effects of Reynolds number on the performance of a high pressure-ratio turbocharger compressor were investigated by both experiments and numerical simulation. The experimental results show that the pressure ratio and the efficiency of the compressor respectively decrease by 7.9% and 6.9% when Reynolds number drops from 9.86×10 5 to 2.96×10 5 . The numerical simulation predicts a similar trend as the experimental results although it underestimates the deterioration of the performance under low Reynolds number conditions. According to simulation results, the boundary layer thickness increases at the inducer, which decreases the throat area and leads to smaller choke mass flow rate. The increments of the boundary thickness are relatively small at the rear part of the impeller. The boundary layer separation flow is severe. The interaction between boundary layer separation flows and leakage flows causes the high loss region at the rear part of the impeller passage under low Reynolds number condition.展开更多
基金Project(P2014-15)supported by the State Key Laboratory of Materials Processing and Die&Mould Technology,Huazhong University of Science and Technology,ChinaProject supported by the Beijing Laboratory of Metallic Materials and Processing for Modern Transportation,China
文摘The forming defects, including thinning, rupture, wrinkling and springback, usually arising in producing a side-door impact beam, were investigated by trial and numerical simulation. A temperature-related constitutive model specific to the temperature range from 350 °C to 500 °C was established and used for the numerical simulation. The trial and numerical simulation were conducted to clarify the quantitative characteristics of forming defects and to analyze the effects of process parameters on the forming defects. Results show that the rupture situation is ameliorated and the springback is eliminated in the aluminum alloy hot stamping. The wrinkling severity decreases with increasing blank holder force (BHF), but the BHF greater than 15 kN causes the rupture at the deepest drawing position of workpiece. The forming defects are avoided with lubricant in the feasible ranges of process parameters: the BHF of 3 to 5 kN and the stamping speed of 50 to 200 mm/s.
基金Project(50934011) supported by the National Natural Science Foundation of ChinaProject(20080430085) supported by the China Postdoctoral Science Foundation
文摘A VBHF(Variable Blank Holder Force) optimization strategy was employed to determine the optimal time-variable and spatial-variable BHF trajectories,aiming at improving the formability of automobile panels with aluminum alloy sheet.The strategy was implemented based on adaptive simulation to calculate the critical wrinkling BHF for each segmented binder of the Numisheet' 05 deck lid in a single round of simulation.The thickness comparison of the stamped part under optimal VBHF and constant BHF shows that the variance of the four sections is decreased by 70%,44%,64% and 61%,respectively,which indicates significant improvement in thickness distribution and variation control.The investigation through strain path comparison reveals the fundamental reason of formability improvement.The study proves the applicability of the new VBHF optimization strategy to complex parts with aluminum alloy sheet.
基金Supported by the National Natural Science Foundation of China (10472024).
文摘The idea that the collapse proceeds from the outer boundary of the cavity cloud towards its center for the ultrasonic cavitation proposed by Hasson and Morch in 1980s is further developed for calculating the collapse pressure and boundaries of cavity cloud at the collapse stage of bubbles for hydraulic cavitation flow in Venturi in present research. The numerical simulation is carried out based on Gilmore's eouations of bubble dynamics, which take account of the compressibility of fluid besides the viscosity and interfacial tension. The collapse of the cavity cloud is considered to proceed layer by layer from the outer cloud towards its inner part. The simulation results indicate that thepredicted boundaries of the cavity cloudat the collapse stage agree.well with the exPerimental ones.It is also found that the maximum collapse pressure of the cavity cloud is several times as high as the collapse pressure of outside boundary, and it is located at a point in the axis, where the cavity cloud disappears completely. This means that a cavity cloud has higher collapse pressure or strength than that of a single bubble due to the interactions of the bubbles. The effects of operation and structural parameters on the collapse pressure are also analyzed in detail.
基金Foundation item: Supported by the National Natural Science Foundation of China (Nos.51309040, 51379033, 51209027, 51309025), Open Research Fund of State Key Laboratory of Ocean Engineering (Shanghai Jiao Tong University) (Grant No.1402), and Fundamental Research Fund for the Central Universities (DMU3132015089).
文摘In this paper, a numerical study of flow in the turbulence boundary layer with adverse and pressure gradients (APGs) is conducted by using Reynolds-averaged Navier-Stokes (RANS) equations. This research chooses six typical turbulence models, which are critical to the computing precision, and to evaluating the issue of APGs. Local frictional resistance coefficient is compared between numerical and experimental results. The same comparisons of dimensionless averaged velocity profiles are also performed. It is found that results generated by Wilcox (2006) k-co are most close to the experimental data. Meanwhile, turbulent quantities such as turbulent kinetic energy and Reynolds-stress are also studied.
基金supported by the National Natural Science Foundation of China(Grant Nos.10932005 and 11202115)
文摘Competition of multiple Gortler modes in hypersonic boundary layer flows are investigated with the local and marching methods. The wall-layer mode (mode W) and the trapped-layer mode (mode T) both occur in the compressible boundary layer where there exists a temperature adjustment layer near the upper edge. The mode T has the largest growth rate at a lower Gortler number while the mode W dominates at larger G/Srtler numbers. These two modes are both responsible for the flow transition in the hypersonic flows especially when Gortler number is in the high value range in which the crossover of these two modes takes place. Such high Gortler numbers are virtually far beyond the neutral regime. The nonparallel base flows, therefore, cease to influence the stability behavior of the Gortler modes. The effects of the Mach number on the multiple Gortler modes are studied within a chosen Mach number of 0.95, 2, 4 and 6. When the flow Mach number is sufficiently large, e.g., Ma ≥4, the growth rate crossover of the mode T and mode W occurs both in the conventional G-β map as well as on the route downstream for a fixed wavelength disturbance. Four particular regions (Region T, T-W, W-T and W) around the crossover point are highlighted with the marching analysis and the result matches that of the local analysis. The initial disturbance of a normal mode maintains the shape in its corresponding dominating region while a shape-transformation occurs outside this region.
基金supported by the National Natural Science Foundation of China(Grant No.11202014)
文摘In order to develop a wall function boundary condition for high-speed flows so as to reduce the grid-dependence of the simula- tion for the skin friction and heat flux, a research was performed to improve the compressible wall function boundary condition proposed by Nichols. Values of parameters in the velocity law-of-the-wall were revised according to numerical experiments and the expression of temperature law-of-the-wall was modified based on theoretical analysis and numerical simulation. Be- sides, the formula of the heat conduction term in near-wall region was derived so that the coupling between the wall function boundary condition and CFD code was realized more accurately. Whereafter, the application study of the modified wall func- tion was carried out. The numerical case of supersonic turbulent boundary layer on a flat plate illustrated that the modified wall function produces reasonable results of skin friction and heat flux, and profiles of velocity, temperature and turbulent eddy viscosity for coarse grids with the initial wall spacing of y+〈400, and that the modifications to the original wall function can obviously improve the simulation precision. As for the application of separation flows, it was found from the numerical cases of supersonic cavity flow and hypersonic axisymmetric compression comer that the compressible velocity law-of-the-wall originally established based on the fully-developed attached turbulent boundary layer approximately holds in the near-wall re- gion inside the separation flows, which ensures that reliable skin friction and heat flux can be given by the wall function inside the separation flows, while for the region near separation and reattachment points, the wall function gives results with a rela- tively large error, because the velocity law-of-the-wall used in the wall function takes on obvious deviation from the real ve- locity profiles near the separation and reattachment points.
基金supported by the National Natural Science Foundation of China (No.50976011)Fundamental Research Funds for the Central Universities of China (No. 2009JBM090)
文摘Gas diffusion layer(GDL) plays a great important role in proton exchange membrane fuel cell(PEMFC).Water transport mechanism in GDL is still not clear.In the present study,an ex-situ transparent setup is built to visualize the transport phenomena and to measure the threshold pressure of water in GDL at different temperatures.It is found that the relationship between the breakthrough pressure and the temperature is nearly linear(i.e.the pressure decreases linearly with the increase of temperature).To avoid the problems faced by the continuum models,the pore network model is developed to simulate the liquid water transport through the carbon paper.A uniform pressure boundary condition is used in simulation and the results are similar to the ones obtained in the experiment.The reason is that the contact angle and surface tension coefficient of water in GDLs change accordingly with the change of temperature.
基金supported by the National Basic Research Program of China (Grant No. 2009CB724103)
文摘The nonlinear evolution of a finite-amplitude disturbance in a 3-D supersonic boundary layer over a cone was investigated recently by Liu et al. using direct numerical simulation (DNS). It was found that certain small-scale 3-D disturbances amplified rapidly. These disturbances exhibit the characteristics of second modes, and the most amplified components have a well- defined spanwise wavelength, indicating a clear selectivity of the amplification. In the case of a cone, the three-dimensionality of the base flow and the disturbances themselves may be responsible for the rapid amplification. In order to ascertain which of these two effects are essential, in this study we carried out DNS of the nonlinear evolution of a spanwise localized disturbance (wavepacket) in a flat-plate boundary layer. A similar amplification of small-scale disturbances was observed, suggesting that the direct reason for the rapid amplification is the three-dimensionality of the disturbances rather than the three-dimensional nature of the base flow, even though the latter does alter the spanwise distribution of the disturbance. The rapid growth of 3-D waves may be attributed to the secondary instability mechanism. Further simulations were performed for a wavepacket of first modes in a supersonic boundary layer and of Tollmien-Schlichting (T-S) waves in an incompressible boundary layer. The re- suits show that the amplifying components are in the band centered at zero spanwise wavenumber rather than at a finite spanwise wavenumber. It is therefore concluded that the rapid growth of 3-D disturbances in a band centered at a preferred large spanwise wavenumber is the main characteristic of nonlinear evolution of second mode disturbances in supersonic boundary layers.
文摘In this paper, we have discussed the food movement in stomach with thermal bound- ary conditions. Eyring-Prandtl fluid model is considered. Formulation of the considered phenomena have been developed for both fixed and moving frame of references. Regular perturbation is used to find the solution of stream function, temperature profile and pressure gradient. Analysis has been carried out for velocity, "stream function, temper- ature, pressure gradient and heat transfer". Appearance of pressure gradient is quite complicated so to get the expression for pressure rise we have used numerical integra- tion. It is perceived that the velocity close to the channel walls is not same in outlook of the Eyrin^Prandtl fluid parameter taken as fl and Hartman number M. The velocity decreases by increasing β and M.
基金supported by the grant GACR 101/08/0623 of the Czech Science Foundation and by the research project MSM4977751303 of the Ministry of Education,Youth and Sports of the Czech Republic
文摘The first part of this study is focused on the numerical modelling and experimental investigation of transonic flow through a 2D model of the male rotor-housing gap in a dry screw compressor.Numerical simulations of the clearance flow are performed with the help of the in-house compressible Navier-Stokes solver.Experimental measurements based on the Schlieren method in Toepler configuration are carried out.The objective of the second part of the study is to derive the analytical solution of gas microflow development in a gap between two parallel plates.The microflow is assumed to be laminar,incompressible and the velocity slip boundary conditions are considered at the walls.The constant velocity profile is prescribed at the inlet.For the mathematical description of the problem,the Oseen equation is used.The analytical results are compared with the numerical ones obtained using the developed incompressible Navier-Stokes solver including the slip flow boundary conditions.
基金supported by the National Natural Science Foundation of China (Grant No. 51176087)
文摘The effects of Reynolds number on the performance of a high pressure-ratio turbocharger compressor were investigated by both experiments and numerical simulation. The experimental results show that the pressure ratio and the efficiency of the compressor respectively decrease by 7.9% and 6.9% when Reynolds number drops from 9.86×10 5 to 2.96×10 5 . The numerical simulation predicts a similar trend as the experimental results although it underestimates the deterioration of the performance under low Reynolds number conditions. According to simulation results, the boundary layer thickness increases at the inducer, which decreases the throat area and leads to smaller choke mass flow rate. The increments of the boundary thickness are relatively small at the rear part of the impeller. The boundary layer separation flow is severe. The interaction between boundary layer separation flows and leakage flows causes the high loss region at the rear part of the impeller passage under low Reynolds number condition.