Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were ex...Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.展开更多
In this paper,to simulate the three dimensional turbulent flow in suddenly expanded rectangular duct numerically,the SIMPLEC algorithm is employed to solve the incompressible Navier-Stckes equation with k-εturbulenc...In this paper,to simulate the three dimensional turbulent flow in suddenly expanded rectangular duct numerically,the SIMPLEC algorithm is employed to solve the incompressible Navier-Stckes equation with k-εturbulence model.The numerical resulis show well the three dimensional turbulent flow field in the rectangular duct behind the sudden expansion cross-section. and agree.fairly well with the experimental result of the length of the main circumfluence.The numerical method of this paper can be applied to numerical analysis of this kind of turbulent flow.展开更多
By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical so...By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical solution of the modified wave action equation and eikonal equation, which can describe the wave refraction and diffraction. The hydrodynamic model is driven by the wave-induced radiation stresses and affected by the wave turbulence. The numerical implementation of the module has used the finite-volume schemes on unstructured grid, which provides great flexibility for modeling the waves and currents in the complex actual nearshore, and ensures the conservation of energy propagation. The applicability of the proposed model is evaluated in calculating the cases of wave set-up, longshore currents, undertow on a sloping beach, rip currents and meandering longshore currents on a tri-cuspate beach. The results indicate that it is necessary to introduce the depth-dependent radiation stresses into the numerical simulation of wave-induced currents, and comparisons show that the present model makes better prediction on the wave procedure as well as both horizontal and vertical structures in the wave-induced current field.展开更多
In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solut...In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.展开更多
The inlcrunnual variation of the vertical distribution of ozone in the tropical stratosphere and its quasi—biennial oscillation (QBO) is analyzed using HALOE data. The results are compared with the wind QBO. A numeri...The inlcrunnual variation of the vertical distribution of ozone in the tropical stratosphere and its quasi—biennial oscillation (QBO) is analyzed using HALOE data. The results are compared with the wind QBO. A numerical experiment is carried out to study the effects of wind QBO on the distribution, and variation of ozone in the stratosphere by using (he NCAR interactive chemical, dynamical, and radiative two—dimensional model (SOCRATES). Data analysis shows that the location of the maximum ozone mixing ratio in the stratosphere changes in the meridional and vertical directions, and assumes a quasi—biennial period. The meridional and vertical motion of the maximum mixing ratio leads to a QBO of column ozone and its hemispheric asymmetry. The QBO of the location of the maximum is closely connected with the zonal wind QBO. The data analysis also shows that in the tropical region, the phase of the QBO for ozone density changes many times with height. Numerical simulation shows that the meridional circulation induced by the wind QBO includes three pairs of cells in the stratosphere, which have hemispheric symmetry. The transport of ozone by the induced meridional circulation in various latitudes and heights is the main dynamic cause for the ozone QBO. Cells of the induced circulation in the middle stratosphere (25-35 km) play an important role in producing the ozone QBO.展开更多
Based on the 2D horizontal plane numerical model, a quasi-3D numerical model is established for coastal regions of shallow water. The characteristics of this model are that the velocity profiles;can be obtained at the...Based on the 2D horizontal plane numerical model, a quasi-3D numerical model is established for coastal regions of shallow water. The characteristics of this model are that the velocity profiles;can be obtained at the same time when the equations of the value of difference between the horizontal current velocity and its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent with the results of the 2D, model. The circulating flow in the rectangular area induced by wind is simulated and applied to the tidal flow field of the radial sandbanks in the South Yellow Sea. The computational results from this quasi-3D model are in good agreement with analytical results and observed data. The solution of the finite difference equations has been found to be stable, and the model is simple, effective and practical.展开更多
The dimension variation of B610E steel during the quenching and tempering process was simulated by ABAQUS,the coefficient of heat transfer during quenching was verified by the buried thermocouple test and the trend of...The dimension variation of B610E steel during the quenching and tempering process was simulated by ABAQUS,the coefficient of heat transfer during quenching was verified by the buried thermocouple test and the trend of the dimension variation during the quenching process was also calculated by ABAQUS. It was shown by the comparison of simulated results and industrially measured results that only the simulation of thickness and width variation was accurate, while length variation needed further simulation. Besides, the dimension variation trend was identical with the measured results.展开更多
According to Cubic law and incompressible fluid law of mass conservation, the seepage character of the fracture surface was simulated with the simulation method of fractal theory and random Brown function. Furthermore...According to Cubic law and incompressible fluid law of mass conservation, the seepage character of the fracture surface was simulated with the simulation method of fractal theory and random Brown function. Furthermore, the permeability coefficient of the single fracture was obtained. In order to test the stability of the method, 500 simulations were conducted on each different fractal dimension. The simulated permeability coefficient was analyzed in probability density distribution and probability cumulative distribution statistics. Statistics showed that the discrete degree of the permeability coefficient increases with the increase of the fractal dimension. And the calculation result has better stability when the fractal dimension value is relatively small. According to the Bayes theory, the characteristic index of the permeability coefficient on fractal dimension P（Dfi｜ Ri） is established. The index, P（Dfi｜ Ri）, shows that when the simulated permeability coefficient is relatively large, it can clearly represent the fractal dimension of the structure surface, the probability is 82%. The calculated results of the characteristic index verify the feasibility of the method.展开更多
Two-dimensional unsteady incompressible viscous flow around a rolling cylinder with ship-like section is numerically simulated by employing the computational scheme previously developed by the authors, in which the co...Two-dimensional unsteady incompressible viscous flow around a rolling cylinder with ship-like section is numerically simulated by employing the computational scheme previously developed by the authors, in which the continuity and momentum equations are satisfied simultaneously at each time step for oscillating flow. The numerical results show that the motion of vortices around a rolling ship hull is cyclical. It is found that the location of the vortices is very similar to the existing experimental result. Using these simulation results, we can calculate the roll damping of ships including viscous effects.展开更多
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.展开更多
The dimensional accuracy of machined parts is strongly influenced by the thermal behavior of machine tools (MT). Minimizing this influence represents a key objective for any modern manufacturing industry. Thermally in...The dimensional accuracy of machined parts is strongly influenced by the thermal behavior of machine tools (MT). Minimizing this influence represents a key objective for any modern manufacturing industry. Thermally induced positioning error compensation remains the most effective and practical method in this context. However, the efficiency of the compensation process depends on the quality of the model used to predict the thermal errors. The model should consistently reflect the relationships between temperature distribution in the MT structure and thermally induced positioning errors. A judicious choice of the number and location of temperature sensitive points to represent heat distribution is a key factor for robust thermal error modeling. Therefore, in this paper, the temperature sensitive points are selected following a structured thermomechanical analysis carried out to evaluate the effects of various temperature gradients on MT structure deformation intensity. The MT thermal behavior is first modeled using finite element method and validated by various experimentally measured temperature fields using temperature sensors and thermal imaging. MT Thermal behavior validation shows a maximum error of less than 10% when comparing the numerical estimations with the experimental results even under changing operation conditions. The numerical model is used through several series of simulations carried out using varied working condition to explore possible relationships between temperature distribution and thermal deformation characteristics to select the most appropriate temperature sensitive points that will be considered for building an empirical prediction model for thermal errors as function of MT thermal state. Validation tests achieved using an artificial neural network based simplified model confirmed the efficiency of the proposed temperature sensitive points allowing the prediction of the thermally induced errors with an accuracy greater than 90%.展开更多
The principle of surface wave plasma discharge in a rectangular cavity is introduced simply based on surface plasmon polariton theory. The distribution of surface-wave electric field at the interface of the plasma-die...The principle of surface wave plasma discharge in a rectangular cavity is introduced simply based on surface plasmon polariton theory. The distribution of surface-wave electric field at the interface of the plasma-dielectric slab is investigated by using the three-dimensional finite-difference time-domain method （3D-FDTD） with different slotantenna structures. And the experimental image of discharge with a novel slot antenna array and the simulation of the electric field with this slot antenna array are both displayed. Combined with the distribution of surface wave excitation and experimental results, the numerical simulation performed by using 3D-FDTD is shown to be a useful tool in the computer-aided antenna design for large area planar-type surface-wave plasma sources.展开更多
Numerical simulation of oil migration and accumulation is to describe the history of oil migration and accumulation in basin evolution. It is of great value in the exploration of oil resources and their rational evalu...Numerical simulation of oil migration and accumulation is to describe the history of oil migration and accumulation in basin evolution. It is of great value in the exploration of oil resources and their rational evaluation. In this paper, from such actual conditions as the effects of mechanics of fluids in porous media and 3-dimensional geology characteristics, a kind of modified method of second order splitting-up implicit interactive scheme is pur forward. For the famous hydraulic experiment of secondary migration-accumulation, the numerical simulation test has been done, and both the computational and experimental results are basically identical. For the actual problem of Dongying hollow of Shengli Petroleum Oil Field, the numerical simulation test and the actual conditions are basically coincident. Thus the well-known problem has been solved.展开更多
The characteristics of three-dimensional （3-D） tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, ma...The characteristics of three-dimensional （3-D） tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, σ 3-D numerical model oil an unstructured triangular grid has been degeloped. The σ coordinate transforination, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simulations in the Oujing Estuary with submerged dyke and moveable boundary problems. The model is evaluated with the in situ data, and the results show that the calculated water elevations at 19 stations and currents at 19 profiler stations are in good agreement with measured data both in magnitude and phase. This numerical model is applied to the 3-D tidal circulation simulations of experiments in stopping flow transport through the South Branch of the Oujiang Estuary, and the feasibility to cutoff the flow in the South Branch of the Oujiang Estuary is demonstrated by numerical simulation experiments. The developed numerical model simulated the 3-D tidal current circulations in complicated coastal and estuarine waters very well.展开更多
Based on the low inductance technology and parallel-plate transmission principle,an experimental apparatus of small-scale slapper initiating primary high explosives driven by electrical explosion is designed and estab...Based on the low inductance technology and parallel-plate transmission principle,an experimental apparatus of small-scale slapper initiating primary high explosives driven by electrical explosion is designed and established.The problem of instantaneously and continuously measuring the velocity of the small-scale slapper is successfully solved by using the technique of laser interference.Compared with the results published on the literatures at home and abroad,data of the experimental and the numerical simulation shown in this paper are more proper to reflect the physical process of electrical explosion driving slapper.One-dimensional numerical simulation of electrical explosion driving slapper is done using the hydrodynamic code.The experimental results are consistent with the computed ones by introducing a power correction factor.In the end,the introduced power correction factor is discussed.展开更多
The analytical solution of the convection diffusion equation is considered by two-dimensional Fourier transform and the inverse Fourier transform. To get the numerical solution, the Crank-Nicolson finite difference me...The analytical solution of the convection diffusion equation is considered by two-dimensional Fourier transform and the inverse Fourier transform. To get the numerical solution, the Crank-Nicolson finite difference method is constructed, which is second-order accurate in time and space. Numerical simulation shows excellent agreement with the analytical solution. The dynamic visualization of the simulating results is realized on ArcGIS platform. This work provides a quick and intuitive decision-making basis for water resources protection, especially in dealing with water pollution emergencies.展开更多
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.展开更多
The characteristics of collisional radio-frequency (rf) sheath dynamics over an elec-trode with a cylindrical hole is simulated by means of a self-consistent model which consists of two-dimensional time-dependent flui...The characteristics of collisional radio-frequency (rf) sheath dynamics over an elec-trode with a cylindrical hole is simulated by means of a self-consistent model which consists of two-dimensional time-dependent fluid equations coupled with Poisson equation. In addition, an equivalent-circuit model is coupled to the fluid equations in order to self-consistently determine re-lationship between the instantaneous potential at the rf-biased electrode and the sheath thickness. Two-dimensional profiles of the potential, the ion fluid velocity, and the distributions of the ion and electron densities within the sheath are computed under various discharge conditions, such as the discharge powers and the gas pressures. The results show that the existence of the cylindrical hole on the electrode significantly affects the sheath structure and generates a potential trap in the horizontal direction, which is particularly strong when the sheath thickness is comparable to the depth of the hole. Moreover, it is found that the collisional effects have a significant influence on the sheath characteristics.展开更多
The characteristics of a collisional dual frequency （DF） sheath near an electrode with a cylindrical hole are studied by utilizing a two-dimensional model which includes time-dependent fluid equations coupled with t...The characteristics of a collisional dual frequency （DF） sheath near an electrode with a cylindrical hole are studied by utilizing a two-dimensional model which includes time-dependent fluid equations coupled with the Poisson equation and an equivalent-circuit model, The effects of the gas pressure on the two-dimensional profiles of the potential, electric field, ion fluid velocity in a DF sheath are investigated. The simulation results show that the cylindrical hole on the electrode has a significant influence on the DF sheath structure, i.e., the sheath profile tends to wrap around the contour of the hole feature. Moreover, it is shown that the structure of the DF sheath is different from that of a single frequency （SF） sheath because the profile of the DF sheath is modulated by the combination of the high and low frequency sources. In addition the characteristics of the DF sheath are obviously affected by the collisional effects in the DF sheath.展开更多
基金Supported by the China National Science and Technology Major Project(2016ZX05033002,2016ZX05033001).
文摘Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.
文摘In this paper,to simulate the three dimensional turbulent flow in suddenly expanded rectangular duct numerically,the SIMPLEC algorithm is employed to solve the incompressible Navier-Stckes equation with k-εturbulence model.The numerical resulis show well the three dimensional turbulent flow field in the rectangular duct behind the sudden expansion cross-section. and agree.fairly well with the experimental result of the length of the main circumfluence.The numerical method of this paper can be applied to numerical analysis of this kind of turbulent flow.
基金financially supported by the the National Natural Science Foundation of China(Grant No.51709054)the Public Science and Technology Research Funds Projects of Ocean(Grant Nos.201405025 and 201505019)
文摘By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical solution of the modified wave action equation and eikonal equation, which can describe the wave refraction and diffraction. The hydrodynamic model is driven by the wave-induced radiation stresses and affected by the wave turbulence. The numerical implementation of the module has used the finite-volume schemes on unstructured grid, which provides great flexibility for modeling the waves and currents in the complex actual nearshore, and ensures the conservation of energy propagation. The applicability of the proposed model is evaluated in calculating the cases of wave set-up, longshore currents, undertow on a sloping beach, rip currents and meandering longshore currents on a tri-cuspate beach. The results indicate that it is necessary to introduce the depth-dependent radiation stresses into the numerical simulation of wave-induced currents, and comparisons show that the present model makes better prediction on the wave procedure as well as both horizontal and vertical structures in the wave-induced current field.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10732010,10972010,and 11332002)
文摘In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.
文摘The inlcrunnual variation of the vertical distribution of ozone in the tropical stratosphere and its quasi—biennial oscillation (QBO) is analyzed using HALOE data. The results are compared with the wind QBO. A numerical experiment is carried out to study the effects of wind QBO on the distribution, and variation of ozone in the stratosphere by using (he NCAR interactive chemical, dynamical, and radiative two—dimensional model (SOCRATES). Data analysis shows that the location of the maximum ozone mixing ratio in the stratosphere changes in the meridional and vertical directions, and assumes a quasi—biennial period. The meridional and vertical motion of the maximum mixing ratio leads to a QBO of column ozone and its hemispheric asymmetry. The QBO of the location of the maximum is closely connected with the zonal wind QBO. The data analysis also shows that in the tropical region, the phase of the QBO for ozone density changes many times with height. Numerical simulation shows that the meridional circulation induced by the wind QBO includes three pairs of cells in the stratosphere, which have hemispheric symmetry. The transport of ozone by the induced meridional circulation in various latitudes and heights is the main dynamic cause for the ozone QBO. Cells of the induced circulation in the middle stratosphere (25-35 km) play an important role in producing the ozone QBO.
基金National Natural Science Foundation of China(Grant No.49236120)
文摘Based on the 2D horizontal plane numerical model, a quasi-3D numerical model is established for coastal regions of shallow water. The characteristics of this model are that the velocity profiles;can be obtained at the same time when the equations of the value of difference between the horizontal current velocity and its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent with the results of the 2D, model. The circulating flow in the rectangular area induced by wind is simulated and applied to the tidal flow field of the radial sandbanks in the South Yellow Sea. The computational results from this quasi-3D model are in good agreement with analytical results and observed data. The solution of the finite difference equations has been found to be stable, and the model is simple, effective and practical.
文摘The dimension variation of B610E steel during the quenching and tempering process was simulated by ABAQUS,the coefficient of heat transfer during quenching was verified by the buried thermocouple test and the trend of the dimension variation during the quenching process was also calculated by ABAQUS. It was shown by the comparison of simulated results and industrially measured results that only the simulation of thickness and width variation was accurate, while length variation needed further simulation. Besides, the dimension variation trend was identical with the measured results.
基金Project(50934006) supported by the National Natural Science Foundation of ChinaProject(CX2012B070) supported by Hunan Provincial Innovation Foundation for Postgraduate,ChinaProject(1343-76140000024) Supported by Academic New Artist Ministry of Education Doctoral Post Graduate in 2012,China
文摘According to Cubic law and incompressible fluid law of mass conservation, the seepage character of the fracture surface was simulated with the simulation method of fractal theory and random Brown function. Furthermore, the permeability coefficient of the single fracture was obtained. In order to test the stability of the method, 500 simulations were conducted on each different fractal dimension. The simulated permeability coefficient was analyzed in probability density distribution and probability cumulative distribution statistics. Statistics showed that the discrete degree of the permeability coefficient increases with the increase of the fractal dimension. And the calculation result has better stability when the fractal dimension value is relatively small. According to the Bayes theory, the characteristic index of the permeability coefficient on fractal dimension P（Dfi｜ Ri） is established. The index, P（Dfi｜ Ri）, shows that when the simulated permeability coefficient is relatively large, it can clearly represent the fractal dimension of the structure surface, the probability is 82%. The calculated results of the characteristic index verify the feasibility of the method.
基金This project is supported by Fok Ying Tung Education Foundation and National Natural Science Foundation of China
文摘Two-dimensional unsteady incompressible viscous flow around a rolling cylinder with ship-like section is numerically simulated by employing the computational scheme previously developed by the authors, in which the continuity and momentum equations are satisfied simultaneously at each time step for oscillating flow. The numerical results show that the motion of vortices around a rolling ship hull is cyclical. It is found that the location of the vortices is very similar to the existing experimental result. Using these simulation results, we can calculate the roll damping of ships including viscous effects.
基金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.
文摘The dimensional accuracy of machined parts is strongly influenced by the thermal behavior of machine tools (MT). Minimizing this influence represents a key objective for any modern manufacturing industry. Thermally induced positioning error compensation remains the most effective and practical method in this context. However, the efficiency of the compensation process depends on the quality of the model used to predict the thermal errors. The model should consistently reflect the relationships between temperature distribution in the MT structure and thermally induced positioning errors. A judicious choice of the number and location of temperature sensitive points to represent heat distribution is a key factor for robust thermal error modeling. Therefore, in this paper, the temperature sensitive points are selected following a structured thermomechanical analysis carried out to evaluate the effects of various temperature gradients on MT structure deformation intensity. The MT thermal behavior is first modeled using finite element method and validated by various experimentally measured temperature fields using temperature sensors and thermal imaging. MT Thermal behavior validation shows a maximum error of less than 10% when comparing the numerical estimations with the experimental results even under changing operation conditions. The numerical model is used through several series of simulations carried out using varied working condition to explore possible relationships between temperature distribution and thermal deformation characteristics to select the most appropriate temperature sensitive points that will be considered for building an empirical prediction model for thermal errors as function of MT thermal state. Validation tests achieved using an artificial neural network based simplified model confirmed the efficiency of the proposed temperature sensitive points allowing the prediction of the thermally induced errors with an accuracy greater than 90%.
基金supported by the Foundation for Returned Scholars,the Ministry of Education of China
文摘The principle of surface wave plasma discharge in a rectangular cavity is introduced simply based on surface plasmon polariton theory. The distribution of surface-wave electric field at the interface of the plasma-dielectric slab is investigated by using the three-dimensional finite-difference time-domain method （3D-FDTD） with different slotantenna structures. And the experimental image of discharge with a novel slot antenna array and the simulation of the electric field with this slot antenna array are both displayed. Combined with the distribution of surface wave excitation and experimental results, the numerical simulation performed by using 3D-FDTD is shown to be a useful tool in the computer-aided antenna design for large area planar-type surface-wave plasma sources.
文摘Numerical simulation of oil migration and accumulation is to describe the history of oil migration and accumulation in basin evolution. It is of great value in the exploration of oil resources and their rational evaluation. In this paper, from such actual conditions as the effects of mechanics of fluids in porous media and 3-dimensional geology characteristics, a kind of modified method of second order splitting-up implicit interactive scheme is pur forward. For the famous hydraulic experiment of secondary migration-accumulation, the numerical simulation test has been done, and both the computational and experimental results are basically identical. For the actual problem of Dongying hollow of Shengli Petroleum Oil Field, the numerical simulation test and the actual conditions are basically coincident. Thus the well-known problem has been solved.
基金The Natural Science Foundation of Tianjin, China under contract No.08JCZDZT00200
文摘The characteristics of three-dimensional （3-D） tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, σ 3-D numerical model oil an unstructured triangular grid has been degeloped. The σ coordinate transforination, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simulations in the Oujing Estuary with submerged dyke and moveable boundary problems. The model is evaluated with the in situ data, and the results show that the calculated water elevations at 19 stations and currents at 19 profiler stations are in good agreement with measured data both in magnitude and phase. This numerical model is applied to the 3-D tidal circulation simulations of experiments in stopping flow transport through the South Branch of the Oujiang Estuary, and the feasibility to cutoff the flow in the South Branch of the Oujiang Estuary is demonstrated by numerical simulation experiments. The developed numerical model simulated the 3-D tidal current circulations in complicated coastal and estuarine waters very well.
基金Sponsored by the Foundation of Institute of Fluid Physics of China Academy of Engineering Physics
文摘Based on the low inductance technology and parallel-plate transmission principle,an experimental apparatus of small-scale slapper initiating primary high explosives driven by electrical explosion is designed and established.The problem of instantaneously and continuously measuring the velocity of the small-scale slapper is successfully solved by using the technique of laser interference.Compared with the results published on the literatures at home and abroad,data of the experimental and the numerical simulation shown in this paper are more proper to reflect the physical process of electrical explosion driving slapper.One-dimensional numerical simulation of electrical explosion driving slapper is done using the hydrodynamic code.The experimental results are consistent with the computed ones by introducing a power correction factor.In the end,the introduced power correction factor is discussed.
文摘The analytical solution of the convection diffusion equation is considered by two-dimensional Fourier transform and the inverse Fourier transform. To get the numerical solution, the Crank-Nicolson finite difference method is constructed, which is second-order accurate in time and space. Numerical simulation shows excellent agreement with the analytical solution. The dynamic visualization of the simulating results is realized on ArcGIS platform. This work provides a quick and intuitive decision-making basis for water resources protection, especially in dealing with water pollution emergencies.
文摘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.
基金The project supported by the National Natural Science Foundation of China(No.19975008 and 10275009)and by MOEC (Ministry of Education,China)Grant for Cross-Century Excellent ScholarSupport from the Natural Sciences and Engineering Research Council of C
文摘The characteristics of collisional radio-frequency (rf) sheath dynamics over an elec-trode with a cylindrical hole is simulated by means of a self-consistent model which consists of two-dimensional time-dependent fluid equations coupled with Poisson equation. In addition, an equivalent-circuit model is coupled to the fluid equations in order to self-consistently determine re-lationship between the instantaneous potential at the rf-biased electrode and the sheath thickness. Two-dimensional profiles of the potential, the ion fluid velocity, and the distributions of the ion and electron densities within the sheath are computed under various discharge conditions, such as the discharge powers and the gas pressures. The results show that the existence of the cylindrical hole on the electrode significantly affects the sheath structure and generates a potential trap in the horizontal direction, which is particularly strong when the sheath thickness is comparable to the depth of the hole. Moreover, it is found that the collisional effects have a significant influence on the sheath characteristics.
基金supported by National Natural Science Foundation of China (Nos.10635010,10572035)
文摘The characteristics of a collisional dual frequency （DF） sheath near an electrode with a cylindrical hole are studied by utilizing a two-dimensional model which includes time-dependent fluid equations coupled with the Poisson equation and an equivalent-circuit model, The effects of the gas pressure on the two-dimensional profiles of the potential, electric field, ion fluid velocity in a DF sheath are investigated. The simulation results show that the cylindrical hole on the electrode has a significant influence on the DF sheath structure, i.e., the sheath profile tends to wrap around the contour of the hole feature. Moreover, it is shown that the structure of the DF sheath is different from that of a single frequency （SF） sheath because the profile of the DF sheath is modulated by the combination of the high and low frequency sources. In addition the characteristics of the DF sheath are obviously affected by the collisional effects in the DF sheath.
