Flow in tidal rivers periodically propagates upstream or downstream under tidal influence. Hydrodynamic models based on the Saint-Venant equations （the SVN model） are extensively used to model tidal rivers. A force-...Flow in tidal rivers periodically propagates upstream or downstream under tidal influence. Hydrodynamic models based on the Saint-Venant equations （the SVN model） are extensively used to model tidal rivers. A force-corrected term expressed as the combination of flow velocity and the change rate of the tidal fevel was developed to represent tidal effects in the SVN model. A momentum equation incorporating with the corrected term was derived based on Newton＇s second law. By combing the modified momentum equation with the continuity equation, an improved SVN model for tidal rivers （the ISVN model） was constructed. The simulation of a tidal reach of the Qiantang River shows that the ISVN model performs better than the SVN model. It indicates that the corrected force derived for tidal effects is reasonable; the ISVN model provides an appropriate enhancement of the SVN model for flow simulation of tidal rivers.展开更多
The metallurgical structure and composition of ingots which depend critically on the fluid motion within the molten pool during the vacuum arc remelting(VAR)process have important effect on the subsequent mechanical p...The metallurgical structure and composition of ingots which depend critically on the fluid motion within the molten pool during the vacuum arc remelting(VAR)process have important effect on the subsequent mechanical processes like forging,rolling and welding.In order to determine the fluid motion of molten pool,a 2D finite element model is established using ANSYS10.0 software,combined with the turbulent fluid flow and heat transfer.The fluid motion caused by thermo buoyancy forces is investigated at different VAR processes in the present study.The results indicate that the fluid flows symmetrically along the axis of the molten pool and clockwisely along the circle at the right pool’s profile.It is also shown that the maximum velocity increases with increasing melting rate and a direct proportional relationship exists.展开更多
Accurate prediction of the aerodynamic response of a compressor under inlet distortion is crucial for next-generation civil aircraft,such as Boundary Layer Ingestion(BLI)silent aircraft.Therefore,research on the Body ...Accurate prediction of the aerodynamic response of a compressor under inlet distortion is crucial for next-generation civil aircraft,such as Boundary Layer Ingestion(BLI)silent aircraft.Therefore,research on the Body Force(BF)model plays a significant role in achieving this objective.However,distorted inlet airflow can lead to varying operating conditions across different spatial locations of the compressor,which may cause some regions to operate outside the stability boundary.Consequently,the accuracy of BF model simulations might be compromised.To address this issue,this paper proposes a numerical simulation strategy for acquiring the steady axisymmetric three-dimensional flow field of a compressor operating at low mass flow rates,which is known as the Underlying Axisymmetric Pressure Rise Characteristic(UAPRC).The proposed simulation accounts for two different rotor speeds of a transonic compressor and identifies initial positions in the flow field where deterioration occurs based on prior experimental investigations.Moreover,simulation results are incorporated into the BF model to replicate hub instability observed in experiments.Obtained results demonstrate that this strategy provides valid predictions of the UAPRC of the compressor,thereby addressing the limitations associated with the BF model.展开更多
Based on the Fourier–Chebyshev spectral method, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with direct numerical simulation(DNS...Based on the Fourier–Chebyshev spectral method, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with direct numerical simulation(DNS) methods for different Reynolds numbers. A formula is derived to express the relation between fluctuating velocities and the friction drag coefficient. With the application of electromagnetic force, the in-depth relations among the fluctuating velocities near the wall, Reynolds stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The fluctuating velocities along the streamwise and normal directions are suppressed significantly,while the fluctuating velocity along the spanwise direction is enhanced dramatically due to the spanwise electromagnetic force. However, the values of Reynolds stress depend on the fluctuating velocities along the streamwise and normal directions rather than that along the spanwise direction. Therefore, the significant effect of drag reduction is obtained. Moreover,the maximum drag reduction is weakened due to the decay of control effect for fluctuating velocities as the Reynolds number increases.展开更多
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.展开更多
Precise function expression of the flow area for the sloping U-shape notch orifice versus the spool stroke was derived. The computational fluid dynamics was used to analyze the flow features of the sloping U-shape not...Precise function expression of the flow area for the sloping U-shape notch orifice versus the spool stroke was derived. The computational fluid dynamics was used to analyze the flow features of the sloping U-shape notch on the spool, such as mass flow rates, flow coefficients, effiux angles and steady state flow forces under different operating conditions. At last, the reliability of the mathematical model of the flow area for the sloping U-shape notch orifice on the spool was demonstrated by the comparison between the orifice area curve derived and the corresponding experimental data provided by the test. It is presented that the bottom arc of sloping U-shape notch （ABU） should not be omitted when it is required to accurately calculate the orifice area of ABU. Although the theoretical flow area of plain bottom sloping U-shape notch （PBU） is larger than that of ABU at the same opening, the simulated mass flow and experimental flow area of ABU are both larger than these of PBU at the same opening, while the simulated flow force of PBU is larger than that of ABU at the same opening. Therefore, it should be prior to adapt the ABU when designing the spool with proportional character.展开更多
With the structure of two air gaps and two rotors, the electromagnetic continuously variable transmission（EMCVT） is a novel power-split continuously variable transmission（CVT）. There are two kinds of power flowing...With the structure of two air gaps and two rotors, the electromagnetic continuously variable transmission（EMCVT） is a novel power-split continuously variable transmission（CVT）. There are two kinds of power flowing through the EMCVT, one is mechanical power and the other is electric power. In the mean time, there are three power ports in the EMCVT, one is the outer rotor named mechanical power port and the other two are the inner rotor and the stator named electric power ports. The mechanical power port is connected to the driving wheels through the final gear and the electric ports are connected to the batteries through the transducers. The two kinds of power are coupled on the outer rotor of the EMCVT. The EMCVT can be equipped on the conventional vehicle being regarded as the CVT and it also can be equipped on the hybrid electric vehicle（HEV） as the multi-energy sources assembly. The power flows of these two kinds of applications are analysed. The back electromotive force（EMF） equations are illatively studied and so the dynamic mathematic model is theorized. In order to certify the feasibility of the above theories, three simulations are carried out in allusion to the above two kinds of mentioned applications of the EMCVT and a five speed automatic transmission（AT） vehicle. The simulation results illustrate that the efficiency of the EMCVT vehicles is higher than that of the AT vehicle owed to the optimized operation area of the engine. Hence the fuel consumption of the EMCVT vehicles is knock-down.展开更多
A direct numerical simulation（DNS） is performed to investigate the control effect and mechanism of turbulent channel flow with the distribution of spanwise Lorentz force. A sinusoidal distribution of constant spanwi...A direct numerical simulation（DNS） is performed to investigate the control effect and mechanism of turbulent channel flow with the distribution of spanwise Lorentz force. A sinusoidal distribution of constant spanwise Lorentz force is selected, of which the control effects, such as flow characters, mean Reynolds stress, and drag reductions, at different parameters of amplitude A and wave number k_x are discussed. The results indicate that the control effects vary with the parameter A and k_x. With the increase of A, the drag reduction rate D_r first increases and then decreases rapidly at low k_x,and slowly at high k_x. The low drag reduction（or even drag increase） is due to a weak suppression or even the enhancements of the random velocity fluctuation and mean Reynolds stress. The efficient drag reduction is due to the quasi-streamwise vortex structure induced by Lorentz force, which contributes to suppressing the random velocity fluctuation and mean Reynolds stress, and the negative vorticity improves the distribution of streamwise velocity. Therefore, the optimal control effect with a drag reduction of up to 58% can be obtained.展开更多
To investigate the movement mechanism of debris flow, a two-dimensional, two-phase, depthintegrated model is introduced. The model uses Mohr-Coulomb plasticity for the solid rheology, and the fluid stress is modeled a...To investigate the movement mechanism of debris flow, a two-dimensional, two-phase, depthintegrated model is introduced. The model uses Mohr-Coulomb plasticity for the solid rheology, and the fluid stress is modeled as a Newtonian fluid. The interaction between solid and liquid phases, which plays a major role in debris flow movement, is assumed to consist of drag and buoyancy forces. The applicability of drag force formulas is discussed. Considering the complex interaction between debris flow and the bed surface, a combined friction boundary condition is imposed on the bottom, and this is also discussed. To solve the complex model equations, a numerical method with second-order accuracy based on the finite volume method is proposed. Several numerical experiments are performed to verify the feasibilities of model and numerical schemes. Numerical results demonstrate that different solid volume fractions substantially affect debris flow movement.展开更多
In real-time hybrid simulation（RTHS）, it is difficult if not impossible to completely erase the error in restoring force due to actuator response delay using existing displacement-based compensation methods. This pa...In real-time hybrid simulation（RTHS）, it is difficult if not impossible to completely erase the error in restoring force due to actuator response delay using existing displacement-based compensation methods. This paper proposes a new force correction method based on online discrete tangent stiffness estimation（online DTSE） to provide accurate online estimation of the instantaneous stiffness of the physical substructure. Following the discrete curve parameter recognition theory, the online DTSE method estimates the instantaneous stiffness mainly through adaptively building a fuzzy segment with the latest measurements, constructing several strict bounding lines of the segment and calculating the slope of the strict bounding lines, which significantly improves the calculation efficiency and accuracy for the instantaneous stiffness estimation. The results of both computational simulation and real-time hybrid simulation show that：（1） the online DTSE method has high calculation efficiency, of which the relatively short computation time will not interrupt RTHS; and（2） the online DTSE method provides better estimation for the instantaneous stiffness, compared with other existing estimation methods. Due to the quick and accurate estimation of instantaneous stiffness, the online DTSE method therefore provides a promising technique to correct restoring forces in RTHS.展开更多
Based on the commercial computational fluid dynamics software CFX-4.3, electrolyte flow fields in a 156 kA pre-baked anode aluminum electrolysis cell were investigated in three different cases where the electrolyte me...Based on the commercial computational fluid dynamics software CFX-4.3, electrolyte flow fields in a 156 kA pre-baked anode aluminum electrolysis cell were investigated in three different cases where the electrolyte melt was driven by different kinds of force, i.e. electromagnetic force only, the anode gas drag force only and both of the former two forces. The results show that when electromagnetic force was introduced only, most of the electrolyte moves horizontally; when anode gas drag force was introduced only, the electrolyte flows mainly around each anode with small circulation; when electromagnetic force and anode gas drag force were both introduced together, the structure of the electrolyte flow fields and the velocity of electrolyte are similar to that of the case where only anode gas drag force is used. The electrolyte flow fields are mainly determined by the anode gas drag force.展开更多
In the course of studying on aerodynamic change and its effect on steering stability and controllability of an automobile in passing, because of multi interaction streams, it is difficult to use traditional methods, s...In the course of studying on aerodynamic change and its effect on steering stability and controllability of an automobile in passing, because of multi interaction streams, it is difficult to use traditional methods, such as wind tunnel test and road test. If the passing process of an automobile is divided into many time segments, so as to avoid the use of moving mesh which takes large calculation resource and CPU processing time in calculating, the segments are simulated with computational fluid dynamics （CFD） method, then the approximate computational results about external flow field will be obtained. On the basis of the idea, the change of external flow field of wagon-shaped car at the moment of passing is simulated through solving three-dimensional, steady and uncompressible N-S equations with finite volume method. Numerical simulation analysis of side force coefficient, stream lines, body surface pressure distribution of wagon-shaped car are presented and a preliminary discussion of aerodynamic characteristics of correlative situations is obtained. Finally, the C3 -x/l curve of side force coefficient（C3） of car following relative distance （x/l） between cars is obtained. By comparison, the curve is coincident well with the experimental data, which shows creditability of numerical simulation methods presented.展开更多
A three-dimensional mathematical and physical model coupling with the heat transfer and the flow of molten metal in the centrifugal casting of the high speed steel roll was established by using CFD software FLUENT. It...A three-dimensional mathematical and physical model coupling with the heat transfer and the flow of molten metal in the centrifugal casting of the high speed steel roll was established by using CFD software FLUENT. It can be used to analyze the distribution of the temperature filed and the flow filed in the centrifugal casting under the gravity, the electromagnetic stirring force and the centrifugal force. Some experiments were carried out to verify the above analysis results. The effects of the electromagnetic force on the centrifugal casting process are discussed. The results showed that under the 0.15 T electromagnetic field intensity, both the absolute pressure of metal flow to mold wall and the metal flow velocity on the same location have some differences between the electromagnetic centrifugal casting and the centrifugal casting. Numerical results for understanding the electromagnetic stirring of the centrifugal casting process have a guiding significance.展开更多
Numerical simulations and experimental research are both carried out to investigate the controlled effect of spanwise oscillating Lorentz force on a turbulent channel flow. The variations of the streaks and the skin f...Numerical simulations and experimental research are both carried out to investigate the controlled effect of spanwise oscillating Lorentz force on a turbulent channel flow. The variations of the streaks and the skin friction drag are obtained through the PIV system and the drag measurement system, respectively. The flow field in the near-wall region is shown through direct numerical simulations utilizing spectral method. The experimental results are consistent with the numerical simulation results qualitatively, and both the results indicate that the streaks are tilted into the spanwise direction and the drag reduction utilizing spanwise oscillating Lorentz forces can be realized. The numerical simulation results reveal more detail of the drag reduction mechanism which can be explained, since the spanwise vorticity generated from the interaction between the induced Stokes layer and intrinsic turbulent flow in the near-wall region can make the longitudinal vortices tilt and oscillate, and leads to turbulence suppression and drag reduction.展开更多
A formulation of the skin-friction drag related to the Reynolds shear stress in a turbulent channel flow is derived. A direct numerical simulation （DNS） of the turbulent control is performed by imposing the spatiall...A formulation of the skin-friction drag related to the Reynolds shear stress in a turbulent channel flow is derived. A direct numerical simulation （DNS） of the turbulent control is performed by imposing the spatially oscillating spanwise Lorentz force. Under the action of the Lorentz force with several proper control parameters, only the periodi- cally well-organized streamwise vortices are finally observed in the near-wall region. The Reynolds shear stress decreases dramatically, especially in the near-wall area, resulting in a drag reduction.展开更多
The ejector-powered engine simulator(EPES)system is an important piece of equipment in conducting an influence test of the intake and jet flow in low-speed wind tunnels.In this work,through the analysis of the structu...The ejector-powered engine simulator(EPES)system is an important piece of equipment in conducting an influence test of the intake and jet flow in low-speed wind tunnels.In this work,through the analysis of the structure and principle of EPES,three parts of the internal flow force were obtained,namely,the additional resistance before the inlet,the internal flow force in the inlet and the thrust produced by the ejector.On the assumption of one-dimensional isentropic adiabatic flow,the theoretical formulae for calculating the forces were derived according to the measured total pressure,static pressure and total temperature of the internal flow section.Subsequently,a calibration tank was used to calibrate the EPES system.On the basis of the characteristics of the EPES system,the process and method of its calibration were designed in detail,and the model installation interface of the calibration tank was reformed.By applying this method,the repeatability accuracy of the inlet flow rate calibration coefficient was less than0.05%,whereas that of the exhaust flow rate and velocity was less than 0.1%.Upon the application of the calibration coefficients to the correction of the wind tunnel experiment data,the results showed good agreement with the numerical simulation results in terms of regularity and magnitude before stall,which validates the reasonableness and feasibility of the calibration method.Analysis of the calibration data also demonstrated the consistency in the variation law and trend between the theoretical calculation and actual measurement of internal flow force,further reflecting the rationality and feasibility of the theoretical calculation.Nevertheless,the numerical difference was large and further widened with a higher ejection flow rate mainly because of the accuracy of flow measurement and the inhomogeneity of internal flow.The thrust deflection angle of EPES is an important factor in correcting this issue.In particular,the thrust deflection angle becomes larger with small ejection flow and becomes smaller with an increase in flow rate,essentially exhibiting a general change of less than 10°.展开更多
Nowadays,more and more attention has been paid to improve the performance of the nozzle flapper servo valve.As a core part of nozzle flapper servo valve,the armature assembly is affected by electromagnetic force,jet f...Nowadays,more and more attention has been paid to improve the performance of the nozzle flapper servo valve.As a core part of nozzle flapper servo valve,the armature assembly is affected by electromagnetic force,jet force and feedback force at the same time.Due to the complex structure of the pilot stage flow field and the high jet pressure,the prediction of the jet force has always been difficult in modeling the transient motion of the servo valve.Whereupon,a numerical simulation method based on the flow-solid interaction(FSI)is applied to observe the variation of the jet force when the flapper is moving.Different parameters are employed to seek a suitable numerical simulation model which can balance the accuracy and computational cost.By comparing with the experiment results,the effectiveness of numerical simulation method in predicting the variation of the jet force and cavitation is verified.By this numerical simulation model,the distribution of flow field and the force on the flapper predicted by the moving and fixed flapper are compared.The results show that more dynamic details are achieved by the transient simulation.By analyzing the numerical simulation results of different inlet pressures and flapper vibration frequencies,the relationship between the movement of the flapper,the flow field distribution,the jet force and the inlet pressure is established,which provides a theoretical basis for the subsequent modeling of the armature assembly.展开更多
In the design of rock sheds for the mitigation of risk due to rapid and long landslides, a crucial role is played by the evaluation of the impact force exerted by the flowing mass on the rock sheds. This paper is focu...In the design of rock sheds for the mitigation of risk due to rapid and long landslides, a crucial role is played by the evaluation of the impact force exerted by the flowing mass on the rock sheds. This paper is focused on the influencing factors of the impact force of dry granular flow onto rock shed and in particular on the evaluation of the maximum impact force. The coupled DEM-FEM model calibrated with small-scale physical experiment is used to simulate the movement of dry granular flow coupled with impact forces on the rock-shed. Based on the numerical results, three key stages were identified of impact process, namely startup streams slippery, impact and pile-up. The maximum impact force increases linearly with bulk density, and the maximum impact force exhibits a power law dependence on the impact height and slop angle respectively. The sensitivities of bulk density, impact height, and slope angle on the maximum impact force are: 1.0, 0.496, and 2.32 respectively in the benchmark model. The parameters with high sensitivity should be given priority in the design of the rock shed. The results obtained from this study are useful for facilitating design of shed against dry granular flow.展开更多
This paper discusses numerical results from three-dimensional large eddy simulations of an oscillating cylinder under prescribed movements in uniform flow. Six cases, namely pure in-line, pure cross-flow and two group...This paper discusses numerical results from three-dimensional large eddy simulations of an oscillating cylinder under prescribed movements in uniform flow. Six cases, namely pure in-line, pure cross-flow and two groups of ＇Figure of Eight＇ oscillation patterns are under investigation at Reynolds number Re = 24000. The ＇ Figure of Eight＇ pattern in each group is with identical shape but oppusite orbital directions. The numerical results on hydrodynamic forces, higher order force components, and vortex shedding modes are extensively studied and compared with the measured experimental data. It is found that the fluid force in phase with the velocity, which represents the energy transfer between the fluid and the cylinder, has opposite sign and different magnitude due to the opposite orbital direction. Higher order force components in cross-flow direction are found to occur at odd nmnber times of the oscillating frequency, while even nmbers dominate the higher order force components in in-llne direction. The 2C and 2T vortex shedding modes are well reproduced due to the opposite orbital direction effect. Comparisons between numerical and experimental results indicate that the present numerical model could be a rational tool for the identification of hydrodynamic coefficients which are normally applied in empirical models to predict the vortex-induced vibrations of slender marine structures.展开更多
基金supported by the National Key Technologies R&D Program of China for the Eleventh Five-Year Plan Period (Grant No. 2008BAB29B08-02)the Program for the Ministry of Education and State Administration of Foreign Experts Affairs of China (Grant No. B08408)
文摘Flow in tidal rivers periodically propagates upstream or downstream under tidal influence. Hydrodynamic models based on the Saint-Venant equations （the SVN model） are extensively used to model tidal rivers. A force-corrected term expressed as the combination of flow velocity and the change rate of the tidal fevel was developed to represent tidal effects in the SVN model. A momentum equation incorporating with the corrected term was derived based on Newton＇s second law. By combing the modified momentum equation with the continuity equation, an improved SVN model for tidal rivers （the ISVN model） was constructed. The simulation of a tidal reach of the Qiantang River shows that the ISVN model performs better than the SVN model. It indicates that the corrected force derived for tidal effects is reasonable; the ISVN model provides an appropriate enhancement of the SVN model for flow simulation of tidal rivers.
基金the National Basic Research Program (973) of China(No.2007CB613802)
文摘The metallurgical structure and composition of ingots which depend critically on the fluid motion within the molten pool during the vacuum arc remelting(VAR)process have important effect on the subsequent mechanical processes like forging,rolling and welding.In order to determine the fluid motion of molten pool,a 2D finite element model is established using ANSYS10.0 software,combined with the turbulent fluid flow and heat transfer.The fluid motion caused by thermo buoyancy forces is investigated at different VAR processes in the present study.The results indicate that the fluid flows symmetrically along the axis of the molten pool and clockwisely along the circle at the right pool’s profile.It is also shown that the maximum velocity increases with increasing melting rate and a direct proportional relationship exists.
基金the National Natural Science Foundation of China(Nos.52322603 and 51976005)the Science Center for Gas Turbine Project,China(Nos.P2022-B-II-004-001 and P2023-B-II-001-001)the Fundamental Research Funds for the Central Universities,and Beijing Nova Program,China(Nos.20220484074 and 20230484479).
文摘Accurate prediction of the aerodynamic response of a compressor under inlet distortion is crucial for next-generation civil aircraft,such as Boundary Layer Ingestion(BLI)silent aircraft.Therefore,research on the Body Force(BF)model plays a significant role in achieving this objective.However,distorted inlet airflow can lead to varying operating conditions across different spatial locations of the compressor,which may cause some regions to operate outside the stability boundary.Consequently,the accuracy of BF model simulations might be compromised.To address this issue,this paper proposes a numerical simulation strategy for acquiring the steady axisymmetric three-dimensional flow field of a compressor operating at low mass flow rates,which is known as the Underlying Axisymmetric Pressure Rise Characteristic(UAPRC).The proposed simulation accounts for two different rotor speeds of a transonic compressor and identifies initial positions in the flow field where deterioration occurs based on prior experimental investigations.Moreover,simulation results are incorporated into the BF model to replicate hub instability observed in experiments.Obtained results demonstrate that this strategy provides valid predictions of the UAPRC of the compressor,thereby addressing the limitations associated with the BF model.
基金Project supported by the National Natural Science Foundation of China(Grant No.11672135)a Foundation for the Author of National Excellent Doctoral Dissertation of China(Grant No.201461)
文摘Based on the Fourier–Chebyshev spectral method, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with direct numerical simulation(DNS) methods for different Reynolds numbers. A formula is derived to express the relation between fluctuating velocities and the friction drag coefficient. With the application of electromagnetic force, the in-depth relations among the fluctuating velocities near the wall, Reynolds stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The fluctuating velocities along the streamwise and normal directions are suppressed significantly,while the fluctuating velocity along the spanwise direction is enhanced dramatically due to the spanwise electromagnetic force. However, the values of Reynolds stress depend on the fluctuating velocities along the streamwise and normal directions rather than that along the spanwise direction. Therefore, the significant effect of drag reduction is obtained. Moreover,the maximum drag reduction is weakened due to the decay of control effect for fluctuating velocities as the Reynolds number increases.
基金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.
基金Project(51004085)supported by the National Natural Science Foundation of China
文摘Precise function expression of the flow area for the sloping U-shape notch orifice versus the spool stroke was derived. The computational fluid dynamics was used to analyze the flow features of the sloping U-shape notch on the spool, such as mass flow rates, flow coefficients, effiux angles and steady state flow forces under different operating conditions. At last, the reliability of the mathematical model of the flow area for the sloping U-shape notch orifice on the spool was demonstrated by the comparison between the orifice area curve derived and the corresponding experimental data provided by the test. It is presented that the bottom arc of sloping U-shape notch （ABU） should not be omitted when it is required to accurately calculate the orifice area of ABU. Although the theoretical flow area of plain bottom sloping U-shape notch （PBU） is larger than that of ABU at the same opening, the simulated mass flow and experimental flow area of ABU are both larger than these of PBU at the same opening, while the simulated flow force of PBU is larger than that of ABU at the same opening. Therefore, it should be prior to adapt the ABU when designing the spool with proportional character.
基金supported by National Natural Science Foundation of China(No.50605020)Guangdong Provincial Science and Technology Project of China(No.2006A10501001).
文摘With the structure of two air gaps and two rotors, the electromagnetic continuously variable transmission（EMCVT） is a novel power-split continuously variable transmission（CVT）. There are two kinds of power flowing through the EMCVT, one is mechanical power and the other is electric power. In the mean time, there are three power ports in the EMCVT, one is the outer rotor named mechanical power port and the other two are the inner rotor and the stator named electric power ports. The mechanical power port is connected to the driving wheels through the final gear and the electric ports are connected to the batteries through the transducers. The two kinds of power are coupled on the outer rotor of the EMCVT. The EMCVT can be equipped on the conventional vehicle being regarded as the CVT and it also can be equipped on the hybrid electric vehicle（HEV） as the multi-energy sources assembly. The power flows of these two kinds of applications are analysed. The back electromotive force（EMF） equations are illatively studied and so the dynamic mathematic model is theorized. In order to certify the feasibility of the above theories, three simulations are carried out in allusion to the above two kinds of mentioned applications of the EMCVT and a five speed automatic transmission（AT） vehicle. The simulation results illustrate that the efficiency of the EMCVT vehicles is higher than that of the AT vehicle owed to the optimized operation area of the engine. Hence the fuel consumption of the EMCVT vehicles is knock-down.
基金supported by the National Natural Science Foundation of China(Grant Nos.11672135 and 11202102)the Fundamental Research Funds for the Central Universities,China(Grant No.30916011347)a Foundation for the Author of National Excellent Doctoral Dissertation,China(Grant No.201461)
文摘A direct numerical simulation（DNS） is performed to investigate the control effect and mechanism of turbulent channel flow with the distribution of spanwise Lorentz force. A sinusoidal distribution of constant spanwise Lorentz force is selected, of which the control effects, such as flow characters, mean Reynolds stress, and drag reductions, at different parameters of amplitude A and wave number k_x are discussed. The results indicate that the control effects vary with the parameter A and k_x. With the increase of A, the drag reduction rate D_r first increases and then decreases rapidly at low k_x,and slowly at high k_x. The low drag reduction（or even drag increase） is due to a weak suppression or even the enhancements of the random velocity fluctuation and mean Reynolds stress. The efficient drag reduction is due to the quasi-streamwise vortex structure induced by Lorentz force, which contributes to suppressing the random velocity fluctuation and mean Reynolds stress, and the negative vorticity improves the distribution of streamwise velocity. Therefore, the optimal control effect with a drag reduction of up to 58% can be obtained.
基金Financial support from the NSFC-ICIMOD(41661144041)NSFC(Grant No.41772312)+1 种基金Key Research and Development Program(2017SZ0041)Sichuan Province Science and Technology Support Project(2016SZ0067)
文摘To investigate the movement mechanism of debris flow, a two-dimensional, two-phase, depthintegrated model is introduced. The model uses Mohr-Coulomb plasticity for the solid rheology, and the fluid stress is modeled as a Newtonian fluid. The interaction between solid and liquid phases, which plays a major role in debris flow movement, is assumed to consist of drag and buoyancy forces. The applicability of drag force formulas is discussed. Considering the complex interaction between debris flow and the bed surface, a combined friction boundary condition is imposed on the bottom, and this is also discussed. To solve the complex model equations, a numerical method with second-order accuracy based on the finite volume method is proposed. Several numerical experiments are performed to verify the feasibilities of model and numerical schemes. Numerical results demonstrate that different solid volume fractions substantially affect debris flow movement.
基金Priority Academic Program Development of Jiangsu Higher Education Institutions under Grant No.1105007002National Natural Science Foundation of China under Grant No.51378107 and No.51678147
文摘In real-time hybrid simulation（RTHS）, it is difficult if not impossible to completely erase the error in restoring force due to actuator response delay using existing displacement-based compensation methods. This paper proposes a new force correction method based on online discrete tangent stiffness estimation（online DTSE） to provide accurate online estimation of the instantaneous stiffness of the physical substructure. Following the discrete curve parameter recognition theory, the online DTSE method estimates the instantaneous stiffness mainly through adaptively building a fuzzy segment with the latest measurements, constructing several strict bounding lines of the segment and calculating the slope of the strict bounding lines, which significantly improves the calculation efficiency and accuracy for the instantaneous stiffness estimation. The results of both computational simulation and real-time hybrid simulation show that：（1） the online DTSE method has high calculation efficiency, of which the relatively short computation time will not interrupt RTHS; and（2） the online DTSE method provides better estimation for the instantaneous stiffness, compared with other existing estimation methods. Due to the quick and accurate estimation of instantaneous stiffness, the online DTSE method therefore provides a promising technique to correct restoring forces in RTHS.
基金Project (G1999064903) supported by the National Key Fundamental Research and Development Programof China
文摘Based on the commercial computational fluid dynamics software CFX-4.3, electrolyte flow fields in a 156 kA pre-baked anode aluminum electrolysis cell were investigated in three different cases where the electrolyte melt was driven by different kinds of force, i.e. electromagnetic force only, the anode gas drag force only and both of the former two forces. The results show that when electromagnetic force was introduced only, most of the electrolyte moves horizontally; when anode gas drag force was introduced only, the electrolyte flows mainly around each anode with small circulation; when electromagnetic force and anode gas drag force were both introduced together, the structure of the electrolyte flow fields and the velocity of electrolyte are similar to that of the case where only anode gas drag force is used. The electrolyte flow fields are mainly determined by the anode gas drag force.
基金National Natural Science Foundation of China(No. 50275052).
文摘In the course of studying on aerodynamic change and its effect on steering stability and controllability of an automobile in passing, because of multi interaction streams, it is difficult to use traditional methods, such as wind tunnel test and road test. If the passing process of an automobile is divided into many time segments, so as to avoid the use of moving mesh which takes large calculation resource and CPU processing time in calculating, the segments are simulated with computational fluid dynamics （CFD） method, then the approximate computational results about external flow field will be obtained. On the basis of the idea, the change of external flow field of wagon-shaped car at the moment of passing is simulated through solving three-dimensional, steady and uncompressible N-S equations with finite volume method. Numerical simulation analysis of side force coefficient, stream lines, body surface pressure distribution of wagon-shaped car are presented and a preliminary discussion of aerodynamic characteristics of correlative situations is obtained. Finally, the C3 -x/l curve of side force coefficient（C3） of car following relative distance （x/l） between cars is obtained. By comparison, the curve is coincident well with the experimental data, which shows creditability of numerical simulation methods presented.
文摘A three-dimensional mathematical and physical model coupling with the heat transfer and the flow of molten metal in the centrifugal casting of the high speed steel roll was established by using CFD software FLUENT. It can be used to analyze the distribution of the temperature filed and the flow filed in the centrifugal casting under the gravity, the electromagnetic stirring force and the centrifugal force. Some experiments were carried out to verify the above analysis results. The effects of the electromagnetic force on the centrifugal casting process are discussed. The results showed that under the 0.15 T electromagnetic field intensity, both the absolute pressure of metal flow to mold wall and the metal flow velocity on the same location have some differences between the electromagnetic centrifugal casting and the centrifugal casting. Numerical results for understanding the electromagnetic stirring of the centrifugal casting process have a guiding significance.
文摘Numerical simulations and experimental research are both carried out to investigate the controlled effect of spanwise oscillating Lorentz force on a turbulent channel flow. The variations of the streaks and the skin friction drag are obtained through the PIV system and the drag measurement system, respectively. The flow field in the near-wall region is shown through direct numerical simulations utilizing spectral method. The experimental results are consistent with the numerical simulation results qualitatively, and both the results indicate that the streaks are tilted into the spanwise direction and the drag reduction utilizing spanwise oscillating Lorentz forces can be realized. The numerical simulation results reveal more detail of the drag reduction mechanism which can be explained, since the spanwise vorticity generated from the interaction between the induced Stokes layer and intrinsic turbulent flow in the near-wall region can make the longitudinal vortices tilt and oscillate, and leads to turbulence suppression and drag reduction.
基金Supported by the National Natural Science Foundation of China(Nos.11172140 and 11372356)the Open Project of State Key Laboratory of Explosion Science and Technology in Beijing Institute of Technology(No.KFJJ13-3M)
文摘A formulation of the skin-friction drag related to the Reynolds shear stress in a turbulent channel flow is derived. A direct numerical simulation （DNS） of the turbulent control is performed by imposing the spatially oscillating spanwise Lorentz force. Under the action of the Lorentz force with several proper control parameters, only the periodi- cally well-organized streamwise vortices are finally observed in the near-wall region. The Reynolds shear stress decreases dramatically, especially in the near-wall area, resulting in a drag reduction.
基金supported by the funda-mental research the Funds of China Aerodynamics Research and Development Center
文摘The ejector-powered engine simulator(EPES)system is an important piece of equipment in conducting an influence test of the intake and jet flow in low-speed wind tunnels.In this work,through the analysis of the structure and principle of EPES,three parts of the internal flow force were obtained,namely,the additional resistance before the inlet,the internal flow force in the inlet and the thrust produced by the ejector.On the assumption of one-dimensional isentropic adiabatic flow,the theoretical formulae for calculating the forces were derived according to the measured total pressure,static pressure and total temperature of the internal flow section.Subsequently,a calibration tank was used to calibrate the EPES system.On the basis of the characteristics of the EPES system,the process and method of its calibration were designed in detail,and the model installation interface of the calibration tank was reformed.By applying this method,the repeatability accuracy of the inlet flow rate calibration coefficient was less than0.05%,whereas that of the exhaust flow rate and velocity was less than 0.1%.Upon the application of the calibration coefficients to the correction of the wind tunnel experiment data,the results showed good agreement with the numerical simulation results in terms of regularity and magnitude before stall,which validates the reasonableness and feasibility of the calibration method.Analysis of the calibration data also demonstrated the consistency in the variation law and trend between the theoretical calculation and actual measurement of internal flow force,further reflecting the rationality and feasibility of the theoretical calculation.Nevertheless,the numerical difference was large and further widened with a higher ejection flow rate mainly because of the accuracy of flow measurement and the inhomogeneity of internal flow.The thrust deflection angle of EPES is an important factor in correcting this issue.In particular,the thrust deflection angle becomes larger with small ejection flow and becomes smaller with an increase in flow rate,essentially exhibiting a general change of less than 10°.
基金Supported by the National Natural Science Foundation of China(51675119)。
文摘Nowadays,more and more attention has been paid to improve the performance of the nozzle flapper servo valve.As a core part of nozzle flapper servo valve,the armature assembly is affected by electromagnetic force,jet force and feedback force at the same time.Due to the complex structure of the pilot stage flow field and the high jet pressure,the prediction of the jet force has always been difficult in modeling the transient motion of the servo valve.Whereupon,a numerical simulation method based on the flow-solid interaction(FSI)is applied to observe the variation of the jet force when the flapper is moving.Different parameters are employed to seek a suitable numerical simulation model which can balance the accuracy and computational cost.By comparing with the experiment results,the effectiveness of numerical simulation method in predicting the variation of the jet force and cavitation is verified.By this numerical simulation model,the distribution of flow field and the force on the flapper predicted by the moving and fixed flapper are compared.The results show that more dynamic details are achieved by the transient simulation.By analyzing the numerical simulation results of different inlet pressures and flapper vibration frequencies,the relationship between the movement of the flapper,the flow field distribution,the jet force and the inlet pressure is established,which provides a theoretical basis for the subsequent modeling of the armature assembly.
文摘In the design of rock sheds for the mitigation of risk due to rapid and long landslides, a crucial role is played by the evaluation of the impact force exerted by the flowing mass on the rock sheds. This paper is focused on the influencing factors of the impact force of dry granular flow onto rock shed and in particular on the evaluation of the maximum impact force. The coupled DEM-FEM model calibrated with small-scale physical experiment is used to simulate the movement of dry granular flow coupled with impact forces on the rock-shed. Based on the numerical results, three key stages were identified of impact process, namely startup streams slippery, impact and pile-up. The maximum impact force increases linearly with bulk density, and the maximum impact force exhibits a power law dependence on the impact height and slop angle respectively. The sensitivities of bulk density, impact height, and slope angle on the maximum impact force are: 1.0, 0.496, and 2.32 respectively in the benchmark model. The parameters with high sensitivity should be given priority in the design of the rock shed. The results obtained from this study are useful for facilitating design of shed against dry granular flow.
文摘This paper discusses numerical results from three-dimensional large eddy simulations of an oscillating cylinder under prescribed movements in uniform flow. Six cases, namely pure in-line, pure cross-flow and two groups of ＇Figure of Eight＇ oscillation patterns are under investigation at Reynolds number Re = 24000. The ＇ Figure of Eight＇ pattern in each group is with identical shape but oppusite orbital directions. The numerical results on hydrodynamic forces, higher order force components, and vortex shedding modes are extensively studied and compared with the measured experimental data. It is found that the fluid force in phase with the velocity, which represents the energy transfer between the fluid and the cylinder, has opposite sign and different magnitude due to the opposite orbital direction. Higher order force components in cross-flow direction are found to occur at odd nmnber times of the oscillating frequency, while even nmbers dominate the higher order force components in in-llne direction. The 2C and 2T vortex shedding modes are well reproduced due to the opposite orbital direction effect. Comparisons between numerical and experimental results indicate that the present numerical model could be a rational tool for the identification of hydrodynamic coefficients which are normally applied in empirical models to predict the vortex-induced vibrations of slender marine structures.
