Modeling of unsteady aerodynamic loads at high angles of attack using a small amount of experimental or simulation data to construct predictive models for unknown states can greatly improve the efficiency of aircraft ...Modeling of unsteady aerodynamic loads at high angles of attack using a small amount of experimental or simulation data to construct predictive models for unknown states can greatly improve the efficiency of aircraft unsteady aerodynamic design and flight dynamics analysis.In this paper,aiming at the problems of poor generalization of traditional aerodynamic models and intelligent models,an intelligent aerodynamic modeling method based on gated neural units is proposed.The time memory characteristics of the gated neural unit is fully utilized,thus the nonlinear flow field characterization ability of the learning and training process is enhanced,and the generalization ability of the whole prediction model is improved.The prediction and verification of the model are carried out under the maneuvering flight condition of NACA0015 airfoil.The results show that the model has good adaptability.In the interpolation prediction,the maximum prediction error of the lift and drag coefficients and the moment coefficient does not exceed 10%,which can basically represent the variation characteristics of the entire flow field.In the construction of extrapolation models,the training model based on the strong nonlinear data has good accuracy for weak nonlinear prediction.Furthermore,the error is larger,even exceeding 20%,which indicates that the extrapolation and generalization capabilities need to be further optimized by integrating physical models.Compared with the conventional state space equation model,the proposed method can improve the extrapolation accuracy and efficiency by 78%and 60%,respectively,which demonstrates the applied potential of this method in aerodynamic modeling.展开更多
Pressure fluctuation due to rotor-stator interaction in turbomachinery is unavoidable,inducing strong vibration in the equipment and shortening its lifecycle.The investigation of optimization methods for an industrial...Pressure fluctuation due to rotor-stator interaction in turbomachinery is unavoidable,inducing strong vibration in the equipment and shortening its lifecycle.The investigation of optimization methods for an industrial centrifugal pump was carried out to reduce the intensity of pressure fluctuation to extend the lifecycle of these devices.Considering the time-consuming transient simulation of unsteady pressure,a novel optimization strategy was proposed by discretizing design variables and genetic algorithm.Four highly related design parameters were chosen,and 40 transient sample cases were generated and simulated using an automatic program.70%of them were used for training the surrogate model,and the others were for verifying the accuracy of the surrogate model.Furthermore,a modified discrete genetic algorithm(MDGA)was proposed to reduce the optimization cost owing to transient numerical simulation.For the benchmark test,the proposed MDGA showed a great advantage over the original genetic algorithm regarding searching speed and effectively dealt with the discrete variables by dramatically increasing the convergence rate.After optimization,the performance and stability of the inline pump were improved.The efficiency increased by more than 2.2%,and the pressure fluctuation intensity decreased by more than 20%under design condition.This research proposed an optimization method for reducing discrete transient characteristics in centrifugal pumps.展开更多
We propose the meshfree-based physics-informed neural networks for solving the unsteady Oseen equations.Firstly,based on the ideas of meshfree and small sample learning,we only randomly select a small number of spatio...We propose the meshfree-based physics-informed neural networks for solving the unsteady Oseen equations.Firstly,based on the ideas of meshfree and small sample learning,we only randomly select a small number of spatiotemporal points to train the neural network instead of forming a mesh.Specifically,we optimize the neural network by minimizing the loss function to satisfy the differential operators,initial condition and boundary condition.Then,we prove the convergence of the loss function and the convergence of the neural network.In addition,the feasibility and effectiveness of the method are verified by the results of numerical experiments,and the theoretical derivation is verified by the relative error between the neural network solution and the analytical solution.展开更多
The computational cost of unsteady adjoint equations remains high in adjoint-based unsteady aerodynamic op-timization.In this letter,the solution of unsteady adjoint equations is accelerated by dynamic mode decomposi-...The computational cost of unsteady adjoint equations remains high in adjoint-based unsteady aerodynamic op-timization.In this letter,the solution of unsteady adjoint equations is accelerated by dynamic mode decomposi-tion(DMD).The pseudo-time marching of every real-time step is approximated as an infinite-dimensional linear dynamical system.Thereafter,DMD is utilized to analyze the adjoint vectors sampled from these pseudo-time marching.First-order zero frequency mode is selected to accelerate the pseudo-time marching of unsteady adjoint equations in every real-time step.Through flow past a stationary circular cylinder and an unsteady aerodynamic shape optimization example,the efficiency of solving unsteady adjoint equations is significantly improved.Re-sults show that one hundred adjoint vectors contains enough information about the pseudo-time dynamics,and the adjoint dominant mode can be precisely predicted only by five snapshots produced from the adjoint vectors,which indicates DMD analysis for pseudo-time marching of unsteady adjoint equations is efficient.展开更多
This article presents an investigation into the flow and heat transfer characteristics of an impermeable stretching sheet subjected to Magnetohydrodynamic Casson fluid. The study considers the influence of slip veloci...This article presents an investigation into the flow and heat transfer characteristics of an impermeable stretching sheet subjected to Magnetohydrodynamic Casson fluid. The study considers the influence of slip velocity, thermal radiation conditions, and heat flux. The investigation is conducted employing a robust numerical method that accounts for the impact of thermal radiation. This category of fluid is apt for characterizing the movement of blood within an industrial artery, where the flow can be regulated by a material designed to manage it. The resolution of the ensuing system of ordinary differential equations (ODEs), representing the described problem, is accomplished through the application of the finite difference method. The examination of flow and heat transfer characteristics, including aspects such as unsteadiness, radiation parameter, slip velocity, Casson parameter, and Prandtl number, is explored and visually presented through tables and graphs to illustrate their impact. On the stretching sheet, calculations, and descriptions of the local skin-friction coefficient and the local Nusselt number are conducted. In conclusion, the findings indicate that the proposed method serves as a straightforward and efficient tool for exploring the solutions of fluid models of this kind.展开更多
The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic inter...The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic interaction between the blade rows. The results show that, compared to the compressor with unurved rotors, the compressor under scrutiny acquires remarkable increases in efficiency with significantly reduced amplitudes of the time-dependent fluctuation. The amplitude of the pressure fluctuation around the stator leading edge decreases at both endwalls, but increases at the mid-span in the curved rotors. The pressure fluctuation near the stator leading edge, therefore, becomes more uniform in the radial direction of this compressor. Except for the leading edge area, the pressure fluctuatinn amplitude declines remarkably in the tip region of stator surface downstream of the curved rotor, but hardly changes in the middle and at the hub.展开更多
The characteristics and mechanism of unsteady aerodynamic heating of a transient hypersonic boundary layer caused by a sudden change in surface temperature are studied. The complete time history of wall heat flux is p...The characteristics and mechanism of unsteady aerodynamic heating of a transient hypersonic boundary layer caused by a sudden change in surface temperature are studied. The complete time history of wall heat flux is presented with both analytical and numerical approaches. With the analytical method, the unsteady compressible boundary layer equation is solved. In the neighborhood of the initial and final steady states, the transient responses can be expressed with a steady-state solution plus a perturbation series. By combining these two solutions, a complete solution in the entire time domain is achieved. In the region in which the analytical approach is applicable, numerical results are in good agreement with the analytical results, showing reliability of the methods. The result shows two distinct features of the unsteady response. In a short period just after a sudden increase in the wall temperature, the direction of the wall heat flux is reverted, and a new inflexion near the wall occurs in the profile of the thermal boundary layer. This is a typical unsteady characteristic. However, these unsteady responses only exist in a very short period in hypersonic flows, meaning that, in a long-term aerodynamic heating process considering only unsteady surface temperature, the unsteady characteristics of the flow can be ignored, and the traditional quasi-steady aerodynamic heating prediction methods are still valid.展开更多
Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous f...Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous flows. In the area far away from the body, the traditional cloud of isotropic points is used, while in the adjacent area, the cloud of anisotropic points is distributed. In this way, the point spacing normal to the wall can be small enough for simulating the boundary layer, and meanwhile, the total number of points in the computational do- main can be controlled due to large spacing in other tangential direction through the anisotropic way. A fast mov- ing technique of clouds of points at each time-step is presented based on the attenuation law of disturbed motion for unsteady flows involving moving boundaries. In the mentioned cloud of points, a uniform weighted least- square curve fit method is utilized to discretize the spatial derivatives of the Navier-Stokes equations. The pro- posed gridless method, coupled with a dual time-stepping method and the Spalart-Allmaras turbulence model, is implemented for the Navier-Stokes equations. The computational results of unsteady viscous flows around a NLR7301 airfoil with an oscillating flap and a pitching NACA0012 airfoil are presented in a good agreement with the available experimental data.展开更多
The pressure fluctuation in the flow passage of both impeller and casing is addressed on design condition. The initial conditions for the unsteady turbulent simulation are resulted from the steady calculations, and th...The pressure fluctuation in the flow passage of both impeller and casing is addressed on design condition. The initial conditions for the unsteady turbulent simulation are resulted from the steady calculations, and the three dimensional unsteady turbulent simulation concerning the rotor-stator interaction is executed by a Navier-Stoke solver embedded with k -ε turbulence model and with appropriate moving interface boundary conditions. Detecting points are distributed in the flow passage in different radial and circumferential positions to capture the static pressure fluctuation character for one cycle of the impeller. The time-domain spectrums show that the static pressure curves are periodic and have five peaks and five valleys. With the radius increasing, the pressure fluctuation peak-to-peak values in the impeller are increasing, and reach the maximum value on the interface. In the casing flow passage, those values are about 7% of local static pressure except some ones near the tongue. The values become decreasingly in the diffuser pipe. The frequency spectrums transformed by fast Fourier transform (FFT) show that the dominant frequency is approximate with the blade passing frequency, and the pressure fluctuations in impeller passage have high frequency content while those in casing ones have no such information.展开更多
The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized secon...The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus.The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates.In addition,the solutions of the shear stresses at the plates are also determined.展开更多
Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of...Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.展开更多
Unsteady supersonic base flows around three afterbodies, cylindrical (Cy), boattailed (BT) and three-step (MS), are investigated in this paper. Reynolds-averaged Navier-Stokes (RANS) and two RANS/LES (large-e...Unsteady supersonic base flows around three afterbodies, cylindrical (Cy), boattailed (BT) and three-step (MS), are investigated in this paper. Reynolds-averaged Navier-Stokes (RANS) and two RANS/LES (large-eddy simulation) hybrid methods, detached eddy simulation (DES) and delayed-DES (DDES), are used to predict the base flow characteristics around the baseline Cy afterbody. All the RANS and hybrid methods are based on the two-equation SST (shear-stress transport) model with compressible corrections (CC). According to the comparison of measurements, both DES and DDES can produce more satisfactory results than RANS. RANS can only present the "stable" flow pat- terns, while the hybrid methods can demonstrate unsteady flow structures. DDES and DES results are little different from one another although the latter exhibits better agreement with the experiment. DES is taken to investigate the 5° BT and three-step afterbodies. The mean flow data and the instantaneous turbulent coherent structures are compared against available measurements.展开更多
The fractional calculus is used in the constitutive relationship model of viscoelastic fluid. A generalized Maxwell model with fractional calculus is considered. Based on the flow conditions described, two flow cases ...The fractional calculus is used in the constitutive relationship model of viscoelastic fluid. A generalized Maxwell model with fractional calculus is considered. Based on the flow conditions described, two flow cases are solved and the exact solutions are obtained by using the Weber transform and the Laplace transform for fractional calculus.展开更多
Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate f...Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate factorization algorithm and internal Newton iterations. An integral boundary layer method based on the dissipation integral is used to account for viscous effects. The computational results about unsteady transonic forces on wings, bodies and control surfaces are in agreement with experimental data.展开更多
Cavitation has a significant e ect on the flow fields and structural behaviors of a centrifugal pump. In this study, the unsteady flow and structural behaviors of a centrifugal pump are investigated numerically under ...Cavitation has a significant e ect on the flow fields and structural behaviors of a centrifugal pump. In this study, the unsteady flow and structural behaviors of a centrifugal pump are investigated numerically under di erent cavitation conditions. A strong two-way coupling fluid-structure interaction simulation is applied to obtain interior views of the e ects of cavitating bubbles on the flow and structural dynamics of a pump. The renormalization-group k-ε turbulence model and the Zwart–Gerbe–Belamri cavitation model are solved for the fluid side, while a transient structural dynamic analysis is employed for the structure side. The di erent cavitation states are mapped in the head-net positive suction head(H-NPSH) curves and flow field features inside the impeller are fully revealed. Results indicate that cavitating bubbles grow and expand rapidly with decreasing NPSH. In addition, the pressure fluctuations, both in the impeller and volute, are quantitatively analyzed and associated with the cavitation states. It is shown that influence of the cavitation on the flow field is critical, specifically in the super-cavitation state. The e ect of cavitation on the unsteady radial force and blade loads is also discussed. The results indicate that the averaged radial force increased from 8.5 N to 54.4 N in the transition progress from an onset cavitation state to a super-cavitation state. Furthermore, the structural behaviors, including blade deformation, stress, and natural frequencies, corresponding to the cavitation states are discussed. A large volume of cavitating bubbles weakens the fluid forces on the blade and decreases the natural frequencies of the rotor system. This study could enhance the understanding of the e ects of cavitation on pump flow and structural behaviors.展开更多
An analysis of unsteady boundary layer flow and heat transfer over an exponentially shrinking porous sheet filled with a copper-water nanofluid is presented.Water is treated as a base fluid.In the investigation,non-un...An analysis of unsteady boundary layer flow and heat transfer over an exponentially shrinking porous sheet filled with a copper-water nanofluid is presented.Water is treated as a base fluid.In the investigation,non-uniform mass suction through the porous sheet is considered.Using Keller-box method the transformed equations are solved numerically.The results of skin friction coefficient,the local Nusselt number as well as the velocity and temperature profiles are presented for different flow parameters.The results showed that the dual non-similar solutions exist only when certain amount of mass suction is applied through the porous sheet for various unsteady parameters and nanoparticle volume fractions.The ranges of suction where dual non-similar solution exists,become larger when values of unsteady parameter as well as nanoparticle volume fraction increase.So,due to unsteadiness of flow dynamics and the presence of nanoparticles in flow field,the requirement of mass suction for existence of solution of boundary layer flow past an exponentially shrinking sheet is less.Furthermore,the velocity boundary layer thickness decreases and thermal boundary layer thickness increases with increasing of nanoparticle volume fraction in both non-similar solutions.Whereas,for stronger mass suction,the velocity boundary layer thickness becomes thinner for the first solution and the effect is opposite in the case of second solution.The temperature inside the boundary layer increases with nanoparticle volume fraction and decreases with mass suction.So,for the unsteadiness and for the presence of nanoparticles,the flow separation is delayed to some extent.展开更多
In this paper, 2-D computational analyses were conducted for unsteady high Reynolds number flows around a smooth circular cylinder in the supercritical and upper-transition flow regimes, i.e. 8.21×104〈Re〈1.54...In this paper, 2-D computational analyses were conducted for unsteady high Reynolds number flows around a smooth circular cylinder in the supercritical and upper-transition flow regimes, i.e. 8.21×104〈Re〈1.54×106. The calculations were performed by means of solving the 2-D Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations with a k-ε turbulence model. The calculated results, produced flow structure drag and lift coefficients, as well as Strouhal numbers. The findings were in good agreement with previous published data, which also supplied us with a good understanding of the flow across cylinders of different high Reynolds numbers. Meanwhile, an effective measure was presented to control the lift force on a cylinder, which points the way to decrease the vortex induced vibration of marine structure in future.展开更多
In this study, an advanced Lagrangian vortex- boundary element method is applied to simulate the unsteady impeller-diffuser interactions in a diffuser pump not only for design but also for off-design considerations. I...In this study, an advanced Lagrangian vortex- boundary element method is applied to simulate the unsteady impeller-diffuser interactions in a diffuser pump not only for design but also for off-design considerations. In velocity calculations based on the Biot-Savart law we do not have to grid large portions of the flow field and the calculation points are concentrated in the regions where vorticity is present. Lagrangian representation of the evolving vorticity field is well suited to moving boundaries. An integral pressure equation shows that the pressure distribution can be estimated directly from the instantaneous velocity and vorticity field. The numerical results are compared with the experimental data and the comparisons show that the method used in this study can provide us insight into the complicated unsteady impeller-diffuser interaction phenomena in a diffuser pump.展开更多
Numerical simulations of unsteady cavitating flow around a NACA66-mod hydrofoil were performed using the partially-averaged Navier-Stokes method with different values of the resolution control parameters (f k =1.0-0....Numerical simulations of unsteady cavitating flow around a NACA66-mod hydrofoil were performed using the partially-averaged Navier-Stokes method with different values of the resolution control parameters (f k =1.0-0.2,fε =1).With decreasing fk,the predicted cavitating flow becomes unsteady as the time-averaged turbulent viscosity at the rear part of the attached cavity is gradually reduced.For fk =0.9 and 0.8,the cavity becomes unstable and its length dramatically expands and shrinks,but the calculation fails to predict the vapor cloud shedding behavior observed experimentally.With smaller fk less than 0.7,the cloud shedding behavior is simulated numerically and the predicted cavity shedding frequency increases.With fk =0.2,the whole cavitating flow evolution can be reasonably reproduced including the cavity growth/destabilization observed previously.The reentrant flow along the suction surface of the hydrofoil is the main trigger to cause the vapor cloud shedding.The wall pressure along the hydrofoil surface oscillates greatly due to the dynamic cavity shedding.Comparing the simulations and experiments,it is confirmed that for the PANS method,resolution control parameters of fk =0.2and fε =1 are recommended for numerical simulations of unsteady cavitating flows.Thus,the present study shows that the PANS method is an effective approach for predicting unsteady cavitating flow over hydrofoils.展开更多
Aerodynamic forces and power requirements in forward flight in a bumblebee (Bombus terrestris) were studied using the method of computational fluid dynamics. Actual wing kinematic data of free flight were used in th...Aerodynamic forces and power requirements in forward flight in a bumblebee (Bombus terrestris) were studied using the method of computational fluid dynamics. Actual wing kinematic data of free flight were used in the study (the speed ranges from 0 m/s to 4.5 m/s; advance ratio ranges from 0-0.66). The bumblebee employs the delayed stall mechanism and the fast pitching-up rotation mechanism to produce vertical force and thrust. The leading-edge vortex does not shed in the translatory phase of the half-strokes and is much more concentrated than that of the fruit fly in a previous study. At hovering and low-speed flight, the vertical force is produced by both the half-strokes and is contributed by wing lift; at medium and high speeds, the vertical force is mainly produced during the downstroke and is contributed by both wing lift and wing drag. At all speeds the thrust is mainly produced in the upstroke and is contributed by wing drag. The power requirement at low to medium speeds is not very different from that of hovering and is relatively large at the highest speed (advance ratio 0.66), i.e. the power curve is Jshaped. Except at the highest flight speed, storing energy elastically can save power up to 20%-30%. At the highest speed, because of the large increase of aerodynamic torque and the slight decrease of inertial torque (due to the smaller stroke amplitude and stroke frequency used), the power requirement is dominated by aerodynamic power and the effect of elastic storage of energy on power requirement is limited.展开更多
基金supported in part by the National Natural Science Foundation of China (No. 12202363)。
文摘Modeling of unsteady aerodynamic loads at high angles of attack using a small amount of experimental or simulation data to construct predictive models for unknown states can greatly improve the efficiency of aircraft unsteady aerodynamic design and flight dynamics analysis.In this paper,aiming at the problems of poor generalization of traditional aerodynamic models and intelligent models,an intelligent aerodynamic modeling method based on gated neural units is proposed.The time memory characteristics of the gated neural unit is fully utilized,thus the nonlinear flow field characterization ability of the learning and training process is enhanced,and the generalization ability of the whole prediction model is improved.The prediction and verification of the model are carried out under the maneuvering flight condition of NACA0015 airfoil.The results show that the model has good adaptability.In the interpolation prediction,the maximum prediction error of the lift and drag coefficients and the moment coefficient does not exceed 10%,which can basically represent the variation characteristics of the entire flow field.In the construction of extrapolation models,the training model based on the strong nonlinear data has good accuracy for weak nonlinear prediction.Furthermore,the error is larger,even exceeding 20%,which indicates that the extrapolation and generalization capabilities need to be further optimized by integrating physical models.Compared with the conventional state space equation model,the proposed method can improve the extrapolation accuracy and efficiency by 78%and 60%,respectively,which demonstrates the applied potential of this method in aerodynamic modeling.
基金Supported by National Key Research and Development Program of China(Grant No.2022YFC3202901)Natural Science Foundation of China(Grant No.51879121)+1 种基金Jiangsu Provincial Primary Research&Development Plan(Grant No.BE2019009-1)China Scholarship Council(Grant No.202108690020).
文摘Pressure fluctuation due to rotor-stator interaction in turbomachinery is unavoidable,inducing strong vibration in the equipment and shortening its lifecycle.The investigation of optimization methods for an industrial centrifugal pump was carried out to reduce the intensity of pressure fluctuation to extend the lifecycle of these devices.Considering the time-consuming transient simulation of unsteady pressure,a novel optimization strategy was proposed by discretizing design variables and genetic algorithm.Four highly related design parameters were chosen,and 40 transient sample cases were generated and simulated using an automatic program.70%of them were used for training the surrogate model,and the others were for verifying the accuracy of the surrogate model.Furthermore,a modified discrete genetic algorithm(MDGA)was proposed to reduce the optimization cost owing to transient numerical simulation.For the benchmark test,the proposed MDGA showed a great advantage over the original genetic algorithm regarding searching speed and effectively dealt with the discrete variables by dramatically increasing the convergence rate.After optimization,the performance and stability of the inline pump were improved.The efficiency increased by more than 2.2%,and the pressure fluctuation intensity decreased by more than 20%under design condition.This research proposed an optimization method for reducing discrete transient characteristics in centrifugal pumps.
基金Project supported in part by the National Natural Science Foundation of China(Grant No.11771259)Shaanxi Provincial Joint Laboratory of Artificial Intelligence(GrantNo.2022JCSYS05)+1 种基金Innovative Team Project of Shaanxi Provincial Department of Education(Grant No.21JP013)Shaanxi Provincial Social Science Fund Annual Project(Grant No.2022D332)。
文摘We propose the meshfree-based physics-informed neural networks for solving the unsteady Oseen equations.Firstly,based on the ideas of meshfree and small sample learning,we only randomly select a small number of spatiotemporal points to train the neural network instead of forming a mesh.Specifically,we optimize the neural network by minimizing the loss function to satisfy the differential operators,initial condition and boundary condition.Then,we prove the convergence of the loss function and the convergence of the neural network.In addition,the feasibility and effectiveness of the method are verified by the results of numerical experiments,and the theoretical derivation is verified by the relative error between the neural network solution and the analytical solution.
基金the Natural Science Foundation of Jiangsu Province(Grants No.BK20230202)Basic Science(Natural Science)Re-search Project of Colleges and Universities in Jiangsu Province(Grant No.22KJB130005)+3 种基金Changzhou Science and Technology Project(Grant No.CJ20220242)for financial supportJiaqing Kou would like to thank the support of the Alexander von Humboldt Foundation(Ref 3.5-CHN-1227287-HFST-P)Wenkai Yang would like to thank the support of the National Natural Science Foundation of China(Grant No.52205335)supported by Changzhou Sci&Tech Pro-gram(Grant No.CM20223013).
文摘The computational cost of unsteady adjoint equations remains high in adjoint-based unsteady aerodynamic op-timization.In this letter,the solution of unsteady adjoint equations is accelerated by dynamic mode decomposi-tion(DMD).The pseudo-time marching of every real-time step is approximated as an infinite-dimensional linear dynamical system.Thereafter,DMD is utilized to analyze the adjoint vectors sampled from these pseudo-time marching.First-order zero frequency mode is selected to accelerate the pseudo-time marching of unsteady adjoint equations in every real-time step.Through flow past a stationary circular cylinder and an unsteady aerodynamic shape optimization example,the efficiency of solving unsteady adjoint equations is significantly improved.Re-sults show that one hundred adjoint vectors contains enough information about the pseudo-time dynamics,and the adjoint dominant mode can be precisely predicted only by five snapshots produced from the adjoint vectors,which indicates DMD analysis for pseudo-time marching of unsteady adjoint equations is efficient.
文摘This article presents an investigation into the flow and heat transfer characteristics of an impermeable stretching sheet subjected to Magnetohydrodynamic Casson fluid. The study considers the influence of slip velocity, thermal radiation conditions, and heat flux. The investigation is conducted employing a robust numerical method that accounts for the impact of thermal radiation. This category of fluid is apt for characterizing the movement of blood within an industrial artery, where the flow can be regulated by a material designed to manage it. The resolution of the ensuing system of ordinary differential equations (ODEs), representing the described problem, is accomplished through the application of the finite difference method. The examination of flow and heat transfer characteristics, including aspects such as unsteadiness, radiation parameter, slip velocity, Casson parameter, and Prandtl number, is explored and visually presented through tables and graphs to illustrate their impact. On the stretching sheet, calculations, and descriptions of the local skin-friction coefficient and the local Nusselt number are conducted. In conclusion, the findings indicate that the proposed method serves as a straightforward and efficient tool for exploring the solutions of fluid models of this kind.
基金National Natural Science Foundation of China (506460210) Chinese Specialized Research Fund for the Doctoral Program of Higher Education (20060213007)Development Program for Outstanding Young Teachers in Harbin Institute of Technology (HITQNJS.2006.046)
文摘The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic interaction between the blade rows. The results show that, compared to the compressor with unurved rotors, the compressor under scrutiny acquires remarkable increases in efficiency with significantly reduced amplitudes of the time-dependent fluctuation. The amplitude of the pressure fluctuation around the stator leading edge decreases at both endwalls, but increases at the mid-span in the curved rotors. The pressure fluctuation near the stator leading edge, therefore, becomes more uniform in the radial direction of this compressor. Except for the leading edge area, the pressure fluctuatinn amplitude declines remarkably in the tip region of stator surface downstream of the curved rotor, but hardly changes in the middle and at the hub.
基金supported by the National Natural Science Foundation of China (No. 90716011)
文摘The characteristics and mechanism of unsteady aerodynamic heating of a transient hypersonic boundary layer caused by a sudden change in surface temperature are studied. The complete time history of wall heat flux is presented with both analytical and numerical approaches. With the analytical method, the unsteady compressible boundary layer equation is solved. In the neighborhood of the initial and final steady states, the transient responses can be expressed with a steady-state solution plus a perturbation series. By combining these two solutions, a complete solution in the entire time domain is achieved. In the region in which the analytical approach is applicable, numerical results are in good agreement with the analytical results, showing reliability of the methods. The result shows two distinct features of the unsteady response. In a short period just after a sudden increase in the wall temperature, the direction of the wall heat flux is reverted, and a new inflexion near the wall occurs in the profile of the thermal boundary layer. This is a typical unsteady characteristic. However, these unsteady responses only exist in a very short period in hypersonic flows, meaning that, in a long-term aerodynamic heating process considering only unsteady surface temperature, the unsteady characteristics of the flow can be ignored, and the traditional quasi-steady aerodynamic heating prediction methods are still valid.
基金Supported by the National Natural Science Foundation of China(10372043,11172134)the Fundingof Jiangsu Innovation Program for Graduate Education(CXZZ11-0192)~~
文摘Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous flows. In the area far away from the body, the traditional cloud of isotropic points is used, while in the adjacent area, the cloud of anisotropic points is distributed. In this way, the point spacing normal to the wall can be small enough for simulating the boundary layer, and meanwhile, the total number of points in the computational do- main can be controlled due to large spacing in other tangential direction through the anisotropic way. A fast mov- ing technique of clouds of points at each time-step is presented based on the attenuation law of disturbed motion for unsteady flows involving moving boundaries. In the mentioned cloud of points, a uniform weighted least- square curve fit method is utilized to discretize the spatial derivatives of the Navier-Stokes equations. The pro- posed gridless method, coupled with a dual time-stepping method and the Spalart-Allmaras turbulence model, is implemented for the Navier-Stokes equations. The computational results of unsteady viscous flows around a NLR7301 airfoil with an oscillating flap and a pitching NACA0012 airfoil are presented in a good agreement with the available experimental data.
基金supported by National Outstanding Young Scientists Funds of China (Grand No.50825902)
文摘The pressure fluctuation in the flow passage of both impeller and casing is addressed on design condition. The initial conditions for the unsteady turbulent simulation are resulted from the steady calculations, and the three dimensional unsteady turbulent simulation concerning the rotor-stator interaction is executed by a Navier-Stoke solver embedded with k -ε turbulence model and with appropriate moving interface boundary conditions. Detecting points are distributed in the flow passage in different radial and circumferential positions to capture the static pressure fluctuation character for one cycle of the impeller. The time-domain spectrums show that the static pressure curves are periodic and have five peaks and five valleys. With the radius increasing, the pressure fluctuation peak-to-peak values in the impeller are increasing, and reach the maximum value on the interface. In the casing flow passage, those values are about 7% of local static pressure except some ones near the tongue. The values become decreasingly in the diffuser pipe. The frequency spectrums transformed by fast Fourier transform (FFT) show that the dominant frequency is approximate with the blade passing frequency, and the pressure fluctuations in impeller passage have high frequency content while those in casing ones have no such information.
基金The project supported by the National Natural Science Foundation of China (10372007,10002003) and CNPC Innovation Fund
文摘The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus.The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates.In addition,the solutions of the shear stresses at the plates are also determined.
基金Project(51375498)supported by the National Natural Science Foundation of China
文摘Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.
基金supported by the National Natural Science Foundation of China (10502030 and 90505005)Innovation and Support Foundation of Chinese Astronautics
文摘Unsteady supersonic base flows around three afterbodies, cylindrical (Cy), boattailed (BT) and three-step (MS), are investigated in this paper. Reynolds-averaged Navier-Stokes (RANS) and two RANS/LES (large-eddy simulation) hybrid methods, detached eddy simulation (DES) and delayed-DES (DDES), are used to predict the base flow characteristics around the baseline Cy afterbody. All the RANS and hybrid methods are based on the two-equation SST (shear-stress transport) model with compressible corrections (CC). According to the comparison of measurements, both DES and DDES can produce more satisfactory results than RANS. RANS can only present the "stable" flow pat- terns, while the hybrid methods can demonstrate unsteady flow structures. DDES and DES results are little different from one another although the latter exhibits better agreement with the experiment. DES is taken to investigate the 5° BT and three-step afterbodies. The mean flow data and the instantaneous turbulent coherent structures are compared against available measurements.
基金The project supported by the National Natural Science Foundation of China (10272067, 10426024)the Doctoral Program Foundation of the Education Ministry of China (20030422046)the Natural Science Foundation of Shandong University at Weihai.
文摘The fractional calculus is used in the constitutive relationship model of viscoelastic fluid. A generalized Maxwell model with fractional calculus is considered. Based on the flow conditions described, two flow cases are solved and the exact solutions are obtained by using the Weber transform and the Laplace transform for fractional calculus.
基金Aeronautical Science Foundation of China (99A52007)
文摘Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate factorization algorithm and internal Newton iterations. An integral boundary layer method based on the dissipation integral is used to account for viscous effects. The computational results about unsteady transonic forces on wings, bodies and control surfaces are in agreement with experimental data.
基金Supported by National Natural Science Foundation of China(Grant Nos.51609212,51606167)China Postdoctoral Science Foundation(Grant No.2016M590546)Zhejiang Provincial Natural Science Foundation(Grant No.2016C31043)
文摘Cavitation has a significant e ect on the flow fields and structural behaviors of a centrifugal pump. In this study, the unsteady flow and structural behaviors of a centrifugal pump are investigated numerically under di erent cavitation conditions. A strong two-way coupling fluid-structure interaction simulation is applied to obtain interior views of the e ects of cavitating bubbles on the flow and structural dynamics of a pump. The renormalization-group k-ε turbulence model and the Zwart–Gerbe–Belamri cavitation model are solved for the fluid side, while a transient structural dynamic analysis is employed for the structure side. The di erent cavitation states are mapped in the head-net positive suction head(H-NPSH) curves and flow field features inside the impeller are fully revealed. Results indicate that cavitating bubbles grow and expand rapidly with decreasing NPSH. In addition, the pressure fluctuations, both in the impeller and volute, are quantitatively analyzed and associated with the cavitation states. It is shown that influence of the cavitation on the flow field is critical, specifically in the super-cavitation state. The e ect of cavitation on the unsteady radial force and blade loads is also discussed. The results indicate that the averaged radial force increased from 8.5 N to 54.4 N in the transition progress from an onset cavitation state to a super-cavitation state. Furthermore, the structural behaviors, including blade deformation, stress, and natural frequencies, corresponding to the cavitation states are discussed. A large volume of cavitating bubbles weakens the fluid forces on the blade and decreases the natural frequencies of the rotor system. This study could enhance the understanding of the e ects of cavitation on pump flow and structural behaviors.
基金the National Board for Higher Mathematics (NBHM),Department of Atomic Energy,Government of India for the financial support in pursuing this workthe financial support from MOHE and the Research Management Center-UTM through FRGS and RUG vote number 4F109 and 02H80 for this research
文摘An analysis of unsteady boundary layer flow and heat transfer over an exponentially shrinking porous sheet filled with a copper-water nanofluid is presented.Water is treated as a base fluid.In the investigation,non-uniform mass suction through the porous sheet is considered.Using Keller-box method the transformed equations are solved numerically.The results of skin friction coefficient,the local Nusselt number as well as the velocity and temperature profiles are presented for different flow parameters.The results showed that the dual non-similar solutions exist only when certain amount of mass suction is applied through the porous sheet for various unsteady parameters and nanoparticle volume fractions.The ranges of suction where dual non-similar solution exists,become larger when values of unsteady parameter as well as nanoparticle volume fraction increase.So,due to unsteadiness of flow dynamics and the presence of nanoparticles in flow field,the requirement of mass suction for existence of solution of boundary layer flow past an exponentially shrinking sheet is less.Furthermore,the velocity boundary layer thickness decreases and thermal boundary layer thickness increases with increasing of nanoparticle volume fraction in both non-similar solutions.Whereas,for stronger mass suction,the velocity boundary layer thickness becomes thinner for the first solution and the effect is opposite in the case of second solution.The temperature inside the boundary layer increases with nanoparticle volume fraction and decreases with mass suction.So,for the unsteadiness and for the presence of nanoparticles,the flow separation is delayed to some extent.
基金Foundation item: Supported by Supported by the National Natural Science Foundation of China (Grant No. 51009070).
文摘In this paper, 2-D computational analyses were conducted for unsteady high Reynolds number flows around a smooth circular cylinder in the supercritical and upper-transition flow regimes, i.e. 8.21×104〈Re〈1.54×106. The calculations were performed by means of solving the 2-D Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations with a k-ε turbulence model. The calculated results, produced flow structure drag and lift coefficients, as well as Strouhal numbers. The findings were in good agreement with previous published data, which also supplied us with a good understanding of the flow across cylinders of different high Reynolds numbers. Meanwhile, an effective measure was presented to control the lift force on a cylinder, which points the way to decrease the vortex induced vibration of marine structure in future.
文摘In this study, an advanced Lagrangian vortex- boundary element method is applied to simulate the unsteady impeller-diffuser interactions in a diffuser pump not only for design but also for off-design considerations. In velocity calculations based on the Biot-Savart law we do not have to grid large portions of the flow field and the calculation points are concentrated in the regions where vorticity is present. Lagrangian representation of the evolving vorticity field is well suited to moving boundaries. An integral pressure equation shows that the pressure distribution can be estimated directly from the instantaneous velocity and vorticity field. The numerical results are compared with the experimental data and the comparisons show that the method used in this study can provide us insight into the complicated unsteady impeller-diffuser interaction phenomena in a diffuser pump.
基金Supported by the National Natural Science Fund of China under Grant Nos 51179091 and 50976061, the Major NationalScientific Instrument and Equipment Development Project (No 2011YQ07004901), and the China Postdoctoral Science Foundation(No 2011M500314).
文摘Numerical simulations of unsteady cavitating flow around a NACA66-mod hydrofoil were performed using the partially-averaged Navier-Stokes method with different values of the resolution control parameters (f k =1.0-0.2,fε =1).With decreasing fk,the predicted cavitating flow becomes unsteady as the time-averaged turbulent viscosity at the rear part of the attached cavity is gradually reduced.For fk =0.9 and 0.8,the cavity becomes unstable and its length dramatically expands and shrinks,but the calculation fails to predict the vapor cloud shedding behavior observed experimentally.With smaller fk less than 0.7,the cloud shedding behavior is simulated numerically and the predicted cavity shedding frequency increases.With fk =0.2,the whole cavitating flow evolution can be reasonably reproduced including the cavity growth/destabilization observed previously.The reentrant flow along the suction surface of the hydrofoil is the main trigger to cause the vapor cloud shedding.The wall pressure along the hydrofoil surface oscillates greatly due to the dynamic cavity shedding.Comparing the simulations and experiments,it is confirmed that for the PANS method,resolution control parameters of fk =0.2and fε =1 are recommended for numerical simulations of unsteady cavitating flows.Thus,the present study shows that the PANS method is an effective approach for predicting unsteady cavitating flow over hydrofoils.
基金The project supported by the National Natural Science Foundation of China(10232010)the National Aeronautic Science fund of China(03A51049)
文摘Aerodynamic forces and power requirements in forward flight in a bumblebee (Bombus terrestris) were studied using the method of computational fluid dynamics. Actual wing kinematic data of free flight were used in the study (the speed ranges from 0 m/s to 4.5 m/s; advance ratio ranges from 0-0.66). The bumblebee employs the delayed stall mechanism and the fast pitching-up rotation mechanism to produce vertical force and thrust. The leading-edge vortex does not shed in the translatory phase of the half-strokes and is much more concentrated than that of the fruit fly in a previous study. At hovering and low-speed flight, the vertical force is produced by both the half-strokes and is contributed by wing lift; at medium and high speeds, the vertical force is mainly produced during the downstroke and is contributed by both wing lift and wing drag. At all speeds the thrust is mainly produced in the upstroke and is contributed by wing drag. The power requirement at low to medium speeds is not very different from that of hovering and is relatively large at the highest speed (advance ratio 0.66), i.e. the power curve is Jshaped. Except at the highest flight speed, storing energy elastically can save power up to 20%-30%. At the highest speed, because of the large increase of aerodynamic torque and the slight decrease of inertial torque (due to the smaller stroke amplitude and stroke frequency used), the power requirement is dominated by aerodynamic power and the effect of elastic storage of energy on power requirement is limited.