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.展开更多
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.展开更多
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.展开更多
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.展开更多
Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panM...Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panMARE code-to simulate the unsteady flow behavior of a full OWT with various combinations of aerodynamic and hydrodynamic loads in the time domain. This code is implemented to simulate potential flows for different applications and is based on a three-dimensional first-order panel method. Three different OWT configurations consisting of a generic 5 MW NREL rotor with three different types of foundations (Monopile, Tripod, and Jacket) are investigated. These three configurations are analyzed using the RANSE solver which is carried out using ANSYS CFX for validating the corresponding results. The simulations are performed under the same environmental atmospheric wind shear and rotor angular velocity, and the wave properties are wave height of 4 m and wave period of 7.16 s. In the present work, wave environmental effects were investigated firstly for the two solvers, and good agreement is achieved. Moreover, pressure distribution in each OWT case is presented, including detailed information about local flow fields. The time history of the forces at inflow direction and its moments around the mudline at each OWT part are presented in a dimensionless form with respect to the mean value of the last three loads and the moment amplitudes obtained from the BEM code, where the contribution of rotor force is lower in the tripod case and higher in the jacket case and the calculated hydrodynamic load that effect on jacket foundation type is lower than other two cases.展开更多
As far as the accuracy of calculating unsteady temperature field is concerned, it is very important to find the accurate physical parameters such as specific heat, thermal conductivity, latent heat of phase transforma...As far as the accuracy of calculating unsteady temperature field is concerned, it is very important to find the accurate physical parameters such as specific heat, thermal conductivity, latent heat of phase transformation and surface heat flux. The model for calculating H and Q is established in this paper. The measurement methods and data processing for physical parameters such as volume specific heat C, thermal conductivity k, volume latent heat of phase transformation c1 and surface heat flux are introduced The physical parameters of 1Cr18Ni9Ti and 45 steels and the surface heat flux for 1 Cr18Ni9Ti probe cooled in water,10% NaCl water and oil with different temperatures are measured, respectively. These data show that the probability of absolute error less than 2* C between the calculated and measured values in temperature field calculation reaches above 80% if using the above physical parameters, which provides a reliable technology basis for precise calculation of temperature field.展开更多
An analysis is presented for an unsteady boundary layer stagnation-point flow of a Newtonian fluid and the heat transfer towards a stretching sheet taking non-conventional partial slip conditions at the sheet.The self...An analysis is presented for an unsteady boundary layer stagnation-point flow of a Newtonian fluid and the heat transfer towards a stretching sheet taking non-conventional partial slip conditions at the sheet.The self-similar equations are obtained using similarity transformations and solved numerically by the shooting method.Effects of the parameters involved in the equations,especially velocity slip and thermal slip parameters on the velocity and temperature profiles,are analyzed extensively.It is revealed that due to the velocity and thermal slip parameters,the rate of heat transfer from the sheet and the wall skin friction change significantly.展开更多
This paper investigates the unsteady stagnation-point flow and heat transfer over a moving plate with mass transfer,which is also an exact solution to the unsteady Navier-Stokes(NS)equations.The boundary layer energy ...This paper investigates the unsteady stagnation-point flow and heat transfer over a moving plate with mass transfer,which is also an exact solution to the unsteady Navier-Stokes(NS)equations.The boundary layer energy equation is solved with the closed form solutions for prescribed wall temperature and prescribed wall heat flux conditions.The wall temperature and heat flux have power dependence on both time and spatial distance.The solution domain,the velocity distribution,the flow field,and the temperature distribution in the fluids are studied for different controlling parameters.These parameters include the Prandtl number,the mass transfer parameter at the wall,the wall moving parameter,the time power index,and the spatial power index.It is found that two solution branches exist for certain combinations of the controlling parameters for the flow and heat transfer problems.The heat transfer solutions are given by the confluent hypergeometric function of the first kind,which can be simplified into the incomplete gamma functions for special conditions.The wall heat flux and temperature profiles show very complicated variation behaviors.The wall heat flux can have multiple poles under certain given controlling parameters,and the temperature can have significant oscillations with overshoot and negative values in the boundary layers.The relationship between the number of poles in the wall heat flux and the number of zero-crossing points is identified.The difference in the results of the prescribed wall temperature case and the prescribed wall heat flux case is analyzed.Results given in this paper provide a rare closed form analytical solution to the entire unsteady NS equations,which can be used as a benchmark problem for numerical code validation.展开更多
The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller(HALLER, G. Exact theory ...The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller(HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257–311(2004)). By analyzing the distribution of the finite-time Lyapunov exponent(FTLE) along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE.This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.展开更多
The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic diff...The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time varia-ble. The image function of water head H can be solved by BEM. We derived the boundary integral equation ofthe transformed variable H and the discretization form of it, so that there is no need to discretize the bounda-ries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.展开更多
The effects of leading-edge blowing-suction on the vortex flow pastan airfoil at high incidence are investigated numerically by solving the Navier-Stokesequations. The results indicate that the frequency of the flowfi...The effects of leading-edge blowing-suction on the vortex flow pastan airfoil at high incidence are investigated numerically by solving the Navier-Stokesequations. The results indicate that the frequency of the flowfield excited by theperiodic blowing-suction locks into the forcing frequency, which is half of the dominantfrequency for the flow past a fixed airfoil without injection. In that case, a well-developed primary leading-edge vortex occupies the upper surface of the airfoil andthe largest lift augmentation is obtained.展开更多
In this paper, the techniques to manage and control the flow over airfoils by using the external unsteady excitations are investigated. The mechanisms of these excitation effects are also explored. The principal goal ...In this paper, the techniques to manage and control the flow over airfoils by using the external unsteady excitations are investigated. The mechanisms of these excitation effects are also explored. The principal goal of this study is to gain a better understanding and to find the possible ways for enhancing the aerodynamic efficients. The experimental investigations are carried out in two low-speed wind tunnels. The test models are two dimensional airfoils with different section geometries. Four means of excitations have been used in these experiments. (1) The pitch oscillation of the airfoil high-angle-of-attack situation. (2) The moving surface effects of the airfoil with a leading edge rotating cylinder. (3) Oscillating leading edge flaperon. (4) Small oscillating spoiler located near the leading edge of the airfoil. The lift, drag and pitch moment coefficients are measured in these experiments. But, we will put the emphasis only on the "dynamic amplifying effects" on aerodynamic lift in this paper.展开更多
A new unsteady three-dimensional convective-diffusive mathematical model for the transportation of macromolecules and water across the arterial wall was proposed . After the formation of leaky junctions due to the mit...A new unsteady three-dimensional convective-diffusive mathematical model for the transportation of macromolecules and water across the arterial wall was proposed . After the formation of leaky junctions due to the mitosis of endothelial cell of the arterial wall, the macromolecular transport happens surrounding the leaky cells. The arterial wall was divided into four layers: the endothelial layer, the subendothelial intima, the internal elastic lamina and the media for the convenience of research. The time-dependent concentration growth,the effect of the shape of endothelial cell and the effect of physiological parameters were analyzed. The analytical solution of velocity field and pressure field of water flow across the arterial wall were obtained; and concentration distribution of three macromolecules ; LDL,HRP and Albumin, were calculated with numerical simulation method. The new theory predicts, the maximum and distribution areas of time dependent concentration with round shape endothelial cell are both larger than that with ellipse-shape endothelial cell. The model also predicts the concentration growth is much alike that of a two-dimensional model and it shows that the concentration reaches its peak at the leaky junction where atherosclerotic formation frequently occurs and falls down rapidly in a limited area beginning from its earlier time growth to the state when macromolecular transfer approaches steadily. These predictions of the new model are in agreement with the experimental observation for the growth and concentration distribution of LDL and Albumin.展开更多
Unsteady aerodynamic characteristics of a seagull wing in level flight are investigated using a boundary element method.Anew no-penetration boundary condition is imposed on the surface of the wing by considering its d...Unsteady aerodynamic characteristics of a seagull wing in level flight are investigated using a boundary element method.Anew no-penetration boundary condition is imposed on the surface of the wing by considering its deformation.The geometry andkinematics of the seagull wing are reproduced using the functions and data in the previously published literature.The proposedmethod is validated by comparing the computed results with the published data in the literature.The unsteady aerodynamicscharacteristics of the seagull wing are investigated by changing flapping frequency and advance ratio.It is found that the peakvalues of aerodynamic coefficients increase with the flapping frequency.The thrust and drag generations are complicatedfunctions of frequency and wing stroke motions.The lift is inversely proportional to the advance ratio.The effects of severalflapping modes on the lift and induced drag(or thrust)generation are also investigated.Among three single modes(flapping,folding and lead & lag),flapping generates the largest lift and can produce thrust alone.For three combined modes,both flapping/foldingand flapping/lead & lag can produce lift and thrust larger than the flapping-alone mode can.Folding is shown toincrease thrust when combined with flapping,whereas lead & lag has an effect of increasing the lift when also combined withflapping.When three modes are combined together,the bird can obtain the largest lift among the investigated modes.Eventhough the proposed method is limited to the inviscid flow assumption,it is believed that this method can be used to the designof flapping micro aerial vehicle.展开更多
Weak axisymmetrical fountains resulting from the injection of a densefluid upwards into a large container of homogeneous fluid of lower density has beenstudied numerically in this paper using a time-accurate finite vo...Weak axisymmetrical fountains resulting from the injection of a densefluid upwards into a large container of homogeneous fluid of lower density has beenstudied numerically in this paper using a time-accurate finite volume scheme.Thebehaviour of fountains for both the uniform and parabolic profiles of the dischargevelocity at the fountain source has been investigated.The evolution of transient foun-tain flow has been analysed and two distinct stages of evolution have been identified.The time trace of the position of the fountain front has been presented and the initial,temporary and final fountain height and fountain width have been determined.展开更多
In this paper, the improved Background Oriented Schlieren technique called CBOS (Colored Background Oriented Schlieren) is described and used to reconstruct the density fields of three-dimensional flows. The Backgroun...In this paper, the improved Background Oriented Schlieren technique called CBOS (Colored Background Oriented Schlieren) is described and used to reconstruct the density fields of three-dimensional flows. The Background Oriented Schlieren technique (BOS) allows the measurement of the light deflection caused by density gradients in a compressible flow. For this purpose the distortion of the image of a background pattern observed through the flow is used. In order to increase the performance of the conventional Background Oriented Schlieren technique, the monochromatic background is replaced by a colored dot pattern. The different colors are treated separately using suitable correlation algorithms. Therefore, the precision and the spatial resolution can be highly increased. Furthermore a special arrangement of the different colored dot patterns in the background allows astigmatism in the region with high density gradients to be overcome. For the first time an algebraic reconstruction technique (ART) is then used to reconstruct the density field of unsteady flows around a spike-tipped model from CBOS measurements. The obtained images reveal the interaction between the free-stream flow and the high-pressure region in front of the model, which leads to large-scale instabilities in the flow.展开更多
基金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.
基金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.
基金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.
基金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.
文摘Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panMARE code-to simulate the unsteady flow behavior of a full OWT with various combinations of aerodynamic and hydrodynamic loads in the time domain. This code is implemented to simulate potential flows for different applications and is based on a three-dimensional first-order panel method. Three different OWT configurations consisting of a generic 5 MW NREL rotor with three different types of foundations (Monopile, Tripod, and Jacket) are investigated. These three configurations are analyzed using the RANSE solver which is carried out using ANSYS CFX for validating the corresponding results. The simulations are performed under the same environmental atmospheric wind shear and rotor angular velocity, and the wave properties are wave height of 4 m and wave period of 7.16 s. In the present work, wave environmental effects were investigated firstly for the two solvers, and good agreement is achieved. Moreover, pressure distribution in each OWT case is presented, including detailed information about local flow fields. The time history of the forces at inflow direction and its moments around the mudline at each OWT part are presented in a dimensionless form with respect to the mean value of the last three loads and the moment amplitudes obtained from the BEM code, where the contribution of rotor force is lower in the tripod case and higher in the jacket case and the calculated hydrodynamic load that effect on jacket foundation type is lower than other two cases.
文摘As far as the accuracy of calculating unsteady temperature field is concerned, it is very important to find the accurate physical parameters such as specific heat, thermal conductivity, latent heat of phase transformation and surface heat flux. The model for calculating H and Q is established in this paper. The measurement methods and data processing for physical parameters such as volume specific heat C, thermal conductivity k, volume latent heat of phase transformation c1 and surface heat flux are introduced The physical parameters of 1Cr18Ni9Ti and 45 steels and the surface heat flux for 1 Cr18Ni9Ti probe cooled in water,10% NaCl water and oil with different temperatures are measured, respectively. These data show that the probability of absolute error less than 2* C between the calculated and measured values in temperature field calculation reaches above 80% if using the above physical parameters, which provides a reliable technology basis for precise calculation of temperature field.
文摘An analysis is presented for an unsteady boundary layer stagnation-point flow of a Newtonian fluid and the heat transfer towards a stretching sheet taking non-conventional partial slip conditions at the sheet.The self-similar equations are obtained using similarity transformations and solved numerically by the shooting method.Effects of the parameters involved in the equations,especially velocity slip and thermal slip parameters on the velocity and temperature profiles,are analyzed extensively.It is revealed that due to the velocity and thermal slip parameters,the rate of heat transfer from the sheet and the wall skin friction change significantly.
文摘This paper investigates the unsteady stagnation-point flow and heat transfer over a moving plate with mass transfer,which is also an exact solution to the unsteady Navier-Stokes(NS)equations.The boundary layer energy equation is solved with the closed form solutions for prescribed wall temperature and prescribed wall heat flux conditions.The wall temperature and heat flux have power dependence on both time and spatial distance.The solution domain,the velocity distribution,the flow field,and the temperature distribution in the fluids are studied for different controlling parameters.These parameters include the Prandtl number,the mass transfer parameter at the wall,the wall moving parameter,the time power index,and the spatial power index.It is found that two solution branches exist for certain combinations of the controlling parameters for the flow and heat transfer problems.The heat transfer solutions are given by the confluent hypergeometric function of the first kind,which can be simplified into the incomplete gamma functions for special conditions.The wall heat flux and temperature profiles show very complicated variation behaviors.The wall heat flux can have multiple poles under certain given controlling parameters,and the temperature can have significant oscillations with overshoot and negative values in the boundary layers.The relationship between the number of poles in the wall heat flux and the number of zero-crossing points is identified.The difference in the results of the prescribed wall temperature case and the prescribed wall heat flux case is analyzed.Results given in this paper provide a rare closed form analytical solution to the entire unsteady NS equations,which can be used as a benchmark problem for numerical code validation.
基金supported by the National Natural Science Foundation of China(Nos.11372340 and 11732016)
文摘The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller(HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257–311(2004)). By analyzing the distribution of the finite-time Lyapunov exponent(FTLE) along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE.This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.
基金supported by the National Natural Science Foundation of China
文摘The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time varia-ble. The image function of water head H can be solved by BEM. We derived the boundary integral equation ofthe transformed variable H and the discretization form of it, so that there is no need to discretize the bounda-ries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.
基金The project supported by the National Defence Research Fund of China
文摘The effects of leading-edge blowing-suction on the vortex flow pastan airfoil at high incidence are investigated numerically by solving the Navier-Stokesequations. The results indicate that the frequency of the flowfield excited by theperiodic blowing-suction locks into the forcing frequency, which is half of the dominantfrequency for the flow past a fixed airfoil without injection. In that case, a well-developed primary leading-edge vortex occupies the upper surface of the airfoil andthe largest lift augmentation is obtained.
文摘In this paper, the techniques to manage and control the flow over airfoils by using the external unsteady excitations are investigated. The mechanisms of these excitation effects are also explored. The principal goal of this study is to gain a better understanding and to find the possible ways for enhancing the aerodynamic efficients. The experimental investigations are carried out in two low-speed wind tunnels. The test models are two dimensional airfoils with different section geometries. Four means of excitations have been used in these experiments. (1) The pitch oscillation of the airfoil high-angle-of-attack situation. (2) The moving surface effects of the airfoil with a leading edge rotating cylinder. (3) Oscillating leading edge flaperon. (4) Small oscillating spoiler located near the leading edge of the airfoil. The lift, drag and pitch moment coefficients are measured in these experiments. But, we will put the emphasis only on the "dynamic amplifying effects" on aerodynamic lift in this paper.
文摘A new unsteady three-dimensional convective-diffusive mathematical model for the transportation of macromolecules and water across the arterial wall was proposed . After the formation of leaky junctions due to the mitosis of endothelial cell of the arterial wall, the macromolecular transport happens surrounding the leaky cells. The arterial wall was divided into four layers: the endothelial layer, the subendothelial intima, the internal elastic lamina and the media for the convenience of research. The time-dependent concentration growth,the effect of the shape of endothelial cell and the effect of physiological parameters were analyzed. The analytical solution of velocity field and pressure field of water flow across the arterial wall were obtained; and concentration distribution of three macromolecules ; LDL,HRP and Albumin, were calculated with numerical simulation method. The new theory predicts, the maximum and distribution areas of time dependent concentration with round shape endothelial cell are both larger than that with ellipse-shape endothelial cell. The model also predicts the concentration growth is much alike that of a two-dimensional model and it shows that the concentration reaches its peak at the leaky junction where atherosclerotic formation frequently occurs and falls down rapidly in a limited area beginning from its earlier time growth to the state when macromolecular transfer approaches steadily. These predictions of the new model are in agreement with the experimental observation for the growth and concentration distribution of LDL and Albumin.
基金supported by a grant from the Academic Research Program of Chungju National University,2006supported by the Korea Research Foundation Grant funded by the korean Govemment through the Ministry of Education and Human Resources Development,Basic Research Promotion Fund(KRF-2007-331-D00081)
文摘Unsteady aerodynamic characteristics of a seagull wing in level flight are investigated using a boundary element method.Anew no-penetration boundary condition is imposed on the surface of the wing by considering its deformation.The geometry andkinematics of the seagull wing are reproduced using the functions and data in the previously published literature.The proposedmethod is validated by comparing the computed results with the published data in the literature.The unsteady aerodynamicscharacteristics of the seagull wing are investigated by changing flapping frequency and advance ratio.It is found that the peakvalues of aerodynamic coefficients increase with the flapping frequency.The thrust and drag generations are complicatedfunctions of frequency and wing stroke motions.The lift is inversely proportional to the advance ratio.The effects of severalflapping modes on the lift and induced drag(or thrust)generation are also investigated.Among three single modes(flapping,folding and lead & lag),flapping generates the largest lift and can produce thrust alone.For three combined modes,both flapping/foldingand flapping/lead & lag can produce lift and thrust larger than the flapping-alone mode can.Folding is shown toincrease thrust when combined with flapping,whereas lead & lag has an effect of increasing the lift when also combined withflapping.When three modes are combined together,the bird can obtain the largest lift among the investigated modes.Eventhough the proposed method is limited to the inviscid flow assumption,it is believed that this method can be used to the designof flapping micro aerial vehicle.
文摘Weak axisymmetrical fountains resulting from the injection of a densefluid upwards into a large container of homogeneous fluid of lower density has beenstudied numerically in this paper using a time-accurate finite volume scheme.Thebehaviour of fountains for both the uniform and parabolic profiles of the dischargevelocity at the fountain source has been investigated.The evolution of transient foun-tain flow has been analysed and two distinct stages of evolution have been identified.The time trace of the position of the fountain front has been presented and the initial,temporary and final fountain height and fountain width have been determined.
文摘In this paper, the improved Background Oriented Schlieren technique called CBOS (Colored Background Oriented Schlieren) is described and used to reconstruct the density fields of three-dimensional flows. The Background Oriented Schlieren technique (BOS) allows the measurement of the light deflection caused by density gradients in a compressible flow. For this purpose the distortion of the image of a background pattern observed through the flow is used. In order to increase the performance of the conventional Background Oriented Schlieren technique, the monochromatic background is replaced by a colored dot pattern. The different colors are treated separately using suitable correlation algorithms. Therefore, the precision and the spatial resolution can be highly increased. Furthermore a special arrangement of the different colored dot patterns in the background allows astigmatism in the region with high density gradients to be overcome. For the first time an algebraic reconstruction technique (ART) is then used to reconstruct the density field of unsteady flows around a spike-tipped model from CBOS measurements. The obtained images reveal the interaction between the free-stream flow and the high-pressure region in front of the model, which leads to large-scale instabilities in the flow.
