Why the stall of an airfoil can be significantly delayed by its pitching-up motion？ Various attempts have been proposed to answer this question over the past half century, but none is satisfactory. In this letter we ...Why the stall of an airfoil can be significantly delayed by its pitching-up motion？ Various attempts have been proposed to answer this question over the past half century, but none is satisfactory. In this letter we prove that a chain of vorticity-dynamics processes at accelerating boundary is fully responsible for the causal mechanism underlying this peculiar phenomenon. The local flow behavior is explained by a simple potential-flow model.展开更多
The dissipation function in turbulent plane Poiseuille flows(PPFs) and plane Couette flows(PCFs) subject to spanwise rotations is analyzed. It is found that, in the PCFs without system rotations, the mean part is cons...The dissipation function in turbulent plane Poiseuille flows(PPFs) and plane Couette flows(PCFs) subject to spanwise rotations is analyzed. It is found that, in the PCFs without system rotations, the mean part is constant while the fluctuation part follows a logarithmic law, resulting in a similar logarithmic skin friction law as PPFs.However, if the flow system rotates in the spanwise direction, no obvious dependence on the rotation number can be evaluated. In the PPFs with rotations, the dissipation function shows an increase with the rotation number, while in the PCFs with rotations,when the rotation number increases, the dissipation function first decreases and then increases.展开更多
In this paper, we consider L<sup>∞</sup> estimates of eigenfunction, or more generally, the L<sup>∞</sup> estimates of equation -Δu=fu. We use heat flow to give a new proof of the L<sup&...In this paper, we consider L<sup>∞</sup> estimates of eigenfunction, or more generally, the L<sup>∞</sup> estimates of equation -Δu=fu. We use heat flow to give a new proof of the L<sup>∞</sup> estimates for such type equations.展开更多
Direct numerical simulation of rotating plane Couette flow(RPCF) at Re_w= 1300 and Ro = 0.02 was performed with different mesh resolutions and different sizes of computation domain. Our results showed that a grid reso...Direct numerical simulation of rotating plane Couette flow(RPCF) at Re_w= 1300 and Ro = 0.02 was performed with different mesh resolutions and different sizes of computation domain. Our results showed that a grid resolution in wall units with ?x^+=8.51, ?z^+= 4.26, ?y^+|_(min)= 0.0873 and ?y^+|_(max)= 3.89 is fine enough to simulate the problem at the present parameters. The streamwise length Lxand spanwise length Lzof the computational box have different impacts on the flow statistics, where the statistics were converged if Lxis longer than 8πh, while no converged results were obtained for different Lz. More importantly,our results with very long simulation time showed that a state transition would happen if L_x≥ 8πh, from a state with four pairs of roll cells to a state with three pairs of roll cells with L_z= 6πh. Each state could survive for more than 1500 h/U_w, and the flow statistics were different.展开更多
This paper presents an extension work of the hybrid scheme proposed by Wang et al.[J.Comput.Phys.229(2010)169-180]for numerical simulation of sub-sonic isotropic turbulence to supersonic turbulence regime.The scheme s...This paper presents an extension work of the hybrid scheme proposed by Wang et al.[J.Comput.Phys.229(2010)169-180]for numerical simulation of sub-sonic isotropic turbulence to supersonic turbulence regime.The scheme still utilizes an 8th-order compact scheme with built-in hyperviscosity for smooth regions and a 7th-order WENO scheme for highly compression regions,but now both in their con-servation formulations and for the latter with the Roe type characteristic-wise recon-struction.To enhance the robustness of the WENO scheme without compromising its high-resolution and accuracy,the recursive-order-reduction procedure is adopted,where a new type of reconstruction-failure-detection criterion is constructed from the idea of positivity-preserving.In addition,a new form of cooling function is proposed,which is proved also to be positivity-preserving.With a combination of these techniques,the new scheme not only inherits the good properties of the original one but also extends largely the computable range of turbulent Mach number,which has been further confirmed by numerical results.展开更多
Compressible flow past a circular cylinder at an inflow Reynolds number of 2×105 is numerically investigated by using a constrained large-eddy simulation(CLES)technique.Numerical simulation with adiabatic wall bo...Compressible flow past a circular cylinder at an inflow Reynolds number of 2×105 is numerically investigated by using a constrained large-eddy simulation(CLES)technique.Numerical simulation with adiabatic wall boundary condition and at a free-stream Mach number of 0.75 is conducted to validate and verify the performance of the present CLES method in predicting separated flows.Some typical and characteristic physical quantities,such as the drag coefficient,the root-mean-square lift fluctuations,the Strouhal number,the pressure and skin friction distributions around the cylinder,etc.are calculated and compared with previously reported experimental data,finer-grid large-eddy simulation(LES)data and those obtained in the present LES and detached-eddy simulation(DES)on coarse grids.It turns out that CLES is superior to DES in predicting such separated flow and that CLES can mimic the intricate shock wave dynamics quite well.Then,the effects of Mach number on the flow patterns and parameters such as the pressure,skin friction and drag coefficients,and the cylinder surface temperature are studied,with Mach number varying from 0.1 to 0.95.Nonmonotonic behaviors of the pressure and skin friction distributions are observed with increasing Mach number and the minimum mean separation angle occurs at a subcritical Mach number of between 0.3 and 0.5.Additionally,the wall temperature effects on the thermodynamic and aerodynamic quantities are explored in a series of simulations using isothermal wall boundary conditions at three different wall temperatures.It is found that the flow separates earlier from the cylinder surface with a longer recirculation length in the wake and a higher pressure coefficient at the rear stagnation point for higher wall temperature.Moreover,the influences of different thermal wall boundary conditions on the flow field are gradually magnified from the front stagnation point to the rear stagnation point.It is inferred that the CLES approach in its current version is a useful and effective tool for simulating wall-bounded compressible turbulent flows with massive separations.展开更多
In this paper, we perform a numerical simulation of the cavitating flow around an underwater hemispherical-head slender body running at a high speed. For the first time, the slip boundary condition is introduced into ...In this paper, we perform a numerical simulation of the cavitating flow around an underwater hemispherical-head slender body running at a high speed. For the first time, the slip boundary condition is introduced into this problem, and we find that the slip boundary condition has a big influence on the cavitation in the flow-separation zone. By simulating the cavitating flow under different cavitation numbers, we demonstrate that the slip boundary condition can effectively reduce the intensity of cavitation, as represented by the length of cavitation bubbles. The present paper provides a new method for utilization of new surface materials to control the cavitation on the underwater moving objects.展开更多
Using the incompressible isotropic turbulent fields obtained from direct numerical simulation and large-eddy simulation,we studied the statistics of oscillation structures based on local zero-crossings and their relat...Using the incompressible isotropic turbulent fields obtained from direct numerical simulation and large-eddy simulation,we studied the statistics of oscillation structures based on local zero-crossings and their relation with inertial-range intermittency for transverse velocity and passive scalar.Our results show that for both the velocity and passive scalar,the local oscillation structures are statistically scaleinvariant at high Reynolds number,and the inertial-range intermittency of the overall flow region is determined by the most intermittent structures characterized by one local zero-crossing.Local flow patterns conditioned on the oscillation structures are characterized by the joint probability density function of the invariants of the filtered velocity gradient tensor at inertial range.We demonstrate that the most intermittent regions for longitudinal velocity tend to lay at the saddle area,while those for the transverse velocity tend to locate at the vortex-dominated area.The connection between the ramp-cliff structures in passive scalar field and the corresponding saddle regions in the velocity field is also verified by the approach of oscillation structure classification.展开更多
基金supported by the National Natural Science Foundation of China(10921202,11221062,11521091,and 11472016)
文摘Why the stall of an airfoil can be significantly delayed by its pitching-up motion？ Various attempts have been proposed to answer this question over the past half century, but none is satisfactory. In this letter we prove that a chain of vorticity-dynamics processes at accelerating boundary is fully responsible for the causal mechanism underlying this peculiar phenomenon. The local flow behavior is explained by a simple potential-flow model.
基金Project supported by the National Natural Science Foundation of China(Nos.11772297 and11822208)
文摘The dissipation function in turbulent plane Poiseuille flows(PPFs) and plane Couette flows(PCFs) subject to spanwise rotations is analyzed. It is found that, in the PCFs without system rotations, the mean part is constant while the fluctuation part follows a logarithmic law, resulting in a similar logarithmic skin friction law as PPFs.However, if the flow system rotates in the spanwise direction, no obvious dependence on the rotation number can be evaluated. In the PPFs with rotations, the dissipation function shows an increase with the rotation number, while in the PCFs with rotations,when the rotation number increases, the dissipation function first decreases and then increases.
文摘In this paper, we consider L<sup>∞</sup> estimates of eigenfunction, or more generally, the L<sup>∞</sup> estimates of equation -Δu=fu. We use heat flow to give a new proof of the L<sup>∞</sup> estimates for such type equations.
基金supported by the National Natural Science Foundation of China(Grant Nos.12225204 and 91752201)Department of Science and Technology of Guangdong Province(Grant No.2020B1212030001)+2 种基金the Key-Area Research and Development Program of Guangdong Province(Grant No.2021B0101190003)the Shenzhen Science and Technology Program(Grant No.KQTD20180411143441009)Numerical simulations have been supported by the Center for Computational Science and Engineering of Southern University of Science and Technology.
基金supported by the National Natural Science Foundation of China(Grant Nos.11822208,11772297,11672123,and 91752201)
文摘Direct numerical simulation of rotating plane Couette flow(RPCF) at Re_w= 1300 and Ro = 0.02 was performed with different mesh resolutions and different sizes of computation domain. Our results showed that a grid resolution in wall units with ?x^+=8.51, ?z^+= 4.26, ?y^+|_(min)= 0.0873 and ?y^+|_(max)= 3.89 is fine enough to simulate the problem at the present parameters. The streamwise length Lxand spanwise length Lzof the computational box have different impacts on the flow statistics, where the statistics were converged if Lxis longer than 8πh, while no converged results were obtained for different Lz. More importantly,our results with very long simulation time showed that a state transition would happen if L_x≥ 8πh, from a state with four pairs of roll cells to a state with three pairs of roll cells with L_z= 6πh. Each state could survive for more than 1500 h/U_w, and the flow statistics were different.
基金supported by National Natural Science Foundation of China(Grant Nos.11702127,11521091,91752202)Science Challenge Project(No.TZ2016001).
文摘This paper presents an extension work of the hybrid scheme proposed by Wang et al.[J.Comput.Phys.229(2010)169-180]for numerical simulation of sub-sonic isotropic turbulence to supersonic turbulence regime.The scheme still utilizes an 8th-order compact scheme with built-in hyperviscosity for smooth regions and a 7th-order WENO scheme for highly compression regions,but now both in their con-servation formulations and for the latter with the Roe type characteristic-wise recon-struction.To enhance the robustness of the WENO scheme without compromising its high-resolution and accuracy,the recursive-order-reduction procedure is adopted,where a new type of reconstruction-failure-detection criterion is constructed from the idea of positivity-preserving.In addition,a new form of cooling function is proposed,which is proved also to be positivity-preserving.With a combination of these techniques,the new scheme not only inherits the good properties of the original one but also extends largely the computable range of turbulent Mach number,which has been further confirmed by numerical results.
基金support from National Natural Science Foundation of China(Grants No.91130001 and No.11221061)supported by the National Basic Research Program of China(Grant No.2009CB724101)support from National Science Foundation for Postdoctoral Scientists of China(Grant No.2012M520109).
文摘Compressible flow past a circular cylinder at an inflow Reynolds number of 2×105 is numerically investigated by using a constrained large-eddy simulation(CLES)technique.Numerical simulation with adiabatic wall boundary condition and at a free-stream Mach number of 0.75 is conducted to validate and verify the performance of the present CLES method in predicting separated flows.Some typical and characteristic physical quantities,such as the drag coefficient,the root-mean-square lift fluctuations,the Strouhal number,the pressure and skin friction distributions around the cylinder,etc.are calculated and compared with previously reported experimental data,finer-grid large-eddy simulation(LES)data and those obtained in the present LES and detached-eddy simulation(DES)on coarse grids.It turns out that CLES is superior to DES in predicting such separated flow and that CLES can mimic the intricate shock wave dynamics quite well.Then,the effects of Mach number on the flow patterns and parameters such as the pressure,skin friction and drag coefficients,and the cylinder surface temperature are studied,with Mach number varying from 0.1 to 0.95.Nonmonotonic behaviors of the pressure and skin friction distributions are observed with increasing Mach number and the minimum mean separation angle occurs at a subcritical Mach number of between 0.3 and 0.5.Additionally,the wall temperature effects on the thermodynamic and aerodynamic quantities are explored in a series of simulations using isothermal wall boundary conditions at three different wall temperatures.It is found that the flow separates earlier from the cylinder surface with a longer recirculation length in the wake and a higher pressure coefficient at the rear stagnation point for higher wall temperature.Moreover,the influences of different thermal wall boundary conditions on the flow field are gradually magnified from the front stagnation point to the rear stagnation point.It is inferred that the CLES approach in its current version is a useful and effective tool for simulating wall-bounded compressible turbulent flows with massive separations.
基金the National Natural Science Foundation of China(Grant No.11172001)
文摘In this paper, we perform a numerical simulation of the cavitating flow around an underwater hemispherical-head slender body running at a high speed. For the first time, the slip boundary condition is introduced into this problem, and we find that the slip boundary condition has a big influence on the cavitation in the flow-separation zone. By simulating the cavitating flow under different cavitation numbers, we demonstrate that the slip boundary condition can effectively reduce the intensity of cavitation, as represented by the length of cavitation bubbles. The present paper provides a new method for utilization of new surface materials to control the cavitation on the underwater moving objects.
文摘Using the incompressible isotropic turbulent fields obtained from direct numerical simulation and large-eddy simulation,we studied the statistics of oscillation structures based on local zero-crossings and their relation with inertial-range intermittency for transverse velocity and passive scalar.Our results show that for both the velocity and passive scalar,the local oscillation structures are statistically scaleinvariant at high Reynolds number,and the inertial-range intermittency of the overall flow region is determined by the most intermittent structures characterized by one local zero-crossing.Local flow patterns conditioned on the oscillation structures are characterized by the joint probability density function of the invariants of the filtered velocity gradient tensor at inertial range.We demonstrate that the most intermittent regions for longitudinal velocity tend to lay at the saddle area,while those for the transverse velocity tend to locate at the vortex-dominated area.The connection between the ramp-cliff structures in passive scalar field and the corresponding saddle regions in the velocity field is also verified by the approach of oscillation structure classification.
