A novel numerical model based on the image Green function and first-order Taylor expansion boundary element method(TEBEM), which can improve the accuracy of the hydrodynamic simulation for the non-smooth body, was dev...A novel numerical model based on the image Green function and first-order Taylor expansion boundary element method(TEBEM), which can improve the accuracy of the hydrodynamic simulation for the non-smooth body, was developed to calculate the side wall effects on first-order motion responses and second-order drift loads upon offshore structures in the wave tank. This model was confirmed by comparing it to the results from experiments on hydrodynamic coefficients, namely the first-order motion response and second-order drift load upon a hemisphere, prolate spheroid, and box-shaped barge in the wave tank. Then,the hydrodynamics of the KVLCC2 model were also calculated in two wave tanks with different widths. It was concluded that this model can predict the hydrodynamics for offshore structures effectively, and the side wall has a significant impact on the firstorder quantities and second-order drift loads, which satisfied the resonant rule.展开更多
According to the Liutex-shear decomposition,vorticity can be decomposed into a rotational part,i.e.,the Liutex vector,and a residual shear part.With this decomposition,the vorticity transport equation can be used to f...According to the Liutex-shear decomposition,vorticity can be decomposed into a rotational part,i.e.,the Liutex vector,and a residual shear part.With this decomposition,the vorticity transport equation can be used to formulate a governing equation for Liutex easily for two-dimensional incompressible flows with a source term depending on the residual shear.The dynamics of Liutex-identified structures is then studied in a Taylor-Green vortex flow and a flow past a cylinder at Reynolds number of 200.It is revealed that such boundaries exist outside which the shear has trivial impact on the evolution of Liutex and inside which enhancing and weakening effects of shear on Liutex can be observed.In addition,there is a strong dissipation effect upon Liutex on these boundaries.Based on the interaction mechanism between Liutex and shear,we argue that the vortex boundaries can be identified by these highly dissipative boundaries.In contrast,traditional methods use iso-surfaces of arbitrarily selected thresholds to represent vortex boundaries.The current method of identifying vortex boundaries based on the Liutex-shear interaction has a clearer theoretical base and avoids the arbitrary selection of thresholds.Extensions to three-dimensional incompressible flows can be made in future following the same procedure but with a slightly more complex vorticity transport equation which includes the velocity gradient induced stretching or tilting term.展开更多
The accurate simulation of turbulence and the implementation of corresponding turbulence models are both critical to the understanding of the complex physics behind turbulent flows in a variety of science and engineer...The accurate simulation of turbulence and the implementation of corresponding turbulence models are both critical to the understanding of the complex physics behind turbulent flows in a variety of science and engineering applications.Despite the tremendous increase in the computing power of central processing units(CPUs),direct numerical simulation of highly turbulent flows is still not feasible due to the need for resolving the smallest length scale,and today’s CPUs cannot keep pace with demand.The recent development of graphics processing units(GPU)has led to the general improvement in the performance of various algorithms.This study investigates the applicability of GPU technology in the context of fast-Fourier transform(FFT)-based pseudo-spectral methods for DNS of turbulent flows for the Taylor–Green vortex problem.They are implemented on a single GPU and a speedup of unto 31x is obtained in comparison to a single CPU.展开更多
Numerical instability may occur when simulating high Reynolds number flows by the lattice Boltzmann method(LBM).The multiple-relaxation-time(MRT)model of the LBM can improve the accuracy and stability,but is still sub...Numerical instability may occur when simulating high Reynolds number flows by the lattice Boltzmann method(LBM).The multiple-relaxation-time(MRT)model of the LBM can improve the accuracy and stability,but is still subject to numerical instability when simulating flows with large single-grid Reynolds number(Reynolds number/grid number).The viscosity counteracting approach proposed recently is a method of enhancing the stability of the LBM.However,its effectiveness was only verified in the single-relaxation-time model of the LBM(SRT-LBM).This paper aims to propose the viscosity counteracting approach for the multiple-relaxationtime model(MRT-LBM)and analyze its numerical characteristics.The verification is conducted by simulating some benchmark cases:the two-dimensional(2D)lid-driven cavity flow,Poiseuille flow,Taylor-Green vortex flow and Couette flow,and threedimensional(3D)rectangular jet.Qualitative and Quantitative comparisons show that the viscosity counteracting approach for the MRT-LBMhas better accuracy and stability than that for the SRT-LBM.展开更多
基金the National Natural Science Foundation of China (Grant No.51709064)the Numerical Tank Project sponsored by the Ministry of Industry and Information Technology (MIIT)of P.R.China.
文摘A novel numerical model based on the image Green function and first-order Taylor expansion boundary element method(TEBEM), which can improve the accuracy of the hydrodynamic simulation for the non-smooth body, was developed to calculate the side wall effects on first-order motion responses and second-order drift loads upon offshore structures in the wave tank. This model was confirmed by comparing it to the results from experiments on hydrodynamic coefficients, namely the first-order motion response and second-order drift load upon a hemisphere, prolate spheroid, and box-shaped barge in the wave tank. Then,the hydrodynamics of the KVLCC2 model were also calculated in two wave tanks with different widths. It was concluded that this model can predict the hydrodynamics for offshore structures effectively, and the side wall has a significant impact on the firstorder quantities and second-order drift loads, which satisfied the resonant rule.
基金supported by the Jiangsu Shuangchuang Project(Grant No.JSSCTD202209)the National Science Foundation of the Jiangsu Higher Education Institutions of China(Grant No.22KJB130011)the Supercomputing Center in Yancheng(Grant No.FW(W)20221001).
文摘According to the Liutex-shear decomposition,vorticity can be decomposed into a rotational part,i.e.,the Liutex vector,and a residual shear part.With this decomposition,the vorticity transport equation can be used to formulate a governing equation for Liutex easily for two-dimensional incompressible flows with a source term depending on the residual shear.The dynamics of Liutex-identified structures is then studied in a Taylor-Green vortex flow and a flow past a cylinder at Reynolds number of 200.It is revealed that such boundaries exist outside which the shear has trivial impact on the evolution of Liutex and inside which enhancing and weakening effects of shear on Liutex can be observed.In addition,there is a strong dissipation effect upon Liutex on these boundaries.Based on the interaction mechanism between Liutex and shear,we argue that the vortex boundaries can be identified by these highly dissipative boundaries.In contrast,traditional methods use iso-surfaces of arbitrarily selected thresholds to represent vortex boundaries.The current method of identifying vortex boundaries based on the Liutex-shear interaction has a clearer theoretical base and avoids the arbitrary selection of thresholds.Extensions to three-dimensional incompressible flows can be made in future following the same procedure but with a slightly more complex vorticity transport equation which includes the velocity gradient induced stretching or tilting term.
文摘The accurate simulation of turbulence and the implementation of corresponding turbulence models are both critical to the understanding of the complex physics behind turbulent flows in a variety of science and engineering applications.Despite the tremendous increase in the computing power of central processing units(CPUs),direct numerical simulation of highly turbulent flows is still not feasible due to the need for resolving the smallest length scale,and today’s CPUs cannot keep pace with demand.The recent development of graphics processing units(GPU)has led to the general improvement in the performance of various algorithms.This study investigates the applicability of GPU technology in the context of fast-Fourier transform(FFT)-based pseudo-spectral methods for DNS of turbulent flows for the Taylor–Green vortex problem.They are implemented on a single GPU and a speedup of unto 31x is obtained in comparison to a single CPU.
基金supported by the National Natural Science Foundation of China(NSFC,Grant Numbers 10572106,10872153 and 11172219)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20130141110013).
文摘Numerical instability may occur when simulating high Reynolds number flows by the lattice Boltzmann method(LBM).The multiple-relaxation-time(MRT)model of the LBM can improve the accuracy and stability,but is still subject to numerical instability when simulating flows with large single-grid Reynolds number(Reynolds number/grid number).The viscosity counteracting approach proposed recently is a method of enhancing the stability of the LBM.However,its effectiveness was only verified in the single-relaxation-time model of the LBM(SRT-LBM).This paper aims to propose the viscosity counteracting approach for the multiple-relaxationtime model(MRT-LBM)and analyze its numerical characteristics.The verification is conducted by simulating some benchmark cases:the two-dimensional(2D)lid-driven cavity flow,Poiseuille flow,Taylor-Green vortex flow and Couette flow,and threedimensional(3D)rectangular jet.Qualitative and Quantitative comparisons show that the viscosity counteracting approach for the MRT-LBMhas better accuracy and stability than that for the SRT-LBM.