With the continuous increase of aeroengine flight ceiling(>20 km),the thin atmosphere at high altitudes and the size effect all cause the compressor component inlet Reynolds number to decrease rapidly to a critical...With the continuous increase of aeroengine flight ceiling(>20 km),the thin atmosphere at high altitudes and the size effect all cause the compressor component inlet Reynolds number to decrease rapidly to a critical value(approximately 2.0×10^(5)),and the significant transition process on the blade/endwall surface leads to the sharp degradation of compressor performance,which seriously affects the engine fuel consumption and working stability at high altitudes.In this paper,the research progress on the internal flow mechanism and flow control methods of axial compressors at low Reynolds numbers is reviewed from the aspects of quantification and prediction of performance variation,flow loss mechanism related to separation and transition,efficient transition control and flow field organization.The development trend of the low-Reynolds-number effect of axial flow compressors is noted,and the difficulties and application prospects of aerodynamic design and efficient flow control methods for compressors under low Reynolds numbers at high altitudes are discussed.展开更多
Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of un...Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.展开更多
The flow ripple caused by an axial piston pump may lead to pipe vibrations and lower hydraulic component reliability,which are of particular concern in hydraulic systems.The valve plate of the pump is considered the p...The flow ripple caused by an axial piston pump may lead to pipe vibrations and lower hydraulic component reliability,which are of particular concern in hydraulic systems.The valve plate of the pump is considered the part most related to flow ripple,and its structural design is an important topic.In this study,an analytical model for the axial piston pump flow ripple was established and verified using a numerical analysis with computational fluid dynamics(CFD)calculations.Moreover,a parametric analysis of the valve plate was performed to investigate the critical parameters and their ranges.A fast optimization method,the rotation vector optimization method(RVOM),was proposed for the valve plate design and compared with the currently used optimization methods to prove its efficiency.As a constant-pressure pump works in different states of swashplate angle,outlet pressure,and pump speed,an optimization principle for the entire working status was proposed to achieve the overall reduction performance.A test rig for an aircraft hydraulic pump was established,and validation experiments were conducted.It was determined that the optimized pump could achieve reduction at multiple working statuses,and the largest pressure pulsation reduction ratios for the typical speed and speed sweep tests reached 64.7%and 71.7%,respectively.The model and method proposed in this study are proven to be effective and accurate.展开更多
Vortex-induced vibration(VIV)of an underwater manipulator in pulsating flow presents a notable engineering problem in precise control due to the velocity variation in the flow.This study investigates the VIV response ...Vortex-induced vibration(VIV)of an underwater manipulator in pulsating flow presents a notable engineering problem in precise control due to the velocity variation in the flow.This study investigates the VIV response of an underwater manipulator subjected to pulsating flow,focusing on how different postures affect the behavior of the system.The effects of pulsating parameters and manipulator arrangement on the hydrodynamic coefficient,vibration response,motion trajectory,and vortex shedding behaviors were analyzed.Results indicated that the cross flow vibration displacement in pulsating flow increased by 32.14%compared to uniform flow,inducing a shift in the motion trajectory from a crescent shape to a sideward vase shape.In the absence of interference between the upper and lower arms,the lift coefficient of the manipulator substantially increased with rising pulsating frequency,reaching a maximum increment of 67.0%.This increase in the lift coefficient led to a 67.05%rise in the vibration frequency of the manipulator in the in-line direction.As the pulsating amplitude increased,the drag coefficient of the underwater manipulator rose by 36.79%,but the vibration frequency in the cross-flow direction decreased by 56.26%.Additionally,when the upper and lower arms remained in a state of mutual interference,the cross-flow vibration amplitudes of the upper and lower arms were approximately 1.84 and 4.82 times higher in a circular-elliptical arrangement compared to an elliptical-circular arrangement,respectively.Consequently,the flow field shifted from a P+S pattern to a disordered pattern,disrupting the regularity of the motion trajectory.展开更多
The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled wi...The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled with flow field. The microsegregation, concentration maximum value, boundary thickness of concentration near upstream dendrite and normal to flow dendrite, and downstream dendrite were studied quantitatively in the case of forced lamia flow. The simulation results show that solute field and flow field interact complexly. Compared with melt without flow, in front of upstream dendrite tip, the concentration boundary thickness is the lowest and the concentration maximum value is the smallest for melt with flow. However, in front of downstream dendrite tip, the results are just the opposite. The zone of poor Cu in upstream dendrite where is the most severely microsegregation and shrinkage cavity is wider and the concentration is lower for melt with flow than that without flow.展开更多
A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic gr...A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.展开更多
The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three ...The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.展开更多
The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, ...The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma：0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.展开更多
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.展开更多
Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks....Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks. The first one is designed on the path-turn incidencerelationship, and it is similar to the computational procedure of link flows. It applies to thetraffic assignment algorithms that can provide detailed path structures. The second utilizes thelink-turn incidence relationship and the conservation of flow on links, a law deriving from thisrelationship. It is actually an improved version of Dial's logit assignment algorithm. The proposedapproaches can avoid the shortcomings both of the estimation methods, e. g. Furness's model andFrator's model, and of the network-expanding method in precision, stability and computation scale.Finally, they are validated by numerical examples.展开更多
From the viewpoint of systems energy conservation, the influences of material flow on its energy consumption in a steel manufacturing process is an important subject. The quantitative analysis of the relationship betw...From the viewpoint of systems energy conservation, the influences of material flow on its energy consumption in a steel manufacturing process is an important subject. The quantitative analysis of the relationship between material flow and the energy intensity is useful to save energy in steel industry. Based on the concept of standard material flow diagram, all possible situations of ferric material flow in steel manufacturing process are analyzed. The expressions of the influence of material flow deviated from standard material flow diagram on energy consumption are put forward.展开更多
In this study, the effects of magnetic field and nanoparticle on the Jeffery- Hamel flow are studied using a powerful analytical method called the Adomian decomposition method （ADM）. The traditional Navier-Stokes eq...In this study, the effects of magnetic field and nanoparticle on the Jeffery- Hamel flow are studied using a powerful analytical method called the Adomian decomposition method （ADM）. The traditional Navier-Stokes equation of fluid mechanics and Maxwell＇s electromagnetism governing equations are reduced to nonlinear ordinary differential equations to model the problem. The obtained results are well agreed with that of the Runge-Kutta method. The present plots confirm that the method has high accuracy for different a, Ha, and Re numbers. The flow field inside the divergent channel is studied for various values of Hartmann ：number and angle of channel. The effect of nanoparticle volume fraction in the absence of magnetic field is investigated.展开更多
The motivation of this study is to investigate the turbulence-chemistry interactions by using probability density function （PDF） method. A consistent hybrid Reynolds Averaged Navier-Stokes （RANS）/PDF method is use...The motivation of this study is to investigate the turbulence-chemistry interactions by using probability density function （PDF） method. A consistent hybrid Reynolds Averaged Navier-Stokes （RANS）/PDF method is used to simulate the turbulent non-reacting and reacting flows. The joint fluctuating velocity-frequency-composition PDF equation coupled with the Reynolds averaged density, momentum and energy equations are solved on unstructured meshes by the Lagrangian Monte Carlo （MC） method combined with the finite volume （FV） method. The simulation of the axisymmetric bluff body stabilized non-reacting flow fields is presented in this paper. The calculated length of the recirculation zone is in good agreement with the experimental data. Moreover, the significant change of the flow pattern with the increase of the jet-to-coflow momentum flux ratio is well predicted. In addition, comparisons are made between the joint PDF model and two different Reynolds stress models.展开更多
Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on e...Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.展开更多
The main goal of this paper is to investigate natural convective heat transfer and flow characteristics of non-Newtonian nanofluid streaming between two infinite vertical flat plates in the presence of magnetic field ...The main goal of this paper is to investigate natural convective heat transfer and flow characteristics of non-Newtonian nanofluid streaming between two infinite vertical flat plates in the presence of magnetic field and thermal radiation.Initially,a similarity transformation is used to convert momentum and energy conservation equations in partial differential forms into non-linear ordinary differential equations (ODE) applying meaningful boundary conditions.In order to obtain the non-linear ODEs analytically,Galerkin method (GM) is employed.Subsequently,the ODEs are also solved by a reliable numerical solution.In order to test the accuracy,precision and reliability of the analytical method,results of the analytical analysis are compared with the numerical results.With respect to the comparisons,fairly good compatibilities with insignificant errors are observed.Eventually,the impacts of effective parameters including magnetic and radiation parameters and nanofluid volume fraction on the velocity,skin friction coefficient and Nusselt number distributions are comprehensively described.Based on the results,it is revealed that with increasing the role of magnetic force,velocity profile,skin friction coefficient and thermal performance descend.Radiation parameter has insignificant influence on velocity profile while it obviously has augmentative and decreasing effects on skin friction and Nusselt number,respectively.展开更多
The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte ...The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.展开更多
Velocities of solid phase and liquid phase in debris flow are one key problem to research on impact and abrasion mechanism of banks and control structures under action of debris flow. Debris flow was simplified as two...Velocities of solid phase and liquid phase in debris flow are one key problem to research on impact and abrasion mechanism of banks and control structures under action of debris flow. Debris flow was simplified as two-phase liquid composed of solid phase with the same diameter particles and liquid phase with the same mechanical features. Assume debris flow was one-dimension two-phase liquid moving to one direction, then general equations of velocities of solid phase and liquid phase were founded in two-phase theory. Methods to calculate average pressures, volume forces and surface forces of debris flow control volume were established. Specially, surface forces were ascertained using Bingham＇s rheology equation of liquid phase and Bagnold＇s testing results about interaction between particles of solid phase. Proportional coefficient of velocities between liquid phase and solid phase was put forward, meanwhile, divergent coefficient between theoretical velocity and real velocity of solid phase was provided too. To state succinctly before, method to calculate velocities of solid phase and liquid phase was obtained through solution to general equations. The method is suitable for both viscous debris flow and thin debris flow. Additionally, velocities every phase can be identified through analyzing deposits in-situ after occurring of debris flow. It is obvious from engineering case the result in the method is consistent to that in real-time field observation.展开更多
Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have on...Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.展开更多
Debris flows are the one type of natural disaster that is most closely associated with hu- man activities. Debris flows are characterized as being widely distributed and frequently activated. Rainfall is an important ...Debris flows are the one type of natural disaster that is most closely associated with hu- man activities. Debris flows are characterized as being widely distributed and frequently activated. Rainfall is an important component of debris flows and is the most active factor when debris flows oc- cur. Rainfall also determines the temporal and spatial distribution characteristics of the hazards. A reasonable rainfall threshold target is essential to ensuring the accuracy of debris flow pre-warning. Such a threshold is important for the study of the mechanisms of debris flow formation, predicting the characteristics of future activities and the design of prevention and engineering control measures. Most mountainous areas have little data regarding rainfall and hazards, especially in debris flow forming re- gions. Therefore, both the traditional demonstration method and frequency calculated method cannot satisfy the debris flow pre-warning requirements. This study presents the characteristics of pre-warning regions, included the rainfall, hydrologic and topographic conditions. An analogous area with abundant data and the same conditions as the pre-warning region was selected, and the rainfall threshold was calculated by proxy. This method resolved the problem of debris flow pre-warning in ar- eas lacking data and provided a new approach for debris flow pre-warning in mountainous areas.展开更多
基金co-supported by the National Natural Science Foundation of China(No.52306053)the Science Center for Gas Turbine Project,China(No.P2022-B-Ⅱ-005-001)the National Science and Technology Major Project of China(No.2017-Ⅱ-0010-0024)。
文摘With the continuous increase of aeroengine flight ceiling(>20 km),the thin atmosphere at high altitudes and the size effect all cause the compressor component inlet Reynolds number to decrease rapidly to a critical value(approximately 2.0×10^(5)),and the significant transition process on the blade/endwall surface leads to the sharp degradation of compressor performance,which seriously affects the engine fuel consumption and working stability at high altitudes.In this paper,the research progress on the internal flow mechanism and flow control methods of axial compressors at low Reynolds numbers is reviewed from the aspects of quantification and prediction of performance variation,flow loss mechanism related to separation and transition,efficient transition control and flow field organization.The development trend of the low-Reynolds-number effect of axial flow compressors is noted,and the difficulties and application prospects of aerodynamic design and efficient flow control methods for compressors under low Reynolds numbers at high altitudes are discussed.
基金supported by the National Natural Science Foundation of China(No.92252201)the Fundamental Research Funds for the Central Universitiesthe Academic Excellence Foundation of Beihang University(BUAA)for PhD Students。
文摘Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.
基金Supported by National Natural Science Foundation of China(Grant No.51975025)National Key Research and Development Program of China(Grant No.2019YFB2004500)。
文摘The flow ripple caused by an axial piston pump may lead to pipe vibrations and lower hydraulic component reliability,which are of particular concern in hydraulic systems.The valve plate of the pump is considered the part most related to flow ripple,and its structural design is an important topic.In this study,an analytical model for the axial piston pump flow ripple was established and verified using a numerical analysis with computational fluid dynamics(CFD)calculations.Moreover,a parametric analysis of the valve plate was performed to investigate the critical parameters and their ranges.A fast optimization method,the rotation vector optimization method(RVOM),was proposed for the valve plate design and compared with the currently used optimization methods to prove its efficiency.As a constant-pressure pump works in different states of swashplate angle,outlet pressure,and pump speed,an optimization principle for the entire working status was proposed to achieve the overall reduction performance.A test rig for an aircraft hydraulic pump was established,and validation experiments were conducted.It was determined that the optimized pump could achieve reduction at multiple working statuses,and the largest pressure pulsation reduction ratios for the typical speed and speed sweep tests reached 64.7%and 71.7%,respectively.The model and method proposed in this study are proven to be effective and accurate.
基金Supported by the National Natural Science Foundation of China(No.51905211)A Project of the“20 Regulations for New Universities”Funding Program of Jinan(No.202228116).
文摘Vortex-induced vibration(VIV)of an underwater manipulator in pulsating flow presents a notable engineering problem in precise control due to the velocity variation in the flow.This study investigates the VIV response of an underwater manipulator subjected to pulsating flow,focusing on how different postures affect the behavior of the system.The effects of pulsating parameters and manipulator arrangement on the hydrodynamic coefficient,vibration response,motion trajectory,and vortex shedding behaviors were analyzed.Results indicated that the cross flow vibration displacement in pulsating flow increased by 32.14%compared to uniform flow,inducing a shift in the motion trajectory from a crescent shape to a sideward vase shape.In the absence of interference between the upper and lower arms,the lift coefficient of the manipulator substantially increased with rising pulsating frequency,reaching a maximum increment of 67.0%.This increase in the lift coefficient led to a 67.05%rise in the vibration frequency of the manipulator in the in-line direction.As the pulsating amplitude increased,the drag coefficient of the underwater manipulator rose by 36.79%,but the vibration frequency in the cross-flow direction decreased by 56.26%.Additionally,when the upper and lower arms remained in a state of mutual interference,the cross-flow vibration amplitudes of the upper and lower arms were approximately 1.84 and 4.82 times higher in a circular-elliptical arrangement compared to an elliptical-circular arrangement,respectively.Consequently,the flow field shifted from a P+S pattern to a disordered pattern,disrupting the regularity of the motion trajectory.
基金Project (10964004) supported by the National Natural Science Foundation of ChinaProject (20070731001) supported by Research Fund for the Doctoral Program of ChinaProject (096RJZA104) supported by the Natural Science Foundation of Gansu Province,China
文摘The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled with flow field. The microsegregation, concentration maximum value, boundary thickness of concentration near upstream dendrite and normal to flow dendrite, and downstream dendrite were studied quantitatively in the case of forced lamia flow. The simulation results show that solute field and flow field interact complexly. Compared with melt without flow, in front of upstream dendrite tip, the concentration boundary thickness is the lowest and the concentration maximum value is the smallest for melt with flow. However, in front of downstream dendrite tip, the results are just the opposite. The zone of poor Cu in upstream dendrite where is the most severely microsegregation and shrinkage cavity is wider and the concentration is lower for melt with flow than that without flow.
基金Projects(51161011,11364024)supported by the National Natural Science Foundation of ChinaProject(1204GKCA065)supported by the Key Technology R&D Program of Gansu Province,China+1 种基金Project(201210)supported by the Fundamental Research Funds for the Universities of Gansu Province,ChinaProject(J201304)supported by the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China
文摘A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.
文摘The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.
基金Supported by the Aeronautical Science Foundation of China(20061453020)Foundation for Basic Research of Northwestern Polytechnical University(03)~~
文摘The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma：0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.
基金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.
文摘Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks. The first one is designed on the path-turn incidencerelationship, and it is similar to the computational procedure of link flows. It applies to thetraffic assignment algorithms that can provide detailed path structures. The second utilizes thelink-turn incidence relationship and the conservation of flow on links, a law deriving from thisrelationship. It is actually an improved version of Dial's logit assignment algorithm. The proposedapproaches can avoid the shortcomings both of the estimation methods, e. g. Furness's model andFrator's model, and of the network-expanding method in precision, stability and computation scale.Finally, they are validated by numerical examples.
基金Item Sponsored by National Basic Research Programof China (200002600)
文摘From the viewpoint of systems energy conservation, the influences of material flow on its energy consumption in a steel manufacturing process is an important subject. The quantitative analysis of the relationship between material flow and the energy intensity is useful to save energy in steel industry. Based on the concept of standard material flow diagram, all possible situations of ferric material flow in steel manufacturing process are analyzed. The expressions of the influence of material flow deviated from standard material flow diagram on energy consumption are put forward.
文摘In this study, the effects of magnetic field and nanoparticle on the Jeffery- Hamel flow are studied using a powerful analytical method called the Adomian decomposition method （ADM）. The traditional Navier-Stokes equation of fluid mechanics and Maxwell＇s electromagnetism governing equations are reduced to nonlinear ordinary differential equations to model the problem. The obtained results are well agreed with that of the Runge-Kutta method. The present plots confirm that the method has high accuracy for different a, Ha, and Re numbers. The flow field inside the divergent channel is studied for various values of Hartmann ：number and angle of channel. The effect of nanoparticle volume fraction in the absence of magnetic field is investigated.
基金The project supported by the National Natural Science Foundation of China (50506028)Action Scheme for Invigorating Education Towards the twenty-first century
文摘The motivation of this study is to investigate the turbulence-chemistry interactions by using probability density function （PDF） method. A consistent hybrid Reynolds Averaged Navier-Stokes （RANS）/PDF method is used to simulate the turbulent non-reacting and reacting flows. The joint fluctuating velocity-frequency-composition PDF equation coupled with the Reynolds averaged density, momentum and energy equations are solved on unstructured meshes by the Lagrangian Monte Carlo （MC） method combined with the finite volume （FV） method. The simulation of the axisymmetric bluff body stabilized non-reacting flow fields is presented in this paper. The calculated length of the recirculation zone is in good agreement with the experimental data. Moreover, the significant change of the flow pattern with the increase of the jet-to-coflow momentum flux ratio is well predicted. In addition, comparisons are made between the joint PDF model and two different Reynolds stress models.
基金Iranian Offshore Oil Company (IOOC) for financial support of this work
文摘Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.
文摘The main goal of this paper is to investigate natural convective heat transfer and flow characteristics of non-Newtonian nanofluid streaming between two infinite vertical flat plates in the presence of magnetic field and thermal radiation.Initially,a similarity transformation is used to convert momentum and energy conservation equations in partial differential forms into non-linear ordinary differential equations (ODE) applying meaningful boundary conditions.In order to obtain the non-linear ODEs analytically,Galerkin method (GM) is employed.Subsequently,the ODEs are also solved by a reliable numerical solution.In order to test the accuracy,precision and reliability of the analytical method,results of the analytical analysis are compared with the numerical results.With respect to the comparisons,fairly good compatibilities with insignificant errors are observed.Eventually,the impacts of effective parameters including magnetic and radiation parameters and nanofluid volume fraction on the velocity,skin friction coefficient and Nusselt number distributions are comprehensively described.Based on the results,it is revealed that with increasing the role of magnetic force,velocity profile,skin friction coefficient and thermal performance descend.Radiation parameter has insignificant influence on velocity profile while it obviously has augmentative and decreasing effects on skin friction and Nusselt number,respectively.
文摘The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.
基金Project supported by the Talent Fund of the Ministry of Communication of China(No.95050508) the Fund of Western Communication of China(No.200332822047) the Key Science Fund of the Ministry of Communication of China(No.95060233)
文摘Velocities of solid phase and liquid phase in debris flow are one key problem to research on impact and abrasion mechanism of banks and control structures under action of debris flow. Debris flow was simplified as two-phase liquid composed of solid phase with the same diameter particles and liquid phase with the same mechanical features. Assume debris flow was one-dimension two-phase liquid moving to one direction, then general equations of velocities of solid phase and liquid phase were founded in two-phase theory. Methods to calculate average pressures, volume forces and surface forces of debris flow control volume were established. Specially, surface forces were ascertained using Bingham＇s rheology equation of liquid phase and Bagnold＇s testing results about interaction between particles of solid phase. Proportional coefficient of velocities between liquid phase and solid phase was put forward, meanwhile, divergent coefficient between theoretical velocity and real velocity of solid phase was provided too. To state succinctly before, method to calculate velocities of solid phase and liquid phase was obtained through solution to general equations. The method is suitable for both viscous debris flow and thin debris flow. Additionally, velocities every phase can be identified through analyzing deposits in-situ after occurring of debris flow. It is obvious from engineering case the result in the method is consistent to that in real-time field observation.
基金supported by National Science Foundation of China (10576015)
文摘Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.
基金supported by the National Natural Science Foundation of China(Nos.40830742 and 40901007)
文摘Debris flows are the one type of natural disaster that is most closely associated with hu- man activities. Debris flows are characterized as being widely distributed and frequently activated. Rainfall is an important component of debris flows and is the most active factor when debris flows oc- cur. Rainfall also determines the temporal and spatial distribution characteristics of the hazards. A reasonable rainfall threshold target is essential to ensuring the accuracy of debris flow pre-warning. Such a threshold is important for the study of the mechanisms of debris flow formation, predicting the characteristics of future activities and the design of prevention and engineering control measures. Most mountainous areas have little data regarding rainfall and hazards, especially in debris flow forming re- gions. Therefore, both the traditional demonstration method and frequency calculated method cannot satisfy the debris flow pre-warning requirements. This study presents the characteristics of pre-warning regions, included the rainfall, hydrologic and topographic conditions. An analogous area with abundant data and the same conditions as the pre-warning region was selected, and the rainfall threshold was calculated by proxy. This method resolved the problem of debris flow pre-warning in ar- eas lacking data and provided a new approach for debris flow pre-warning in mountainous areas.
