The sealing performance of magnetic fluid is related to the magnetic fluid itself. Many factors can influence the magnetic field and the seal pressure differences of magnetic fluid seals, such as the sealing gap, the ...The sealing performance of magnetic fluid is related to the magnetic fluid itself. Many factors can influence the magnetic field and the seal pressure differences of magnetic fluid seals, such as the sealing gap, the shaft eccentricity, the shaft diameter, the volume of the magnetic fluid and the centrifugal force. These factors are analyzed by numerical computation. When the seal material and structure are the same, the seal pressure difference is directly proportional to the magnetic field intensity and the saturation magnetization of the magnetic fluid. The sealing performance of the magnetic fluid will reduce with the increase of the sealing gap and shaft eccentricity. The sealing performance will increase with the volume of the magnetic fluid and decrease with the increase of the shaft diameter taking gravity into account. The increase of the shaft diameter is the same as the reduction of the volume of the magnetic fluid. The magnetic fluid cross-section can change because of the centrifugal force. Some improvements can reducc the influence of the centrifugal force. The centrifugal force can be utilized to improve the sealing performance.展开更多
In this paper, two novel numerical computation methods are introduced which have been recently developed at Research Institute for Applied Mechanics ( R/AM ), Kyushu University, for strongly nonlinear wave-body inte...In this paper, two novel numerical computation methods are introduced which have been recently developed at Research Institute for Applied Mechanics ( R/AM ), Kyushu University, for strongly nonlinear wave-body interaction problems, such as ship motions in rough seas and resulting green-water impact on deck. The first method is the CIP-based Cartesian grid method, in which the free surface flow is treated as a multi-phase flow which is solved using a Cartesian grid. The second method is the MPS method, which is a so-called particle method and hence no grid is used. The features and calculation procedures of these numerical methods are described. One validation computation against a newly conducted experiment on a dam break problem, which is also described in this paper, is presented.展开更多
The rubber circular plate is considered as a kind of membrane. Based on the character that there exists no bending moment inside a membrane, the geometric behavior of the rubber circular plate in expanding state was d...The rubber circular plate is considered as a kind of membrane. Based on the character that there exists no bending moment inside a membrane, the geometric behavior of the rubber circular plate in expanding state was described with the aid of a group of mathematic method. The relationship between deflection and load was attained by means of calculating stress and strain inside the curved surface of rubber plate. Meantime, based on Hencky method, the relationship between deflection and load was attained and considered as the Hencky solution. The different results given rise by the two different resolving methods were compared. The deviation results from the Hencky method was discussed, and a kind of correcting method was put forward.展开更多
To study the grouting reinforcement mechanism in jointed rock slope, first, the theoretical deduction was done to calculate the critical length of slipping if the slope angle is larger than that of joint inclination; ...To study the grouting reinforcement mechanism in jointed rock slope, first, the theoretical deduction was done to calculate the critical length of slipping if the slope angle is larger than that of joint inclination; Second, the numerical calculation model was founded by FLAG^3D, so as to find the stress and deformation responses of rock mass in the state before and after grouting, the analysis results show that the range between the boundary of critical slipping block and the joint plane that passes the slope toe is the effective grouting area (EGA). After excavation, large deformation occurs along the joint plane. After grouting, the displacements of rock particles become uniform and continuous, and large deformations along the joint plane are controlled; the dynamic displacement can re- flect the deformation response of slope during excavation in the state before and after grouting, as well as the shear location of potential slip plane. After grouting, the dynamic displacement of each monitoring point reaches the peak value with very few time steps, which indicate that the parameters of the joint plane, such as strength and stiffness, are improved; the stress field becomes uniform. Tensile area reduces gradually; whole stability of the slope and its ability to resist tensile and shear stress are improved greatly.展开更多
Nowadays,there are many studies conducted in the field of marine hydrodynamics which focus on two vessels traveling and floating in sufficiently close proximity to experience significant interactions.The hydrodynamic ...Nowadays,there are many studies conducted in the field of marine hydrodynamics which focus on two vessels traveling and floating in sufficiently close proximity to experience significant interactions.The hydrodynamic behavior of parallel moving ships in waves is an interesting and important topic of late.A numerical investigation has been carried out for the prediction of wave exciting forces and motion responses of parallel moving ships in regular waves.The numerical solution was based on 3D distribution technique and using the linear wave theory to determine the exciting forces and ship's motion.The speed effects have been considered in the Green function for more realistic results.The numerical computations of wave exciting forces and motion responses were carried out for a Mariner and Series 60 for the purpose of discovering different Froude numbers and different separation distances in head sea conditions.Based on the numerical computations,it was revealed that the sway,roll and yaw have a significant effect due to hydrodynamic interaction.展开更多
As a new type of structure which has never been built, submerged floating tunnel was studied mainly by numerical simulations. To further study the seismic response of a submerged floating tunnel, the first model exper...As a new type of structure which has never been built, submerged floating tunnel was studied mainly by numerical simulations. To further study the seismic response of a submerged floating tunnel, the first model experiment of submerged floating tunnel (SFT) under the earthquake was carried out on the unique underwater shaking table in China. The experimental results show that vertical excitation induces larger response than horizontal and different inclination degrees of the tether also cause different seismic responses. Subsequently, based on the fluid-structure interaction theory, the corresponding numerical model is established. And comparing the numerical results with the experimental results, those of shaking table test. Numerical model adopted is effective for it is shown that the numerical results are basically identical with dynamic response of SFT.展开更多
Experiments involving a sonar platform with a sound absorption wedge were carried out for the purpose of obtaining the low frequency acoustic characteristics. Acoustic characteristics of a sonar platform model with a ...Experiments involving a sonar platform with a sound absorption wedge were carried out for the purpose of obtaining the low frequency acoustic characteristics. Acoustic characteristics of a sonar platform model with a sound absorption wedge were measured, and the effects of different wedge laid areas on platform acoustic characteristic were tested. Vibration acceleration and self-noise caused by model vibration were measured in four conditions: 0%, 36%, 60%, and 100% of wedge laid area when the sonar platform was under a single frequency excitation force. An experiment was performed to validate a corresponding numerical calculation. The numerical vibration characteristics of platform area were calculated by the finite element method, and self-noise caused by the vibration in it was predicted by an experiential formula. The conclusions prove that the numerical calculation method can partially replace the experimental process for obtaining vibration and sound characteristics.展开更多
We determi:ned a suitable gate road layout in slice mining in an ultra-thick unstable coal seam, using theoretical anallysis and numerical calculations. Based on plasticity theory in terms of limiting equilibrium, th...We determi:ned a suitable gate road layout in slice mining in an ultra-thick unstable coal seam, using theoretical anallysis and numerical calculations. Based on plasticity theory in terms of limiting equilibrium, the width of chain pillar in the upper slice was calculated to be 18 m. The stress distribution in the chain pillar after the upper slice was mined out was described with numerical simulation. The extent of the effect of stress on the upper chain pillar on the lower solid coal was obtained on the basis of an elastic solution of a distributed force loaded on a half-plane. Three layout designs for lower gate roads were pro- posed and a stability factor was introduced to analyze the stability of the lower pillar with numerical calculation. Gate road translation was determined as the most suitable layout method, which maximizes the extraction rate on the basis of the pillar stability.展开更多
A numerical model to predict the temperature field in a hermetic reciprocating compressor for household refrigeration appliances is presented in this work. The model combines a high resolution three-dimensional heat c...A numerical model to predict the temperature field in a hermetic reciprocating compressor for household refrigeration appliances is presented in this work. The model combines a high resolution three-dimensional heat conduction formulation of the compressor's solid parts, a three-dimensional CFD (computational fluid dynamics) approach for the gas line domain and lumped formulations of the shell gas and the lubrication oil. Heat transfer coefficients are determined by applying CFD to the gas line side and correlations from the literature on the shell gas and oil side, respectively. The valve in the gas line simulation is modelled as a parallel moving fiat plate. By means of an iterative loop the temperature field of the solid parts acts as boundary condition for the CFD calculation of the gas line which returns a cycle averaged quantity of heat to the solid parts. Using an iteration method which is based on the temperature deviation between two iteration steps, the total number of iterations and consequently the computational time can be reduced. The loop is continued until a steady-state temperature field is obtained. Calculated temperatures of the solid parts are verified by temperature measurements of a calorimeter test bench.展开更多
The numerical method for computing the live load distribution coefficients in bridge decks is presented. The grillage analogy for representation of bridge decks is adopted in determining the general behavior under tra...The numerical method for computing the live load distribution coefficients in bridge decks is presented. The grillage analogy for representation of bridge decks is adopted in determining the general behavior under traffic loads. The principles of Maxwell's reciprocal theorem are developed in computing live load distribution coefficients and their influence lines. The presented method uses the approach developed in traditional methods of transversal live load distribution but bridge decks are modeled more realistic with the help of well-established grillage analogy. Simple numerical programs for grillage analysis can be used and no special software is needed. While computing the distribution coefficients for a bridge deck the rest of the analysis can be performed with habitual procedures of structural mechanics.展开更多
Recently, microreactors have become available to be fabricated and used safely. The performance of these microreactors depends on the behavior of the multiphase flow hydrodynamics. Gas-liquid flow through T-junction m...Recently, microreactors have become available to be fabricated and used safely. The performance of these microreactors depends on the behavior of the multiphase flow hydrodynamics. Gas-liquid flow through T-junction microchannel reactor is simulated numerically using VOF (volume of fluid) method. 2-D (Two-dimensional) and 3-D (three-dimensional) models of the T-junction microchannel reactor were introduced to the simulations. Both 2-D and 3-D simulations for nitrogen-water flow were performed in the FLUENT (Fluent. Inc.) computational fluid dynamics package. The third direction effect has been studied by comparing the results of the 2-D and 3-D simulations with the published experimental data. Also, the bubble slug length was calculated for the 2-D and 3-D simulations. Furthermore, the hydrodynamics of the flow was studied for the 2-D and 3-D simulations, and compared with other experimental data. The pressure drop, mean bubble velocity, the velocity distribution and the void fraction were calculated and found to be in good agreement with published data.展开更多
The main aim of this paper is determine the centrifugal fan with forward curved blades aerodynamic characteristics based on numerical modeling. Three variants of geometry were investigated. The first, basic "A&qu...The main aim of this paper is determine the centrifugal fan with forward curved blades aerodynamic characteristics based on numerical modeling. Three variants of geometry were investigated. The first, basic "A" variant contains 12 blades. The geometry of second "B" variant contains 12 blades and 12 semi-blades with optimal length[1]. The third, control variant "C" contains 24 blades without semi-blades. Numerical calculations were performed by CFD Ansys. Another aim of this paper is to compare results of the numerical simulation with results of approximate numerical procedure. Applied approximate numerical procedure [2] is designated to determine characteristics of the turbulent flow in the bladed space of a centrifugal-flow fan impeller. This numerical method is an extension of the hydro-dynamical cascade theory for incompressible and inviscid fluid flow. Paper also partially compares results from the numerical simulation and results from the experimental investigation. Acoustic phenomena observed during experiment, during numerical simulation manifested as deterioration of the calculation stability, residuals oscillation and thus also as a flow field oscillation. Pressure pulsations are evaluated by using frequency analysis for each variant and working condition.展开更多
In order to expand the study on flow instability of supercritical circulating fluidized bed(CFB) boiler,a new numerical computational model considering the heat storage of the tube wall metal was presented in this pap...In order to expand the study on flow instability of supercritical circulating fluidized bed(CFB) boiler,a new numerical computational model considering the heat storage of the tube wall metal was presented in this paper.The lumped parameter method was proposed for wall temperature calculation and the single channel model was adopted for the analysis of flow instability.Based on the time-domain method,a new numerical computational program suitable for the analysis of flow instability in the water wall of supercritical CFB boiler with annular furnace was established.To verify the code,calculation results were respectively compared with data of commercial software.According to the comparisons,the new code was proved to be reasonable and accurate for practical engineering application in analysis of flow instability.Based on the new program,the flow instability of supercritical CFB boiler with annular furnace was simulated by time-domain method.When 1.2 times heat load disturbance was applied on the loop,results showed that the inlet flow rate,outlet flow rate and wall temperature fluctuated with time eventually remained at constant values,suggesting that the hydrodynamic flow was stable.The results also showed that in the case of considering the heat storage,the flow in the water wall is easier to return to stable state than without considering heat storage.展开更多
With the rapid development of the computational fluid dynamics(CFD),a parameter-free upwind scheme capable of simulating all speeds accurately and efficiently is in high demand.To achieve this goal,we present a new up...With the rapid development of the computational fluid dynamics(CFD),a parameter-free upwind scheme capable of simulating all speeds accurately and efficiently is in high demand.To achieve this goal,we present a new upwind scheme called AUSMPWM in this paper.This scheme computes the numerical mass flux as the AUSMPW+and computes the interfacial sound speed in a different way.Also,it computes the pressure flux by limiting the dissipation if the Mach number is less than 1.Series of numerical experiments show that AUSMPWM can satisfy the following attractive properties independent of any tuning coefficient:(1)Robustness against the shock anomaly and high discontinuity’s resolution;(2)high accuracy on hypersonic heating prediction and capability to give smooth reproductions of heating profiles;(3)low dissipation at low speeds;and(4)strong grid,reconstruction scheme,and Mach number independence in low speeds’simulations.These properties suggest that AUSMPWM is promising to be widely used to accurately and efficiently simulate flows of all speeds.展开更多
Factorization of polynomials is one of the foundations of symbolic computation.Its applications arise in numerous branches of mathematics and other sciences.However,the present advanced programming languages such as C...Factorization of polynomials is one of the foundations of symbolic computation.Its applications arise in numerous branches of mathematics and other sciences.However,the present advanced programming languages such as C++ and J++,do not support symbolic computation directly.Hence,it leads to difficulties in applying factorization in engineering fields.In this paper,the authors present an algorithm which use numerical method to obtain exact factors of a bivariate polynomial with rational coefficients.The proposed method can be directly implemented in efficient programming language such C++ together with the GNU Multiple-Precision Library.In addition,the numerical computation part often only requires double precision and is easily parallelizable.展开更多
This paper analyzes the correctly-expanded supersonic jet from a convergent-divergent axisymmetric nozzle by using numerical simulation of turbulent flow.And the calculated density distributions in this flow are compa...This paper analyzes the correctly-expanded supersonic jet from a convergent-divergent axisymmetric nozzle by using numerical simulation of turbulent flow.And the calculated density distributions in this flow are compared with the present experimental data using rainbow schlieren deflectometry.The value of the density from the experimental data agrees well with the results calculated by this simulation.Therefore,the present method of the measurement using rainbow schlieren deflectometry is useful for the measurement of the density of the correctly-expanded supersonic jet.展开更多
The Hellmann potential, which is a superposition of an attractive Coulomb potential -air and a Yutmwa potential b e-δr /r , is often used to compute bound-state normalizations and energy levels of neutral atoms. By u...The Hellmann potential, which is a superposition of an attractive Coulomb potential -air and a Yutmwa potential b e-δr /r , is often used to compute bound-state normalizations and energy levels of neutral atoms. By using the generalized parametric Nikiforov-Uvarov (NU) method, we have obtained the approximate analytical solutions of the radial Schr6dinger equation (SE) for the Hellmann potential. The energy eigenvalues and corresponding eigenfunctions are calculated in closed forms. Some numerical results are presented, which show good agreement with a numerical amplitude phase method and also those previously obtained by other methods. As a particular case, we find the energy levels of the pure Coulomb potential.展开更多
基金Projects 10172031 supported by National Natural Science Foundation of China and HIT. MD2002.15 by the Multidiscipline Scientific Research Founda-tion of Harbin Institute of Technology
文摘The sealing performance of magnetic fluid is related to the magnetic fluid itself. Many factors can influence the magnetic field and the seal pressure differences of magnetic fluid seals, such as the sealing gap, the shaft eccentricity, the shaft diameter, the volume of the magnetic fluid and the centrifugal force. These factors are analyzed by numerical computation. When the seal material and structure are the same, the seal pressure difference is directly proportional to the magnetic field intensity and the saturation magnetization of the magnetic fluid. The sealing performance of the magnetic fluid will reduce with the increase of the sealing gap and shaft eccentricity. The sealing performance will increase with the volume of the magnetic fluid and decrease with the increase of the shaft diameter taking gravity into account. The increase of the shaft diameter is the same as the reduction of the volume of the magnetic fluid. The magnetic fluid cross-section can change because of the centrifugal force. Some improvements can reducc the influence of the centrifugal force. The centrifugal force can be utilized to improve the sealing performance.
文摘In this paper, two novel numerical computation methods are introduced which have been recently developed at Research Institute for Applied Mechanics ( R/AM ), Kyushu University, for strongly nonlinear wave-body interaction problems, such as ship motions in rough seas and resulting green-water impact on deck. The first method is the CIP-based Cartesian grid method, in which the free surface flow is treated as a multi-phase flow which is solved using a Cartesian grid. The second method is the MPS method, which is a so-called particle method and hence no grid is used. The features and calculation procedures of these numerical methods are described. One validation computation against a newly conducted experiment on a dam break problem, which is also described in this paper, is presented.
文摘The rubber circular plate is considered as a kind of membrane. Based on the character that there exists no bending moment inside a membrane, the geometric behavior of the rubber circular plate in expanding state was described with the aid of a group of mathematic method. The relationship between deflection and load was attained by means of calculating stress and strain inside the curved surface of rubber plate. Meantime, based on Hencky method, the relationship between deflection and load was attained and considered as the Hencky solution. The different results given rise by the two different resolving methods were compared. The deviation results from the Hencky method was discussed, and a kind of correcting method was put forward.
基金Supported by the National Natural Science Foundation of China (50099620, 40804027)
文摘To study the grouting reinforcement mechanism in jointed rock slope, first, the theoretical deduction was done to calculate the critical length of slipping if the slope angle is larger than that of joint inclination; Second, the numerical calculation model was founded by FLAG^3D, so as to find the stress and deformation responses of rock mass in the state before and after grouting, the analysis results show that the range between the boundary of critical slipping block and the joint plane that passes the slope toe is the effective grouting area (EGA). After excavation, large deformation occurs along the joint plane. After grouting, the displacements of rock particles become uniform and continuous, and large deformations along the joint plane are controlled; the dynamic displacement can re- flect the deformation response of slope during excavation in the state before and after grouting, as well as the shear location of potential slip plane. After grouting, the dynamic displacement of each monitoring point reaches the peak value with very few time steps, which indicate that the parameters of the joint plane, such as strength and stiffness, are improved; the stress field becomes uniform. Tensile area reduces gradually; whole stability of the slope and its ability to resist tensile and shear stress are improved greatly.
基金support of JASSO to conduct this research work during the author’s stay at Japan
文摘Nowadays,there are many studies conducted in the field of marine hydrodynamics which focus on two vessels traveling and floating in sufficiently close proximity to experience significant interactions.The hydrodynamic behavior of parallel moving ships in waves is an interesting and important topic of late.A numerical investigation has been carried out for the prediction of wave exciting forces and motion responses of parallel moving ships in regular waves.The numerical solution was based on 3D distribution technique and using the linear wave theory to determine the exciting forces and ship's motion.The speed effects have been considered in the Green function for more realistic results.The numerical computations of wave exciting forces and motion responses were carried out for a Mariner and Series 60 for the purpose of discovering different Froude numbers and different separation distances in head sea conditions.Based on the numerical computations,it was revealed that the sway,roll and yaw have a significant effect due to hydrodynamic interaction.
基金Projects(51108224,51179026) supported by the National Natural Science Foundation of China
文摘As a new type of structure which has never been built, submerged floating tunnel was studied mainly by numerical simulations. To further study the seismic response of a submerged floating tunnel, the first model experiment of submerged floating tunnel (SFT) under the earthquake was carried out on the unique underwater shaking table in China. The experimental results show that vertical excitation induces larger response than horizontal and different inclination degrees of the tether also cause different seismic responses. Subsequently, based on the fluid-structure interaction theory, the corresponding numerical model is established. And comparing the numerical results with the experimental results, those of shaking table test. Numerical model adopted is effective for it is shown that the numerical results are basically identical with dynamic response of SFT.
基金Supported by Harbin Talents of Science and Technology Innovation Special Fund(2011RFQXG021)
文摘Experiments involving a sonar platform with a sound absorption wedge were carried out for the purpose of obtaining the low frequency acoustic characteristics. Acoustic characteristics of a sonar platform model with a sound absorption wedge were measured, and the effects of different wedge laid areas on platform acoustic characteristic were tested. Vibration acceleration and self-noise caused by model vibration were measured in four conditions: 0%, 36%, 60%, and 100% of wedge laid area when the sonar platform was under a single frequency excitation force. An experiment was performed to validate a corresponding numerical calculation. The numerical vibration characteristics of platform area were calculated by the finite element method, and self-noise caused by the vibration in it was predicted by an experiential formula. The conclusions prove that the numerical calculation method can partially replace the experimental process for obtaining vibration and sound characteristics.
基金provided by the Research Fund of the Fundamental Research Funds for the Central Universities of China University of Mining & Technology (No. 2010ZDP02B02)the State Key Laboratory of Coal Resources and Mine Safety (No.SKLCRSM08X2)+2 种基金the Jiangsu "333"High Qualified Talentsthe National Natural Science Foundation of China (Nos. 50904063 and51004101)the Scientific Research Foundation of China University of Mining & Technology (Nos. 2008A003 and 2009A001)
文摘We determi:ned a suitable gate road layout in slice mining in an ultra-thick unstable coal seam, using theoretical anallysis and numerical calculations. Based on plasticity theory in terms of limiting equilibrium, the width of chain pillar in the upper slice was calculated to be 18 m. The stress distribution in the chain pillar after the upper slice was mined out was described with numerical simulation. The extent of the effect of stress on the upper chain pillar on the lower solid coal was obtained on the basis of an elastic solution of a distributed force loaded on a half-plane. Three layout designs for lower gate roads were pro- posed and a stability factor was introduced to analyze the stability of the lower pillar with numerical calculation. Gate road translation was determined as the most suitable layout method, which maximizes the extraction rate on the basis of the pillar stability.
文摘A numerical model to predict the temperature field in a hermetic reciprocating compressor for household refrigeration appliances is presented in this work. The model combines a high resolution three-dimensional heat conduction formulation of the compressor's solid parts, a three-dimensional CFD (computational fluid dynamics) approach for the gas line domain and lumped formulations of the shell gas and the lubrication oil. Heat transfer coefficients are determined by applying CFD to the gas line side and correlations from the literature on the shell gas and oil side, respectively. The valve in the gas line simulation is modelled as a parallel moving fiat plate. By means of an iterative loop the temperature field of the solid parts acts as boundary condition for the CFD calculation of the gas line which returns a cycle averaged quantity of heat to the solid parts. Using an iteration method which is based on the temperature deviation between two iteration steps, the total number of iterations and consequently the computational time can be reduced. The loop is continued until a steady-state temperature field is obtained. Calculated temperatures of the solid parts are verified by temperature measurements of a calorimeter test bench.
文摘The numerical method for computing the live load distribution coefficients in bridge decks is presented. The grillage analogy for representation of bridge decks is adopted in determining the general behavior under traffic loads. The principles of Maxwell's reciprocal theorem are developed in computing live load distribution coefficients and their influence lines. The presented method uses the approach developed in traditional methods of transversal live load distribution but bridge decks are modeled more realistic with the help of well-established grillage analogy. Simple numerical programs for grillage analysis can be used and no special software is needed. While computing the distribution coefficients for a bridge deck the rest of the analysis can be performed with habitual procedures of structural mechanics.
文摘Recently, microreactors have become available to be fabricated and used safely. The performance of these microreactors depends on the behavior of the multiphase flow hydrodynamics. Gas-liquid flow through T-junction microchannel reactor is simulated numerically using VOF (volume of fluid) method. 2-D (Two-dimensional) and 3-D (three-dimensional) models of the T-junction microchannel reactor were introduced to the simulations. Both 2-D and 3-D simulations for nitrogen-water flow were performed in the FLUENT (Fluent. Inc.) computational fluid dynamics package. The third direction effect has been studied by comparing the results of the 2-D and 3-D simulations with the published experimental data. Also, the bubble slug length was calculated for the 2-D and 3-D simulations. Furthermore, the hydrodynamics of the flow was studied for the 2-D and 3-D simulations, and compared with other experimental data. The pressure drop, mean bubble velocity, the velocity distribution and the void fraction were calculated and found to be in good agreement with published data.
基金supported by Technology Agency of the Czech Republic TAR CK TEO 01020036
文摘The main aim of this paper is determine the centrifugal fan with forward curved blades aerodynamic characteristics based on numerical modeling. Three variants of geometry were investigated. The first, basic "A" variant contains 12 blades. The geometry of second "B" variant contains 12 blades and 12 semi-blades with optimal length[1]. The third, control variant "C" contains 24 blades without semi-blades. Numerical calculations were performed by CFD Ansys. Another aim of this paper is to compare results of the numerical simulation with results of approximate numerical procedure. Applied approximate numerical procedure [2] is designated to determine characteristics of the turbulent flow in the bladed space of a centrifugal-flow fan impeller. This numerical method is an extension of the hydro-dynamical cascade theory for incompressible and inviscid fluid flow. Paper also partially compares results from the numerical simulation and results from the experimental investigation. Acoustic phenomena observed during experiment, during numerical simulation manifested as deterioration of the calculation stability, residuals oscillation and thus also as a flow field oscillation. Pressure pulsations are evaluated by using frequency analysis for each variant and working condition.
基金supported by the "Strategic Priority Research Program" of the Chinese Academy of Sciences,Grant No.XDA07030100the National Key Technology R&D Program of China during the 12th Five-Year Plan Period No.2015BAA03B01-01
文摘In order to expand the study on flow instability of supercritical circulating fluidized bed(CFB) boiler,a new numerical computational model considering the heat storage of the tube wall metal was presented in this paper.The lumped parameter method was proposed for wall temperature calculation and the single channel model was adopted for the analysis of flow instability.Based on the time-domain method,a new numerical computational program suitable for the analysis of flow instability in the water wall of supercritical CFB boiler with annular furnace was established.To verify the code,calculation results were respectively compared with data of commercial software.According to the comparisons,the new code was proved to be reasonable and accurate for practical engineering application in analysis of flow instability.Based on the new program,the flow instability of supercritical CFB boiler with annular furnace was simulated by time-domain method.When 1.2 times heat load disturbance was applied on the loop,results showed that the inlet flow rate,outlet flow rate and wall temperature fluctuated with time eventually remained at constant values,suggesting that the hydrodynamic flow was stable.The results also showed that in the case of considering the heat storage,the flow in the water wall is easier to return to stable state than without considering heat storage.
基金supported by the National Basic Research Program of China("973"Project)(Grant No.2009CB724104)
文摘With the rapid development of the computational fluid dynamics(CFD),a parameter-free upwind scheme capable of simulating all speeds accurately and efficiently is in high demand.To achieve this goal,we present a new upwind scheme called AUSMPWM in this paper.This scheme computes the numerical mass flux as the AUSMPW+and computes the interfacial sound speed in a different way.Also,it computes the pressure flux by limiting the dissipation if the Mach number is less than 1.Series of numerical experiments show that AUSMPWM can satisfy the following attractive properties independent of any tuning coefficient:(1)Robustness against the shock anomaly and high discontinuity’s resolution;(2)high accuracy on hypersonic heating prediction and capability to give smooth reproductions of heating profiles;(3)low dissipation at low speeds;and(4)strong grid,reconstruction scheme,and Mach number independence in low speeds’simulations.These properties suggest that AUSMPWM is promising to be widely used to accurately and efficiently simulate flows of all speeds.
基金partly supported by the National Natural Science Foundation of China under Grant Nos.91118001 and 11170153the National Key Basic Research Project of China under Grant No.2011CB302400Chongqing Science and Technology Commission Project under Grant No.cstc2013jjys40001
文摘Factorization of polynomials is one of the foundations of symbolic computation.Its applications arise in numerous branches of mathematics and other sciences.However,the present advanced programming languages such as C++ and J++,do not support symbolic computation directly.Hence,it leads to difficulties in applying factorization in engineering fields.In this paper,the authors present an algorithm which use numerical method to obtain exact factors of a bivariate polynomial with rational coefficients.The proposed method can be directly implemented in efficient programming language such C++ together with the GNU Multiple-Precision Library.In addition,the numerical computation part often only requires double precision and is easily parallelizable.
文摘This paper analyzes the correctly-expanded supersonic jet from a convergent-divergent axisymmetric nozzle by using numerical simulation of turbulent flow.And the calculated density distributions in this flow are compared with the present experimental data using rainbow schlieren deflectometry.The value of the density from the experimental data agrees well with the results calculated by this simulation.Therefore,the present method of the measurement using rainbow schlieren deflectometry is useful for the measurement of the density of the correctly-expanded supersonic jet.
文摘The Hellmann potential, which is a superposition of an attractive Coulomb potential -air and a Yutmwa potential b e-δr /r , is often used to compute bound-state normalizations and energy levels of neutral atoms. By using the generalized parametric Nikiforov-Uvarov (NU) method, we have obtained the approximate analytical solutions of the radial Schr6dinger equation (SE) for the Hellmann potential. The energy eigenvalues and corresponding eigenfunctions are calculated in closed forms. Some numerical results are presented, which show good agreement with a numerical amplitude phase method and also those previously obtained by other methods. As a particular case, we find the energy levels of the pure Coulomb potential.