Organic solid waste(OSW)contains many renewable materials.The pyrolysis and gasification of OSW can realize resource utilization,and its products can be used for methanation reaction to produce synthetic natural gas i...Organic solid waste(OSW)contains many renewable materials.The pyrolysis and gasification of OSW can realize resource utilization,and its products can be used for methanation reaction to produce synthetic natural gas in the specific reactor.In order to understand the dynamic characteristics of the reactor,a three-dimensional numerical model has been established by the method of Computational Fluid Dynamics(CFD).Along the height of the reactor,the particle distribution in the bed becomes thinner and the mean solid volume fraction decreases from 4.18%to 0.37%.Meanwhile,the pressure fluctuation range decreased from 398.76 Pa at the entrance to a much lower value of 74.47 Pa at the exit.In this simulation,three parameters of gas inlet velocity,operating temperature and solid particle diameter are changed to explore their influences on gas-solid multiphase flow.The results show that gas velocity has a great influence on particle distribution.When the gas inlet velocity decreases from 6.51 to 1.98 m/s,the minimum height that particles can reach decreases from 169 to 100 mm.Additionally,as the operating temperature increases,the particle holdup inside the reactor changes from 0.843%to 0.700%.This indicates that the particle residence time reduces,which is not conducive to the follow-up reaction.Moreover,with the increase of particle size,the fluctuation range of the pressure at the bottom of the reactor increases,and its standard deviation increases from 55.34 to 1266.37 Pa.展开更多
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.展开更多
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.展开更多
Based on the particle-in-cell technology and the secondary electron emission theory, a three-dimensional simulation method for multipactor is presented in this paper. By combining the finite difference time domain met...Based on the particle-in-cell technology and the secondary electron emission theory, a three-dimensional simulation method for multipactor is presented in this paper. By combining the finite difference time domain method and the panicle tracing method, such an algorithm is self-consistent and accurate since the interaction between electromagnetic fields and particles is properly modeled. In the time domain aspect, the generation of multipactor can be easily visualized, which makes it possible to gain a deeper insight into the physical mechanism of this effect. In addition to the classic secondary electron emission model, the measured practical secondary electron yield is used, which increases the accuracy of the algorithm. In order to validate the method, the impedance transformer and ridge waveguide filter are studied. By analyzing the evolution of the secondaries obtained by our method, multipactor thresholds of these components are estimated, which show good agreement with the experimental results. Furthermore, the most sensitive positions where multipactor occurs are determined from the phase focusing phenomenon, which is very meaningful for multipactor analysis and design.展开更多
This paper presents a three-dimensional particle-in-cell (PIC) simulation of a Ka-band relativistic Cherenkov source with a slow wave structure (SWS) consisting of metal photonic band gap (PBG) structures. In th...This paper presents a three-dimensional particle-in-cell (PIC) simulation of a Ka-band relativistic Cherenkov source with a slow wave structure (SWS) consisting of metal photonic band gap (PBG) structures. In the simulation, a perfect match layer boundary is employed to absorb passing band modes supported by the PBG lattice with an artificial metal boundary. The simulated axial field distributions in the cross section and surface of the SWS demonstrate that the device operates in the vicinity of the π point of a TM01-1ike mode. The Fourier transformation spectra of the axial fields as functions of time and space show that only a single frequency appears at 36.27 GHz, which is in good agreement with that of the intersection of the dispersion curve with the slow space charge wave generated on the beam. The simulation results demonstrate that the SWS has good mode selectivity.展开更多
The fabrication process dependent effects on single event effects (SEEs) are investigated in a commercial silicon- germanium heterojunction bipolar transistor (SiGe HBT) using three-dimensional (3D) TCAD simulat...The fabrication process dependent effects on single event effects (SEEs) are investigated in a commercial silicon- germanium heterojunction bipolar transistor (SiGe HBT) using three-dimensional (3D) TCAD simulations. The influences of device structure and doping concentration on SEEs are discussed via analysis of current transient and charge collection induced by ions strike. The results show that the SEEs representation of current transient is different from representation of the charge collection for the same process parameters. To be specific, the area of C/S junction is the key parameter that affects charge collection of SEE. Both current transient and charge collection are dependent on the doping of collector and substrate. The base doping slightly influences transient currents of base, emitter, and collector terminals. However, the SEEs of SiGe HBT are hardly affected by the doping of epitaxial base and the content of Ge.展开更多
This paper describes the method of constructing three dimensional model of bioactive ceramic to simulate the degrading process in vivo. Based on the construction of three dimensional simulation model and realizing it ...This paper describes the method of constructing three dimensional model of bioactive ceramic to simulate the degrading process in vivo. Based on the construction of three dimensional simulation model and realizing it by computer programming, it proves that the three dimensional stimulant model can correctly represent the real degrading process, and thus verifies its feasibility.展开更多
Particles and fields represent two major modeling paradigms in pure and applied science at all. In this paper a methodology and some of the results for three-dimensional (3D) simulations that include both field and pa...Particles and fields represent two major modeling paradigms in pure and applied science at all. In this paper a methodology and some of the results for three-dimensional (3D) simulations that include both field and particle abstractions are presented. Electromagnetic field calculations used here are based on the discrete differential form representation of the finite elements method, while the Monte Carlo method makes foundation of the particle part of the simulations. The first example is the simulation of the feature profile evolution during SiO2 etching enhanced by Ar + /CF4 non-equilibrium plasma based on the sparse field method for solving level set equations. Second example is devoted to the design of a spiral inflector which is one of the key devices of the axial injection system of the VINCY Cyclotron.展开更多
The difference of sintering crunodes of metal powders and fibers is discussed. The mathematical model of the surface diffusion described by the difference in mean curvature is defined as a Hamilton-Jacobi-type equatio...The difference of sintering crunodes of metal powders and fibers is discussed. The mathematical model of the surface diffusion described by the difference in mean curvature is defined as a Hamilton-Jacobi-type equation, and the model is numerically solved by the level set method. The three-dimensional numerical simulations of two metal powders and fibers(the fiber angle is 0° or 90°) are implemented by this mathematical model, respectively. The numerical simulation results accord with the experimental ones. The sintering neck growth trends of metal powders and metal fibers are similar. The sintering neck radius of metal fibers is larger than that of metal powders. The difference of the neck radius is caused by the difference of geometric structure which makes an important influence on the curvature affecting the migration rate of atoms.展开更多
The geometry of joints has a significant influence on the mechanical properties of rocks.To simplify the curved joint shapes in rocks,the joint shape is usually treated as straight lines or planes in most laboratory e...The geometry of joints has a significant influence on the mechanical properties of rocks.To simplify the curved joint shapes in rocks,the joint shape is usually treated as straight lines or planes in most laboratory experiments and numerical simulations.In this study,the computerized tomography (CT) scanning and photogrammetry were employed to obtain the internal and surface joint structures of a limestone sample,respectively.To describe the joint geometry,the edge detection algorithms and a three-dimensional (3D) matrix mapping method were applied to reconstruct CT-based and photogrammetry-based jointed rock models.For comparison tests,the numerical uniaxial compression tests were conducted on an intact rock sample and a sample with a joint simplified to a plane using the parallel computing method.The results indicate that the mechanical characteristics and failure process of jointed rocks are significantly affected by the geometry of joints.The presence of joints reduces the uniaxial compressive strength (UCS),elastic modulus,and released acoustic emission (AE) energy of rocks by 37%–67%,21%–24%,and 52%–90%,respectively.Compared to the simplified joint sample,the proposed photogrammetry-based numerical model makes the most of the limited geometry information of joints.The UCS,accumulative released AE energy,and elastic modulus of the photogrammetry-based sample were found to be very close to those of the CT-based sample.The UCS value of the simplified joint sample (i.e.38.5 MPa) is much lower than that of the CT-based sample (i.e.72.3 MPa).Additionally,the accumulative released AE energy observed in the simplified joint sample is 3.899 times lower than that observed in the CT-based sample.CT scanning provides a reliable means to visualize the joints in rocks,which can be used to verify the reliability of photogrammetry techniques.The application of the photogrammetry-based sample enables detailed analysis for estimating the mechanical properties of jointed rocks.展开更多
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.展开更多
Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different ...Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different computing tools have to be developed so as to solve particular fields at different scales and for different processes.Therefore,the integration of different types of software is inevitable.However,it is difficult to perform the transfer of the meshes and simulated results among software packages because of the lack of shared data formats or encrypted data formats.An image processing based method for three-dimensional model reconstruction for numerical simulation was proposed,which presents a solution to the integration problem by a series of slice or projection images obtained by the post-processing modules of the numerical simulation software.By means of mapping image pixels to meshes of either finite difference or finite element models,the geometry contour can be extracted to export the stereolithography model.The values of results,represented by color,can be deduced and assigned to the meshes.All the models with data can be directly or indirectly integrated into other software as a continued or new numerical simulation.The three-dimensional reconstruction method has been validated in numerical simulation of castings and case studies were provided in this study.展开更多
Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed bounda...Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.展开更多
Aerodynamic characteristic of the beetle, Trypoxylus dichotornus, which has a pair of elytra (forewings) and hind wings, is numerically investigated. Based on the experimental results of wing kinematics, two-dimensi...Aerodynamic characteristic of the beetle, Trypoxylus dichotornus, which has a pair of elytra (forewings) and hind wings, is numerically investigated. Based on the experimental results of wing kinematics, two-dimensional (2D) and three-dimensional (3D) computational fluid dynamic simulations were carried out to reveal aerodynamic performance of the hind wing. The roles of the spiral Leading Edge Vortex (LEV) and the spanwise flow were clarified by comparing 2D and 3D simulations. Mainly due to pitching down of chord line during downstroke in highly inclined stroke plane, relatively high averaged thrust was produced in the free forward flight of the beetle. The effects of the local corrugation and the camber variation were also investigated for the beetle's hind wings. Our results show that the camber variation plays a significant role in improving both lift and thrust in the flapping. On the other hand, the local corrugation pattern has no significant effect on the aerodynamic force due to large angle of attack during flapping.展开更多
Carbon dioxide(CO2),the main gas emitted from fossil burning,is the primary contributor to global warming.Circulating fluidized bed reactor(CFBR)is proved as an energy-efficient method for post-combustion CO2 capture....Carbon dioxide(CO2),the main gas emitted from fossil burning,is the primary contributor to global warming.Circulating fluidized bed reactor(CFBR)is proved as an energy-efficient method for post-combustion CO2 capture.The numerical simulation by computational fluid dynamics(CFD)is believed as a promising tool to study CO2 adsorption process in CFBR.Although three-dimensional(3D)simulations were proved to have better predicting performance with the experimental results,two-dimensional(2D)simulations have been widely reported for qualitative and quantitative studies on gas-solid behavior in CFBR for its higher computational efficiency recently.However,the discrepancies between 2D and 3D simulations have rarely been evaluated by detailed study.Considering that the differences between the 2D and 3D simulations will vary substantially with the changes of independent operating conditions,it is beneficial to lower computational costs to clarify the effects of dimensionality on the numerical CO2 adsorption runs under various operating conditions.In this work,the comparative analysis for CO2 adsorption in 2D and 3D simulations was conducted to enlighten the effects of dimensionality on the hydrodynamics and reaction behaviors,in which the separation rate,species distribution and hydrodynamic characteristics were comparatively studied for both model frames.With both accuracy and computational costs considered,the viable suggestions were provided in selecting appropriate model frame for the studies on optimization of operating conditions,which directly affect the capture and energy efficiencies of cyclic CO2 capture process in CFBR.展开更多
The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameteri...The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.展开更多
Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging ...Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.展开更多
基金Funding Statement:This work was supported by the National Key Research and Development Program of China[Grant No.2019YFC1906802].
文摘Organic solid waste(OSW)contains many renewable materials.The pyrolysis and gasification of OSW can realize resource utilization,and its products can be used for methanation reaction to produce synthetic natural gas in the specific reactor.In order to understand the dynamic characteristics of the reactor,a three-dimensional numerical model has been established by the method of Computational Fluid Dynamics(CFD).Along the height of the reactor,the particle distribution in the bed becomes thinner and the mean solid volume fraction decreases from 4.18%to 0.37%.Meanwhile,the pressure fluctuation range decreased from 398.76 Pa at the entrance to a much lower value of 74.47 Pa at the exit.In this simulation,three parameters of gas inlet velocity,operating temperature and solid particle diameter are changed to explore their influences on gas-solid multiphase flow.The results show that gas velocity has a great influence on particle distribution.When the gas inlet velocity decreases from 6.51 to 1.98 m/s,the minimum height that particles can reach decreases from 169 to 100 mm.Additionally,as the operating temperature increases,the particle holdup inside the reactor changes from 0.843%to 0.700%.This indicates that the particle residence time reduces,which is not conducive to the follow-up reaction.Moreover,with the increase of particle size,the fluctuation range of the pressure at the bottom of the reactor increases,and its standard deviation increases from 55.34 to 1266.37 Pa.
基金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.
基金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 Key Laboratory Foundation,China(Grant No.9140C530103110C5301)
文摘Based on the particle-in-cell technology and the secondary electron emission theory, a three-dimensional simulation method for multipactor is presented in this paper. By combining the finite difference time domain method and the panicle tracing method, such an algorithm is self-consistent and accurate since the interaction between electromagnetic fields and particles is properly modeled. In the time domain aspect, the generation of multipactor can be easily visualized, which makes it possible to gain a deeper insight into the physical mechanism of this effect. In addition to the classic secondary electron emission model, the measured practical secondary electron yield is used, which increases the accuracy of the algorithm. In order to validate the method, the impedance transformer and ridge waveguide filter are studied. By analyzing the evolution of the secondaries obtained by our method, multipactor thresholds of these components are estimated, which show good agreement with the experimental results. Furthermore, the most sensitive positions where multipactor occurs are determined from the phase focusing phenomenon, which is very meaningful for multipactor analysis and design.
基金Project supported by the National Key Basic Research Program of China (Grant No 2007CB31040)the National Natural Science Foundation of China (Grant No 60571020)
文摘This paper presents a three-dimensional particle-in-cell (PIC) simulation of a Ka-band relativistic Cherenkov source with a slow wave structure (SWS) consisting of metal photonic band gap (PBG) structures. In the simulation, a perfect match layer boundary is employed to absorb passing band modes supported by the PBG lattice with an artificial metal boundary. The simulated axial field distributions in the cross section and surface of the SWS demonstrate that the device operates in the vicinity of the π point of a TM01-1ike mode. The Fourier transformation spectra of the axial fields as functions of time and space show that only a single frequency appears at 36.27 GHz, which is in good agreement with that of the intersection of the dispersion curve with the slow space charge wave generated on the beam. The simulation results demonstrate that the SWS has good mode selectivity.
基金supported by the National Natural Science Foundation of China(Grant Nos.61274106,11175138,and 61601352)
文摘The fabrication process dependent effects on single event effects (SEEs) are investigated in a commercial silicon- germanium heterojunction bipolar transistor (SiGe HBT) using three-dimensional (3D) TCAD simulations. The influences of device structure and doping concentration on SEEs are discussed via analysis of current transient and charge collection induced by ions strike. The results show that the SEEs representation of current transient is different from representation of the charge collection for the same process parameters. To be specific, the area of C/S junction is the key parameter that affects charge collection of SEE. Both current transient and charge collection are dependent on the doping of collector and substrate. The base doping slightly influences transient currents of base, emitter, and collector terminals. However, the SEEs of SiGe HBT are hardly affected by the doping of epitaxial base and the content of Ge.
文摘This paper describes the method of constructing three dimensional model of bioactive ceramic to simulate the degrading process in vivo. Based on the construction of three dimensional simulation model and realizing it by computer programming, it proves that the three dimensional stimulant model can correctly represent the real degrading process, and thus verifies its feasibility.
基金supported by O171037,Ⅲ 41011 and Ⅲ45006 Projects of Ministry of Education and Science,Serbia.
文摘Particles and fields represent two major modeling paradigms in pure and applied science at all. In this paper a methodology and some of the results for three-dimensional (3D) simulations that include both field and particle abstractions are presented. Electromagnetic field calculations used here are based on the discrete differential form representation of the finite elements method, while the Monte Carlo method makes foundation of the particle part of the simulations. The first example is the simulation of the feature profile evolution during SiO2 etching enhanced by Ar + /CF4 non-equilibrium plasma based on the sparse field method for solving level set equations. Second example is devoted to the design of a spiral inflector which is one of the key devices of the axial injection system of the VINCY Cyclotron.
基金Projects(51174236,51134003)supported by the National Natural Science Foundation of ChinaProject(2011CB606306)supported by the National Basic Research Program of ChinaProject(PMM-SKL-4-2012)supported by the Opening Project of State Key Laboratory of Porous Metal Materials(Northwest Institute for Nonferrous Metal Research),China
文摘The difference of sintering crunodes of metal powders and fibers is discussed. The mathematical model of the surface diffusion described by the difference in mean curvature is defined as a Hamilton-Jacobi-type equation, and the model is numerically solved by the level set method. The three-dimensional numerical simulations of two metal powders and fibers(the fiber angle is 0° or 90°) are implemented by this mathematical model, respectively. The numerical simulation results accord with the experimental ones. The sintering neck growth trends of metal powders and metal fibers are similar. The sintering neck radius of metal fibers is larger than that of metal powders. The difference of the neck radius is caused by the difference of geometric structure which makes an important influence on the curvature affecting the migration rate of atoms.
基金supported by the National Natural Science Foundation of China(Grant Nos.42277150,41977219)Henan Provincial Science and Technology Research Project(Grant No.222102320271).
文摘The geometry of joints has a significant influence on the mechanical properties of rocks.To simplify the curved joint shapes in rocks,the joint shape is usually treated as straight lines or planes in most laboratory experiments and numerical simulations.In this study,the computerized tomography (CT) scanning and photogrammetry were employed to obtain the internal and surface joint structures of a limestone sample,respectively.To describe the joint geometry,the edge detection algorithms and a three-dimensional (3D) matrix mapping method were applied to reconstruct CT-based and photogrammetry-based jointed rock models.For comparison tests,the numerical uniaxial compression tests were conducted on an intact rock sample and a sample with a joint simplified to a plane using the parallel computing method.The results indicate that the mechanical characteristics and failure process of jointed rocks are significantly affected by the geometry of joints.The presence of joints reduces the uniaxial compressive strength (UCS),elastic modulus,and released acoustic emission (AE) energy of rocks by 37%–67%,21%–24%,and 52%–90%,respectively.Compared to the simplified joint sample,the proposed photogrammetry-based numerical model makes the most of the limited geometry information of joints.The UCS,accumulative released AE energy,and elastic modulus of the photogrammetry-based sample were found to be very close to those of the CT-based sample.The UCS value of the simplified joint sample (i.e.38.5 MPa) is much lower than that of the CT-based sample (i.e.72.3 MPa).Additionally,the accumulative released AE energy observed in the simplified joint sample is 3.899 times lower than that observed in the CT-based sample.CT scanning provides a reliable means to visualize the joints in rocks,which can be used to verify the reliability of photogrammetry techniques.The application of the photogrammetry-based sample enables detailed analysis for estimating the mechanical properties of jointed rocks.
基金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.
基金funded by National Key R&D Program of China(No.2021YFB3401200)the National Natural Science Foundation of China(No.51875308)the Beijing Nature Sciences Fund-Haidian Originality Cooperation Project(L212002).
文摘Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different computing tools have to be developed so as to solve particular fields at different scales and for different processes.Therefore,the integration of different types of software is inevitable.However,it is difficult to perform the transfer of the meshes and simulated results among software packages because of the lack of shared data formats or encrypted data formats.An image processing based method for three-dimensional model reconstruction for numerical simulation was proposed,which presents a solution to the integration problem by a series of slice or projection images obtained by the post-processing modules of the numerical simulation software.By means of mapping image pixels to meshes of either finite difference or finite element models,the geometry contour can be extracted to export the stereolithography model.The values of results,represented by color,can be deduced and assigned to the meshes.All the models with data can be directly or indirectly integrated into other software as a continued or new numerical simulation.The three-dimensional reconstruction method has been validated in numerical simulation of castings and case studies were provided in this study.
基金financially supported by the National Key R&D Program of China (Grant No.2022YFB2603000)the National Natural Science Foundation of China (Grant Nos.51779172 and 52179076)。
文摘Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.
文摘Aerodynamic characteristic of the beetle, Trypoxylus dichotornus, which has a pair of elytra (forewings) and hind wings, is numerically investigated. Based on the experimental results of wing kinematics, two-dimensional (2D) and three-dimensional (3D) computational fluid dynamic simulations were carried out to reveal aerodynamic performance of the hind wing. The roles of the spiral Leading Edge Vortex (LEV) and the spanwise flow were clarified by comparing 2D and 3D simulations. Mainly due to pitching down of chord line during downstroke in highly inclined stroke plane, relatively high averaged thrust was produced in the free forward flight of the beetle. The effects of the local corrugation and the camber variation were also investigated for the beetle's hind wings. Our results show that the camber variation plays a significant role in improving both lift and thrust in the flapping. On the other hand, the local corrugation pattern has no significant effect on the aerodynamic force due to large angle of attack during flapping.
基金supported by the National Natural Science Foundation of China(21506181,21506179)Natural Science Foundation of Hunan Province(2020JJ3033,2019JJ40281,2018SK2027,2018RS3088,2019SK2112)+1 种基金Research Foundation of Education Bureau of Hunan Province(18B088)Hunan Key Laboratory of Environment Friendly Chemical Process Integration and Hunan 2011 Collaborative Innovation Center of Chemical Engineering&Technology with Environmental Benignity and Effective Resource Utilization,State Key Laboratory of High-efficiency Utilization of Coal and Green Chemical Engineering(2020-KF-11).
文摘Carbon dioxide(CO2),the main gas emitted from fossil burning,is the primary contributor to global warming.Circulating fluidized bed reactor(CFBR)is proved as an energy-efficient method for post-combustion CO2 capture.The numerical simulation by computational fluid dynamics(CFD)is believed as a promising tool to study CO2 adsorption process in CFBR.Although three-dimensional(3D)simulations were proved to have better predicting performance with the experimental results,two-dimensional(2D)simulations have been widely reported for qualitative and quantitative studies on gas-solid behavior in CFBR for its higher computational efficiency recently.However,the discrepancies between 2D and 3D simulations have rarely been evaluated by detailed study.Considering that the differences between the 2D and 3D simulations will vary substantially with the changes of independent operating conditions,it is beneficial to lower computational costs to clarify the effects of dimensionality on the numerical CO2 adsorption runs under various operating conditions.In this work,the comparative analysis for CO2 adsorption in 2D and 3D simulations was conducted to enlighten the effects of dimensionality on the hydrodynamics and reaction behaviors,in which the separation rate,species distribution and hydrodynamic characteristics were comparatively studied for both model frames.With both accuracy and computational costs considered,the viable suggestions were provided in selecting appropriate model frame for the studies on optimization of operating conditions,which directly affect the capture and energy efficiencies of cyclic CO2 capture process in CFBR.
基金supported by the National Natural Science Foundation of China(Grant Nos.42175099,42027804,42075073)the Innovative Project of Postgraduates in Jiangsu Province in 2023(Grant No.KYCX23_1319)+3 种基金supported by the National Natural Science Foundation of China(Grant No.42205080)the Natural Science Foundation of Sichuan(Grant No.2023YFS0442)the Research Fund of Civil Aviation Flight University of China(Grant No.J2022-037)supported by the National Key Scientific and Technological Infrastructure project“Earth System Science Numerical Simulator Facility”(Earth Lab)。
文摘The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.
基金supported by the National Natural Science Foundation of China(NNSFC)grants 42074202,42274196Strategic Priority Research Program of Chinese Academy of Sciences grant XDB41000000ISSI-BJ International Team Interaction between magnetic reconnection and turbulence:From the Sun to the Earth。
文摘Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.