The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three ...The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.展开更多
The present study examines the results of the researches related to the gob bulking factor carried out at home and abroad.A mathematical function of a three-dimensional gob bulking factor is described based on a three...The present study examines the results of the researches related to the gob bulking factor carried out at home and abroad.A mathematical function of a three-dimensional gob bulking factor is described based on a three-dimensional gob model.The method of taking value for interstice and permeability ratios is also proposed.The law of air leakage of fully mechanized top coal is researched in this study.The results show that the speed of air flow near the upper and lower crossheadings is higher than that in the central section of the gob at the same distance from the working face.When the amount of air at the working face exceeds a critical amount,the width of the spontaneous combustion zone in the upper and lower crossheadings is also larger than that in the central section.In this situation,the key is preventing the coal left in the upper and lower crossheadings from self-igniting.Reducing the amount of air at the working face can decrease the width of the spontaneous combustion zone,especially the width near the upper and lower crossheadings.This also moves the spontaneous combustion zone in the direction of the working face.It can prevent the coal in the gob from self-igniting by making the coal left in the crossheadings to be inert and by effectively controlling the amount of air at the working face.展开更多
The three-dimensional(3 D) structures of pores directly affect the CH4 flow.Therefore,it is very important to analyze the3 D spatial structure of pores and to simulate the CH4 flow with the connected pores as the carr...The three-dimensional(3 D) structures of pores directly affect the CH4 flow.Therefore,it is very important to analyze the3 D spatial structure of pores and to simulate the CH4 flow with the connected pores as the carrier.The result shows that the equivalent radius of pores and throats are 1-16 μm and 1.03-8.9 μm,respectively,and the throat length is 3.28-231.25 μm.The coordination number of pores concentrates around three,and the intersection point between the connectivity function and the X-axis is 3-4 μm,which indicate the macro-pores have good connectivity.During the single-channel flow,the pressure decreases along the direction of CH4 flow,and the flow velocity of CH4 decreases from the pore center to the wall.Under the dual-channel and the multi-channel flows,the pressure also decreases along the CH4 flow direction,while the velocity increases.The mean flow pressure gradually decreases with the increase of the distance from the inlet slice.The change of mean flow pressure is relatively stable in the direction horizontal to the bedding plane,while it is relatively large in the direction perpendicular to the bedding plane.The mean flow velocity in the direction horizontal to the bedding plane(Y-axis) is the largest,followed by that in the direction horizontal to the bedding plane(X-axis),and the mean flow velocity in the direction perpendicular to the bedding plane is the smallest.展开更多
A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-co...A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept (EDC) model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.展开更多
Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale...Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.展开更多
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.展开更多
An optimized device structure for reducing the RESET current of phase-change random access memory （PCRAM） with blade-type like （BTL） phase change layer is proposed. The electrical thermal analysis of the BTL cell ...An optimized device structure for reducing the RESET current of phase-change random access memory （PCRAM） with blade-type like （BTL） phase change layer is proposed. The electrical thermal analysis of the BTL cell and the blade heater contactor structure by three-dimensional finite element modeling are compared with each other during RESET operation. The simulation results show that the programming region of the phase change layer in the BTL cell is much smaller, and thermal electrical distributions of the BTL cell are more concentrated on the TiN/GST interface. The results indicate that the BTL cell has the superiorities of increasing the heating efficiency, decreasing the power consumption and reducing the RESET current from 0.67mA to 0.32mA. Therefore, the BTL cell will be appropriate for high performance PCRAM device with lower power consumption and lower RESET current.展开更多
A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fractu...A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.展开更多
A method of constructing three-dimensional process model for the punching cartridge cases is presented based on DEFORM simulation analysis. Using DEFORM software,the finite element simulation models for the punching a...A method of constructing three-dimensional process model for the punching cartridge cases is presented based on DEFORM simulation analysis. Using DEFORM software,the finite element simulation models for the punching and forming process of cartridge cases are established,and the corresponding simulation result model of each intermediate procedure is obtained by continuously performing the forming process simulation. The simulation model cannot annotate size and process information due to poor interface between DEFORM software and CAD software. Thus,a 3D annotation module is developed with secondary development technology of UG NX software. Consequently,the final process model with dimension and process information is obtained. Then,with the current 3D process management system,the 3D punching and forming process design of cartridge cases can be completed further. An example is also provided to illustrate that the relative error between the simulation process model and the physical model is less than 2%,which proves the validity and reliability of the proposed method in this study.展开更多
The temperature field, velocity field, as well as species distribution in three-dimensional space are successfully calculated by establishing three-dimensional geometry model and solving plasma jet-substrate interacti...The temperature field, velocity field, as well as species distribution in three-dimensional space are successfully calculated by establishing three-dimensional geometry model and solving plasma jet-substrate interaction equations, optimized particle trajecory models, as well as particle-particle heat transfer equations in space. Under typical working conditions, the flying trajectories and distribution of ZrO2 ceramic particles and Ni metal particles are also simulated. Results show that, the plasma jet becomes wider near the substrate, and the stochastic trajectory model is preferable to simulate the turbulent diffusion effect of particles. In addition, Ni metal particles penetrate relatively more deeply than ZrO2 ceramic particles due to larger density.展开更多
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.展开更多
Direct numerical simulations（DNSs） of purely elastic turbulence in rectilinear shear flows in a three-dimensional（3D） parallel plate channel were carried out,by which numerical databases were established.Based on ...Direct numerical simulations（DNSs） of purely elastic turbulence in rectilinear shear flows in a three-dimensional（3D） parallel plate channel were carried out,by which numerical databases were established.Based on the numerical databases,the present paper analyzed the structural and statistical characteristics of the elastic turbulence including flow patterns,the wall effect on the turbulent kinetic energy spectrum,and the local relationship between the flow motion and the microstructures＇ behavior.Moreover,to address the underlying physical mechanism of elastic turbulence,its generation was presented in terms of the global energy budget.The results showed that the flow structures in elastic turbulence were 3D with spatial scales on the order of the geometrical characteristic length,and vortex tubes were more likely to be embedded in the regions where the polymers were strongly stretched.In addition,the patterns of microstructures＇ elongation behave like a filament.From the results of the turbulent kinetic energy budget,it was found that the continuous energy releasing from the polymers into the main flow was the main source of the generation and maintenance of the elastic turbulent status.展开更多
In this paper, we investigate the performance of the bulk fin field effect transistor （FinFET） through a three- dimensional （3D） full band Monte Carlo simulator with quantum correction. Several scattering mechanis...In this paper, we investigate the performance of the bulk fin field effect transistor （FinFET） through a three- dimensional （3D） full band Monte Carlo simulator with quantum correction. Several scattering mechanisms, such as the acoustic and optical phonon scattering, the ionized impurity scattering, the impact ionization scattering and the surface roughness scattering are considered in our simulator. The effects of the substrate bias and the surface roughness scattering near the Si/SiO2 interface on the performance of bulk FinFET are mainly discussed in our work. Our results show that the on-current of bulk FinFET is sensitive to the surface roughness and that we can reduce the substrate leakage current by modulating the substrate bias voltage.展开更多
Three-dimensional simulations of ferroelectric hysteresis and butterfly loops are carried out based on solving the time dependent Ginzburg-Landau equations using a finite volume method. The influence of externally mec...Three-dimensional simulations of ferroelectric hysteresis and butterfly loops are carried out based on solving the time dependent Ginzburg-Landau equations using a finite volume method. The influence of externally mechanical loadings with a tensile strain and a compressive strain on the hysteresis and butterfly loops is studied numerically. Different from the traditional finite element and finite difference methods, the finite volume method is applicable to simulate the ferroelectric phase transitions and properties of ferroelectric materials even for more realistic and physical problems.展开更多
Supersonic axisymmetric jet flow over a missile afterbody containing exhaust jet is simulated using the second order accurate positive schemes method developed for solving the axisymmetric Euler equations based on the...Supersonic axisymmetric jet flow over a missile afterbody containing exhaust jet is simulated using the second order accurate positive schemes method developed for solving the axisymmetric Euler equations based on the 2-D conservation laws.Comparisons between the numerical results and the experimental measurements show excellent agreements.The computed results are in good agreement with the numerical solutions obtained by using third order accurate RKDG finite element method.The results show larger gradient at discontinuous points compared with those obtained by second order accurate TVD schemes.It indicates that the presented method is efficient and reliable for solving the axisymmetric jet with external freestream flows,and shows that the method captures shocks well without numerical noise.展开更多
The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simu...The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simulate the subgrid eddy viscosity, and the eddy-break-up (EBU) combustion subgrid-scale model is used to determine the chemical reaction rate. A two-step turbulent combustion subgrid-scale model is employed for calculating carbon monoxide CO concentration, and the NO subgrid-scale pollutant formation model for the evaluation of the rate of NO formation. The heat flux model is applied to the prediction of radiant heat transfer. The gas phase is solved with the SIMPLE algorithm and a hybrid scheme in the staggered grid system. The liquid phase equations are solved in a Lagrangian frame in reference of the particle-source-in-cell (PSIC) algorithm. From simulation results, the exchange of mass, moment and energy between gas and particle fields for the reacting flow in the afterburner with a V-gutter flame holder can be obtained. By the comparison of experimental and simulation results, profile temperature and pollutant of the outlet are quite in agreement with experimental data. Results show that the LES approach for predicting the two-phase instantaneous reacting flow and pollutant emissions in the afterburner is feasible.展开更多
文摘The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.
基金Financial support for this work was obtained from the National Natural Science Foundation of China(No.51074059)
文摘The present study examines the results of the researches related to the gob bulking factor carried out at home and abroad.A mathematical function of a three-dimensional gob bulking factor is described based on a three-dimensional gob model.The method of taking value for interstice and permeability ratios is also proposed.The law of air leakage of fully mechanized top coal is researched in this study.The results show that the speed of air flow near the upper and lower crossheadings is higher than that in the central section of the gob at the same distance from the working face.When the amount of air at the working face exceeds a critical amount,the width of the spontaneous combustion zone in the upper and lower crossheadings is also larger than that in the central section.In this situation,the key is preventing the coal left in the upper and lower crossheadings from self-igniting.Reducing the amount of air at the working face can decrease the width of the spontaneous combustion zone,especially the width near the upper and lower crossheadings.This also moves the spontaneous combustion zone in the direction of the working face.It can prevent the coal in the gob from self-igniting by making the coal left in the crossheadings to be inert and by effectively controlling the amount of air at the working face.
基金financially supported by the National Key Research and Development Plan(No.2018YFB0605601)the National Natural Science Foundation of China(No.41972168)。
文摘The three-dimensional(3 D) structures of pores directly affect the CH4 flow.Therefore,it is very important to analyze the3 D spatial structure of pores and to simulate the CH4 flow with the connected pores as the carrier.The result shows that the equivalent radius of pores and throats are 1-16 μm and 1.03-8.9 μm,respectively,and the throat length is 3.28-231.25 μm.The coordination number of pores concentrates around three,and the intersection point between the connectivity function and the X-axis is 3-4 μm,which indicate the macro-pores have good connectivity.During the single-channel flow,the pressure decreases along the direction of CH4 flow,and the flow velocity of CH4 decreases from the pore center to the wall.Under the dual-channel and the multi-channel flows,the pressure also decreases along the CH4 flow direction,while the velocity increases.The mean flow pressure gradually decreases with the increase of the distance from the inlet slice.The change of mean flow pressure is relatively stable in the direction horizontal to the bedding plane,while it is relatively large in the direction perpendicular to the bedding plane.The mean flow velocity in the direction horizontal to the bedding plane(Y-axis) is the largest,followed by that in the direction horizontal to the bedding plane(X-axis),and the mean flow velocity in the direction perpendicular to the bedding plane is the smallest.
文摘A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept (EDC) model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.
基金Supported by the National Natural Science Foundation of China(52374043)Key Program of the National Natural Science Foundation of China(52234003).
文摘Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.
基金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.
基金Supported by the Strategic Priority Research Program of the Chinese Academy of Sciences under Grant No XDA09020402the National Integrate Circuit Research Program of China under Grant No 2009ZX02023-003+1 种基金the National Natural Science Foundation of China under Grant Nos 61261160500,61376006,61401444 and 61504157the Science and Technology Council of Shanghai under Grant Nos 14DZ2294900,15DZ2270900 and 14ZR1447500
文摘An optimized device structure for reducing the RESET current of phase-change random access memory （PCRAM） with blade-type like （BTL） phase change layer is proposed. The electrical thermal analysis of the BTL cell and the blade heater contactor structure by three-dimensional finite element modeling are compared with each other during RESET operation. The simulation results show that the programming region of the phase change layer in the BTL cell is much smaller, and thermal electrical distributions of the BTL cell are more concentrated on the TiN/GST interface. The results indicate that the BTL cell has the superiorities of increasing the heating efficiency, decreasing the power consumption and reducing the RESET current from 0.67mA to 0.32mA. Therefore, the BTL cell will be appropriate for high performance PCRAM device with lower power consumption and lower RESET current.
基金Supported by the National Natural Science Foundation of China(51827804)CNPC Strategic Cooperation Science and Technology Major Project(ZLZX2020-01-05)Open Fund of State Key Laboratory of Rock Mechanics and Engineering(SKLGME021024).
文摘A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.
基金Supported by the National Defense Basic Scientific Research Project(A1020131011)
文摘A method of constructing three-dimensional process model for the punching cartridge cases is presented based on DEFORM simulation analysis. Using DEFORM software,the finite element simulation models for the punching and forming process of cartridge cases are established,and the corresponding simulation result model of each intermediate procedure is obtained by continuously performing the forming process simulation. The simulation model cannot annotate size and process information due to poor interface between DEFORM software and CAD software. Thus,a 3D annotation module is developed with secondary development technology of UG NX software. Consequently,the final process model with dimension and process information is obtained. Then,with the current 3D process management system,the 3D punching and forming process design of cartridge cases can be completed further. An example is also provided to illustrate that the relative error between the simulation process model and the physical model is less than 2%,which proves the validity and reliability of the proposed method in this study.
基金Sponsored by the Ministerial Level Foundation ( A12050914)the Excellent Young Teacher Foundation of Beijing Institute of Technology(1040012040101)
文摘The temperature field, velocity field, as well as species distribution in three-dimensional space are successfully calculated by establishing three-dimensional geometry model and solving plasma jet-substrate interaction equations, optimized particle trajecory models, as well as particle-particle heat transfer equations in space. Under typical working conditions, the flying trajectories and distribution of ZrO2 ceramic particles and Ni metal particles are also simulated. Results show that, the plasma jet becomes wider near the substrate, and the stochastic trajectory model is preferable to simulate the turbulent diffusion effect of particles. In addition, Ni metal particles penetrate relatively more deeply than ZrO2 ceramic particles due to larger density.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51276046 and 51506037)the Foundation for Innovative Research Groups of the National Natural Science Foundation of China(Grant No.51421063)+2 种基金the China Postdoctoral Science Foundation(Grant No.2016M591526)the Heilongjiang Postdoctoral Fund,China(Grant No.LBH-Z15063)the China Postdoctoral International Exchange Program
文摘Direct numerical simulations（DNSs） of purely elastic turbulence in rectilinear shear flows in a three-dimensional（3D） parallel plate channel were carried out,by which numerical databases were established.Based on the numerical databases,the present paper analyzed the structural and statistical characteristics of the elastic turbulence including flow patterns,the wall effect on the turbulent kinetic energy spectrum,and the local relationship between the flow motion and the microstructures＇ behavior.Moreover,to address the underlying physical mechanism of elastic turbulence,its generation was presented in terms of the global energy budget.The results showed that the flow structures in elastic turbulence were 3D with spatial scales on the order of the geometrical characteristic length,and vortex tubes were more likely to be embedded in the regions where the polymers were strongly stretched.In addition,the patterns of microstructures＇ elongation behave like a filament.From the results of the turbulent kinetic energy budget,it was found that the continuous energy releasing from the polymers into the main flow was the main source of the generation and maintenance of the elastic turbulent status.
基金Project supported by the National Basic Research Program of China (Grant No. 2011CBA00604)
文摘In this paper, we investigate the performance of the bulk fin field effect transistor （FinFET） through a three- dimensional （3D） full band Monte Carlo simulator with quantum correction. Several scattering mechanisms, such as the acoustic and optical phonon scattering, the ionized impurity scattering, the impact ionization scattering and the surface roughness scattering are considered in our simulator. The effects of the substrate bias and the surface roughness scattering near the Si/SiO2 interface on the performance of bulk FinFET are mainly discussed in our work. Our results show that the on-current of bulk FinFET is sensitive to the surface roughness and that we can reduce the substrate leakage current by modulating the substrate bias voltage.
基金Supported by the Research Starting Funds for Imported Talents of Ningxia University under Grant No BQD2012011
文摘Three-dimensional simulations of ferroelectric hysteresis and butterfly loops are carried out based on solving the time dependent Ginzburg-Landau equations using a finite volume method. The influence of externally mechanical loadings with a tensile strain and a compressive strain on the hysteresis and butterfly loops is studied numerically. Different from the traditional finite element and finite difference methods, the finite volume method is applicable to simulate the ferroelectric phase transitions and properties of ferroelectric materials even for more realistic and physical problems.
基金Supported by the National Natural Defense Basic Scientific Research Program of China(A262006-1288)the Key Disciplines Program of Shanghai Municipal Commission of Education(J50501)~~
文摘Supersonic axisymmetric jet flow over a missile afterbody containing exhaust jet is simulated using the second order accurate positive schemes method developed for solving the axisymmetric Euler equations based on the 2-D conservation laws.Comparisons between the numerical results and the experimental measurements show excellent agreements.The computed results are in good agreement with the numerical solutions obtained by using third order accurate RKDG finite element method.The results show larger gradient at discontinuous points compared with those obtained by second order accurate TVD schemes.It indicates that the presented method is efficient and reliable for solving the axisymmetric jet with external freestream flows,and shows that the method captures shocks well without numerical noise.
文摘The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simulate the subgrid eddy viscosity, and the eddy-break-up (EBU) combustion subgrid-scale model is used to determine the chemical reaction rate. A two-step turbulent combustion subgrid-scale model is employed for calculating carbon monoxide CO concentration, and the NO subgrid-scale pollutant formation model for the evaluation of the rate of NO formation. The heat flux model is applied to the prediction of radiant heat transfer. The gas phase is solved with the SIMPLE algorithm and a hybrid scheme in the staggered grid system. The liquid phase equations are solved in a Lagrangian frame in reference of the particle-source-in-cell (PSIC) algorithm. From simulation results, the exchange of mass, moment and energy between gas and particle fields for the reacting flow in the afterburner with a V-gutter flame holder can be obtained. By the comparison of experimental and simulation results, profile temperature and pollutant of the outlet are quite in agreement with experimental data. Results show that the LES approach for predicting the two-phase instantaneous reacting flow and pollutant emissions in the afterburner is feasible.
