Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate ...Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.展开更多
A calculation scheme, which combines a horizontal upwind finite element method with vertical implicit differences, is used to establish a three-dimensional mathematical model of tidal motion and sediment transport in...A calculation scheme, which combines a horizontal upwind finite element method with vertical implicit differences, is used to establish a three-dimensional mathematical model of tidal motion and sediment transport in tidal current. Compared with those of the relative theoretical formula, the results are satisfactory. The model mentioned above has been applied to the water area of the Lianzhou Bay, Guangxi Province. On the basis of the analysis and comparison with the field data, it shows clearly that the model calculation results are reasonable.展开更多
This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method (FEM) by recursive application of the one-dimensional (1D) element energy pro...This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method (FEM) by recursive application of the one-dimensional (1D) element energy projection (EEP) technique. The main idea is to conceptually treat multi-dimensional problems as generalized 1D problems, based on which the concepts of generalized 1D FEM and its consequent EEP formulae have been developed in a unified manner. Equipped with these concepts, multi-dimensional problems can be recursively discretized in one dimension at each step, until a fully discretized standard finite element (FE) model is reached. This conceptual dimension-by- dimension (D-by-D) discretization procedure is entirely equivalent to a full FE discretization. As a reverse D-by-D recovery procedure, by using the unified EEP formulae together with proper extraction of the generalized nodal solutions, super-convergent displacements and first derivatives for two-dimensional (2D) and three-dimensional (3D) problems can be obtained over the domain. Numerical examples of 3D Poisson's equation and elasticity problem are given to verify the feasibility and effectiveness of the proposed strategy.展开更多
In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images ...In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome (OSAHS), three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.展开更多
Exact solutions of three-dimensional(3D)crack problems are much less in number than those of two-dimensional ones,especially for multi-field coupling media exhibiting a certain kind of material anisotropy.An exact3Dth...Exact solutions of three-dimensional(3D)crack problems are much less in number than those of two-dimensional ones,especially for multi-field coupling media exhibiting a certain kind of material anisotropy.An exact3Dthermoelastic solution has been reported for a uniformly heated penny-shaped crack in an infinite magnetoelectric space,with impermeable electromagnetic conditions assumed on the crack faces.Exact 3Dsolutions for the penny-shaped crack subjected to uniform or point temperature load are further presented here when the crack faces are electrically and magnetically permeable.The solutions,obtained by the potential theory method,are exact in the sense that all field variables are explicitly derived and expressed in terms of elementary functions.Along with the previously reported solution,the limits or bounds of the stress intensity factor at the crack-tip for a practical crack can be identified.展开更多
The effects of low-level PES on the overall accuracy of the final surface constructed by using hierarchical construction were investigated with the constructions of a number of global surfaces for the H3 system at UHF...The effects of low-level PES on the overall accuracy of the final surface constructed by using hierarchical construction were investigated with the constructions of a number of global surfaces for the H3 system at UHF (UMP2, DFT-B3LYP, UCCSD(T))/vtz, and UCCSD(T)/avqz levels of theory. The total reaction probabilities for the H+H2 reaction calculated on these surfaces revealed that the accuracy of UCCSD(T)/avqz surface is very close to the well-known BKMP2 surface, while the UCCSD(T)/vtz PES has a slightly higher barrier. In contrast, the low-level theories (UHF, UMP2, DFT-B3LYP) with vtz basis set can only provide a qualitative description of this simplest reaction despite the fact that they are widely used to study reactions in complex systems. On the other hand, although these theories are not accurate on describing the reaction, they can be used to provide the low- level PESs for hierarchical construction of the UCCSD(T)/avqz PES with the number of UCCSD(T)/avqz energies substantially reduced.展开更多
基金supported by the National Natural Science Foundation of China (10872096)the Open Fund of State Key Laboratory of Explosion Science and Technology, Beijing University of Science and Technology (KFJJ09-13)
文摘Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.
文摘A calculation scheme, which combines a horizontal upwind finite element method with vertical implicit differences, is used to establish a three-dimensional mathematical model of tidal motion and sediment transport in tidal current. Compared with those of the relative theoretical formula, the results are satisfactory. The model mentioned above has been applied to the water area of the Lianzhou Bay, Guangxi Province. On the basis of the analysis and comparison with the field data, it shows clearly that the model calculation results are reasonable.
基金supported by the National Natural Science Foundation of China(Nos.51378293 and 51078199)
文摘This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method (FEM) by recursive application of the one-dimensional (1D) element energy projection (EEP) technique. The main idea is to conceptually treat multi-dimensional problems as generalized 1D problems, based on which the concepts of generalized 1D FEM and its consequent EEP formulae have been developed in a unified manner. Equipped with these concepts, multi-dimensional problems can be recursively discretized in one dimension at each step, until a fully discretized standard finite element (FE) model is reached. This conceptual dimension-by- dimension (D-by-D) discretization procedure is entirely equivalent to a full FE discretization. As a reverse D-by-D recovery procedure, by using the unified EEP formulae together with proper extraction of the generalized nodal solutions, super-convergent displacements and first derivatives for two-dimensional (2D) and three-dimensional (3D) problems can be obtained over the domain. Numerical examples of 3D Poisson's equation and elasticity problem are given to verify the feasibility and effectiveness of the proposed strategy.
基金The project supported by the National Natural Science Foundation of China(10672036,10472025 and 10421002)the Natural Science Foundation of Liaoning Province(20032109)
文摘In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome (OSAHS), three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.
基金This work was supported by the National Natural Sci- ence Foundation of China (11321202) and the Specialized Research Fund for the Doctoral Program of Higher Educa- tion (2013010 1110120).
文摘Exact solutions of three-dimensional(3D)crack problems are much less in number than those of two-dimensional ones,especially for multi-field coupling media exhibiting a certain kind of material anisotropy.An exact3Dthermoelastic solution has been reported for a uniformly heated penny-shaped crack in an infinite magnetoelectric space,with impermeable electromagnetic conditions assumed on the crack faces.Exact 3Dsolutions for the penny-shaped crack subjected to uniform or point temperature load are further presented here when the crack faces are electrically and magnetically permeable.The solutions,obtained by the potential theory method,are exact in the sense that all field variables are explicitly derived and expressed in terms of elementary functions.Along with the previously reported solution,the limits or bounds of the stress intensity factor at the crack-tip for a practical crack can be identified.
文摘The effects of low-level PES on the overall accuracy of the final surface constructed by using hierarchical construction were investigated with the constructions of a number of global surfaces for the H3 system at UHF (UMP2, DFT-B3LYP, UCCSD(T))/vtz, and UCCSD(T)/avqz levels of theory. The total reaction probabilities for the H+H2 reaction calculated on these surfaces revealed that the accuracy of UCCSD(T)/avqz surface is very close to the well-known BKMP2 surface, while the UCCSD(T)/vtz PES has a slightly higher barrier. In contrast, the low-level theories (UHF, UMP2, DFT-B3LYP) with vtz basis set can only provide a qualitative description of this simplest reaction despite the fact that they are widely used to study reactions in complex systems. On the other hand, although these theories are not accurate on describing the reaction, they can be used to provide the low- level PESs for hierarchical construction of the UCCSD(T)/avqz PES with the number of UCCSD(T)/avqz energies substantially reduced.